Performance analysis of an optically pumped magnetometer in Earth's magnetic field

TL;DR

This paper experimentally analyzes how the orientation of optically pumped magnetometers affects their performance in Earth's magnetic field, focusing on a mode promising femtotesla resolution and demonstrating robustness against heading errors.

Contribution

It provides an experimental evaluation of the LSD-Mz regime's key parameters and shows that using balanced laser beams reduces heading errors and maintains sensitivity across orientations.

Findings

Heading error is minimal with balanced laser beams.

Field sensitivity remains stable over a wide orientation range.

The regime achieves femtotesla-level resolution in Earth's magnetic field.

Abstract

We experimentally investigate the influence of the orientation of optically pumped magnetometers in Earth's magnetic field. We focus our analysis to an operational mode that promises femtotesla field resolu-tions at such field strengths. For this so-called light-shift dispersed Mz(LSD-Mz) regime, we focus on the key parameters defining its performance. That are the reconstructed Larmor frequency, the transfer function between output signal and magnetic field amplitude as well as the shot noise limited field resolution. We demonstrate that due to the use of two well balanced laser beams for optical pumping with different helicities the heading error as well as the field sensitivity of a detector both are only weakly influenced by the heading in a large orientation angle range.

Equations

$$I=I(\nu )+I_{\_}dc=I_{\_}0\frac{\Delta {\nu }^2}{4(\nu -{\nu }_{\_}0{)}^2+\Delta {\nu }^2}+I_{\_}dc,$$

$$B_{\_}sn=\frac{\sqrt{2e\left(I_{\_}{dc}^{{\sigma }_{\_}+}({\nu }_{\_}L)+I_{\_}{dc}^{{\sigma }_{\_}-}({\nu }_{\_}L)\right)}}{\gamma s},$$

$$H=\frac{h{\nu }_{\_}0}{2}{\sigma }_{\_}z+h\Omega \cos 2\pi \nu t\left[\cos \alpha {\sigma }_{\_}x+\sin \alpha {\sigma }_{\_}z\right].$$

$$H^{\prime }=\frac{h{\nu }_{\_}0}{2}{\sigma }_{\_}z+h\Omega \cos 2\pi \nu t\sin \alpha \left[e^{ig\sin 2\pi \nu t}{\sigma }_{\_}++e^{-ig\sin 2\pi \nu t}{\sigma }_{\_}-\right].$$

$$⟨{\sigma }_{\_}z⟩=-\frac{{\Gamma }_{\_}r{\Gamma }_{\_}{\phi }^{\prime }}{{\Gamma }_{\_}r{\Gamma }_{\_}{\phi }^{\prime }+{\Omega }^2{\cos }^2\alpha }.$$

$$\begin{array}{rl}{\Gamma }_{\_}r& =\frac{{\gamma }_{\_}r}{2}+\frac{{\Omega }_{\_}p\left|\cos \alpha \right|}{2}\left(1-\cos \alpha \right),\\ {\Gamma }_{\_}e& =\frac{{\gamma }_{\_}r}{2}+\frac{{\Omega }_{\_}p\left|\cos \alpha \right|}{2}\left(1+\cos \alpha \right).\end{array}$$

$$I(\omega )=I_{\_}0⟨{\sigma }_{\_}z⟩=I_{\_}0\frac{({\Gamma }_{\_}e-{\Gamma }_{\_}r){\Gamma }_{\_}{\phi }^{\prime }}{({\Gamma }_{\_}r+{\Gamma }_{\_}e){\Gamma }_{\_}{\phi }^{\prime }+c^2(\alpha ){\Omega }^2},$$

$$I_{\_}dc=I_{\_}0\frac{{\Gamma }_{\_}e-{\Gamma }_{\_}r}{{\Gamma }_{\_}r+{\Gamma }_{\_}e}=I_{\_}0\frac{{\cos }^2\alpha }{\left|\cos \alpha \right|+p_{\_}1},$$

$$I(\nu )=I_{\_}dc\left(1-\frac{c^2(\alpha ){\Omega }^2}{\left({\Gamma }_{\_}e+{\Gamma }_{\_}r\right){\Gamma }_{\_}{\phi }^{\prime }+c^2(\alpha ){\Omega }^2}\right).$$

$$i_{\_}0=I(\nu ={\nu }_{\_}0)/I_{\_}dc=\frac{c^2(\alpha )}{c^2(\alpha )+p_{\_}2\left(p_{\_}1+\cos \alpha \right)},$$

$$\Delta \nu =2{\Gamma }_{\_}\phi \sqrt{1+\frac{c^2(\alpha )}{p_{\_}2(p_{\_}1+\left|\cos \alpha \right|)}+p_{\_}3}.$$

$$\Delta {\nu }_{\_}{B}_{\_}1=\frac{c^2{\Gamma }_{\_}\phi }{p_{\_}2\left(p_{\_}1+\left|\cos \alpha \right|\right))}=\frac{c^2{\Omega }^2}{{\Gamma }_{\_}\phi \left(\gamma +{\Omega }_{\_}p\left|\cos \alpha \right|\right)}$$

$$\vec{D}\cdot \vec{E}=\frac{E_{\_}0}{2\sqrt{2}}{e}^{i\left(\vec{k}\vec{r}-2\pi {\nu }_{\_}Lt\right)}\left(\left[\cos \alpha \pm 1\right]D_{\_}++\left[\cos \alpha \mp 1\right]D_{\_}-+2\sin \alpha D_{\_}z\right).$$

$${\Omega }_{\_}L=\frac{1}{h}\sqrt{\frac{I_{\_}L}{4cn{ϵ}_{\_}0}}⟨J=1/2\left|\left|\vec{D}\right|\right|J^{\prime }=1/2⟩.$$

$$\Delta {\nu }_{\_}laser=\frac{{\Omega }_{\_}{L}^2}{2{\Gamma }_{\_}\phi {\Gamma }_{\_}opt}=3.5,$$

$$\frac{s}{2}=\left|\frac{dI}{d\nu }\right|=\frac{8I_{\_}0\Delta {\nu }^2(\nu -{\nu }_{\_}0)}{{\left(4{\left(\nu -{\nu }_{\_}0\right)}^2+\Delta {\nu }^2\right)}^2}$$

$$B_{\_}sn=\sqrt{\frac{e}{I_{\_}dc}}\frac{1}{8\gamma i_{\_}0}\frac{{\left(4{\nu }_{\_}{LS}^2+\Delta {\nu }^2\right)}^2}{\Delta {\nu }^2{\nu }_{\_}LS},$$

Full text

I Introduction

Optically pumped magnetometers (OPMs) are, in principle, scalar-type quantum sensors for magnetic fields based on the Zeeman effect. That is the shift of energy levels due to the interaction of atoms with an external magnetic field Cohen-Tannoudji et al. (1977). Usually alkali vapors in paraffin-coated glass cells are used as sensing element. Because the energy shift by the Zeeman interaction is based on the scalar product of the measured external magnetic field $\vec{B}_{\_}0$ and the magnetic moment of the atom, such a magnetometer measures the absolute value of the field. This fact makes them interesting for the realization of absolute field sensors Leger et al. (2009); Nabighian et al. (2010); Korth et al. (2016).

On the other hand, the alkali vapor is usually polarized by optical pumping with a circular laser beam to create large signal amplitudes. With the pump beam an additional direction dependence is introduced leading besides dead zones also to unwanted effects usually summarized under the label ”heading error”Budker and Kimball (2013). These effects can be understood from the change in the atom-light coupling that is in dipole approximation given by the scalar product of the laser’s electric field and the atoms electric dipole moment Oelsner et al. (2019). Suppression of the heading error is intensively studied Jensen et al. (2009); Bao et al. (2018); Hu et al. (2018) and it has important consequences in the application of OPMs Colombo et al. (2016); Weis et al. (2017).

The most sensitive magnetic field sensors are based on superconducting quantum interference devices Clarke and Braginski (2004); Schmelz and Stolz (2017). They have been optimized to allow sub-femtotesla field gradient resolution even in Earth’s magnetic field Stolz et al. (2001) and today, besides others, are used for geomagnetic and archeological explorations Chwala et al. (2001); Foley (2004); Stolz (2015). Still the requirements due to cryogenic liquids is demanding and for low temperature superconducting sensors the costs of liquid helium are high. In this context, the newly introduced operational modes which enable shot-noise-limited field resolutions of OPMs on the femtotesla scale in Earth’s magnetic field strengths are promising alternatives. Namely, those are the light narrowing (LN) Scholtes et al. (2011); Guo et al. (2019); Fu et al. (2019) and the LSD-Mz Schultze et al. (2017) mode.

In this work, we analyze the performance of magnetometers based on the LSD-Mz mode and experimentally investigate their characteristics as a function of their orientation in an external magnetic field. The paper is organized as follows: We describe the experimental setup before we discuss a theoretical description of the measured values. Afterwards we present experimental data and relate it to our theoretical description. We demonstrate the influence of different rotational axis to the results. Finally, our model allows an estimate of the influence of the heading direction on the field resolution which we use to conclude on the usability of OPMs based on this mode in Earth’s magnetic field.

II Methods

II.a Experimental description

Because we are interested in the characteristics of a high resolution OPM as a function of its orientation relative to the Earth’s magnetic field, our experimental setup is based on a vapor cell design suitable for the LSD-Mz regime. Namely, the used micro-fabricated buffer gas cell contains two active volumes connected to the same reservoir. It is created by ultrasonic milling the desired structure into a 4 mm thick silicon wafer. After filling the reservoir region with droplets of diluted cesium azide and a drying step, the cell is closed by anodic bonding with Borofloat glass plates. Finally, decomposing the azide to cesium and nitrogen as buffer gas with excimer laser irradiation finalizes the fabrication. A detailed description of the cell fabrication can be found in Woetzel et al. (2011). A photograph of the used cell arrangement is shown in Fig. 1 a).

The cell is mounted at the center of a rotational table and accompanied by optic elements for the preparation and detection of the required circularly polarized laser beams, see Fig. 1 b). A schematic of the optical setup is shown in Fig. 1 c). The laser light is guided to the experiment by a polarization maintaining fiber. The beam preparation (BP) is done in a compact setup. Therein a lens collimates the beam to a diameter of about 4 mm and a linear polarizer (LP) is used to define the polarization direction. The beam is then split by the use of a polarizing beam splitter (PBS). A right angle deflecting prism (DP) and two quarter-lambda plates ($\lambda$/4) create two parallel beams with right and left handed circular polarization, respectively. An additional linear polarizer (LP) in the straight beam is used for balancing their intensities. The so prepared beams are passed through the active volume of the two cells, before they are slightly focused with lenses (L) onto two photodiodes (PD) for recording the transmitted intensities. Unwanted light on the diodes, as for example from the heating laser, is blocked with an absorbing bandpass filter F.

The pump laser (PL) with a wavelength of 895 nm is stabilized to the Doppler free absorption line for the $F=3\rightarrow F^{\prime}$ = 4 D${}_{\_}1$ transition of the Cs vapor of an additional paraffin coated glass cell. A 978 nm heat laser is guided to the setup by a fiber and an optical setup, containing collimating lens and deflecting prism (not show in Fig. 1). The sensor is heated to a temperature of roughly 100 degree Celsius.

The whole optical installation is mounted inside of a three layer mu-metal shielding (see Fig. 1b) and Ref. Schultze et al. (2010)). Also, a three axis Helmholtz-coil system is included allowing the application of arbitrary magnetic fields. In all discussed experiments we apply a static magnetic field $\vec{B}_{\_}0$ of about $50$ $\mu$T along the z-axis that corresponds to the rotational axis of the cylinder-shaped shields. The normal of the rotational table and thus the rotational axis is labeled by y. Two additional Helmhotz coil configurations are mounted around the cell allowing for the application of magnetic rf-fields ($\vec{B}_{\_}1$) perpendicular to the laser light for the detection of the magnetic resonance. We respectively denote the amplitude and the frequency of the $B_{\_}1$ field by $\Omega$ and $\nu$. Note, the $\vec{B}_{\_}1$ field created by the y-coil is perpendicular to both, laser light direction $\vec{k}$ and magnetic field $\vec{B}_{\_}0$ for each rotation angle $\alpha$. But using the second rf-coil leads to an angle modification between $\vec{B}_{\_}0$ and $\vec{B}_{\_}1$ from perpendicular to parallel configuration during rotation. Nevertheless, we label the latter as ”x-coil” according to its initial orientation. Our setup thus resembles the real situation of an OPM moved in Earth’s magnetic field.

The measurements are carried out as follows: At first the heading angle is adjusted from outside of the shielding by a cable pull to an angle $\alpha$ between $\vec{B}_{\_}0$ and the propagation direction of the laser light $\vec{k}$. The latter is given by $\sin\alpha\vec{e}_{\_}x+\cos\alpha\vec{e}_{\_}z$. Here, the $\vec{e}_{\_}i$ denote the unit vectors in i-direction. Then, consecutively the $B_{\_}1$ field is swept through the magnetic resonance using the $x$ as well as $y$ coil. Thus, for each angle $\alpha$ recordings are made with a magnetic rf-field applied in $y$ and $\cos\alpha\vec{e}_{\_}x+\sin\alpha\vec{e}_{\_}z$ direction. Because at an angle $\alpha=0$ both configurations are equivalent, we adjusted the current feed to the different $B_{\_}1$ coils to produce magnetic resonances with the same height and width at these angles to correct for the slightly different coil constants.

In the starting configuration ($\alpha=0$) the LSD-Mz mode requires two detuned circularly polarized laser beams both of which are oriented in parallel to the external magnetic field $B_{\_}0$ Schultze et al. (2017). When using buffer gas cells, the population is pumped to dark states at $m_{\_}F=\pm 4$ of the $F=4$ ground states, respectively for $\sigma_{\_}+$ and $\sigma_{\_}-$ light. By the magnetic rf-field the transition between the Zeeman-split levels can be driven, which is observable by an increase in laser absorbtion. The substraction of the signals for both helicities results in a steep linear measurement curve around the actual Larmor frequency $\gamma B_{\_}0$. Here $\gamma=3.5$ Hz/nT denotes the gyromagnetic ratio of cesium. This substraction also reduces the common noise present on the two signal, as for example intensity noise from the laser. For illustration, corresponding measurements of the magnetic resonances are presented in Fig. 2 for an angle of $\alpha=0$ using the y-coil.

By fitting this measurement curves of the photocurrent as a function of rf-frequency $I(\nu)$ with a Lorentzian function

[TABLE]

we extract the DC-photocurrents $I_{\_}{dc}$, the amplitudes $I_{\_}0$ and widths (FWHM) $\Delta\nu$ of the magnetic resonances, as well as the resonant frequencies $\nu_{\_}0$. The latter is given by the light shifted Larmor frequency $\nu_{\_}0=\nu_{\_}L\pm\nu_{\_}{LS}$, where $\nu_{\_}{LS}$ denotes the frequency shift due to the AC-Stark effect. Although the helicity of the two laser beams as well as their relative intensities were balanced to less than 1 % prior to recording the experimental data, a remarkable deviation in both, the dc-current and the resonance amplitude, is observed. Therefore, in all measurements including the one presented in Fig. 2 an additional electronic balancing through amplifying the $\sigma_{\_}+$ signal by an additional factor of $a=1.61$ is introduced. The reason for the observed deviation, that is dependent on the applied magnetic field requires further investigations.

From the steepness of the difference signal $s=d(I_{\_}{\sigma_{\_}+}-I_{\_}{\sigma_{\_}-})/d\nu$ at the Larmor frequency we calculate the theoretical shot-noise limited sensitivity as

[TABLE]

where the photocurrents $I_{\_}{dc}$ are taken at the Larmor frequency and $e$ is the elementary charge. The heading characteristics of the fitting parameters of Eq. 1 as well as of the sensitivity are in the following discussed.

II.b Theoretical description

Two theoretical concepts are required for the description of our experimental results. On the one hand, the measured resonance frequencies are strongly modified by the light shift due to the intense off-resonant pumping. This effect has been extensively studied Oelsner et al. (2019); Jensen et al. (2009); Hu et al. (2018); Mathur et al. (1968) and we use the method presented in Ref. Oelsner et al. (2019) for a description of our data.

On the other hand, the magnetic resonance itself can be described by Bloch equations. Here, we firstly analyze the Hamiltonian describing the magnetic resonance of a two-level system

[TABLE]

for the x-coil. Please note that the same Hamiltonian is valid for the y-coil keeping $\cos\alpha=1$ and $\sin\alpha=0$. In Hamiltonian (3), $\Omega=\gamma B_{\_}1$ is introduced as the amplitude of the $B_{\_}1$ field in frequency units, $\sigma_{\_}x$, $\sigma_{\_}y$, as well as $\sigma_{\_}z$ denote Pauli matrices, and $h$ the Planck constant. In order to proceed, we aim to discuss the system in a rotating frame that removes the time dependent diagonal term. That can be achieved by a unitary transformation $U_{\_}1=e^{iG(t)\sigma_{\_}z}$ where we choose $G(t)=\frac{\Omega}{\nu}\sin 2\pi\nu t\cos\alpha$. As constructed, $U_{\_}1$ commutes with $\sigma_{\_}z$ and the term $i\hbar\dot{U}_{\_}1U_{\_}1^{\dagger}$ removes the diagonal coupling Shevchenko et al. (2014). The transformation therefore results in

[TABLE]

Here, we used $\sigma_{\_}\pm=0.5\left(\sigma_{\_}x\pm i\sigma_{\_}y\right)$ and the abbreviation $g=\frac{\Omega}{\nu}\cos\alpha$. Although we could continue with the Jacobi-Anger expansion and finally find the stationary terms in a frame rotating with $\nu$, we note that the ratio $\Omega/\nu$ is small in our experimental realization. Thus, we can summarize the lower line to $\sigma_{\_}x$ and solve the Bloch equations in a frame rotating with $\nu$ around the z-axis in rotating wave approximation. The modification of the stationary result for the z-component of the magnetization compared to the usual result is limited to a reduction of the effective rf-amplitude with a factor $\cos\alpha$ as

[TABLE]

Above we introduced the rates $\Gamma_{\_}r$ and $\Gamma_{\_}\varphi$ as respected inverse relaxation and coherence times $T_{\_}1$ and $T_{\_}2$, as well as $\Gamma_{\_}\varphi^{\prime}=(\Gamma_{\_}\varphi^{2}+\delta^{2})/\Gamma_{\_}\varphi$. The detuning of the $B_{\_}1$ frequency from resonance is included by $\delta=\nu_{\_}0-\nu$.

Still, this modification for the x-coil alone fails to accurately explain our measured data. What is missing in Eq. (5) is the modification of population transfer by optical pumping, or in other words the change in the dc-photocurrent that depends on the population difference. To include the optical transition and the decay from the atoms excited state into the two-level model, we modify the relaxation and excitation dissipative dynamics, similar as in Refs. Dréau et al. (2011) and Shi et al. (2018). Namely, we set for the respective excitation $\Gamma_{\_}e$ and decay rate $\Gamma_{\_}r$

[TABLE]

Here we assumed a repopulation rate of the ground state levels $\gamma_{\_}r/2$ and an effective pumping rate between the two levels $\Omega_{\_}p\propto\Omega_{\_}L$. We note that the equations above are written for $\sigma_{\_}+$ polarization of the light. Still, they take the same form for $\sigma_{\_}-$ light. This other helicity only differs by a change of the signs inside the brackets. On the other hand, because the $m_{\_}F=-3$ level has higher energy compared to the $m_{\_}F=-4$ in two level approximation, the result for $\sigma_{\_}-$ would be equivalent to Eq. 6. The rates resulting from the optical pumping in Eq. 6 include contributions of the scalar product between the atoms’s dipole moment and the electric field ($\propto 1\pm\cos\alpha$) (compare to Ref. Oelsner et al. (2019)) as well as an additional factor of $\left|\cos\alpha\right|$. Also, we did not include the linear polarized component scaling with $\sin\alpha$ that, in principle, can also contribute to a population difference. This statement is especially true close to a perpendicular pumping orientation $\vec{B}_{\_}0\perp\vec{k}$. The rates defined in (6) enter to the observable photocurrent. It relates to the expectation value of $\sigma_{\_}z$ as

[TABLE]

Here, the function $c(\alpha)$ is constantly one for the x-coil and $\cos\alpha$ for the y-coil. As constructed, our model accurately describes the dc-photocurrent as shown in Section III. Namely, if $\delta\rightarrow\pm\infty$ or equivalently $\Omega\rightarrow 0$, we find

[TABLE]

with the dimensionless fitting parameter $p_{\_}1=\gamma_{\_}r/\Omega_{\_}p$. The above function changes with decreasing optical pumping from a $\cos$ to a $\cos^{2}$ angular dependence while at the same time the amplitude is decreased. This reduction of the dc-current summarizes the fact that for an effective pumping to the dark states the rate of equalization of population $\gamma_{\_}r$ should be significantly smaller than the optical pumping $\propto\Omega_{\_}p$.

To bring Eq. 7 into a similar form as Eq. 1 we separate the dc-current to find the Lorentzian-like resonance function as

[TABLE]

The resonance amplitude is accordingly the value of the second term in above equation at resonance, namely at $\nu=\nu_{\_}0$. It is

[TABLE]

where we introduced a second dimensionless fitting parameter $p_{\_}2=\Gamma_{\_}\varphi\Omega_{\_}p/\Omega^{2}$.

Finally, the full width at half maximum (FWHM) of the resonance curve of Eq. (7) is given by

[TABLE]

This relation describes the quadratic addition of the natural linewidth $\Gamma_{\_}\varphi$ and the power broadening by the $B_{\_}1$ field

[TABLE]

of the magnetic resonance line. To fit our experimental results, it is necessary to introduce an additional constant broadening term $p_{\_}3$ that we account to the broadening of the magnetic resonance line by the laser power as we will demonstrate below. We introduce this solely to the resonance width, since a saturation in the optical transition does not yield the same in the magnetic resonance and thus will not influence to the amplitude of the measured Lorentzian-shaped signal.

III Results

III.a Resonant frequency - Light shift

In a first step we analyze the angular dependence of the magnetic resonance center frequencies as plotted in Fig. 3. We added to the figure also the mean value of the measured center frequencies found for the two different circular polarizations.

In general, the magnetic resonances are strongly shifted by the ac-Stark shift due to the strong off-resonant pumping. The angle between the laser beam direction, represented by its $k$-vector, and the magnetic field $\vec{B}_{\_}0$ influences to the atom-light interaction by the transition dipole moment Oelsner et al. (2019)

[TABLE]

Therein $E_{\_}0$ is the amplitude of the electric field, $\vec{k}$ the k-vector of the laser, $\vec{r}$ the atoms position, $\nu_{\_}L$ the frequency of the laser beam, and the $D_{\_}i$ are components of the dipole operator. In words, above equation states that with modifying the angle $\alpha$ the weight of pumping to different excited states is strongly influenced. Namely, $D_{\_}+$, $D_{\_}-$, and $D_{\_}z$ components couple to excited states with increased, decreased, as well as unchanged magnetic quantum number $m_{\_}F^{\prime}=m_{\_}F+\left\{+1,-1,0\right\}$, respectively.

The strongest characteristic results from the vector light shift that can be interpreted as virtual magnetic field added in direction of the lights angular momentum Mathur et al. (1968). It is strongest close to angles of zero and $\pm\pi$ and reduces to zero close to perpendicular orientation. Also as expected, it changes its sign for a change in the laser’s helicity. A more detailed discussion of the light shift can be found in Ref. Oelsner et al. (2019). We use the findings therein to calculate the expected light-shifted transition frequencies between Zeeman levels with highest and lowest magnetic quantum numbers to their neighboring states. These correspond to the respective blue and orange solid lines plotted in Fig. 3 and can be identified with the transitions probed when pumping to the dark states $m_{\_}F=\pm 4$, where the plus and minus sign are to be used for $\sigma_{\_}+$ and $\sigma_{\_}-$ light, respectively. For the curves we used a magnetic field strength of $B_{\_}0=49.664$ $\mu$T, detunings of $\delta_{\_}{F^{\prime}=4}=-8$ GHz and $\delta_{\_}{F^{\prime}=3}=-9.2$ GHz from the respective optical transitions $F=4\rightarrow F^{\prime}=4$ and $F=4\rightarrow F^{\prime}=3$, a linewidth of the optical transitions $\Gamma_{\_}{opt}\approx 4$ GHz, and optical driving amplitudes in frequency units $\Omega_{\_}L=3.15$ MHz and 3.45 MHz for $\sigma_{\_}+$ and $\sigma_{\_}-$ polarized beam, respectively. The last values corresponds to the on-resonance Rabi frequencies Oelsner et al. (2013) introduced by the pumping beams.

We achieve a good correspondence between our model and the experimental results. It is best close to angles of zero and $\pm\pi$, where additionally to the strong light shift also a very good pumping to the dark states is achieved. Close to angles of $\pm\pi/2$, not only the observed light shift is reduced but also the population is distributed between several ground state levels. This effect reduces the magnetic resonance amplitude and, additionally, enables probing more ground state transitions both reducing the agreement between theory and experiment.

Nevertheless the model allows the reconstruction of the laser intensities $I_{\_}L$ from

[TABLE]

Here the vacuum permittivity $\epsilon_{\_}0$, speed of light $c$, reduced transition dipole moment $\left\langle\!J=1/2\!\left|\!\left|\vec{D}\right|\!\right|\!J^{\prime}=1/2\!\right\rangle$ Steck (2010), as well as the refractive index $n$ enter to the equation. We found the two intensities to be 64 and 76 W/m2. Since their ratio does not correspond to the scaling factor $a$, we expect some influence of slightly varying beam profiles or non-perfect circular polarization. If we assume the laser power distributed equally to a circular beam with a diameter of 4 mm we find a total value of about 1 mW which fits well to separate measurements.

III.b The dc-photocurrent

A second basic characteristic is found for the dc-photocurrent as introduced in Eq. 1. It is plotted as a function of the heading angle in Fig. 4.

As demonstrated in the figure, we observe a dependence that roughly follows a $\left|\cos\alpha\right|$ or $\left|\cos\alpha\right|^{2}$ function. Furthermore, the orientation of the $B_{\_}1$ coils has no influence to the photocurrent away from the magnetic resonance, as expected. Although we used an electronic balancing, we still observe a slightly higher dc-photocurrent for the $\sigma_{\_}-$ beams compared to the $\sigma_{\_}+$ helicity of about 2 % in maximum that is best demonstrated from the difference signal.

The explanation of the shape of the measurement curves is found in the interaction of the dipole moment $\vec{D}$ with the laser’s electric field $\vec{E}$ as given in Eq. (13). We included this effect by modified relaxation and excitation rates depending on the effective circularly polarized laser intensity (6). The measured photocurrent is increased when the pumping of the atoms to dark states is more efficient because they cannot absorb anymore light. A large dark state population is related to large differences in the ratio of light pumping to larger and smaller magnetic quantum numbers. They are respectively proportional to $D_{\_}+=D_{\_}x+iD_{\_}y$ and $D_{\_}-=D_{\_}x-iD_{\_}y$.

Our model (Eq. 8) predicts a change from a $\cos^{2}\alpha$ dependence to one proportional to $\left|\cos\alpha\right|$ when the effective pumping amplitude $\Omega_{\_}L$ is increased compared to the relaxation of the polarization given by a rate $\gamma_{\_}r$. That allows for a more accurate fit of the photocurrent as a function of the orientation angle $\alpha$ as demonstrated by the solid lines in Fig. 4. For this curve, we estimate the ratio of relaxation to optical pumping rate to be $p_{\_}1\approx 0.3$ and the amplitudes $I_{\_}{0}=0.49$ mA. Keeping in mind the electronic amplification of the $\sigma_{\_}+$ channel, the variable $p_{\_}1$ is multiplied by $1/1.6$ for this channel making it necessary also to adjust the amplitude $I_{\_}0$ when fitting the curve measured for this helicity to 0.56 mA. Additionally, in consistency with the experiment, our model predicts a smaller photocurrent with decreasing the optical pumping amplitude $\propto\Omega_{\_}p$ that usually is not included in Bloch equations.

III.c Magnetic resonance amplitude

In contrast to the very similar angular dependencies of the dc-photocurrent for the two different $B_{\_}1$ coils, a clear discrepancy is found in the normalized amplitude of the magnetic resonance signal $i_{\_}0=I_{\_}0/I_{\_}{dc}$, as demonstrated in Fig. 5.

This mentioned discrepancy is well explained by the additional modification of the effective $B_{\_}1$ field amplitude when modifying its direction compared to $\vec{B}_{\_}0$. We included this in our calculation by the factor $c$ that is constantly 1, when using the $y$-coil, and $\left|\cos\alpha\right|$, when the $x$-coil is used. With Eq. 5 we fit the normalized amplitudes presented in Fig. 5 and found very good agreement between experiment and theory with a factor $p_{\_}2=0.12$. Note, for $\sigma_{\_}+$ we additionally had to adjust the value of $p_{\_}2$ by a factor of $1.6$. The dependence we observed for the resonance amplitude driven by the x-coil is strongly influenced by the $\cos$ function given for the effective $B_{\_}1$-field amplitude $\Omega$. Thus the amplitude values is maximal at $\alpha=0$ and $\pi$ and tends to zero in vicinity of $\alpha=\pm\pi/2$, where effectively no $B_{\_}1$ field remains.

In contrast, when using the $y$-coil, our theory predicts an increase in resonance amplitude towards one for angles close to $\alpha=\pm\pi/2$. Except in close vicinity of these angles this general tendency is also found in our experiment. From our model, it is clear that this increase is connected to a reduction of the optical pumping to the dark state: A given strength of the $B_{\_}1$-field corresponds to a certain rate of population shifting back from the dark state to absorbing states. In other words, the $B_{\_}1$ field introduces a Rabi oscillation whose frequency depends on the field amplitude. Thus for larger powers, the shifting of population to the absorbing state is faster. If this rate is smaller than the possible population transfer that can be achieved by the pumping laser, the measured photocurrent in resonance does not reach down to zero. Therefore the amplitude $I_{\_}0/I_{\_}{dc}$ is smaller than one. When rotating towards $\pm\pi/2$ and using the $y$-coil, the effective $B_{\_}1$ field amplitude stays constant. Because at the same time the optical pumping to the dark state gets less effective, the photocurrent amplitude is increased towards one. That means that all the atoms that are optically pumped contribute to the resonance with the $B_{\_}1$ field. Still, our model does not accurately reproduce the extracted values in close vicinity of $\alpha=\pm\pi/2$, where the experimentally observed amplitudes drop to zero. This deviation probably results from not considering the linearly polarized pumping at these angles, that not only leads to optical alignment but also takes the role of the $B_{\_}1$ field in redistributing population.

III.d Resonance width

A fitting of the experimentally observed resonance widths with the same parameters $p_{\_}1$ and $p_{\_}2$ turned out to be unsuccessful. Thus we found it necessary to introduce an additional fitting parameter $p_{\_}3$ to Eq. 11. Theoretical curves with parameters $p_{\_}3=3.5$ and $\Gamma_{\_}\varphi=350$ Hz are plotted together with the experimental data points in Fig. 6.

Our adjusted model again fits nicely to the experimental results in the case of the $x$-coil as demonstrated by the solid lines’ correspondence to the data presented as diamonds in Fig. 6. We observe a reduction of the power broadening introduced by the rf-field. Therefore, we assume that close to angles of $\alpha=\pm\pi/2$ the minimal possible magnetic resonance width is achieved for this certain temperature and laser power.

In contrast, applying a constant effective magnetic rf-field by the $y$-coil, leads to a strong increase of the resonance width in the experiment. We account this to a large redistribution of population between Zeeman states resulting in an overlay of several ground state transitions and the linear optical pumping remaining at these angles. Since our model is restricted to two-levels and we neglected the linear pumping, we fail to catch the magnitude of this increase in magnetic resonance width. Still, we note that the qualitative behavior is accurately reproduced by the theoretical model.

As already mentioned above, we account the factor $p_{\_}3$ to a power broadening due to the strong laser power. To justify this assumption we can estimate the expected laser power broadening similar to Eq. 12 by (compare e.g. Ref. Loudon (2000))

[TABLE]

where we used the values for $\Omega_{\_}L$ of the $\sigma_{\_}+$ beam and $\Gamma_{\_}{opt}$ as noted above. Since this value is in agreement with the fitting parameter, we identify the strong power of the detuned laser as one important source of broadening of the magnetic resonance.

III.e Shot-noise-limited sensitivity

Finally, as a key parameter of the proposed magnetometer we extracted the shot-noise limited sensitivity of the LSD-Mz signal as explained in Sec. II.a. The respective experimental values as a function of the angle $\alpha$ for the two rf-coil-orientations are presented in Fig.7.

Our experiment shows a minimal shot-noise limited resolution for our experimental parameters of 20 fT/$\sqrt{\text{Hz}}$. This is already close to the optimal conditions Schultze et al. (2017). This optimal resolution is found in the case of parallel alignment of the laser direction to the magnetic field. This results from the strong dependence of the sensors performance on the pumping to the dark state and achievement of a reasonable light shift. Both of these requirements are connected to the atom-light interaction and thus are reduced towards perpendicular orientation. There, both the magnetic resonance signal as well as the light shift are reduced resulting in a loss of sensitivity.

We observe a strong dependence of the amplitude and width of the magnetic resonance signals on the angle that additionally is slightly different for the two used orientations of the $B_{\_}1$ coils. Still, because the width of the magnetic resonance as well as the light shift are both dominated by the laser intensity, the shot noise limited resolution is quite stable for angles of $\pm 20^{\circ}{}$ around optimal orientation.

As shown above, our model allows for a description of the orientation dependence of all the key parameters of the magnetic resonance. Thus it enables us to estimate the shot noise limited sensitivity as a function of $\alpha$ in the case of perfect balancing of the channels. To do so, we calculate the steepness of a single resonance curve as the derivative of (1)

[TABLE]

and evaluate it at $\nu=\nu_{\_}0+\nu_{\_}{LS}$. There the maximal steepness of the LSD-Mz signal is achieved and corresponds to twice the value of a single magnetic resonance $s(\nu_{\_}0+\nu_{\_}{LS})=2\left|\frac{dI}{d\nu}\right|(\nu_{\_}0+\nu_{\_}{LS})$. Substituting into the equation for the shot-noise-limited resolution (2) results in

[TABLE]

with the parameters $I_{\_}{dc}$, $i_{\_}0$, and $\Delta\nu$ defined by the respective equations (8), (10), and (11). The theoretical estimated sensitivity is added to Fig. 7. In the calculation we used the experimental parameters of the starting angle and for $\sigma_{\_}+$ light and assumed perfect channel balancing. In general a slightly better sensitivity is expected from our calculation but the qualitative shape fits very well to the experiment. Note, the two additional peaks around angles of $\pm\pi/2$ result from the overlay of the two resonance curves for the $\sigma_{\_}\pm$ beams meaning $\omega_{\_}{LS}=0$, as can be seen in the plotted resonance frequencies of Fig. 3. This overlay results in a difference signal constantly equal to zero and thus $s\rightarrow 0$. Still, due to the lack of contrast at perpendicular configuration, this effect is not visible in the experimental data. Finally, we note that although a perfect balancing has only little influence to the sensitivity, it makes the sensor more robust against heading error, concerning the reconstructed absolute magnetic field, see green data points of Fig. 3.

Discussion of results and conclusion

We experimentally analyzed the performance of an OPM operated in the LSD-Mz mode at Earth magnetic field strengths as a function of the heading of the sensor. We found that the reconstructed Larmor frequency for all heading angles corresponds accurately to the magnetic field when the two channels with different helicities are well balanced. The shot noise limited resolution strongly depends on the orientation in the external magnetic field. We demonstrated that this is related to the modified atom-light coupling responsible for both, a reduction of the photocurrent due to a reduction of spin polarization and a smaller light shift. The strong optical pumping leads to a strong contribution of power broadening to the magnetic resonance widths. This is added to the power broadening induced by the $B_{\_}1$ field. The orientation of the latter has a strong influence to the amplitude and width of the resonance signal. Nevertheless, in terms of shot noise limited resolution the two coils produce similar dependencies on the orientation angle.

Our experimental results can be explained in frame of light-shift calculations as well as Bloch equations. It was necessary to include the modified optical pumping into the latter by introducing heading dependent relaxation and excitation rates to a two-level model. This allows us to qualitatively and to a large extent also quantitatively describe our experimental results. We note, that such a model is not restricted to the description of the LSD-Mz mode, discussed in this work.

Finally, we conclude that a sensor operated in the LSD-Mz regime should be roughly $\pm 20^{\circ}{}$ aligned to the direction of the measured magnetic field vector to achieve a good field sensitivity. Notably, for well balanced beams the actual error in the Larmor frequency induced by light shift to each individual magnetic resonance is canceled by their subtraction. That gives a reasonable flexibility in adjusting the heading of an OPM sensor in an external magnetic field which allows their operation in Earth field strengths.

List of abbreviations

[TABLE]

Acknowledgements

We acknowledge the support of the clean room staff under supervision of U. Hübner in the Center of Micro- and Nanotechnologies at Leibniz IPHT in the fabrication of the cesium vapor cell used in the experiment. We thank the Leibniz Association’s Open Access Publishing Fund for financial support of this publication.

Funding

The authors acknowledge the financial support by the Federal Ministry of Education and Research (BMBF) of Germany under Grant No. 033R130E (DESMEX). The project 2017 FE 9128, funded by the Free State of Thuringia, was co-financed by European Union funds under the European Regional Development Fund (ERDF). This work was conducted using the infrastructure supported by the Free State of Thuringia under Grant No. 2015 FGI 0008 and co-financed by European Union funds under the European Regional Development Fund (ERDF). This work has received funding from the German Research Foundation (DFG) under Grant No. SCHU 2845/2-1; AOBJ 621093.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The data generated and analysed during the current study is available from the corresponding author on reasonable request.

Author’s contributions

RIJ, VS, and RS planned the experiment. RIJ, VS, and GO created the experimental setup and carried out the measurements. GO made the data analysis and theoretical calculations. All authors participated in interpretation of results and writing the manuscript. RS supervised the project.