EPJ Quantum Technology
Performance analysis of an optically pumped magnetometer in Earth's magnetic field
Gregor Oelsner ORCID: orcid.org/0000-0003-2002-19791,
Volkmar Schultze1,
Rob IJsselsteijn2 &
Ronny Stolz1
EPJ Quantum Technology volume 6, Article number: 6 (2019) Cite this article
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 resolutions at such field strengths. For this so-called light-shift dispersed \(M_{z}\) (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.
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 [1]. 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 [2,3,4].
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"[5]. 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 [6]. Suppression of the heading error is intensively studied [7,8,9] and it has important consequences in the application of OPMs [10, 11].
The most sensitive magnetic field sensors are based on superconducting quantum interference devices [12, 13]. They have been optimized to allow sub-femtotesla field gradient resolution even in Earth's magnetic field [14] and today, besides others, are used for geomagnetic and archeological explorations [15,16,17]. 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) [18,19,20] and the LSD-Mz [21] 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.
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. The diameter of the two cylindrical cells is 11 mm and their distance 16 mm. 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. For the used cell the buffer gas pressure is about 200 mbar as necessary for the LSD-Mz mode. A detailed description of the cell fabrication can be found in [22]. A photograph of the used cell arrangement is shown in Fig. 1(b).
(a) The experimental setup (see also Ref. [6]) is mounted on a rotational table inside of a Helmholtz coil system and a mu-metal shielding. The orientation of the setup concerning the created magnetic field is adjusted by a cable pull. (b) The vapor cell used for our experimental investigations features two circularly shaped active volumes and a rectangular reservoir. (c) The schematic drawing of the setup viewed from the side for rotating angles of \(\alpha = 0\) includes the optics for generating two circularly polarized laser beams. The axis of rotation is parallel to the y-axis. The pump (PL) and heat laser (HL) are colored in red and blue, respectively. The whole optical setup can be rotated inside the magnetic field \(\vec{B}_{0}\) that is fixed along the z-axis
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(a). 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 (KL) 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 centers of the active volumes 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 placed close to the laser on the optical table (not shown in Fig. 1). 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 for a suitable optical density of cesium atoms.
The whole optical installation is mounted inside of a three layer mu-metal shielding (see Fig. 1(a) and Ref. [23]). 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 \mbox{T}\) along the horizontal z-axis that corresponds to the symmetry axis of the cylinder-shaped shields. The normal of the rotational table and thus the vertical 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 Ω and ν. Note, the \(\vec{B}_{1}\) field created by the y-coil is perpendicular to both, laser light direction k⃗ and magnetic field \(\vec{B}_{0}\) for each rotation angle α. 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 α between \(\vec{B}_{0}\) and the propagation direction of the laser light 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 α 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}\) [21]. 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 \(\nu _{L} = \gamma B_{0}\). Here \(\gamma = 3.5\text{ Hz/nT}\) denotes the gyromagnetic ratio of cesium. This substraction also reduces the common noise present on the two signals, 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.
Measured magnetic resonance for both helicities as well as difference signal. The latter yields a large transfer function exploited in the LSD-Mz mode. This example curves are measured at an angle \(\alpha = 0\) with the y-coil
By fitting this measurement curves of the photocurrent I as a function of rf-frequency ν with a Lorentzian function
$$ I = I(\nu ) + I_{dc} = I_{0} \frac{\Delta \nu ^{2}}{4(\nu -\nu _{0})^{2}+ \Delta \nu ^{2}} + I_{dc}, $$
we extract the DC-photocurrents \(I_{dc}\), the amplitudes \(I_{0}\) and widths (FWHM) Δν 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
$$ B_{sn} = \frac{\sqrt{2e (I_{dc}^{\sigma _{+}}(\nu _{L})+ I_{dc} ^{\sigma _{-}}(\nu _{L}) )}}{\gamma s} , $$
by error propagation of the two channels' shot noises where the photocurrents \(I_{dc}\) are taken at the Larmor frequency and e is the elementary charge. This value sets the theoretical lower limit for the field resolution of an OPM operated in the LSD-Mz mode. Please note that a possible sensor operated in the \(M_{z}\)-mode suffers from low-frequency noise due to the direct measurement of the photo-current. However, by the use of two laser beams from the same source, noise contributions such as the amplitude noise of the laser can be reduced to a large extent as we will show elsewhere. The heading characteristics of the fitting parameters of Eq. (1) as well as of the sensitivity are in the following discussed.
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 [6, 7, 9, 24] and we use the method presented in Ref. [6] 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
$$ H= \frac{h \nu _{0}}{2}\sigma _{z} + h\varOmega \cos 2\pi \nu t [ \cos \alpha \sigma _{x} + \sin \alpha \sigma _{z} ] $$
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), \(\varOmega = \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{\varOmega }{\nu }\sin 2\pi \nu t \cos \alpha \). As constructed, \(U_{1}\) commutes with \(\sigma _{z}\) and the term \(i\hbar \dot{U}_{1} U_{1}^{\dagger }\) removes the diagonal coupling [25]. The transformation therefore results in
$$ H^{\prime }= \frac{h \nu _{0}}{2} \sigma _{z} + h \varOmega \cos 2\pi \nu t \cos \alpha \bigl[ e^{ig\sin 2\pi \nu t} \sigma _{+} + e^{-ig \sin 2\pi \nu t} \sigma _{-} \bigr]. $$
Here, we used \(\sigma _{\pm }= 0.5 ( \sigma _{x} \pm i \sigma _{y} )\) and the abbreviation \(g = \frac{\varOmega }{\nu }\cos \alpha \). Although we could continue with the Jacobi–Anger expansion and finally find the stationary terms in a frame rotating with ν, we note that the ratio \(\varOmega /\nu \) is small in our experimental realization. Thus, we can summarize the terms in the brackets to \(\sigma _{x}\) and solve the Bloch equations in a frame rotating with ν 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α as
$$ \langle \sigma _{z} \rangle = -\frac{\varGamma _{r} \varGamma _{\varphi }^{\prime }}{\varGamma _{r} \varGamma _{\varphi }^{\prime }+ \varOmega ^{2} \cos ^{2}\alpha }. $$
Above we introduced the rates \(\varGamma _{r}\) and \(\varGamma _{\varphi }\) as respected inverse relaxation and coherence times \(T_{1}\) and \(T_{2}\), as well as \(\varGamma _{\varphi }^{\prime }= (\varGamma _{\varphi } ^{2}+\delta ^{2})/\varGamma _{\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. [26] and [27]. Namely, we set for the respective excitation \(\varGamma _{e}\) and decay rate \(\varGamma _{r}\)
$$ \begin{aligned} &\varGamma _{r} = \frac{\gamma _{r}}{2} + \frac{\varOmega _{p} \vert \cos \alpha \vert }{2} (1-\cos \alpha ), \\ &\varGamma _{e} = \frac{\gamma _{r}}{2} + \frac{\varOmega _{p} \vert \cos \alpha \vert }{2} (1+\cos \alpha ). \end{aligned} $$
Here we assumed a repopulation rate of the ground state levels \(\gamma _{r}/2\) and an effective pumping rate between the two levels \(\varOmega _{p}\) proportional to the laser field \(\varOmega _{L}\) (see below). 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. [6]) as well as an additional factor of \(\vert \cos \alpha \vert \). Also, we did not include the linear polarized component scaling with sinα 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
$$ I(\nu ) = I_{0} \langle \sigma _{z} \rangle = I_{0}\frac{( \varGamma _{e}-\varGamma _{r}) \varGamma _{\varphi }^{\prime }}{(\varGamma _{r}+\varGamma _{e}) \varGamma _{\varphi }^{\prime }+ c^{2}(\alpha )\varOmega ^{2} }. $$
Here, the function \(c(\alpha )\) is constantly one for the x-coil and cosα for the y-coil. As constructed, our model accurately describes the dc-photocurrent as shown in Sect. 3. Namely, if \(\delta \rightarrow \pm \infty \) or equivalently \(\varOmega \rightarrow 0\), we find
$$ I_{dc} = I_{0} \frac{\varGamma _{e}-\varGamma _{r}}{\varGamma _{r}+\varGamma _{e} }= I _{0}\frac{\cos ^{2} \alpha }{ \vert \cos \alpha \vert + p_{1}}, $$
with the dimensionless fitting parameter \(p_{1} = \gamma _{r}/\varOmega _{p}\). The above function changes with decreasing optical pumping from a cos to a cos2 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 \varOmega _{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
$$ I(\nu ) = I_{dc} \biggl( 1- \frac{c^{2}(\alpha )\varOmega ^{2}}{ (\varGamma _{e}+\varGamma _{r} )\varGamma _{\varphi }^{\prime }+ c^{2}(\alpha ) \varOmega ^{2}} \biggr). $$
The resonance amplitude is accordingly the value of the second term in above equation at resonance, namely at \(\nu = \nu _{0}\). It is
$$ i_{0} = I(\nu = \nu _{0})/I_{dc} = \frac{c^{2}(\alpha )}{c^{2}(\alpha )+p_{2} ( p_{1}+\cos \alpha )}, $$
where we introduced a second dimensionless fitting parameter \(p_{2} = \varGamma _{\varphi }\varOmega _{p} / \varOmega ^{2}\).
Finally, the full width at half maximum (FWHM) of the resonance curve of Eq. (7) is given by
$$ \Delta \nu = 2\varGamma _{\varphi }\sqrt{1+\frac{c^{2}(\alpha )}{p_{2}(p _{1}+ \vert \cos \alpha \vert )}+ p_{3}}. $$
This relation describes the quadratic addition of the natural linewidth \(\varGamma _{\varphi }\) and the power broadening by the \(B_{1}\) field
$$ \Delta \nu _{B_{1}} = \frac{c^{2} \varGamma _{\varphi }}{p_{2} (p_{1}+ \vert \cos \alpha \vert ))} = \frac{c^{2}\varOmega ^{2}}{ \varGamma _{\varphi } ( \gamma + \varOmega _{p} \vert \cos \alpha \vert )} $$
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.
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.
Extracted center frequencies of the magnetic resonances for both circular polarizations of the laser as a function of the heading angle α. Also the mean frequencies of the two helicities are given. The latter roughly corresponds to the Larmor frequency measured by the LSD-Mz mode. The diamonds and dots correspond to a \(B_{1}\) field that is applied respectively in plane of rotation and perpendicular to it. Away from \(\pm \pi /2\) the two coils produce the same result. Additionally, theoretical curves describing the expected light shift are added as solid orange and blue lines
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 [6]
$$ \vec{D}\cdot \vec{E} = \frac{E_{0}}{2\sqrt{2}} e^{i ( \vec{k} \vec{r}-2\pi \nu _{L} t )} \bigl( [ \cos \alpha \pm 1 ] D_{+} + [ \cos \alpha \mp 1 ] D_{-} +2 \sin \alpha D_{z} \bigr). $$
Therein \(E_{0}\) is the amplitude of the electric field, k⃗ the k-vector of the laser, 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 α 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} + \{ +1,-1,0 \} \), 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 [24]. It is strongest close to angles of zero and ±π 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. [6]. 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 \mbox{T}\), detunings of \(\delta _{F^{\prime }= 4} = -8\text{ GHz}\) and \(\delta _{F^{\prime }= 3} = -9.2\text{ 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 \(\varGamma _{\mathrm{opt}} \approx 4\text{ GHz}\), and optical driving amplitudes in frequency units \(\varOmega _{L} = 3.15\text{ MHz}\) and 3.45 MHz for \(\sigma _{+}\) and \(\sigma _{-}\) polarized beam, respectively. The last values corresponds to the on-resonance Rabi frequencies [28] 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 ±π, 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
$$ \varOmega _{L} = \frac{1}{h} \sqrt{\frac{I_{L}}{4cn\epsilon _{0}}} \bigl\langle J=1/2 \bigl\Vert \vec{D} \bigr\Vert J^{\prime }= 1/2 \bigr\rangle . $$
Here the vacuum permittivity \(\epsilon _{0}\), speed of light c, reduced transition dipole moment \(\langle J=1/2 \Vert \vec{D} \Vert J^{\prime }= 1/2 \rangle \) [29], 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.
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.
Photocurrent as a function of the heading angle. The diamonds and dots correspond to the data recorded with the x and y coil producing the \(B_{1}\) field, respectively. In the picture the difference Δ of \(\sigma _{-}\) minus \(\sigma _{+}\) values are added to visualize the balancing. The solid lines are calculated according to Eq. (8). The dotted gray lines correspond to the functions \(0.38~\mathrm{mA} \times \vert \cos \alpha \vert \) and \(0.38~\mathrm{mA} \times \cos ^{2} \alpha \) for a better visual comparison of the results to these functions
As demonstrated in the figure, we observe a dependence that roughly follows a \(\vert \cos \alpha \vert \) or \(\vert \cos \alpha \vert ^{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 D⃗ with the laser's electric field 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} + i D_{y}\) and \(D_{-} = D_{x} - i D_{y}\).
Our model (Eq. (8)) predicts a change from a cos2α dependence to one proportional to \(\vert \cos \alpha \vert \) when the effective pumping amplitude \(\varOmega _{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 α 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~\mbox{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 \varOmega _{p}\) that usually is not included in Bloch equations.
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.
Normalized magnetic resonance height \(I_{0}/I_{dc}\) as a function of the heading angle. Again the diamonds and dots correspond to the x and y orientation of the \(B_{1}\) field, respectively. The calculation of the solid and dashed theoretical lines for the respective x and y coil is explained in the main text. Their colors are adjusted to the corresponding data points
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 \(\vert \cos \alpha \vert \), 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 normalized resonance amplitude driven by the x-coil is strongly influenced by the cos function given for the effective \(B_{1}\)-field amplitude Ω. Namely, the resonance height is maximal at \(\alpha = 0\) and π 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 the 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.
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 \(\varGamma _{\varphi }= 350\text{ Hz}\) are plotted together with the experimental data points in Fig. 6.
FWHM of the magnetic resonances as a function of the angle α. Again diamonds and dots correspond to the use of x and y coil for applying the rf-field, respectively. In green, we added the corresponding difference between \(\sigma _{+}\) and \(\sigma _{-}\) result. The solid lines correspond to calculation results following Eq. (11) in the same color as the corresponding data points
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. [30])
$$ \Delta \nu _{\mathrm{laser}} = \frac{\varOmega _{L}^{2}}{2\varGamma _{\varphi }\varGamma _{\mathrm{opt}}} = 3.5, $$
where we used the values for \(\varOmega _{L}\) of the \(\sigma _{+}\) beam and \(\varGamma _{\mathrm{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.
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 Sect. 2.1. The respective experimental values as a function of the angle α for the two rf-coil-orientations are presented in Fig. 7. Our experiment shows a minimal shot-noise limited resolution for the used parameters of \(20~\mbox{fT/}\sqrt{\text{Hz}}\) already close to the optimal conditions [21]. 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.
Shot-noise-limited resolution as a function of the heading angle. The black diamonds and dots correspond to the values extracted from the experimental curves for the x and y coil, respectively. The calculated, expected dependencies in case of optimal signal balancing from the theoretical model are added as solid orange and green line
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 ±20∘ 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 α 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)
$$ \frac{s}{2} = \biggl\vert \frac{dI}{d\nu } \biggr\vert = \frac{8 I_{0} \Delta \nu ^{2} (\nu -\nu _{0})}{ ( 4 (\nu - \nu _{0} )^{2} + \Delta \nu ^{2} )^{2}} $$
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 \vert \frac{dI}{d \nu } \vert (\nu _{0} + \nu _{LS})\). Substituting into the equation for the shot-noise-limited resolution (2) results in
$$ B_{sn} = \sqrt{\frac{e}{I_{dc}}}\frac{1}{8\gamma i_{0} } \frac{ ( 4\nu _{LS}^{2} + \Delta \nu ^{2} )^{2}}{\Delta \nu ^{2} \nu _{LS}}, $$
with the parameters \(I_{dc}\), \(i_{0}\), and Δν 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 and 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 ±20∘ 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.
beam splitter
DP:
deflecting prism
bandpass filter
FWHM:
full width half maximum
heat laser
KL:
collimating lens
LN:
light narrowing
linear polarizer
LSD-Mz:
Light-shift dispersed \(M_{z}\)
OPM:
optically-pumped magnetometer
polarizing beam splitter
photo diode
pump laser
rf:
Cohen-Tannoudji C, Diu B, Laloe F. Quantum mechanics. Volume I. New York: Wiley; 1977.
Leger J-M, Bertrand F, Jager T, Prado ML, Fratter I, Lalaurie J-C. Swarm absolute scalar and vector magnetometer based on helium 4 optical pumping. Proc Chem. 2009;1(1):634–7. https://doi.org/10.1016/j.proche.2009.07.158.
Nabighian MN, Grauch VJS, Hansen RO, LaFehr TR, Li Y, Peirce JW, Phillips JD, Ruder ME. 9. Magnetic exploration methods. In: Geophysics today. Tulsa: Society of Exploration Geophysicists; 2010. p. 183–213. https://doi.org/10.1190/1.9781560802273.ch9.
Korth H, Strohbehn K, Tejada F, Andreou AG, Kitching J, Knappe S, Lehtonen SJ, London SM, Kafel M. Miniature atomic scalar magnetometer for space based on the rubidium isotope 87 rb. J Geophys Res Space Phys. 2016;121(8):7870–80. https://doi.org/10.1002/2016ja022389.
Budker D, Kimball DFJ, editors. Optical magnetometry. New York: Cambridge University Press; 2013.
Oelsner G, Schultze V, IJsselsteijn R, Wittkämper F, Stolz R. Sources of heading errors in optically pumped magnetometers operated in the Earth's magnetic field. Phys Rev A. 2019;99(1):013420. https://doi.org/10.1103/physreva.99.013420.
Jensen K, Acosta VM, Higbie J, Ledbetter M, Rochester S, Budker D. Cancellation of nonlinear Zeeman shifts with light shifts. Phys Rev A. 2009;79(2):023406. https://doi.org/10.1103/PhysRevA.79.023406.
Bao G, Wickenbrock A, Rochester S, Zhang W, Budker D. Suppression of the nonlinear Zeeman effect and heading error in Earth-field-range alkali-vapor magnetometers. Phys Rev Lett. 2018;120(3):033202. https://doi.org/10.1103/PhysRevLett.120.033202.
Hu Q-Q, Freier C, Sun Y, Leykauf B, Schkolnik V, Yang J, Krutzik M, Peters A. Observation of vector and tensor light shifts in rb87 using near-resonant, stimulated Raman spectroscopy. Phys Rev A. 2018;97(1):013424. https://doi.org/10.1103/physreva.97.013424.
Colombo S, Dolgovskiy V, Scholtes T, Grujić ZD, Lebedev V, Weis A. Orientational dependence of optically detected magnetic resonance signals in laser-driven atomic magnetometers. Appl Phys B. 2016;123(1):35. https://doi.org/10.1007/s00340-016-6604-8.
Weis A, Bison G, Grujić ZD. In: Grosz A, Haji-Sheikh MJ, Mukhopadhyay SC, editors. Magnetic resonance based atomic magnetometers. Cham: Springer; 2017. p. 361–424. https://doi.org/10.1007/978-3-319-34070-8_13.
Clarke J, Braginski AI, editors. The SQUID handbook: fundamentals and technology of SQUIDs and SQUID systems. Weinheim: Wiley-VCH; 2004.
Schmelz M, Stolz R. In: Grosz A, Haji-Sheikh MJ, Mukhopadhyay SC, editors. Superconducting quantum interference devices (SQUID) magnetometers. Cham: Springer; 2017. p. 279–311. https://doi.org/10.1007/978-3-319-34070-8_13.
Stolz R, Zakosarenko VM, Fritzsch L, Oukhanski N, Meyer H-G. Long baseline thin film SQUID gradiometers. IEEE Trans Appl Supercond. 2001;11(1):1257–60. https://doi.org/10.1109/77.919578.
Chwala A, Stolz R, IJsselsteijn R, Schultze V, Ukhansky N, Meyer H-G, Schüler T. SQUID gradiometers for archaeometry. Supercond Sci Technol. 2001;14(12):1111–4. https://doi.org/10.1088/0953-2048/14/12/327.
Foley CP. In: Clarke J, Braginski AI, editors. SQUID system issues. vol. 7. Weinheim: Wiley-VCH; 2004. p. 251.
Stolz R. In: Seidel P, editor. Geophysical exploration. Weinheim: Wiley-VCH; 2015. p. 1020. Chap. 9.3.4.
Scholtes T, Schultze V, IJsselsteijn R, Woetzel S, Meyer H-G. Light-narrowed optically pumped m x magnetometer with a miniaturized cs cell. Phys Rev A. 2011;84(4):043416. https://doi.org/10.1103/physreva.84.043416.
Guo Y, Wan S, Sun X, Qin J. Compact, high-sensitivity atomic magnetometer utilizing the light-narrowing effect and in-phase excitation. Appl Opt. 2019;58(4):734. https://doi.org/10.1364/ao.58.000734.
Fu Y, Liu Y, Yuan J. Light narrowing of cesium magnetic-resonance lines in a radio-frequency atomic magnetometer. AIP Adv. 2019;9(1):015304. https://doi.org/10.1063/1.5043231.
Schultze V, Schillig B, IJsselsteijn R, Scholtes T, Woetzel S, Stolz R. An optically pumped magnetometer working in the light-shift dispersed mz mode. Sensors. 2017;17(3):561. https://doi.org/10.3390/s17030561.
Woetzel S, Schultze V, IJsselsteijn R, Schulz T, Anders S, Stolz R, Meyer H-G. Microfabricated atomic vapor cell arrays for magnetic field measurements. Rev Sci Instrum. 2011;82(3):033111. https://doi.org/10.1063/1.3559304.
Schultze V, IJsselsteijn R, Meyer H-G. Noise reduction in optically pumped magnetometer assemblies. Appl Phys B. 2010;100(4):717–24. https://doi.org/10.1007/s00340-010-4084-9.
Mathur BS, Tang H, Happer W. Light shifts in the alkali atoms. Phys Rev. 1968;171(1):11–9. https://doi.org/10.1103/physrev.171.11.
Shevchenko SN, Oelsner G, Greenberg YS, Macha P, Karpov DS, Grajcar M, Hübner U, Omelyanchouk AN, Il'ichev E. Amplification and attenuation of a probe signal by doubly dressed states. Phys Rev B. 2014;89(18):184504. https://doi.org/10.1103/physrevb.89.184504.
Dréau A, Lesik M, Rondin L, Spinicelli P, Arcizet O, Roch J-F, Jacques V. Avoiding power broadening in optically detected magnetic resonance of single NV defects for enhanced dc magnetic field sensitivity. Phys Rev B. 2011;84(19):195204. https://doi.org/10.1103/physrevb.84.195204.
Shi Y, Scholtes T, Grujić ZD, Lebedev V, Dolgovskiy V, Weis A. Quantitative study of optical pumping in the presence of spin-exchange relaxation. Phys Rev A. 2018;97(1):013419. https://doi.org/10.1103/physreva.97.013419.
Oelsner G, Macha P, Astafiev OV, Il'ichev E, Grajcar M, Hübner U, Ivanov BI, Neilinger P, Meyer H-G. Dressed-state amplification by a single superconducting qubit. Phys Rev Lett. 2013;110(5):053602. https://doi.org/10.1103/PhysRevLett.110.053602.
Steck DA. Cesium D line data. https://steck.us/alkalidata/cesiumnumbers.1.6.pdf.
Loudon R. The quantum theory of light. Oxford science publications. New York: Oxford University Press; 2000. p. 72–9.
We acknowledge the support of the clean room staff under supervision of U. Hübner at the Competence Center Micro- and Nanotechnologies of the Leibniz IPHT in the fabrication of the cesium vapor cell used in the experiment.
The data generated and analysed during the current study is available from the corresponding author on reasonable request.
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. The publication of this article was funded by the Open Access Fund of the Leibniz Association.
Leibniz Institute of Photonic Technology, Jena, Germany
Gregor Oelsner
, Volkmar Schultze
& Ronny Stolz
Supracon AG, Jena, Germany
Rob IJsselsteijn
Search for Gregor Oelsner in:
Search for Volkmar Schultze in:
Search for Rob IJsselsteijn in:
Search for Ronny Stolz in:
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. All authors read and approved the final manuscript.
Correspondence to Gregor Oelsner.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Oelsner, G., Schultze, V., IJsselsteijn, R. et al. Performance analysis of an optically pumped magnetometer in Earth's magnetic field. EPJ Quantum Technol. 6, 6 (2019) doi:10.1140/epjqt/s40507-019-0076-9
DOI: https://doi.org/10.1140/epjqt/s40507-019-0076-9
Optically pumped magnetometers
Heading error
Optical pumping
Bloch equations