circularly polarized, yielding atomic orientation, a configuration thoroughly studied in the 1970's [13–16].

In parallel, these last years, there is a renewed interest for optically pumped magnetometers based on linearly polarized light, yielding atomic alignment [17–21]. We recently worked on PRM based on linearly polarized light. For metastable helium-4 [22–24], this configuration is advantageous, showing better sensitivities in a multi-axis configuration, and no vector light-shifts [25].

Both for orientation and alignment-based PRM two approximations are usual for obtaining closed-form expressions for the photodetection signals:

* The three step approach consists in modelling the magnetometer dynamics as if it happened in three sequential steps: (i) state preparation by optical pumping, (ii) evolution under a magnetic field and relaxation, and (iii) measurement of the system state. This approach introduces a steady-state polarization resulting from the pumping. This assumption is accurate as far as the pumping and probing light have sufficiently low intensities [26–28].

* The rotating wave approximation (RWA) requires the RF frequencies to be much larger than other characteristic frequencies [16, 25], including the total spin relaxation rate and the precession rates corresponding to each component of the magnetic field e. g. $ω_x = γB_x$, with γ the gyromagnetic ratio of the atomic species.

The conditions for these approximations are hardly fulfilled in practical magnetometers. The first one because it is well known that the best sensitivities are obtained with pumping rates much larger than the relaxation rates in the dark [14]. The second one because the slower RF field frequency $Ω$ is usually just a few times larger than the relaxation rate. Indeed, it has to comply with $Ω \ll ω$, while the faster RF field frequency $ω$ remains low enough so that the ratio $γB_1/ω$ is near the optimum of sensitivity [16, 25]. For instance, with helium-4, choosing $ω = 2π 100 kHz$ yields an optimal $B_1$ of 3.5 μT, which is difficult to generate without increasing noise around DC.

It is important to have an exact picture of the consequences of the breakdown of these approximations. Indeed, optically pumped magnetometers are often used for applications requiring good accuracy in addition to ultra-low noise. Therefore the presence of offsets in the magnetic field readings, or the emergence of additional sources of noise, could seriously compromise the interest of OPMs. From the same perspective, pumping with linearly polarized light has the advantage of inducing no vector light-shifts, which could result in field offsets or increased noise [29]. However linearly-polarized light detuned from an optical transition causes tensor light-shifts (TLS) [30–32]. It is not clear if such TLS could also translate into systematic errors or increased noise for an alignment-based magnetometer operating in very low magnetic fields.

In this paper we aim to study these three second-order effects in order to clarify how they translate to metrological properties of optically pumped magnetometers based on atomic alignment.

We will start (in Sect. 2) by briefly recalling the Hanle effect and PRM schemes that we study. In Sect. 3 we show that the three-step approach results in inaccurate predictions for photodetection signals under high pump intensities, and how to refine it. In Sect. 4 we calculate the first correction terms to RWA and discuss on their impact on magnetic measurements. In Sect. 5 we analyze the impact of TLS on magnetic field measurements. Finally in Sect. 6 we discuss the results and conclude.
## 2 Overview of the schemes used for studying the Hanle effect and the parametric resonance magnetometer based on alignment

We will focus here on the metastable state of helium-4 ($F = 1$) as magnetically-sensitive species. The spin polarization of this state can be prepared by optical pumping towards the $F' = 0$ state $2^3P_0$ [33]. Both atomic orientation and alignment can be obtained in this way. We will focus here on alignment only, which has been shown to display interesting advantages [25]. Both these polarizations decay, notably due to the collisional relaxation of metastable state, at rates $\Gamma_e$ of a few kHz in 5–10 mm sized cells [34].

The setup for observation of the Hanle effect is represented in Fig. 1(a): a linearly polarized light beam is sent through a helium-4 cell along the z axis with the polarization electric field along the x axis. When the magnetic field transverse to x is not null, the alignment created by the pumping is partially suppressed and optical absorption increases [35]. The photodetection signals are even-symmetric in magnetic field, which is not convenient for measuring the magnetic field. However, like in the case of atomic orientation, dispersive (i.e. odd-symmetric) photodetection signals can be obtained by adding an independent probe beam [36].

Parametric resonance magnetometer scheme is shown in Fig. 1(b). Compared to the previous scheme, two RF fields are added along the z and y axes, both orthogonal to the linear polarization of the pump, with angular frequencies $\omega$ and $\Omega$ which comply with $\omega \gg \Omega \gg \Gamma$. The photodetection signal contains harmonics and inter-harmonics of $\omega$ and $\Omega$. The first harmonic of each RF is sensitive to the $B_0$ component parallel to the axis of the corresponding RF, and their first inter-harmonics ($\omega \pm \Omega$) are sensitive to the third component of the magnetic field with lower sensitivity. It is usual to operate PRM

<figure><img src="image_1.png" /><figcaption>**Figure 1** Null-field alignment-based magnetometry schemes. (a) Hanle effect on atomic alignment: a cell containing helium, with the metastable state populated by a high-frequency discharge, is subject to a linearly-polarized light-beam, with light electric field $\vec{E}_p$. This beam is photodetected after the cell. (b) Parametric resonance magnetometer based on atomic alignment: in the addition to the elements of the former scheme, this one includes two RF fields set along the axes orthogonal to $\vec{E}_p$. (c) Experimental setup which allows an implementation of both schemes</figcaption></figure>
in a closed-loop mode, where these three signals are used to cancel the measured field by applying opposite fields on compensation coils [37].

In both configurations, our experimental setup is composed of the elements shown in Fig. 1(c). The optical beam is generated by a fiber-coupled DFB laser diode (QD Laser QLD1061-8330) tuned to D0 line of helium-4 (at 1083.206 nm wavelength in the vacuum) with a linewidth <5 MHz and no wavelength stabilisation. Its intensity can be varied thanks to a variable fiber optical attenuator. The beam is collimated with a 5.7 mm $1/e^2$ waist by a 3D-printed collimator. Although we use polarization-maintaining fibers, a linear polarizer is placed before the cell to improve the polarization extinction ratio. This beam is sent through a pyrex cell of 12 mm length and 10 mm diameter containing ultra-pure helium-4. Due to practical issues of availability, the setup for observing Hanle effect includes an aluminium mirror at the end of the cell, and the light, with 5° incidence angle, undergoes a double pass in the helium cell. The setup for parametric resonance does not include such a mirror, and the light undergoes a single pass in the cell. Under appropriate high-frequency excitation, the density of helium atoms excited to the metastable state is on the order of $10^{11}$ cm$^{-3}$. The optical beam is photodetected by an InGaAs photodiode, connected to a home-made 25 kΩ transimpedance amplifier. Both DC and RF magnetic fields applied to the atoms are generated with the same tri-axial coils, also 3D printed, of 5 cm diameter. The whole setup is enclosed in a three-layer cylindrical magnetic shielding made of mu-metal, of 12 cm inner diameter.

### 3 Strong pumping toward aligned states: anisotropic relaxation and its dressing

The three-step approach is a useful method for describing the evolution of optically pumped systems. This approach is however known to be limited to low light intensities [26–28]. For depopulation pumping of a spin 1/2 system with circularly polarized light, it was shown long ago that this approach remains valid for higher light intensities if an additional relaxation term $\gamma_p$ due to the pumping, and proportional to light intensity $I$ is added to the natural relaxation $\Gamma_e$.

Some authors have discussed if a similar approach may be valid for a spin > 1/2 pumped by linearly polarized light towards an aligned state [38]. We study here if adding well-chosen relaxation terms to the three-step approach model could yield accurate predictions for the photodetection signals of helium-4 interacting with linearly polarized light. We first study the simplest scheme which is the Hanle effect. Then we consider the case of parametric resonance magnetometers.

#### 3.1 Hanle effect

We will consider here the situation sketched in Fig. 1(a). Let us recall the multipole decomposition of the density matrix on the irreducible tensor operators $\hat{T}_q^{(k)}$ [39]:

$$ \hat{\rho} = \sum_{k,q} m_q^{(k)} \hat{T}_q^{(k)\dagger} \qquad (1) $$

with $m_q^{(k)} = \langle \hat{T}_q^{(k)} \rangle$. The order $k = 1$ corresponds to orientation and $k = 2$ to alignment. The latter can be described by a column matrix $M$ with five components $m_{-2}^{(2)}$, $m_{-1}^{(2)}$, $m_0^{(2)}$, $m_1^{(2)}$, $m_2^{(2)}$.
Within the three-step approach, the evolution of *M* is given by [25, 35]:

$$ \frac{d}{dt} M = \mathbb{H}(B_0) \cdot M - \Gamma(M - M_{ss}) \qquad (2) $$

with $\mathbb{H}(B_0)$ the matrix describing the magnetic evolution, the total relaxation rate $\Gamma = \Gamma_e + \gamma_p$ and the pumping steady-state is $M_{ss} = \gamma_p/\Gamma(0, 0, 1/\sqrt{6}, 0, 0)^t$. The optical pumping rate $\gamma_p = \eta I$ with $I$ the light intensity and $\eta$ a constant. The evolution of *M* causes the light absorption to vary. For an input intensity $I_0$, the intensity at the output of the cell will be $I_0 - \Delta I$ with [40]:

$$ \Delta I = I_0 \alpha_r \left[ \frac{m_0^{(0)}}{\sqrt{3}} + \frac{m_0^{(2)}}{\sqrt{6}} - \Re(m_2^{(2)}) \right], \qquad (3) $$

where $\alpha_r$ a constant which depends on the optical transition, and the quantization axis is set along the propagation axis $z$.

In order to probe this approach we have recorded a set of Hanle-effect resonance curves for growing light intensities from 35 µW/cm² up to 1.96 mW/cm². For each light power, a quasi-static ramp of the $B_z$ component of the magnetic field is applied to the cell, and the photodetection signals are recorded with a National Instruments DAQmx board at the output of the transimpedance amplifier. These Hanle resonance curves are shown in Fig. 2.

This data is compared to the theoretical predictions resulting from Eqs. (2) and (3). The only free parameters of this model are the proportionality constants $\eta$ and $\alpha_r$ which are common for the whole data set: we have fitted them on the resonance curve corresponding to the lowest light intensity. The resulting theoretical curves show a fair agreement for the lowest light intensities, which is progressively degraded for larger light intensities. A qualitative explanation for these disagreements is the following. In this model the effective relaxation $\gamma_p$ does only depend on the light intensity, and not on the system alignment *M*. This is only accurate as far as the system is far from a fully aligned state, i.e. as far as optical power remains very low.

For larger optical powers, the complete equations describing optical pumping [41, 42], suggest that it may be possible to refine the previous model by considering the anisotropic nature of the effective relaxation brought by the pumping. The equation of evolution thus becomes:

$$ \frac{d}{dt} M = \mathbb{H}(B_0) \cdot M - \mathbb{R} \cdot M + \gamma_p M_{ss}, \qquad (4) $$

where $\mathbb{R} = \mathbb{R}^{(2)} + \Gamma_e \mathbb{I}$ with $\mathbb{I}$ the identity matrix and $\mathbb{R}^{(2)}$ is a 5 × 5 matrix which accounts for relaxation due to the optical pumping. From the complete equations of optical pumping [42], this latter is found to be for $D_0$ line of helium-4

$$ \mathbb{R}^{(2)} = \begin{pmatrix} 0 & 0 & 0 & 0 & 0 \\ 0 & \frac{3}{2}\gamma_p & 0 & 0 & 0 \\ 0 & 0 & 2\gamma_p & 0 & 0 \\ 0 & 0 & 0 & \frac{3}{2}\gamma_p & 0 \\ 0 & 0 & 0 & 0 & 0 \end{pmatrix} \qquad (5) $$
<figure><img src="image_2.png" /><figcaption>**Figure 2** Comparison of the three-step approach and the complete model for predicting Hanle resonances shapes at various optical powers. The five first plots show the photodetected optical power $P_{pd}$ as a function of the magnetic field, in units of the magnetic linewidth. Experimental data (in black) is compared to the theoretical expectations resulting from the three-step approach comprising only with isotropic relaxation (blue line) and complete model including anisotropic relaxation resulting from the pumping (red line). The bottom right plot shows the dependence of the full-width at half maximum of the curves as a function of the optical power for the experiment and the two theoretical models. The optical powers are given with 10% systematic uncertainty coming from optical power calibration</figcaption></figure>

and the pumping steady-state $M_{ss} = (0, 0, \sqrt{2/3}, 0, 0)^t$, both with quantization axis parallel to the light polarization. The pumping rate is:

$$ \gamma_p = I_0 \frac{2\pi^2 r_e c f}{\hbar\omega} \Re[\mathcal{V}(\omega - \omega_0)] \qquad (6) $$

with $r_e$ the electron classical radius, $f$ the oscillator strength of the transition, $\omega$ and $\omega_0$ the angular frequencies of the light and of the transition respectively, $\Re$ is the real part and $\mathcal{V}$ the complex Voigt profile function. Since the same beam is used as pump and probe, it is possible to write $\alpha_r$ in terms of $\gamma_p$: $\alpha_r = 3\alpha\gamma_p$ with $\alpha = \omega\ln\hbar/I_0$, where $l$ is the inner cell length and $n_m$ the density of atoms in the metastable state [40].

A first measurement at 9 µW optical power allows to deduce the metastable density $n_m$ and a second one at 37 µW, the ratio between the effective pumping rate in the double-pass configuration and the $\gamma_p$ given by Eq. (6). Except for these independently-measured parameters, the model contains no free parameters at all. We have calculated the resonance curves from Eq. (4) for the optical powers used in the experiments described above (Fig. 2). The agreement with experimental data is much better than with the previous
model. Figure 2 also shows the widths of the Hanle resonances. For growing powers the experimentally measured widths become narrower than the three step approach prediction. Including an anisotropic relaxation term as explained above allows a very good agreement between the predicted widths and those observed.

## 3.2 Parametric resonance

We now consider the case of practical interest which is the PRM. Thanks to the dressed atom formalism, PRM equations of motion can be mapped to the ones of the Hanle effect of the atom dressed by the two RF fields [16, 25].

For optical powers larger than those described by the three step approach, this equivalence remains valid if, in addition to the usual conditions, all the elements of the $\mathbb{R}$ matrix to remain much lower than the RF frequencies $\omega$ and $\Omega$. If this condition is fulfilled, the dressed anisotropic relaxation matrix $\mathcal{R}$ can be found to be [43]:

$$ \mathcal{R}_{q,q'} = \sum_{r=-\infty}^{+\infty} J_{r,q} \mathbb{R}_{q,q'} J_{r,q'} = J_{0,|q-q'|} \mathbb{R}_{q,q'}, \quad (7) $$

where the second equality comes from Graf's theorem [44]. The PRM equations can then be solved as usual.

For single-RF parametric resonance (the scheme of Fig. 1(b) with $B_2 = 0$), the absorption signal at odd harmonics of the RF frequency $\omega$ is sensitive to the $B_z$ component of the field. From Eq. (5), (7) and [25] one can find around null field:

$$ \Delta I = \frac{48I_0\alpha}{2\Gamma_e + 3\gamma_p} \times \frac{\gamma B_z \Gamma_e \gamma_p^2 J_{0,2} J_{p,2}}{(2\Gamma_e + \gamma_p)(2\Gamma_e + 3\gamma_p) - 3\gamma_p^2 J_{0,2}^2} \sin(p\omega t), \quad (8) $$

where $p$ is an odd integer and $J_{q,q'} = J_q(q'\gamma B_1/\omega)$ with $J_q(x)$ the Bessel function of first kind and order $q$.

To probe if this description of PRM is accurate we have acquired a set of parametric resonance dispersive response curves at different optical powers, in the presence of a single RF field. To do so we slowly varied the $B_z$ component of the magnetic field in the presence of a RF field along the $z$ axis, of 20 kHz frequency and $B_1 = \omega/(2\gamma)$. The photodetection signal at the output of the transimpedance amplifier was sent to the input of a Zurich MFLI lock-in amplifier. The output of this lock-in is then acquired with the same National DAQmx board used before.

The resulting resonance curves are plotted with crosses in Fig. 3(a). The theoretical predictions with no free parameters are superposed to this data set, and show an excellent agreement with the experimental data in the central null-field region. The agreement becomes slightly worse for high light intensities and larger magnetic field. We believe this may be due to slight imperfections on the cancellation of the transverse components of the magnetic field.

The Fig. 3(b) displays the slopes of the linear low-field portion of the Fig. 3(a). The agreement is very good for light power up to 1 mW (i.e. an intensity of 3.9 mW/cm²), which is well above the light power which optimizes the magnetometer signal [14]. For larger powers a disagreement is clearly visible, probably due to the incipient saturation of the optical transition. Indeed the saturation intensity in our system is 7.42 mW/cm², corresponding to an integrated power of 2.1 mW.
<figure><img src="image_3.png" /><figcaption>**Figure 3** Parametric resonance magnetometer response: comparison between the theoretical prediction (with no free parameters) and the experimental data. (a) Superposition of experimental measurements (color disks) for different optical powers (given in the inset, with 10% systematic uncertainty coming from optical power calibration) and theoretical predictions (solid lines). (b) Slopes of the dispersive response near the null magnetic field, as a function of the optical power: experimental data (blue crosses) and theoretical predictions (blue solid line)</figcaption></figure>

The two RF case can be addressed in the very same way, dressing the atom with the faster RF field first, then rotating the quantization axis from z to y and finally dressing the atom with the slower RF field. The resulting relaxation matrix is $\mathcal{R} = \Gamma_v \mathbb{I} + \gamma_p \mathbb{D}$, where $\mathbb{D}$ is

$$ \begin{pmatrix} \frac{3}{4}(1+J_{0,2}) & 0 & \sqrt{\frac{3}{32}(1+J_{0,2})\mathcal{J}_{0,2}} & 0 & 0 \\ 0 & \frac{3}{8}(3-J_{0,2}) & 0 & \frac{3}{8}(1+J_{0,2})\mathcal{J}_{0,2} & 0 \\ \sqrt{\frac{3}{32}(1+J_{0,2})\mathcal{J}_{0,2}} & 0 & \frac{1}{4}(5-3J_{0,2}) & 0 & \sqrt{\frac{3}{32}(1+J_{0,2})\mathcal{J}_{0,2}} \\ 0 & \frac{3}{8}(1+J_{0,2})\mathcal{J}_{0,2} & 0 & \frac{3}{8}(3-J_{0,2}) & 0 \\ 0 & 0 & \sqrt{\frac{3}{32}(1+J_{0,2})\mathcal{J}_{0,2}} & 0 & \frac{3}{4}(1+J_{0,2}) \end{pmatrix}. \quad (9) $$

## 4 Breakdown of the RWA

The dressed atom formalism relies on a RWA which requires $\omega \gg \Omega$, $\gamma B_1$, $\gamma B_0$, $\Gamma$ and $\Omega \gg \gamma B_0$, $\Gamma$. In practice, this latter condition may not be well achieved. In the case of orientation-based PRM, the breakdown of RWA is somehow beneficial because it yields a resonance at the interharmonics $\omega \pm \Omega$ which is sensitive to the third component of the magnetic field, which otherwise cannot be measured [15, 16].

We study here the case of alignment-based PRM. The calculation can be performed in a very similar way as done for orientation [16]: after dressing the atom with the fast RF field, one obtains a set of differential equations for the dressed alignments $\bar{M} =$
$(\overline{m}_{-2}^{(2)}, \overline{m}_{-1}^{(2)}, \overline{m}_{0}^{(2)}, \overline{m}_{1}^{(2)}, \overline{m}_{2}^{(2)})^t$ as explained in [25]. For instance the equation for $\overline{m}_{2}^{(2)}$ is:

$$ \frac{d}{dt} \overline{m}_{2}^{(2)} + [\Gamma - 2i\overline{\omega}_{y} - 2i\overline{\Omega}_{1} \cos(\Omega\tau)] \overline{m}_{2}^{(2)} = (i\overline{\omega}_{+} \overline{m}_{1}^{(2)} + \Gamma \overline{m}_{2,ss}^{(2)}), \quad (10) $$

where we kept the same notations used in [25] i.e. $\overline{\omega}_{x,y} = \gamma J_{0,2}B_{x,y}$, $\overline{\omega}_{z} = \gamma B_{z}$, $\overline{\Omega}_{1} = \gamma J_{0,2}B_{1}$, $\overline{\omega}_{\pm} = \omega_{z} \pm i\omega_{x}$, with $J_{q,q'} = J_{q}(q'\gamma B_{1}/\omega)$ and $J_{q}(x)$ is the Bessel function of first kind of order $q$ and argument $x$.

These equations can be written in an integral form, for instance Eq. (10) yields:

$$ \overline{m}_{2}^{(2)} = \exp\left[-\Gamma t + 2i\overline{\omega}_{y}t + 2i\frac{\overline{\Omega}_{1}}{\Omega}\sin(\Omega t)\right] \times \int_{-\infty}^{t} \exp\left[\Gamma\tau - 2i\overline{\omega}_{y}\tau - 2i\frac{\overline{\Omega}_{1}}{\Omega}\sin(\Omega\tau)\right] \times (i\overline{\omega}_{+}\overline{m}_{1}^{(2)} + \Gamma\overline{m}_{2,ss}^{(2)}) d\tau. \quad (11) $$

Developing these equations to the first order in $\gamma B_{x}/\Omega$, $\gamma B_{y}/\Omega$, $\gamma B_{z}/\Omega$ and $\Gamma/\Omega$, and replacing the dressed alignments $\overline{M}$ appearing in the integrals with the zero-order solutions of [25] we obtain:

$$ \begin{aligned} \overline{m}_{0}^{(2)} &= \overline{m}_{0,ss}^{(2)} + \sqrt{6} \left[ \left( \frac{\overline{\omega}_{z}}{\Omega} \Re \overline{m}_{1}^{(2)} + \frac{\overline{\omega}_{x}}{\Omega} \Im \overline{m}_{1}^{(2)} \right) C_{0} \right. \\ & \left. + \left( \frac{\overline{\omega}_{x}}{\Omega} \Re \overline{m}_{1}^{(2)} - \frac{\overline{\omega}_{z}}{\Omega} \Im \overline{m}_{1}^{(2)} \right) N_{0} \right], \\ \Im \overline{m}_{1}^{(2)} &= \left( \frac{\Gamma}{\Omega} \overline{m}_{1,ss}^{(2)} - \sqrt{\frac{3}{2}} \frac{\overline{\omega}_{x}}{\Omega} \overline{m}_{0}^{(2)} \right) C_{1} + \left( -\frac{\overline{\omega}_{x}}{\Omega} \Re \overline{m}_{2}^{(2)} + \frac{\overline{\omega}_{z}}{\Omega} \Im \overline{m}_{2}^{(2)} \right) C_{3} \\ & \qquad + \left( \frac{\overline{\omega}_{z}}{\Omega} \Re \overline{m}_{2}^{(2)} + \frac{\overline{\omega}_{x}}{\Omega} \Im \overline{m}_{2}^{(2)} \right) N_{1} - \sqrt{\frac{3}{2}} \frac{\overline{\omega}_{z}}{\Omega} \overline{m}_{0}^{(2)} N_{2}, \\ \Re \overline{m}_{2}^{(2)} &= \left( \frac{\overline{\omega}_{x}}{\Omega} \Re \overline{m}_{1}^{(2)} + \frac{\overline{\omega}_{z}}{\Omega} \Im \overline{m}_{1}^{(2)} \right) N_{3} - \frac{\Gamma}{\Omega} \overline{m}_{2,ss}^{(2)} N_{4} \\ & \qquad + \left( \frac{\overline{\omega}_{x}}{\Omega} \Im \overline{m}_{1}^{(2)} - \frac{\overline{\omega}_{z}}{\Omega} \Re \overline{m}_{1}^{(2)} \right) C_{2}, \end{aligned} \quad (12) $$

where $\Re$ and $\Im$ are the real and imaginary parts, $\overline{m}_{0}^{(2)}$ are the doubly dressed alignments and the expressions of $N_{i}$ and $C_{i}$ are:

$$ \begin{aligned} N_0 &= \sum_{q>0} 2 \sin 2q\Omega t \frac{\mathcal{J}_{2q,1}}{2q}, \\ N_1 &= \sum_{q,r>0} 2 \sin 2r\Omega t \frac{\mathcal{J}_{q,1}(\mathcal{J}_{q-2r,1} - \mathcal{J}_{q-2r,1})}{q}, \\ N_2 &= \sum_{q,r>0} 2 \sin 2r\Omega t \frac{\mathcal{J}_{q,1}(\mathcal{J}_{q+2r,1} - \mathcal{J}_{q-2r,1})}{q}, \\ N_3 &= \sum_{q,r>0} 2 \sin 2r\Omega t \frac{\mathcal{J}_{q,1}(\mathcal{J}_{q+2r,2} - \mathcal{J}_{q-2r,2})}{q}, \\ N_4 &= \sum_{q,r>0} 2 \sin 2r\Omega t \frac{\mathcal{J}_{q,2}(\mathcal{J}_{q+2r,2} - \mathcal{J}_{q-2r,2})}{q}, \end{aligned} \quad (13) $$
$$C_0 = \sum_{q>0} 2\cos(2q-1)\Omega t \frac{\mathcal{J}_{2q-1,1}}{2q-1},$$

$$C_1 = \sum_{q,r>0} 2\cos(2r-1)\Omega t \frac{\mathcal{J}_{q,1}(\mathcal{J}_{q+2r-1,1} + \mathcal{J}_{q-2r+1,1})}{q},$$

$$C_2 = \sum_{q,r>0} 2\cos(2r-1)\Omega t \frac{\mathcal{J}_{q,1}(\mathcal{J}_{q+2r-1,2} + \mathcal{J}_{q-2r+1,2})}{q},$$

$$C_3 = \sum_{q,r>0} 2\cos(2r-1)\Omega t \frac{\mathcal{J}_{q,1}(\mathcal{J}_{-q+2r-1,1} + \mathcal{J}_{-q-2r+1,1})}{q}$$

with $\mathcal{J}_{q,q'} = J_q(q'\gamma\bar{B}_2/\Omega)$ with $\bar{B}_2 = J_{0,1}B_2$.

The variation of light intensity contains the following terms:

$$3\left(\frac{\overline{m}_0^{(2)}}{\sqrt{6}} - \frac{\Re\overline{m}_2^{(2)}}{3}\right) \cos \frac{2\omega_1 \sin \omega t}{\omega} - \frac{\overline{m}_0^{(2)}}{\sqrt{6}} + \Re\overline{m}_2^{(2)} - 2\sin \frac{2\omega_1 \sin \omega t}{\omega} \Im\overline{m}_1^{(2)}. \quad (14)$$

Putting together Eqs. (12), (13) and (14) in a symbolic calculation software like Mathematica we have studied the non-secular terms which appear at the three frequencies $\omega$, $\Omega$, $\omega \pm \Omega$ which are used for magnetic field measurement.

The only term which is constant in magnetic field, and thus could cause a fixed offset, is $\Re\overline{m}_2^{(2)}$. Expanding Eq. (14) with Jacobi-Anger relations [44] shows however that all the resulting terms are at even harmonics of $\omega$.

All other terms have at least first-order dependence on the magnetic field. We have analyzed if these terms could bring spurious sensitivities to components other than $B_z$, $B_y$ and $B_x$ for the three frequencies $\omega$, $\Omega$ and $\omega \pm \Omega$ respectively. This is not the case. Indeed one can notice in Eqs. (13) that for all the ones with a $\sin(p\Omega t)$ time-dependence $p$ is even and all the ones with a $\cos(p\Omega t)$ time-dependence $p$ is odd. This is exactly the opposite for secular terms [25]. Therefore the terms resulting from RWA breakdown are in quadrature with the secular terms. An important example is the interharmonic term which allows the measurement of the $B_x$ component of the magnetic field. The secular sensitivity to $B_x$ has a time dependence $\sin(\omega t)\sin(\Omega t)$ [25]. The new term which appears here has a time dependence $\sin(\omega t)\cos(\Omega t)$.

We can therefore conclude that no systematic errors in magnetic field readings should result from the first order corrections of RWA.

## 5 Tensor light-shifts

The absence of vector light-shifts is often cited as an advantage of optical pumping with linearly instead of circularly polarized light [20, 25, 45]. However, a detuning between the linearly-polarized light and the atomic transition leads to a TLS. The effect of such a shift can be modelled by adding to the Hamiltonian a term [30, 31]:

$$\hat{H}_{\text{TLS}} = -\sqrt{6}Q_2\Delta E \hat{T}_0^{(2)} \quad (15)$$

with quantization axis along light polarization. $Q_2 = 1$ for $D_0$ and

$$\Delta E = I_0 \frac{\pi^2 r_e c f}{\hbar \omega} \Im[\mathcal{V}(\omega - \omega_0)]. \quad (16)$$
$\hat{H}_{TLS}$ is formally equivalent to the Hamiltonian of quadrupole effects [46]. For a scalar magnetometer based on atomic alignment [38] the $\hat{H}_{TLS}$ can be treated as a perturbation, like it is usually done for quadrupole effects [47]. At first order it causes a broadening of the magnetic resonance line, and to second order a shift of the center of the resonance. The resulting systematic error on scalar magnetometers of the ESA Swarm mission [20, 48, 49] is around 1 pT/mW in the worst case. This is negligible as compared to other systematic errors [50].

Our goal here is to study if there is a systematic error arising from TLS for the PRM. From Eq. (15) Ehrenfest theorem yields:

$$ \begin{cases} \frac{d}{dt} m_{\pm 1}^{(1)} = -i \langle [\hat{T}_{\pm 1}^{(1)}, \hat{H}_{TLS}] \rangle = \mp 3 Q_2 i \Delta E m_{\pm 1}^{(2)}, \\ \frac{d}{dt} m_{\pm 1}^{(2)} = -i \langle [\hat{T}_{\pm 1}^{(2)}, \hat{H}_{TLS}] \rangle = \mp 3 Q_2 i \Delta E m_{\pm 1}^{(1)}. \end{cases} \qquad (17) $$

These terms describe alignment-to-orientation conversion. Adding these terms in the differential equations of Hanle effect presented in [25], and calculating the new steady-state yields cumbersome equations, which can be developed to first-order in each magnetic field component, yielding the alignment:

$$ \begin{align} m_0^{(2)} &= \sqrt{\frac{1}{6} \frac{\gamma_p}{\Gamma_e + 2\gamma_p}}, \nonumber \\ m_1^{(2)} &= \omega_y \frac{2\Gamma_e + 3\gamma_p}{\Gamma_e + 2\gamma_p} \frac{2\gamma_p}{(2\Gamma_e + 3\gamma_p)^2 + 36\Delta E^2}, \tag{18} \\ m_2^{(2)} &= \gamma_p \frac{(2\Gamma_e + 3\gamma_p)(4i\omega_z + 2\Gamma_e + 3\gamma_p) + 36\Delta E^2}{(\Gamma_e + 2\gamma_p)[(2\Gamma_e + 3\gamma_p)^2 + 36\Delta E^2]}. \nonumber \end{align} $$

Comparing these solutions to the one in the absence of TLS ($\Delta E = 0$) one sees that TLS causes no field offset, since $\Delta E$ is only involved in the factors dividing the magnetic field. The resonance curve is therefore only widened in the presence of a non null TLS.

The effect of a TLS on PRM can be deduced from the above considerations. Indeed PRM can be modelled as the Hanle effect of the dressed atom. The tensor light-shift rate $\Delta E$ is dressed in the same way as the pumping rate $\gamma_p$. Once dressed, the problem is formally equivalent to the Hanle effect problem considered above. The photodetection signal de-modulated at the RF frequency has an expression analog to Eq. (18) showing a widening but no shift in the presence of a TLS.

# 6 Discussion and conclusion

We have studied above three effects which are disregarded in the usual analysis of parametric resonance magnetometers.

For optimizing the signal-to-noise ratio of PRM it is usual to raise the pump intensity well above the limits of the three-step approach. This is a situation which has been broadly addressed in the field of non-linear magneto-optical rotation (NMOR) [51]. The most usual NMOR scheme is based on the Hanle effect observed in birefringence [27]. This system is usually modelled by calculating the evolution of the whole density matrix dynamics, including the excited state. This approach has been used with great success also for the non-linear Hanle effect observed in fluorescence [52].
Our starting point was the Hanle effect, but observed in the pump absorption signal. We probed a simpler model, which consists in refining the three-step approach by introducing an anisotropic relaxation resulting from optical pumping. This model, which contains no free parameters, is in good agreement with a set of experimental acquisitions on a broad range of optical powers, as long as optical saturation is not reached. This latter can indeed be let aside, since for our helium-4 magnetometers the best sensitivities are achieved for optical intensities lamost an order of magnitude lower than optical saturation.

A former work [53] on Hanle effect deals with a somehow complementary situation where an anisotropic relaxation is probed with low light intensities, yielding isotropic pump-induced relaxation.

From our model of the Hanle effect we calculated the PRM photodetection signals by studying how anisotropic relaxation can be treated in the dressed-atom formalism for RF. We also found an expression for the doubly-dressed anisotropic relaxation term, which should allow to perform further checks of the accuracy of this model on two-RF PRM.

A second aspect we studied is the breakdown of RWA in PRM. For double-resonance magnetometers it is well known that such breakdown causes a systematic error: the Bloch-Siegert shift [21]. For PRM based on orientation no such systematic errors appear, and these effects are even beneficial, since they provide an additional sensitivity to the third axis of the magnetic field [15]. We have shown that for alignment-based PRM neither systematic errors are expected to result from this effect.

Finally we studied TLS which appear when the linearly-polarized pumping light is not well tuned to the transition, We have shown that these TLS cause a widening of the resonance curves, but no shift of their center, in contrast with vector light-shifts. If this widening varies with time it can cause a measurement error when the PRM is operated in open loop, relying on an initial calibration of the response curve to deduce the field from the optical signals. However it is usual to operate helium-4 PRM in closed-loop mode [37], where a compensation field continuously cancels the three components of the magnetic field. In such a configuration, a slope variation does not cause an error in the measurement.

As a conclusion, a reasonably good model of the photodetection signals can be obtained by refining the previous model of alignment-based PRM with the three mechanisms treated here. None of these mechanisms has been found to cause systematic errors. Therefore, the origin and ultimate order of magnitude of the systematic errors in PRM remains an open question.

**Acknowledgements**

We acknowledge H. Gilles, G. Le Gal and R. Garcés Malonda for interesting discussions, F. Alcouffe for cell filling, S. Gidon and W. Fourcault for help with optical setups. This research work was supported by the French ANR via Carnot funding.

**Abbreviations**

OPM, Optically Pumped Magnetometer; PRM, Parametric Resonance Magnetometer; RWA, Rotating Wave Approximation; TLS, Tensor Light Shift; NMOR, Non-linear Magneto Optical Rotation.

**Availability of data and materials**

The data and models are available upon request to the corresponding author.

**Competing interests**

The authors declare that they have no competing interests.

**Authors' contributions**

FB developed the theoretical models with the help of APL. FB made the experiment. APL wrote the manuscript with the help of FB. All authors read and approved the final manuscript.
**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 31 December 2019 Accepted: 14 April 2020 Published online: 07 May 2020

**References**

1. Kominis IK, Kornack TW, Allred JC, Romalis MV. A subfemtotesla multichannel atomic magnetometer. Nature. 2003;422(6932):596-9.

2. Knappe S, Alem O, Sheng D, Kitching J. Microfabricated optically-pumped magnetometers for biomagnetic applications. J Phys Conf Ser. 2016;723:012055.

3. 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.

4. Xia H, Ben-Amar Baranga A, Hoffman D, Romalis MV. Magnetoencephalography with an atomic magnetometer. Appl Phys Lett. 2006;89(21):211104.

5. Boto E, Holmes N, Leggett J, Roberts G, Shah V, Meyer SS, Muñoz LD, Mullinger KJ, Tierney TM, Bestmann S, Barnes GR, Bowtell R, Brookes MJ. Moving magnetoencephalography towards real-world applications with a wearable system. Nature. 2018;555(7698):657-61.

6. Weis A, Bison G, Castagna N, Cook S, Hofer A, Kasprzak M, Knowles P, Schenker J-L. Mapping the cardiomagnetic field with 19 room temperature second-order gradiometers. In: 17th international conference on biomagnetism advances in biomagnetism—Biomag2010. Berlin: Springer; 2010. p. 58–61.

7. Broser PJ, Knappe S, Kajal D, Noury N, Alem O, Shah V, Braun C. Optically pumped magnetometers for magneto-myography to study the innervation of the hand. IEEE Trans Neural Syst Rehabil Eng. 2018;26(11):2226-30.

8. Dolgovskiy V, Fescenko I, Sekiguchi N, Colombo S, Lebedev V, Zhang J, Weis A. A magnetic source imaging camera. Appl Phys Lett. 2016;109(2):023505.

9. Allred JC, Lyman RN, Kornack TW, Romalis MV. High-sensitivity atomic magnetometer unaffected by spin-exchange relaxation. Phys Rev Lett. 2002;89(13):130801.

10. Wyllie R, Kauer M, Wakai RT, Walker TG. Optical magnetometer array for fetal magnetocardiography. Opt Lett. 2012;37(12):2247-9.

11. Sulai IA, DeLand ZJ, Bulatowicz MD, Wahl CP, Wakai RT, Walker TG. Characterizing atomic magnetic gradiometers for fetal magnetocardiography. Rev Sci Instrum. 2019;90(8):085003.

12. Alem O, Mhaskar R, Jiménez-Martínez R, Sheng D, LeBlanc J, Trahms L, Sander T, Kitching J, Knappe S. Magnetic field imaging with microfabricated optically-pumped magnetometers. Opt Express. 2017;25(7):7849-58.

13. Cohen-Tannoudji C, Dupont-Roc J, Haroche S, Laloë F. Diverses résonances de croisement de niveaux sur des atomes pompés optiquement en champ nul ii. Applications à la mesure de champs faibles. Rev Phys Appl. 1970;5(1):102-8.

14. Cohen-Tannoudji C, Dupont-Roc J, Haroche S, Laloë F. Diverses résonances de croisement de niveaux sur des atomes pompés optiquement en champ nul. i. Théorie. Rev Phys Appl. 1970;5(1):95-101.

15. Dupont-Roc J. Détermination par des méthodes optiques des trois composantes d'un champ magnétique très faible. Rev Phys Appl. 1970;5(6):853-64.

16. Dupont-Roc J. Etude théorique de diverses résonances observables en champ nul sur des atomes habillés par des photons de radiofréquence. J Phys. 1971;32(2-3):135-44.

17. Yudin VI, Taichenachev AV, Dudin YO, Velichansky VL, Zibrov AS, Zibrov SA. Vector magnetometry based on electromagnetically induced transparency in linearly polarized light. Phys Rev A. 2010;82(3):033807.

18. Bevilacqua G, Breschi E. Magneto-optic spectroscopy with linearly polarized modulated light: theory and experiment. Phys Rev A. 2014;89(6):062507.

19. Ingleby SJ, O'Dwyer C, Griffin PF, Arnold AS, Riis E. Vector magnetometry exploiting phase-geometry effects in a double-resonance alignment magnetometer. Phys Rev Appl. 2018;10(3):034035.

20. Lieb G, Jager T, Palacios-Laloy A, Gilles H. All-optical isotropic scalar 4He magnetometer based on atomic alignment. Rev Sci Instrum. 2019;90(7):075104.

21. Sudyka JI, Pustelny S, Gawlik W. Limitations of rotating-wave approximation in magnetic resonance: characterization and elimination of the Bloch-Siegert shift in magneto-optics. New J Phys. 2019;21(2):023024.

22. Slocum RE, Reilly FN. Low field helium magnetometer for space applications. IEEE Trans Nucl Sci. 1963;10(1):165-71.

23. Slocum RE. Zero-field level-crossing resonances in optically pumped 2³S₁ He 4. Phys Rev Lett. 1972;29(25):1642-5.

24. Slocum RE, Marton B. Measurement of weak magnetic fields using zero-field parametric resonance in optically pumped He4. IEEE Trans Magn. 1973;9(3):221-6.

25. Beato F, Belorizky E, Labyt E, Le Prado M, Palacios-Laloy A. Theory of a 4He parametric-resonance magnetometer based on atomic alignment. Phys Rev A. 2018;98(5):053431.

26. Kanorsky SI, Weis A, Wurster J, Hänsch TW. Quantitative investigation of the resonant nonlinear Faraday effect under conditions of optical hyperfine pumping. Phys Rev A. 1993;47(2):1220-6.

27. Weis A, Wurster J, Kanorsky SI. Quantitative interpretation of the nonlinear Faraday effect as a Hanle effect of a light-induced birefringence. J Opt Soc Am B. 1993;10(4):716-24.

28. Budker D, Gawlik W, Kimball DF, Rochester SM, Yashchuk VV, Weis A. Resonant nonlinear magneto-optical effects in atoms. Rev Mod Phys. 2002;74(4):1153-201.

29. Sulai IA, Wyllie R, Kauer M, Smetana GS, Wakai RT, Walker TG. Diffusive suppression of ac-Stark shifts in atomic magnetometers. Opt Lett. 2013;38(6):974-6.

30. Happer W, Mathur BS. Effective operator formalism in optical pumping. Phys Rev. 1967;163(1):12-25.

31. Cohen-Tannoudji C, Dupont-Roc J. Experimental study of Zeeman light shifts in weak magnetic fields. Phys Rev A. 1972;5(2):968-84.

32. Chalupczak W, Wojciechowski A, Pustelny S, Gawlik W. Competition between the tensor light shift and nonlinear Zeeman effect. Phys Rev A. 2010;82(2):023417.

33. Colegrove FD, Franken PA. Optical pumping of helium in the 3S₁ metastable state. Phys Rev. 1960;119(2):680-90.

34. Rutkowski J, Fourcault W, Bertrand F, Rossini U, Gétin S, Le Calvez S, Jager T, Herth E, Gorecki C, Le Prado M, Léger JM, Morales S. Towards a miniature atomic scalar magnetometer using a liquid crystal polarization rotator. Sens Actuators A, Phys. 2014;216:386-93.
35. Breschi E, Weis A. Ground-state Hanle effect based on atomic alignment. Phys Rev A. 2012;86(5):053427.

36. Le Gal G, Lieb G, Beato F, Jager T, Gilles H, Palacios-Laloy A. Dual-axis Hanle magnetometer based on atomic alignment with a single optical access. Phys Rev Appl. 2019;12(6):064010.

37. Labyt E, Corsi M, Fourcault W, Palacios Laloy A, Bertrand F, Lenouvel F, Cauffet G, Prado ML, Berger F, Morales S. Magnetoencephalography with optically pumped 4He magnetometers at ambient temperature. IEEE Trans Med Imaging. 2019;38(1):90–8.

38. Weis A, Bison G, Pazgalev AS. Theory of double resonance magnetometers based on atomic alignment. Phys Rev A. 2006;74(3):033401.

39. Blum K. Density matrix theory and applications. Springer series on atomic, optical, and plasma physics. vol. 64. Berlin: Springer; 2012.

40. Laloë F, Leduc M, Minguzzi P. Relations entre l'état angulaire d'une vapeur atomique soumise au pompage optique et ses propriétés d'absorption et de dispersion. J Phys. 1969;30(2-3):277-88.

41. Barrat JP, Cohen-Tannoudji C. Étude du pompage optique dans le formalisme de la matrice densité. J Phys Radium. 1961;22(6):329-36.

42. Omont A. Irreducible components of the density matrix. Application to optical pumping. Prog Quantum Electron. 1977;5:69-138.

43. Haroche S. Etude théorique et expérimentale des propriétés physiques d'atomes habillés par des photons de radiofréquence. PhD thesis. Univ. Paris VI; 1971.

44. Abramowitz M, Stegun IA, editors. Handbook of mathematical functions: with formulas, graphs, and mathematical tables. 0009-revised ed. New York: Dover; 1965.

45. Leger J-M, Bertrand F, Jager T, Fratter I. Spaceborn scalar magnetometers for Earth's field studies. In: 62nd international astronautical congress 2011, IAC 2011. vol. 4; 2011. p. 2671-6.

46. Peck SK, Lane N, Ang DG, Hunter LR. Using tensor light shifts to measure and cancel a cell's quadrupolar frequency shift. Phys Rev A. 2016;93(2):023426.

47. Volk CH, Mark JG, Grover BC. Spin dephasing of Kr83. Phys Rev A. 1979;20(6):2381-8.

48. Guttin C, Leger JM, Stoeckel F. An isotropic earth field scalar magnetometer using optically pumped helium 4. J Phys IV. 1994;4(C4):655-9.

49. Leger J-M, Bertrand F, Jager T, Le Prado M, Fratter I, Lalaurie J-C. Swarm absolute scalar and vector magnetometer based on helium 4 optical pumping. Proc Chem. 2009;1(1):634.

50. Jager T, Léger J-M, Bertrand F, Fratter I, Lalaurie J-C. Swarm absolute scalar magnetometer accuracy: analyses and measurement results. In: 2010 IEEE sensors. 2010. p. 2392-5.

51. Auzinsh M, Budker D, Rochester S. Optically polarized atoms: understanding light-atom interactions. Oxford: Oxford University Press; 2010.

52. Papoyan AV, Auzinsh M, Bergmann K. Nonlinear Hanle effect in Cs vapor under strong laser excitation. Eur Phys J D. 2002;21(1):63-71.

53. Castagna N, Weis A. Measurement of longitudinal and transverse spin relaxation rates using the ground-state Hanle effect. Phys Rev A. 2011;84(5):053421.

**Submit your manuscript to a SpringerOpen**
**journal and benefit from:**

► Convenient online submission

► Rigorous peer review

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article

Submit your next manuscript at ► springeropen.com