## Investigating the hyperfine systematic error and relative phase in low spin-polarization alkali FID magnetometers

D. P. Hewatt, 1, 2, 3, ∗ M. Ellmeier, 2 C. Kiehl, 1, 3 T. S. Menon, 1, 3 J. W. Pollock, 2 C. A. Regal, 1, 3 and S. Knappe 2, 4, 5

1 Department of Physics, University of Colorado, Boulder, Colorado 80309, USA

2 Paul M. Rady Department of Mechanical Engineering,

University of Colorado, Boulder, Colorado 80309, USA

3 JILA, National Institute of Standards and Technology and University of Colorado, Boulder, Colorado 80309, USA

4 FieldLine Medical, Boulder Colorado 80301, USA

5 FieldLine Industries, Boulder Colorado 80301, USA

(Dated: March 17, 2025)

Alkali-metal optically pumped magnetometers that utilize vapor cells containing high buffer-gas pressure and operate at elevated temperatures are prone to inaccuracies arising from the overlap of the average F = I +1 2 and / F = I -1 2 ground-state Zeeman resonances. We employ density-matrix / simulations and experiments to investigate how this hyperfine systematic error varies with spin polarization in a 87 Rb free-induction-decay (FID) magnetometer. At low spin polarizations P ≤ 0 5, . this effect causes single-frequency magnetic-field extraction techniques to exhibit inaccuracies up to approximately 3.5 nT. Density-matrix simulations reveal that this bias can be traced to the relative amplitude and phase between the F = I ± 1 2 hyperfine ground-state manifolds in the FID spin / precession signal. We show that this systematic error can be mitigated using either a doublefrequency fitting model that accounts for the relative amplitude and phase or synchronous-pulse pumping that minimizes the F = 1 contribution to the FID signal. Theoretical simulations predict accuracies within 0.5 nT for both techniques across a wide range of spin polarizations, suggesting a sevenfold enhancement over single-frequency extraction methods. Our experiments validate this, showcasing a variation in the extracted field below 1 nT, a 3.5-fold improvement compared to singlefrequency extraction methods. Furthermore, the mitigation techniques demonstrate agreement of the extracted magnetic field within 1.5 nT.

## I. INTRODUCTION

Scalar optically pumped magnetometers (OPMs), known for their high sensitivity, robustness, and ability to be miniaturized, have found widespread applications in fields such as fundamental physics [1-3], space magnetometry [4, 5], and biomedical imaging [6-9]. The accuracy of alkali-metal OPMs is limited by frequency shifts of the unresolved Zeeman sublevels, caused by effects such as the nonlinear Zeeman (NLZ) shift (Fig. 1(a)). Microfabricated vapor cells often contain buffer gas to reduce depolarizing wall collisions and operate at elevated temperatures (80 ◦ C-120 ◦ C) to increase atomic density, enhancing sensitivity. The addition of buffer gas broadens and shifts the optical resonances hindering hyperfineselective pumping, while elevated temperatures increase spin-exchange collisions, broadening the magnetic resonance. Since the F = I ± 1 / 2 manifolds have different hyperfine Land´ e g -factors, there is a frequency difference between the average Zeeman resonances of the groundstate hyperfine manifolds [10, 11]. In regimes of high buffer gas and elevated temperatures, the average Zeeman resonances of the F = I ± 1 / 2 hyperfine manifolds become spectrally unresolved, introducing an additional source of inaccuracy we call the hyperfine systematic error.

Techniques like analytical correction terms [10], double-pass configurations [12], light-polarization modulation [13], double-modulated synchronous pumping [14], and spin-locking [15] can mitigate NLZ effects. While these approaches effectively minimize NLZ shifts, they do not account for the hyperfine systematic error, limiting their accuracy in regimes of low initial atomic spin polarization P ≤ 0 5, . where this error becomes nonnegligible. Techniques utilizing higher-order polarization modes are inherently immune to the NLZ effect and hyperfine systematic error [16, 17]. However, the transition amplitudes of these higher-order modes are significantly suppressed in vapor cells with high buffer-gas pressure, limiting the sensitivity [18, 19]. To address the hyperfine systematic error, Lee et al. proposed averaging precession signals between two orthogonal probe beams, but this approach requires a second probe laser, increasing sensor complexity [10].

∗ Dawson.Hewatt@colorado.edu

Our analysis reveals that unresolved Zeeman resonances between the F = I ± 1 2 manifolds cause tech-/ niques that use a single frequency to extract the magnetic field to produce inaccuracies exceeding 1 nT for atomic spin polarizations below 0.5 (Fig. 2(a)) [10, 20]. Lowpower consumption systems, such as compact sensors, often lack access to the high-powered lasers required for generating high spin polarizations in optically dense media used in microfabricated cells. Therefore, these magnetometers are susceptible to these accuracy constraints. This systematic error is a barrier to applications like navigation, fundamental physics, and space magnetometry,

where the demand for accuracies less than or equal to 1 nT is commonplace [21-23].

There are many OPM detection methods including Bell-Bloom [24, 25], M x and M z [3, 26], coherent population trapping (CPT) [27-29], Rabi [30-32], spin-exchange-relaxation-free regime [33, 34], and freeinduction-decay regime (FID) [35-37]. Hyperfine measurement approaches such as CPT and microwave interrogation techniques can overcome these systematic inaccuracies because they can resolve selected transitions [32, 38, 39]. However, these hyperfine techniques exhibit limited sensitivity when compared to methods that measure Zeeman resonances. Spin-exchange-relaxation-free OPMs are well known for displaying sub-fT/(Hz) 1 / 2 sensitivities, but can only operate at near-zero fields [4042]. In contrast, FID magnetometers can operate at geomagnetic fields (approximately 50 µ T) and have demonstrated fT/(Hz) 1 / 2 -level sensitivities [10, 20, 22, 43]. Unlike driven magnetometers, FID sensors do not require prior knowledge of the magnetic field to operate. This can improve the magnetic field slew rate by eliminating the need for an additional control loop [35]. Additionally, the pulsed operation used in FID magnetometers minimizes light shifts by off-resonantly probing the atomic state.

In an FID measurement scheme, the atoms are optically pumped, and then freely precess in the dark, sampling the six Zeeman splittings f i shown in Fig. 1(a). Although there are various ways to implement an FID magnetometer, a schematic of the apparatus used in this study is shown in Fig. 1(b). The magnetometer consists of a 3 × × 3 2 mm 3 micro-fabricated vapor cell containing enriched 87 Rb and heated to roughly 108 ◦ C via resistive heaters. Approximately 1000 Torr of N 2 buffer gas is added to the cell to limit wall collisions. Circularly polarized laser light, resonant with the 87 Rb D 1 transition, pumps the atoms into one of the stretched states | F = 2, m F = ± ⟩ 2 . Free-induction-decay magnetometers commonly employ one of two primary pumping methods to create an initial spin polarization: single-pulse pumping with one uninterrupted pulse of pulse length T p or synchronous-pulse pumping with many short pulses, amplitude modulated (or frequency modulated) at the Larmor frequency (Fig. 1(d)) [35]. A balanced polarimeter measures the precession of the atomic spin ensemble by detecting the Faraday rotation of a far-detuned, linearly polarized probe beam (Fig. 1(b)). Here we use a collinear pump-probe geometry. The average Larmor frequency f L of the spin precession is extracted to determine the magnetic field.

As depicted in Fig. 1(c), optical pumping affects the relative amplitudes of the various frequency components within FID signals and consequently the accuracy of the magnetometer. Higher optical pumping rates R OP generate larger initial spin polarizations. As the spin polarization decreases, the amplitude of the F = I -1/2 frequency component increases, shifting the average Larmor frequency. Furthermore, we also find that the rela- tive phase between the F = I ± 1 / 2 Zeeman resonances depends on optical pumping and can affect the extracted field.

In this paper we use simulation and experiment to investigate how the hyperfine systematic error affects the accuracy of magnetic-field extraction as a function of initial spin polarization in a 87 Rb FID magnetometer. At low spin polarizations P ≤ 0 5, we demonstrate . that implementing a double-frequency fitting model more accurately extracts the magnetic field. Conversely, synchronous-pulse pumping can enhance the initial spin polarization and consequently the accuracy of singlefrequency extraction techniques. Theoretically, we show that both double-frequency fitting and synchronous-pulse pumping result in an enhancement of the FID accuracy by a factor of 7 compared to single-frequency fitting of FID signals produced by single-pulse pumping. Experimentally, we observe agreement within 1.5 nT between both the double-frequency and synchronous-pulse pumping techniques.

## II. THEORETICAL MODEL

We describe the ground-state dynamics of this system using an 8 × 8 density matrix, ρ , in the total atomic spin basis basis | F , m F ⟩ z , where the z indicates that the basis is chosen in the ˆ direction of the laboratory frame. z The Hamiltonian that characterizes the hyperfine and Zeeman interactions of the 87 Rb atom is given by

$$\underset { \text{rel/v} } { \text{ning} } \underset { \text{rel/v} } { H } = & R _ { z } ( \alpha ) R _ { y } ( \beta ) \{ \mathcal { M } \left [ A _ { \text{hfs} } \mathbf I \cdot \mathbf S + \mu _ { B } ( g _ { s } S _ { z } + g _ { i } I _ { z } ) B \right ] \mathcal { M } ^ { \dagger } \\ \underset { \text{larg/r} } { \text{ gas} } \underset { \text{al/r} } { - I _ { 2 } \Delta E \} R _ { y } ( \beta ) ^ { \dagger } R _ { z } ( \alpha ) ^ { \dagger }, }$$

(1)

where A hfs is the hyperfine splitting constant, I is the nuclear spin, S is the electron spin, B is the dc magnetic field strength, µ B is the Bohr magneton, g s and g i are the electron and nuclear g -factors, respectively, as defined by [11], I 2 is the F = 2 identity operator, and ∆ E = h (6 834 . GHz) is the hyperfine splitting energy. Because the Hamiltonian is not generally diagonal in the | F , m F ⟩ z basis unless the magnetic field is in the ˆ direcz tion, we define M as the operator that diagonalizes the hyperfine and Zeeman terms [32]. This serves to maintain the NLZ effects during the rotating-wave approximation (RWA). The first term in Eq. (1) describes the hyperfine splitting of the 87 Rb atoms and the second term describes the Zeeman interaction between the atoms and the magnetic field. The third term in this Hamiltonian eliminates the high-energy scale caused by the hyperfine splitting, thereby enhancing computational efficiency. We rotate the system into the direction of the magnetic field by applying the rotation operators R z ( α ) = e -iF z α and R y ( β ) = e -iF y β (Fig. 1(b)). In 87 Rb, the Zeeman effect causes splitting of the m F levels with a gyromagnetic ratio of approximately 7 Hz/nT. The energy eigenvalues of the F = I ± 1 / 2 sublevels are given by the Breit-Rabi

FIG. 1. (a) Level diagram for the 87 Rb ground state. Free-induction-decay magnetometers measure the Zeeman resonance frequencies f i shown in green. (b) Magnetometer apparatus and example FID signal. Here, HWP is the half-wave plate; PBS is a Wollaston polarizing beam splitter, and BP is a balanced polarimeter. The magnetic field B is specified by a magnitude B azimuthal angle α , and polar angle β . The bottom graph shows an FID signal (black) and single-frequency fit (red). (c) Normalized Fourier transforms of simulated FID signals, with black solid and gray dashed lines representing high (100 R OP ) and low ( R OP ) pumping rates, respectively. These simulations neglect spin exchange so that the F = I ± 1 / 2 resonances could be resolved. The green vertical lines are the individual Zeeman transitions at 50 µ T, while the purple vertical lines are the average Larmor frequencies for the F = I ± 1 2 hyperfine manifolds. / (d) FID pumping schemes for single-pulse pumping (top) and synchronous-pulse pumping (bottom). Single-pulse pumping uses one pump pulse with a pumping rate R OP and pulse length T p . Synchronous-pulse pumping modulates the amplitude of the pump laser at the average Larmor frequency f L . Note that the time axes for these graphs are not to scale. After pumping, the atoms precess in the dark for a time T free .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00898v2-figure-08b69079-d7bd-44ea-9d97-de761ae34aa2.jpg)

formula [44]

$$E _ { m _ { F } } = \frac { \Delta E } { 2 } \left ( \frac { - 1 } { 2 I + 1 } + \frac { 2 g _ { I } m _ { F } x } { g _ { s } - g _ { I } } \pm \sqrt { x ^ { 2 } + \frac { 4 x m _ { F } } { 2 I + 1 } + 1 } \right ), \quad \begin{array} { c c c } \text{The r} \\ \text{the n} \\ \text{the n} \\ ( 2 ) } \quad \text{gas at} \\ \text{where} \quad \varkappa \, - \, \int \nolimits ( \alpha \, - \, \partial \, R / \Lambda \, E \, \text{ and } I \, - \, 2 / 2 \, \text{ } \text{the $e$} \, \Pi \, \text{ $the cc} \end{array}$$

where x = µ B ( g s -g I ) B/ ∆ E and I = 3 2. / The FID signal contains six frequency components given by the difference between adjacent Zeeman sublevels (Fig. 1(a)).

The spin dynamics are approximated by the density matrix equation

$$\frac { d \rho } { d t } = & \frac { [ H, \rho ] } { i \hslash } + R _ { \text{SE} } [ \varphi ( 1 + 4 \langle \text{S} \rangle \cdot \text{S} ) - \rho ] + R _ { \text{SD} } ( \varphi - \rho ) \quad \text{est $SI$} \\ & + D \nabla ^ { 2 } \rho + R _ { \text{OP} } [ \varphi ( 1 + 2 \text{s} \cdot \text{S} ) - \rho ], \quad \text{$R_{\text{W} } =$} \\ \text{where to} = \alpha / 4 + \text{S} \cdot \text{s} \text{ is the river} \text{ nuclear art of the} \quad \text{film}$$

where φ = ρ /4 + S · ρ S is the purely nuclear part of the density matrix [45]. The first term is the time evolution of the Hamiltonian [Eq. (1)]. The second term describes the effect of spin-exchange collisions between alkali-metal atoms at a rate R SE . The third term expresses the effect of spin-destruction collisions between Rb with itself and

N 2 at a rate of R SD , where R SD = R Rb-Rb SD + R Rb-N 2 SD . The rates of these mechanisms are given by the definition, R SE (SD) ≡ n Rb (N ) 2 σ SE (SD) v r [33]. Here n Rb (N ) 2 is the number density of the alkali-metal atoms and buffergas atoms, respectively, σ SE (SD) is the cross section of the collision, and v r is the relative velocity of the atoms given by v r = √ 8 k B T/πM , where k B is the Boltzmann constant, T is the temperature, and M is the reduced mass of the colliding pair. The fourth term accounts for depolarizing wall collisions. We assume that the lowest spin-diffusion decay mode dominates the wall collision term such that we can make the approximation D ∇ 2 ρ → R W 2 ( ρ e -ρ ) where the wall-collision rate is R W = Dπ 2 /( l x + l 2 y + l 2 z ), D is the diffusion constant, l i is the cell length along the i th axis, and ρ e is the equilibrium density matrix containing equal populations in all states [46]. The last term describes optical pumping along a single axis in the limit of high buffer-gas pressure. This term is characterized by a pumping rate R OP and a mean photon spin s = sˆ z ∈ -( 1 1) [45]. , Spin-exchange

collisions are the dominant source of decoherence. Table I shows the decoherence rates used in this paper. Since we applied the RWA to the Hamiltonian, we must also apply the RWA to the decoherence and pumping terms by numerically averaging such that the high-frequency counterrotating terms are eliminated [31].

TABLE I. Decoherence rates for various collisions used in this paper. The cross sections are taken from ref. [31]. All rates were calculated using a cell temperature of T = 108 ◦ C, buffer-gas pressure of P N 2 = 1000 Torr, and inner cell dimensions of 3 × 3 × 2 mm . 3

| Collision    | Cross section                    | Rate (s - 1 )     |
|--------------|----------------------------------|-------------------|
| (Rb-Rb) SE   | σ SE = 1.90 × 10 - 18 m 2        | R SE = 6600       |
| (Rb-Rb) SD   | σ Rb-Rb SD = 1.77 × 10 - 21 m 2  | R Rb-Rb SD = 6    |
| (Rb-N 2 ) SD | σ Rb-N 2 SD = 1.44 × 10 - 26 m 2 | R Rb-N 2 SD = 226 |
| Wall         | D = 1.69 × 10 - 5 m 2 s - 1      | R W = 79          |

The simulation of a single FID signal is divided into two steps: pumping and free precession. During the pumping step, we solve Eq. (3) in 100-ns time steps using the equilibrium density matrix ρ e as the initial state and a nonzero pumping rate R OP for 10 µ s. After the pumping step, the pumped density matrix is used as the initial state to solve Eq. (3) with R OP = 0. We simulate the free precession in 100-ns time steps for 1 ms. Finally, we compute the expectation value of Faraday rotation for a far-detuned D 2 probe beam using the Faraday rotation operator [31]

$$\theta _ { F } = \frac { c r _ { e } f _ { D _ { 2 } } n _ { \text{Rb} } l } { 2 ( 2 I + 1 ) } \mathcal { L } ( \nu - \nu _ { D _ { 2 } } ) \begin{cases} - m _ { F } & ( F = 1 ) \\ + m _ { F } & ( F = 2 ). \end{cases} \quad ( 4 ) \quad \text{single} \\ \text{SFF} ) _ { \text{tin} \sigma }$$

Here, c is the speed of light, r e is the classical electron radius, f D 2 ≈ 2/3 is the D 2 oscillator strength, l = 2 mm is the propagation distance through the atomic media, and L ( ν -ν D 2 ) = ( ν D 2 -ν / ) [( ν -ν D 2 ) 2 +(Γ D 2 / 2) 2 ] is the value of the line shape with FWHM Γ D 2 and optical detuning ν -ν D 2 . For simplicity, we set the constants equal to 1. The Faraday rotation expectation value is ⟨ θ F ⟩ = Tr( ρθ F ).

## III. FITTING METHODS

There are many ways of extracting the frequency information from an FID signal such as frequency counters [10], demodulation and instantaneous phase retrieval techniques [36, 47], and least-squares fitting [22, 35, 48]. Fitting is often desirable for its flexibility, enabling the utilization of different models. Specifically, we focus on single- and double-frequency time-domain fitting models. Other frequency extraction methods that only account for one frequency, such as demodulation, exhibit behavior similar to the single-frequency fitting model discussed here.

We fit the Faraday rotation signals to the model

$$S _ { y } ^ { 2 } V ( t ) = \sum _ { i = 1 } ^ { N } V _ { i } \cos ( 2 \pi f _ { L, i } t + \phi _ { i } ) e ^ { - t / \tau _ { i } } + V _ { \text{dc} } ( t ), \quad ( 5 )$$

where V i , f L,i , ϕ i , and τ i are the ac amplitude, Larmor frequency, initial phase, and transverse spin-relaxation time for the i th frequency component, respectively. Here, V dc ( ) t = V 0 + V off e t/T off is the dc offset containing an overall offset, V 0 , and an exponentially decaying offset characterized by an amplitude, V off , and time constant, T off . The single-frequency fitting model includes one term in the sum ( N = 1). Given that the separation of the average Zeeman resonance frequencies f L,i in the F = 1 and F = 2 manifolds is much larger than the separation between adjacent Zeeman resonance frequencies (1.4 kHz versus 36 Hz at 50 µ T), we can more accurately model low spin-polarization FID signals with a double-frequency fitting model ( N = 2). This model includes a term for each hyperfine manifold. Additionally, we constrain the two frequencies in the double-frequency fitting model to the average Larmor frequency of each hyperfine manifold using the Breit-Rabi formula. This sets the magnetic field as a free parameter of the fit and prevents over-fitting.

## IV. SIMULATED PUMP-POWER DEPENDENCE

To investigate the accuracy of magnetic-field extraction from FID signals at varying spin polarizations, we simulate and analyze FID signals using three methods: single-pulse pumping with single-frequency fitting (SPSFF), single-pulse pumping with double-frequency fitting (SP-DFF), and synchronous-pulse pumping with single-frequency fitting (SY-SFF). We then compare the extracted-magnetic-field values to the known magneticfield strength. The obtained spin polarization depends on the mean photon spin s , pumping rate ( R OP ), and magnetic field B . Given that a high mean photon spin ( | s | ≥ 0 9) can be readily achieved using a quarter wave-. plate, we simulate FID signals as a function of pumping rate while keeping the mean photon spin constant. This allows us to compare our numerical simulations with experimental measurements as a function of pumping power because we can relate the pumping rate ( R OP ) to the pumping power, P pump , using the equation

$$P _ { \text{pump} } & = \frac { E _ { \text{ph} } A _ { b } R _ { \text{OP} } } { \sigma _ { \text{abs} } } \\ & = \frac { E _ { \text{ph} } A _ { b } \Gamma _ { D _ { 1 } } } { 2 c r _ { e } f _ { D _ { 1 } } } R _ { \text{OP} } \quad ( \text{on resonance} ),$$

where E ph is the energy of the pumping photons, A b is the area of the pumping beam written in terms of the beam waist, Γ D 1 is the FWHM of the D 1 absorption resonance, σ abs is the absorption cross section of the pumping light, and f D 1 is the D 1 oscillator strength [49, 50].

FIG. 2. (a) Simulated magnetic-field accuracy for SP-SFF (blue solid line), SP-DFF (green dash-dotted line), and SY-SFF (purple dashed line) as a function of pumping rate, R OP (bottom axis) and pump power (top axis). The pump power was calculated using Eq. (6) with a beam waist of 1.5 mm and FWHM of 21.3 GHz. The inset plots the same accuracy data as a function of initial atomic spin polarization P = 2 ⟨ S z ⟩ . (b) Spin polarization versus pumping rate for single-pulse pumping (cyan solid line) and synchronous-pulse pumping (orange dashed line). (c) Simulated relative phase ϕ rel versus pumping rate for single-pulse pumping (solid cyan) and synchronous-pulse pumping (dashed orange). The insets show the Faraday rotation θ F of the F = 2 (red) and F = 1 (blue) manifolds at R OP = 0.3 Ms -1 (left) and R OP = 20.1 Ms -1 (right). The Faraday rotation signals were calculated by splitting the Faraday rotation operator into submatrices for each hyperfine manifold and calculating the expectation value for each hyperfine manifold separately.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00898v2-figure-f05621f6-f701-4a05-b011-4ab4d22528ae.jpg)

The second line of this equation holds when the pumping light is resonant with the D 1 absorption line.

proximately 3.5 nT from the true field. As the pumping rate increases, the extracted magnetic field converges toward the true magnetic field.

This systematic error is caused by non-negligible contributions from the F = 1 state. The separation between the average Larmor frequencies for the F = 1 and F = 2 manifolds, | f L, 2 -f L, 1 | , is proportional to the magnetic field and is approximately 1.4 kHz at 50 µ T. Therefore, coherent Larmor precession of atoms in the F = 1 manifold shifts the average frequency extracted by a single-frequency fit. As the pumping rate increases, the F = 1 Larmor signal becomes negligible, making single-frequency extraction techniques accurate approximations. Consequently, the relative amplitude of the F = 1 Larmor signal is responsible for the saturation of the SP-SFF accuracy in Fig. 2(a).

These simulations assume a magnetic field of 50 µ T in the ˆ direction with pumping rates ranging from 0.3 to x 20 Ms -1 and a mean photon spin s = 0 9ˆ. . z By orienting a magnetic field perpendicular to the pump-probe laser axis, we minimize inaccuracies from the NLZ effect, allowing the hyperfine systematic error to dominate because the perpendicular field causes the Zeeman resonances to be symmetric, effectively canceling out the NLZ effect. For single-pulse pumping, we use T p = 10 µ s and for synchronous-pulse pumping, we use 200 pulses, each 30 ns long, square-wave modulated between 0 and R OP at the average F = 2 Larmor frequency f L, 2 (Fig. 1(c)). After pumping, we simulate free precession for T free = 1 ms (Fig. 1(d)).

## A. SP-SFF Method

The simulated accuracies for each magnetic-field extraction method are shown in Fig. 2(a). The SPSFF exhibits the largest inaccuracy as the pumping rate decreases. This method exhibits a distinctive pattern consistent with experimental findings (Fig. 5). At lower pumping rates, the single-frequency fitting method demonstrates a dip in accuracy, with a deviation of ap-

However, the relative amplitude alone does not explain the behavior of the SP-SFF accuracy at low pumping rates. A complete understanding necessitates considering the relative phase ϕ rel between the F = 1 and F = 2 Larmor precession signals. Because the atoms in the F = 1 and F = 2 manifolds precess in opposite directions, it is commonly assumed that the relative phase is π . Using an analytical model described in Appendix A, we can extract the relative phase from the density matrices used in the simulations of Fig. 2(a). The phases for each Larmor transition within a manifold are approx-

imately equal. However, the relative phase between different manifolds varies depending on the initial density matrix after pumping, which is influenced by the pumping rate and method (Fig. 2(c)). This relative phase arises from the steady-state orientation of the average F = I ± 1 / 2 atomic spins after optical pumping, a concept further explored in Sec. V. As the pumping rate goes to zero, ϕ rel → 0 and as the pumping rates goes to infinity, ϕ rel → π .

The insets of Fig. 2(c) depict the average F = 1 (blue) and F = 2 (red) Faraday rotation signals after optical pumping. These Faraday rotation signals were calculated by splitting the Faraday rotation operator into 3 × 3 and 5 × 5 sub-matrices for each hyperfine manifold and computing their expectation values separately. From these insets, we see that at low pumping rates ( R OP = 0.3 Ms -1 ), the F = I ± 1 / 2 Faraday rotation signals are approximately in phase. Conversely, at high pumping rates ( R OP = 20 Ms -1 ), these signals approach a phase difference of nearly 0 8 . π .

The shifting relative phase alters the interference between the average F = 1 and F = 2 Larmor precession signals, thereby affecting the accuracy of single-frequency extraction. Further detailed in Appendix B, this varying relative phase is responsible for the observed dip in the SP-SFF accuracy at low pumping rates (Fig. 2(a)). It is worth mentioning that two orthogonal probe beams can observe this relative phase and thereby distinguish between the F = 1 and F = 2 signals [10].

## B. SP-DFF Method

In regimes of low spin polarization, i.e. low pumping power and/or small mean photon spin, the systematic inaccuracy associated with the relative amplitude and phase between the average F = 1 and F = 2 precession can be mitigated by using the double-frequency fitting model described in Sec. III. The accuracy associated with SP-DFF is shown by the green dash-dotted line in Fig. 2(a). The simulated accuracy remains under about 0.5 nT within the investigated pumping rates. Comparing the worst accuracies, the SP-DFF method has an accuracy enhancement of roughly 7 (3.5 nT/0.5 nT) compared to the SP-SFF method. In contrast to singlefrequency fitting, the double-frequency fitting model provides a more accurate estimation of the magnetic field at low spin polarizations by accounting for the relative amplitude and phase of both F = 1 and F = 2 Larmor precession signals. Since the frequency separation between adjacent Zeeman transitions (36 Hz) is much less than the frequency separation between the average Zeeman transitions for the two hyperfine manifolds (1.4 kHz), extending the model to include all six Larmor frequencies present in a 87 Rb FID signal does not significantly enhance the accuracy.

## C. SY-SFF Method

The hyperfine systematic error can also be mitigated using synchronous-pulse pumping. The purple dashed line in Fig. 2(a) displays the accuracy associated with SY-SFF. Since synchronous-pulse pumping is a resonant process, high spin polarizations can be achieved with relatively low pumping rates compared to single-pulse pumping (Fig. 2(b)). Additionally, the relative phase, ϕ rel , approaches π at lower pumping rates compared to single-pulse pumping (Fig. 2(c)). This resonant process minimizes the F = 1 contribution to the Larmor precession signal, making single-frequency extraction methods more accurate. The simulated accuracy using SY-SFF is under 0.5 nT for the pumping rates investigated here. Similar to SP-DFF, this method also has an accuracy enhancement of roughly 7 (3.5 nT/0.5 nT) relative to SPSFF when comparing the most inaccurate data points. The peak of ϕ rel at 0.7 Ms -1 in Fig. 2(c) is not well understood and warrants further investigation.

Despite quickly saturating to the highest possible spin polarization of P = 0 9 (limited here by the mean pho-. ton spin | s | = 0 9), . the magnetic field extracted using the SY-SFF technique has a bias of approximately 0.3 nT. Similarly, when the SP-SFF technique also reaches the maximum possible spin polarization P = 0 9, the ex-. tracted field error is approximately 0.3 nT. This is due to imperfect spin polarization P &lt; 1 and inaccuracies in the average gyromagnetic ratio. If we use a perfect spin polarization of P = 1, then the inaccuracies due to the hyperfine systematic error are suppressed to the picotesla level. In this work, we define the gyromagnetic ratio as γ = f L, 2 /B = 6.99579 Hz/nT, where f L, 2 is the average F = 2 Larmor frequency. Furthermore, the inset of Fig. 2(a) shows that the accuracy at the same spin polarization is different for the SY-SFF and SP-SFF techniques. This occurs because different initial distributions can result in the same initial atomic spin polarization, and these distributions will have nonidentical amplitudes for each Larmor frequency component, leading to differences in accuracy.

Both the SP-DFF and the SY-SFF techniques successfully mitigate the F = 1 systematic error, with their respective accuracies consistently agreeing within 1 nT of each other throughout all the explored pumping rates. The advantage of the SP-DFF method is that it does not require prior knowledge of the magnetic field, whereas the SY-SFF method necessitates either prior field information or the implementation of a feedback loop to maintain resonance between the synchronous pulses and the Larmor precession. The SY-SFF method's advantage lies in its highly efficient pumping which enhances the accuracy of single-frequency extraction techniques such as demodulation or frequency counting. Although not explored in this study, 'enhanced optical pumping' (EOP) can generate high spin polarizations with relatively low optical power by applying a strong bias field (on the order of millitesla) along the laser axis during optical pumping

FIG. 3. Expectation values for (a) electron, (b) nuclear, and (c) total atomic spins during optical pumping. The red (blue) vectors represent the F = 2 ( F = 1) spin component. This simulation was performed with R OP = 5 Ms -1 at a 50 µ T magnetic field. The black vector indicates the laser axis. The horizontal axis represents the spin component along the y axis, while the vertical axis represents the spin component along the z axis. The darker vector represents the spin after 0 2 . µ s of optical pumping and the lighter vector represents the spin after 1 4 . µ s of optical pumping.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00898v2-figure-9c79f173-1c2e-4236-ada5-5ac9e8b17747.jpg)

[20]. Like the SY-SFF method, EOP also reduces the hyperfine systematic error since the F = 1 contributions are negligible.

It is important to note that these methods do not mitigate NLZ shifts within a hyperfine manifold. Therefore, they should be paired with another technique, such as the analytical correction outlined in [10], to improve magnetometer accuracy across all field directions and spin polarizations.

manifold, they are antiparallel [49]. The pumping light interacts strongly with the valence electron but not directly with the nucleus, because it is shielded by inner electron shells. However, the nucleus indirectly interacts with the pumping light via the hyperfine coupling with the electron.

## V. RELATIVE PHASE EFFECT

Grasping the factors that influence the relative phase is crucial, as it affects the accuracy of an FID magnetometer. The relative phase results from the initial orientation of the average F = I ± 1 / 2 spins after pumping. Therefore, understanding the relative phase necessitates understanding optical pumping. In this section we focus on examining the relative phase after single-pulse pumping and then extrapolate these insights to the case of synchronous-pulse pumping.

Since the 50 µ T magnetic field is oriented in the ˆ dix rection, the spins precess in the yz plane (Fig. 1). As before, we use a mean photon spin of s = 0.9ˆ, meanz ing that the pumping axis is oriented in the ˆ direction z resulting in σ + transitions. We calculate the electron, nuclear, and total atomic spins in the | F, m F ⟩ z basis for both the F = 1 and F = 2 manifolds by splitting the spin operators into 3 × 3 and 5 × 5 submatrices and calculating their expectation values. Figures 3(a), 3(b), and 3(c) depict the electron, nuclear, and total atomic spins, respectively, during single-pulse pumping.

Due to opposite signs of the F = 1 and F = 2 g -factors, the respective spins rotate in opposite directions. Furthermore, in the F = 2 manifold, the average electron and nuclear spins are parallel, while in the F = 1

Since the electron strongly interacts with the pumping light, the electron spin is initially pumped parallel to the pumping axis (+ˆ) for both hyperfine manifolds, as obz served at 0 2 . µ s after the pumping light is switched on in Fig. 3(a). In the F = 2 manifold, the nuclear spin is also initially pumped parallel to the pumping axis because it must be parallel to the electron spin. Conversely, in the F = 1 manifold, the nuclear spin is initially pumped antiparallel to the pumping axis because it must be antiparallel to the electron spin [0 2 . µ s vector in Fig. 3(b)]. Because the nucleus stores most of the angular momentum ( I = 3 2 versus / S = 1 2), the average total atomic / spin mirrors the behavior of the average nuclear spin. Consequently, the F = 1 total atomic spin is initially pumped antiparallel to the pumping axis, polarizing the F = 1 manifold towards the | 1 , - ⟩ 1 z state (Fig. 3(c)).

Eventually, the atoms reach a steady state between the precession caused by the magnetic field and the optical pumping. For the 5 Ms -1 pumping rate used in these simulations, the spins reach their steady state after approximately 1 4 . µ s of pumping. Interestingly, repopulation pumping causes the average F = 1 spin to rotate more than the average F = 2 spin (Fig. 3). Initially pumped antiparallel to the pumping axis, into the | 1 , - ⟩ 1 z state, repopulation rotates the average F = 1 spin towards the pumping axis, moving population into the | 1 1 , ⟩ z state.

In the presence of buffer gas, collisional mixing quickly depolarizes electron spins in the excited state, leading to an equal probability of decaying into either groundstate hyperfine manifold [49]. However, these collisions

FIG. 4. Simulated total atomic spins after single-pulse pumping for 10 µ s with σ + photons at various pumping rates with a 50 µ T field oriented in the ˆ direction. x The blue vectors are the total atomic spin for the F = 1 manifold and the red vectors are the total atomic spin for the F = 2 manifold. The black arrow represents the pumping laser axis. The green vectors indicate the F = 1 atomic spin as measured by the probe laser. Here the z component of the F = 1 atomic spin is multiplied by -1 to account for the Faraday rotation as described by Eq. (4). The angle between the F = 2 vector (red) and the Faraday rotated F = 1 vector (green) represents the relative phase, ϕ rel .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00898v2-figure-7647f307-f598-4dfa-9ab9-1ae54a4694bc.jpg)

are sudden with respect to the nuclear spins such that the atoms preserve their nuclear spin upon deexcitation [45]. Therefore, the atoms do not have equal probability of decaying into each hyperfine sublevel but are restricted to the sublevels that preserve the nuclear angular momentum. After an interval of optical pumping, a majority of the population lies in the F = 2 manifold, biased towards the positive m F sublevels which carry positive nuclear angular momentum. After these atoms are excited, half of them decay into the F = 1 manifold, mainly repopulating the | 1 1 , ⟩ z and | 1 0 , ⟩ z states, flipping the direction of the average F = 1 spin. Combined with the precession caused by the magnetic field, the repopulation effect causes the F = 1 spin vectors to rotate more during optical pumping than the F = 2 spin vectors such that the F = 1 and F = 2 spins reach a steady state nearly parallel to each other (Fig. 3(c)).

As shown in Fig. 4, a higher pumping rate results in the steady-state total atomic spins approaching closer to alignment with the pumping axis. The probe laser, aligned with the ˆ direction, z exclusively detects the z component of atomic spins. Since the Faraday rotation operator is proportional to -m F in the F = 1 manifold [Eq. (4)], we illustrate the effect of Faraday rotation by multiplying the z component of the F = 1 atomic spin by -1 (green vector). The angle between the F = 2 atomic spin vector (red) and the F = 1 atomic spin vec- tor as seen by the probe laser (green) matches the relative phase extracted from the density matrix.

In the infinite pumping rate limit, the steady-state spins perfectly align with the pumping axis (Fig. 4). Consequently, the Faraday rotation has a phase shift of π between the F = 1 and F = 2 Larmor signals. Conversely, at very low pumping rates, the spins take longer to reach a steady state, potentially exceeding the duration of the pumping pulse. Given the comparable magnitude of the F = 1 and F = 2 g -factors, the spins will oscillate in opposite directions with the average F = 2 spin starting parallel and the average F = 1 spin starting antiparallel to the pumping axis such that the relative phase of the respective Larmor signals is always approximately zero. In between these limits, the average F = 1 and F = 2 spins reach a steady state such that the relative phase lies within 0 ≤ ϕ rel ≤ π .

The physics is the same for synchronous-pulse pumping schemes. The main difference is that the spins are allowed to precess an integer number of Larmor periods between pump pulses. Since synchronous-pulse pumping is more efficient than single-pulse pumping, the initial atomic spins are oriented closer to the pumping axis [35, 43]. Therefore, the relative phase converges to π at much lower pumping rates compared to single-pulse pumping (Fig. 2(c)).

## VI. FITTING EXPERIMENTAL DATA

To confirm the theory laid out in the previous sections, we measured FID signals at various optical pumping powers from 28 to 465 mW while keeping the probe power fixed at 5 mW. The pumping laser frequency was resonant with the D 1 transition, while the probe laser frequency was approximately -80 GHz detuned from the D 2 transition. A simplified diagram of our experimental setup is shown in Fig. 1(b). The measurements were conducted in a magnetic field of approximately 20 4 . µ T, oriented in the ˆ direction. x We measured 100 FID signals at each pumping power for both single-pulse pumping and synchronous-pulse pumping. In the single-pulse pumping scheme, we used T p = 10 µ s, while the synchronous-pulse pumping scheme used 200 30-ns pulses modulated at the average Larmor frequency.

Figure 5(a) displays the mean magnetic field extracted from the FID signals across a range of pumping powers using the SP-SFF, SP-DFF, and SY-SFF methods. Figure 5(b) displays the simulated accuracy with the same parameters as before but at approximately the same magnetic field as the experimental data. The extracted magnetic field associated with the SP-SFF method (blue solid line) has the same distinctive pattern observed in the simulated data of Fig. 5(b). The simulations revealed that the dip in the extracted field at low pumping powers can be attributed to the relative phase between the F = 1 and F = 2 Larmor precession signals, while the saturation of the extracted magnetic field at high pumping

FIG. 5. (a) Experimentally measured pump-power dependence of the extracted magnetic field. The signals generated via single-pulse pumping data are fit to the single-frequency model (blue solid line with squares) and the double-frequency model (green dash-dotted line with circles). The synchronously pumped data are fit to the single-frequency model (purple dashed line with circles). The error bars are the standard deviations of the magnetic field outputted from the fits. The inset of this graph plots the average relative phase between the F = 2 and F = 1 manifolds extracted from the double-frequency fits. (b) Simulated magnetic-field accuracy at a bias field of 20 435 . µ T. The blue solid line is the SP-SFF accuracy, the green dash-dotted line is the SP-DFF accuracy, and the purple dashed line is the SY-SFF accuracy.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00898v2-figure-69018fae-5c98-41bd-a08c-dcbf8ecb0785.jpg)

powers can be attributed to the diminishing amplitude of the F = 1 signal. Indicated by the blue bar on the right side of Fig. 5(a), the extracted magnetic field using SP-SFF spans approximately 3 nT, aligning closely with our simulation.

to the SP-SFF method.

By fitting the same single-pulse pumped data with the double-frequency fitting model, we can minimize the systematic error caused by the F = 1 contributions to the FID signal, as demonstrated by the green dash-dotted line in Fig. 5. Utilizing the SP-DFF technique, the extracted magnetic field demonstrates a diminished dependence on pumping power, spanning about 0.9 nT, highlighted by the green bar on the graph's right side. Comparing the ranges of the extracted field, this is an enhancement of approximately 3.3 (3 nT/0.9 nT) compared to the SP-SFF method. Just as in the simulation, this is attributed to the double-frequency fitting model accounting for the relative amplitude and phase of the F = 1 Larmor signal. Moreover, the relative phase extracted from the fit parameters of the double-frequency fits, depicted in the inset of Fig. 5, mirrors the relative phase behavior observed in the simulations shown in Fig. 2. At low pumping powers, the relative phase approaches zero and as the pumping power increases the relative phase increases to π .

The purple dashed line in Fig. 5 shows the extracted magnetic field using the SY-SFF method. Synchronouspulse pumping efficiently drives the atomic ensemble into the stretched state, minimizing F = 1 contributions even at lower pumping rates (Fig. 2(b)). This reduces the dependence on pumping power. The magnetic field extracted using the SY-SFF technique spans 0.8 nT, as indicated by the purple bar on the right side of the graph, which is an enhancement of 3.75 (3 nT/0.8 nT) relative

Significantly, both SP-DFF and SY-SFF magneticfield-extraction techniques exhibit agreement, demonstrating discrepancies of less than 1.5 nT between each other at all investigated pumping powers. This agreement is apparent when inspecting the overlapping range bars for both data sets. This, coupled with the outcomes of our simulations, strongly suggests that the magnetic field extracted through the SP-DFF and SY-SFF techniques offers superior accuracy across an expanded range of spin polarizations compared to the magnetic field extracted utilizing the SP-SFF technique. Notably, the extracted magnetic fields from both the SP-DFF and SYSFF techniques deviate from each other at the lowest pumping powers. This might be explained by a decreasing signal-to-noise ratio or by field drifts between measurements which were not included in the simulations. While we cannot assert absolute accuracy for the magnetometer as there may be other systematic effects like light shifts, nonlinear Zeeman shifts, and magnetic-field gradients, these improved magnetic-field-extraction techniques have effectively mitigated inaccuracies caused by the hyperfine systematic error to the ± 1 nT level.

## VII. CONCLUSION

We studied the accuracy of extracting the magnetic field from FID signals over a wide range of spin polarizations. As the spin polarization changes, the relative amplitudes and phases between the F = 1 and F = 2 Zeeman frequency components change, ultimately affecting the accuracy of magnetic-field extraction. Simulations reveal that F = 1 contributions introduce inaccu-

racies of approximately 3.5 nT when employing the SPSFF method. Theoretically, these inaccuracies can be reduced to under 0.5 nT through the utilization of either SP-DFF or SY-SFF techniques. Measurements of the extracted magnetic field closely match the simulations, displaying the same characteristic patterns. The variation of the magnetic field extracted from the SP-SFF measurements is approximately 3 nT, while the variations of the magnetic field extracted from the SP-DFF and SY-SFF methods are both under 1 nT. Furthermore, the magnetic fields measured using the SP-DFF and SY-SFF techniques agree to within 1.5 nT at all pumping powers investigated.

The SP-DFF technique is well suited for accurately extracting the magnetic field from FID signals in the lowspin-polarization regime P ≤ 0 5. . Furthermore, it does not require feedback loops or prior knowledge of the magnetic field to achieve high accuracies. In contrast, the SY-SFF scheme generates higher spin polarizations at lower pumping rates, making single-frequency extraction techniques such as demodulation, instantaneous-phase retrieval, or frequency counters more accurate. However, this enhancement necessitates prior knowledge of the magnetic field to guess the pulse-modulation frequency. While the techniques introduced in this paper do not directly address absolute accuracy or the NLZ effect, they effectively mitigate inaccuracies linked to the overlap of F = 1 and F = 2 Zeeman resonances. Compact sensors utilizing microfabricated cells with high buffer-gas pressure, operating at elevated temperatures and technically constrained by moderate optical power levels, are particularly susceptible to this systematic inaccuracy. The techniques described in this paper can significantly improve the accuracy of such sensors under limited optical power conditions.

## ACKNOWLEDGMENTS

We wish to acknowledge helpful discussions with Ricardo Jim´ enez-Mart´ ınez, Georg Bison, and Tobias Thiele. This work was supported by the DARPA SAVaNT program through grant W911NF-21-1-0127, and NSF QLCI Award OMA-2016244.

## Appendix A: Analytical Model for Extracting Phases

We can write an analytical model that does not include decoherence to help us understand the physics of the system and approximate parameters such as the phase of each Larmor transition. We start by diagonalizing the Hamiltonian from Eq. (4), ˜ H → ˜ H d . Next we convert the Faraday rotation operator and the input density matrix into the basis that diagonalizes the Hamiltonian, θ F → θ F,d = U θ † F U and ρ 0 → ρ 0 ,d = U ρ U † 0 , where U is the unitary matrix of eigenvectors that di-

FIG. 6. Faraday rotation signals without decoherence simulated using the full density-matrix model (red solid line) and the analytical model (black dashed line). The consistency between the analytical model and the density matrix highlights the analytical model's ability to accurately predict the phases of each Larmor transition.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00898v2-figure-771686b3-a554-45b5-89c4-1b43ae783ca5.jpg)

agonalizes the Hamiltonian. The Schr¨dinger o equation can be solved without decoherence by ρ 0 ,d (t) = e -iH t/ d ℏ ρ 0 ,d iH t/ d ℏ . The Faraday rotation can be calculated using θ F ∝ Tr[ θ F,d ρ 0 ,d ( )]. t Finally, we find that the Faraday rotation can be calculated at any time based on the input density matrix using the equation

̸

$$\text{tical} & & \theta _ { F } \infty & \sum _ { i = 1 } ^ { 8 } F _ { d } ^ { i i } \rho _ { 0, d } ^ { i i } \\ & & + \sum _ { i \neq j = 1 } ^ { 8 } 2 | F _ { d } ^ { i j } \rho _ { 0, d } ^ { j i } | \cos [ 2 \pi ( f _ { i } - f _ { j } ) t + \arg ( F _ { d } ^ { i j } \rho _ { 0, d } ^ { j i } ) ]. \\ \text{with} & & \text{has} & & \text{The first term accounts for the dc offset, while the second}$$

The first term accounts for the dc offset, while the second term represents the precession. The argument of the cosine consists of two terms: The first term is the frequency and the second term is the phase. Figure 6 shows the excellent agreement between the analytical model and density-matrix simulations. For this density-matrix simulation, we included decoherence during pumping, but not during free precession.

The agreement observed here demonstrates the analytical model's ability to accurately predict the initial phases of the simulated FID signals, even when decoherence effects are included during pumping. Furthermore, the phase of each transition, calculated using arg( F ij d ρ ji 0 ,d ), are approximately equal within a hyperfine manifold. Consequently, we can precisely determine relative phases by computing the difference between the absolute values of the average F = 1 and F = 2 phases.

Appendix B: Separating the Effects of Relative Amplitudes and Phases on the Accuracy

FIG. 7. Accuracy of single-frequency magnetic field extraction versus relative phase, ϕ rel , at various magnetic fields. The F = 2 and F = 1 amplitudes and time constants are fixed while the F = 2 and F = 1 phases are varied from π/ 2 to 0 and π/ 2 to π , respectively, in 100 evenly spaced steps.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00898v2-figure-c8519e25-37f9-496c-9c97-8c3832ce8306.jpg)

To show that the interference caused by the relative phase between the hyperfine manifolds contributes to the fitting inaccuracies, we fit fake signals generated using the double-frequency fitting equation to the singlefrequency fitting model. Using the Breit-Rabi formula, the magnetic field constrains the frequencies of each manifold. The dc offset is set to 0. To approximate an FID signal, the F = 2 and F = 1 amplitudes are fixed to 1.5 and 0.25, respectively, while the F = 2 and F = 1 time constants are fixed to 1 and 0.25 ms, respectively. We vary the F = 2 phase from π/ 2 to 0 and the F = 1 phase from π/ 2 to π in 100 evenly spaced steps. Figure 7 displays the magnetic-field-extraction accuracy using the single-frequency fitting equation for various magnetic fields. This graph shows that the accuracy of fitting a sig-

- [1] C. Abel, S. Afach, N. J. Ayres, C. A. Baker, G. Ban, G. Bison, K. Bodek, V. Bondar, M. Burghoff, and E. e. a. Chanel, Measurement of the permanent electric dipole moment of the neutron, Phys. Rev. Lett. 124 , 081803 (2020).
- [2] C. Abel, N. J. Ayres, G. Ban, G. Bison, K. Bodek, V. Bondar, M. Daum, M. Fairbairn, V. V. Flambaum, and P. e. a. Geltenbort, Search for axionlike dark matter through nuclear spin precession in electric and magnetic fields, Phys. Rev. X 7 , 041034 (2017).
- [3] S. Groeger, G. Bison, P. E. Knowles, R. Wynands, and A. Weis, Laser-pumped cesium magnetometers for high-

nal made up of two frequencies varies with the relative phase.

FIG. 8. Magnetic-field accuracy of the single-frequency fit as the amplitudes are varied. The horizontal axis of this plot is the F = 2 amplitude A F =2 minus the F = 1 amplitude A F =1 . Here the time constants and phases are fixed, while the F = 2 and F = 1 amplitudes are varied from 0.25 to 2 and from 0.15 to 0, respectively, in 100 evenly spaced steps.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00898v2-figure-6d824393-007d-4ef8-8aa7-349ad3c6c7f3.jpg)

Since we keep the amplitudes fixed in the rudimentary simulation of Fig. 7, we do not see the saturation of the single-frequency magnetic field accuracy in the high spin-polarization limit as described in Sec. IV. We can separate the saturation from the interference by fixing the F = 2 and F = 1 time constants to 1 and 0.25 ms and the phases to 0 and π , respectively, while varying the F = 2 amplitude A F =2 from 0.25 to 2, and the F = 1 amplitude A F =1 from 0.15 to 0 in 100 evenly spaced steps. We see from Fig. 8 that as the F = 2 population starts to dominate the signal, the magnetic field extracted from the single-frequency fit saturates to the true value. These two rudimentary simulations confirm that the accuracy saturation is linked to the relative amplitudes while the dip is linked to the relative phase between hyperfine manifolds.

resolution medical and fundamental research, Sensors and Actuators A: Physical 129 , 1 (2006), EMSA 2004.

- [4] H. Korth, K. Strohbehn, F. Tejada, A. G. Andreou, J. Kitching, S. Knappe, S. J. Lehtonen, S. M. London, and M. Kafel, Miniature atomic scalar magnetometer for space based on the Rubidium isotope 87 Rb, Journal of Geophysical Research: Space Physics 121 , 7870 (2016).
- [5] J. S. Bennett, B. E. Vyhnalek, H. Greenall, E. M. Bridge, F. Gotardo, S. Forstner, G. I. Harris, F. A. Miranda, and W. P. Bowen, Precision magnetometers for aerospace applications: A review, Sensors 21 , 5568 (2021).
- [6] T. M. Tierney, N. Holmes, S. Mellor, J. D. L´ opez,

- G. Roberts, R. M. Hill, E. Boto, J. Leggett, V. Shah, and M. J. B. et. al., Optically pumped magnetometers: From quantum origins to multi-channel magnetoencephalography, NeuroImage 199 , 598 (2019).
- [7] E. Boto, S. S. Meyer, V. Shah, O. Alem, S. Knappe, P. Kruger, T. M. Fromhold, M. Lim, P. M. Glover, and P. G. M. et. al., A new generation of magnetoencephalography: Room temperature measurements using opticallypumped magnetometers, NeuroImage 149 , 404 (2017).
- [8] G. Bison, R. Wynands, and A. Weis, Dynamical mapping of the human cardiomagnetic field with a roomtemperature, laser-optical sensor, Opt. Express 11 , 904 (2003).
- [9] T. H. Sander, J. Preusser, R. Mhaskar, J. Kitching, L. Trahms, and S. Knappe, Magnetoencephalography with a chip-scale atomic magnetometer, Biomed. Opt. Express 3 , 981 (2012).
- [10] W. Lee, V. G. Lucivero, M. V. Romalis, M. E. Limes, E. L. Foley, and T. W. Kornack, Heading errors in alloptical alkali-metal-vapor magnetometers in geomagnetic fields, Phys. Rev. A 103 , 063103 (2021).
- [11] D. Steck, Rubidium 87 D line data, available online at http://steck.us/alkalidata (2021).
- [12] Y. Rosenzweig, D. Tokar, I. Shcerback, M. Givon, and R. Folman, Heading error compensation in a portable optical magnetometer using a double-pass single beam configuration (2023), arXiv:2307.13982 [physics.app-ph].
- [13] G. Oelsner, V. Schultze, R. IJsselsteijn, F. Wittk¨mper, a and R. Stolz, Sources of heading errors in optically pumped magnetometers operated in the Earth's magnetic field, Phys. Rev. A 99 , 013420 (2019).
- [14] S. J. Seltzer, P. J. Meares, and M. V. Romalis, Synchronous optical pumping of quantum revival beats for atomic magnetometry, Phys. Rev. A 75 , 051407(R) (2007).
- [15] G. Bao, A. Wickenbrock, S. Rochester, W. Zhang, and D. Budker, Suppression of the nonlinear Zeeman effect and heading error in Earth-field-range alkali-vapor magnetometers, Phys. Rev. Lett. 120 , 033202 (2018).
- [16] R. Zhang, D. Kanta, A. Wickenbrock, H. Guo, and D. Budker, Heading-error-free optical atomic magnetometry in the Earth-field range, Phys. Rev. Lett. 130 , 153601 (2023).
- [17] V. M. Acosta, M. Auzinsh, W. Gawlik, P. Grisins, J. M. Higbie, D. F. J. Kimball, L. Krzemien, M. P. Ledbetter, S. Pustelny, and S. M. R. et. al., Production and detection of atomic hexadecapole at Earth's magnetic field, Opt. Express 16 , 11423 (2008).
- [18] A. Nagel, Precision Spectroscopy of Coherent Dark States in Thermal Cesium Vapour , Ph.D. thesis, University of Bonn (1999).
- [19] L. M. Rushton, L. Elson, A. Meraki, and K. Jensen, Alignment-based optically pumped magnetometer using a buffer-gas cell, Physical Review Applied 19 , 064047 (2023).
- [20] D. Hunter, M. S. Mrozowski, A. McWilliam, S. J. Ingleby, T. E. Dyer, P. F. Griffin, and E. Riis, Optical pumping enhancement of a free-induction-decay magnetometer, JOSA B 40 , 2664 (2023).
- [21] A. Canciani and J. Raquet, Absolute positioning using the Earth's magnetic anomaly field, Navigation 63 , 111 (2016).
- [22] Z. D. Gruji´, P. A. Koss, G. Bison, and A. Weis, A senc sitive and accurate atomic magnetometer based on free
- spin precession, The European Physical Journal D 69 , 10.1140/epjd/e2015-50875-3 (2015).
- [23] I. Fratter, J.-M. L´ger, F. Bertrand, T. Jager, G. Hulot, e L. Brocco, and P. Vigneron, Swarm absolute scalar magnetometers first in-orbit results, Acta Astronautica 121 , 76 (2016).
- [24] W. E. Bell and A. L. Bloom, Optically driven spin precession, Phys. Rev. Lett. 6 , 280 (1961).
- [25] V. Gerginov, S. Krzyzewski, and S. Knappe, Pulsed operation of a miniature scalar optically pumped magnetometer, J. Opt. Soc. Am. B 34 , 1429 (2017).
- [26] P. D. D. Schwindt, B. Lindseth, S. Knappe, V. Shah, J. Kitching, and L.-A. Liew, Chip-scale atomic magnetometer with improved sensitivity by use of the M x technique, Applied Physics Letters 90 , 10.1063/1.2709532 (2007), 081102.
- [27] M. Ellmeier, C. Amtmann, A. Pollinger, W. Magnes, C. Hagen, A. Betzler, I. Jernej, M. Ag´, L. Windholz, and u R. Lammegger, Frequency shift compensation for single and dual laser beam pass sensors of a coherent population trapping resonance based coupled dark state magnetometer, Measurement: Sensors 25 , 100606 (2023).
- [28] M. St¨hler, R. Wynands, S. Knappe, J. Kitching, L. Holla berg, A. Taichenachev, and V. Yudin, Coherent population trapping resonances in thermal 85 Rb vapor: D 1 versus D 2 line excitation, Opt. Lett. 27 , 1472 (2002).
- [29] P. D. D. Schwindt, S. Knappe, V. Shah, L. Hollberg, J. Kitching, L.-A. Liew, and J. Moreland, Chip-scale atomic magnetometer, Applied Physics Letters 85 , 6409 (2004).
- [30] T. Thiele, Y. Lin, M. O. Brown, and C. A. Regal, Self-calibrating vector atomic magnetometry through microwave polarization reconstruction, Phys. Rev. Lett. 121 , 153202 (2018).
- [31] C. Kiehl, D. Wagner, T.-W. Hsu, S. Knappe, C. A. Regal, and T. Thiele, Coherence of Rabi oscillations with spin exchange, Phys. Rev. Res. 5 , L012002 (2023).
- [32] C. Kiehl, T. S. Menon, D. P. Hewatt, S. Knappe, T. Thiele, and C. A. Regal, Correcting heading errors in optically pumped magnetometers through microwave interrogation, Phys. Rev. Appl. 22 , 014005 (2024).
- [33] J. C. Allred, R. N. Lyman, T. W. Kornack, and M. V. Romalis, High-sensitivity atomic magnetometer unaffected by spin-exchange relaxation, Phys. Rev. Lett. 89 , 130801 (2002).
- [34] M. P. Ledbetter, I. M. Savukov, V. M. Acosta, D. Budker, and M. V. Romalis, Spin-exchange-relaxation-free magnetometry with Cs vapor, Phys. Rev. A 77 , 033408 (2008).
- [35] D. Hunter, S. Piccolomo, J. D. Pritchard, N. L. Brockie, T. E. Dyer, and E. Riis, Free-induction-decay magnetometer based on a microfabricated Cs vapor cell, Phys. Rev. Appl. 10 , 014002 (2018).
- [36] S. Afach, G. Ban, G. Bison, K. Bodek, Z. Chowdhuri, Z. D. Gruji´, c L. Hayen, V. H´ elaine, M. Kasprzak, and K. K. et. al., Highly stable atomic vector magnetometer based on free spin precession, Opt. Express 23 , 22108 (2015).
- [37] A. Fabricant, I. Novikova, and G. Bison, How to build a magnetometer with thermal atomic vapor: A tutorial, New Journal of Physics 25 , 025001 (2023).
- [38] A. Pollinger, R. Lammegger, W. Magnes, C. Hagen, M. Ellmeier, I. Jernej, M. Leichtfried, C. K¨ urbisch, R. Maierhofer, and R. W. et. al., Coupled dark state mag-

- netometer for the China Seismo-Electromagnetic Satellite, Measurement Science and Technology 29 , 095103 (2018).
- [39] S.-Q. Liang, G.-Q. Yang, Y.-F. Xu, Q. Lin, Z.-H. Liu, and Z.-X. Chen, Simultaneously improving the sensitivity and absolute accuracy of CPT magnetometer, Opt. Express 22 , 6837 (2014).
- [40] I. K. Kominis, T. W. Kornack, J. C. Allred, and M. V. Romalis, A subfemtotesla multichannel atomic magnetometer, Nature 422 , 596 (2003).
- [41] H. B. Dang, A. C. Maloof, and M. V. Romalis, Ultrahigh sensitivity magnetic field and magnetization measurements with an atomic magnetometer, Applied Physics Letters 97 , 151110 (2010).
- [42] D. Sheng, S. Li, N. Dural, and M. V. Romalis, Subfemtotesla scalar atomic magnetometry using multipass cells, Phys. Rev. Lett. 110 , 160802 (2013).
- [43] V. Gerginov, M. Pomponio, and S. Knappe, Scalar magnetometry below 100 fT/Hz 1 / 2 in a microfabricated cell, IEEE Sensors Journal 20 , 12684 (2020).
- [44] G. Breit and I. I. Rabi, Measurement of nuclear spin, Phys. Rev. 38 , 2082 (1931).
- [45] S. Appelt, A. B.-A. Baranga, C. J. Erickson, M. V. Romalis, A. R. Young, and W. Happer, Theory of spinexchange optical pumping of 3 He and 129 Xe, Phys. Rev. A 58 , 1412 (1998).
- [46] Y.-Y. Jau, New studies of optical pumping, spin resonances, and spin exchange in mixtures of inert gases and alkali-metal vapors , Ph.D. thesis, Princeton University (2005).
- [47] N. Wilson, C. Perrella, R. Anderson, A. Luiten, and P. Light, Wide-bandwidth atomic magnetometry via instantaneous-phase retrieval, Physical Review Research 2 , 013213 (2020).
- [48] D. Hunter, T. E. Dyer, and E. Riis, Accurate optically pumped magnetometer based on Ramsey-style interrogation, Opt. Lett. 47 , 1230 (2022).
- [49] S. J. Seltzer, Developments in alkali-metal atomic magnetometry , Ph.D. thesis, Princeton University (2008).
- [50] M. Hern´ andez Ruiz, Y. Ma, H. Medhat, C. Mazzinghi, V. G. Lucivero, and M. W. Mitchell, Cavity-enhanced detection of spin polarization in a microfabricated atomic vapor cell, Phys. Rev. Appl. 21 , 064014 (2024).
D. P. Hewatt M. Ellmeier C. Kiehl T. S. Menon J. W. Pollock C. A. Regal S. Knappe
Investigating the hyperfine systematic error and relative phase in low spin-polarization alkali FID magnetometers
Alkali-metal optically-pumped magnetometers are prone to inaccuracies arising
from the overlap of the average F = I + 1/2 and F = I - 1/2 ground-state Zeeman
resonances. We employ density-matrix simulations and experiments to investigate
how this hyperfine systematic error varies with spin polarization in a
$^{87}$Rb free-induction-decay (FID) magnetometer. At low spin polarizations,
($P \leq 0.5$), this effect causes single-frequency magnetic-field extraction
techniques to exhibit inaccuracies up to approximately 3.5 nT. Density-matrix
simulations reveal that this bias can be traced to the relative amplitude and
phase between the F = I $\pm$ 1/2 hyperfine ground-state manifolds in the FID
spin precession signal. We show that this systematic error can be mitigated
using either a double-frequency fitting model that accounts for the relative
amplitude and phase or synchronous-pulse pumping, that minimizes the F = 1
contribution to the FID signal. Theoretical simulations predict accuracies
within 0.5 nT for both techniques across a wide range of spin polarizations,
suggesting a sevenfold enhancement over single-frequency extraction methods.
Our experiments validate this, showcasing a variation in the extracted field
below 1 nT, a 3.5-fold improvement compared to single-frequency extraction
methods. Furthermore, the mitigation techniques demonstrate agreement of the
extracted magnetic field within 1.5 nT.