Sensitive and stable vector magnetometer for operation in zero and finite fields

We report on an optically pumped magnetometer that uses multiple laser beams to pump and probe spin-polarized Cs atoms. The selected sensor geometry allows for operation in finite magnetic fields as well as close to zero field. In finite fields the magnetometer employs free spin precession signals to determine the field modulus and direction as described in a separate publication. This publication focuses on the magnetometer operation close to zero field, which is based on a ground state Hanle resonance. The four laser beams permit the simultaneous measurement of two orthogonal magnetic field components in a differential detection scheme that greatly suppresses technical laser power noise. Sensitivities better than 54 fT/Hz1/2 could be demonstrated simultaneously for both measurement channels in a well shielded environment. A minimum Allan deviation, limited by residual field fluctuations, of better than 40 fT was observed for integration times of 2s. The magnetometer achieves high sensitivity and stability in offset fields as well as close to zero field and is, thus, a universal tool for low frequency magnetic field measurements. © 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement


Introduction
Optically pumped magnetometeres (OPM) [1,2] are among the most sensitive magnetic field sensors known today [3]. In such sensors, atomic spin polarization is created using optical pumping [4,5] and its change due to the interaction with the magnetic field is monitored by optical means. Most commonly, the sensors monitor vector spin polarization (spin orientation), which is created with circularly polarized pumping light [6] and precesses in the magnetic field with the Larmor frequency ω L = γ |B|. In this study we focus on 133 Cs atoms in the ground state with the overall atomic spin F = 4. The precession of that spin is characterized by the gyromagnetic ratio γ ≈ 3.5 Hz/nT and the effective transverse spin coherence time T 2 with its corresponding spin relaxation rate Γ = 1/T 2 .
The large variety of OPM sensors can be classified according to the magnetic field strength |B| they are designed for. The first class are zero-field sensors which operate at field strengths |B| Γ/γ. Their Larmor precession is much slower than the spin relaxation and their operation is based on the ground state Hanle effect [7]. The very sensitive SERF magnetometers [3,8] are a sub class of zero-field sensors. Zero-field sensors are usually sensitive to a single magnetic field component given by the geometry of the sensor. Using modulation techniques that periodically change one or several magnetic field components, it is possible to gain sensitivity to more than one component of the magnetic field vector [9].
A second class of optically pumped magnetometers is designed to operate in a finite magnetic fields |B| Γ/γ [10,11]. Under those conditions the Larmor precession is much faster than the spin-relaxation rate, which means that many Larmor periods can be observed during the coherence time using magnetic resonance techniques. Magnetometers in this regime include the classical M x and M z [11] designs that use oscillating magnetic fields to manipulate the precessing The modulation output of one lia is used to synchronize the other lia and to generate the modulation field B M . The in-phase output of lia 1..4 is digitized by the data acquisition (DAQ) PC, which also controls a function generator that can produce a ramping field B R . pd1 measures the light power transmitted along the +xz direction while pd2 measures along −xz, pd3 along +yz, and pd4 along −yz.
spins as well as all-optical variants where all necessary interactions are mediated by light [12]. Such magnetometers primarily measure the modulus of the magnetic field and can in addition be sensitive to the direction of the field [12,13].
The magnetometer setup used in this study is capable of delivering highly sensitive and stable readings in zero field as well as in finite fields. The performance in finite fields is published in [13], where we show that the magnetometer can be used for sensitive vector field measurements with 80 fT magnitude resolution and better than 8 µrad angular resolution for a few seconds integration time. In the work reported here we show that the same module can achieve similar performance in zero field when operated in Hanle mode.

Experimental setup
The experimental setup uses four laser beams, which all contribute to optical pumping and probing of the Cs spin. We describe the ensemble average F of the Cs spin as the magnetization M ∝ F of the medium. Figure 1 shows the magnetometer module with the integrated optics that create the four laser beams. The magnetometer module was placed inside the Biomagnetic Shielded Room (BMSR-2) [14] at PTB Berlin. Inside the shield, a system of Helmholtz-coils with 1.40 m and 1.60 m diameter could provide homogeneous magnetic fields along the x y and z axes, respectively. The active medium of the magnetometer consisted of Cs vapor at room temperature contained in an evacuated glass cell. The inside surface of the cell was coated with paraffin to prevent spin depolarization during wall collisions [15]. The light from a external-cavity diode laser, resonant with the F=4 to F=3 Cs D 1 transition (895 nm), was guided to the magnetometer head via a polarization-maintaining single-mode fiber. Close to the magnetometer module, a non-magnetic beam-splitting setup distributed the light to four multi-mode fibers with 0.4 mm core diameter and approx. 1 m length. The large core diameter was chosen because non-magnetic (plastic) setups do not permit coupling into single-mode fibers. The fibers were well fixed in order to minimize the conversion of mechanical vibrations to light intensity changes. Each of the four light beams in the magnetometer head was formed by an optical module (om in Fig. 1) that linearly polarizes and collimates the light from each multi-mode fiber. A second polarizer in the optical module ensured that each beam is circularly polarized with the same helicity when it enters the Cs cell. The power in each beam was individually adjusted to be approximately 1 µW by rotating the two polarizers in the optical module with respect to each other. All k-vectors of the beams had a polar angle of θ k =45°with respect to the z axis. Beams 1 and 2 lay in the xz-plane, while beams 3 and 4 lay in the yz-plane. The light power, transmitted through the Cs cell in each beam, was measured by a non-magnetic photo-diode (pd1 to 4) and a corresponding transimpendance amplifier (pa1 to 4). The experimental setup described so far in this section is identical to the one used in [13] to measure vector components of a magnetic field with a magnitude of |B| ≈ 1 µT, which was mainly oriented along z. In this study, we operated close to zero magnetic field, which was achieved by degaussing [16] the BMSR-2 and without the need for compensation coils. The remanent field in BMSR-2 B 0 = B 0x , B 0y , B 0z did not exceed the limit |B 0 | < 400 pT during the measurements.
As described in section 3, the magnetometer performance could be significantly improved by using phase-sensitive detection, which was implemented using four lock-in amplifiers referenced to the same modulation frequency ω M . The sinusoidal reference output of one lock-in amplifier was used to generate a modulated magnetic field B M = A M sin(ω M t) along z using the corresponding Helmholtz-coil. The modulation frequency and amplitude was experimentally optimized and chosen to be ω M = 2π · 10.217 Hz and A M = 2 nT. The reference phase of all lock-in amplifiers was chosen to be identical and optimized to yield the largest in-phase signal, which was recorded by the data acquisition system. Additionally, the output signals of the DC-coupled pre-amplifiers (pa 1..4) could be recorded as a measure of the DC light power in each beam.
Magnetic test signals were generated by an arbitrary waveform generator that supplied the Helmholtz-coils along x or y via a precision resistor R t . Slow ramps of the magnetic field B R along x or y were used to measure the line shape of the Hanle resonance (see Fig. 2 and 3). A simple calibration of the magnetometer response was gained by applying a square-wave to those coils with a known magnetic field step.

Signal modeling and data analysis
In this section we derive the line shapes of the Hanle resonances for the particular experimental geometry of our four-beam magnetometer module in order to compare it to experimental data. We start with the evolution of the magnetization vector M = (M x , M y , M z ) that is described by the Bloch equations where Γ is the effective relaxation rate, which includes the broadening due to optical pumping and assumes isotropic relaxation. The combined optical pumping of all four laser beams results in a source of magnetization along the z-direction, which is modeled as a relaxation towards M = (0, 0, M 0 ). In the following expressions, we parametrize spin relaxation using the magnetic line width δB = Γ/γ. The atoms' absorption coefficient for circularly polarized light contains a term proportional to the projection of the magnetization M onto the k vector of the considered beam. Consequently, the four beams probe different projections of M, which for the steady-state solutions ( M = 0) of Eq. (1) are given by Here M +xz and M −xz are the projections of the magnetization probed by the laser beams along the +xz-and the −xz-directions respectively. Correspondingly, M +yz and M −yz are the projections of the magnetization onto the the two beams in the yz-plane. In order to verify the validity of those line shapes, a slow ramp of the x field component B x = B 0x + B R was applied in the range from −35.5 nT ≤ B R ≤ 35.5 nT, while the transmitted light power of the four laser beams was recorded. The top row of Fig. 2 shows the preamplifier voltages, which are proportional to the transmitted light power recorded for each of the four laser beams. These signals are in principle usable for sensitive magnetometry in the regions, where they show a significant slope s = dM/dB with respect to the applied field. The most sensitive part of the slope for beams 3 and 4 are, however, not centered with respect to B R = 0 which can be corrected using the lock-in technique introduced in Sec. 2. In this mode of operation, a modulation field B m = A m sin(ω m t) was applied along the z-direction. The oscillating field leads to a corresponding modulation of the detected signals with an amplitude that depends characteristically on the field components in x-and y-directions. The lock-in amplifiers are used to demodulate the signal and generate an output that is proportional to the signal component that oscillates in phase with the modulation. The top row of Fig. 2 shows the modulated signals at the input to the lock-in amplifiers, while the bottom row shows the output demodulated signals. The most sensitive part of the slope in the demodulated signals for beams 1 and 2 are now centered respect to B R = 0 and the two signals from the same plane have very similar shapes except for a different sign. Subtracting the two mirrored signals further increases the signal amplitude and suppresses correlated noise. This differential detection scheme is discussed later and leads to the data presented in Fig. 3.
To model the output demodulated signals of the lock-in amplifier, we substitute B z → B z + B m in (Eq. 2), expand the expressions into a Taylor series and select the terms with sin(ω M t) up to the third order. We then take the difference between the expressions obtained for M +xz and M −xz (beam 1 and beam 2); and M +yz and M −yz (beam 3 and beam 4), respectively. This results in the line shapes of the differential output signals f xz and f yz , measured by beams in xz-and yzplanes: Here D = B 2 x + B 2 y + B 2 z + δB 2 is the typical resonance denominator. The terms proportional to A 3 M disappear if the Taylor series is expanded only to linear order. Close to zero field (|B| δB), the magnetometer signals simplify to linear zero crossings This approximation shows that the f xz -signal is mostly sensitive to the B x component of the field, and the f yz signal to the B y component. To first order both signals are not sensitive to B z . In order to compensate for possible misalignments between the laser beams that define the xzand yz-planes and the the coils producing the ramping field, the anticipated line shapes of the measured signals are parameterized by where the amplitudes A 1 , A 2 , A 3 and A 4 take into account mutual contributions of the differential signals.
Since the functions f are dimensionless, the amplitudes contain the overall laser power and scaling factors such as amplifier gains. They also permit to model differences in laser power in the different beams. The experimental line shapes are shown in Fig. 3 along with a simultaneous fit of the two functions from Eq. (5). One can conclude, that experimental data are well reproduced by the model predictions. In this fit, the quantities A 1...4 , B x...z , and δB were set as free parameters. This permits to determine values of the transverse components of the residual magnetic field B 0x and B 0y , which yields: B 0x = 0.329(32) nT, B 0y = −0.366(24) nT, δB = 3.483(43) nT. The given uncertainties are statistical 1σ errors determined by the fitting algorithm.
The field component B x = B 0x + B R was ramped very slowly (60 s ramping time) in order to record the data shown in Fig. 3, since the magnetometer has a limited bandwidth. The effective bandwidth is given by the time constant of the low-pass filter in the lock-in amplifier (T c = 100 ms with 24 db/octave), which has to be long enough to ensure an artefact-free demodulation of the signals. Thus, it has to be longer than a period of the modulation frequency, which itself has a lower bound given by the coherence time of the spin polarization. Recording a magnetic field ramp like in Fig. 3 is useful to verify the validity of the data model and to determine the linewidth δB. Once this is done, the field range (|B| δB) in which the signal can be aproximated by a linear zero-crossing (see Eq. 4) can be determined. In this range, the measured output values of the lock-in amplifier can be directly converted to magnetic field readings without the need for ramping fields. Given the fit results from above, the fields are approximately an order of magnitude smaller than the linewidth, which justifies this approximation. Since the linewidth δB can change with temperature and laser power, we calibrated the slope before and after taking long time seies by applying a known magnetic field step. In the next section we show an analysis of the magnetometers sensitivity and stability using recordings close to zero field and without the use of ramping fields.

Sensitivity and long term stability
In order to determine the sensitivity and stability of the magnetometer, we recorded the lock-in amplifier output for long periods of time (overnight). We prepared the measurement by carefully degaussing the BMSR-2 and tried to keep the field inside the room as constant as possible. The signal corresponded to the point B R = 0 in Fig. 3, where small changes of the field are linearly translated to signal changes according to Eq. (4). The resulting time series of magnetic field values for the x and y direction are depicted in Fig. 4 and show a typical behavior for the BMSR-2. The field in the x-direction is mainly influenced by the door and needs two hours after the door was closed to stabilize. Apart from that, field changes are dominated by the magnetic disturbances caused by the Berlin subway lines. Usually, the subway stops operating for a few Here several bins were averaged in order to reduce the scatter. The red curve is computed using a segment of the B x time series that shows particularly low noise (see text). The gray curve shows the original spectrum ρ(B x ) before low-pass filtering.
hours in the night and the best noise performance is obtained during that time. Unfortunately, construction works in the subway tunnels were carried out during our data taking, which caused disturbances almost constantly throughout the night. As indicated in Fig. 4, the following noise and stability analysis is based on a time series that covers the time from 11:42 pm till 08:02 am on December 15, 2016, which does not include the two hour period needed to let the door settle. Figure 5 shows the noise density spectrum of the measured magnetic field components B x and B y . In addition to the two spectra for the entire time series, the spectrum of a 3.4 h period starting at 1:09 am (red curve) of the B x time series is shown, which was selected for particularly low noise probably due to reduced subway activity during that time. All spectra show a white noise region extending down to ≈ 0.1 Hz. The selected segment shows the lowest white noise density averaged over a frequency span from 0.5 Hz to 1.1 Hz of 35 fT/Hz 1/2 . The noise shown in the time series that cover the whole night is slightly higher with ρ(B x ) = 44 fT/Hz 1/2 and ρ(B y ) = 54 fT/Hz 1/2 . The difference between the two channels is due to slightly different probe powers and the non-uniformity of the noise in BMSR-2.
For frequencies above 2 Hz, all original spectra showed an increase in noise due to vibrational resonances in the shielding room and an imperfect suppression of noise in the lock-in amplifier (see gray curve in Fig. 5). The imperfection was due to a reduced lock-in time constant of T c = 30 ms for this measurement which was chosen in order to maximize the usable bandwidth. In order to suppress that noise, all time series were filtered with a steep FIR low-pass filter [17] with a flat pass-band from DC to 1.9 Hz and an equivalent noise bandwidth of 2.3 Hz. The long time stability of the readings for B x and B y was characterized using Allan deviations (ADEV) of the same times series as in the spectra. Figure 6 shows the dependence of the ADEV on the integration time (σ − τ plot). All ADEV-plots were calculated with overlapping samples (see Sec. 5.2.2 in [18]) in order to reduce the scatter. The use of the low-pass filter induces correlations in the samples which lead to a rising ADEV at small integration times τ < 0.2 s. The white noise region for the measurements of B x (black and red curves in Fig. 5) are slightly sloped, while the one for B y is almost flat. As a consequence, only the ADEV for B y shows a segment between τ = 0.2 s and τ = 0.5 s that follows the expected τ −1/2 scaling. The ADEV for longer integration times is dominated by the instability of the magnetic field in the BMSR-2 and shows typical values for this shield. We expect that the curves would follow the statistical scaling for longer integration times if the experiment were repeated in a more stable field. Under the given conditions we find a minimum ADEV of 32 fT for τ = 1 s for the selected data set and slightly worse values of about 40 fT for the entire time series.

Conclusions
In summary, the results from our experiment demonstrate highly sensitive simultaneous measurements of the magnetic field vector components along two orthogonal directions. The measurements are based on a ground-state Hanle resonance, that is well described by the derived algebraic expressions for the line shapes. We find that the sensor geometry with four laser beams allows for efficient magnetometer operation in offset fields [13] as well in zero-field conditions. The use of two laser beams in one plane permits a differential detection scheme, which greatly suppresses the influence of common mode power fluctuations. This is particularly useful in the zero-field case, where the detection is based on signal components with frequencies of about 10 Hz which can be strongly affected by technical noise of the laser power and frequency. A lock-in detection scheme further suppresses instabilities in the zero-field case by shifting the detection frequency away from DC. We found a high stability of the readings with an almost flat part of the noise spectrum extending down to 0.1 Hz. In contrast to most SERF magnetometers, the four beam magnetometer is thus particularly well suited to detect low-frequency magnetic signals. The instability that leads to a rising Allan deviation (Fig. 6) for integration times τ > 1 s is typical for BMSR-2 and thus probably not induced by the magnetometer.
The Cs medium used is optimized for room-temperature operation with a narrow linewidth of about 3 nT. Using a SERF medium with the four beam pumping and readout scheme would result in a larger detection bandwidth and potentially higher sensitivity. However, the magnetometer cannot measure in an offset field if a medium in the SERF regime is used. The ability to switch from zero-field to finite-field operation without need to change anything in the sensor head allows for a very flexible use of the sensor. The device and methods are thus very useful as a general-purpose magnetometer and especially in applications that require operation in both magnetic field regimes. A possible application of this sensor is its use in precision experiment like the search for an electric dipole moment (EDM) of the neutron [19,20]. In normal operation such experiment use a coil producing a field of 1 µT is powered inside a strong magnetic shield. Without the need to change the magnetometer and open the shield the magnetometer may now be used to investigate with the coil powered off the properties of the shield, i.e. characterizing the quality of degaussing [16].
The magnetometer is well suited for all applications that require detecting small changes in the magnetic field. In precision experiments magnetic field changes that correlate with changes in the experimental conditions are often a source of systematic uncertainties. In EDM searches such uncertainties, for example, can be caused by magnetic field changes correlated with the reversal of the electric field (see Sec. 8B in [20]). The four beam magnetometer permits to study such correlations with or without offset field and thus gains additional information about the possible cause of the correlation. The noise spectrum in Fig. 5 shows that the highest sensitivity is obtained for detection frequencies between 1 Hz and 2 Hz. This fact is also reflected by the part of the ADEV in Fig. 6 that scales approximately statistically e.g. proportional to τ −1/2 . A detection frequency of 1 Hz corresponds to averaging the time series over segments of τ = 0.5 s. The curve ρ(B selected x ) in Fig. 6, which is least effected by external noise, shows that the magnetometer can compare one segment of length τ = 0.5 s to the next one with a statistical uncertainty of 35 fT.
The corresponding statistical uncertainty for τ = 0.5 s in a finite magnetic field is 110 fT according to the recordings presented in [13]. This measurement was taken in a pulsed mode of operation that is intrinsically not as sensitive as the continuous recording used in zero field. Continuous magnetic resonance techniques [11] can be adapted to the four beam setup and should further improve the performance in finite field mode.
In both modes of operation the ADEV curves demonstrate a high stability which permits to compare events that are separated in time by hundreds of seconds with a statistical uncertainty of better than 300 fT. This is important, for example, for studying correlations due to experimental parameters that cannot be quickly modulated like the reversal of a strong electric field.