Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Evaluation of optical parameters for a microminiature Rb vapor cell in a dual-beam SERF magnetometer

Open Access Open Access

Abstract

In the spin-exchange relaxation-free (SERF) magnetometer of a perpendicular pump-probe configuration, the pump and probe beam characteristics significantly affect the performance. In this paper, an efficient evaluation of optical parameters to improve the sensitivity of a miniature magnetometer has been presented. We have determined the pump light’s optimal intensity and wavelength through theoretical analysis and the zero-field resonance experiments. Chirp signals are applied to measure the optical rotations at different probe intensities and frequencies. Through theoretical and experimental analysis of noise source characterization under different beam intensities and wavelengths, we demonstrate that dual-beam magnetometer performance is mainly limited by photon shot noise. Based on the optimum pump and probe beam parameters, we demonstrate magnetic field sensitivity of 6.3 ${\textrm {fT/}}\sqrt {\textrm {Hz}}$ in an 87Rb vapor cell filled with nitrogen gas, with an active measurement volume of 3 × 3 × 3 mm3.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

The dominant technology for high-sensitive magnetic field measurements have been liquid-helium cooled superconducting quantum interference device (SQUID) magnetometer for nearly half a century due to their exquisite advantages, i.e., the high sensitivity of approximate 150 ${\textrm {aT/}}\sqrt {\textrm {Hz}}$ in the white noise regime between about 20 ${\textrm {kHz}}$ and a large bandwidth over 2.5 ${\textrm {MHz}}$ [1]. These sensitive magnetometers have already been used in biomagnetic signal detection, such as magnetoencephalography (MEG) [2] and magnetocardiography (MCG) [3]. Although the SQUID technology has been fully developed and SQUID-based magnetometers have been commercially available, the prerequisite cryogenics and the continual supply of expensive liquid helium limit their reconfiguration and increase the operating cost. Over the last 20 years, atomic magnetometers based on the interaction of resonant magneto-optics with atomic vapor show the capacity of high sensitivity in magnetic field measurements. Especially the descendant, spin-exchange relaxation-free (SERF) magnetometers have advantages in the detection of the weak magnetic field with the level of subfemtotesla [4,5], which eliminates the relaxation due to spin-exchange collision and achieves superior sensitivity performance [6,7]. Except for precisely fundamental measurements [8,9], the sensitivity and bandwidth of SERF magnetometer are sufficiently high to observe MCG [1012] and MEG [13,14] signals. In recent years, a commercial 15 ${\textrm {fT/}}\sqrt {\textrm {Hz}}$ in 3-100 ${\textrm {Hz}}$ band QuSpin Zero-Field Magnetometer (QZFM) was released, which is compact, low cost, and easy to operate [15]. Combined with the inherent advantages of being free from cryogenic cooling and allowing flexible positioning of the sensor, tens to hundreds of microfabricated SERF magnetometers can be easily placed around a subject to create a spatial map of the magnetic field at a low cost [16,17].

The application above requires the magnetometer to be miniaturized enough to obtain sufficient spatial resolution and ultra sensitivity to detect the weak magnetic signals. Considering spatial resolution consideration, we make a microfabricated vapor cell with an inner volume of $3 \times 3 \times 3\;{\textrm {mm}^3}$. But in early articles, the sensitivity of magnetometers is below 8 ${\textrm {fT/}}\sqrt {\textrm {Hz}}$ [15,18] with the same measurement volume, which is much larger than the theoretical limited sensitivity, about 1 ${\textrm {fT/}}\sqrt {\textrm {Hz}}$ [19]. So it is worthwhile to analyze the factor restraining the sensitivity of magnetometers. With respect to improve sensitivity, there are two approaches to be considered. One is narrowing the magnetic linewidth and increasing transverse relaxation time $T_{2}$. In our previous work [20], we studied the single-beam SERF magnetometer with $^{87}$Rb vapor cells under different nitrogen gas pressures, and achieved a high magnetic sensitivity of 8.8 ${\textrm {fT/}}\sqrt {\textrm {Hz}}$. The other is improving the noise-to-signal ratio, maximizing magnetometer response, and suppressing the system noise. A differential method to measure the optical rotation has been employed to reduce the noise from the fluctuation of probe beam. Furthermore, the performance of magnetometer is affected by additional types of noise including electronic and magnetic noise. The pioneering work by Romalis $et al.$ [19,21] analyzed the quantum noise, derived the noise limit, and presented a theoretical optimization of the magnetometer. Krzyzewski $et al.$ [22] proposed the test methods of various types of noise and experimentally evaluated noise arising from the laser intensity, laser frequency, and other sources of single-beam SERF magnetometer, but they did not give the analysis and test conclusions of the optimum intensity. Xing $et al.$ [23] presented a probe noise characteristic model demonstrating the relationship between the intensity and the signal frequency but they did not analyze the combined effect caused by the probe wavelength and intensity. Moreover, the photo-elastic modulator modulated system technique is unsuitable for the microfabricated structure in this paper.

In this paper, an efficient approach to improve the sensitivity of a miniature magnetometer has been presented. We analyze and measure the different noise sources of a microfabricated vapor cell with the power spectral density and propose a method combining multiple noise sources using the noise floor to optimize the probe beam characteristics. The different noise sources are characterized by converting noise spectrum to magnetic field units using a magnetometer signal response curve. We modify the shot noise and electronic noise equation to analyze the combined effect caused by the probe wavelength and intensity. We demonstrate the optimal pump and probe beam characterization based on the theoretical and experimental results. This method is helpful in noise analysis and could promote the sensitivity of the SERF magnetometer.

2. Theory analysis

The Bloch equation is the theoretical basis of the spin dynamics of the dual-beam SERF magnetometer, describing the fundamental physical processes of optical pumping, spin relaxation, and spin precession [7,19].

$$\frac{d}{{dt}}{\mathbf{P}} = \frac{1}{q}({\gamma ^{\textrm{e}}}{\mathbf{P}} \times {\mathbf{B}} + {R_{\textrm{OP}}}\left( {s{\mathbf{\hat z}} - {\mathbf{P}}} \right){\textrm{ + }}{R_{\textrm{PR}}}(s{\mathbf{\hat x}} - {\mathbf{P}}) - {R_{\textrm{rel}}}{\mathbf{P}}),$$
where $\gamma ^{\mathrm {e}}$ is gyromagnetic ratio of the bare electron, $q$ is nuclear slow-down factor. $\boldsymbol {P}$ is atomic polarization, $\boldsymbol {B}$ is magnetic-field vector. Photon spin operator $s=\pm 1$ for the circularly-polarized pump beam, $s=0$ for the linearly polarized probe beam. $R_{\textrm {OP}}$ and ${R_{\textrm {PR}}}$ are the optical pumping rate of pump and probe beam, respectively. $R_{\textrm {rel}}$ is relaxation rate.
$$R = \frac{{I{r_e}fc}}{{hv}}\frac{{{\Gamma _{\textrm{D}}}/2}}{{\Delta {v^2} + {{({\Gamma _{\textrm{D}}}/2)}^2}}},$$

The pumping rate for pump and probe beam can be expressed above [19], where $I$ is the intensity of the laser beam, ${{r_e}}$ is the classical electron radius, $f$ is the oscillator strength, ${\Delta v}$ is the frequency detuning from the corresponding transition. ${{\Gamma _{\textrm {D}}}}$ is the Full widths at half maximum of Lorentzian curve due to pressure broadening and natural lifetime. For convenience, we define ${R_{\textrm {PR}}} = \eta {I_{\textrm {probe}}}$ where the $\eta$ is the optical pumping factor which is influenced by the detuning.

We solve the Eq. (1) in steady-state. After the magnetic field in the shield is nulled and only applied small magnetic field of $B_{y}$ in $y$ axis, we can simplify the solution as:

$${P_x} = \frac{{{R_{\textrm{OP}}}}}{{\Gamma ({\Gamma ^2} + {\Omega ^2})}}(\Gamma {\Omega _y} + {\Omega _x}{\Omega _z}) \approx \frac{{{R_{\textrm{OP}}}{\Omega _y}}}{{{\Gamma ^2}}},$$
where $\Omega _{x}=\gamma ^{\mathrm {e}} B_{x}$, $\Omega _{y}=\gamma ^{\mathrm {e}} B_{y}$, $\Omega _{z} = \gamma ^{\mathrm {e}} B_{z}$, $\Omega = {\Omega _{\textrm {x}}}{\textrm { + }}{\Omega _{\textrm {y}}}{\textrm { + }}{\Omega _{\textrm {z}}}$, $\Gamma =R_{\textrm {OP}}+R_{\textrm {rel}}+R_{\textrm {PR}}$.

Linearly polarized probe light on ${\textrm {D}_2}$ off-resonance frequency along the $x$ axis is used to detect the optical rotation $\theta$, given by the Eq. (4). After the cell the probe beam passes through a polarizing beam splitter setting at ${45^ \circ }$ to the initial polarization.

$$\theta = \frac{1}{4}l{r_e}nc{f_{{\textrm{D2}}}}\frac{{\Delta v}}{{\Delta {v^2} + {{({\Gamma _{\textrm{D}}}/2)}^2}}}\frac{{{R_{{\textrm{OP}}}}{\Omega _{\textrm{y}}}}}{{{\Gamma ^2}}} = G(\Delta {v_{{\textrm{PR}}}},{I_{{\textrm{PR}}}}){B_{\textrm{y}}},$$
where $n$ is rubidium vapor density, $l$ is optical path length, $G$ is the magnetometer response, ${f_{{\textrm {D2}}}}$ is the oscillator strength. In small-angle approximation, combining the equations above and considering the attenuation of probe intensity, the output voltage signal can be written as :
$$V = 2K{I_{{\textrm{PR}}}}{e^{ - {\textrm{O}}{{\textrm{D}}_{{\textrm{PR}}}}}}G{B_{\textrm{y}}},$$
where $\textrm { OD }$ is the optical depth. $K$ is the conversion factor between light intensity and voltage of photodiode amplifier.

The performance of magnetometer is mainly limited by quantum noise, electronic noise, and magnetic noise. With a magnetic field $B$ to be measured, we evaluate those noise sources in magnetic field units for simplification as in Eq. (6).

$$\delta B = \sqrt {\delta B_{{\textrm{psn}}}^2 + \delta B_{{\textrm{spn}}}^2 + \delta B_{{\textrm{m}}}^2 + \delta B_{{\textrm{e}}}^2},$$

The fundamental quantum noise mainly contain spin projection noise and photon shot noise which can be expressed as $\delta {B_{{\textrm {spn}}}} = \sqrt {2{\Gamma ^3}} /({R_{{\textrm {OP}}}}{\gamma ^{\textrm {e}}}\sqrt {nV} )$ and $\delta {B_{{\textrm {psn}}}} = {k_{{\textrm {psn}}}}/(G\sqrt {{I_{{\textrm {PR}}}}{e^{ - {\textrm {OD}}}}} )$ respectively [19,21,24], $k_{{\textrm {psn}}}$ is the ratio of the shot noise unrelated with the probe intensity. The electronic noise is mainly from the instruments as PD, PDA and DAQs, according to Eq. (5), the equivalent magnetic noise can rearranged as $\delta {B_{\textrm {e}}} = {k_{\textrm {e}}}/({I_{{\textrm {PR}}}}{e^{ - {\textrm {OD}}}}G)$. Magnetic noise ${\delta {B_{\textrm {m}}}}$ is mainly from the magnetic shields and nulling coils. Combining the equations above, we can calculate the optimum optical parameters to ensure that the magnetometer achieves the highest sensitivity as the Fig. 1 shows. The simulation sensitivity is analyzed by the intensity and detuning of the probe beam, which could be calculated from Eq. (6) for probe beam center wavelength $\lambda$ = 780.255 ${\textrm {nm}}$, relaxation rate of the cell ${R_{\textrm {rel}}}$ = 100 Hz, full width at half maximum (FWHM) of pressure broadening $\Gamma _{{\textrm {D}}}$ = 24.5 GHz, and the rubidium vapor density $n$ = $1 \times {10^{14}}\;{{\textrm {cm}}^{ - 3}}$.

 figure: Fig. 1.

Fig. 1. The simulation sensitivity as a function of probe intensity and detuning are calculated from Eq. (6) for ${R_{\textrm {rel}}}$ = 100 Hz, $\Gamma _{{\textrm {D}}}$ = 24.5 GHz, rubidium vapor density $n$ = $1 \times {10^{14}}\;{{\textrm {cm}}^{ - 3}}$.

Download Full Size | PDF

3. Experimental setup

A schematic of the experimental setup is shown in Fig. 2, which is similar to the one used in our earlier studies [20]. Selecting borosilicate with the advantages of a low melted point and low cost as cell glass material, we made a microminiature alkali cell with an inner volume $3 \times 3 \times 3 {\textrm {m}}{{\textrm {m}}^3}$, containing a droplet of $^{87}\textrm {Rb}$ metal. In order to meet the needs of long-term testing, this cell is filled with N$_2$ buffer gas with a pressure of 1.4 amg. FWHM $\Gamma _{\textrm {D}}$ of this cell is 24.5 GHz, the center wavelength of ${{\textrm {D}}_1}$ line is 794.9984 nm, and the center wavelength of ${\textrm {D}_2}$ line is 780.2561 nm. Optical depth is 6.38 for pump beam which is measured by spectral profiles at 423 K. The cell is heated by a wire heater using high-frequency (200 kHz) alternating current as high as 150 $^{\circ }$C inside a boron nitride ceramic oven, and the temperature, which is monitored by a PT1000 resistor, is stable to 0.01$^{\circ }$C and vary less than 1$^{\circ }$C across the cell. The rubidium vapor density $n$ is about $1 \times {10^{14}}\;{{\textrm {cm}}^{ - 3}}$. The magnetic shield consists of four layers of $\mu$-metal and an inner layer of MnZn ferrite is constructed to attenuate the earth’s magnetic field and maintain the SERF regime. A triaxial uniform field Coils’ system is established to further compensate the residual magnetic field and provide the external magnetic field [25]. The radial coil is designed with a quasi-elliptic function fitting method in a simple structure with the constant of 74.55 ${\textrm {nT/mA}}$ [26] and the axial coil is designed with a Taylor expansion method with the constant of 92.32 ${\textrm {nT/mA}}$ [27]. As a design result, the absolute maximum relative error in the cubic target region (with a side length of $R$/2, where the radius of the ferrite $R$ = 45.5 mm) is optimized to 0.023$\%$. The coil constant is optimized considering the coupling effect between the coil and the ferrite material.

 figure: Fig. 2.

Fig. 2. Schematic of the experimental setup for measurements. A circularly polarized laser beam at $^{87}\textrm {Rb}\;{\textrm {D}_1}$ resonance line optical pumps the rubidium atoms along the $z$ direction, and the transmitted light intensity is monitored by PD0. And a linearly polarized probe beam which is blue detuned from the ${{\textrm {D}}_2}$ line, propagates along the $x$ direction, passing sequentially through the Rb cell, a polarizing beam splitter, and two subtraction photodiodes (PD+ and PD-). LP denotes linearly polarized plate; QWP is a quarter-wave plate; PD is a photodiode; PBS is a polarized beam splitter; DAQ is the data acquisition system. Shutters are used to turn on/off the probe and pump beam when measuring the system and probe noise.

Download Full Size | PDF

A circularly polarized laser beam at $^{87}\textrm {Rb}\;{\textrm {D}_1}$ line optical pumps the atoms along the $z$ direction. And a linearly polarized probe beam which is tuned to the ${\textrm {D}_2}$ line, propagates along the $x$ direction, passing sequentially through the Rb cell, a polarizing beam splitter, and two subtraction photodiodes (PD+ and PD-). The diameter $d$ of the pump and probe beam is 2.7 mm, and the linewidth of both is on the order of 1 MHz. The polarization beam splitter and noise eater combination are used to control the probe and pump laser intensity which changes less than 0.1$\%$ during the experiment. The extinction ratio of two beams passing through the fibers is larger than 40 dB. In the experiment, all data are collected by DAQs for 100 s with the sampling rate of 1000 Hz and analyzed with power spectral density of averaging spectral resolution of 1 Hz.

4. Results and discussion

In order to evaluate the noise characteristics under optimized condition, we optimize the pump intensity and wavelength first. On-resonance circular wavelength used to pump the alkali atoms is 794.9984 nm, which avoids the creation of a virtual magnetic field by near-resonant light. Although the AC-Stark shift can be compensated by the reverse magnetic field generated by the coil, gradient magnetic field is difficult to cancel due to attenuation of near-resonance light intensity especially in a cell with a large $\textrm { OD }$ parameter [28]. After determining the wavelength of the pump light, we optimize the optimal pump intensity. The fundamental sensitivity limit of the SERF magnetometer is given by the quadratic sum of the spin-projection noise, the photon-shot noise, electronic noise, and magnetic noise. We see that if the quantum limit is dominated by the spin-projection noise, the fundamental sensitivity is minimized by setting ${R_\textrm {OP}} = 2({R_\textrm {rel}}{\textrm { + }}{R_\textrm {PR}})$, where the polarizability is 0.67. If the photon shot noise, electronic noise dominates, then the sensitivity limit is optimized by setting ${R_\textrm {OP}} = {R_\textrm {rel}}{\textrm { + }}{R_\textrm {PR}}$, where the polarizability is 0.5, which in general result in the largest signal strength described in the Eq. (3). Under the experimental conditions of this paper. We simulate $\delta {B _\textrm {psn}}{\textrm { < }} 1\textrm {fT}/\sqrt {\textrm {Hz}}$ and $\delta {B _\textrm {spn}}{\textrm { > }} 4\textrm {fT}/\sqrt {\textrm {Hz}}$. The ${R_\textrm {OP}} = {R_\textrm {rel}}{\textrm { + }}{R_\textrm {PR}}$ is optimal choice for the following experiments. After nulling the magnetic fields, we sweep the amplitude of ${B_y}$ from −80 nT to 80 nT. Figure 3(a) shows the optical rotation signal of the cell at different incident pump beam intensity, while the probe beam intensity remain 8.73 ${\textrm {mW/c}}{{\textrm {m}}^2}$ and about 100 GHz detuning from the center wavelength. Combining the Eq. (3) with Eq. (5), we can simplify the solution as:

$$V = k\frac{{\eta {I_{{\textrm{pump}}}}}}{{{{(\eta {I_{{\textrm{pump}}}} + {R_{{\textrm{rel}}}}{\textrm{ + }}{R_{{\textrm{PR}}}})}^2}}},$$
where the slope $k$ is the scale factors, the data for different pumping intensity are well described by a dispersion curve. As we can see, pump intensity influence the peak-to-peak ${B_y}$ amplitude and the slope $k$ at the center portion. The slope $k$ at zero-field resonance, fitted by linear center portion when $\left | {{B_y}} \right | {\textrm { < }} 0.1\;{\textrm {nT}}$, vary from different pump power, which is shown in Fig. 3(b). The slope $k$ represents the DC-mode magnetometer’s sensitivity, and gets maximal value when the pumping rate equals spin-relaxation rate, namely ${R_{{\textrm {OP}}}} = {R_{{\textrm {rel}}}}{\textrm { + }}{R_{{\textrm {PR}}}}$. The black dots are measured data, while the red solid curve is fitted by Eq. (7). The optimized pump power can be ascertained under this condition.

 figure: Fig. 3.

Fig. 3. (a): Optical rotation signal of probe beam as a function of $B_{y}$ at different pump beam intensity. The probe intensity is 7${\textrm {mW/c}}{{\textrm {m}}^2}$. (b): The slope $k$ of the center of magnetometer response by sweeping ${B_y}$ amplitude at different pump beam intensity. The black dots are experimental data, and the red solid curve is fitted by Eq. (7).

Download Full Size | PDF

Using the method above, we can also calculate the optical pumping rate of the probe beam, which will be used in the following calculation. The half-width at half maximum (HWHM) of the Lorentzian curve represents the whole relaxation rate $\Gamma$. When we change the intensity of the probe beam, the slope represents the optical pumping factor $\eta$ as shown in Fig. 4(a). Figure 4(b) shows the optical pumping factor as a function of wavelength, which is discussed in Eq. (2). The red Lorentzian curve fits the black experimental data very well, but the width of the fitted Lorentzian curve is not consistent with that of pressure broadening of the cell. The reason is that the intensity of the probe light decays along the direction of the probe beam in the alkali cell. The smaller the detune is, the higher the attenuation of the probe intensity is. thereby, the actual measured value $\eta$ departs from theoretical value seriously. For reasonability, we use Lorentzian lineshape for fitting in this experiment.

 figure: Fig. 4.

Fig. 4. (a): The optical pumping rate of probe beam as a function of intensity at different wavelengths. The pump intensity is 7.43 ${\textrm {mW/c}}{{\textrm {m}}^2}$. (b): The optical pumping factor of probe beam under different wavelengths.

Download Full Size | PDF

As shown in Eq. (6), the fundamental noise are related to the factor $G$ except for the constant value of magnetic noise ${\delta {B_{\textrm {m}}}}$. At the condition ${R_{{\textrm {OP}}}} = {R_{{\textrm {rel}}}}{\textrm { + }}{R_{{\textrm {PR}}}}$ when magnetometer has the best sensitivity, the magnetometer response $G$ can be simplified as:

$$G = \frac{\kappa D(v) }{{{R_{{\textrm{rel}}}} + {R_{{\textrm{PR}}}}}},$$
where the factor $\kappa = l{r_e}nc{f_{{\textrm {D2}}}}{\gamma ^{\textrm {e}}}/16$, $D(v)$=$\Delta v/(\Delta {v^2} + {({\Gamma _D}/2)^2})$. Obviously, the coefficient $G$ is affected by the intensity of the probe beam due to the change of the pumping rate of probe beam ${R_{{\textrm {PR}}}}$, it is also affected by the wavelength of probe beam as a result of absorption cross-sectional area changed. The frequency response is used to calculate the sensitivity of the magnetometer. A chirp signal in the $x$ direction, varied from 1 Hz to 200 Hz in 50 s period with the amplitude of 100 ${{\textrm {pT}}_{{\textrm {rms}}}}$, is used to calibrate the frequency response under different probe beam intensities [29]. The frequency reponse can be described by a first-order low-pass filter with a fit to $A/{({(f - B)^2} + f_{3{\textrm {dB}}}^2)^{1/2}}$ where $f$ is the frequency and ${f_{3{\textrm {dB}}}}$ is the 3 - dB bandwidth [30]. As Fig. 5(a) shows, the magnetometer response is inversely proportional to the probe intensity, which is due to the pumping rate of probe beam ${R_{{\textrm {PR}}}}$. Then, we apply a 100 ${{\textrm {pT}}_{{\textrm {rms}}}}$ magnetic field with a frequency of 173.7 Hz to the magnetometer during the process of scanning the wavelength. The demodulation signal through the lock-in amplifier divided by incident intensity of probe beam is the magnetometer response signal $G{e^{ - {\textrm {OD}}}}$, taking into account absorption of the probe beam. These signals recorded as a function of probe wavelength at different probe beam intensities are shown in Fig. 5(b).

 figure: Fig. 5.

Fig. 5. (a): The magnetometer frequency response $G$ under different probe beam intensity at the wavelength of 780.183 nm. (b): The magnetometer response as a function of probe wavelength at different probe beam intensities, taking into account absorption of the probe beam.

Download Full Size | PDF

As shown by the photon shot noise, $\delta {B_{{\textrm {psn}}}} = {k_{{\textrm {psn}}}}/(G\sqrt {{I_{{\textrm {PR}}}}{e^{ - {\textrm {OD}}}}} )$, although increasing the intensity of probe beam can reduce the noise, as discussed above, it will cause the factor $G$ to become smaller so as to offset the noise reduction effect of increasing the intensity. Combined with the Eq. (8), the $\delta {B_{{\textrm {psn}}}}$ expands to:

$$\delta {B_{{\textrm{psn}}}} = \frac{{{k_{{\textrm{psn}}}}}}{{\kappa D(v)\sqrt {{e^{ - {\textrm{OD}}}}} }}\left( {\frac{{{R_{{\textrm{rel}}}}}}{{\sqrt {{I_{{\textrm{PR}}}}} }} + \eta \sqrt {{I_{{\textrm{PR}}}}} } \right),$$
Setting $\partial \delta B_{{\textrm {psn}}}^{}/\partial {I_{{\textrm {PR}}}} = 0$, we find the $\delta {B_{{\textrm {psn}}}} = 2{k_{{\textrm {psn}}}}\sqrt {{R_{{\textrm {rel}}}}\eta } /(\kappa D(v)\sqrt {{e^{ - {\textrm {OD}}}}} )$ is minimum for ${I_{{\textrm {PR}}}} = {R_{{\textrm {rel}}}}/\eta$. Accroding to the Eq. (9), the intensity characteristics of the photon shot noise is shown in Fig. 6. The optimal probe intensity corresponding to minimal shot noise becomes larger with the increasing detuning. At low detuning state, the optimal intensity is small. However, the shot noise increases rapidly with the increasing probe intensity after the optimum value due to the larger factor $\eta$. With the detuning increasing, not only does the optimum light intensity become larger, but the sensitivity changes more slowly as the probe intensity increases. In conjunction with Fig. 1, it can be seen that in the case of large detuning, the total noise of the dual-beam magnetometers is consistent with the photon shot noise, which also demonstrates that the magnetometer noise is mainly affected by the photon shot noise. We can optimize the photon shot noise $\delta {B_{{\textrm {psn}}}}$ by reducing the ${R_{{\textrm {rel}}}}$, for example, increasing the volume of cell, selecting Helium as buffer gas or encapsulating the cell with potassium.

 figure: Fig. 6.

Fig. 6. The simulation results of photon shot noise as a function of probe intensity at different wavelengths from Eq. (9)

Download Full Size | PDF

Through the formulation of electronic noise, it can be seen that the electronic noise is inversely proportional to the intensity of the probe beam. This attenuation coefficient of inverse proportional function is mainly affected by the detuning of the probe wavelength. As Fig. 7(a) shows, we measure the parameter $2K{I_{{\textrm {PR}}}}{e^{ - \textrm {OD}}}G$ as a function of probe wavelength at different probe beam. This quasi dispersion curve, in which the center wavelength is 780.255 nm, has two peaks, about 780.12 nm and 780.42 nm. The center zero area is due to the thick optical depth. As shown in Fig. 7(b), the total noise will be affected by electronic noise in the condition of slight detuning and weak intensity of the probe intensity. For example, the electronic noise is significant when the wavelength is 780.20 nm due to the small signal amplitude. When the probe intensity increases, the electronic noise approaches a small constant below 0.5 ${\textrm {fT/}}\sqrt {{\textrm {Hz}}}$ and can be neglected.

 figure: Fig. 7.

Fig. 7. (a): The signal amplitude of magnetometer as a function of probe wavelength at different probe intensities. (b): The electronic noise as a function of probe intensity at different wavelength.

Download Full Size | PDF

As discussed above, we should analyze the different noise sources ${\delta {B_{{\textrm {spn}}}}}$, ${\delta {B_{{\textrm {psn}}}}}$, ${\delta {B_{\textrm {e}}}}$ and ${\delta {B_{\textrm {m}}}}$. Assuming the magnetometer working under optimal conditions, the spin-projection noise ${\delta {B_{{\textrm {spn}}}}}$ is the level of 0.73 ${\textrm {fT/}}\sqrt {{\textrm {Hz}}}$. For the magnetic noise ${\delta {B_{\textrm {m}}}}$ consists of the coils’ noise and magnetic shield’s noise in our experiment condition, the magnetic coils’ noise floor which caused by the circuit voltage uncertainty is measured by applying a calibrating sine signal of 100 ${{\textrm {pT}}_{{\textrm {rms}}}}$ from the function generator and calculated using coil calibration constant. The magnetic shield consists of four layers of $\mu$-metal and an inner layer of MnZn ferrite with a calculated noise floor of 2 ${\textrm {fT/}}\sqrt {{\textrm {Hz}}}$ at 1 Hz [31,32]. The magnetic noise is a constant no matter what the performance of the magnetometer is. Due to the magnetometer frequency response, the electronic noise ${\delta {B_{\textrm {e}}}}$ measured by blocking off the pump and probe beam is different under different probe beam intensity. In our experiment, optical rotation signal is measured by a pair of difference PDs, which can estimate the noise from the fluctuation of the laser power. However, there are some offsets that remain between the two PDs, the optical rotation signal might be slightly influenced. We monitor the probe incident intensity to evaluate the probe intensity noise, including noises from the probe beam, fiber, quantum, and electronic devices. By blocking the pump beam, we can get the probe shot noise which contains probe intensity noise and photon shot noise, and define the photon shot noise by comparing the probe intensity noise with probe shot noise. Figure 8 shows the magnetic noise spectrum of probe beam, coils, electronic system, when the pump beam intensity is 7 ${\textrm {mW/c}}{{\textrm {m}}^2}$ and probe beam intensity is 11 ${\textrm {mW/c}}{{\textrm {m}}^2}$ with the wavelength of 780.12 nm. As we can see, the electronic and coils noise is about 0.1-0.2 ${\textrm {fT/}}\sqrt {{\textrm {Hz}}}$, and the probe intensity noise is about 1.5 ${\textrm {fT/}}\sqrt {{\textrm {Hz}}}$, while the probe shot noise is about 6 ${\textrm {fT/}}\sqrt {{\textrm {Hz}}}$. Compared to other noise sources, the performance of the magnetometer is mainly limited by the photon shot noise.

 figure: Fig. 8.

Fig. 8. Magnetic noise spectrum of probe photon shot noise, probe intensity noise, electronic noise, coils noise, and theoretical ferrite noise.

Download Full Size | PDF

To measure the sensitivity of the magnetometer and analyze the influence of probe beam, a small sine calibrating magnetic field at 30.5 Hz with an amplitude of 100 ${{\textrm {pT}}_{{\textrm {rms}}}}$ is applied at $x$-direction. We choose several probe beam wavelengths (780.02 nm, 780.065 nm, 780.12 nm, 780.165 nm, and 780.2 nm) to evaluate the magnetometer’s noise density spectrum under several different probe beam intensities as Fig. 9 shows. The five noise floors are too close to tell the difference under the relatively small probe intensity. For brief verification, we average the noise floor from 10 Hz to 100 Hz noise density spectrum except the calibrating and power frequency while the average range contains the magnetometer bandwidth.

 figure: Fig. 9.

Fig. 9. The noise density spectrum of the magnetometer at several different probe intensities while the probe beam wavelength is 780.165 nm. Calibrating magnetic field is 30.5 Hz, and the power frequency is 50 Hz.

Download Full Size | PDF

To analyze the noise of a dual-beam magnetometer, we have tested the sensitivity of the magnetometer under different probe intensities at different probe wavelengths plotted in Fig. 10. By comprehensively considering the several noises, we have fitted the data using Eq. (6). As shown in Fig. 10, the dots are experimental data and the solid line are fitted with good by the parameters ${R_{{\textrm {rel}}}}{\textrm { = 103 Hz}}$, ${k_{{\textrm {psn}}}} = 0.65 \pm 0.05\;{\textrm {fT/}}\sqrt {{\textrm {Hz}}} \sqrt {{\textrm {mW/c}}{{\textrm {m}}^{\textrm {2}}}}$, and ${k_{\textrm {e}}} = 0.025 \pm 0.005\;{\textrm {fT/}}\sqrt {{\textrm {Hz}}} {\textrm {(mW/c}}{{\textrm {m}}^{\textrm {2}}})$. The five fitting curves at different probe wavelengths are similar to the simulation in Fig. 1. Under different wavelengths, the sensitivities all decrease at first and increase subsequently. Meanwhile, when the detuning of the probe beam is small, the sensitivity is poor and the suitable light intensity range is also very small. When the detuning is bigger, the sensitivity gets better and the suitable light intensity range becomes bigger. For small detuning, e.g., 780.20 nm, the total noise is mainly affected by electronic and photon shot noise under low intensity. Although the electronic noise decreases with the probe intensity increasing, the photon shot noise increases rapidly due to the large factor $\eta$. The total noise floor is always large at different probe intensities. When the detuning gets larger, e.g., 780.12 nm and 780.165 nm, the electronic noise gradually decreases and does not affect total noise. The factor $\eta$ is relatively small, so the photon shot noise grows gradually with probe intensity. When the detuning is large, e.g., 780.02 nm and 780.065 nm, the trend of total noise with probe intensity is almost the same as the photon shot noise.

 figure: Fig. 10.

Fig. 10. The sensitivities under different probe intensities at different probe wavelengths. The results averaged by 10 Hz to 100 Hz noise density spectrum except the calibrating and power frequency. The dots are experimental data, and the solid lines are calculated from Eq. (6).

Download Full Size | PDF

In the case of large detuning, optimum sensitivity can be achieved at large probe intensity, so the transimpedance amplifier is required to have the ability of a large measuring range. Sometimes the magnetometer output has bias voltage due to the position of the polarization axis of the half-wave plate, and the transimpedance amplifier will overload with the probe intensity increasing. In the case of the larger probe intensity, the photodiode’s photoelectric conversion coefficient will turn in the nonlinear region, resulting in output error. Through the above analysis, the optimal parameter is about 50 - 70 GHz detuning from center wavelength and about 10 - 15 ${\textrm {mW/c}}{{\textrm {m}}^2}$ probe beam intensity, the best sensitivity is about 6.3 ${\textrm {fT/}}\sqrt {{\textrm {Hz}}}$.

5. Conclusion

In this paper, we present a theoretical and experimental analysis of noise sources characterized from a microfabricated cell with an active measurement volume of $3 \times 3 \times 3\;{{\textrm {mm}}^3}$ in a perpendicular pump-probe configuration. Then we use several noises of the magnetometer, including spin projection noise, photon shot noise, electronic noise, and magnetic noise to evaluate the probe beam characteristics and find that the optical pumping rate of the probe beam significantly affects the magnetometer performance. Compared to other noise sources, the performance of the magnetometer is mainly limited by the probe noise at the optimal pumping condition. The optimal probe intensity range increases with more detuning from the resonance frequency. However, limited by the performance of transimpedance amplifiers and photodetectors, the optimal probe intensity is not suitable for large detuning. The experimental results show that the optimal parameter is about 50 - 70 GHz detuning from center wavelength and about 10 - 15 ${\textrm {mW/c}}{{\textrm {m}}^2}$ probe beam intensity with a noise floor of 6.3 ${\textrm {fT/}}\sqrt {{\textrm {Hz}}}$. In order to further improve the sensitivity, the photon shot noise can be suppressed by reducing relaxation ${R_{{\textrm {rel}}}}$, such as replacing rubidium with potassium or filling the cell with Helium gas in the same volume. Fabricating a big size cell can also settle this matter, but this is not friendly to miniaturized integration of magnetometers. There is another sophisticated method of spin squeezing to enhance the magnetometer by reducing the photon noise [33,34]. We theoretically calculate the total noise and fit it with the experimental data in a good arrangement. Based on the results, we demonstrate the optimal pump and probe beam characterizations, which could promote the sensitivity in the SERF magnetometer.

Funding

Beijing Municipal Natural Science Foundation (1222025); Key Research and Development Program of Zhejiang, China under Grant (2020C01037); National Natural Science Foundation of China (62003022); National Key Research and Development Program of China under Grant (2018YFB2002405).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. J.-H. Storm, P. Hömmen, D. Drung, and R. Körber, “An ultra-sensitive and wideband magnetometer based on a superconducting quantum interference device,” Appl. Phys. Lett. 110(7), 072603 (2017). [CrossRef]  

2. R. Hari and R. Salmelin, “Magnetoencephalography: From squids to neuroscience: Neuroimage 20th anniversary special edition,” NeuroImage 61(2), 386–396 (2012). [CrossRef]  

3. H. Koch, “Squid magnetocardiography: status and perspectives,” IEEE Trans. Appl. Supercond. 11(1), 49–59 (2001). [CrossRef]  

4. W. C. Griffith, S. Knappe, and J. Kitching, “Femtotesla atomic magnetometry in a microfabricated vapor cell,” Opt. Express 18(26), 27167–27172 (2010). [CrossRef]  

5. J. Tang, Y. Zhai, L. Cao, Y. Zhang, L. Li, B. Zhao, B. Zhou, B. Han, and G. Liu, “High-sensitivity operation of a single-beam atomic magnetometer for three-axis magnetic field measurement,” Opt. Express 29(10), 15641–15652 (2021). [CrossRef]  

6. W. Happer and A. Tam, “Effect of rapid spin exchange on the magnetic-resonance spectrum of alkali vapors,” Phys. Rev. A 16(5), 1877–1891 (1977). [CrossRef]  

7. J. Allred, R. Lyman, T. Kornack, and M. V. Romalis, “High-sensitivity atomic magnetometer unaffected by spin-exchange relaxation,” Phys. Rev. Lett. 89(13), 130801 (2002). [CrossRef]  

8. I. Savukov and M. V. Romalis, “Nmr detection with an atomic magnetometer,” Phys. Rev. Lett. 94(12), 123001 (2005). [CrossRef]  

9. D. Bear, R. Stoner, R. Walsworth, V. A. Kosteleckỳ, and C. D. Lane, “Limit on lorentz and cpt violation of the neutron using a two-species noble-gas maser,” Phys. Rev. Lett. 85(24), 5038–5041 (2000). [CrossRef]  

10. G. Bison, R. Wynands, and A. Weis, “Dynamical mapping of the human cardiomagnetic field with a room-temperature, laser-optical sensor,” Opt. Express 11(8), 904–909 (2003). [CrossRef]  

11. G. Bison, N. Castagna, A. Hofer, P. Knowles, J.-L. Schenker, M. Kasprzak, H. Saudan, and A. Weis, “A room temperature 19-channel magnetic field mapping device for cardiac signals,” Appl. Phys. Lett. 95(17), 173701 (2009). [CrossRef]  

12. R. Wyllie, M. Kauer, G. Smetana, R. Wakai, and T. Walker, “Magnetocardiography with a modular spin-exchange relaxation-free atomic magnetometer array,” Phys. Med. Biol. 57(9), 2619–2632 (2012). [CrossRef]  

13. H. Xia, A. Ben-Amar Baranga, D. Hoffman, and M. Romalis, “Magnetoencephalography with an atomic magnetometer,” Appl. Phys. Lett. 89(21), 211104 (2006). [CrossRef]  

14. T. Sander, J. Preusser, R. Mhaskar, J. Kitching, L. Trahms, and S. Knappe, “Magnetoencephalography with a chip-scale atomic magnetometer,” Biomed. Opt. Express 3(5), 981–990 (2012). [CrossRef]  

15. I. Savukov, Y. J. Kim, V. Shah, and M. G. Boshier, “High-sensitivity operation of single-beam optically pumped magnetometer in a kHz frequency range,” Meas. Sci. Technol. 28(3), 035104 (2017). [CrossRef]  

16. G. Zhang, S. Huang, F. Xu, Z. Hu, and Q. Lin, “Multi-channel spin exchange relaxation free magnetometer towards two-dimensional vector magnetoencephalography,” Opt. Express 27(2), 597–607 (2019). [CrossRef]  

17. J. Iivanainen, R. Zetter, M. Grön, K. Hakkarainen, and L. Parkkonen, “On-scalp meg system utilizing an actively shielded array of optically-pumped magnetometers,” NeuroImage 194, 244–258 (2019). [CrossRef]  

18. V. K. Shah and R. T. Wakai, “A compact, high performance atomic magnetometer for biomedical applications,” Phys. Med. Biol. 58(22), 8153–8161 (2013). [CrossRef]  

19. M. Ledbetter, I. Savukov, V. Acosta, D. Budker, and M. Romalis, “Spin-exchange-relaxation-free magnetometry with cs vapor,” Phys. Rev. A 77(3), 033408 (2008). [CrossRef]  

20. Y. Wang, G. Jin, J. Tang, W. Zhou, B. Han, B. Zhou, and T. Shi, “Optimized gas pressure of an rb vapor cell in a single-beam serf magnetometer,” Opt. Express 30(1), 336–348 (2022). [CrossRef]  

21. I. M. Savukov, S. Seltzer, M. Romalis, and K. Sauer, “Tunable atomic magnetometer for detection of radio-frequency magnetic fields,” Phys. Rev. Lett. 95(6), 063004 (2005). [CrossRef]  

22. S. Krzyzewski, A. Perry, V. Gerginov, and S. Knappe, “Characterization of noise sources in a microfabricated single-beam zero-field optically-pumped magnetometer,” J. Appl. Phys. 126(4), 044504 (2019). [CrossRef]  

23. B. Xing, C. Sun, Z. Liu, J. Zhao, J. Lu, B. Han, and M. Ding, “Probe noise characteristics of the spin-exchange relaxation-free (serf) magnetometer,” Opt. Express 29(4), 5055–5067 (2021). [CrossRef]  

24. S. Smullin, I. Savukov, G. Vasilakis, R. Ghosh, and M. Romalis, “Low-noise high-density alkali-metal scalar magnetometer,” Phys. Rev. A 80(3), 033420 (2009). [CrossRef]  

25. S. Zou, W. Quan, X. Chen, J. Fang, and H. Zhang, “Determination of residual magnetic field based on optically-detected free spin precession of hyperpolarized 21ne,” IEEE Trans. Instrum. Meas. 71, 1–7 (2022). [CrossRef]  

26. W. Zhou, B. Han, J. Wang, D. Wu, F. Zhao, and N. Li, “Design of uniform magnetic field coil by quasi-elliptic function fitting method with multiple optimizations in miniature atomic sensors,” IEEE Trans. Ind. Electron. (posted 13 November 2021, in press).

27. W. Wu, B. Zhou, Z. Liu, J. Wang, H. Pang, L. Chen, W. Quan, and G. Liu, “Design of highly uniform magnetic field coils based on a particle swarm optimization algorithm,” IEEE Access 7, 125310–125322 (2019). [CrossRef]  

28. I. Sulai, R. Wyllie, M. Kauer, G. Smetana, R. Wakai, and T. Walker, “Diffusive suppression of ac-stark shifts in atomic magnetometers,” Opt. Lett. 38(6), 974–976 (2013). [CrossRef]  

29. J. Lu, Z. Qian, and J. Fang, “A fast determination method for transverse relaxation of spin-exchange-relaxation-free magnetometer,” Rev. Sci. Instrum. 86(4), 043104 (2015). [CrossRef]  

30. I. Kominis, T. Kornack, J. Allred, and M. V. Romalis, “A subfemtotesla multichannel atomic magnetometer,” Nature 422(6932), 596–599 (2003). [CrossRef]  

31. T. Kornack, S. Smullin, S.-K. Lee, and M. Romalis, “A low-noise ferrite magnetic shield,” Appl. Phys. Lett. 90(22), 223501 (2007). [CrossRef]  

32. S.-K. Lee and M. V. Romalis, “Calculation of magnetic field noise from high-permeability magnetic shields and conducting objects with simple geometry,” J. Appl. Phys. 103(8), 084904 (2008). [CrossRef]  

33. J. Geremia, J. K. Stockton, A. C. Doherty, and H. Mabuchi, “Quantum kalman filtering and the heisenberg limit in atomic magnetometry,” Phys. Rev. Lett. 91(25), 250801 (2003). [CrossRef]  

34. R. Sewell, M. Koschorreck, M. Napolitano, B. Dubost, N. Behbood, and M. Mitchell, “Magnetic sensitivity beyond the projection noise limit by spin squeezing,” Phys. Rev. Lett. 109(25), 253605 (2012). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (10)

Fig. 1.
Fig. 1. The simulation sensitivity as a function of probe intensity and detuning are calculated from Eq. (6) for ${R_{\textrm {rel}}}$ = 100 Hz, $\Gamma _{{\textrm {D}}}$ = 24.5 GHz, rubidium vapor density $n$ = $1 \times {10^{14}}\;{{\textrm {cm}}^{ - 3}}$.
Fig. 2.
Fig. 2. Schematic of the experimental setup for measurements. A circularly polarized laser beam at $^{87}\textrm {Rb}\;{\textrm {D}_1}$ resonance line optical pumps the rubidium atoms along the $z$ direction, and the transmitted light intensity is monitored by PD0. And a linearly polarized probe beam which is blue detuned from the ${{\textrm {D}}_2}$ line, propagates along the $x$ direction, passing sequentially through the Rb cell, a polarizing beam splitter, and two subtraction photodiodes (PD+ and PD-). LP denotes linearly polarized plate; QWP is a quarter-wave plate; PD is a photodiode; PBS is a polarized beam splitter; DAQ is the data acquisition system. Shutters are used to turn on/off the probe and pump beam when measuring the system and probe noise.
Fig. 3.
Fig. 3. (a): Optical rotation signal of probe beam as a function of $B_{y}$ at different pump beam intensity. The probe intensity is 7${\textrm {mW/c}}{{\textrm {m}}^2}$. (b): The slope $k$ of the center of magnetometer response by sweeping ${B_y}$ amplitude at different pump beam intensity. The black dots are experimental data, and the red solid curve is fitted by Eq. (7).
Fig. 4.
Fig. 4. (a): The optical pumping rate of probe beam as a function of intensity at different wavelengths. The pump intensity is 7.43 ${\textrm {mW/c}}{{\textrm {m}}^2}$. (b): The optical pumping factor of probe beam under different wavelengths.
Fig. 5.
Fig. 5. (a): The magnetometer frequency response $G$ under different probe beam intensity at the wavelength of 780.183 nm. (b): The magnetometer response as a function of probe wavelength at different probe beam intensities, taking into account absorption of the probe beam.
Fig. 6.
Fig. 6. The simulation results of photon shot noise as a function of probe intensity at different wavelengths from Eq. (9)
Fig. 7.
Fig. 7. (a): The signal amplitude of magnetometer as a function of probe wavelength at different probe intensities. (b): The electronic noise as a function of probe intensity at different wavelength.
Fig. 8.
Fig. 8. Magnetic noise spectrum of probe photon shot noise, probe intensity noise, electronic noise, coils noise, and theoretical ferrite noise.
Fig. 9.
Fig. 9. The noise density spectrum of the magnetometer at several different probe intensities while the probe beam wavelength is 780.165 nm. Calibrating magnetic field is 30.5 Hz, and the power frequency is 50 Hz.
Fig. 10.
Fig. 10. The sensitivities under different probe intensities at different probe wavelengths. The results averaged by 10 Hz to 100 Hz noise density spectrum except the calibrating and power frequency. The dots are experimental data, and the solid lines are calculated from Eq. (6).

Equations (9)

Equations on this page are rendered with MathJax. Learn more.

d d t P = 1 q ( γ e P × B + R OP ( s z ^ P )  +  R PR ( s x ^ P ) R rel P ) ,
R = I r e f c h v Γ D / 2 Δ v 2 + ( Γ D / 2 ) 2 ,
P x = R OP Γ ( Γ 2 + Ω 2 ) ( Γ Ω y + Ω x Ω z ) R OP Ω y Γ 2 ,
θ = 1 4 l r e n c f D2 Δ v Δ v 2 + ( Γ D / 2 ) 2 R OP Ω y Γ 2 = G ( Δ v PR , I PR ) B y ,
V = 2 K I PR e O D PR G B y ,
δ B = δ B psn 2 + δ B spn 2 + δ B m 2 + δ B e 2 ,
V = k η I pump ( η I pump + R rel  +  R PR ) 2 ,
G = κ D ( v ) R rel + R PR ,
δ B psn = k psn κ D ( v ) e OD ( R rel I PR + η I PR ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.