Abstract
Typical weak signal search experiments rely on resonant effects, where the resonance frequency is scanned over a broad range, resulting in significant time consumption. In this study, we demonstrate an accelerated strategy that surpasses the typical resonance-bandwidth limited scan step without compromising sensitivity. We apply this method to an alkali-noble-gas spin system, achieving an approximately 30-fold increase in scanning step size. Additionally, we obtain an ultrahigh sensitivity of 1.29 fTââ âHzâ1/2 at around 5 Hz, corresponding to an energy resolution of approximately 1.8âÃâ10â23eVââ âHzâ1/2, which is among the highest quantum energy resolutions reported. Furthermore, we use this sensor to search for axion-like particles, setting stringent constraints on axion-like particles (ALPs) in the 4.5â15.5 Hz Compton-frequency range coupling to neutrons and protons, improving on previous limits by several-fold. This accelerated strategy has potential applications in other resonant search experiments.
Similar content being viewed by others
Introduction
Precision measurements utilizing resonant effects, such as those found in optical cavities1, mechanical resonators2, and superconducting circuit resonances3, are widely employed in detecting weak signals. Ultralight bosonic particles, including axions and axion-like particles (ALPs), are viable dark matter (DM) candidates that exist in galaxies as oscillating fields. These particles can couple with various standard model particles, producing weak signals that motivate both astrophysical and laboratory experiments4,5,6,7,8. In our Milky Way galaxy, these particles are typically nonrelativistic with a velocity dispersion of approximately 220âkmââ âsâ1, leading to a sharp energy spectrum with a spread of about 10â6 9. Consequently, experimental searches for ALPs often rely on resonant effects10 to identify a distinct bump-like signal in the noise spectrum. However, when applying these techniques to ultralight bosonic dark matter, the resonance frequency must often be scanned due to the unknown mass of the dark matter. This necessitates tuning the resonant frequency over a wide range, which is both time-consuming and labor-intensive7. Resonant experiments for dark matter detection include those utilizing cavities11,12,13,14,15,16,17,18,19,20,21,22 and nuclear magnetic resonance (NMR) techniques23,24,25,26,27,28. In cavity experiments, the resonance frequency corresponds to the cavity resonant frequency, while in NMR experiments, it corresponds to the Larmor frequency of nuclear spins.
Our current ChangE-NMR experiment employs an ultrasensitive atomic magnetometer, integral to precision measurements and the exploration of physics beyond the Standard Model23,24,25,26,27,28. The magnetometer is based on hybrid alkali-noble-gas spin ensembles and operates close to the spin-exchange relaxation-free (SERF) regime. The advantages of our alkali-noble-gas spin sensor over other ultrasensitive alkali-only magnetometer and K-3He magnetometer are the fine noise suppression and the utilization of 21Ne spins, whose gyromagnetic ratio is about three orders of magnitude lower than alkali spin and about one order of magnitude lower than 3He respectively. Typically, comagnetometers such as the self-compensation (SC) comagnetometer29,30 and clock-comparison (CC) comagnetometer31 are used to suppress magnetic noise. However, ultrasensitive SC comagnetometers are limited to a relatively low- frequency range due to the low resonance frequency of noble-gas nuclear spin; the sensitivity of CC comagnetometers are not high enough for our ALP searches. Thus, instead of working in the comagnetometer regime, we use another approach to work in the NMR mode by developing a ferrite magnetic shield with sub-fT magnetic noise and successfully reducing magnetic noise from heating 32.
By sweeping the external magnetic field to scan the resonant frequency, we record a time series of voltage readings from the probe detector. Employing Fourier transformation, we derive the power spectral density (PSD) of the voltage within the frequency domain near the resonant frequency, Fig. 1. The PSD typically exhibits a flat profile without the manual injection of fabricated signals27,28, as the background noise is predominantly white. Nevertheless, our experiment reveals a distinct peaked PSD. We categorize these two situations as flat and peaked regimes and provide analytical explanations for their unique characteristics. A similar feature was discussed in refs. 33,34 and used for DM search projections 35,36,37,38,39,40,41,42. In the peaked regime, the sensitivity is limited by amplifiable noise, and the high amplification factor does not enhance the sensitivity. However, we discovered that, with the high amplification factor, one can significantly extend the bandwidth while maintaining similar sensitivity levels, between one and two orders of magnitude. We have developed and experimentally verified an expression for this expanded bandwidth. This approach is applicable to various types of resonance experiments, effectively accelerating the scanning process (see refs. 43,44 for different approaches).
We operate in the NMR mode corresponding to the peaked regime. As a result, this regime results in an \({{{{{{{\mathcal{O}}}}}}}}(30)\) increase in the full width at half maximum (FWHM) of the sensitivity peak compared to the width of the spin resonance, thus accelerating scan speeds by the same factor, see the scan step boost factor in Fig. 2. We incorporate the stochastic nature of ultralight axion DM into the sensitivity calculation, extending it to the scenario with the magnetometer sensitive to the ALP-field gradient (see below) along two transverse axes. With a scan step of 0.25âHz, much larger than the spin-resonance width of about 0.006âHz, our investigations yield improved limits for axion-nucleon couplings within the 5â15âHz frequency range, see Fig. 3.
We achieve energy resolution of about 1.8âÃâ10â23âeVââ âHzâ1/2, which is about two orders of magnitude better than the most sensitive alkali-spin magnetometer45, and about three times better than the most sensitive K-3He comagnetometer31. Exploring sub-eV (<1017âHz) parameter space for ultralight bosonic dark matter remains meaningful and crucial due to the unknown DM mass, necessitating comprehensive investigation. With our current setup, we achieve optimal sensitivity around \({{{{{{{\mathcal{O}}}}}}}}(10)\)âHz, offering a valuable window into the parameter space. However, the vastness of this space presents exciting opportunities for future exploration and discovery.
Results
Properties of two regimes
Given the sharp spectral peak of the resonant signal, we confine our analysis to a frequency band centered around the peak. Within this band, the noises are flat and can be characterized as random Gaussian white noise. We classify the noise into two categories: non-amplifiable noise Nna (e.g., photon shot noise) and amplifiable noise Na (e.g., magnetic noise).
We analyze the signal assuming amplification at the resonant frequency f0 by the amplification factor q and Lorentzian amplification function,
where Πrepresents the width of the function. Due to the independence between noise sources, the total noise variance Σ(f) is
where Σna and Σa denote the variances for non-amplifiable and amplifiable noise, respectively.
The variance of the signal takes the form:
where gs signifies the signal strength. For simplicity, we assume the measurement spans a coherence time46,47,48\({\tau }_{a} \sim 2\hslash /({m}_{a}{v}_{{{{{{{{\rm{vir}}}}}}}}}^{2})\), with ma and vvir the mass and virial velocity of the ALP, ensuring that the signal is confined to a single frequency bin within the data. Consequently, confidence-level limits can be established at 95%, which yields the expression
with c0â=â5.9949, the experimental premises of qââ«â1 and Σnaââ«âΣa.
For a single measurement at a resonant frequency f0, the best sensitivity \({g}_{{{{{{{{\rm{s}}}}}}}}}^{\min }\) is achieved at fâ=âf0 and follows a Lorentzian profile. Additionally, we can establish the FWHM \({\Gamma }_{{{{{{{{\rm{FWHM}}}}}}}}}^{{{{{{{{\rm{sen}}}}}}}}}\) to characterize the sensitivity bandwidth, requiring \({g}_{{{{{{{{\rm{s}}}}}}}}}({f}_{0}+{\Gamma }_{{{{{{{{\rm{FWHM}}}}}}}}}^{{{{{{{{\rm{sen}}}}}}}}}/2)=2{g}_{{{{{{{{\rm{s}}}}}}}}}^{\min }\). Therefore, we have
with κââ¡âq2Σa/Σna, reflecting the level of amplification.
In Fig. 1, we maintain constant noise variances, Σa and Σna, while employing two amplification factors, q1,2 with q1ââªâq2, to illustrate the distinct regimes. Monte-Carlo simulations are utilized to generate PSD data for both regimes in the top panel, with corresponding sensitivities presented in the bottom panel. Note that experimentally, the value of q can be controlled via changing the degree of polarization of the noble gas and its coherence time as discussed below.
In the flat regime, κââªâ1, the peak sensitivity follows \({g}_{{{{{{{{\rm{s}}}}}}}}}^{\min }\,\approx \,\sqrt{{c}_{0}{\Sigma }_{{{{{{{{\rm{na}}}}}}}}}}/q\), which can still be improved by increasing the amplification factor or decreasing the non-amplifiable noise, and it is not affected by the amplifiable noise. The single scan of the experiment can effectively cover the region of \({\Gamma }_{{{{{{{{\rm{FWHM}}}}}}}}}^{{{{{{{{\rm{sen}}}}}}}}}\,\approx \,2\sqrt{3}\Gamma\) which is not improved by increasing κ. The characteristic of this regime is a flat PSD near the resonant frequency, which is common in the previous cavity experiments (e.g., ADMX50,51,52, HAYSTAC53, CAST-CAPP22 with flat excess PSD, even for superconducting radiofrequency (SRF) cavity54,55,56) and magnetometer experiments (e.g., ref. 27, NASDUCK28). In general, for the flat regime, the scan step is taken to be the largest physical width, the cavity resonant linewidth for cavity, the nuclear linewidth for NMR, and DM linewidth for SRF cavity.
Conversely, in the peaked regime (right panel of Fig. 1), κââ«â1, the peak sensitivity is approximately \({g}_{{{{{{{{\rm{s}}}}}}}}}^{\min }\,\approx \,\sqrt{{c}_{0}{\Sigma }_{{{{{{{{\rm{a}}}}}}}}}}\), entirely dictated by the amplifiable noise Σa. Reducing the non-amplifiable noise Σna and increasing amplification can no longer enhance sensitivity. However, adjusting these parameters can augment the spectral width \({\Gamma }_{{{{{{{{\rm{FWHM}}}}}}}}}^{{{{{{{{\rm{sen}}}}}}}}}\,\approx \,2\sqrt{3}\Gamma \sqrt{\kappa }\) by a factor of \(\sqrt{\kappa }\), see also refs. 33,34.
ALP search principle
We seek ALP DM that couples to nuclei through the gradient interaction given by the Hamiltonian
where ÏN represents the nucleus spin, encompassing both neutron and proton contributions, and gaNN signifies the effective axion-nucleus coupling. Consequently, the axion-wind operator induces an anomalous magnetic field baââ¡âgaNNâa/γN, with γN denoting the nuclear gyromagnetic ratio.
We have employed the alkali-noble-gas magnetometer with hybrid spins of 21Ne-Rb-K which are contained in a glass vapor cell, as shown in Fig. 4a. Through spin-exchange optical pumping, the pump light polarizes the alkali atoms, which subsequently polarize the noble-gas nuclear spins via spin-exchange interactions with alkali electronic spins. To enhance the polarization homogeneity and efficiency, hybrid optical pumping is employed, following refs. 57,58. When we apply a bias field along the pump-light direction, the noble-gas nuclear spins are resonant at the Larmor frequency. The precession of polarized 21Ne nuclear spins occurs under the influence of ALP fields. When the ALP-field frequency approaches the resonance frequency, the precession of 21Ne nuclear spins is maximized. While traditional methods employ pick-up coils or extra magnetometers to read the dipole magnetic field due to the precession of noble-gas nuclear spins, we utilize in situ alkali spins for this purpose. The Fermi-contact interaction between overlapping alkali and noble-gas spins significantly augments this dipole magnetic field with an amplification factor q. Therefore, the precession of noble-gas spins due to the ALP field is amplified and read out by in situ alkali spins. In this Fermi-contact enhanced NMR, the response to the ALP field is amplified by two factors: the resonance of noble-gas spins and the Fermi-contact enhancement. Finally, the precession of alkali spins is measured with linearly polarized probe light via optical rotation.
We employ magnetic fields for calibrating the sensor response to ALP fields. The response to the ALP field, generally different from that to a magnetic field59,60,61, is determined by relating the responses to magnetic and pseudo-magnetic fields. Due to the pump-light alignment along the z direction, the symmetry of the x and y components in the Bloch equation causes a 90-degree phase difference between their respective response phases. When nuclear spin precession frequency significantly exceeds its relaxation rate, the nuclear spin projection shifts between x and y axis rapidly, making the projections nearly equal, see Methods 62,63,64,65.
However, the lineshape of amplification for the usual magnetic field is not entirely symmetric due to the well-known Fano resonance66,67, as illustrated in Fig. 4b. When both resonant (21Ne) and continuous-background (alkalies) transitions coexist in a single system, their interplay results in Fano resonances. The amplification function, incorporating the Fano effect, is described by
with the resonant frequency f0 and the relaxation rate \({R}_{2}^{{{{{{{{\rm{n}}}}}}}}}={({T}_{2{{{{{{{\rm{n}}}}}}}}})}^{-1}=2\pi \Gamma\). The amplification factor q is primarily determined by the nuclear spin coherence time T2n and the effective magnetic field \({B}_{0}^{{{{{{{{\rm{n}}}}}}}}}\), represented as \(q=\frac{1}{2}{\gamma }_{{{{{{{{\rm{Ne}}}}}}}}}{B}_{0}^{{{{{{{{\rm{n}}}}}}}}}{T}_{2{{{{{{{\rm{n}}}}}}}}}\) in Methods. The magnitude of q determines the degree of asymmetry, with a higher q indicating a more pronounced resonance. As q approaches infinity, the Fano lineshape transforms into a symmetric Lorentz lineshape in Eq. (1) as ηF(f)âââ1â+âη(f).
In Fig. 4b, we set the resonance frequency to 4.98âHz and measure the q. We fit the measured q with Eq. (7) and find it to be in good agreement with the experimental parameters \({B}_{0}^{{{{{{{{\rm{n}}}}}}}}}\simeq 250\)ânT and T2nâââ31âs. In Fig. 4c, we show the noise spectrum with a Fano resonance seen. The Fano resonance has minimal impact on sensitivity at the resonance frequency, but significantly affects the non-resonance region. For a usual magnetic field, Fano lineshape accurately describes the response of the system. However, since we only consider the ALP field coupling to the nuclear spins, the response is more accurately described by a Lorentzian shape. A detailed analytical derivation and simulations for both usual magnetic fields and ALP fields are provided in Methods.
The measured sensitivity is shown in Fig. 4c. The energy resolution at 5.25âHz reaches 1.8âÃâ10â23âeVââ âHzâ1/2. The best sensitivity is achieved at the gain peak seen in Fig. 4b, where the amplification exceeds 100-fold. Here, the sensitivity is limited by the (amplified) magnetic noise. However, the Fano effect arising from the interference between resonant (21Ne) and continuous-background (alkali) transitions, renders a narrow region on the left side of the resonant peak insensitive to signals (see also ref. 67). Analyzing the averaged sensitivity curve reveals that beyond the resonant region, the averaged sensitivity hovers around ~1âÃâ10â22âeVââ âHzâ1/2, influenced by non-amplifiable noise, probe background noise (shown with dot-dashed line), which is close to the theoretical value of photon shot noise. The sensitivity achieved is below the spin projection noise and photon shot noise of the alkali-metal magnetometer utilizing Fermi-contact amplification68.
ALP search results
We conducted 49 individual scans within a step interval of approximately 0.25âHz, chosen according to Eq. (5), with a cumulative measurement duration of 1440âhours, over the frequency range of 5â15âHz. Employing frequency-domain analysis outlined in ref. 69, we searched for ALP signals, taking into consideration the stochastic effects47. It is important to note a distinction from ref. 69: while their investigation involves a sole sensitive axis, our study incorporates both the x and y axes as sensitive directions in the NMR mode. Consequently, we have adapted the algorithm accordingly.
When compared to the previous simple estimation, the Fano effect can significantly modify the response of the system to magnetic fields, while having less impact on the response to ALP. In Methods, we account for this effect and calculate the sensitivity bandwidth and the best-sensitivity frequency,
Because Σnaââ«âΣa we have q2ââ«âκ, thus Eq. (8) reverts to Eq. (5).
In the experiment, we have two apparatus, denoted as A and B in Methods, collecting the data. Apparatus A has \(\left\vert q\right\vert \sim 40-50\) and κ in the range of 10â1000, while apparatus B has \(\left\vert q\right\vert \sim 100\) and κ ranging between 400 and 1000. We obtain the values of q,âÎ from the calibration process and extract κ from the data. Therefore, we can compute ÎFWHM and \({\tilde{f}}_{0}\) and compare them with the experiment results. In Fig. 5, we observe that the theory and experimental results are in good agreement for ÎFWHM, differing by only 7%. Meanwhile, the actual PSD and sensitivity are shown, confirming the shift of \({\tilde{f}}_{0}\) and the enhancement of ÎFWHM.
After computing the sensitivity, we evaluate the ALP-field nucleon couplings as illustrated in Fig. 3. The coupling constants gann for neutrons and gapp for protons are derived following70. The sensitivity of every scan spans over a range of a few tenth of a hertz. Our sensitivity surpasses our prior ChangE-HSR study by a few-fold in both couplings. Additionally, our results exceed the astrophysical bounds within the frequency range of 4.6â6.6âHz and around 7.5âHz.
The goal of our investigation is the detection of possible ALP-dark-matter signals in our dataset. We undertake a post-analysis to evaluate the significance of the best-fit signal in comparison to the background-only model and to statistically assess data points surpassing the 5Ï threshold69,70. We found several candidates surpass the 5Ï threshold of the background-only hypothesis, with consideration for the look-elsewhere effect. However, each of these candidate frequencies can be tested by multiple adjacent individual scans. Notably, these 5Ï peaks were not consistently present in other individual scans, implying potential transient background fluctuations. Consequently, we conclude the absence of ALP DM signals in our experiments.
In summary, we conducted a search for ultralight ALP DM using a magnetometer operating in the Fermi-contact enhanced NMR mode. Our findings revealed that the experiment operates in the âpeakedâ regime, dominated by magnetic noise. In this regime, we optimize the scanning speed by choosing a step size exceeding the intrinsic NMR linewidth by approximately a factor of 30, as determined by the amplification parameter and the noise level in the system. We analytically elucidated the properties of this peaked regime and found that the sensitivity bandwidth inferred from the parameters of the system agrees with experimental results at a 7% level. Furthermore, we established stringent constraints on the coupling of ALP to both neutrons and protons, surpassing previous searches by a few-fold. Our laboratory results surpass the astrophysical ones derived from neutron-star cooling in the frequency range of 4.6â6.6âHz and around 7.5âHz.
Discussion
This search is rooted in the NMR regime, which diverges significantly from the previous self-compensation (SC) regime30,32 and the hybrid-spin-resonance (HSR) regime70, despite their common basis in the K-Rb-21Ne spin ensembles. In the SC regime, the magnetization of noble-gas nuclear spins compensates for low-frequency magnetic noise under a certain external bias field, where the value equals the sum of the Fermi-contact interaction fields of alkali spins and noble-gas nuclear spins. There is no physical spin resonance in the SC regime. In contrast, in the HSR regime, the alkali spins and noble-gas nuclear spins are strongly coupled by setting the bias field to a specific value, which decreases the total magnetic field experienced by the alkali spins, distinguishing it from the SC regime. Actually, in the HSR regime, the dynamic performance of noble-gas nuclear spins is significantly amplified by several orders of magnitude through coupling to the highly dynamic alkali spins, resulting in resonance between noble-gas nuclear spins and alkali spins, where their Zeeman splittings are similar to each other. However, in the NMR regime, the bias magnetic field is tuned such that the Zeeman splitting of the nuclear spin matches the dark matter mass, enabling resonant absorption of the dark matter particles. To explore a wide range of dark matter masses, the bias field has to be scanned. In this experiment, the resonance frequencies of noble-gas nuclear spins and alkali spins are widely separated by several kHz, unlike the HSR regime. Moreover, there is no compensation for magnetic noise compared to the SC regime. Furthermore, in the NMR regime, the resonance is asymmetric, with a profile significantly different from the symmetric Lorentzian shape, attributed to the well-known Fano resonance. Therefore, the NMR regime search differs from searches conducted under the SC and HSR regimes.
In resonant dark matter search experiments, achieving both high sensitivity and broad coverage is often challenging. While maintaining world-leading sensitivity can be difficult, achieving a wide search range is also not straightforward. In cases where the measurement sensitivity is limited by amplifiable noise, increasing the response amplitude (usually denoted as the amplification factor q) does not necessarily improve the sensitivity. However, we found that increasing the amplification factor can extend the measurement bandwidth. We formulated expressions for the expanded bandwidths, determined by the amplification factor and the ratio of amplifiable noise to non-amplifiable noise. The validity of the bandwidth-expansion method was verified using ultrasensitive hybrid spin-based sensors. This approach can be applied also to other types of resonance experiments to broaden the bandwidth. The analysis will guide the design of future, higher-sensitivity experiments.
The ADMX and CAPP experiments have achieved unprecedented sensitivity in probing the Quantum chromodynamics(QCD) axion region through their investigations of the axion-photon coupling71,72. Similarly, the JEDI Collaboration introduced a storage ring for pioneering axion dark matter research73. These experiments may target vastly disparate frequencies, resulting in varying coherence times for ALP-dark matter detection. For instance, our experiment targets frequencies around 10âHz, whereas cavity experiments often operate near 1âGHz. This discrepancy leads to eight orders of magnitude differences in the coherent time. Additionally, the scan step is influenced by experiment quality factors. For experiments with quality factors lower than the dark matterâs ~106, the scan step typically equals the measurement width of the system. Conversely, experiments with higher quality factors (e.g., superconducting cavities) often adopt scan steps equivalent to the dark matter signal width. This difference affects the number of data points measured and, in conjunction with total measurement time, dictates search coverage.
Methods
In the Methods section, we begin by providing a detailed description of the experimental setup and its fundamental underlying principles. This includes a discussion on the manifestation of the Fano effect within our experiment and its impact on our results. Subsequently, we delve into the stochastic nature of ultralight axion DM and its implications. Finally, we explore the statistical methods employed to calculate signal sensitivity and outline the post-analysis process for discriminating potential axion DM candidates.
Experiment setup
The experimental setup for the axion DM search is illustrated in Fig. 6. A spherical vapor cell with a diameter of 12âmm containing alkali-metal atoms and noble-gas atoms is at the center of the apparatus. The cell is heated to 468âK with several flexible heating films driven by AC current. The films are arranged to compensate for the magnetic field produced by the heating current. The cell and oven are inside a polyetheretherketone (PEEK) chamber, and a water-cooling tube is wrapped around the PEEK chamber to prevent the ferrite outside the chamber from heating up and generating thermo-magnetic noise. The geomagnetic environment around the cell is shielded by a five-layer magnetic shield and magnetic compensation coils wrapped on the PEEK chamber. Atomic spins are polarized with circularly polarized pump light tuned to the D1 line of K atoms propagating along the z-axis. Alkali-atom polarization is transferred to the noble-gas nuclei in spin-exchange collisions. The precession information is read out with a linearly polarized probe light propagating along x, which is detuned towards longer wavelengths by 0.5ânm from the D1 line of Rb atoms. In the presence of a signal, which is the rotation of the nuclear spins, the polarization direction of the linearly polarized probe light is rotated upon passing through the cell, and the optical rotation angle can be read by detecting the transmitted light intensity. The probe light polarization is modulated with a photoelastic modulator (PEM) at a frequency of 50âkHz to suppress low-frequency noise. The transmitted intensity detected with a photodiode is demodulated with a lock-in amplifier. In order to suppress low-frequency interference caused by air convection, the cell, the magnetic shield system, and the optical system are placed in a water-cooling box (not shown in the figure).
The responses to the magnetic field and ALP field
The spin dynamics of the Fermi-contact enhanced spin resonator could be described by the coupled Bloch equations29,61,74:
where γe and γNe are gyromagnetic ratios of electrons and 21Ne nucleons, respectively. Q is the slowing-down factor arising from spin-exchange collisions and hyperfine interaction, Pe and Pn are the spin polarizations of alkali electrons and 21Ne nucleons, respectively. B and Ω are external magnetic field and inertial rotation. βe and βNe are ALP-fields couplings to alkali electron spins and noble-gas nuclear spins, respectively. Here we focus on the βNe. The Fermi-contact interaction between alkali atoms and 21Ne atoms can be described by an effective magnetic field \(\lambda {M}_{0}^{{{{{{{{\rm{e}}}}}}}},{{{{{{{\rm{n}}}}}}}}}{{{{{{{{\bf{P}}}}}}}}}^{{{{{{{{\rm{e}}}}}}}},{{{{{{{\rm{n}}}}}}}}}\), where \({M}_{0}^{{{{{{{{\rm{n}}}}}}}}}\) (\({M}_{0}^{{{{{{{{\rm{e}}}}}}}}}\)) is the maximum magnetization of 21Ne nucleon (alkali electron). For a uniformly spin-polarized spherical cell, λâ=â8Ïκ0/3, where κ0 is the enhancement factor. \({P}_{{{{{{{{\rm{z}}}}}}}}0}^{{{{{{{{\rm{e}}}}}}}}}\) and \({P}_{{{{{{{{\rm{z}}}}}}}}0}^{{{{{{{{\rm{n}}}}}}}}}\) are steady electronic and nuclear spin polarization. \({R}_{1}^{{{{{{{{\rm{e}}}}}}}}}\) and \({R}_{2}^{{{{{{{{\rm{e}}}}}}}}}\) are the longitudinal and transverse relaxation rates for alkali electron spin, respectively, and \({R}_{1}^{{{{{{{{\rm{n}}}}}}}}}\) and \({R}_{2}^{{{{{{{{\rm{n}}}}}}}}}\) are the longitudinal and transverse relaxation rates for the 21Ne nucleon spin.
The response of the nuclear spin to the ALP field in the rotating coordinate system can be obtained as
Here, \({\tilde{P}}_{\perp }^{{{{{{{{\rm{n}}}}}}}}}={\tilde{P}}_{{{{{{{{\rm{x}}}}}}}}}^{{{{{{{{\rm{n}}}}}}}}}+{{{{{{{\rm{i}}}}}}}}{\tilde{P}}_{{{{{{{{\rm{y}}}}}}}}}^{{{{{{{{\rm{n}}}}}}}}}\), \({\tilde{\beta }}_{\perp }^{{{{{{{{\rm{Ne}}}}}}}}}={\tilde{\beta }}_{{{{{{{{\rm{x}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}+{{{{{{{\rm{i}}}}}}}}{\tilde{\beta }}_{{{{{{{{\rm{y}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}},{\omega }_{{{{{{{{\rm{n}}}}}}}}}={\gamma }_{{{{{{{{\rm{Ne}}}}}}}}}\left({B}_{{{{{{{{\rm{z}}}}}}}}}+\lambda {M}_{0}^{{{{{{{{\rm{e}}}}}}}}}{P}_{{{{{{{{\rm{z}}}}}}}}0}^{{{{{{{{\rm{e}}}}}}}}}\right)\) is the Larmor frequency of nuclear spin, and Ï is the frequency of the ALP field. \({R}_{2}^{{{{{{{{\rm{n}}}}}}}}}={({T}_{2{{{{{{{\rm{n}}}}}}}}})}^{-1}\). Then the effective nuclear magnetic field arising in response to the ALP field can be expressed as \({\tilde{B}}_{\perp }^{{{{{{{{\rm{eff}}}}}}}}}=\lambda {M}_{0}^{{{{{{{{\rm{n}}}}}}}}}{\tilde{P}}_{\perp }^{{{{{{{{\rm{n}}}}}}}}}\). Thus, the response of electronic spin induced by the nuclear response magnetic field can be obtained as
with \({\omega }_{{{{{{{{\rm{e}}}}}}}}}=\frac{{\gamma }_{{{{{{{{\rm{e}}}}}}}}}}{Q}\left({B}_{{{{{{{{\rm{z}}}}}}}}}+\lambda {M}_{0}^{{{{{{{{\rm{n}}}}}}}}}{P}_{{{{{{{{\rm{z}}}}}}}}0}^{{{{{{{{\rm{n}}}}}}}}}\right)\). It can be seen that the coefficients of \({\beta }_{{{{{{{{\rm{x}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}\) and \({\beta }_{{{{{{{{\rm{y}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}\) only differ by one imaginary unit, so the x and y axes are symmetric, and the responses of nuclear and electronic spins to the ALP fields in the x and y directions differ by 90 degrees of phase.
The ALP field is equivalent to an oscillating magnetic field which couples to nuclear spins, so the possibility of using the oscillating magnetic field to calibrate the ALP field is analyzed. Magnetic or non-magnetic perturbations can be expressed in a complex form. Based on the symmetry of the x and y axes, the transverse input can be expressed as Bâ¥â=âBxâ+âiBy and \({\beta }_{\perp }^{{{{{{{{\rm{Ne}}}}}}}}}={\beta }_{{{{{{{{\rm{x}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}+{{{{{{{\rm{i}}}}}}}}{\beta }_{{{{{{{{\rm{y}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}\), respectively. Thus, the response of the sensor corresponding to the two parts of the signal can be calculated separately as \({M}_{+}^{{{{{{{{\rm{e}}}}}}}}}{{{{{{{{\rm{e}}}}}}}}}^{{{{{{{{\rm{i}}}}}}}}\omega t}\) and \({M}_{-}^{{{{{{{{\rm{e}}}}}}}}}{{{{{{{{\rm{e}}}}}}}}}^{-{{{{{{{\rm{i}}}}}}}}\omega t}\), where
Here, \({\omega }_{{{{{{{{\rm{en}}}}}}}}}=\frac{{\gamma }_{{{{{{{{\rm{e}}}}}}}}}}{Q}\lambda {M}_{0}^{{{{{{{{\rm{n}}}}}}}}}{P}_{{{{{{{{\rm{z}}}}}}}}0}^{{{{{{{{\rm{e}}}}}}}}},{\omega }_{{{{{{{{\rm{ne}}}}}}}}}={\gamma }_{{{{{{{{\rm{Ne}}}}}}}}}\lambda {M}_{0}^{{{{{{{{\rm{e}}}}}}}}}{P}_{{{{{{{{\rm{z}}}}}}}}0}^{{{{{{{{\rm{n}}}}}}}}},{\omega }_{{{{{{{{\rm{nn}}}}}}}}}={\gamma }_{{{{{{{{\rm{Ne}}}}}}}}}\lambda {M}_{0}^{{{{{{{{\rm{n}}}}}}}}}{P}_{{{{{{{{\rm{z}}}}}}}}0}^{{{{{{{{\rm{n}}}}}}}}}\). When the ALP field or the magnetic field is resonant with the nuclear spins, the coefficients of responses to the ALP field and magnetic field are approximately equal (when \({R}_{2}^{{{{{{{{\rm{n}}}}}}}}}\,\ll \,{\omega }_{{{{{{{{\rm{nn}}}}}}}}}\)). Therefore, the ALP field can be calibrated with an oscillating magnetic field. Since the sensor detects the information of the x-axis electron spin polarization, the x-direction information can eventually be read by separating the real part as \({P}_{{{{{{{{\rm{x}}}}}}}}}^{{{{{{{{\rm{e}}}}}}}}}={{{{{{{\rm{Re}}}}}}}}\left({M}_{+}^{{{{{{{{\rm{e}}}}}}}}}{{{{{{{{\rm{e}}}}}}}}}^{{{{{{{{\rm{i}}}}}}}}\omega t}+{M}_{-}^{{{{{{{{\rm{e}}}}}}}}}{{{{{{{{\rm{e}}}}}}}}}^{-{{{{{{{\rm{i}}}}}}}}\omega t}\right)\). According to the solution, the response to the usual magnetic field and ALP field can be simulated as shown in Fig. 7. It can be seen that the nuclear spin-resonance frequency is 16.38âHz when a 5000ânT bias magnetic field is applied, and the response of the nuclear spin to the usual magnetic field as well as the ALP field is almost equal at the resonance. It is noteworthy that the response to the usual magnetic field is asymmetric because both nuclear and electronic spins respond to the usual magnetic field67. However, the response to the ALP field is more symmetric because there is only discrete resonance (due to 21Ne) in the system. Thus, the response to the usual magnetic field is fitted by the Fano lineshape and the response to the ALP field is calibrated by the Lorentz lineshape in the small range around the resonance peak.
FANO effect and its response
The magnetic field response can be used to calibrate the ALP field. In practice, we find that the exact lineshape of usual magnetic field amplification as a function of the test frequency is asymmetric with the profile significantly distorted compared to the symmetric Lorentzian shape, consistent with the simulation in Fig. 7. This asymmetric lineshape is caused by the well-known Fano resonance. Specifically, Rb and 21Ne spins both experience the measured oscillating usual magnetic field and both generate signatures at the frequency of the measured field, and the asymmetric shape is caused by the interference of the response of the two spin ensembles as shown in Fig. 8. The resonance frequencies of noble-gas and alkali-metal atoms are different, similar to pendulums with two different oscillation periods. The resonance between the ALP field and noble-gas atoms causes a significant oscillation of the ânoble-gas pendulumâ. Subsequently, due to the Fermi-contact coupling between the noble gas and the alkali-metal atoms (shown as a spring), the alkali-metal atoms are pulled to oscillate.
The Rb spins and noble gas 21Ne spins both react to the oscillating usual magnetic field and generate their own responses. Unlike the case in the self-compensation comagnetometer, where the resonance frequencies of noble-gas atoms and alkali atoms are close to each other in the near-zero field condition, here, the bias field is large enough that noble-gas atoms and alkali atoms are decoupled. In the 21Ne-Rb-K experiment, due to the long coherence time of the 21Ne nuclear spins (T2n ~40âs), the 21Ne spins have a sharp resonance line and, therefore, can be considered as having a discrete system. The Rb-K spins have a shorter coherence time (T2e ~1âms), and thus can be considered as having a continuous spectrum. Assuming the apparatus is sensitive in \(\hat{{{{{{{{\bf{y}}}}}}}}}\) direction, we have the result of output effective magnetic field
with By being the measured magnetic field. Therefore the profile of amplification reads
where f0 is the resonant frequency and the relations \({R}_{2}^{{{{{{{{\rm{n}}}}}}}}}={({T}_{2{{{{{{{\rm{n}}}}}}}}})}^{-1}=2\pi \Gamma\) and amplification factor \(q=\frac{1}{2}{\gamma }_{{{{{{{{\rm{Ne}}}}}}}}}{B}_{0}^{{{{{{{{\rm{n}}}}}}}}}{T}_{2{{{{{{{\rm{n}}}}}}}}}\) have been used.
ALP field \({\beta }_{{{{{{{{\rm{y}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}\) couples only to noble-gas nuclear spins. The resulting effective magnetic field \({\beta }_{{{{{{{{\rm{y}}}}}}}}}^{{{{{{{{\rm{eff}}}}}}}}}\) is similar to Eq. (14) without the \({\beta }_{{{{{{{{\rm{y}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}\cos \omega t\) term. The amplification of the ALP field reads
which is the same as the Lorentzian lineshape.
One can easily see that in the large-q limit the amplification profile reduces to the Lorentzian shape. Using Bloch equations, we can obtain the responses to the usual magnetic field and ALP field coupling to noble-gas nuclear spins. We use the usual magnetic field to calibrate the response of the system. The measured response to the usual magnetic field is fitted with Eq. (15). And then we substitute the measured parameters into the response equation for ALP field coupling to noble-gas nuclear spins based on Eq. (16).
Stochastic properties of the signal
In this section, we implement the gradient axion stochastic spectrum derived from ref. 69, while making minor adjustments to account for the new response. In our NMR setup, the comagnetometer is simultaneously sensitive to the anomalous magnetic fields along both \(\hat{x}\) and \(\hat{y}\) directions, which is different from ref. 69 with only one sensitive axis. Therefore, the calculation details need to be adapted accordingly.
The gradient of the axion field within the Galactic halo is assumed to be
where ÏDMâ=â0.4âGeVââ âcmâ362 represents the local DM density. Here, Ïp denotes the axion energy associated with momentum p, αp follows standard Rayleigh distribution as a random variable, and Φp is a random phase uniformly distributed within \(\left[0,2\pi \right]\). The summation runs over all the infinitesimal momentum space Ωp domain.
The DM Maxwell-Boltzmann distribution in the Galactic coordinate system is expressed as69
where vE is the Earth velocity and v represent the DM velocity in the lab frame, thus vEâ+âv corresponds to the axion DM velocity in the rest frame of the Galaxy. The parameter \(\sqrt{2}{\sigma }_{{{{{{{{\rm{v}}}}}}}}}\approx 220\,\,{{\mbox{km}}}/{{\mbox{s}}}\,\) represents the virial velocity.
By choosing the Cartesian coordinate system \(\left\{\hat{u},\hat{{{{{{{{\mathscr{\nu }}}}}}}}},\hat{s}\right\}\), where \(\hat{\nu }\) is parallel to vE and \(\hat{u}\) and \(\hat{s}\) form the other orthonormal basis vectors, and applying certain mathematical transformations, we can derive the expression as follows:
where j is the velocity index and \({f}_{j}\equiv {(2\pi {\sigma }_{{{{{{{{\rm{v}}}}}}}}}^{2})}^{-3/2}\exp \left[-\frac{{v}_{j}^{2}+{v}_{E}^{2}}{2{\sigma }_{{{{{{{{\rm{v}}}}}}}}}^{2}}\right]\). The detailed derivation and the complex definition of ϵi,j are provided in ref. 69.
A bias magnetic field is applied along the laboratory \(\hat{{{{{{{{\rm{z}}}}}}}}}\) axis (outward from the ground), and the response signal of the alkali metal ensemble is detected along the \(\hat{{{{{{{{\rm{x}}}}}}}}}\) axis (east-westward). The response of the alkali atom \({P}_{x}^{e}\) receives contributions from the gradient axion field along the \(\hat{{{{{{{{\rm{x}}}}}}}}}-\hat{{{{{{{{\rm{y}}}}}}}}}\) plane. Based on the argument below Eq. (12), there exists a Ï/2 phase difference between the responses \({{{{{{{{\bf{b}}}}}}}}}_{0}\cdot \hat{{{{{{{{\rm{x}}}}}}}}}\) and \({{{{{{{{\bf{b}}}}}}}}}_{0}\cdot \hat{{{{{{{{\rm{y}}}}}}}}}\), which needs to be addressed in the calculations of ALP signal response variance.
For the case of ALP DM, we have \({b}_{0}^{{{{{{{{\rm{x}}}}}}}}}=\overrightarrow{\nabla }a\cdot {\hat{{{{{{{{\bf{m}}}}}}}}}}_{{{{{{{{\rm{x}}}}}}}}}\) and \({b}_{0}^{{{{{{{{\rm{y}}}}}}}}}=\overrightarrow{\nabla }a\cdot {\hat{{{{{{{{\bf{m}}}}}}}}}}_{{{{{{{{\rm{y}}}}}}}}}\). The sensitive axes \({\hat{{{{{{{{\bf{m}}}}}}}}}}_{{{{{{{{\rm{x}}}}}}}}}\) and \({\hat{{{{{{{{\bf{m}}}}}}}}}}_{{{{{{{{\rm{y}}}}}}}}}\) in the Cartesian coordinate system \(\left\{u,v,s\right\}\) during the experiment can be described as follows:
where \(i\in \left\{u,v,s\right\},{\omega }_{{{{{{{{\rm{E}}}}}}}}}=2\pi \times 1.16\,\mu \,{{\mbox{Hz}}}\,\) is the Earthâs sidereal frequency, and the coefficients \({C}_{i}^{(1)},{C}_{i}^{(2)},{D}_{i}^{(1)},{D}_{i}^{(2)},{\theta }_{i}^{(1)},\) and \({\theta }_{i}^{(2)}\) are determined through coordinate transformations from the laboratory coordinate system to the celestial coordinate system, and subsequently to the Galactic coordinate system. These coefficients are calculated before each individual scanning.
The calculation for the axion-nucleon coupling sensitivity
We conduct a profile analysis using the log-likelihood ratio (LLR) method based on the response power spectrum (V2). In this experiment, we are in the peaked regime, indicating that the power spectrum exhibits a peak at the resonant frequency even without injecting a magnetic field in the xââây plane. This phenomenon arises because the amplified magnetic noise after amplification far surpasses the photon shot noise, which cannot be amplified. Therefore, we cannot assume that the background noise is frequency-independent white noise within the narrow LLR frequency range (as assumed in refs. 27,28,69,70). Instead, we model the background noise as a sum of two zero-mean Gaussian components: one being the non-amplifiable photon shot noise with a variance of Σna, which remains constant in frequency, and the other being the amplifiable magnetic noise with a variance of Σaââ âηF(f)2, following the Fano profile. The axion signal is modeled as a multivariate Gaussian distribution28,69, with a variance of ΣA(fk, fr)â=âgcali(fk)gcali(fr)ÏA(fk, fr).
The calibration function gcali(f) is the transfer function that converts magnetic field strength in Tesla to the apparatus response in Voltage at a given frequency f, which includes the amplification function for the signal. We obtain this function by injecting magnetic fields of known amplitude and varying frequency along the sensitive axes, subsequently measuring the apparatus response. We fit the obtained response using a Fano lineshape given by Eq. (15). From the fit parameters q and Î, we can then derive the calibration function for the ALP field, which is modeled by a Lorentzian lineshape. The complete response of the ALP field is depicted in Fig. 7. For a narrow band near the resonant frequency (e.g., f0â±â0.25âHz), the shape is accurately described by the Lorentzian function. Lastly, the aforementioned calibration procedure is conducted before each measurement.
The covariance matrix ÏA(fk, fr) is defined by Eqs. B33 and B34 in ref. 69:
where μNeâ=â2ÏâÃâ3.36MHzââ âTâ1 is the 21Ne gyromagnetic ratio and
and \({\tilde{\beta }}_{k}\) is the discrete Fourier transform of obtained data time series βn
The definitions of Eik(v) and Fik(v) are slightly modified to account for two sensitive axes and their Ï/2 relative phase in the response of \({P}_{{{{{{{{\rm{x}}}}}}}}}^{{{{{{{{\rm{e}}}}}}}}}\):
with \(\omega ={m}_{a}(1+\frac{1}{2}{v}^{2})\).
Since both the signal and the background are multivariate Gaussian random variables, we can construct the likelihood function as follows:
with vector \({{{{{{{\bf{d}}}}}}}}=\left\{{A}_{k},{B}_{k}\right\}\) and variance matrix \(\Sigma ={\Sigma }_{{{{{{{{\rm{A}}}}}}}}}({g}_{{{{{{{{\rm{aNN}}}}}}}}})+{\Sigma }_{{{{{{{{\rm{na}}}}}}}}}\cdot {\mathbb{1}}+\,{{\mbox{diag}}}\,({\Sigma }_{{{{{{{{\rm{a}}}}}}}}}\cdot {\eta }_{{{{{{{{\rm{F}}}}}}}}}^{2}(f))\)
To set a quantitative limit, we use the LLR defined below
where we treat both Σa and Σna as nuisance parameters. \({\tilde{\Sigma }}_{{{{{{{{\rm{a}}}}}}}}}\) and \({\tilde{\Sigma }}_{{{{{{{{\rm{na}}}}}}}}}\) maximizes the likelihood L for a fixed nucleus coupling gaNN and \({\hat{\Sigma }}_{{{{{{{{\rm{a}}}}}}}}},{\hat{\Sigma }}_{{{{{{{{\rm{na}}}}}}}}}\) and \(\hat{g}\) are unconditional estimators that maximize likelihood L for given d. In order to establish an upper limit, we set LLRâ=â0 if \(\hat{g} \, > \,{g}_{{{{{{{{\rm{aNN}}}}}}}}}\). The defined LLR can be asymptotically shown to follow a half-Ï2 distribution with the cumulative distribution function63
thus the 95% CL upper limits are determined by finding the value of gaNN where LLR(gaNN)â=â2.7055. The effective coupling between axions and the nucleus N is given by gaNNâ=âξngannâ+âξpgapp, where \({\xi }_{{{{{{{{\rm{n}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}=0.58\) and \({\xi }_{{{{{{{{\rm{p}}}}}}}}}^{{{{{{{{\rm{Ne}}}}}}}}}=0.04\) are the spin-polarization fractions for neutrons and protons in 21Ne64,65. When setting limits on the coupling between ALP and neutrons, we set gappâ=â0. Similarly, when setting the limits on the coupling of protons, we set gannâ=â0.
The shape of the sensitivity curve
For measuring times smaller or equal to the coherent time \({\tau }_{a} \sim 2\hslash /({m}_{a}{v}_{{{{{{{{\rm{vir}}}}}}}}}^{2}) \sim 165\,\,{{\mbox{hours}}}\,\left(\frac{\,{{\mbox{Hz}}}\,}{{f}_{a}}\right)\), the signal predominantly falls within a single bin. As discussed in the main text, in this scenario, the lineshape of the signal limit can be estimated.
In each single bin, both the signal and the background are treated as the Gaussian variables in the frequency domain. To establish the upper limit for the signal strength, we make the assumption that the data consists only of background with a certain variance. Therefore, the data power spectrum follows a Ï2 distribution with two degrees of freedom because the variance is complex. The corresponding cumulative distribution function follows an exponential distribution49,
Therefore, the 95% confidence-level limit on the variance contributed by the signal is y95%â=â5.99âÏ2. Consequently, the signal coupling at the frequency bin f should satisfy:
Σa and Σna denote the variances of the amplifiable and non-amplifiable background components, e.g., the magnetic noise and the photon shot noise respectively. c1 represents the signal coefficient independent of f. ηF(f) and η(f) correspond to the Fano profile and Lorentzian profiles, respectively. In this paper, we assume that the axion DM-electron coupling is negligible, resulting in no interference between alkali metal and noble gas responses. Therefore, the signal profile can be described by a Lorentzian shape at the leading order, which is symmetric near the resonant frequency. As illustrated in Fig. 7, the response to the ALP field does show a Lorentzian shape and only deviates in the region far from the resonance, which is different from the usual magnetic field. Since we are primarily interested in a narrow band near the resonant frequency, e.g., f0â±â0.25âHz, the Lorentzian profile is a good description for the signal response. Following Eq. (31), we can derive the sensitivity as a function of frequency gaNN(f) and its FWHM for the signal as
where Îââ¡â(fâââf0)/Î,âq represents the amplification factor at the Larmor frequency, and κââ¡âq2Σa/Σna is the amplification level.
In practice, we obtain values for q and Î from the calibration process, and κ from the experimental data. In Fig. 9, we present the amplification factor q and the amplification level κ for the 49 individual scans, represented by red and blue dots for apparatus A and B, respectively. Apparatus A has an amplification factor of approximately âqâââ50, while Apparatus B has a value of around 100. The amplification factor q is negative because the dip, attributed to the Fano effect, appears on the left side of the resonant frequency, see Fig. 7. As for the amplification level κ, Apparatus B shows a range between 400â1000, whereas Apparatus A mostly falls within the range of 100 and 1000. Since we have observed in the experiment that the photon shot noise is much larger than the magnetic noise (Σnaââ«âΣa), we can make the approximation q2ââ«âκ, as supported by Fig. 9. Using this approximation, the sensitivity \({g}_{{{{{{{{\rm{aNN}}}}}}}}}^{95 \% }(f)\) and its bandwidth ÎFWHM in Eqs. (32) and (33), will revert to the simple Lorentzian expression.
Due to the Fano effect, there is another difference from the simple Lorentzian case: the best sensitivity for \({g}_{{{{{{{{\rm{aNN}}}}}}}}}^{95 \% }(f)\) is not achieved at the resonant frequency f0. Following Eq. (33), we can determine that the best sensitivity \({\tilde{g}}_{{{{{{{{\rm{aNN}}}}}}}}}^{\min }\equiv {g}_{{{{{{{{\rm{aNN}}}}}}}}}^{95 \% }({\tilde{f}}_{0})\) and its corresponding frequency \({\tilde{f}}_{0}\) are
In the second line, the best sensitivity \({\tilde{g}}_{{{{{{{{\rm{aNN}}}}}}}}}^{\min }\) is very close to the Lorentzian value when q2ââ«âκ. The frequency \({\tilde{f}}_{0}\) deviates from its Lorentzian value f0, which is significant and can be observed.
In Fig. 10, we present a summary plot of the sensitivity curves for all 49 individual scans, each smoothed with a ±250 bins moving average. Each sensitivity curve exhibits a peak centered around its resonant frequency, but the peaks vary in width, corresponding to their respective κ values. Their widths are clearly around 0.25âHz, which are much broader than the typical nuclear spin-resonance width for 21Ne. Specifically, there are four points for Apparatus A with κ ~10, which exactly match the sharp peaks in Fig. 10. These points at around 6, 6.5, and 7.5âHz are better than the nearby scans, due to the suppression of the experimental environmental noises.
In addition, our experiment operates within a frequency range of 5â15âHz. We maintained stable experimental conditions throughout the measurements, allowing us to estimate the measurement error of the transfer function gcali at key frequency points (5, 10.75, 11.5, and 13.75âHz). Assuming the error trend follows a linear behavior across these points, we extrapolate it to estimate the systematic error for the entire range.
For each key frequency point, we performed three independent measurements of gcali. The systematic error of our measurements is defined as the error of gcali, and we use it to calculate the standard deviation of neutron coupling gann. In Fig. 11, we show the sensitivity with the error bands reflecting systematic uncertainties from the calibration function.
The scan strategy
As discussed in refs. 43,44,69, when the transverse spin-relaxation time is significantly shorter than both the ALP DM coherent time Ïa and the measurement time T (as is the case in this experiment), the 95% upper limits follow a power law of gâââTâ1/2 and Tâ1/4 for TââªâÏa and Tââ«âÏa, respectively. Since our experiment operates at low frequencies, approximately on the order of \({{{{{{{\mathcal{O}}}}}}}}(10)\) Hz, the coherent time can extend up to around ~30âh. Therefore, we choose the optimal measurement time as TâââÏa for each run, with a cumulative measurement duration of 1440âh for 49 individual scans.
As explained in the previous chapter, given that we are in the magnetic-dominated regime, widening the scan step can be considerably beneficial, potentially leading to a reduction in the number of scans needed. In practice, we opt for a frequency step of 0.25âHz to scan the ALP frequency range between 5â15âHz.
The post-analysis of possible candidates
Since the axion mass is not known as a prior, one has to test all the frequencies in a given range. We explore if there are possible axion DM signals in the data that exclude the background-only hypothesis. For a quantitative analysis of possible axion candidates, we define the following test statistics to test the significance of a best-fit signal compared to the background-only model69,
where \({\tilde{\Sigma }}_{{{{{{{{\rm{a}}}}}}}}}\) and \({\tilde{\Sigma }}_{{{{{{{{\rm{na}}}}}}}}}\) are maximum LLR estimators for gaNNâ=â0, while \(\hat{g},{\hat{\Sigma }}_{{{{{{{{\rm{na}}}}}}}}}\) and \({\hat{\Sigma }}_{{{{{{{{\rm{a}}}}}}}}}\) maximize the LLP across all parameter spaces.
Since we need to test a large number of ALP candidate frequencies within a given frequency range, we have to consider the look-elsewhere effect. In the case of NF independent tests, the global p value, denoted as pglobal, is given by the formula:
where p is the p value of a single measurement. For sufficiently small p, we can approximate it as pglobalâ=â1âââNFp. Utilizing Eq. (29), we can calculate the corresponding qglobal after accounting for the look-elsewhere effect:
To determine the number of independent tests, NF, for the frequencies, we employ the method outlined in ref. 69. NF is given by:
where F represents the test frequency interval for ALP. The function α(v) is determined through Monte-Carlo simulations. Using the above equations, the 5Ï threshold corresponds to qglobal within the range [41.33, 47.36] for 49 individual measurements.
After analyzing all 49 individual scans, several frequency candidates surpass the threshold of the background-only hypothesis, LLRdiscoveryâ=âqglobal, signifying a 5Ï significance with consideration for the look-elsewhere effect. However, each of these candidate frequencies can be tested by multiple adjacent individual scans. Upon closer examination, we observed that these 5Ï anomalous peaks were not consistently present in other individual scans, suggesting they may be transient background fluctuations. Therefore, we conclude that there is no ALP signal in our experiments.
Data availability
All relevant data were available from the corresponding author upon request.
References
Takahashi, H., Kassa, E., Christoforou, C. & Keller, M. Strong coupling of a single ion to an optical cavity. Phys. Rev. Lett. 124, 013602 (2020).
Forn-DÃaz, P., Lamata, L., Rico, E., Kono, J. & Solano, E. Ultrastrong coupling regimes of light-matter interaction. Rev. Mod. Phys. 91, 025005 (2019).
Atalaya, J., Opremcak, A., Nersisyan, A., Lee, K. & Korotkov, A. N. Measurement of small photon numbers in circuit qed resonators. Phys. Rev. Lett. 132, 203601 (2024).
Raffelt, G. G. Astrophysical methods to constrain axions and other novel particle phenomena. Phys. Rept. 198, 1â113 (1990).
Graham, P. W., Irastorza, I. G., Lamoreaux, S. K., Lindner, A. & van Bibber, K. A. Experimental searches for the axion and axion-like particles. Ann. Rev. Nucl. Part. Sci. 65, 485â514 (2015).
Safronova, M. S. et al. Search for new physics with atoms and molecules. Rev. Mod. Phys. 90, 025008 (2018).
Chadha-Day, F., Ellis, J. & Marsh, D. J. E. Axion dark matter: what is it and why now? Sci. Adv. 8, abj3618 (2022).
Semertzidis, Y. K. & Youn, S. Axion dark matter: how to see it? Sci. Adv. 8, abm9928 (2022).
Kimball, D. F. J. & van Bibber, K. (eds.) The Search for Ultralight Bosonic Dark Matter (Springer, 2022).
Sikivie, P. Experimental tests of the invisible axion. Phys. Rev. Lett. 51, 1415â1417 (1983).
De Panfilis, S. et al. Limits on the abundance and coupling of cosmic axions at 4.5â<âmaâ<â5.0 μeV. Phys. Rev. Lett. 59, 839â842 (1987).
Hagmann, C. et al. Results from a search for cosmic axions. Phys. Rev. D 42, 1297â1300 (1990).
CAST Collaboration. New CAST limit on the axion-photon interaction. Nat. Phys. 13, 584â590 (2017).
Brubaker, B. M., Zhong, L., Lamoreaux, S.K., Lehnert, K. W. & van Bibber, K. A. Haystac axion search analysis procedure. Phys. Rev. D 96, 123008 (2017).
Boutan, C. et al. Piezoelectrically tuned multimode cavity search for axion dark matter. Phys. Rev. Lett. 121, 261302 (2018).
Lee, S. et al. Axion dark matter search around 6.7μeV. Phys. Rev. Lett. 124, 101802 (2020).
Jeong, J. et al. Search for invisible axion dark matter with a multiple-cell haloscope. Phys. Rev. Lett. 125, 221302 (2020).
Bartram, C. et al. Axion dark matter experiment: run 1b analysis details. Phys. Rev. D 103, 032002 (2021).
Bartram, C. et al. Dark matter axion search using a josephson traveling wave parametric amplifier. Rev. Sci. Instrum. 94, 044703 (2023).
Backes, K. M. et al. A quantum enhanced search for dark matter axions. Nature 590, 238â242 (2021).
Kwon, O. et al. First results from an axion haloscope at CAPP around 10.7μeV. Phys. Rev. Lett. 126, 191802 (2021).
Adair, C. M. et al. Search for dark matter axions with CAST-CAPP. Nat. Commun. 13, 6180 (2022).
Budker, D., Graham, P. W., Ledbetter, M., Rajendran, S. & Sushkov, A. Proposal for a cosmic axion spin precession experiment (CASPEr). Phys. Rev. X 4, 021030 (2014).
Graham, P. W. & Rajendran, S. New observables for direct detection of axion dark matter. Phys. Rev. D 88, 035023 (2013).
Wu, T. et al. Search for axionlike dark matter with a liquid-state nuclear spin comagnetometer. Phys. Rev. Lett. 122, 191302 (2019).
Garcon, A. et al. Constraints on bosonic dark matter from ultralow-field nuclear magnetic resonance. Sci. Adv. 5, eaax4539 (2019).
Jiang, M., Su, H., Garcon, A., Peng, X. & Budker, D. Search for axion-like dark matter with spin-based amplifiers. Nat. Phys. 17, 1402â1407 (2021).
Bloch, I. M. et al. New constraints on axion-like dark matter using a Floquet quantum detector. Sci. Adv. 8, abl8919 (2022).
Kornack, T., Ghosh, R. & Romalis, M. Nuclear spin gyroscope based on an atomic comagnetometer. Phys. Rev. Lett. 95, 230801 (2005).
Wei, K. et al. Ultrasensitive atomic comagnetometer with enhanced nuclear spin coherence. Phys. Rev. Lett. 130, 063201 (2023).
Terrano, W. A. & Romalis, M. V. Comagnetometer probes of dark matter and new physics. Quantum Sci. Technol. 7, 014001 (2022).
Wei, K. et al. Constraints on exotic spin-velocity-dependent interactions. Nat. Commun. 13, 7387 (2022).
Chaudhuri, S., Irwin, K., Graham, P. W. & Mardon, J. Optimal impedance matching and quantum limits of electromagnetic axion and hidden-photon dark matter searches. Preprint at arXiv:1803.01627 (2018).
Chaudhuri, S., Irwin, K. D., Graham, P. W. & Mardon, J. Optimal electromagnetic searches for axion and hidden-photon dark matter. Preprint at arXiv:1904.05806 (2019).
Liu, H., Elwood, B. D., Evans, M. & Thaler, J. Searching for axion dark matter with birefringent cavities. Phys. Rev. D 100, 023548 (2019).
Berlin, A. et al. Axion dark matter detection by superconducting resonant frequency conversion. JHEP 07, 088 (2020).
Berlin, A., DâAgnolo, R. T., Ellis, S. A. R. & Zhou, K. Heterodyne broadband detection of axion dark matter. Phys. Rev. D 104, L111701 (2021).
Brady, A. J. et al. Entangled sensor-networks for dark-matter searches. PRX Quantum 3, 030333 (2022).
Brouwer, L. et al. Proposal for a definitive search for GUT-scale QCD axions. Phys. Rev. D 106, 112003 (2022).
Brouwer, L. et al. Projected sensitivity of DMRadio-m3: a search for the QCD axion below 1 μeV. Phys. Rev. D 106, 103008 (2022).
Brady, A. J. et al. Entanglement-enhanced optomechanical sensor array with application to dark matter searches. Commun. Phys. 6, 237 (2023).
Benabou, J. N., Foster, J. W., Kahn, Y., Safdi, B. R. & Salemi, C. P. Lumped-element axion dark matter detection beyond the magnetoquasistatic limit. Phys. Rev. D 108, 035009 (2023).
Dror, J. A., Gori, S., Leedom, J. M. & Rodd, N. L. Sensitivity of spin-precession axion experiments. Phys. Rev. Lett. 130, 181801 (2023).
Zhang, Y. et al. Frequency-scanning considerations in axionlike dark matter spin- precession experiments. Ann. Phys. 536, 2300223 (2023).
Kornack, T., Smullin, S., Lee, S.-K. & Romalis, M. A low-noise ferrite magnetic shield. Appl. Phys. Lett. 90, 223501 (2007).
Schive, H.-Y., Chiueh, T. & Broadhurst, T. Cosmic structure as the quantum interference of a coherent dark wave. Nat. Phys. 10, 496â499 (2014).
Centers, G. P. et al. Stochastic fluctuations of bosonic dark matter. Nat. Commun. 12, 7321 (2021).
Gramolin, A. V. et al. Spectral signatures of axionlike dark matter. Phys. Rev. D 105, 035029 (2022).
Walck, C. Hand-book on Statistical Distributions for Experimentalists (1996).
Asztalos, S. J. et al. Design and performance of the ADMX SQUID-based microwave receiver. Nucl. Instrum. Meth. A 656, 39â44 (2011).
Du, N. et al. A search for invisible axion dark matter with the axion dark matter experiment. Phys. Rev. Lett. 120, 151301 (2018).
Khatiwada, R. et al. Axion dark matter experiment: detailed design and operations. Rev. Sci. Instrum. 92, 124502 (2021).
Brubaker, B. M., Zhong, L., Lamoreaux, S. K., Lehnert, K. W. & van Bibber, K. A. HAYSTAC axion search analysis procedure. Phys. Rev. D 96, 123008 (2017).
Alesini, D. et al. Search for galactic axions with a high-Q dielectric cavity. Phys. Rev. D 106, 052007 (2022).
Cervantes, R. et al. Deepest sensitivity to wavelike dark photon dark matter with SRF cavities. United States. https://www.osti.gov/servlets/purl/1882569 (2022).
Tang, Z. et al. SRF cavity searches for dark photon dark matter: first scan results. Preprint at arXiv:2305.09711 (2023).
Babcock, E. et al. Hybrid spin-exchange optical pumping of 3-He. Phys. Rev. Lett. 91, 123003 (2003).
Romalis, M. Hybrid optical pumping of optically dense alkali-metal vapor without quenching gas. Phys. Rev. Lett. 105, 243001 (2010).
Kornack, T. & Romalis, M. Dynamics of two overlapping spin ensembles interacting by spin exchange. Phys. Rev. Lett. 89, 253002 (2002).
Ghosh, R. K. & Romalis, M. V. Measurement of spin-exchange and relaxation parameters for polarizing ne 21 with k and rb. Phys. Rev. A 81, 043415 (2010).
Padniuk, M. et al. Response of atomic spin-based sensors to magnetic and nonmagnetic perturbations. Sci. Rep. 12, 324 (2022).
de Salas, P. F. & Widmark, A. Dark matter local density determination: recent observations and future prospects. Rept. Prog. Phys. 84, 104901 (2021).
Cowan, G., Cranmer, K., Gross, E. & Vitells, O. Asymptotic formulae for likelihood-based tests of new physics. Eur. Phys. J. C 71, 1554 (2011).
Brown, B. A., Bertsch, G. F., Robledo, L. M., Romalis, M. V. & Zelevinsky, V. Nuclear matrix elements for tests of local Lorentz invariance violation. Phys. Rev. Lett. 119, 192504 (2017).
Almasi, A., Lee, J., Winarto, H., Smiciklas, M. & Romalis, M. V. New limits on anomalous spin-spin interactions. Phys. Rev. Lett. 125, 201802 (2020).
Miroshnichenko, A. E., Flach, S. & Kivshar, Y. S. Fano resonances in nanoscale structures. Rev. Modern Phys. 82, 2257 (2010).
Jiang, M. et al. Enhanced quantum sensing with amplification and deamplification. Preprint at arXiv:2309.00177 (2023).
Ledbetter, M. P., Savukov, I. M., Acosta, V. M., Budker, D. & Romalis, M. V. Spin-exchange-relaxation-free magnetometry with Cs vapor. Phys. Rev. A 77, 033408 (2008).
Lee, J., Lisanti, M., Terrano, W. A. & Romalis, M. Laboratory constraints on the neutron-spin coupling of feV-scale axions. Phys. Rev. X 13, 011050 (2023).
Wei, K. et al. Dark matter search with a strongly-coupled hybrid spin system. Preprint at arXiv: 2306.08039 (2023).
Bartram, C. et al. Search for invisible axion dark matter in the 3.3â4.2 μeV mass range. Phys. Rev. Lett. 127, 261803 (2021).
Yang, B., Yoon, H., Ahn, M., Lee, Y. & Yoo, J. Extended axion dark matter search using the CAPP18T haloscope. Phys. Rev. Lett. 131, 081801 (2023).
Karanth, S. et al. First search for axionlike particles in a storage ring using a polarized deuteron beam. Phys. Rev. X 13, 031004 (2023).
Ji, W. et al. New experimental limits on exotic spin-spin-velocity-dependent interactions by using smco5 spin sources. Phys. Rev. Lett. 121, 261803 (2018).
The data pertaining to axion-nucleon coupling limits from ChangE-NMR. https://github.com/JiaLiu2023/Axion_Data_ChangE_NMR (2023).
Buschmann, M., Dessert, C., Foster, J. W., Long, A. J. & Safdi, B. R. Upper limit on the QCD axion mass from isolated neutron star cooling. Phys. Rev. Lett. 128, 091102 (2022).
Acknowledgements
The work of K.W. is supported by NSFC under Grants No. 62203030 and 61925301 for Distinguished Young Scholars and in part by the Innovation Program for Quantum Science and Technology under Grant 2021ZD0300401. The work of J.L. is supported by NSFC under Grant No. 12075005 and 12235001. The work of X.-P.W. is supported by the National Science Foundation of China under Grant No. 12005009 and the Fundamental Research Funds for the Central Universities. The work of W.J. and D.B. is supported by the DFG Project ID 390831469: EXC 2118 (PRISMA+ Cluster of Excellence), by the COST Action within the project COSMIC WISPers (Grant No. CA21106), and by the QuantERA project LEMAQUME (DFG Project No. 500314265).
Author information
Authors and Affiliations
Contributions
K.W., W.J., J.Liu, X-P.W., and D.B. proposed the experimental concept and devised the experimental protocols, and wrote the manuscript. Z.X., K.W., X.Heng, and X.Huang designed and performed experiments, analyzed the data, and wrote the manuscript. X.M., Y.H., T.A., J.Liao, and J.Liu analyzed the data, added theoretical discussion, and edited the manuscript. All authors contributed with discussions and checking the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Communications Physics thanks the anonymous reviewers for their contribution to the peer review of this work.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Xu, Z., Ma, X., Wei, K. et al. Constraining ultralight dark matter through an accelerated resonant search. Commun Phys 7, 226 (2024). https://doi.org/10.1038/s42005-024-01713-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s42005-024-01713-7