Detection of zero-field resonance

Abstract

The present invention relates to a method of detecting a zero-field resonance (ZFR), comprising the steps of i) directing a pump light beam in a first direction x through a vapor cell comprising gaseous atoms, ii) applying a magnetic field component to the vapor cell parallel or antiparallel to the first direction x with different magnitudes: formula (I), iii) detecting light of the pump light beam directed through the vapor cell for the different magnitudes: formula (I) of the magnetic field component applied parallel or antiparallel to the first direction x, iv) determining the maximum optical transmission Tmax and the minimum optical transmission Tmin based on the detected light, v) determining the difference between the maximum optical transmission and the minimum optical transmission D = Tmax - Tmin, vi) determining the magnitude of the x-component of an ambient magnetic field: formula (II) by determining the particular magnitude: formula (III) of the x-component of the applied magnetic field for which the minimum optical transmission Tmin is determined, vii) repeating at least steps i) to v) while applying magnetic field components to the vapor cell in second and third directions y and z perpendicular to the first direction x with different magnitudes: formula (IV) and formula (V) and viii) determining the magnitudes: formula (VI) and formula (VII) of the ambient magnetic field components in the second and third directions y and z by determining the particular magnitudes of the y- and z-components of the applied magnetic field: formula (VIII) and formula (IX) for which the difference between the maximum optical transmission and the minimum optical transmission D is minimized.

Description

[0001] Detection of Zero-Field Resonance

[0002] Field of Invention

[0003] The present invention relates to a method and an apparatus for detecting zero-field resonance and, in particular, magnetometers employing the zero-field resonance effect.

[0004] Background of the invention

[0005] Zero-field resonance (ZFR) is a phenomenon in which an atomic medium is largely transparent to laser light at zero or near-zero magnetic field and is largely opaque to laser light at other values of the magnetic field. The transparency range might be, for example, -100 nT to +100 nT in any of the three vector components of the magnetic field. For comparison, the field of the earth is about 50 pT or 50,000 times larger than the transparency range.

[0006] Employing the zero-field resonance effect has a wide and growing range of applications including the detection of biomagnetic signals, neuroscience tools, geophysical exploration and prospecting, navigation and space applications as well as microwave communication, in general.

[0007] Particularly, the zero-field resonance effect can be employed by highly sensitive magnetometers that measure the direction and magnitude of an external magnetic field through the induced changes in the atomic spin polarization of an ensemble of atoms. In particular, optically pumped magnetometers (OPMs) operate based on the zero-field resonance effect. An optically pumped magnetometer is a device that uses lasers and an atomic medium to obtain an optical signal that indicates the magnetic field of interest. Typically, an OPM operating based on ZFR incorporates magnetic field coils arranged around the atomic medium and uses active feedback (biasing by means of the magnetic coils) to maintain the magnetic field at the sensed location within the transparency range. Thus, in the active feedback application the OPM must generate a field that cancels the ambient field or a residual ambient magnetic field still present after some magnetic field shielding which is not known beforehand but should be determined by the OPM. Therefore, the magnetic field applied by the magnetic coils is to be scanned over a range of values in the attempt of finding the correct strengths of the individual magnetic field components to compensate for the ambient field (thus “finding the zero” of the field). This search is complicated, particularly, since it has to be performed in all three spatial dimensions.

[0008] US 2016 / 0223627 A1 proposes a solution for the search problem based on the application of a bias magnetic field in the direction of the laser pump beam to simultaneously increase both the height and width of the ZFR. Subsequently, the ZFR is detected and a magnetic field is scanned / swept in directions perpendicular to that of the pump light beam (transversal directions of the magnetic field) and the transversal magnetic field components generated by one or more magnetic coils are adjusted to minimize the width and maximize the height of the zero-field resonance (transmission). However, this method requires precise adjustment of the applied magnetic field components relative to the ambient magnetic field and the ZFR width. It also involves detecting small changes in relatively large measurement quantities. The procedure, therefore, is rather time-consuming.

[0009] It is, therefore, an object of the present invention to provide a technique for employing the zero-field resonance effect with high accuracy and at high operation speed, particularly, in the context of magnetometer applications.

[0010] Description of the Invention

[0011] The present invention addresses the above-mentioned object by providing a method of detecting a zero-field resonance (ZFR), comprising the step of i) directing a pump light beam in a first direction x through a vapor cell comprising gaseous atoms. The gaseous atoms are or comprise one of alkali atoms, rubidium atoms, cesium atoms, potassium atoms, sodium atoms and helium atoms. The vapor cell may, additionally, contain a buffer gas provided to slow down the rate at which the atoms collide with the inner walls of the vapor cell such that randomization of the spins of the atoms caused by such collisions can be suppressed. The method further comprises the steps of ii) applying (by means of a magnetic field coil) a (longitudinal) magnetic field component to the vapor cell parallel or antiparallel to the first direction x with different magnitudes B*pp(i = 1 , N, N being an integer indicating the number of different magnitudes), iii) detecting light of the pump light beam directed through the vapor cell for the different magnitudes B*ppof the magnetic field component applied parallel or antiparallel to the first direction x (i.e., sweeping / scanning the applied longitudinal magnetic field component) and iv) determining (for the different N magnitudes of the applied longitudinal magnetic field) the maximum optical transmission Tmax and the minimum optical transmission Tmin based on the detected light. Then, in step v) the difference between the maximum optical transmission and the minimum optical transmission D = Tmax - Tmin (i.e., the depth of the optical transmission dip) is determined.

[0012] The method further comprises vi) determining the magnitude of the x-component of an ambient magnetic field B®xtby determining the particular magnitude B®pplof the x-compo- nent of the applied magnetic field for which the minimum optical transmission Tmin is determined. Here and in the following the ambient magnetic field is an external naturally occurring magnetic field or a residual external naturally occurring magnetic field that remains after shielding the vapor cell by some magnetic shielding. This magnitude of the x-component of an ambient magnetic field B®xtmay be considered to be equal to the particular one of the magnitudes of the applied longitudinal magnetic field component for which the minimum transmission Tmin is determined.

[0013] Tmax, Tmin and, therefore, D depend on (transversal) magnetic field components applied to the vapor cell in second and third directions y and z perpendicular to the first direction x. Therefore, the method, furthermore, comprises vii) repeating at least steps i) to v) (with possibly varying N) while applying (transversal) magnetic field components to the vapor cell in second and third directions y and z perpendicular to the first direction x with different magnitudes B*ppand B*pp, and viii) determining the magnitudes B®xtand B®xtof the ambient magnetic field components in the second and third directions y and z by determining the particular magnitudes of the y- and z-components of the applied magnetic field B®ppand B®ppfor which the difference between the maximum optical transmission and the minimum optical transmission D is minimized. The magnitudes B®xtand B®xtof the ambient magnetic field components are considered to be equal to or are determined based on the particular magnitudes ByPPand B*ppof the applied transversal magnetic field components.

[0014] It is noted that when carrying out step ii) the first time before the repetition cycle vii) the y- and z-components of the applied magnetic field may be chosen to be zero, for example. Further, when repeating step ii) during the repetition cycle vii) the number of different magnitudes N of the applied longitudinal magnetic field component may be varied with respect to step ii) carried out for the first time before the repetition cycle vii). According to a particular example, N=1 during the repetition cycle vii) and the only applied longitudinal magnetic field component has the particular magnitude B^ppdetermined in step vi).

[0015] It is noted that step vi) of determining the magnitude of the x-component of an ambient magnetic field B®xtby determining the particular magnitude B^ppof the x-component of the applied magnetic field for which the minimum optical transmission Tmin is determined may also be repeated through the cycle of step vii) and, for example, the x-component of an ambient magnetic field B®xtmay be determined as an average of the particular magnitudes B^ppof the x-component of the applied magnetic field for which the minimum optical transmission Tmin is determined for different magnitudes of the y- and z-components of the applied magnetic field B®xtand B®xt.

[0016] This method allows for a fast and reliable determination of the ambient magnetic field and, thus, the ZFR. The sweeping / scanning of the applied longitudinal magnetic field component while observing the transmission allowing for determining optical transmission dips with depths D is essential for the present invention. The inventive method is based on a minimization of D such that all nonzero values of D are significant and, consequently, optimization can be performed speedily without high precision measurements that rely on small differences between relatively large finite values rather than differences with respect to zero (see also detailed description below) as needed in the prior art. Scan and overall processing speed can be increased as compared to the art. Particularly, according to the inventive method there is no need for a computationally expensive fitting of the width of the ZFR or any time expensive modulation / demodulation processing. Further, accuracy of the determined ambient magnetic field is increased by the inventive method as compared to the art. According to an implementation, the magnitudes and BappOf the applied transversal magnetic field components in the second and third directions y and z are varied in step vii) based on at least one of a raster scan, spiral search starting from a starting magnitude and moving outward in a spiral pattern, Lissajous search with magnitudes oscillating sinusoidally at different frequencies, a Rose curve search by oscillating the transversal applied magnetic field magnitudes J + (Bapp)2sinusoidally while linearly scanning the angle arctan(Bzpp / Bapp) and a sample of representative points. These kinds of searches allow for a speedy and reliable adjustment of the applied transversal magnetic field components to the actual transversal magnetic field components of the ambient magnetic field in order to determine the ZFR.

[0017] According to an implementation, in step viii) the magnitudes B®xtand Bzxtof the transversal ambient magnetic field components are determined by weighted averages of the magnitudes Bappand Bappof the applied transversal magnetic field components wherein the weights are functions of the determined optical transmissions. This weighting procedure allows for a fast and reliable determination of the magnitudes B®xtand Bzxtof the transversal ambient magnetic field components. For example, the magnitudes B®xtand Bzxtare determined by wherein Tmaxis the maximum value of the determined transmission Tt, x is a power-law exponent, v is the threshold fraction and i = 1 , .., n, n being an integer indicating the number of different magnitudes. According to an implementation, the method further comprises generating the pump light beam by a laser device. Alternatively, the pump light beam may be generated by a suitable vapor lamp.

[0018] The above-mentioned object is also addressed by providing an apparatus for detecting a zero-field resonance, comprising a vapor cell comprising gaseous atoms (for example, being or comprising one of alkali atoms, rubidium atoms, cesium atoms, potassium atoms, sodium atoms and helium atoms), a light source configured to direct a pump light beam in a first direction x through the vapor cell, magnetic coils arranged and configured to apply a magnetic field to the (interior of the) vapor cell, a detector device configured to detect light of the pump light beam directed through the vapor cell, and a processing and controlling means. The processing and controlling means is configured to perform

[0019] I) directing a pump light beam in a first direction x through a vapor cell comprising gaseous atoms;

[0020] II) applying a (longitudinal) magnetic field component to the vapor cell parallel or antiparallel to the first direction x with different magnitudes B^pp(i = 1 , .., N, N being an integer indicating the number of different magnitudes);

[0021] III) detecting light of the pump light beam directed through the vapor cell for the different magnitudes B^ppof the magnetic field component applied parallel or antiparallel to the first direction x;

[0022] IV) determining the maximum optical transmission Tmax and the minimum optical transmission Tmin based on the detected light;

[0023] V) determining the difference between the maximum optical transmission and the minimum optical transmission D = Tmax- Tmin;

[0024] VI) determining the magnitude of the x-component of an ambient magnetic field B®xtby determining the particular magnitude Bppof the x-component of the applied magnetic field for which the minimum optical transmission Tmin is determined; VII) repeating at least I) to V) while applying (transversal) magnetic field components to the vapor cell in second and third directions y and z perpendicular to the first direction x with different magnitudes ByPPand B^pp; and

[0025] VIII) determining the magnitudes B®xtand B®xtof the ambient magnetic field components in the second and third directions y and z by determining the particular magnitudes of the y- and z-components of the applied magnetic field ByPPand Bzppfor which the difference between the maximum optical transmission and the minimum optical transmission D is minimized.

[0026] According to an implementation, the magnitudes ByPPand B^ppof the magnetic field components in the second and third directions y and z are varied in step VII) performed by the processing and controlling means based on at least one of a raster scan, spiral search starting from a starting magnitude and moving outward in a spiral pattern, Lissajous search with magnitudes oscillating sinusoidally at different frequencies and Rose curve search by oscillating the transversal applied magnetic field magnitudes J + (Bapp)2sinusoidally while linearly scanning the angle arctan(Bapp / Bapp).

[0027] According to an implementation, in step VIII) performed by the processing and controlling means the magnitudes B®xtand B®xtof the ambient transversal magnetic field components are determined by weighted averages of the magnitudes Bappand Bappof the applied transversal magnetic field components wherein the weights are functions of the determined optical transmissions. The magnitudes B®xtand B®xtof the ambient transversal magnetic field components may be determined by wherein Tmaxis the maximum value of the determined transmission Tt, x is the power-law exponent, v is the threshold fraction and i = 1 , n, n being an integer indicating the number of different magnitudes.

[0028] According to an implementation, the light source is a laser device configured for generating the pump (laser) light beam.

[0029] According to an implementation, the apparatus further comprises a magnetic shielding configured for substantially shielding an external magnetic field from the vapor cell. The components of the ambient magnetic field determined as described above, in this case, are the components of the remaining external magnetic field remaining after magnetic shielding.

[0030] The apparatus and any implementation thereof provide respective similar advantages as discussed above with respect to the provided method and implementations thereof. The provided method and any implementation thereof may be performed by means of the provided apparatus and any implementation thereof.

[0031] Furthermore, it is provided a zero-field resonance optically pumped magnetometer comprising the apparatus according to any of the above-described examples. Very accurate measurement of 3D magnetic fields can be provided by such a zero-field resonance optically pumped magnetometer.

[0032] Additional features and advantages of the present invention will be described with reference to the drawings. In the description, reference is made to the accompanying figures that are meant to illustrate preferred embodiments of the invention. It is understood that such embodiments do not represent the full scope of the invention.

[0033] Brief description of the drawings

[0034] Figure 1 schematically illustrates components and operation of an apparatus for detecting a zero-field resonance of an exemplary embodiment of the invention wherein a laser beam optically pumps spins of an atomic species and the optical transmission of the laser beam is detected while sweeping an applied magnetic field component orientated parallel or antiparallel to the direction of propagation of the laser beam.

[0035] Figures 2a and b represent a flow chart illustrating a method of detecting a zero-field resonance in accordance with an exemplary embodiment of the invention.

[0036] Figure 3 illustrates a comparison of an approach of the prior art based on minimizing the width and maximize the height of the zero-field resonance (a) and an approach based on minimizing the transmission depth according to an embodiment (b).

[0037] Figure 4 illustrates the effect of small variations of applied transversal magnetic field components on depths of optical transmission dips.

[0038] Description of specific embodiments of the invention

[0039] The present invention provides a method of detecting a zero-field resonance (ZFR) and an apparatus for detecting a zero-field resonance based on scanning / sweeping an applied magnetic field component orientated parallel or anti-parallel to the direction of propagation of a laser beam pumping spins of an atomic species contained in a vapor cell.

[0040] Figure 1 schematically illustrates components of an apparatus 10 for detecting a zero-field resonance according to an embodiment. The apparatus 10 comprises a laser device 1. The laser device 1 may be a Vertical Cavity Surface Emitting Laser (VCSEL) with an emission wavelength of about 795 nm (D1 line of87Rb), for example.

[0041] The laser device 1 is used to generate a pump light (laser) beam in order to optically pump (the spins of) atoms contained in a vapor cell 2. For example, the pump light beam may generated by the laser device 1 at a power level of the order of pW. In general, the wavelength of the light emitted by the laser device 1 has to be adapted to the wavelength of a D1 or D2 line atomic transition of the atoms. The atomic medium contained in the vapor cell 2 is largely transparent to laser light at zero or near-zero magnetic field and is largely opaque to laser light at other values of the magnetic field (Zero-field resonance (ZFR) effect). For example, the atomic medium may comprise or consists of an alkali metal vapor. In the example shown in Figure 1 the pump light beam enters the vapor cell 2 in the x- direction (beam propagation direction). The same pump light beam is used for monitoring the optical transmission through the vapor cell 2. For example, the pump light beam originating from the laser device 1 is collimated and guided to the vapor cell over a free space path of about 10 cm. A mirror, half-wave plate, a polarizing beam splitter cube and a quarter wave plate (not shown in Figure 1) may be used to circularly polarize the pump light beam and direct it through the vapor cell 2.

[0042] According to an embodiment, a buffer gas, for example, helium, neon or nitrogen, is provided in the vapor cell 2 in order to slow down the rate at which the atoms (for example, atoms of an alkali metal vapor) collide with the inner walls of the vapor cell 2 such that randomization of the spins of the alkali metal atoms caused by such collisions can be suppressed. The buffer gas pressure in the vapor cell 2 may be, for example, in a range of about 1 to 5 bars at room temperature but the method is equally applicable for any buffer gas density.

[0043] To increase the number density of the atoms of the active atomic medium the vapor cell 2 may be heated, for example, up to a temperature over 100°C, for example, 150 °C to 200 °C. For example, the vapor cell 2 may be located in an oven or a separate heating device may be provided for heating the active atomic medium. For example, a density of the active atomic medium in the range of 107cm-3to 1015cm-3may be considered appropriate. According to a particular embodiment, the vapor cell 2 is filled with isotopically enriched Ru- bidium-87 and a nitrogen buffer gas obtained through the process of Rubidium-azide (RbNs) decomposition caused by irradiation of UV light. Besides rubidium other suitable active atomic media may consist of or comprise cesium, potassium or sodium.

[0044] The vapor cell 2 may be sealed through anodic bonding of borosilicate glass to silicon. Special coatings on the inner walls of the vapor cell 2, such as octadecyltrichlorosilane (OTS) or paraffin, may be provided to reduce spin randomization from wall collisions of the active atomic medium. All components of the housing of the vapor cell 2 should be chosen to be either non-magnetic or to have a very low residual magnetization. The apparatus 10, furthermore, comprises a set of (electrically activated biasing) biplanar magnetic field generating coils 3 with optimized current traces for generating a homogeneous magnetic field (for example, with a strength in the range of a few 10 nT to 100 nT or a few 100 nT or of the order of pT) at the position of the vapor cell 2 (to cancel out / compensate for some (residual) ambient magnetic field) and minimum stray fields around the module. The coils 3 may be configured to provide a magnetic field with a strength in order of an ambient magnetic field (or, after shielding, a residual ambient magnetic field). The vapor cell 2 is positioned in between the biplanar coils 3 and the whole module may be placed inside a four-layer magnetic shield housing giving raise to some residual ambient magnetic field being present at the location of the active atomic medium. The shield housing may be comprised of 3 layers of mu-metal (non-oriented nickel-iron soft ferromagnetic alloy) and one layer of ferrite, but any material that can suitably shield the ambient magnetic noise may equally be used acceptable. Further, only a partially shielding housing or no shield housing at all may be present according to different embodiments and applications.

[0045] The apparatus 10, furthermore, comprises a function generator (processing means) 4 configured for linearly sweeping / scanning (for example, from some predetermined negative value to some predetermined positive value, see sketch above the vapor cell 2 in Figure 1) the applied magnetic field component Bppthat is orientated parallel or antiparallel to the beam propagation direction x. The function generator 4 may control magnitudes and directions of electric currents applied to one of the coils 3 in order to perform the sweeping / scanning of the magnetic field component Bx. Further, the function generator 4 may be configured to alter applied y- and z- components of the magnetic field ByPPand B*ppat the location of the vapor cell 2 (again, by controlling the application of currents to the corresponding ones of the coils 3).

[0046] The apparatus 10, furthermore, comprises a detector device 5 configured for detecting transmission of the pump light beam through the (active atomic medium of) the vapor cell 2. The detector device 5 may be or comprises an Si photodetector and may be supplemented by some (for example, photodiode) signal amplifier. Particularly, the detector device 5 is configured for detecting a ZFR signal defined by the optical transmission T(BX, By, Bz) as a function of the magnetic field at the vapor cell 2 B = Bext- Bappgiven by the difference of the ambient magnetic vector field and the applied magnetic vector field, and, more particularly, for detecting optical transmission dips (see sketch above the detector device 5 in Figure 1 and description below). It is noted that the transmission T as a function of the applied magnetic field T(B^PP, ByPP, B^pp) has the same shape as the transmission T as a function of the magnetic field at the vapor cell 2 T(BX, By, Bz) with a center at Bext.

[0047] Figures 2a and 2b represent a flow chart illustrating a method 20 of detecting a zero-field resonance (ZFR) in accordance with an exemplary embodiment of the invention. For example, the method 20 can be carried out by means of the apparatus 10 illustrated in Figure 1 . Whereas the steps of the method 20 are shown in a particular order for illustration purposes individual steps may be carried out concurrently or even in a different order where considered appropriate.

[0048] In step S21 of the method 20 illustrated in Figure 2a a pump light beam is directed in a first direction x through a vapor cell comprising gaseous atoms. The pump light beam may be generated by the laser device 1 shown in Figure 1. The vapor cell may be the vapor cell 2 shown in Figure 1. The gaseous atoms may comprise one of alkali atoms, rubidium atoms, cesium atoms, potassium atoms, sodium atoms and helium atoms. The vapor cell may further contain a buffer gas.

[0049] A magnetic field component is applied S22 to the vapor cell parallel or antiparallel to the first direction x with different magnitudes Bpp, i = 1 , .., N, N being an integer indicating the number of different magnitudes. This step of sweeping / scanning the longitudinal component of the applied magnetic field (i.e., the component parallel / antiparallel to the direction of propagation of the pump light beam) with N different magnitudes is essential for the present invention.

[0050] The method 20 illustrated in Figure 2a, further, comprises detecting S23 light of the pump light beam directed through the vapor cell for the different magnitudes Bppof the magnetic field component applied parallel or antiparallel to the first direction x and determining S24 the maximum optical transmission Tmax and the minimum optical transmission Tmin based on the detected light. The depth D of the dip of the transmission T is determined by determining S25 the difference between the maximum optical transmission and the minimum optical transmission determined in step S24: D = Tmax - Tmin. It is noted that Tmax and Tmin and, thus, the depth D of the dip are also functions of the transversal components of the applied magnetic field ByPPand B*pp. The magnitude of the x-component of an ambient magnetic field B®xtis determined in step S26 by determining the particular magnitude B®pplof the x-component of the applied magnetic field for which the minimum optical transmission Tmin is obtained. This magnitude B®xtof the x-component of the ambient magnetic field may be considered to be equal to the particular magnitude Bppof the N magnitudes of the magnetic field component applied parallel or antiparallel to the first direction x in step S22 for which the minimum optical transmission Tmin is obtained.

[0051] As already said Tmax and Tmin and, thus, the depth D of the dip are also functions of the transversal components of the applied magnetic field ByPPand B*pp. Therefore, at least steps S21 to S25 are repeated S27 while applying (transversal) magnetic field components ByPPand Bppto the vapor cell in second and third directions y and z perpendicular to the first direction x with different magnitudes ByPPand B^ppand the magnitudes B®xtand B®xtof the ambient magnetic field components in the second and third directions y and z are determined S28 by determining the y- and z-components of the applied magnetic field for which the difference between the maximum optical transmission and the minimum optical transmission D is minimized. These magnitudes By*1and B *1of the ambient magnetic field components in the second and third directions y and z may be considered to be equal to or may be determined based on the particular ones ByPPand Bppof the N magnitudes of the applied transversal magnetic field components for which the difference between the maximum optical transmission and the minimum optical transmission D is minimized. It is noted that when carrying out step S22 the first time before the repetition cycle S27 the y- and z-components of the applied magnetic field may be chosen to be zero. Further, when repeating step S22 during the repetition cycle S27 the number of different magnitudes N of the applied longitudinal magnetic field component may be varied with respect to the first performance of step S22 before the repetition cycle S27. According to a particular example, N=1 during the repetition cycle S27 and the only applied longitudinal magnetic field component has the particular magnitude B^ppdetermined in step S26.

[0052] The magnitudes ByPPand Bppof the magnetic field components in the second and third directions y and z may be varied in step S27 based on at least one of a a. raster scan b. spiral search starting from a starting magnitude and moving outward in a spiral pattern c. Lissajous search with magnitudes oscillating sinusoidally at different frequencies; and d. Rose curve search by oscillating the transversal applied magnetic field

[0053] > magnitudes J + (Bapp)2sinusoidally while linearly scanning the angle arctan(Bapp / Bapp).

[0054] The magnitudes B®xtand B®xtof the ambient transversal magnetic field components may be determined is step S28 by weighted averages of the magnitudes ByPPand Bappof the applied transversal magnetic field components wherein the weights are functions of the determined optical transmissions. For example, the magnitudes B®xtand B®xtof the ambient transversal magnetic field components are determined by

[0055] Bext y W< next > — yni=i B ozaPiP

[0056] Wi, with r11. > vvT1max otherwise wherein Tmaxis the maximum value of the determined transmission Tt, x is a power-law exponent, v is the threshold fraction and i = 1 , n, n being an integer indicating the number of different magnitudes.

[0057] It is noted that step S26 of determining the magnitude of the x-component of an ambient magnetic field B®xtby determining the particular magnitude Bxppof the x-component of the applied magnetic field for which the minimum optical transmission Tmin is determined may also be repeated through the cycle of step S27 and, for example, the x-component of an ambient magnetic field B®xtmay be determined as an average of the particular magnitudes Bxppof the x-component of the applied magnetic field for which the minimum optical transmission Tmin is determined for different magnitudes of the y- and z-components of the applied magnetic field B®ppand B®pp.

[0058] As described above, according to the invention an applied longitudinal magnetic field component is scanned and the depth D of a transmission dip (the difference between the maximum optical transmission and the minimum optical transmission) is determined S25 and the depth D is minimized S28 to adjust ByPPand B*ppto the ambient field to cancel it out at the location of the vapor cell. This procedure is different from the art exemplified by US 2016 / 0223627 A1 as it is illustrated in Figure 3. US 2016 / 0223627 A1 proposes application of a bias magnetic field in the direction of the laser pump beam to simultaneously increase both the height and width of the ZFR and, after detection of the ZFR, a magnetic field is scanned / swept in directions perpendicular to that of the pump light beam (transversal directions of the magnetic field) and these transversal magnetic field components are adjusted to minimize the width and maximize the height of the zero-field resonance (see left column (a) of Figure 3 showing optical transmission versus a transversal magnetic field component). Contrary, according to the invention optical transmission as a function of a scanned longitudinal magnetic field component is swept through zero (see right column (b) of Figure 3 showing optical transmission versus the longitudinal magnetic field component).

[0059] As compared to the art, the inventive approach provides the advantage that optimization has to be with respect to a substantially non-zero value of the depth D of a transmission dip which is possible with a relatively low calculation precision needed as compared to an optimization with respect to relatively large finite values of the height and width of the ZFR wherein the detection of small changes of the relatively large measuring quantities have to be considered. Particularly, significantly shorter averaging times for obtaining reasonably useful measurement results are needed according to the inventive method and, thus, the scan speed and overall processing speed can be increased. Further, accuracy of the determined ambient magnetic field is increased by the inventive method.

[0060] As it is illustrated in Figure 4 even a small variation of the current applied to a coil providing applied transversal magnetic field components has a relatively large effect on the depth of the dip (D) of the optical transmission. In the example shown in Figure 4 the optical transmission measured in V is plotted against the longitudinal scanning time in s (and, thus, variation in the applied longitudinal magnetic field component) for coil currents for the application of the transversal magnetic field components of ly= 0 mA and Iz = 0.000 mA (curve a), ly= 0 mA and Iz = -0.005 mA (curve b), and ly= 0 mA and Iz = -0.100 mA (curve c), respectively. The minimum of D is obtained for ly= 0 mA and Iz = 0.000 mA (curve a). Even a small variation of magnitude of the applied current of 0.005 mA has a significant effect on D close to the zero field (see curve b). Thus, adjustment to the ambient field can be achieved quickly and accurately.

[0061] The inventive method and apparatus, for example, the embodiments illustrated in Figures 1 and 2a and 2b can be advantageously used in a great variety of applications including optically pumped magnetometers (OPMs), for example, in the context of magnetoencephalography (a medical imaging technique that provides functional information about the brain and is used for pre-surgical characterization of epilepsy and for diagnosis of mild traumatic brain injury) or in the context of the characterization of batteries, e.g., for transportation, wherein ZFR-OPMs are used to non-invasively detect and localize leakage currents and diagnose battery health. Other applications include magnetic navigation using the Earth’s field as a reference and magnetic communication for underwater and underground vehicles and installations, geotechnical and minerals exploration techniques, and clinical and pre- clinical magnetic resonance investigations wherein ZFR-OPMs are used together with 'hyperpolarized’ nuclear spins in molecules used, for example, as tracer compounds for metabolism studies in patients.

[0062] All previously discussed embodiments are not intended as limitations but serve as examples illustrating features and advantages of the invention. It is to be understood that some or all of the above-described features can also be combined in different ways.

[0063] Acknowledgement

[0064] This work has been partially funded by the European Union’s Horizon Europe research and innovation programme under grant agreement No. 101099379, by the Secretariat of Digital Policies of the Government of Catalonia - G.A. GOV / 51 / 2022, and has received the support of a fellowship from ”la Caixa” Foundation (ID 100010434), under the fellowship code LCF / BQ / DI24 / 12070013.

Claims

Claims1 . Method (20) of detecting a zero-field resonance, comprising the steps of i) directing (S21) a pump light beam in a first direction x through a vapor cell comprising gaseous atoms; ii) applying (S22) a magnetic field component to the vapor cell parallel or antiparallel to the first direction x with different magnitudes B^pp; iii) detecting (S23) light of the pump light beam directed through the vapor cell for the different magnitudes B®ppof the magnetic field component applied parallel or antiparallel to the first direction x; iv) determining (S24) the maximum optical transmission Tmax and the minimum optical transmission Tmin based on the detected light; v) determining (S25) the difference between the maximum optical transmission and the minimum optical transmission D = Tmax- Tmin; vi) determining (S26) the magnitude of the x-component of an ambient magnetic field B®xtby determining the particular magnitude Bppof the x-component of the applied magnetic field for which the minimum optical transmission Tmin is determined; vii) repeating (S27) at least steps i) to v) while applying magnetic field components to the vapor cell in second and third directions y and z perpendicular to the first direction x with different magnitudes ByPPand B^pp; and viii) determining (S28) the magnitudes B®xtand B®xtof the ambient magnetic field components in the second and third directions y and z by determining the particular magnitudes of the y- and z-components of the applied magnetic field Rappand Bppfor which the difference between the maximum optical transmission and the minimum optical transmission D is minimized.

2. The method (20) according to claim 1, wherein the magnitudes B^ and B^pOf the magnetic field components in the second and third directions y and z are varied in step vii) based on at least one of a a. raster scan b. spiral search starting from a starting magnitude and moving outward in a spiral pattern c. Lissajous search with magnitudes oscillating sinusoidally at different frequencies; and d. Rose curve search.

3. The method (20) according to claim 1 or 2, wherein in step viii) the magnitudes B®xtand B®xtof the ambient transversal magnetic field components are determined (S28) by weighted averages of the magnitudesand B^PP of the applied transversal magnetic field components wherein the weights are functions of the determined optical transmissions.

4. The method (20) according to claim 3, wherein the magnitudes B®xtand B®xtof the ambient transversal magnetic field components are determined byB^ = X?=1B^Wi, withwherein Tmaxis the maximum value of the determined transmission Tt, x is a power-law exponent, v is the threshold fraction and i = 1, .., n, n being an integer indicating the number of different magnitudes.

5. The method (20) according to one of the preceding claims, further comprising generating the pump light beam by a laser.

6. The method (20) according to one of the preceding claims, wherein the gaseous atoms are or comprise one of alkali atoms, rubidium atoms, cesium atoms, potassium atoms, sodium atoms and helium atoms.

7. An apparatus (10) for detecting a zero-field resonance, comprising a vapor cell (2) comprising gaseous atoms; a light source (1) configured to direct a pump light beam in a first direction x through the vapor cell (2); magnetic coils (3) arranged and configured to apply a magnetic field to the vapor cell (2); a detector device (5) configured to detect light of the pump light beam directed through the vapor cell; and a processing and controlling means configured to performI) directing a pump light beam in a first direction x through a vapor cell comprising gaseous atoms;II) applying a magnetic field component to the vapor cell parallel or antiparallel to the first direction x with different magnitudes B^pp;III) detecting light of the pump light beam directed through the vapor cell for the different magnitudes B®ppof the magnetic field component applied parallel or antiparallel to the first direction x;IV) determining the maximum optical transmission Tmax and the minimum optical transmission Tmin based on the detected light;V) determining the difference between the maximum optical transmission and the minimum optical transmission D = Tmax- Tmin;VI) determining the magnitude of the x-component of an ambient magnetic field B®xtby determining the particular magnitude Bppof the x-component of the applied magnetic field for which the minimum optical transmission Tmin is determined;VII) repeating at least I) to V) while applying magnetic field components to the vapor cell in second and third directions y and z perpendicular to the first direction x with different magnitudes ByPPand B^pp; andVIII) determining the magnitudes B®xtand B®xtof the ambient magnetic field components in the second and third directions y and z by determining the particular magnitudes of the y- and z-components of the applied magnetic field ByPPand B®ppfor which the difference between the maximum optical transmission and the minimum optical transmission D is minimized.

8. The apparatus (10) according to claim 7, wherein the magnitudes ByPPand B^ppof the magnetic field components in the second and third directions y and z are varied in step VII) performed by the processing and controlling means based on at least one of a a. raster scan b. spiral search starting from a starting magnitude and moving outward in a spiral pattern c. Lissajous search with magnitudes oscillating sinusoidally at different frequencies; and d. Rose curve search.

9. The apparatus (10) according to claim 7 or 8, wherein in step VIII) performed by the processing and controlling means the magnitudes B®xtand B®xtof the ambient transversal magnetic field components are determined by weighted averages of themagnitudes B^ andof the applied transversal magnetic field components wherein the weights are functions of the determined optical transmissions.

10. The apparatus (10) according to claim 9, wherein the magnitudes B®xtand B®xtof the ambient transversal magnetic field components are determined bywherein Tmaxis the maximum value of the determined transmission Tt, x is the power-law exponent, v is the threshold fraction and i = 1, n, n being an integer indicating the number of different magnitudes.

11. The apparatus (10) according any of the claims 7 to 10, wherein the light source (1) is a laser device (1) configured for generating the pump light beam.

12. The apparatus (10) according to any of the claims 7 to 11, wherein the gaseous atoms are or comprise one of alkali atoms, rubidium atoms, cesium atoms, potassium atoms, sodium atoms and helium atoms.

13. The apparatus (10) according to any of the claims 7 to 12, further comprising a magnetic shielding configured for substantially shielding an external magnetic field from the vapor cell.

14. A zero-field resonance optically pumped magnetometer comprising the apparatus (10) according to any of the claims 7 to 13.

Images