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

In-situ measurement of the density ratio of K-Rb hybrid vapor cell using spin-exchange collision mixing of the K and Rb light shifts

Open Access Open Access

Abstract

We investigate a new method that enables the direct measurement of the density ratio of a K-Rb hybrid vapor cell, using the spin-exchange collision mixing of the K and Rb light shifts. The densities for each alkali metals can be further determined using Raoult’s law. The mixture of the light shifts in both magnetometers and comagnetometers is formulated using Bloch equations and explained by considering the fast spin-exchange interaction. The relationship between the density ratio and the mixed light shifts is both formulated and simulated. The method was performed on several K-Rb magnetometer- and K-Rb-21Ne comagnetometer-cells at different temperatures, pump light powers, and mole fractions of K. The method was further verified by the conventional laser-absorption-spectroscopy method. The new approach has the advantage to measure the density ratio of the optically-thick hybrid alkali atoms, while requiring no additional magnetic field necessary for conventional magnetic-field induced Faraday-rotation techniques. It also has the advantage of in-situ measuring the density ratio under exactly the normal operation of the devices, which means that the errors caused by the heating-effect of the strong pump light and the temperature drift during long-term operation can be real-time monitored.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Optical pumping, to produce strong atomic polarization, is used in many areas, partially including fundamental physics research [1,2], biomagnetism [3], and inertial navigation [4]. Unfortunately, in an optically thick medium, the strong absorption of pump light produces a significant polarization gradient and imposes a limit on the efficiency of optical pumping [5]. Recently, a hybrid spin-exchange optical-pumping (HSEOP) approach that uses a mixture of K and Rb was developed to overcome this difficulty, where the optically pumped low-density alkali metal was used to polarize the other high-density alkali metal via spin-exchange collisions [6,7]. In spin-exchange relaxation free (SERF) magnetometers and comagnetometers, the condition to achieve optimal sensitivity, demanding the optical pumping rate equaling the spin relaxation rate, cannot be maintained throughout the optical thick cell, not even across the intersection of the pump-probe beams [8, 9]. Using HSEOP, the optimal sensitivity can be achieved throughout the vapor by the virtue of uniform polarization. This large sensitive area is useful to yield multi-channel biomagnetic field measurements with fine spatial resolution [10,11]. Besides, in the spin-exchange optical pumping of noble gas such as 3He and 21Ne, HSEOP increases the pumping efficiency by nearly an order of magnitude, compared to the pure alkali-metal pumping method [7].

The main challenge for HSEOP is that both the homogeneity of polarization and the pumping efficiency are determined by the density ratio of the used hybrid alkali metals. The higher the density of the optically pumped alkali metal is, the lower is the polarization homogeneity. On the other hand, the lower the density is of the optically pumped alkali metal, the more pumping power is required. The motivation to accurately determine the density ratio of hybrid alkali metals vapor is to find the optimal density ratio that is necessary to improve both the sensitivity of the magnetometer [10,11] and the pumping efficiency in hyperpolarizing noble gas [7].

At present, the commonly used methods to measure atomic density in vapor cells are laser- absorption-spectroscopy (LAS) [12], Faraday-rotation [13], and measuring-spin-exchange-relaxation [14, 15]. The LAS method, where the density is derived from fitted absorption spectrum, is difficult to use with dense alkali vapor cells. This is because it is not only very challenging to directly measure the extremely weak transmitted light near the resonance, but also hard to indirectly infer the absorption line-shape from measurements in far-detuned wings due to near-complete attenuation of the light within many linewidths of resonance. Even at lower densities, the accuracy is limited by the broadening and distortion of the line-shape. In contrast, the magnetic-field-induced and polarization-induced Faraday-rotation methods can be used to monitor the density of dense vapor [13,16,17]. However, a very strong magnetic field (∼1 T) is needed for the magnetic-field-induced Faraday-rotation method [13], which would magnetize the magnetic field shields in some devices. Besides, the polarization-induced Faraday-rotation method requires an additional absolute measurement of the polarization [16,17]. The recently proposed measuring-spin-exchange-relaxation method measures the spin-exchange relaxation rate of alkali atoms, in the condition of a moderate magnetic field and low spin polarization, to derive the density [14, 15]. Unfortunately, measurement accuracy is affected by the spin polarization, and it cannot be used under exactly the normal operation of many devices with higher spin polarization. Furthermore, this method cannot measure the density in a hybrid alkali cell due to fast spin-exchange collisions.

In this paper, a new method is used to directly measure the density ratio of K-Rb in hybrid alkali vapor based on the mixture of the K and Rb light shifts caused by fast spin-exchange collisions. Firstly, the mixture of light shifts of K and Rb atoms in magnetometer and comagnetometer is described based on Bloch equations, and further explained by fast spin-exchange collisions with an intuitive model. We establish the relationship between the density ratio and the mixed light shifts, and present the measurement procedure of this method. Next, measurements of the density ratio using this method were performed with a K-Rb SERF magnetometer and a K-Rb-21Ne SERF comagnetometer at different temperatures, pump light powers and mole fractions of K. As a further check, the measured density ratios were compared with those measured using the LAS method. In the last section, we discuss limitations of this new method. This new method has the advantage that it does not need additional strong magnetic field and absolute polarization measurement, which are necessary for magnetic-field-induced and polarization-induced Faraday-rotation methods respectively, to directly measure the density ratio of the optically-thick hybrid alkali atoms. It also has the advantage of in-situ measuring the density ratio under exactly normal operations, thus the errors caused by the strong pump light heating-effect and the long-term temperature drift can be real-time monitored by this method. Besides, the cell temperature and alkali atoms densities in the K-Rb cell could be further determined by measuring the density ratio based on our method.

2. Basic principle

The energy levels of an atom can be shifted by a far-off-resonant laser beam, even in the absence of light absorption, which can be classified into scalar, vector, and tensor contributions [18]. Circularly-polarized light generates a vector light shift, while linearly-polarized light produces scalar and tensor shifts. The vector light shift affects the individual Zeeman sublevels separately, like a fictitious magnetic field along the light-propagation direction, while a scalar light shift shifts the Zeeman sublevels collectively. The tensor light-shift is negligible compared to the vector light shift [19]. In this work, we concentrate on the effect of the vector light shift (which we refer to simply as “light shift” from now on) of circularly polarized light. This study is performed on a K-Rb SERF magnetometer and a K-Rb-21Ne SERF comagnetometer. The K-Rb magnetometer consists of a vapor cell filled with a small mixed droplet of K-Rb and 4He gas (buffer gas) as well as N2 gas (quenching gas). The mole fraction of K in the hybrid K-Rb cell is kept small to obtain a relative low K density at high temperature, where the density ratio of K to Rb, Dr = nK/nRb, is typically on the order of 0.01. Here, nK and nRb are the densities of the K and Rb atoms, respectively. The low-density K atoms, whose absorption of pump light along the propagation direction is much smaller than that of a pure K vapor cell at the same temperature, are uniformly polarized along the Z-axis by circularly polarized D1 resonant light of K. The high-density Rb atoms are polarized through spin-exchange collisions with the polarized K atoms to obtain uniform polarization [6]. The main difference between the K-Rb SERF magnetometer and the K-Rb-21Ne SERF comagnetometer is that the 4He buffer gas is replaced by 21Ne gas with a nonzero nuclear magnetic moment. Furthermore, 21Ne nuclear spins are polarized along the Z-axis by spin-exchange collisions with alkali atoms. When the frequency of the circularly polarized pump light is detuned away from the resonance frequency of the K D1 line, light shifts for K and Rb atoms occur, which are referred to as LK and LRb, respectively [20].

The properties of the comagnetometer are described by the Bloch equations [4,20], coupling K electron, Rb electron, 21Ne nuclear polarizations PK, PRb, and Pn with spin-exchange interactions - see Eqs. (1)(3). The K-Rb magnetometer can be described by Eqs. (1) and (2) with the parameters related to 21Ne atoms omitted [10],

PKt=γeQ(PK)(B+LK+λKNeM0nPn)×PK+RpSp+RseRbK(PRbPK)Q(PK)PK{T1K,T2K,T2K},
PRbt=γeQ(PRb)(B+LRb+λRbNeM0nPn)×PRb+RseKRb(PKPRb)Q(PRb)PRb{T1Rb,T2Rb,T2Rb},
Pnt=γn(B+λRbNeM0RbPRb)×PnΩ×Pn+RseRbNe(PRbPn)Pn{T1n,T2n,T2n}.
Here γe, and γn are the gyromagnetic ratios of electron and 21Ne nucleon, while Q(PK) and Q(PRb) are the slow-down factors [21]. B, LK/Rb and Ω are magnetic field, light shifts, and inertia rotation, respectively. The spin-exchange interaction between alkali and 21Ne atoms is represented by an effective magnetic field, λM0P, which is experienced by one spin species due to the average magnetization of the other [9]. Here λ = 8πk/3, and k is the Fermi-contact-shift enhancement factor [9, 22], and M0 is the magnetization corresponding to full polarization, and P is the polarization. The effective magnetic fields of Rb-21Ne pair are denoted as Bee=λRbNeM0RbPRb and Bnn=λRbNeM0nPn. Rp and Sp are pumping rate and photon polarization of pump light, respectively. RseKRb=nKσKRbSEν¯KRb is the spin-exchange rate from K atoms to Rb atoms, and the other way around for RseRbK=nRbσKRbSEν¯KRb. Here, σKRbSE is the spin-exchange cross-section between K and Rb, and ν̄KRb is the relative velocity between K and Rb [23]. Thus, RseKRb=DrRseRbK. Similarly, RseRbNe is the spin-exchange rate from Rb atoms to 21Ne atoms. The longitudinal and transverse relaxation times of K atoms are T1K and T2K. For the Rb spin ensemble, the longitudinal and transverse relaxation times are T1Rb and T2Rb. T1n and T2n are the longitudinal and transverse relaxation times of 21Ne.

When the external magnetic field is relatively small compared to the effective spin-destruction rate and pump rate, the equilibrium polarization of alkali atoms along the Z-axis for both magnetometer and comagnetometer is given by [20]:

PK0ePRb0eDrRpDrRp+1/T1Rb.
Due to the extremely large spin-exchange cross-section between K-Rb, the K-Rb spin-exchange rate exceeds 106 1/s at typical densities. This means that the K and Rb atoms are in spin-temperature equilibrium and have the same polarization [6, 7], which is consistent with the derivation from the Bloch equations shown in Eq. (4).

To extract the excitation signals, the transverse polarization of Rb spins in both magnetometer and comagnetometer, which are induced by either a magnetic field or inertia rotation, are measured using linearly polarized probe light that is detuned towards the red side of the Rb D1 line using optical rotation [4]. When the input excitations are relatively small, the Bloch equations can be linearized, and the quasi-static signals of the magnetometer and comagnetometer are given by [9,20]:

PRbxmag=PRb0eByRtotRb/γe+BxBzRtotRb2/γe2+Bx2+By2+(Bz+Lz)2,
SRbxcomag=kPRbxcomag=kPRb0eRtotRb/γeRtotRb2/γe2+(Lz+δBz)2×(Ωyγn+Ωxγe(Lz+δBz)γnRtotRbByδBzBnn).
Here,
Lz=LzRb+DrLzK,
LzK=recfD1KΦSpγeννD1K(ννD1K)2+(ΓD1K/2),
LzRb=recΦSpγe(fD1Rb1ννD1Rb+12fD2Rb1ννD2Rb),
where δBz=BzBzeeBznn, RtotRb is the transverse relaxation rate of Rb, and k is the scale factor. There are also several constants, re, c, fD1K, fD1Rb, fD2Rb, representing the classical electron radius, the speed of light, and the oscillator strengths of K D1 line, Rb D1 and D2 lines respectively. The photons flux per area per time and frequency of pump light are Φ and ν respectively. Besides, νD1K, νD1Rb, and νD2Rb are the D1 and D2 resonant frequencies of K and Rb, and ΓD1K is the D1 linewidth of K (full-width-at-half-maximum).

According to Eq. (7), the effective light shift experienced by the Rb atoms is a mixture of the light shifts of Rb and K [20]. This can be explained by the rapid spin-exchange collision between K and Rb atoms - see Fig. 1. Initially, the transverse polarizations PxK and PxRb are aligned, and the total transverse polarization is Px. When applying the light shifts LzK and LzRb, PxK and PxRb are rotated by their own light shift during the short period between spin-exchange collision, which leads to PxK and PxRb respectively. Then, the spin-exchange collision can re-orient PxK and PxRb together in a new direction P′x. Therefore, the precession of one spin species around their own light shift can be transferred to the other due to rapid spin-exchange collision. Thus, the K and Rb atoms can experience a mixture of the light shifts that is related to the density ratio.

 figure: Fig. 1

Fig. 1 Vector diagram illustrating the mixture of light shifts due to rapid spin-exchange collision.

Download Full Size | PDF

Near the D1 line of K, the light shift of Rb, LzRb, is nearly independent of the light frequency. The direction and amplitude of the light shift of K, LzK, however, changes significantly with the light frequency. Thus, the effective light shift Lz can become zero by changing LzK to compensate LzRb. According to Eq. (7), at this Lz = 0 point, the density ratio Dr is given by

Dr=LzRbLzK=1.025νD2RbνD1Rb(ννD1Rb)(ννD2Rb)(ννD1K)2+(ΓD1K/2)2ννD1K,
where the resonant frequencies νD1K, νD1Rb, and νD2Rb, and the linewidth ΓD1K can either be calculated using the composition of the cell or derived from the measured absorption curves. After the pump light frequency ν for Lz = 0 was determined, the density ratio of K to Rb, Dr, can be calculated directly by substituting ν into Eq. (10). By scanning the pump light frequency and measuring the corresponding effective light shift, Lz, the frequency ν for Lz = 0 can be obtained.

Moreover, under typical experimental conditions, the frequency shifts of the D1 and D2 resonance frequencies of Rb (νD1Rb,νD2Rb) from the resonance frequencies in vacuum (ν0D1Rb, ν0D2Rb) are usually on the order of several GHz, and the pump light frequency ν is usually detuned several GHz away from the D1 resonance frequency of K in vacuum (ν0D1K). Because these are negligible compared to the frequency differences ννD1Rb and ννD2Rb (on the order of several THz), ννD1Rb and ννD2Rb in Eq. (10) can be approximated by ν0D1Kν0D1Rb and ν0D1Kν0D2Rb respectively. In addition, both νD1K and ΓD1K can be extracted by fitting the dispersion curve for the measured effective light shift with Eqs. (7)(9). Therefore, we only need to scan the pump light frequency ν around the resonance frequency and measure the corresponding Lz to find the ν for Lz = 0. We then fit the measured Lz with Eqs. (7)(9) to obtain νD1K and ΓD1K, and substitute the frequency detune δν=ννD1K and linewidth ΓD1K into Eq. (10) to determine the density ratio Dr.

The effective light shift in the K-Rb magnetometer can be measured using the cross-modulation method, which is usually applied to zero residual magnetic field along the Z-axis [23]. When a small field modulation is applied along the X-axis, with a slow modulation frequency, the response to the field can be simplified due to the small residual magnetic field inside the magnetic shield :

PRbxmagPRb0eBxmodsin(ωxt)(Bz+Bzres+Lz)γe/RtotRbRtotRb2/γe2+(Bz+Bzres+Lz)2,
where Bxmod and ωx are the modulation amplitude and frequency, and Bzres is the residual magnetic field inside the magnetic shielding. As indicated by Eq. (11), the Bz+Bzres+Lz=0 point can be found by adjusting Bz to zero the oscillation of the magnetometer signal at that modulation frequency. Because Bzres is constant and Lz is proportional to the light intensity Iz with light frequency fixed, the Bzres can be determined by changing the light intensity. After the determination of Bzres, the effective light shift Lz=BzBzres at different pumping light frequencies and intensities can be measured based on the cross-modulation method.

The effective light shift in the K-Rb-21Ne comagnetometer can be measured using the Bz magnetic-field-zeroing method [24]. According to Eq. (6), the response of a By square wave modulation is given by:

ΔSxe=kΔPxe=kPzeγeRtotRbRtotRb2/γe2(Lz+δBz)2δBzBnnΔBy,
where ΔBy is the amplitude of the square wave. The response ΔSxe for different δBz forms a Lorentzian-shape curve, and Lz determines the center of the curve. Thus, Lz can be measured by fitting the response ΔSxe for different δBz with Eq. (12).

To further verify the accuracy of this measurement method, the density ratios measured with this method for different conditions were compared to those measured using the conventional LAS method and Raoult’s law. According to Raoult’s law, the density of K is nK = fKnK0, where nK0 is the saturated vapor density for pure K alkali metal, and fK is the mole fraction of K in the mixed alkali metals [6]. Because at a typical high temperature, around 463 K, the LAS method becomes inaccurate due to strong light-absorption, it cannot be used to measure the high density at high temperatures directly. However, the LAS method can be used to measure the K density at lower temperature to determine the mole fraction of K, fK, and then to derive the density of K at high temperature using Raoult’s law. The basic principle of the LAS method is briefly described below [23]. For linear polarized light, transmitting through alkali vapor, the optical depth OD is given as:

OD=ln(I/I0),
OD=nLrecfΓL/2(νν0)2+(ΓL/2)2.
Here, I0 and I are the incident and transmitted light intensities respectively, n is the alkali metal density, L is the diameter of the cell, re is the classical electron radius, c the speed of light, f the oscillator strength, ΓL the linewidth, ν0 the resonance frequency, and ν the light frequency. By scanning the light frequency ν around ν0 and measuring I0 and I, the extracted OD = − ln(I/I0) can be fitted with Eq. (14) to obtain the density.

3. Experimental setup and results

The experimental setup is shown in Fig. 2, which is similar to the one used in our earlier studies [25, 26]. The spherical cell was heated by a 110 kHz AC electrical heater, which was situated within a boron nitride ceramic oven and installed in a PEEK vacuum vessel. The vessel, cooled by the water-cooling tube, was enclosed by a ferrite barrel and a 5-layer cylindrical μ-metal shield to shield ambient magnetic fields and reduce inner magnetic noise. In addition, a three-axis coil was used to compensate the residual magnetic fields. The pump light, produced by an external cavity diode laser, was tuned to the D1 resonance of K, amplified by a tapered amplifier to 1.5 W (Toptica TApro), and expanded to cover the cell to polarize K atoms along the Z-axis. The transverse component of the Rb spin polarization was measured by a linear-polarized probe light, propagating along X-axis from a distributed feedback (DFB) laser, whose wavelength was red-detuned (0.4 nm) from the D1 resonance of Rb. The probe light was modulated by a photo-elastic modulator (Hinds Instruments PEM100) with the amplitude of about 0.08 rad and a frequency of about 50 kHz. It was demodulated with a lock-in amplifier (Zurich Instruments HF2LI). The wavelength of the light was monitored by a wavelength meter with a measurement accuracy of 60 MHz (High Finesse WS7). In this work, four different cells were measured: two spherical cells with a diameter of 25 mm, which were filled with the same 2660 Torr 4He, 50 Torr N2 gases but with a different K-Rb mole fraction to be used in the magnetometer (referred to as cells Mag1 and Mag2). Another two cells with a diameter of 12 mm were filled with the same 2665 Torr 21Ne, 52 Torr N2 gases but with a different K-Rb mole fraction to be used in the comagnetometer (referred to as cells Comag1 and Comag2).

 figure: Fig. 2

Fig. 2 Schematic of the experimental setup. TA, tapered amplifier; BE, beam expander; GT, Glan-Thompson polarizer; PD, photo detector; PBS, polarization beam splitter; PEM, photo elastic modulator.

Download Full Size | PDF

Firstly, we simulated the light shifts in a K-Rb SERF magnetometer. At typical experimental conditions, the cell is filled with 2660 Torr 4He and 50 Torr N2. Thus, the frequency shifts of the K D1 line, the Rb D1 line and the D2 line are 10.67 GHz, 9.06 GHz, and 0.51 GHz respectively. The linewidth of the K D1 line, due to pressure broadening, is 41.28 GHz. In Fig. 3, the light shift for Rb with a pump light power of 200 mW/cm2 is calculated using Eq. (9) and plotted on the left-Y-axis-bottom-X-axis coordinates (black curve). The zero point of the bottom-X-axis is the D1 line of K in vacuum, 770.108 nm. Around the D1 line of K, the LzRb is nearly constant, with a variation below 1%. According to Eqs. (7)(9), the effective light shifts Lz for a density ratio of Dr = 0.01 are simulated for different pump-light powers and drawn using the left-Y-axis-bottom-X-axis coordinates (red, blue and pink curves). These effective light shift curves intersect at Lz = 0 with a detune frequency of δν=ννD1K=5.43GHz. For different density ratios, the frequency detune δν for Lz = 0 is different, and their relationship can be simulated using Eq. (10) and plotted using the right-Y-axis-top-X-axis coordinates (green curve). The smaller the density ratio is, the larger is the frequency detune δν for Lz = 0. Thus, it can be used to measure the density ratio directly.

 figure: Fig. 3

Fig. 3 Simulation of light shifts in the K-Rb magnetometer under different conditions.

Download Full Size | PDF

Next, we measured the density ratio using this method in the K-Rb SERF magnetometer. To study the validity of this method for different pump light intensities, the density ratios of the cell Mag1 with different pump light intensities at 463 K were measured. Using the cross-modulation method from above, the residual magnetic field Bzres is the intercept of the line Bz=BzreskIz with the vertical coordinate, 0.73 nT - see the inset in Fig. 4. Lz at different pump light frequencies and incident intensities were measured - see Fig. 4. After fitting the measured Lz with Eqs. (7)(9), we can determine the light shifts for Rb LzRb, the linewidths ΓD1K and the frequency detunes δν for Lz = 0 point. The fitted light shifts for Rb LzRb are shown in Fig. 5(a) as red squares. They can be fitted with LzRb=0.0039Iz+0.00014 (blue solid line), as predicted by Eq. (9). The slope of the fitted line is slightly smaller than the theoretically calculated line LzRb=0.0072Iz (black dashed line). This observed difference is due to the absorption of pump light along the propagation direction and inhomogeneity of the light intensity. The fitted linewidths ΓD1K and the frequency detunings δν for Lz = 0 are shown in Fig. 5(b) as green circular and blue triangular symbols respectively. The fitted linewidths are larger than the theoretically calculated pressure-broadened linewidth (47.9 GHz) using the filled gas pressure. This can be explained by the power broadening of the high light intensity [20,27]. Substituting the measured ΓD1K and δν into Eq. (10), the density ratios Dr at these light intensities can be calculated respectively - see the red squares in Fig. 5(b). The mean value and the standard deviation of these measured density ratios at different light intensities are 0.0436 and 0.0033, respectively. That means the measured Dr at different light intensities agree with each other, thus this method is valid for different pump light intensities.

 figure: Fig. 4

Fig. 4 The effective light shift as a function of pump light frequency detuning. The point data were obtained at different light intensities, while the corresponding curves are fitted lines. The zero point of the horizontal axis is the D1 line of K in vacuum, 770.108 nm. Inset: The compensating magnetic field Bz as a function of light intensity Iz.

Download Full Size | PDF

 figure: Fig. 5

Fig. 5 (a) The light shift of Rb LzRb as a function of light intensity, (b) Fitted linewidths and the frequency detunes for Lz = 0, plotted against the left Y-axis, while the plotted frequency detunes are scaled up 10-fold for clarity. The calculated density ratios are plotted against the right Y-axis.

Download Full Size | PDF

To further verify this method, the density ratios for the cell Mag1 at different temperatures and identical incident pump light intensity were measured and compared with the results of the LAS method. The mole fraction of K in the cell Mag1 was calibrated using a pure K cell (filled with the same gases) by the LAS method. As shown in the inset of Fig. 6, the density of K, nKp, in a pure K cell was measured to be 1.217×1013 cm−3 using the LAS method with a temperature set to 423 K. This value approaches the saturated vapor density of 1.222×1013 cm−3 at 423 K, which indicates that the temperature of the inner cell is consistent with the set value. At 443 K, the mole fraction for K in the cell Mag1 fK=nK/nKp was measured to be 0.204 ± 0.021. Thus, the density ratios at higher temperatures could be calculated using Raoult’s law with an estimation error about 7.4% - see the black circles and shaded area in Fig. 7(a). With the same pump light intensity (166 mW/cm2) but different temperatures (443 K, 453 K, 463 K, and 473 K), the density ratios were measured using the new method - see Fig. 6. These measured effective light shifts, Lz, were fitted using Eq. (7)(9). The Lz of the higher-temperature cell is smaller than that of the lower-temperature cell, because the Lz is proportional to the pump light intensity according to Eqs. (7)(9) and the pump light intensity in the probe area of the higher-temperature cell is smaller due to the stronger absorption of the identical incident pump light at higher-temperature. The fitted linewidths ΓD1K and frequency detunes δν for Lz = 0 are shown in Fig. 7(a) using blue and green triangles, respectively. The linewidths at different temperatures are consistent with each other, which suggests that the linewidths are not significantly dependent on temperature [28]. Using Eq. (10), the density ratios at these temperatures were calculated and drawn in Fig. 7(a) with red squares. The density ratios measured with our method (red squares) are consistent with those using the LAS method (black circles). The mean relative error of these two methods is within 10%. The density ratios for these two methods all increase with temperature as predicted by Raoult’s law. The density ratios of our method are slightly larger than those of the LAS method, which is attributed to the heating of the cell by the strong pump light under the normal operation of K-Rb magnetometer. The measurement of density ratio by our method was performed under the normal operation of K-Rb magnetometer (with strong pump light), while the measurement of density ratio by the LAS method was performed with a weak linearly polarized light, thus the error caused by the heating-effect of the strong pump light under the normal operation of K-Rb magnetometer was considered in our method compared to the LAS method. To verify the heating effect, we monitored the absorption of the probe light when the pump light was turned on or off. The intensity of the transmitted probe light decreased significantly within several minutes when the pump light was turned on and slowly heated the cell.

 figure: Fig. 6

Fig. 6 The effective light shift as a function of the pump light frequency detuning. The data points were measured at different temperatures, while the corresponding curves are fitted lines. The inset shows the laser absorption spectra of the pure K vapor cell (black triangles) and the cell Mag1 (red squares). The zero point of the horizontal axis is the D2 line of K in vacuum 766.701 nm, while the vertical axis represents the log base e of the transmitted light intensity. The black triangles and red squares are fitted with y = 4119/[(x − 0.3039)2 + 30.722] − 1.811 and y = 1315/[(x + 0.0709)2 + 31.252] − 1.796 based on Eq. (14).

Download Full Size | PDF

 figure: Fig. 7

Fig. 7 Density ratios at different temperatures plotted against the left Y-axis. The fitted linewidths and the frequency detunes for Lz = 0 at different temperatures plotted against the right Y-axis. The shaded areas are the estimation errors of the calculated density ratios. (a): Data for the cell Mag1. The plotted frequency detunes are scaled up 10-fold for clarity. (b) Data for the cell Mag2. The plotted frequency detunes are scaled by a factor of five for clarity.

Download Full Size | PDF

Using the K-Rb SERF magnetometer, another K-Rb cell (cell Mag2) with a smaller mole fraction of K was also used to verify this method. Likewise, the density ratios of the cell Mag2 were measured at different temperatures (443 K, 453 K 463 K and 473 K) with identical pump light intensity using our method. They were then compared with the results of the LAS method. The mole fraction of K in the cell Mag2 was calibrated to be 0.076 ± 0.057 at 423 K, using the LAS method, which is nearly three times smaller than that of the cell Mag1. Using Raoult’s law, the density ratios were calculated at these temperatures with an estimation error about 6.1% - see the black circles and shaded area in Fig. 7(b). Using our method, the effective light shifts Lz were measured at these temperatures and fitted with Eqs. (7)(9). The fitted linewidths and frequency detunes for Lz = 0 are shown in Fig. 7(b) with blue and green triangles. The density ratios calculated with our method for these temperatures are marked as red squares in Fig. 7(b). There is a good agreement between the results of these two methods, with a mean relative error below 5%.

We also verified this method with the K-Rb-21Ne comagnetometer. The effective light shift was measured using the Bz magnetic field zeroing method from above: Scanning δBz, the response signals ΔSxe of a By square wave modulation were measured. The measured ΔSxe were fitted with Eq. (12) to acquire Lz. The measured ΔSxe at different pump light frequencies for cell Comag1 are shown in Fig. 8, and fitted using Eq. (12) to obtain Lz. Similarly, the light shift for Rb, LzRb, the linewidth ΓD1K, and frequency detune δν for Lz = 0 could be acquired by fitting the measured Lz with Eqs. (7)(9) - see Fig. 9(a). The fitted light shift of Rb LzRb is 1.75 nT, which is consistent with the theoretically calculated result (1.77 nT) using pump light of 242 mW/cm2. The fitted linewidth is ΓD1K=34.8GHz, larger than the theoretically calculated pressure broadened linewidth of K, 26.9 GHz. This can be explained by power broadening of the strong pump light, similar to the K-Rb magnetometer. The frequency detune for Lz = 0 is δν = 1.9 GHz. Substituting ΓD1K and δν into Eq. (10), the density ratio is 0.0191 with an operating temperature of 463 K. Similarly, to verify the accuracy of this method, the measured density ratio with our method was compared to the LAS method. As shown in the inset of Fig. 9(a), the density of the K atoms was measured at a lower temperature (443 K) with the LAS method, resulting in the mole fraction fK = 0.095 ± 0.009 for K. Thus, the density ratio is 0.0166 at 463 K with an estimation error about 8.6% according to Raoult’s law. This value is consistent with the result obtained with our method with a relative error of 15.1%. This difference is due to the optical heating-effect and the measurement error of the mole fraction of K both.

 figure: Fig. 8

Fig. 8 The response signals ΔSxe of a By square wave modulation for different δBz and pump light frequencies. The dots are measured signals, while the solid lines are fitted lines.

Download Full Size | PDF

 figure: Fig. 9

Fig. 9 Effective light shift Lz as a function of the pump light frequency detuning. The zero point of the horizontal axis is the K D1 line in vacuum, 770.108 nm. Inset: The absorption curve of K atoms in the cell at 443 K. The zero point of the horizontal axis is the D2 line of K in vacuum (766.701 nm), while the vertical axis is the log base e of the transmitted light intensity. (a) Data for cell Comag1. The black dots in the inset are fitted using y = 284.7/[(x +7.411)2 +14.522]−1.597 based on Eq. (14). (b) Data for cell Comag2. The black dots in the inset are fitted using y = 195.9/[(x +6.511)2 + 16.942] − 1.425 based on Eq. (14).

Download Full Size | PDF

To further verify the accuracy of our method, another K-Rb-21Ne spherical cell (cell Comag2) with a smaller mole fraction of K, was measured with our method and compared with the LAS method. As shown in the inset of Fig. 9(b), the mole fraction of K was calibrated to be fK = 0.034 ± 0.002 using the LAS method at 453 K. This value is about 2.8 times smaller than that of the cell Comag1. The corresponding density ratio at 463 K is Dr = 0.0056 with an estimation error about 5.7% using Raoult’s law. The effective light shift was also measured with our method - see Fig. 9(b). The fitted linewidth is 29.13 GHz, and the fitted light shift of Rb LzRb is 1.76 nT. These values are consistent with the theoretically calculated results. The frequency detune for Lz = 0 is 4.43 GHz, which means the density ratio is Dr = 0.0062 at 463 K. This is consistent with the result measured using the LAS method. The relative error between these two results is 10.7%.

As a further check, we measured the heating-effect of the strong light with different intensities by the LAS method. With the electric heater temperature set to 443 K, the absorptions of the weak linearly polarized the K D2 light (along X-axis) in a K-Rb-21Ne cell were measured with the strong linearly polarized Rb D1 light (along Z-axis) turned off or on in Fig. 10(a). Fitting these absorption curves with Eq. (14), we obtained the densities of K. And then the densities were used to derive the corresponding true in-situ temperatures based on the Rauolt’s law and the saturated vapor density. The temperature increments between the true in-situ temperature and the setting temperature of electric heater are plotted on the left-Y-axis-bottom-X-axis in Fig. 10(b), which increase significantly with the increasing of light intensities. Optical heating via the absorption of light has been studied previously [29,30]. A 200 mW laser light (λ = 915 nm) was used to heat the cell to 97 °C [29], while a cell temperature of 150 °C was achieved with a 140 mW laser light (λ = 1.5 um) absorbed by absorption filters attached to the cell windows [30]. Compared with these previous works, the optical heating effect in our experiment is much smaller. This is because the optical heating efficiency strictly depends on the size and absorption element of the cell, and thermal isolation [29, 30]. The relationship between the saturated vapor density ratio Dr0 = nK0/nRb0 and temperature is simulated and plotted on the right-Y-axis-top-X-axis in Fig. 10(b). In our experiment, the pump light of about 100 mW/cm2 − 700 mW/cm2 causes an increment of cell temperature about 5 − 10 °C resulting in an increment of the saturated vapor density ratio smaller than 4%. Therefore, it would cause an increment of the density ratio smaller than 4% according to DrDr0fK/(1 − fK).

 figure: Fig. 10

Fig. 10 Measurement of the optical heating effect. (a). The absorption curves of K atoms at 443 K with different light intensities. The zero point of the horizontal axis is the D2 line of K in vacuum (766.701 nm), while the vertical axis is the log base e of the transmitted light intensity. The dots are measured data, and solid curves are corresponding fitting curves based on Eq. (14). (b). The cell temperature increments derived from the absorption curves is plotted on the left-Y-axis-bottom-X-axis. The simulated relationship between the saturated vapor density ratio Dr0 and temperature is plotted on the right-Y-axis-top-X-axis.

Download Full Size | PDF

4. Discussions and conclusion

We investigated a new method to directly measure the density ratio of an alkali metal mixture. The mixture of the light shifts was formulated using Bloch equations and explained by the fast spin-exchange interaction. The relationship between the density ratio and the frequency detune of Lz = 0 from the resonance was both formulated and simulated. We described and executed the procedure to measure the effective light shifts, and we found the Lz = 0 point for both the K-Rb SERF magnetometer and the K-Rb-21Ne SERF comagnetometer. The density ratios of different K-Rb magnetometer and K-Rb-21Ne comagnetometer cells were measured using our method for different temperatures, pump light powers, and mole fractions of K. We compared these with the results measured using the conventional LAS method. The measured density ratios of these two methods are consistent. When operating at different temperatures and mole fractions of K, the density ratios measured with our method are all slightly larger than the results measured by the LAS method, which is due to the heating of the cell by the strong pump light and the measurement errors of fK both.

This method is limited, when the density ratio of K to Rb is too small, because when the density of K is too small the light shift of K cannot compensate the light shift of Rb, which means the Lz = 0 point cannot be found. This limitation depends on the gas composition in the cell, and the minimum density ratio that can be measured is about 1/500 for typical operating conditions of the K-Rb SERF magnetometer and the K-Rb-21Ne SERF comagnetometer. However, this problem can be well solved by pumping Rb instead of K. Thus the light shift of Rb can be changed to compensate that of K to find the Lz = 0 point.

This method can be used to directly measure the density ratio of optical-thick hybrid alkali atoms with no need for additional strong magnetic field and absolute polarization measurement, which are necessary for magnetic-field-induced and polarization-induced Faraday-rotation methods respectively. With the advantage of in-situ measuring the density ratio under exactly normal operations, the errors caused by the strong pump light heating-effect and the long-term temperature drift can be real-time monitored by this method. This method can also be used for K-Cs and Rb-Cs hybrid cells. Precise determination of the density ratio is useful to find the optimal density ratio to decrease the gradient of polarization, increase the sensitive area, and improve the sensitivity of both magnetometer and comagnetometer, as well as to increase the pumping efficiency to hyperpolarize noble gases.

Funding

National Natural Science Foundation of China (NSFC) (61227902, 61374210); National Key R&D Program of China (2016YFB0501601).

References

1. D. Budker and M. Romalis, “Optical magnetometry,” Nat. Phys. 3, 227–234 (2007). [CrossRef]  

2. J. Lee, A. Almasi, and M. Romalis, “Improved limits on spin-mass interactions,” Phys. Rev. Lett. 120, 161801 (2018). [CrossRef]   [PubMed]  

3. E. Boto, N. Holmes, J. Leggett, G. Roberts, V. Shah, S. Meyer, L. D. Muñoz, K. Mullinger, T. Tierney, S. Bestmann, G. Barnes, R. Bowtell, and M. Brookes, “Moving magnetoencephalography towards real-world applications with a wearable system,” Nature 555, 657–661 (2018). [CrossRef]   [PubMed]  

4. R. Li, W. Fan, L. Jiang, L. Duan, W. Quan, and J. Fang, “Rotation sensing using a K − Rb −21 Ne comagnetometer,” Phys. Rev. A 94, 032109 (2016). [CrossRef]  

5. M. Romalis, “Hybrid optical pumping of optically dense alkali-metal vapor without quenching gas,” Phys. Rev. Lett. 105, 243001 (2010). [CrossRef]  

6. E. Babcock, I. Nelson, S. Kadlecek, B. Driehuys, L. W. Anderson, F. W. Hersman, and T. G. Walker, “Hybrid spin-exchange optical pumping of 3He,” Phys. Rev. Lett. 91, 123003 (2003). [CrossRef]  

7. W. C. Chen, T. R. Gentile, and T. G. Walker, “Spin-exchange optical pumping of 3He with Rb − K mixtures and pure K,” Phys. Rev. A 75, 013416 (2007). [CrossRef]  

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

9. T. W. Kornack, R. K. Ghosh, and M. Romalis, “Nuclear spin gyroscope based on an atomic comagnetometer,” Phys. Rev. Lett. 95, 230801 (2005). [CrossRef]   [PubMed]  

10. Y. Ito, D. Sato, K. Kamada, and T. Kobayashi, “Optimal densities of alkali metal atoms in an optically pumped K − Rb hybrid atomic magnetometer considering the spatial distribution of spin polarization,” Opt. Express 24, 15391–15402 (2016). [CrossRef]   [PubMed]  

11. K. Nishi, Y. Ito, and T. Kobayashi, “High-sensitivity multi-channel probe beam detector towards meg measurements of small animals with an optically pumped K − Rb hybrid magnetometer,” Opt. Express 26, 1988–1996 (2018). [CrossRef]   [PubMed]  

12. M. W. Millard, P. P. Yaney, B. N. Ganguly, and C. A. DeJoseph, “Diode laser absorption measurements of metastable helium in glow discharges,” Plasma Sources Sci. Technol. 7, 389–394 (1998). [CrossRef]  

13. E. Vliegen, S. Kadlecek, L. W. Anderson, T. G. Walker, C. J. Erickson, and W. Happer, “Faraday rotation density measurements of optically thick alkali metal vapors,” Nucl. Instruments Methods Phys. Res. A 460, 444–450 (2001). [CrossRef]  

14. H. Zhang, S. Zou, X. Chen, M. Ding, G. Shan, Z. Hu, and W. Quan, “On-site monitoring of atomic density number for an all-optical atomic magnetometer based on atomic spin exchange relaxation,” Opt. Express 24, 17234–17241 (2016). [CrossRef]   [PubMed]  

15. H. Yao, H. Zhang, D. Ma, J. Zhao, and M. Ding, “In situ determination of alkali metal density using phase-frequency analysis on atomic magnetometers,” Meas. Sci. Technol. 29, 6 (2018). [CrossRef]  

16. B. Chann, E. Babcock, L. W. Anderson, and T. G. Walker, “Measurements of 3He spin-exchange rates,” Phys. Rev. A 66, 032703 (2002). [CrossRef]  

17. E. D. Babcock, “Spin exchange optical pumping with alkali metal vapors,” Ph.D. thesis, University of WisconsinMadison (2008).

18. B. S. Mathur, H. Tang, and W. Happer, “Light shifts in the alkali atoms,” Phys. Rev. 171, 11–19 (1968). [CrossRef]  

19. E. Zhivun, A. Wickenbrock, B. Patton, and D. Budker, “Alkali-vapor magnetic resonance driven by fictitious radiofrequency fields,” Appl. Phys. Lett. 105, 192406 (2014). [CrossRef]  

20. Y. Chen, W. Quan, L. Duan, Y. Lu, L. Jiang, and J. Fang, “Spin-exchange collision mixing of the K and Rb ac stark shifts,” Phys. Rev. A 94, 052705 (2016). [CrossRef]  

21. I. M. Savukov and M. Romalis, “Effects of spin-exchange collisions in a high-density alkali-metal vapor in low magnetic fields,” Phys. Rev. A 71, 023405 (2005). [CrossRef]  

22. R. K. Ghosh and M. V. Romalis, “Measurement of spin-exchange and relaxation parameters for polarizing 21Ne with K and Rb,” Phys. Rev. A 81, 043415 (2010). [CrossRef]  

23. S. J. Seltzer, “Developments in alkali-metal atomic magnetometry,” Ph.D. thesis, Princeton University (2008).

24. J. M. Brown, “A new limit on Lorentz- and CPT-violating neutron spin interactions using a K −3 He comagnetometer,” Ph.D. thesis, Princeton University (2011).

25. W. Quan, K. Wei, and H. R. Li, “Precision measurement of magnetic field based on the transient process in a K − Rb −21 Ne comagnetometer,” Opt. Express 25, 8470–8483 (2017). [CrossRef]   [PubMed]  

26. W. Ji, Y. Chen, C. Fu, M. Ding, J. Fang, Z. Xiao, K. Wei, and H. Yan, “New experimental limits on exotic spin-spin-velocity-dependent interactions by using smco5 spin sources,” Phys. Rev. Lett. 121, 261803 (2018). [CrossRef]  

27. M. L. Citron, H. R. Gray, C. W. Gabel, and C. R. Stroud, “Experimental study of power broadening in a two-level atom,” Phys. Rev. A 16, 1507 (1977). [CrossRef]  

28. R. Li, Y. Li, L. Jiang, W. Quan, M. Ding, and J. Fang, “Pressure broadening and shift of K D1 and D2 lines in the presence of 3He and 21Ne,” Eur. Phys. J. D 70, 139 (2016). [CrossRef]  

29. J. Preusser, V. Gerginov, S. Knappe, and J. Kitching, “A microfabricated photonic magnetometer,” Proc. IEEE Sensors Conf. (IEEE, 2008) pp. 344–346.

30. R. Mhaskar, S. Knappe, and J. Kitching, “A low-power, high-sensitivity micromachined optical magnetometer,” Appl. Phys. Lett. 101, 241105 (2012). [CrossRef]  

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 Vector diagram illustrating the mixture of light shifts due to rapid spin-exchange collision.
Fig. 2
Fig. 2 Schematic of the experimental setup. TA, tapered amplifier; BE, beam expander; GT, Glan-Thompson polarizer; PD, photo detector; PBS, polarization beam splitter; PEM, photo elastic modulator.
Fig. 3
Fig. 3 Simulation of light shifts in the K-Rb magnetometer under different conditions.
Fig. 4
Fig. 4 The effective light shift as a function of pump light frequency detuning. The point data were obtained at different light intensities, while the corresponding curves are fitted lines. The zero point of the horizontal axis is the D1 line of K in vacuum, 770.108 nm. Inset: The compensating magnetic field Bz as a function of light intensity Iz.
Fig. 5
Fig. 5 (a) The light shift of Rb L z Rb as a function of light intensity, (b) Fitted linewidths and the frequency detunes for Lz = 0, plotted against the left Y-axis, while the plotted frequency detunes are scaled up 10-fold for clarity. The calculated density ratios are plotted against the right Y-axis.
Fig. 6
Fig. 6 The effective light shift as a function of the pump light frequency detuning. The data points were measured at different temperatures, while the corresponding curves are fitted lines. The inset shows the laser absorption spectra of the pure K vapor cell (black triangles) and the cell Mag1 (red squares). The zero point of the horizontal axis is the D2 line of K in vacuum 766.701 nm, while the vertical axis represents the log base e of the transmitted light intensity. The black triangles and red squares are fitted with y = 4119/[(x − 0.3039)2 + 30.722] − 1.811 and y = 1315/[(x + 0.0709)2 + 31.252] − 1.796 based on Eq. (14).
Fig. 7
Fig. 7 Density ratios at different temperatures plotted against the left Y-axis. The fitted linewidths and the frequency detunes for Lz = 0 at different temperatures plotted against the right Y-axis. The shaded areas are the estimation errors of the calculated density ratios. (a): Data for the cell Mag1. The plotted frequency detunes are scaled up 10-fold for clarity. (b) Data for the cell Mag2. The plotted frequency detunes are scaled by a factor of five for clarity.
Fig. 8
Fig. 8 The response signals Δ S x e of a By square wave modulation for different δBz and pump light frequencies. The dots are measured signals, while the solid lines are fitted lines.
Fig. 9
Fig. 9 Effective light shift Lz as a function of the pump light frequency detuning. The zero point of the horizontal axis is the K D1 line in vacuum, 770.108 nm. Inset: The absorption curve of K atoms in the cell at 443 K. The zero point of the horizontal axis is the D2 line of K in vacuum (766.701 nm), while the vertical axis is the log base e of the transmitted light intensity. (a) Data for cell Comag1. The black dots in the inset are fitted using y = 284.7/[(x +7.411)2 +14.522]−1.597 based on Eq. (14). (b) Data for cell Comag2. The black dots in the inset are fitted using y = 195.9/[(x +6.511)2 + 16.942] − 1.425 based on Eq. (14).
Fig. 10
Fig. 10 Measurement of the optical heating effect. (a). The absorption curves of K atoms at 443 K with different light intensities. The zero point of the horizontal axis is the D2 line of K in vacuum (766.701 nm), while the vertical axis is the log base e of the transmitted light intensity. The dots are measured data, and solid curves are corresponding fitting curves based on Eq. (14). (b). The cell temperature increments derived from the absorption curves is plotted on the left-Y-axis-bottom-X-axis. The simulated relationship between the saturated vapor density ratio Dr0 and temperature is plotted on the right-Y-axis-top-X-axis.

Equations (14)

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

P K t = γ e Q ( P K ) ( B + L K + λ K Ne M 0 n P n ) × P K + R p S p + R se Rb K ( P Rb P K ) Q ( P K ) P K { T 1 K , T 2 K , T 2 K } ,
P Rb t = γ e Q ( P Rb ) ( B + L Rb + λ Rb Ne M 0 n P n ) × P Rb + R se K Rb ( P K P Rb ) Q ( P Rb ) P Rb { T 1 Rb , T 2 Rb , T 2 Rb } ,
P n t = γ n ( B + λ Rb Ne M 0 Rb P Rb ) × P n Ω × P n + R se Rb Ne ( P Rb P n ) P n { T 1 n , T 2 n , T 2 n } .
P K 0 e P Rb 0 e D r R p D r R p + 1 / T 1 Rb .
P Rb x mag = P Rb 0 e B y R tot Rb / γ e + B x B z R tot Rb 2 / γ e 2 + B x 2 + B y 2 + ( B z + L z ) 2 ,
S Rb x comag = k P Rb x comag = k P Rb 0 e R tot Rb / γ e R tot Rb 2 / γ e 2 + ( L z + δ B z ) 2 × ( Ω y γ n + Ω x γ e ( L z + δ B z ) γ n R tot Rb B y δ B z B n n ) .
L z = L z Rb + D r L z K ,
L z K = r e c f D 1 K Φ S p γ e ν ν D 1 K ( ν ν D 1 K ) 2 + ( Γ D 1 K / 2 ) ,
L z Rb = r e c Φ S p γ e ( f D 1 Rb 1 ν ν D 1 Rb + 1 2 f D 2 Rb 1 ν ν D 2 Rb ) ,
D r = L z Rb L z K = 1.025 ν D 2 Rb ν D 1 Rb ( ν ν D 1 Rb ) ( ν ν D 2 Rb ) ( ν ν D 1 K ) 2 + ( Γ D 1 K / 2 ) 2 ν ν D 1 K ,
P Rb x mag P Rb 0 e B x mod sin ( ω x t ) ( B z + B z res + L z ) γ e / R tot Rb R tot Rb 2 / γ e 2 + ( B z + B z res + L z ) 2 ,
Δ S x e = k Δ P x e = k P z e γ e R tot Rb R tot Rb 2 / γ e 2 ( L z + δ B z ) 2 δ B z B n n Δ B y ,
OD = ln ( I / I 0 ) ,
OD = n L r e c f Γ L / 2 ( ν ν 0 ) 2 + ( Γ L / 2 ) 2 .
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.