Long spin-1/2 noble gas coherence times in mm-sized anodically bonded batch-fabricated 3He-129Xe-87Rb cells

M. E. Limes    N. Dural    M. V. Romalis Department of Physics, Princeton University, Princeton, New Jersey, 08544, USA    E. L. Foley    T. W. Kornack    A. Nelson    L. R. Grisham Twinleaf LLC, Princeton, New Jersey, 08544, USA
(October 21, 2024)
Abstract

As the only stable spin-1/2 noble gas isotopes, 3He and 129Xe are promising systems for inertial rotation sensing and searches for exotic spin couplings. Spin-1/2 noble gases have intrinsic coherence times on the order of hours to days, which allows for incredibly low frequency error of free-precession measurements. However relaxation in miniature cells is dominated by interactions with the cell wall, which limits the performance of a chip-scale sensor that uses noble gases. While 129Xe wall relaxation times have previously been limited to 10’s of seconds in mm-sized cells, we demonstrate the first anodically bonded batch-fabricated cells with dual 3He-129Xe isotopes and 87Rb in a 6 mm3 volume with 3He and 129Xe T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT coherence times of, respectively, 4 h and 300 s. We use these microfabricated cells in a dual noble-gas comagnetometer and discuss its limits.

pacs:
32.30.Dx, 06.30.Gv,39.90.+d

The research and commercial development of minature high-performance atomic sensors [1, 2, 3, 4, 5] has naturally driven interest in chip-scale systems for hyperpolarization of noble gas nuclear spins [6, 7]. Compact noble-gas NMR systems are used in a variety of studies, including searches of short-range spin-dependent forces [8], inertial rotation sensing [9, 10, 11, 12, 13], magnetometry [14], and microfluidic detection [15, 16]. These applications use vapor cells that contain the readily available 129Xe or 131Xe that are hyperpolarized by spin exchange with optically pumped alkali metals. Noble-gas nuclei are notable for having extraordinarily long T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT coherence times while having high densities, making the quantum projection frequency noise δω=1/T2Nt𝛿𝜔1subscript𝑇2𝑁𝑡\delta\omega=1/\sqrt{T_{2}Nt}italic_δ italic_ω = 1 / square-root start_ARG italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_N italic_t end_ARG well below pHz levels [17], with N𝑁Nitalic_N the number of nuclei sampled and measurement time t𝑡titalic_t. The potential of noble gases for rotation measurements has long been considered [18, 19, 20], but more than one simultaneous measurement is needed to discern rotations of the apparatus from changes in the local magnetic field. Thus the use of a comagnetometer that samples the same magnetic field in the same volume is advantageous for rotation measurements done in a compact and cost-effective package. In principle the spin-1/2 isotopes 129Xe or 3He [21, 22, 23] are the best candidates for systematic-free rotation measurements, as higher-spin isotopes like 131Xe or 21Ne are susceptible to additional frequency shifts and relaxation mechanisms that can be deleterious [24, 25, 26, 27, 28]. Detection of nuclear spins at a quantum shot noise level is difficult to achieve. A popular and sensitive detection method that takes advantage of the wavefunction overlap of the optically addressable electronic state of an alkali metal with a noble gas nuclei [29, 30]. While backaction and relaxation from the alkali on the noble gas must be mitigated [31, 32], this method of detection is convenient because Rb alkali vapor is already present in these systems to hyperpolarize 3He and 129Xe through spin-exchange optical pumping.

Spin-1/2 noble gases intrinsically decohere very slowly through intraspecies collisional processes that form transient or persistent dimers which allow for dipole-dipole and spin-rotation relaxation [33, 34, 35]. Various extrinsic factors that cause nuclear relaxation include diffusion through magnetic field gradients, spin exchange and spin destruction of the noble gases with warm alkali vapors, and relaxation due to the cell wall. Given a particular coating or cell surface, the characteristic wall relaxation time Twallsubscript𝑇wallT_{\text{wall}}italic_T start_POSTSUBSCRIPT wall end_POSTSUBSCRIPT scales linearly with surface-to-volume ratio and follows an Arrhenius temperature dependence. In small cells with a large surface-to-volume ratio, the noble-gas coherence times are often limited by interactions with the cell walls, with past studies showing 129Xe relaxation of 10’s of seconds [11]. We demonstrate batch fabrication of stemless anodically bonded cells that have a well-defined geometry and contain 3He and 129Xe with nuclear-spin coherence times of 4 hours and 300 s, respectively, a significant improvement in the 129Xe coherence time for micro-fabricated cells and the first successful detection of 3He signals in such cells. We then use these cells as a free-precession dual noble-gas comagnetometer in a laboratory-sized prototype for a chip-scale gyroscope, achieving a stability of 0.2 deg/h at 48 h.

Refer to caption
Figure 1: Drilled silicon wafers can be filled with a variety of noble gases and enriched alkali isotopes then anodically bonded and diced to make many uniform batch fabricated vapor cells. We fill cells with 3He-129Xe-87Rb-N2 for use in a dual noble-gas comagnetometer detected by an in situ alkali magnetometer and use long noble-gas relaxation times to study a low-drift chip-scale gyroscope. For detection of the noble gases, we use co-linear probe and pump beams, with pump helicity flipped by an electro-optical modulator (EOM) and probe shuttered with an acousto-optical modulator (AOM), along with a fast 87Rb π𝜋\piitalic_π pulse train applied with small Helmholtz coils.

For this study we make stemless cells shown in Fig. 1 by using a silicon wafer of thickness 2 mm, drilled with an array of 2 mm holes, machined with two 0.2 mm thick SD-2 aluminosilcate glass plates on each side [1, 36]. The wafer and glass are baked under high vacuum, which is crucial to obtain long 129Xe lifetimes. We first bond glass to one side of the wafer, and then introduce 99.9%percent99.999.9\%99.9 % isotopically enriched 87Rb metal into the cells. Then the other side of the wafer is bonded to enclose the cell under a 1400:6.5:80 torr 3He:129Xe:N2 gas mixture in the bonding chamber. The N2subscript𝑁2N_{2}italic_N start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is included for quenching during optical pumping of 87Rb. Heat applied during bonding can cause buffer gas pressure loss in cells of up to 30%percent\%%; we check the 3He and 87Rb densities by measuring the Rb absorption spectrum of the probe laser. A cryogenic storage system reclaims 3He from the bonding chamber, enabling a pathway for fabricating many low-cost devices. Anodically bonded glass on machined silicon also allows for excellent control of the cell geometry. Typical glass blown or optically contacted cells have a glass stem for cell filling and sealing [37, 10], which creates an inherent asymmetry in the cell shape. Some are able to lessen cell stem effects by plugging the stem with alkali metal [38] or a small piece of movable glass [39]. Others have made stemless cells by allowing 3He to diffuse through quartz walls at a high temperature [40]. The cell shape is an important consideration for precision measurements and chip-scale sensors that use noble-gas precession frequency measurements, as it affects the long-range dipolar magnetic interactions between spin-polarized nuclei [41]. Such interactions cause significant frequency shifts in precision comagnetometer measurements [42], where long-term drifts can be correlated to placement and magnitude of noble-gas polarization [43, 44, 45, 46, 47]. To study this dipolar effect, we recently used our cell fabrication process to make a slight wedge in the Si wafer to vary the height-to-diameter aspect ratios across many cells. We demonstrated the existance of an optimal aspect ratio that nulls the dipolar effect, as well as the existence of nuclear J𝐽Jitalic_J-coupling between different noble gases 3He-129Xe [44]. This J𝐽Jitalic_J-coupling in principle also provides a mechanism for relaxation and polarization transfer between the different noble gas species, although this has not yet been directly observed.

We include 87Rb in the cells in order to polarize nuclear spins by spin-exchange with 87Rb that is polarized by optical pumping using a circularly polarized laser beam parallel to a longitudinal bias magnetic field B0=0.5subscript𝐵00.5B_{0}=0.5italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 µT. We also detect noble-gas spin precession with a 87Rb magnetometer [48, 31, 49], which gives a high signal-to-noise ratio due to wavefunction overlap of the 87Rb with the noble gases [50, 51, 52, 53]. Since the cells have only a single optical axis, counter-propagating pump and probe beams are used with fast 87Rb magnetic field π𝜋\piitalic_π pulses to make a pulse-train magnetometer similar to that described in Ref. [31, 44]. For 87Rb magnetometer operation, the polarization of an 795 nm pump laser is alternated between σ+superscript𝜎\sigma^{+}italic_σ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and σsuperscript𝜎\sigma^{-}italic_σ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT light at 13 kHz by an EOM. Synchronous 2 μ𝜇\muitalic_μs long magnetic field pulses along the y𝑦yitalic_y-axis (or any axis in the transverse plane) quickly Rabi flop the 87Rb polarization back and forth along B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. If the pulses are sufficiently fast and short, they can suppress alkali-alkali spin-exchange relaxation as well as precession of 87Rb spins around B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [31]. The 87Rb polarization along B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is detected with a probe beam that is detuned off the D1𝐷1D1italic_D 1 resonance, passing through the cell to a balanced polarimeter. The probe laser is shuttered by an AOM to open only during π𝜋\piitalic_π pulses to mitigate probe broadening. A DC or low frequency Bysubscript𝐵𝑦B_{y}italic_B start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT field parallel to the direction of the fast pulses delays or advances the 87Rb signal zero-crossing during the π𝜋\piitalic_π Rabi flop. This signal is measured by a lock-in amplifier referenced to half the EOM frequency with the lock-in output proportional to Bysubscript𝐵𝑦B_{y}italic_B start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT. The sensitivity of the 87Rb magnetometer is 300 fT/HzHz\sqrt{\text{Hz}}square-root start_ARG Hz end_ARG. We use sub-mW power from Photodigm DBR lasers through the active volume, similar to the powers from commercial VCSEL laser packages.

Refer to caption
Figure 2: Top: Spin precession signal from 129Xe at 73.3C show relaxation times of T2=308subscript𝑇2308T_{2}=308italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 308 s. Bottom: The dependence of 129Xe T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT on the Rb number density. The inset shows that at 120superscript120120^{\circ}120 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPTC 87Rb pump light intensity and deviations of the alkali magnetometer π𝜋\piitalic_π pulse amplitude from optimal conditions shorten 129Xe T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT due to Xe diffusion in a 87Rb polarization gradient.

Although self-relaxation for the noble gases in our conditions yields coherence times of hundreds of hours for both 129Xe and 3He, our system is dominated by wall relaxation and Rb-noble gas interactions. To study 129Xe relaxation mechanisms, we first spin-exchange optically pump the 129Xe along the bias field for 30-50 s, and have negligible 3He polarization. A tipping pulse then places 129Xe transverse to B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and its precession is detected with the pulse-train 87Rb magnetometer. A representative signal from the lock-in amplifier is shown in Fig. 2a. We fit a given data set to the function Aexp(t/T2)sin(ωXet)𝐴𝑡subscript𝑇2subscript𝜔Xe𝑡A\exp(-t/T_{2})\sin(\omega_{\text{Xe}}t)italic_A roman_exp ( - italic_t / italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) roman_sin ( italic_ω start_POSTSUBSCRIPT Xe end_POSTSUBSCRIPT italic_t ) to extract T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, and show 129Xe 1/T21subscript𝑇21/T_{2}1 / italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT as a function of cell temperature and Rb density in Fig. 2b. In this cell we find the Rb absorption FWHM of 23.4 ±1plus-or-minus1\pm 1± 1 GHz corresponding to 1000100010001000 torr 3He [54]. From the slope in Fig. 2b we find a Rb-Xe spin-exchange rate of (7.8±0.7)×1016plus-or-minus7.80.7superscript1016(7.8\pm 0.7)\!\times\!10^{-16}( 7.8 ± 0.7 ) × 10 start_POSTSUPERSCRIPT - 16 end_POSTSUPERSCRIPT cm3/s. This value is in agreement with a calculation of the spin-exchange rate 7.2×10167.2superscript10167.2\times 10^{-16}7.2 × 10 start_POSTSUPERSCRIPT - 16 end_POSTSUPERSCRIPT cm3/s based on previously measured cross-sections [55, 56]. The relaxation rate does not change significantly from the Rb-Xe spin-exchange rate across our temperature range, so we neglect any Arrhenius temperature dependence in our analysis. The inset of Fig. 2b shows the Rb magnetometer can shorten 129Xe T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT by causing a Rb polarization gradient that 129Xe diffuses through. The additional relaxation depends on the 87Rb pump laser power and the accuracy of 87Rb π𝜋\piitalic_π pulses. Increasing the 87Rb π𝜋\piitalic_π pulse repetition rate decreases this relaxation because the Rb gradient is reversed more rapidly. For accurate T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT measurements we use proper π𝜋\piitalic_π pulses and low pump power. The 3He coherence exceeds 129Xe T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT by at least an order of magnitude, thus sufficient for our comagnetometer operation.

Refer to caption
Figure 3: Top: Time-domain 3He-129Xe free-precession comagnetometer signals detected with an in situ 87Rb magnetometer at an operating temperature of 120C. Bottom: Spectral density of comagnetometer in mm-scale cells.
Refer to caption
Figure 4: Our system stability is shown with Allan deviations for 3He-129Xe comagnetometer detected with an in situ 87Rb magnetometer. We compare a continuous-mode pulsed magnetometer operation with a Ramsey-style sequence that has an Rb-noble gas decoupling and dark period separating two detection periods. Both data sets are taken in a batch-fabricated anodically bonded cell at 120C, and have a practical limit of achievable steady-state 3He polarization.

To study the feasibility of these long noble gas lifetime cells in a compact gyroscope we operate the noble gases in a free-precession comagnetometer, similar to measurements in Ref. [31] that were with larger glass-blown cells. For these measurements we hyperpolarize 3He and 129Xe along a bias field B0=0.5subscript𝐵00.5B_{0}=0.5italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 µT that is co-linear to the pump and probe beams. We then apply a tipping pulse that places 3He and 129Xe spins in a plane transverse to B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The 87Rb pulse-train magnetometer detects the precession of both noble gases about B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT as shown in Fig. 3. To find the rotation rate for a given shot with both measurement styles, we first find the effective free-precession frequencies of the noble gases ωHe=γHeB0+Ωsubscript𝜔Hesubscript𝛾Hesubscript𝐵0Ω\omega_{\text{He}}=\gamma_{\text{He}}B_{0}+\Omegaitalic_ω start_POSTSUBSCRIPT He end_POSTSUBSCRIPT = italic_γ start_POSTSUBSCRIPT He end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + roman_Ω and ωXe=γXeB0+Ωsubscript𝜔Xesubscript𝛾Xesubscript𝐵0Ω\omega_{\text{Xe}}=\gamma_{\text{Xe}}B_{0}+\Omegaitalic_ω start_POSTSUBSCRIPT Xe end_POSTSUBSCRIPT = italic_γ start_POSTSUBSCRIPT Xe end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + roman_Ω. We extract the rotation rate ΩΩ\Omegaroman_Ω by taking the ratio fr=ωHe/ωXesubscript𝑓𝑟subscript𝜔Hesubscript𝜔Xef_{r}=\omega_{\text{He}}/\omega_{\text{Xe}}italic_f start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = italic_ω start_POSTSUBSCRIPT He end_POSTSUBSCRIPT / italic_ω start_POSTSUBSCRIPT Xe end_POSTSUBSCRIPT and finding Ω=ωXe(γrfr)/(γr1)Ωsubscript𝜔Xesubscript𝛾𝑟subscript𝑓𝑟subscript𝛾𝑟1\Omega=\omega_{\text{Xe}}(\gamma_{r}-f_{r})/(\gamma_{r}-1)roman_Ω = italic_ω start_POSTSUBSCRIPT Xe end_POSTSUBSCRIPT ( italic_γ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT - italic_f start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ) / ( italic_γ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT - 1 ) with γr=γHe/γXe=2.7540813(3)subscript𝛾𝑟subscript𝛾Hesubscript𝛾Xe2.75408133\gamma_{r}=\gamma_{\text{He}}/\gamma_{\text{Xe}}=2.7540813(3)italic_γ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = italic_γ start_POSTSUBSCRIPT He end_POSTSUBSCRIPT / italic_γ start_POSTSUBSCRIPT Xe end_POSTSUBSCRIPT = 2.7540813 ( 3 ) [57]. After a measurement the noble gas spins are placed back along the bias field by dumping feedback from the lock-in signal [58]. Repeating this sequence, we run the comagnetometer for an extended period of time to determine its sensitivity and long-term stability. In Fig. 4 we compare the Allan deviations for a Ramsey-style sequence with active Rb depolarization and for measurements with the Rb π𝜋\piitalic_π-pulse train magnetometer continuously running. For the Ramsey-style measurement each shot takes about 115 s, with a 25 s inital and 45 s final detection period. We fit the detection period signals to decaying sine waves to determine the phases ϕ1subscriptitalic-ϕ1\phi_{1}italic_ϕ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and ϕ2subscriptitalic-ϕ2\phi_{2}italic_ϕ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT with which the spins enter and leave the dark detection period Tdsubscript𝑇𝑑T_{d}italic_T start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT = 45 sec. Assuming there are no significant changes in the precession frequency, we find the in-the-dark free-precession frequencies of the noble gases by taking ω=Δϕ/Δt=(2nπ+ϕ2ϕ1)/Td𝜔Δitalic-ϕΔ𝑡2𝑛𝜋subscriptitalic-ϕ2subscriptitalic-ϕ1subscript𝑇𝑑\omega=\Delta\phi/\Delta t=(2n\pi+\phi_{2}-\phi_{1})/T_{d}italic_ω = roman_Δ italic_ϕ / roman_Δ italic_t = ( 2 italic_n italic_π + italic_ϕ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_ϕ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) / italic_T start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT, where n𝑛nitalic_n is an integer determined from an estimate of the precession frequency. During the dark period a two-axis decoupling pulse train is used to mitigate noble-gas frequency shifts due to 87Rb back polarization and nulls Bloch-Siegert shifts introduced by the pulse train [31]. For continuous measurements each shot lasts about 50 s, and we fit two decaying sine waves to find ωXesubscript𝜔Xe\omega_{\text{Xe}}italic_ω start_POSTSUBSCRIPT Xe end_POSTSUBSCRIPT and ωHesubscript𝜔He\omega_{\text{He}}italic_ω start_POSTSUBSCRIPT He end_POSTSUBSCRIPT. In both measurement styles we find frequency uncertainty is limited by the 3He frequency error due to low steady-state 3He magnetization at a temperature that allows for long 129Xe T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, as higher temperatures cause greater Rb-Xe relaxation. The Ramsey-style measurement integrates for more than a day, longer than the continous method by about an order of magnitude, with both methods reaching around 0.2 deg/h precision before drifting. The uncertainties achieved here are about an order of magnitude worse than achieved in the same system for glass-blown cells with 100 times larger volume used in Ref. [31]. Using microfabricated cells with higher 3He pressure would improve the performance, as in the current cells the small steady-state magnetization becomes insufficient with practical shot-to-shot pumping periods. To determine expected gyroscope performance we can use Cramér-Rao lower bounds (CRLB) of a decaying sine wave [59, 60, 61] to estimate a lower bound for the rotation error

δΩ=1γr1γr2δωXe2+δωHe2.𝛿Ω1subscript𝛾𝑟1superscriptsubscript𝛾𝑟2𝛿superscriptsubscript𝜔Xe2𝛿superscriptsubscript𝜔He2\delta\Omega=\frac{1}{\gamma_{r}-1}\sqrt{\gamma_{r}^{2}\delta\omega_{\text{Xe}% }^{2}+\delta\omega_{\text{He}}^{2}}.italic_δ roman_Ω = divide start_ARG 1 end_ARG start_ARG italic_γ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT - 1 end_ARG square-root start_ARG italic_γ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_δ italic_ω start_POSTSUBSCRIPT Xe end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_δ italic_ω start_POSTSUBSCRIPT He end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (1)

For a decaying sine wave,

δω2π=12C2πTm3/2SNR[Hz],δϕ=2DTm1/2SNR[rad],formulae-sequence𝛿𝜔2𝜋12𝐶2𝜋superscriptsubscript𝑇𝑚32SNRdelimited-[]Hz𝛿italic-ϕ2𝐷superscriptsubscript𝑇𝑚12SNRdelimited-[]rad\begin{split}\frac{\delta\omega}{2\pi}&=\frac{\sqrt{12C}}{2\pi T_{m}^{3/2}% \text{SNR}}~{}[\text{Hz}],\\ \delta\phi&=\frac{2\sqrt{D}}{T_{m}^{1/2}\text{SNR}}~{}[\text{rad}],\end{split}start_ROW start_CELL divide start_ARG italic_δ italic_ω end_ARG start_ARG 2 italic_π end_ARG end_CELL start_CELL = divide start_ARG square-root start_ARG 12 italic_C end_ARG end_ARG start_ARG 2 italic_π italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT SNR end_ARG [ Hz ] , end_CELL end_ROW start_ROW start_CELL italic_δ italic_ϕ end_CELL start_CELL = divide start_ARG 2 square-root start_ARG italic_D end_ARG end_ARG start_ARG italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT SNR end_ARG [ rad ] , end_CELL end_ROW (2)

where Tmsubscript𝑇𝑚T_{m}italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT is the measurement time and SNR =A/ρabsent𝐴𝜌=A/\rho= italic_A / italic_ρ, where A𝐴Aitalic_A is the initial sine wave amplitude and ρ𝜌\rhoitalic_ρ is the spectral density of white Gaussian noise. C𝐶Citalic_C and D𝐷Ditalic_D are factors depending on the measurement and relaxation time Tmsubscript𝑇𝑚T_{m}italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT and T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT that both approach unity as T2subscript𝑇2T_{2}\rightarrow\inftyitalic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT → ∞. For Ramsey-style detection in the limit of long T2subscript𝑇2T_{2}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, the optimal detection intervals are 1/6161/61 / 6 of the total measurement time and the overall statistical sensitivity is three times worse than for continuous detection of equal time. For finite relaxation one can show that

C=Tm3(e2TmT21)3T23cosh(2TmT2)3(T23+2T2Tm2),D1=e2TmT2Tm(2T2Tm+2Tm2+T22)T22Tm2T23cosh(2TmT2)2(2T2Tm2+T23),D2=T22Tm(e2TmT21)2Tm32T2Tm22T23cosh(2TmT2)2(2T2Tm2+T23),formulae-sequence𝐶superscriptsubscript𝑇𝑚3superscript𝑒2subscript𝑇𝑚subscript𝑇213superscriptsubscript𝑇232subscript𝑇𝑚subscript𝑇23superscriptsubscript𝑇232subscript𝑇2superscriptsubscript𝑇𝑚2formulae-sequencesubscript𝐷1superscript𝑒2subscript𝑇𝑚subscript𝑇2subscript𝑇𝑚2subscript𝑇2subscript𝑇𝑚2superscriptsubscript𝑇𝑚2superscriptsubscript𝑇22superscriptsubscript𝑇22subscript𝑇𝑚2superscriptsubscript𝑇232subscript𝑇𝑚subscript𝑇222subscript𝑇2superscriptsubscript𝑇𝑚2superscriptsubscript𝑇23subscript𝐷2superscriptsubscript𝑇22subscript𝑇𝑚superscript𝑒2subscript𝑇𝑚subscript𝑇212superscriptsubscript𝑇𝑚32subscript𝑇2superscriptsubscript𝑇𝑚22superscriptsubscript𝑇232subscript𝑇𝑚subscript𝑇222subscript𝑇2superscriptsubscript𝑇𝑚2superscriptsubscript𝑇23\begin{split}C&=\frac{{T_{m}}^{3}\left(e^{\frac{2{T_{m}}}{{T_{2}}}}-1\right)}{% 3{T_{2}}^{3}\cosh\left(\frac{2{T_{m}}}{{T_{2}}}\right)-3\left({T_{2}}^{3}+2{T_% {2}}{T_{m}}^{2}\right)},\\ D_{1}&=\frac{e^{\frac{2T_{m}}{T_{2}}}T_{m}\left(-2T_{2}T_{m}+2T_{m}^{2}+T_{2}^% {2}\right)-T_{2}^{2}T_{m}}{2T_{2}^{3}\cosh\left(\frac{2T_{m}}{T_{2}}\right)-2% \left(2T_{2}T_{m}^{2}+T_{2}^{3}\right)},\\ D_{2}&=\frac{T_{2}^{2}T_{m}\left(e^{\frac{2T_{m}}{T_{2}}}-1\right)-2T_{m}^{3}-% 2T_{2}T_{m}^{2}}{2T_{2}^{3}\cosh\left(\frac{2T_{m}}{T_{2}}\right)-2\left(2T_{2% }T_{m}^{2}+T_{2}^{3}\right)},\end{split}start_ROW start_CELL italic_C end_CELL start_CELL = divide start_ARG italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( italic_e start_POSTSUPERSCRIPT divide start_ARG 2 italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG end_POSTSUPERSCRIPT - 1 ) end_ARG start_ARG 3 italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT roman_cosh ( divide start_ARG 2 italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ) - 3 ( italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 2 italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG , end_CELL end_ROW start_ROW start_CELL italic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL = divide start_ARG italic_e start_POSTSUPERSCRIPT divide start_ARG 2 italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ( - 2 italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT + 2 italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) - italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT roman_cosh ( divide start_ARG 2 italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ) - 2 ( 2 italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ) end_ARG , end_CELL end_ROW start_ROW start_CELL italic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL = divide start_ARG italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ( italic_e start_POSTSUPERSCRIPT divide start_ARG 2 italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG end_POSTSUPERSCRIPT - 1 ) - 2 italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - 2 italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT roman_cosh ( divide start_ARG 2 italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ) - 2 ( 2 italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ) end_ARG , end_CELL end_ROW (3)

where D1subscript𝐷1D_{1}italic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is used to for estimate of the trailing phase of the first detection period and D2subscript𝐷2D_{2}italic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT for the leading phase of the second period. Using the signal height, noise, and relaxation times from Fig. 3, and considering detection, in-the-dark, and shot-to-shot times used in Fig. 4 we find a Cramér-Rao frequency uncertainty lower bound that corresponds to an angular random walk of 0.094 deg/hh\sqrt{\text{h}}square-root start_ARG h end_ARG for a Ramsey-style measurement and 0.056 deg/hh\sqrt{\text{h}}square-root start_ARG h end_ARG for continuous operation, about an order of magnitude better than we realized experimentally.

In conclusion, we demonstrate fabrication of miniature stemless anodically bonded vapor cells containing 3He, 129Xe, 87Rb, and N2 with long 129Xe and 3He coherence times. We use these cells in a dual noble-gas comagnetometer detected by an in situ 87Rb pulse-train magnetometer, and study the long-term stability of the system. Though far from fundamental limits, this system has great potential for high-precision, accurate rotation measurements in a small, low-cost package.

This work was funded by DARPA (Defense Advanced Research Project Agency) under contract FA8650-13-1-7326 and NSF Grant No. 1404325.

References

  • Liew et al. [2004] L.-A. Liew, S. Knappe, J. Moreland, H. Robinson, L. Hollberg, and J. Kitching, Microfabricated alkali atom vapor cells, Appl. Phys. Lett. 84, 2694 (2004).
  • Knappe et al. [2006] S. Knappe et al., J. Opt. A: Pure Appl. Opt. 8, S318 (2006).
  • Kitching [2018] J. Kitching, Chip-scale atomic devices, Applied Physics Reviews 5, 031302 (2018)https://doi.org/10.1063/1.5026238 .
  • Limes et al. [2020] M. Limes, E. Foley, T. Kornack, S. Caliga, S. McBride, A. Braun, W. Lee, V. Lucivero, and M. Romalis, Portable magnetometry for detection of biomagnetism in ambient environments, Phys. Rev. Applied 14, 011002 (2020).
  • McGilligan et al. [2020] J. P. McGilligan, K. R. Moore, A. Dellis, G. D. Martinez, E. de Clercq, P. F. Griffin, A. S. Arnold, E. Riis, R. Boudot, and J. Kitching, Laser cooling in a chip-scale platform, Applied Physics Letters 117, 054001 (2020)https://doi.org/10.1063/5.0014658 .
  • Ledbetter et al. [2008] M. Ledbetter, I. Savukov, D. Budker, V. Shah, S. Knappe, J. Kitching, D. Michalak, S. Xu, and A. Pines, Zero-field remote detection of nmr with a microfabricated atomic magnetometer, Proceedings of the National Academy of Sciences 105, 2286 (2008).
  • Jiménez-Martínez et al. [2014] R. Jiménez-Martínez, D. J. Kennedy, M. Rosenbluh, E. A. Donley, S. Knappe, S. J. Seltzer, H. L. Ring, V. S. Bajaj, and J. Kitching, Optical hyperpolarization and nmr detection of 129 xe on a microfluidic chip, Nature communications 5, 1 (2014).
  • Bulatowicz et al. [2013] M. Bulatowicz, R. Griffith, M. Larsen, J. Mirijanian, C. B. Fu, E. Smith, W. M. Snow, H. Yan, and T. G. Walker, Laboratory search for a long-range t𝑡titalic_t-odd, p𝑝pitalic_p-odd interaction from axionlike particles using dual-species nuclear magnetic resonance with polarized Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_Xe and Xe131superscriptXe131{}^{131}\mathrm{Xe}start_FLOATSUPERSCRIPT 131 end_FLOATSUPERSCRIPT roman_Xe gas, Phys. Rev. Lett. 111, 102001 (2013).
  • Donley [2010] E. A. Donley, Nuclear magnetic resonance gyroscopes, in 2010 IEEE Sensors Proc. (2010) p. 17.
  • Larsen and Bulatowicz [2012] M. Larsen and M. Bulatowicz, Nuclear magnetic resonance gyroscope, in 2012 IEEE Int. Freq. Cont. Symp. Proc. (2012).
  • Walker and Larsen [2016] T. Walker and M. Larsen, Chapter eight - spin-exchange-pumped {NMR} gyros, Adv. Atom, Mol. Opt. Phys. 65, 373 (2016).
  • Buchs et al. [2018] G. Buchs, S. Karlen, T. Overstolz, N. Torcheboeuf, E. Onillon, J. Haesler, and D. L. Boiko, AIP Conf. Proc. 1936, 020011 (2018).
  • Thrasher et al. [2019] D. A. Thrasher, S. S. Sorensen, J. Weber, M. Bulatowicz, A. Korver, M. Larsen, and T. G. Walker, Continuous comagnetometry using transversely polarized xe isotopes, Phys. Rev. A 100, 061403 (2019).
  • Bulatowicz and Larsen [2012] M. Bulatowicz and M. Larsen, Compact atomic magnetometer for global navigation, in 2012 IEEE/ION PLANS Proc. (2012) p. 1088.
  • Jiménez-Martínez et al. [2014] R. Jiménez-Martínez et al., Nat. Comm. 5, 3908 (2014).
  • Kennedy et al. [2017] D. Kennedy et al., Sci. Rep. 7, 43994 (2017).
  • Gentile et al. [2017] T. R. Gentile, P. J. Nacher, B. Saam, and T. G. Walker, Optically polarized He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_HeRev. Mod. Phys. 89, 045004 (2017).
  • Kastler [1973] A. Kastler, The hanle effect and its use for the measurements of very small magnetic fields, Nuclear Instruments and Methods 110, 259 (1973).
  • Lam et al. [1983] L. K. Lam, E. Phillips, E. Kanegsberg, and G. W. Kamin, Application Of CW Single-Mode GaAlAs Lasers To Rb-Xe NMR Gyroscopes, in Fiber Optic and Laser Sensors I, Vol. 0412, edited by E. L. Moore and O. G. Ramer, International Society for Optics and Photonics (SPIE, 1983) pp. 272 – 276.
  • Mehring et al. [1992] M. Mehring, A. Appelt, B. Menke, and P. Scheufler, High Precision Navigation 91: Proceedings of the 2nd International Workshop on High Precision Navigation, Stuttgart and Freudenstadt, 1991, edited by Linkwitz and Hangleiter (Dummler Verlag, 1992).
  • Chupp et al. [1987] T. E. Chupp, M. E. Wagshul, K. P. Coulter, A. B. McDonald, and W. Happer, Polarized, high-density, gaseous He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He targets, Phys. Rev. C 36, 2244 (1987).
  • Stoner et al. [1996] R. E. Stoner, M. A. Rosenberry, J. T. Wright, T. E. Chupp, E. R. Oteiza, and R. L. Walsworth, Demonstration of a two species noble gas maser, Phys. Rev. Lett. 77, 3971 (1996).
  • Bear et al. [2000] D. Bear, R. E. Stoner, R. L. 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, 5038 (2000).
  • Chupp and Coulter [1985] T. E. Chupp and K. P. Coulter, Polarization of Ne21superscriptNe21{}^{21}\mathrm{Ne}start_FLOATSUPERSCRIPT 21 end_FLOATSUPERSCRIPT roman_Ne by spin exchange with optically pumped rb vapor, Phys. Rev. Lett. 55, 1074 (1985).
  • Wu et al. [1987] Z. Wu, W. Happer, and J. M. Daniels, Coherent nuclear-spin interactions of adsorbed Xe131superscriptXe131{}^{131}\mathrm{Xe}start_FLOATSUPERSCRIPT 131 end_FLOATSUPERSCRIPT roman_Xe gas with surfaces, Phys. Rev. Lett. 59, 1480 (1987).
  • Chupp et al. [1988] T. E. Chupp, E. R. Oteiza, J. M. Richardson, and T. R. White, Precision frequency measurements with polarized He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He, Ne21superscriptNe21{}^{21}\mathrm{Ne}start_FLOATSUPERSCRIPT 21 end_FLOATSUPERSCRIPT roman_Ne, and Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_Xe atoms, Phys. Rev. A 38, 3998 (1988).
  • Butscher et al. [1994] R. Butscher, G. Wäckerle, and M. Mehring, Nuclear quadrupole interaction of highly polarized gas phase 131xe with a glass surface, The Journal of Chemical Physics 100, 6923 (1994)https://doi.org/10.1063/1.467008 .
  • Donley et al. [2009] E. A. Donley, J. L. Long, T. C. Liebisch, E. R. Hodby, T. A. Fisher, and J. Kitching, Phys. Rev. A 79, 013420 (2009).
  • Walker and Happer [1997] T. G. Walker and W. Happer, Spin-exchange optical pumping of noble-gas nuclei, Rev. Mod. Phys. 69, 629 (1997).
  • Appelt et al. [1998] S. Appelt, A. B.-A. Baranga, C. J. Erickson, M. V. Romalis, A. R. Young, and W. Happer, Theory of spin-exchange optical pumping of He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He and Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_XePhys. Rev. A 58, 1412 (1998).
  • Limes et al. [2018] M. E. Limes, D. Sheng, and M. V. Romalis, Phys. Rev. Lett. 120, 033401 (2018).
  • Korver et al. [2015] A. Korver, D. Thrasher, M. Bulatowicz, and T. G. Walker, Synchronous spin-exchange optical pumping, Phys. Rev. Lett. 115, 253001 (2015).
  • Newbury et al. [1993] N. R. Newbury, A. S. Barton, G. D. Cates, W. Happer, and H. Middleton, Gaseous 33{}^{3}\mathrm{-}^{3}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT - start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPThe magnetic dipolar spin relaxation, Phys. Rev. A 48, 4411 (1993).
  • Chann et al. [2002] B. Chann, I. A. Nelson, L. W. Anderson, B. Driehuys, and T. G. Walker, x129eXesuperscript𝑥129𝑒Xe{}^{129}xe-\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT italic_x italic_e - roman_Xe molecular spin relaxation, Phys. Rev. Lett. 88, 113201 (2002).
  • Anger et al. [2008] B. C. Anger, G. Schrank, A. Schoeck, K. A. Butler, M. S. Solum, R. J. Pugmire, and B. Saam, Gas-phase spin relaxation of Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_XePhys. Rev. A 78, 043406 (2008).
  • Dellis et al. [2016] A. Dellis, V. Shah, E. Donley, S. Knappe, and J. Kitching, Optics Letters 41, 2775 (2016).
  • Eklund et al. [2008] E. J. Eklund, A. M. Shkel, S. Knappe, E. Donley, and J. Kitching, Glass-blown spherical microcells for chip-scale atomic devices, Sensors and Act. A: 143, 175 (2008).
  • Lee et al. [2018] J. Lee, A. Almasi, and M. Romalis, Improved limits on spin-mass interactions, Phys. Rev. Lett. 120, 161801 (2018).
  • Balabas et al. [2010] M. V. Balabas, T. Karaulanov, M. P. Ledbetter, and D. Budker, Polarized alkali-metal vapor with minute-long transverse spin-relaxation time, Phys. Rev. Lett. 105, 070801 (2010).
  • Maul et al. [2016] A. Maul, P. Blümler, W. Heil, A. Nikiel, E. Otten, A. Petrich, and T. Schmidt, Rev. Sci. Instr. 87, 015103 (2016).
  • Stoltz et al. [1995] E. Stoltz, J. Tannenhauser, and P. Nacher, J. Low Temp. Phys. 101, 839 (1995).
  • Allmendinger et al. [2014] F. Allmendinger, W. Heil, S. Karpuk, W. Kilian, A. Scharth, U. Schmidt, A. Schnabel, Y. Sobolev, and K. Tullney, New limit on lorentz-invariance- and cpt𝑐𝑝𝑡cptitalic_c italic_p italic_t-violating neutron spin interactions using a free-spin-precession He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He-Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_Xe comagnetometer, Phys. Rev. Lett. 112, 110801 (2014).
  • Romalis et al. [2014] M. V. Romalis, D. Sheng, B. Saam, and T. G. Walker, Comment on “new limit on lorentz-invariance- and cpt𝑐𝑝𝑡cptitalic_c italic_p italic_t-violating neutron spin interactions using a free-spin-precession He3129Xesuperscript129superscriptHe3Xe{}^{3}\mathrm{He}\text{$-$}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He - start_POSTSUPERSCRIPT 129 end_POSTSUPERSCRIPT roman_Xe comagnetometer”, Phys. Rev. Lett. 113, 188901 (2014).
  • Limes et al. [2019] M. E. Limes, N. Dural, M. V. Romalis, E. L. Foley, T. W. Kornack, A. Nelson, L. R. Grisham, and J. Vaara, Dipolar and scalar He3129Xesuperscript129superscriptHe3Xe{}^{3}\mathrm{He}\text{$-$}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He - start_POSTSUPERSCRIPT 129 end_POSTSUPERSCRIPT roman_Xe frequency shifts in stemless cells, Phys. Rev. A 100, 010501 (2019).
  • Vaara and Romalis [2019] J. Vaara and M. V. Romalis, Calculation of scalar nuclear spin-spin coupling in a noble-gas mixture, Phys. Rev. A 99, 060501 (2019).
  • Terrano et al. [2019] W. A. Terrano, J. Meinel, N. Sachdeva, T. E. Chupp, S. Degenkolb, P. Fierlinger, F. Kuchler, and J. T. Singh, Frequency shifts in noble-gas comagnetometers, Phys. Rev. A 100, 012502 (2019).
  • Sachdeva et al. [2019] N. Sachdeva, I. Fan, E. Babcock, M. Burghoff, T. E. Chupp, S. Degenkolb, P. Fierlinger, S. Haude, E. Kraegeloh, W. Kilian, S. Knappe-Grüneberg, F. Kuchler, T. Liu, M. Marino, J. Meinel, K. Rolfs, Z. Salhi, A. Schnabel, J. T. Singh, S. Stuiber, W. A. Terrano, L. Trahms, and J. Voigt, New limit on the permanent electric dipole moment of Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_Xe using He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He comagnetometry and squid detection, Phys. Rev. Lett. 123, 143003 (2019).
  • Sheng et al. [2014] D. Sheng, A. Kabcenell, and M. V. Romalis, New classes of systematic effects in gas spin comagnetometers, Phys. Rev. Lett. 113, 163002 (2014).
  • Zhivun et al. [2019] E. Zhivun, M. Bulatowicz, A. Hryciuk, and T. Walker, Dual-axis π𝜋\piitalic_π-pulse magnetometer with suppressed spin-exchange relaxation, Phys. Rev. Applied 11, 034040 (2019).
  • Schaefer et al. [1989] S. R. Schaefer, G. D. Cates, T.-R. Chien, D. Gonatas, W. Happer, and T. G. Walker, Frequency shifts of the magnetic-resonance spectrum of mixtures of nuclear spin-polarized noble gases and vapors of spin-polarized alkali-metal atoms, Phys. Rev. A 39, 5613 (1989).
  • Romalis and Cates [1998] M. V. Romalis and G. D. Cates, Accurate He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He polarimetry using the rb zeeman frequency shift due to the Rb3Hesuperscript3RbHe\mathrm{Rb}{-}^{3}\mathrm{He}roman_Rb - start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT roman_He spin-exchange collisions, Phys. Rev. A 58, 3004 (1998).
  • Ma et al. [2011] Z. L. Ma, E. G. Sorte, and B. Saam, Collisional He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He and Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_Xe frequency shifts in rb˘noble-gas mixtures, Phys. Rev. Lett. 106, 193005 (2011).
  • Nahlawi et al. [2019] A. I. Nahlawi, Z. L. Ma, M. S. Conradi, and B. Saam, High-precision determination of the frequency-shift enhancement factor in rb-Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_XePhys. Rev. A 100, 053415 (2019).
  • Romalis et al. [1997] M. V. Romalis, E. Miron, and G. D. Cates, Pressure broadening of rb D1subscript𝐷1{D}_{1}italic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and D2subscript𝐷2{D}_{2}italic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT lines by 3he, 4he, n2, and xe: Line cores and near wings, Phys. Rev. A 56, 4569 (1997).
  • Nelson [2001] I. A. Nelson, Physics of practical spin-exchange optical pumping, Ph.D. thesis, University of Wisconsin–Madison (2001).
  • Schrank et al. [2009] G. Schrank, Z. Ma, A. Schoeck, and B. Saam, Characterization of a low-pressure high-capacity Xe129superscriptXe129{}^{129}\text{X}\text{e}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_X roman_e flow-through polarizer, Phys. Rev. A 80, 063424 (2009).
  • Makulski [2015] W. Makulski, 129xe and 131xe nuclear magnetic dipole moments from gas phase nmr spectra, Magn. Reson. Chem. 53, 273 (2015).
  • Alem et al. [2013] O. Alem, K. L. Sauer, and M. V. Romalis, Spin damping in an rf atomic magnetometer, Phys. Rev. A 87, 013413 (2013).
  • Chupp et al. [1994] T. E. Chupp, R. J. Hoare, R. L. Walsworth, and B. Wu, Spin-exchange-pumped He3superscriptHe3{}^{3}\mathrm{He}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_He and Xe129superscriptXe129{}^{129}\mathrm{Xe}start_FLOATSUPERSCRIPT 129 end_FLOATSUPERSCRIPT roman_Xe zeeman masers, Phys. Rev. Lett. 72, 2363 (1994).
  • Ying-Xian Yao and Pandit [1995] Ying-Xian Yao and S. M. Pandit, Cramer-rao lower bounds for a damped sinusoidal process, IEEE Transactions on Signal Processing 43, 878 (1995).
  • Gemmel et al. [2010] C. Gemmel, W. Heil, S. Karpuk, K. Lenz, C. Ludwig, Y. Sobolev, K. Tullney, M. Burghoff, W. Kilian, S. Knappe-Grüneberg, W. Müller, A. Schnabel, F. Seifert, L. Trahms, and S. Baeßler, Ultra-sensitive magnetometry based on free precession of nuclear spins, The European Physical Journal D 57, 303 (2010).