Magnetic field measurement method and magnetic field measurement device for measuring and offsetting original magnetic field

Abstract

A light source unit irradiates a gas cell disposed in a measurement region with linearly polarized light in which the direction of travel is a z-axis direction and the vibration direction of an electric field is a y-axis direction. A polarimeter detects optical characteristics of light passing through the gas cell. A magnetic field generator applies an artificial magnetic field, capable of varying an x-axis component, a y-axis component, and a z-axis component, to the measurement region. A calculation control unit generates a plurality of artificial magnetic fields, calculates a magnetization value or a value corresponding to the magnetization value on the basis of the detection results of the polarimeter, and calculates an original magnetic field present in the measurement region, using an artificial magnetic field when the magnetization value or the value corresponding to the magnetization value satisfies a condition for external value.

Claims

1. A magnetic field measurement method of measuring a magnetic field of a measurement region in a magnetic field measurement device, the magnetic field measurement device including (i) a light source unit that emits linearly polarized light in which a first direction, a second direction and a third direction are orthogonal to each other, (ii) a medium, disposed in the measurement region, which changes optical characteristics of the linearly polarized light in accordance with a magnetic field, and that is irradiated along the third direction with the linearly polarized light of which a vibration direction of an electric field of the linearly polarized light is the second direction, (iii) an optical detector that detects the optical characteristics, and (iv) a magnetic field generator that applies an artificial magnetic field to the measurement region, the method comprising: generating and applying to the measurement region a plurality of artificial magnetic fields, obtained by changing an artificial magnetic field component in the third direction and one artificial magnetic field component in the first direction and the second direction, in the magnetic field generator, in a state where the other artificial magnetic field component in the first direction and the second direction is set to a fixed value in each of the plurality of artificial magnetic fields; calculating a magnetization value that is a component of a magnetization vector of the medium in the first direction or a value corresponding to the magnetization value based on detection results of the optical detector for each of the plurality of artificial magnetic fields; and calculating an original magnetic field present in the measurement region, based on a selected artificial magnetic field of the plurality of artificial magnetic fields when the magnetization value in the selected artificial magnetic field or the value corresponding to the magnetization value in the selected artificial magnetic field satisfies a condition for external value.

2. The magnetic field measurement method according to claim 1, wherein the calculating of the original magnetic field includes calculating using: a first magnetic field that is an artificial magnetic field when the magnetization value or value corresponding to the magnetization value satisfies a maximum value condition; and a second magnetic field that is an artificial magnetic field when the magnetization value or value corresponding to the magnetization value satisfies a minimum value condition.

3. The magnetic field measurement method according to claim 2, wherein the one artificial magnetic field component is an artificial magnetic field component in the first direction, and the other artificial magnetic field component is an artificial magnetic field component in the second direction, and the calculating of the original magnetic field includes calculating using the following Expressions (4) to (6): $\begin{array}{cc}{C}_x=-\frac{A_{\mathrm{py}}+A_{\mathrm{vy}}}{2}-A_{\mathrm{fx}}& \left(4\right)\\ C_y=\frac{A_{\mathrm{py}}-A_{\mathrm{vy}}}{2}& \left(5\right)\\ C_z=-\frac{A_{\mathrm{pz}}+A_{\mathrm{vz}}}{2}& \left(6\right)\end{array}$ where, A.sub.py is a component of the first magnetic field in the second direction, A.sub.vy is a component of the second magnetic field in the second direction, A.sub.pz is a component of the first magnetic field in the third direction, A.sub.vz is a component of the second magnetic field in the third direction, A.sub.fx is the fixed value, C.sub.x is a component of the original magnetic field in the first direction, C.sub.y is a component of the original magnetic field in the second direction, and C.sub.z is a component of the original magnetic field in the third direction.

4. The magnetic field measurement method according to claim 1, wherein the one artificial magnetic field component is an artificial magnetic field component in the second direction, and the other artificial magnetic field component is an artificial magnetic field component in the first direction, and the calculating of the original magnetic field includes calculating using the following Expressions (1) to (3): $\begin{array}{cc}{C}_x=-\frac{A_{\mathrm{px}}+A_{\mathrm{vx}}}{2}& \left(1\right)\\ C_y=\begin{array}{c}{A}_{\mathrm{px}}-A_{\mathrm{vx}}\\ 2\end{array}-A_{\mathrm{fy}}& \left(2\right)\\ C_z=\begin{array}{c}{A}_{\mathrm{pz}}+A_{\mathrm{vz}}\\ 2\end{array}& \left(3\right)\end{array}$ where, A.sub.px is a component of the first magnetic field in the first direction, A.sub.vx, is a component of the second magnetic field in the first direction, A.sub.pz is a component of the first magnetic field in the third direction, A.sub.vz is a component of the second magnetic field in the third direction, A.sub.fy is the fixed value, C.sub.x is a component of the original magnetic field in the first direction, C.sub.y is a component of the original magnetic field in the second direction, and C.sub.z is a component of the original magnetic field in the third direction.

5. The magnetic field measurement method according to claim 1, wherein the fixed value is zero.

6. The magnetic field measurement method according to claim 1, further comprising: generating a magnetic field of a difference in the original magnetic field with respect to a target magnetic field, in the magnetic field generator; disposing a measuring object in the measurement region; and measuring a magnetic field which is radiated by the measuring object, using the detection results of the optical detector while the magnetic field of the difference is generated.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) Various embodiments of the invention will be described with reference to the accompanying drawings, wherein like numbers reference like elements.

(2) FIG. 1 is a diagram illustrating an entire configuration example of a magnetic field measurement device according to Embodiment 1.

(3) FIG. 2 is a diagram illustrating an outline of an arrangement relation between a light source unit, a gas cell, and a polarimeter.

(4) FIG. 3 is a diagram illustrating a rotation of the plane of polarization.

(5) FIG. 4 is a diagram illustrating a relationship between alignment orientation angles and detection results of probe light.

(6) FIG. 5 is a diagram illustrating a distribution of spin polarization degrees M.sub.x.

(7) FIG. 6 is a diagram illustrating a positional relationship between a maximum point and a minimum point of the distribution of FIG. 5 in a plane B.sub.xB.sub.z.

(8) FIG. 7 is another diagram illustrating a distribution of spin polarization degrees M.sub.x.

(9) FIG. 8 is a diagram illustrating a positional relationship a maximum point and a minimum point of the distribution of FIG. 7 in a plane B.sub.xB.sub.z.

(10) FIG. 9 is a diagram illustrating a relationship between a coordinate system of a magnetic field of a measurement region of Embodiment 1 and a coordinate system of an artificial magnetic field.

(11) FIG. 10 is a flow diagram illustrating a procedure of a magnetic field forming process of Embodiment 1.

(12) FIG. 11 is a diagram illustrating an entire configuration example of a magnetic field measurement device of Embodiment 3.

(13) FIG. 12 is a diagram illustrating a three-dimensional distribution relating to B.sub.x and B.sub.z of a differential value M.sub.x/B.sub.z.

(14) FIG. 13 is a diagram illustrating a two-dimensional distribution relating to B.sub.z of the differential value M.sub.x/B.sub.z when B.sub.x=0.

(15) FIG. 14 is a diagram illustrating a relationship between a coordinate system of a magnetic field of a measurement region of Embodiment 3 and a coordinate system of an artificial magnetic field.

(16) FIG. 15 is a flow diagram illustrating a procedure of a magnetic field forming process of Embodiment 3.

DETAILED DESCRIPTION

(17) Hereinafter, one form for implementing a magnetic field measurement method and a magnetic field measurement device according to various embodiments of the invention will be described. Meanwhile, the invention is not limited by embodiments described below, and forms capable of applying the invention are not limited to the following embodiments. In addition, in the accompanying drawings, the same components are denoted by the same reference numerals and signs.

(18) Embodiment 1

(19) Entire Configuration

(20) FIG. 1 is a diagram illustrating an entire configuration example of a magnetic field measurement device 1 of the present embodiment. In addition, FIG. 2 is a diagram illustrating an outline of an arrangement relation between a light source unit 2, a gas cell 3, and a polarimeter 4 constituting the magnetic field measurement device 1. The magnetic field measurement device 1 of the present embodiment is used in a magnetocardiograph that measures a magnetocardiogram or a magnetoencephalography that measures a magnetoencephalo. The magnetic field measurement device 1 has a so-called one-beam-type magnetic sensor incorporated therein which uses both pump light irradiation and probe light irradiation, as an optically-pumped magnetic sensor, and measures a magnetic field using nonlinear magneto-optical rotation (NMOR). Meanwhile, without being limited to a one-beam type, a so-called two-beam-type configuration may be used in which a light source unit that performs irradiation with pump light and a light source unit that performs irradiation with probe light are separated from each other.

(21) As shown in FIG. 1, the magnetic field measurement device 1 includes the light source unit 2, the gas cell 3, the polarimeter 4 as an optical detector, a calculation control unit 5, and a magnetic field generator 7. In the present embodiment, as shown in FIG. 2, a third direction which is the direction of travel of linearly polarized light (irradiation light) in which irradiation with both the pump light and the probe light is performed by the light source unit 2 is defined as a z-axis direction, a second direction which is the vibration direction of a electric field of the linearly polarized light is defined as a y-axis direction, and a first direction perpendicular to the z-axis direction and the y-axis direction is defined as an x-axis direction. A space in which the light source unit 2, the gas cell 3, and the polarimeter 4 are disposed is denoted as an orthogonal three-axis xyz coordinate space.

(22) The light source unit 2 is constituted by a light source 21 and a polarizing plate 23, and emits linearly polarized light that propagates in the z-axis direction and vibrates along the y-axis direction, as irradiation light. The light source 21 is a laser generating device that generates a laser beam having a frequency corresponding to the transition of a hyperfine structure level of gaseous atoms sealed in the gas cell 3. Specifically, the wavelength of a laser beam is a wavelength equivalent to state transition between hyperfine structure quantum numbers F and F (=F1) of a line D1 of gaseous atoms (such as, for example, cesium, potassium, or rubidium) within the gas cell 3. The polarizing plate 23 is an element that polarizes a laser beam from the light source 21 in a predetermined direction and forms linearly polarized light. The irradiation light emitted from the light source unit 2 is guided by, for example, an optical fiber or the like, and the gas cell 3 is irradiated therewith.

(23) The gas cell 3 is an element made of glass in which alkali metal atoms such as potassium (K), rubidium (Rb), or cesium (Cs) are sealed in a gaseous state. The alkali metal atoms are excited (optically pumped) by the irradiation light (pump light and probe light) from the light source unit 2, and has a property as a medium for rotating the plane of polarization of light passing through the gas cell 3 in accordance with the intensity of a magnetic field. The gas cell 3 is disposed inside a magnetic shield 8 shown by a dashed-two dotted line in FIG. 1. The magnetic shield 8 is used for shielding above a certain magnetism and forming a space in which magnetism is reduced as compared to the outside of the magnetic shield 8, and the area of a subject to be tested as a measuring object such as the heart or the brain is located at a measurement region (peripheral region of the gas cell 3) in which the gas cell 3 is disposed inside the magnetic shield 8 during measurement. As described later, the measurement region can be set to a target magnetic field (for example, zero magnetic field) by the magnetic field generator 7 installed within the magnetic shield 8. The magnetic field measurement device 1 sets the magnetic field of a measurement region to be in a state of the target magnetic field, and then disposes the measurement area of a subject to be tested in the measurement region, to thereby measure a magnetic field which is radiated by the measurement area. Therefore, it can be also said that a magnetic field forming device is included in the magnetic field measurement device 1. Meanwhile, the gaseous atoms within the gas cell 3 may be set to be in a gaseous state during the measurement of a magnetic field, and may be set to be in a gaseous state at all times. In addition, the material of the gas cell 3 may be a material that transmits the irradiation light without being limited to glass, and may be a resin or the like.

(24) The polarimeter 4 is constituted by a polarization splitter 41 and two photodetectors 431 and 433, and is configured to split the irradiation light (probe light) passing through the gas cell 3 into two polarization components perpendicular to each other and detect each light intensity. The polarization splitter 41 is an element that splits the irradiation light from the gas cell 3 into components of an -axis and a -axis (see FIG. 3) perpendicular to each other in each axial direction. One polarization component which is split is guided to the photodetector 431, and the other polarization component is guided to the photodetector 433. The polarization splitter 41 is constituted by, for example, a Wollaston prism, a polarization beam splitter or the like. The photodetectors 431 and 433 receives the polarization components split by the polarization splitter 41, generates signals based on the amount of light received, and outputs the signals to a signal processing unit 51 of the calculation control unit 5.

(25) The calculation control unit 5 is configured using a microprocessor such as a CPU (Central Processing Unit), an ASIC (Application Specific Integrated Circuit), an IC (Integrated Circuit) memory, or the like, and controls an operation of each unit of the device as a whole. The calculation control unit 5 includes the signal processing unit 51, a magnetic field calculation unit 53, and a correction magnetic field setting unit 55, and performs magnetic field measurement in the measurement region, a magnetic field forming process (see FIG. 10) for forming a target magnetic field (for example, zero magnetic field) in the measurement region, or the like. In addition, the calculation control unit 5 includes a storage unit such as flash memory or a hard disk, and readably stores a program in the storage unit when the magnetic field forming process is realized by executing the program.

(26) In addition, an input unit 61 for inputting required information and a display unit 63 for displaying measurement results of a magnetic field or the like are appropriately connected to the calculation control unit 5. The input unit 61 is constituted by input devices of various types of switch such as a button switch, a lever switch, or a dial switch, a touch panel, a keyboard, a mouse, and the like. The display unit 63 is constituted by a display device such as an LCD (Liquid Crystal Display) or an EL display (Electroluminescence display).

(27) The signal processing unit 51 measures the magnetic field of the measurement region by calculating the rotation angle of the plane of polarization which is rotated by the probe light passing through the gas cell 3. The signal processing unit 51 processes signals from the photodetectors 431 and 433, and calculates a square sum W.sub.+ and a square difference W.sub. of the components of the -axis and the -axis in each axial direction in accordance with the following Expressions (7) and (8). E.sub. indicates the light intensity of the component in an -axis direction, and E.sub. indicates the light intensity of the component in a -axis direction.
W.sub.+=E.sub..sup.2+E.sub..sup.2(7)
W.sub.=E.sub..sup.2E.sub..sup.2(8)

(28) The magnetic field calculation unit 53 calculates a magnetic field present in the measurement region, using the value of an artificial magnetic field when the square difference W.sub. calculated by the signal processing unit 51 satisfies specified conditions such as a maximum value condition and a minimum value condition described later as a condition for external value. The magnetic field calculation unit 53 includes an artificial magnetic field control unit 531. The artificial magnetic field control unit 531 applies an artificial magnetic field to the measurement region by controlling the magnetic field generator 7, and the magnetic field calculation unit 53 measures a spin polarization degree M.sub.x in that case. For example, during the measurement of a magnetic field (referred to as an original magnetic field) present in the measurement region when the artificial magnetic field is set to zero, in a state where an artificial magnetic field component (referred to as a y-axis component) in the y-axis direction is set to a fixed value (predetermined value), a plurality of artificial magnetic fields obtained by changing a artificial magnetic field component (referred to as an x-axis component) in the x-axis direction and a artificial magnetic field component (referred to as a z-axis component) in the z-axis direction are sequentially generated in the measurement region, and in a state where one artificial magnetic field is applied, the spin polarization degree M.sub.x in that case is measured. In short, a first artificial magnetic field is generated, and the spin polarization degree M.sub.x when the first artificial magnetic field is applied is measured. Similarly, a second artificial magnetic field is generated, and the spin polarization degree M.sub.x when the second artificial magnetic field is applied is measured. Hereinafter, by repeating the same process, an N-th artificial magnetic field is generated (N is an integer equal to or greater than 2), and the spin polarization degree M.sub.x when the N-th artificial magnetic field is applied is measured. In the first artificial magnetic field to the N-th artificial magnetic field, the y-axis component is a fixed value and is common, but the x-axis component and the z-axis component are scanned, and the spin polarization degree M.sub.x is measured in each state. Meanwhile, the original magnetic field is typically an external magnetic field entering the magnetic shield 8 from the outside, but the original magnetic field itself may be present in a magnetic field to be measured, for example, such as magnetocardiogram or magnetoencephalo which is radiated by a living body.

(29) The correction magnetic field setting unit 55 sets a correction magnetic field for correcting the original magnetic field calculated by the magnetic field calculation unit 53, and forms a target magnetic field in the measurement region by controlling the magnetic field generator 7 to generate the correction magnetic field in the measurement region. When the original magnetic field is an external magnetic field, an example of correction of the original magnetic field is to offset the original magnetic field, and set the target magnetic field to a zero magnetic field.

(30) The magnetic field generator 7 is constituted by a three-axis Helmholtz coil for applying magnetic fields of the x-axis, the y-axis, and the z-axis in each axial direction, and includes a pair of coils disposed in pairs in each axial direction with the gas cell 3 interposed therebetween within the magnetic shield 8, and a current supply portion that supplies a current to these coils. The magnetic field generator 7 can generate a magnetic field in any three-dimensional direction in the measurement region.

(31) Meanwhile, the z-axis direction is a third direction in embodiments of the invention, but as shown in FIG. 2, the irradiation light (pump light) is not necessarily limited to be emitted from the light source unit 2 in the z-axis direction. The irradiation light (pump light) may be incident on the gas cell 3 in the z-axis direction after emission.

(32) Principle

(33) In the magnetic field measurement device 1 configured in this manner, when the gas cell 3 is irradiated with the pump light from the light source unit 2, gaseous atoms within the gas cell 3 are spin-polarized. The probability distribution of magnetic moments generated by this spin polarization when energy transitions from the hyperfine structure quantum number F to F (=F1) forms an ellipsoidal shape extending along the y-axis direction which is the vibration direction of the linearly polarized light. This biased probability distribution is referred to as alignment, and generating the alignment is referred to as optical pumping. When the measurement region is a zero magnetic field, the generated alignment remains along the y-axis direction which is the vibration direction of the pump light. However, when a magnetic field which is not zero is present in the measurement region, the alignment performs a precessional motion using the direction of a magnetic field as a rotational axis. As a result, the plane of polarization of the linearly polarized light rotates at an angle based on the magnetic field of the measurement region, using the z-axis direction which is the direction of travel as a rotational axis.

(34) FIG. 3 is a diagram illustrating a rotation of the plane of polarization. As described above, the alignment performs a precessional motion in accordance with the magnetic field (magnetic field received by the gas cell 3) of the measurement region. An optical pumping action caused by the pump light and a relaxation action caused by the collision or the like of the gaseous atoms with the inner wall of the gas cell 3 are applied, and thus the alignment is set to be in a steady state in an arrangement in which the alignment is rotated by an angle () based on the intensity of a magnetic field with respect to the y-axis, as shown in an ellipsoidal body hatched in FIG. 3.

(35) The probe light passing through the gas cell 3 receives an action of linear dichroism due to this alignment. The linear dichroism refers to a property in which transmissivities of the linearly polarized light different from each other a direction (direction of .sub.p) along the alignment and a direction (direction of .sub.s) perpendicular to the alignment. Specifically, since larger components in the direction perpendicular to the alignment are absorbed than those in the direction along the alignment, the plane of polarization of the probe light is rotated so as to approach the direction along the alignment.

(36) For example, in the present embodiment, the probe light incident on the gas cell 3 is linearly polarized light of a vector E.sub.0 in which the vibration direction of an electric field is the y-axis direction. Due to the alignment, an component of the probe light in a .sub.p direction passes through the gas cell 3 at transmissivity t.sub.p, and an component in a .sub.s direction passes through the gas cell 3 at transmissivity t.sub.s. Since the relation of t.sub.p>t.sub.s is established due to the linear dichroism, the vibration direction of an electric field of the probe light passing through the gas cell 3 is set to a vector E.sub.1. That is, the plane of polarization of the probe light passing through the gas cell 3 is rotated so as to approach the .sub.P direction. Specifically, when a vector of a component along the alignment of the vector E.sub.0 is denoted by E.sub.0P, a vector of a component along a direction perpendicular to the alignment of the vector E.sub.0 and the direction of travel of the linearly polarized light is denoted by E.sub.0S a vector of a component along of the alignment of the vector E.sub.1 is denoted by E.sub.1P, and a vector of a component along a direction perpendicular to the alignment of the vector E.sub.1 and the direction of travel of the linearly polarized light is denoted by E.sub.1s, a relationship between E.sub.1P=t.sub.pE.sub.0P and E.sub.1s=t.sub.sE.sub.0s is established.

(37) As shown in FIG. 3, when an angle (hereinafter, referred to as an alignment orientation angle) formed by the direction (.sub.p direction) along the alignment and the vibration direction (vector E.sub.0 direction) of an electric field of the probe light incident on the gas cell 3 is set to , each component of the vector E.sub.1 in the .sub.p direction and the .sub.s direction is calculated by the following Expression (9) from the above-mentioned relation, and each component of (,) in a coordinate system is calculated by the following Expression (10).

(38) $\begin{array}{cc}{\stackrel{.fwdarw.}{E}}_1=\left(0E_0\right)\left(\begin{array}{cc}\cos & \sin \\ -\sin & \cos \end{array}\right)\left(\begin{array}{cc}{t}_s& 0\\ 0& t_p\end{array}\right)& \left(9\right)\\ \left(E_{}{E}_{}\right)={\stackrel{.fwdarw.}{E}}_1\left(\begin{array}{cc}\cos \left(-\frac{}{4}-\right)& \sin \left(-\frac{}{4}-\right)\\ -\sin \left(-\frac{}{4}-\right)& \cos \left(-\frac{}{4}-\right)\end{array}\right)& \left(10\right)\end{array}$

(39) FIG. 4 is a diagram illustrating a relationship between the alignment orientation angle and the detection result of the probe light. In FIG. 4, when focused on the value of the square difference W.sub., the square difference W.sub. vibrates using 180 degrees as a period with respect to the alignment orientation angle . In the square difference W.sub., the alignment orientation angle substantially linearly changes with respect to the alignment orientation angle in the range of 45 degrees to +45 degrees, and thus high sensitivity is obtained. In addition, since the center of the linear change is 0 degrees, and the range of the linear change is wider than others (square sum W.sub.+ and the like), it is suitable to measure a magnetic field which is generated in the measurement region.

(40) A biomagnetic field such as magnetocardiogram or magnetoencephalo is weak, and the alignment orientation angle is small, and thus the rotation angle of the plane of polarization can be observed with high sensitivity when the square difference W.sub. is used. However, when an external magnetic field is present in the measurement region, as described above, sensitivity decreases due to the influence, and a decrease in the accuracy of measurement is caused. Generally, the measurement of a magnetic field to be measured such as magnetocardiogram or magnetoencephalo is performed under the environment where the infiltration of a magnetic field (external magnetic field) from the outside to the measurement region is suppressed by the magnetic shield 8 (state where the infiltration of an external magnetic field is low), but it is difficult to sufficiently reduce the external magnetic field to an extent of not influencing measurement through the magnetic shield 8, and thus the original magnetic field (that is, external magnetic field entering the measurement region) may be present in the measurement region. In other words, the infiltration of the external magnetic field is not able to be completely shielded by the magnetic shield 8, and thus the original magnetic field which is not zero may be present in the measurement region. A magnetic shield device capable of shielding magnetism completely is large in size, is expensive in cost, and is high in installation cost or operating cost. Consequently, in the present embodiment, the magnetic shield 8 is used, the original magnetic field within the magnetic shield 8 is then measured, and a magnetic field to be measured is measured in a state where the original magnetic field is reduced by the magnetic field generator 7. However, when the external magnetic field is originally low (when the external magnetic field and the original magnetic field are equal to each other and are small) or when the external magnetic field is stabilized, it is also possible to configure the present embodiment without using the magnetic shield 8.

(41) Since the square difference W.sub. is substantially proportional to the value of an x-axis component M.sub.x (hereinafter, denoted by the spin polarization degree M.sub.x) of spin polarization degrees (M.sub.x, M.sub.y, M.sub.z) of the alignment generated within the gas cell 3, the spin polarization degree M.sub.x is measured from the output value of the square difference W.sub.. The spin polarization degrees (M.sub.x, M.sub.y, M.sub.z) are equivalent to the magnetization vector of a medium (gas of alkali metal atoms) within the gas cell 3. That is, the square difference W.sub. is an example of a value corresponding to a magnetization value. Hereinafter, it is examined how the magnetization value (spin polarization degree M.sub.x) is changed by a change in the x-axis component, the y-axis component, and the z-axis component (absolute magnetic flux densities B.sub.x, B.sub.y, and B.sub.z) of the magnetic field of the measurement region, using the spin polarization degree M.sub.x as a magnetization value indicating an component (x-axis component) of the magnetization vector in the first direction.

(42) The temporal development of the spin polarization degrees (M.sub.x, M.sub.y, M.sub.z) of the alignment generated by optical pumping is approximated by Bloch equations represented by the following Expressions (11) to (13). Here, .sub.F indicates a gyromagnetic ratio which is determined by the type of medium gas (alkali metal atom gas) within the gas cell 3. In addition, .sub.0 indicates a relaxation rate of the spin polarization degrees (M.sub.x, M.sub.y, M.sub.z), and .sub.p indicates an optical pumping rate. The relaxation rate .sub.0 of the spin polarization degree and the optical pumping rate .sub.p are represented by the same unit system as that of an angular frequency, and specifically have a unit of radian per second (rad/s). Alternatively, a cut-off angular frequency .sub.C is a sum of the relaxation rate .sub.0 of the spin polarization degree and the optical pumping rate .sub.p (.sub.C=.sub.0+.sub.p).

(43) $\begin{array}{cc}\frac{{\mathrm{dM}}_x}{\mathrm{dt}}={}_F\left(M_y{B}_z-M_z{B}_y\right)-{}_0{M}_x-{}_p{M}_x& \left(11\right)\\ \frac{{\mathrm{dM}}_y}{\mathrm{dt}}={}_F\left(M_z{B}_x-M_x{B}_z\right)-{}_0{M}_y-{}_p\left(M_p-M_y\right)& \left(12\right)\\ \frac{{\mathrm{dM}}_z}{\mathrm{dt}}={}_F\left(M_x{B}_y-M_y{B}_x\right)-{}_0{M}_z-{}_p{M}_z& \left(13\right)\end{array}$

(44) Since the gas cell 3 is irradiated with the pump light and the probe light with steadily constant power, steady solutions of the spin polarization degrees (M.sub.x, M.sub.y, M.sub.z) can be solved by setting the left-hand sides of Expressions (11) to (13) to zero. The solutions are represented by the following Expression (14) to (17).

(45) $\begin{array}{cc}{M}_x=\left(\frac{c}{a}\right)\frac{B_x{B}_y+a{B}_z}{a^2+B_x^2+B_y^2+B_z^2}& \left(14\right)\\ M_y=\left(\frac{c}{a}\right)\frac{a^2+B_y^2}{a^2+B_x^2+B_y^2+B_z^2}& \left(15\right)\\ M_z=\left(\frac{c}{a}\right)\frac{B_y{B}_z-a{B}_x}{a^2+B_x^2+B_y^2+B_z^2}& \left(16\right)\\ a=\frac{{}_0+{}_p}{{}_F}c={}_p{M}_p& \left(17\right)\end{array}$

(46) First, FIG. 5 shows a distribution of the spin polarization degree M.sub.x in which the value of B.sub.y which is the y-axis component of the magnetic field of the measurement region is fixed to B.sub.y=0 [nT], and which is obtained by Expression (14) while changing the respective values of B.sub.x which is the x-axis component and B.sub.z which is the z-axis component. As shown in FIG. 5, a maximum point P and a minimum point V appear one by one in the distribution of the spin polarization degree M.sub.x. In addition, FIG. 6 is a diagram illustrating a positional relationship between the maximum point P and the minimum point V shown in FIG. 5 in a plane B.sub.xB.sub.z. As shown in FIG. 6, a straight line that links the maximum point P to the minimum point V is parallel to a B.sub.z axis when the plane B.sub.xB.sub.z is seen in plan view, and a midpoint F between the respective points P and V passes through the origin of the plane B.sub.xB.sub.z.

(47) Next, FIG. 7 shows a distribution of the spin polarization degree M.sub.x in which the value of B.sub.y which is the y-axis component of the magnetic field of the measurement region is fixed to B.sub.y=100 [nT], and which is obtained while changing the respective values of B.sub.x and B.sub.z in the same way as in FIG. 5. In addition, FIG. 8 is a diagram illustrating a positional relationship between a maximum point P and a minimum point V of the distribution of the spin polarization degree M.sub.x shown in FIG. 7 in the plane B.sub.xB.sub.z. As shown in FIG. 8, when the value of B.sub.y is 100 [nT], a straight line that links the maximum point P to the minimum point V is inclined with respect to the B.sub.z axis, but a characteristic in which the midpoint F between the respective points P and V passes through the origin of the plane B.sub.xB.sub.z is not different from the case of B.sub.y=0 [nT] shown in FIG. 5.

(48) When B.sub.y is any fixed value B.sub.f (B.sub.y=B.sub.f).sub.f a magnetic field P=(P.sub.x, P.sub.y, P.sub.z) in which the spin polarization degree M.sub.x is set to the maximum point P and a magnetic field V=(V.sub.x, V.sub.y, V.sub.z) in which the spin polarization degree M.sub.x is set to a minimum point V are obtained by differentiating M.sub.x of Expression (14) by B.sub.x, and solving each value obtained by differentiating M.sub.x of Expression (14) by B.sub.z in terms of B.sub.x and B.sub.z, setting each value to zero (the following Expression (18)).

(49) $\begin{array}{cc}\frac{M_x}{B_x}=0\frac{M_x}{B_z}=0& \left(18\right)\end{array}$

(50) The solutions are represented by the following Expressions (19) and (20), and can be obtained from the value of B.sub.y and a constant a shown in Expression (17). Here, a is a constant relating to relaxation.
{right arrow over (P)}=(P.sub.xP.sub.yP.sub.z)=(B.sub.yB.sub.ya)=(B.sub.fB.sub.fa)(19)
{right arrow over (V)}=(V.sub.xV.sub.yV.sub.z)=(B.sub.yB.sub.ya)=(B.sub.fB.sub.fa)(20)

(51) In addition, as shown in the following Expression (21), the midpoint F between the respective points P and V indicated by the solutions of Expressions (19) and (20) constantly satisfies the relation of B.sub.x=B.sub.z=0 regardless of the value of B.sub.y. meanwhile, in a case of B.sub.y=B.sub.f=0, the midpoint F is constantly coincident with the origin.

(52) $\begin{array}{cc}\stackrel{.fwdarw.}{F}=\frac{\stackrel{.fwdarw.}{P}+\stackrel{.fwdarw.}{V}}{2}=\left(\frac{B_y-B_y}{2}{B}_y\frac{a-a}{2}\right)=\left(0B_{f}0\right)& \left(21\right)\end{array}$

(53) Consequently, in the magnetic field forming process of the present embodiment, the magnetic field of the measurement region, that is, the magnetization value (spin polarization degree M.sub.x) is measured while the magnetic field generator 7 generates a plurality of artificial magnetic fields A obtained by changing the x-axis component (A.sub.x) and the z-axis component (A.sub.z), in a state where the y-axis component is first set to a fixed value (A.sub.y=A.sub.fy), on the basis of the characteristics (that is, characteristics of the square difference W.sub.) of the spin polarization degree M.sub.x described above. Specifically, the signal processing unit 51 calculates and stores the square difference W.sub. as the spin polarization degree M.sub.x whenever the artificial magnetic field A is applied, on the basis of the detection results of the polarimeter 4. The original magnetic field is then calculated using a first magnetic field (artificial magnetic field in which the magnetic field of the measurement region is set to the maximum point P) when a maximum value condition in which the square difference W.sub. is set to a maximum value is satisfied, and a second magnetic field (artificial magnetic field in which the magnetic field of the measurement region is set to the minimum point V) when a minimum value condition in which the square difference W.sub. is set to a minimum value is satisfied.

(54) FIG. 9 is a diagram illustrating a relationship between a coordinate system of a magnetic field B=(B.sub.x, B.sub.y, B.sub.z) of the measurement region and a coordinate system of an artificial magnetic field A=(A.sub.x, A.sub.y, A.sub.z) generated by the magnetic field generator 7, when the y-axis component of the artificial magnetic field is zero (A.sub.y=0). As shown in FIG. 9, the magnetic field B=(B.sub.x, B.sub.y, B.sub.z) of the measurement region can be calculated by performing vector addition of the artificial magnetic field A=(A.sub.x, A.sub.y, A.sub.z) to an original magnetic field C=(C.sub.x, C.sub.y, C.sub.z). A calculation expression is represented by the following Expression (22).
{right arrow over (B)}={right arrow over (A)}+{right arrow over (C)}=(B.sub.xB.sub.yB.sub.z)=(A.sub.x+C.sub.xA.sub.y+C.sub.yA.sub.z+C.sub.z)(22)

(55) Assuming that an artificial magnetic field when the magnetization value (spin polarization degree M.sub.x) in the measurement region becomes the maximum point P=(P.sub.x, P.sub.y, P.sub.z) is set to the first magnetic field A.sub.p=(A.sub.px, A.sub.py, A.sub.pz), and that an artificial magnetic field when the magnetization value (spin polarization degree M.sub.x) becomes the minimum point V=(V.sub.x, V.sub.y, V.sub.z) is set to the second magnetic field A.sub.v=(A.sub.vx, A.sub.vy, A.sub.vz), these artificial magnetic fields are represented by the following Expressions (23) and (24) from the relation of Expression (22) and a condition in which the y-axis component of the artificial magnetic field is set to a fixed value (A.sub.y=A.sub.fy).

(56) $\begin{array}{cc}\begin{array}{c}\left(P_x{P}_y{P}_z\right)=\left(A_{\mathrm{px}}+C_x{A}_{\mathrm{py}}+C_y{A}_{\mathrm{pz}}+C_z\right)\\ =\left(A_{\mathrm{px}}+C_x{A}_{\mathrm{fy}}+C_y{A}_{\mathrm{pz}}+C_z\right)\end{array}& \left(23\right)\\ \begin{array}{c}\left(V_x{V}_y{V}_z\right)=\left(A_{\mathrm{vx}}+C_x{A}_{\mathrm{vy}}+C_y{A}_{\mathrm{vz}}+C_z\right)\\ =\left(A_{\mathrm{vx}}+C_x{A}_{\mathrm{fy}}+C_y{A}_{\mathrm{vz}}+C_z\right)\end{array}& \left(24\right)\end{array}$

(57) The midpoint F between the maximum point P and the minimum point V is deformed as in the following Expression (25) from Expression (21) and Expressions (23) and (24).

(58) $\begin{array}{cc}\begin{array}{c}\stackrel{.fwdarw.}{F}=\frac{\stackrel{.fwdarw.}{P}+\stackrel{.fwdarw.}{V}}{2}=\left(\frac{A_{\mathrm{px}}+A_{\mathrm{vx}}}{2}+C_x{A}_{\mathrm{fy}}+C_y\frac{A_{\mathrm{pz}}+A_{\mathrm{vz}}}{2}+C_z\right)\\ =\left(0A_{\mathrm{fy}}+C_{y}0\right)\end{array}& \left(25\right)\end{array}$

(59) Therefore, C.sub.x which is the x-axis component of the original magnetic field C can be calculated by the following Expression (26), and C.sub.z which is the z-axis component of the original magnetic field C can be calculated by the following Expression (27).

(60) 0$\begin{array}{cc}{C}_x=-\frac{A_{\mathrm{px}}+A_{\mathrm{vx}}}{2}& \left(26\right)\\ C_z=-\frac{A_{\mathrm{pz}}+A_{\mathrm{vz}}}{2}& \left(27\right)\end{array}$

(61) In addition, when the x-axis component and the y-axis component are equal to each other in the maximum point P, C.sub.y which is the y-axis component of the original magnetic field C can be calculated by the following Expression (28) from the relation of P.sub.x=P.sub.y of Expression (19) and Expression (23).

(62) $\begin{array}{cc}{C}_y=A_{\mathrm{px}}+C_x-A_{\mathrm{fy}}=\frac{A_{\mathrm{px}}-A_{\mathrm{vx}}}{2}-A_{\mathrm{fy}}& \left(28\right)\end{array}$

(63) Meanwhile, the y-axis component A.sub.y of the artificial magnetic field A which is generated in the measurement region may be the fixed value A.sub.fy, and is arbitrary. When the fixed value A.sub.fy is zero, C.sub.y can be calculated according to the following Expression (29) at A.sub.fy=0.

(64) $\begin{array}{cc}{C}_y=\frac{A_{\mathrm{px}}-A_{\mathrm{vx}}}{2}& \left(29\right)\end{array}$

(65) In this manner, the y-axis component of the artificial magnetic field A is set to a fixed value, and a plane A.sub.xA.sub.y of A.sub.y=A.sub.fy is scanned by changing the x-axis component and the z-axis component, thereby allowing the original magnetic field C to be obtained from the artificial magnetic field (first magnetic field A.sub.p) in which the magnetization value (spin polarization degree M.sub.x) becomes maximum and the artificial magnetic field (second magnetic field A.sub.z) in which the magnetization value (spin polarization degree M.sub.x) becomes minimum. As described above, since the original magnetic field C is an external magnetic field entering the magnetic shield 8, the controlling of the magnetic field generator 7 so as to offset the original magnetic field C obtained in this manner during the magnetic measurement of a measuring object such as magnetocardiogram makes it possible to reduce the adverse influence of the external magnetic field and then measure the magnetic field of the measuring object with low noise.

(66) Flow of Processes

(67) FIG. 10 is a flow diagram illustrating a procedure of a magnetic field forming process of Embodiment 1. The magnetic field measurement device 1 performs the magnetic field forming process shown in FIG. 10 before a subject to be tested is loaded into the magnetic shield 8 and the measurement of a biomagnetic field is performed.

(68) As shown in FIG. 10, in the magnetic field forming process, first, the magnetic field calculation unit 53 calculates the original magnetic field C present in the measurement region (step S1). As a specific processing procedure, first, in a state where the relation of A.sub.y=A.sub.fy (A.sub.fy=0 in a preferred example) is established, the artificial magnetic field control unit 531 causes the magnetic field generator 7 to sequentially generate a plurality of artificial magnetic fields A by scanning each value of A.sub.x and A.sub.z within a predetermined range, and calculates the original magnetic field C through Expressions (26) to (28) using the first magnetic field A.sub.p when the square difference W.sub. satisfies the maximum value condition and the second magnetic field A.sub.v when the square difference satisfies the minimum value condition. In this case, the original magnetic field C becomes an external magnetic field entering the magnetic shield 8.

(69) Subsequently, the correction magnetic field setting unit 55 subtracts the original magnetic field C calculated in step S1 from a target magnetic field T=(T.sub.x, T.sub.y, T.sub.z) and sets a correction magnetic field TC (step S3). The correction magnetic field setting unit 55 then generates the set correction magnetic field TC in the magnetic field generator 7. Then, the original magnetic field C is offset by Expression (22) and the target magnetic field T is formed in the measurement region (step S5). In the present embodiment, as an example, the correction magnetic field TC=C is set as the target magnetic field T=(0, 0, 0), and a zero magnetic field is formed in the measurement region. In this manner, in a state where the target magnetic field such as the zero magnetic field is formed in the measurement region (that is, in a state where the correction magnetic field TC is generated in the measurement region), a minute magnetic field (magnetocardiogram or magnetoencephalo) which is emitted from a measuring object (for example, person) disposed in the measurement region is measured. Since the target magnetic field is formed in the measurement region, it is possible to measure a minute magnetic field with a high degree of accuracy.

(70) As described above, according to the present embodiment, it is possible to calculate the original magnetic field C present in the measurement region. In addition, the zero magnetic field is formed by offsetting the original magnetic field C in the measurement region, and the subject to be tested is load thereon into the magnetic shield 8, thereby allowing the measurement of a biomagnetic field to be performed. According to this, it is possible to measure the biomagnetic field with good sensitivity and with a high degree of accuracy.

(71) Embodiment 2

(72) In the aforementioned Embodiment 1, in a state where the y-axis component A.sub.y is set to the fixed value A.sub.fy (for example, zero), a plurality of artificial magnetic fields A are generated by changing the x-axis component A.sub.x and the z-axis component A.sub.z. On the other hand, in the present embodiment, in a state where the x-axis component A.sub.x is set to the fixed value A.sub.fx (for example, zero), a plurality of artificial magnetic fields A may be generated by changing the y-axis component A.sub.y and the z-axis component A.sub.z. The original magnetic field C=(C.sub.x, C.sub.y, C.sub.z) may be calculated by the following Expressions (30) to (32) using the first magnetic field A.sub.p=(A.sub.fx, A.sub.py, A.sub.pz) when the square difference W.sub. satisfies the maximum value condition and the second magnetic field A.sub.v=(A.sub.fx, A.sub.vy, A.sub.vz) when the square difference W.sub. satisfies the minimum value condition in which a minimum value is set.

(73) $\begin{array}{cc}{C}_x=-\frac{A_{\mathrm{py}}+A_{\mathrm{vy}}}{2}-A_{\mathrm{fx}}& \left(30\right)\\ C_y=\frac{A_{\mathrm{py}}-A_{\mathrm{vy}}}{2}& \left(31\right)\\ C_z=-\frac{A_{\mathrm{pz}}+A_{\mathrm{vz}}}{2}& \left(32\right)\end{array}$

(74) A method to obtain Expression (30) to Expression (32) is the same as that in Embodiment 1.

(75) In the magnetic field measurement method (how to calculate the original magnetic field C) of the present embodiment, the same operational effect as that in the above Embodiment 1 is also obtained.

(76) Embodiment 3

(77) Next, a magnetic field measurement device of Embodiment 3 and a magnetic field measurement method using the device will be described with reference to FIGS. 11 to 15.

(78) Entire Configuration

(79) First, the magnetic field measurement device of the present embodiment will be described with reference to FIG. 11.

(80) FIG. 11 is a diagram illustrating an entire configuration example of a magnetic field measurement device of Embodiment 3.

(81) As shown in FIG. 11, a magnetic field measurement device 10 of the present embodiment has substantially the same configuration as that of the magnetic field measurement device 1 of the above Embodiment 1. The magnetic field calculation unit 53 of the calculation control unit 5 includes an amplitude detection unit 533 in addition to the artificial magnetic field control unit 531. Therefore, the same components as those of the magnetic field measurement device 1 of Embodiment 1 are denoted by the same reference numerals and signs, and thus the detailed description thereof will not be given. Hereinafter, components having different configurations will be described.

(82) The artificial magnetic field control unit 531 controls the magnetic field generator 7, and sequentially generates, in the measurement region, artificial magnetic fields of a plurality of combinations obtained by changing an artificial magnetic field component in each axial direction of the x-axis, the y-axis, and the z-axis, that is, artificial magnetic fields obtained by periodically changing the z-axis component by superimposing a predetermined alternating-current component on the artificial magnetic field component (referred to as the z-axis component) in the z-axis direction. The amplitude detection unit 533 detects the ratio of the temporal change of the square difference W.sub. to the temporal change of the z-axis component by extracting an amplitude and a phase from the temporal change of the square difference W.sub. as an example of a value corresponding to the magnetization value calculated by the signal processing unit 51. The amplitude detection unit 533 can be configured using a lock-in amplifier or the like.

(83) How to calculate the original magnetic field C in the calculation control unit 5 of the present embodiment, that is, a magnetic field measurement method is different from that in Embodiment 1 or Embodiment 2.

(84) The magnetic field measurement method of the present embodiment is to obtain the original magnetic field C with a focus on a differential value M.sub.x/B.sub.z obtained by differentiating, by B.sub.z, the spin polarization degree M.sub.x as a magnetization value which is the x-axis component (component in the first direction) of the magnetization vector of the medium (gaseous atoms) in the gas cell 3.

(85) Specifically, the differential value M.sub.x/B.sub.z obtained by differentiating the spin polarization degree M.sub.x of Expression (14) by B.sub.z is represented by the following Expression (33).

(86) $\begin{array}{cc}\frac{M_x}{B_z}=\left(\frac{c}{a}\right)\frac{a\left(a^2+B_x^2+B_y^2-B_z^2\right)-2B_x{B}_y{B}_z}{{\left(a^2+B_x^2+B_y^2+B_z^2\right)}^2}& \left(33\right)\end{array}$

(87) The differential value M.sub.x/B.sub.z represented by Expression (33) shows a change in output (spin polarization degree M.sub.x) with respect to a change in the amount of detection (z component B.sub.z of a magnetic field), and thus means detection sensitivity. That is, in a condition in which the differential value M.sub.x/B.sub.z becomes maximum, the measurement sensitivity of the magnetic field measurement device 10 becomes maximum. Even when B.sub.x or B.sub.y has a certain value, the denominator of Expression (33) increases by the fourth power of these values. On the other hand, the numerator of Expression (33) changes at the square of B.sub.x or B.sub.y, and thus the differential value M.sub.x/B.sub.z becomes maximum when B.sub.x=B.sub.y=0 with respect to B.sub.x and B.sub.y. The differential value M.sub.x/B.sub.z in this case is represented by the following Expression (34).

(88) $\begin{array}{cc}\frac{M_x}{B_z}=\frac{c\left(a^2-B_z^2\right)}{{\left(a^2+B_z^2\right)}^2}& \left(34\right)\end{array}$

(89) FIG. 12 is a diagram illustrating a three-dimensional distribution relating to B.sub.x and B.sub.z of the differential value M.sub.x/B.sub.z, and FIG. 13 is a diagram illustrating a two-dimensional distribution relating to B.sub.z of the differential value M.sub.x/B.sub.z when B.sub.x=0.

(90) When each value of B.sub.x and B.sub.z is changed at B.sub.y=0[nT], a result as shown in FIG. 12 is obtained. Alternatively, when each value of B.sub.z is changed at B.sub.x=B.sub.y=0 [nT], a result as shown in FIG. 13 is obtained. As shown in FIGS. 12 and 13, one maximum value appears in the distribution of the differential value M.sub.x/B.sub.z and the differential value M.sub.x/B.sub.z becomes maximum when B.sub.z=0 in a case of B.sub.x=B.sub.y=0. Therefore, when the measurement region is a zero magnetic field (each value of B.sub.x, B.sub.y, and B.sub.z is zero), a change in the spin polarization degree M.sub.x (that is, square difference W.sub.) with respect to a change in magnetic field along the z-axis direction (direction of travel of the probe light) becomes maximum, and sensitivity becomes maximum. In other words, in order to minimize the original magnetic field such as an external magnetic field in the measurement region, the artificial magnetic field may be adjusted so that the differential value M.sub.x/B.sub.z becomes maximum.

(91) Here, the differential value M.sub.x/B.sub.z (ratio of the temporal change of the spin polarization degree M.sub.x to the temporal change of B.sub.z) can be replaced by the ratio of the temporal change of the square difference W.sub. to the temporal change of B.sub.z. Consequently, in the magnetic field forming process of the present embodiment, artificial magnetic fields of a plurality of combinations obtained by changing the x-axis component, the y-axis component, and the z-axis component, that is, artificial magnetic fields obtained by periodically changing the z-axis component are sequentially generated in the measurement region, and an artificial magnetic field in which the differential value M.sub.x/B.sub.z which is a ratio of the temporal change of the spin polarization degree M.sub.x to the temporal change of the B.sub.z becomes maximum is searched for. In this case, it is preferable that the angular frequency of the z-axis component is set to a value equal to or less than a cut-off angular frequency .sub.C. The cut-off angular frequency .sub.C is a sum of the relaxation rate .sub.0 of the spin polarization degree M.sub.x and the optical pumping rate .sub.p (.sub.C=.sub.0+.sub.p), and is approximately 100 Hz in the present embodiment. That is, it is preferable that the z-axis component of the artificial magnetic field satisfies a relation of <.sub.C=.sub.0+.sub.p, and the angular frequency is equal to or less than 100 Hz in the present embodiment. When the angular frequency of the z-axis component is set to a value equal to or less than the cut-off angular frequency .sub.C, dM.sub.x/dt can be regarded as approximately zero, and the approximation of the left-hand side of Expression (11) to be zero is justified. That is, when the relation of <.sub.C=.sub.0+.sub.p is satisfied, the original magnetic field C present in the measurement region can be accurately measured. However, since the magnetic field measurement device 10 acts like a primary low-pass filter, a gain and a phase in the vicinity of the cut-off angle frequency .sub.C gently decreases with an increase in angular frequency . For this reason, the angular frequency of a periodic function which is superimposed in reality may be made larger by approximately 10 percent than the cut-off angular frequency .sub.C.

(92) As a specific measurement method, the x-axis component A.sub.x and the y-axis component A.sub.y of the artificial magnetic field are set to a fixed magnetic field, and a magnetic field represented by a periodic function is applied to the z-axis component, to thereby measure a ratio of the temporal change of the spin polarization degree M.sub.x to the temporal change of B.sub.z. This is measured at various levels, and an artificial magnetic field in which the differential value M.sub.x/B.sub.z becomes maximum is specified. For example, as a first measurement, A.sub.x=0, A.sub.y=0, and A.sub.z are set to a vibration magnetic field around 0 (as an example, A.sub.z=sin t), and a first differential value M.sub.x/B.sub.z is measured. Next, as a second measurement, A.sub.x=0, A.sub.y=0, and A.sub.z are set to a vibration magnetic field around 1 (as an example, A.sub.z=1+sin t), and a second differential value M.sub.x/B.sub.z is measured. Next, as a third measurement, A.sub.x=0, A.sub.y=0, and A.sub.z are set to a vibration magnetic field around 1 (as an example, A.sub.z=1+sin t), and a third differential value M.sub.x/B.sub.z is measured. Next, as a fourth measurement, A.sub.x=0, A.sub.y=1, and A.sub.z are set to a vibration magnetic field around 0 (as an example, A.sub.z=sin t), and a fourth differential value M.sub.x/B.sub.z is measured. Next, as a fifth measurement, A.sub.x=0, A.sub.y=1, and A.sub.z are set to a vibration magnetic field around 1 (as an example, A.sub.z=1+sin t), and a fifth differential value M.sub.x/B.sub.z is measured. Next, as a sixth measurement, A.sub.x=0, A.sub.y=1, and A.sub.z are set to a vibration magnetic field around 1 (as an example, A.sub.z=1+sin t), and a sixth differential value M.sub.x/B.sub.z is measured. Next, as a seventh measurement, A.sub.x=0, A.sub.y=1, and A.sub.z are set to a vibration magnetic field around 0 (as an example, A.sub.z=sin t), and a seventh differential value M.sub.x/B.sub.z is measured. Next, as an eighth measurement, A.sub.x=0, A.sub.y=1, and A.sub.z are set to a vibration magnetic field around 1 (as an example, A.sub.z=1+sin t), and an eighth differential value M.sub.x/B.sub.z is measured. Next, as a ninth measurement, A.sub.x=0, A.sub.y=1, and A.sub.z are set to a vibration magnetic field around 1 (as an example, A.sub.z=1+sin t), and a ninth differential value M.sub.x/B.sub.z is measured. In this manner, a lot of measurements are repeated centering on A.sub.y=0, an artificial magnetic field in which the differential value M.sub.x/B.sub.z becomes maximum is specified from a lot of differential values M.sub.x/B.sub.z obtained in this manner. In other words, the differential value M.sub.x/B.sub.z is detected by extracting an amplitude from the period change of the square difference W.sub. calculated by the signal processing unit 51, on the basis of the detection results of the polarimeter 4, for each of a plurality of artificial magnetic fields of each combination. The original magnetic field C is calculated using an artificial magnetic field when a condition for external value in which the differential value M.sub.x/B.sub.z is set to a maximum value is satisfied. A process of obtaining the maximum value of the differential value M.sub.x/B.sub.z while changing the artificial magnetic field can be realized using publicly known optimization processing. As described above, in order to search for an artificial magnetic field in which the differential value M.sub.x/B.sub.z becomes maximum, making a search for the y-axis component A.sub.y of the artificial magnetic field from the vicinity of zero is preferable from the aspect of efficiency.

(93) FIG. 14 is a diagram illustrating a relationship between a coordinate system of a magnetic field B=(B.sub.x, B.sub.y, B.sub.z) of the measurement region and a coordinate system of an artificial magnetic field A=(A.sub.x, A.sub.y, A.sub.z) generated by the magnetic field generator 7, in the present embodiment. As shown in FIG. 14, each component of the magnetic field B=(B.sub.x, B.sub.y, B.sub.z) of the measurement region can be calculated by performing vector addition of each component of the artificial magnetic field A=(A.sub.x, A.sub.y, A.sub.z) to the original magnetic field C=(C.sub.x, C.sub.y, C.sub.z). A calculation expression is represented by the following Expression (35).
{right arrow over (B)}={right arrow over (A)}+{right arrow over (C)}=(B.sub.xB.sub.yB.sub.z)=(A.sub.x+C.sub.xA.sub.y+C.sub.yA.sub.z+C.sub.z)(35)
As described above, the magnetic field of the measurement region when the condition for external value in which the differential value M.sub.x/B.sub.z is set to a maximum value is satisfied, that is, when the extracted amplitude becomes maximum is a zero magnetic field (B=0), and thus the original magnetic field C (C.sub.x, C.sub.y, C.sub.z) is represented by the following Expressions (36) to (38), using an artificial magnetic field A.sub.h (A.sub.hx, A.sub.hy, A.sub.hz) when the condition for external value is satisfied, from the relation of Expression (35).
C.sub.x=A.sub.hx(36)
C.sub.y=A.sub.hy(37)
C.sub.z=A.sub.hz(38)
Flow of Processes

(94) FIG. 15 is a flow diagram illustrating a procedure of a magnetic field forming process in the present embodiment. The magnetic field measurement device 10 performs the magnetic field forming process shown in FIG. 15 before a subject to be tested is loaded into the magnetic shield 8 and the measurement of a biomagnetic field is performed.

(95) As shown in FIG. 15, in the magnetic field forming process, first, the magnetic field calculation unit 53 calculates the original magnetic field C present in the measurement region (step S11). As a specific processing procedure, first, the artificial magnetic field control unit 531 sequentially generates the above-mentioned artificial magnetic fields of a plurality of combinations in the measurement region through the magnetic field generator 7, and the amplitude detection unit 533 extracts the amplitude of the square difference W.sub. which is obtained as measurement results every time. The magnetic field calculation unit 53 then calculates the original magnetic field C by Expressions (36) to (38) using the artificial magnetic field A.sub.h when the amplitude of the square difference W.sub. becomes maximum.

(96) Subsequently, the correction magnetic field setting unit 55 sets the correction magnetic field T-C by subtracting the original magnetic field C calculated in step S11 from the target magnetic field T (T.sub.x, T.sub.y, T.sub.z) (step S13). The correction magnetic field setting unit 55 then generates the set correction magnetic field TC in the magnetic field generator 7, and thus forms the target magnetic field T in the measurement region by offsetting the original magnetic field C (step S15). In the present embodiment, as an example, the correction magnetic field TC is set as the target magnetic field T=(0, 0, 0), and a zero magnetic field is formed in the measurement region.

(97) As described above, according to the present embodiment, it is possible to calculate the original magnetic field C present in the measurement region. In addition, the zero magnetic field is formed by offsetting the original magnetic field C in the measurement region, and the subject to be tested is load thereon into the magnetic shield 8, thereby allowing the measurement of a biomagnetic field to be performed. According to this, it is possible to measure the biomagnetic field with good sensitivity and with a high degree of accuracy.

(98) A basic technical idea common to the magnetic field measurement methods of Embodiment 1 to the above Embodiment 3 has a technical feature in which there is provided a magnetic field measurement method of measuring a magnetic field of a measurement region in a magnetic field measurement device including a light source unit that emits linearly polarized light in which a first direction (x-axis direction), a second direction (y-axis direction) and a third direction (z-axis direction) are orthogonal to each other, a gas cell having a medium filled therein, which is disposed in the measurement region, which changes optical characteristics of the linearly polarized light in accordance with a magnetic field, and which is irradiated along the third direction (z-axis direction) with the linearly polarized light of which a vibration direction of an electric field is the second direction (y-axis direction), an optical detector that detects the optical characteristics, and a magnetic field generator that applies an artificial magnetic field to the measurement region, the method including generating a plurality of artificial magnetic fields, obtained by changing an artificial magnetic field component in the third direction (z-axis direction), in the magnetic field generator, calculating a magnetization value which is a component of a magnetization vector of the medium in the first direction or a value corresponding to the magnetization value on the basis of detection results of the optical detector, and calculating an original magnetic field C present in the measurement region, using an artificial magnetic field when the magnetization value or the value corresponding to the magnetization value satisfies a specified condition.

MODIFICATION EXAMPLE

(99) In addition, in the aforementioned embodiment, the gas cell 3 having alkali metal atoms sealed in a gaseous state is used as a medium for generating polarization in a magnetic moment and rotating the plane of polarization of transmitted light in accordance with the intensity of a magnetic field, but a medium other than the alkali metal atoms may be used. For example, a solid-state element such as a diamond providing a lattice defect due to nitrogen may be used as a medium. In addition, the magnetic field measurement method and the magnetic field measurement device according to various embodiments of the invention can also be applied to an atomic oscillator using a small gas cell of a millimeter size, in addition to a magnetic sensor.

(100) In addition, the embodiment in which the target magnetic field is set to a zero magnetic field has been described, but the target magnetic field can be set to any magnetic field other than the zero magnetic field.