Abstract
Axions, originally proposed to solve the strong CP problem of quantum chromodynamics, emerge now as leading candidates of WISP dark matter. The rich phenomenology associated to the light and stable QCD axion can be described as an effective magnetic field that can be experimentally investigated. For the QUAX experiment, dark matter axions are searched by means of their resonant interactions with electronic spins in a magnetized sample. In principle, axion-induced magnetization changes can be detected by embedding a sample in an rf cavity in a static magnetic field. In this work we describe the operation of a prototype ferromagnetic haloscope, with a sensitivity limited by thermal fluctuations and receiver noise. With a preliminary dark matter search, we are able to set an upper limit on the coupling constant of DFSZ axions to electrons \(g_{aee}<4.9\times 10^{-10}\) at 95% C.L. for a mass of \(58\,\upmu \)eV (i. e. 14 GHz). This is the first experimental result with an apparatus exploiting the coupling between cosmological axions and electrons.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
A major fraction of the mass content of the universe is composed of dark matter (DM), i.e. particles not interacting significantly with electromagnetic radiation, with ordinary matter or self-interacting (cold dark matter) [1,2,3]. Up-to-date results [4] show that with respect to the universe critical density the DM fraction is the \(25.8\%\) while the luminous matter fraction is \(5.7\%\), meaning that DM is about five times more abundant than ordinary baryonic matter. This outstanding result triggered theoretical studies aiming to understand the nature of DM, for instance in the form of new particles beyond the Standard Model (SM).
The axion is a good candidate for DM but was not originally introduced to account for this specific issue. To solve the strong CP problem Peccei and Quinn added a new symmetry to the SM [5], which breaks at an extremely high energy scale \(F_a\) producing a pseudo-Goldstone boson, the axion [6]. Among the proposed models, the “invisible axion” model classes KSVZ and DFSZ still hold [7,8,9,10]. For scales \(F_a\sim 10^{12}\,\)GeV, corresponding to typical mass values \(m_a\lesssim 1\,\)meV, large quantities of axions may have been produced in the early universe and could account even for the totality of cold dark matter [11]. Consequently, several detection schemes have been devised during the last decades to search for relic axions. The value of \(F_a\) is not fixed by the theory, however, cosmological considerations and astrophysical observations [12,13,14,15,16,17] provide boundaries on \(F_a\) and suggest a favoured axion mass range \(1\,\upmu \mathrm {eV}<m_a<10\,\mathrm {meV}\). In addition, lattice results on QCD topological susceptibility, based on reliable computations of the axion relic density, indicate a preferred window for the axion mass in the range of tens of \(\upmu \)eV [18,19,20,21,22,23].
Axion model classes can be tested with different experimental techniques [24,25,26,27,28]. Most of these experiments are based on the Primakoff effect, i.e. an axion to photon conversion in a strong static magnetic field [29,30,31,32,33,34,35,36]. In particular, the ADMX experiment reached the cosmologically relevant sensitivity to exclude the axion mass range \(1.9\,\upmu \mathrm {eV}\lesssim m_a\lesssim 3.7\,\upmu \mathrm {eV}\) for the KSVZ model and \(2.66\,\upmu \mathrm {eV}<m_a<2.81\,\upmu \mathrm {eV}\) for the DFSZ model [37], assuming a local DM density of 0.45 GeV/cm\(^3\). On the other hand the axion-fermion coupling, explicitly predicted in different axion models including DFSZ [38,39,40,41,42], allows for designing new detectors that exploit the interaction between axions and fermionic spins [43,44,45]. Among these, the QUAX detector [46] takes advantage of the resonant interaction between relic axions and a magnetized magnetic sample housed in a microwave cavity. In this paper we present results on the operation of a QUAX demonstrator, based on 5 GaYIG (Gallium Yttrium Iron Garnet) 1 mm diameter spheres placed in a 14 GHz resonant cavity. The apparatus is operated at cryogenic temperatures and its sensitivity is limited only by thermal effects. Section 2 describes the proposed detection scheme, Sects. 3 and 4 report on the measurement of an upper limit on the axion interaction with electronic spins, using the small-scale prototype of the final apparatus. Conclusions are eventually drawn in Sect. 5.
2 Axion detection by resonant interaction with electron spin
The first ideas on axion detection via their conversion to magnons, collective excitations of the spins in a ferromagnet, were discussed in Refs. [47,48,49,50]. As the DFSZ axion and other axion models [38,39,40,41,42] does not suppress the coupling between an axion a and an electron \(\psi \) at the tree level, the Lagrangian reads
where \(\hbar \) is the reduced Planck constant, \(\gamma ^\mu \) is the Dirac matrices vector, m is the mass of the electron and c is the speed of light. The second term of Eq. (1) describes the interaction between a and the spin of the fermion, proportional to the dimensionless coupling constant \(g_{aee}\). In the non-relativistic limit, the interaction term can be expressed as a function of the Bohr magneton \(\mu _B\) and of the effective axionic field \(\mathbf {B}_a\)
where \(\hat{\sigma }\) is the Pauli matrices vector and e is the charge of the electron.
Due to the Earth motion through the DM halo of the Galaxy, relic axions can be seen as a wind in an Earth-based laboratory, thus a non zero value of \(\nabla a\) is expected. The DM wind average speed is \(v_a\simeq 220\,\)km/s with a dispersion of about 270 km/s [51]. Axions will interact with an electron spin as an effective magnetic field pointing roughly in the direction of Vega [52, 53]. The effective field frequency \(f_a\) and amplitude \(B_a\) are determined by the mass of the axion \(m_a\) and the coupling constant \(g_{aee}=3\times 10^{-11}(m_a/1\,\mathrm {eV})\)
meaning that \(\mathbf {B}_a\) is an extremely weak effective rf magnetic field with a linewidth of \(\varDelta f_a=7.0\,(m_a/58.5\,\upmu \mathrm {eV})\,\)kHz, due to the dispersion of \(v_a\). The axion occupation number is \(n_a=\varrho _\mathrm {dm}/m_a\), where \(\varrho _\mathrm {dm}=0.45\) GeV/cm\(^3\) is the local DM density [4]. For a reference mass \(m_a=58.5\,\upmu \)eV the mean de Broglie wavelength is \(\lambda _{\nabla a}=0.74\lambda _a=0.75h/(m_av_a)=5.1\,\mathrm {m}\), while the coherence time is \(\tau _{\nabla a }=0.68\tau _a=58\,\upmu \)s.Footnote 1
Placing a sample in a static magnetic field \(\mathbf {B}\), perpendicular to the axion wind, it is possible to tune the Larmor frequency of the electrons \(f_L\) to \(f_a\), for \(\mathbf {B}_0=(0,0,B_0)\), the direction of the electron spin \(\hat{\sigma }\) is along the z-axis. The axionic field \(\mathbf {B}_a\), acting on the spins of matter, deposits in the material an amount of power \(P_\mathrm {in}\)
where \(V_s\) is the volume of the material, M its magnetization, \(n_s\) its spin density, \(\gamma \) is the gyromagnetic ratio of the electron, and \(\tau _\mathrm {min}\) the minimum relaxation time of the system. The absorbed power is then re-emitted in the form of rf radiation, which can be collected and represents our axion signal. In free space \(\tau _\mathrm {min}\) is mainly determined by radiation damping mechanisms (i.e. magnetic dipole emission of the sample) [54,55,56], with values much smaller than the material relaxation times. To avoid the radiation damping issue and thus increase the sensitivity, the magnetic sample is placed inside a resonant cavity. A cavity mode with frequency \(f_c \simeq f_L\) couples to the Kittel mode (uniform spin precession) of the material, hybridization takes place and the single cavity mode splits into two hybrid modes with frequencies \(f_-\) and \(f_+\) (strong coupling regime) [57,58,59,60]. This phenomenon limits the phase space of the dipole emission avoiding radiation damping, and is described by the anti-crossing curve represented in Fig. 1, which also justifies the strong coupling regime approximation. The coupling between the cavity mode and the Kittel mode is
where \(\mu _0\) is vacuum magnetic permeability and \(V_m=\xi V_c\) is the product of the cavity volume \(V_c\) and a mode-dependent form factor \(\xi \). The linewidths of the hybrid modes \(k_{+,-}\) are an average of the linewidth of the cavity \(k_c\) and of the material \(k_m\), i. e. \(k_{+,-}=\frac{1}{2}(k_c+k_m)\equiv k_h\). The calculated power spectral density of an empty cavity and of a cavity with the volume \(V_s\) and \(5V_s\) of material are shown in Fig. 2. The two hybrid modes are more sensitive to the power deposited by the axion field since they are not affected by radiation damping, the minimum relaxation time is \(\tau _\mathrm {min}=\min (\tau _h,\tau _{\nabla a})\), where \(\tau _h=1/k_h\). With an antenna critically coupled to one of the hybrid resonant modes, the extracted power is \(P_\mathrm {out}=P_\mathrm {in}/2\). The scalar product \(\hat{\sigma }\cdot \nabla a\) of Eq. (2) shows that the effect is directional. Due to earth rotation, an earth-based experiment experiences a full daily modulation of the signal, due to the variation of the axion wind direction.
3 The QUAX prototype
To implement the scheme presented in Sect. 2 we use a cylindrical copper cavity TM110 mode with resonance frequency \(f_c\simeq 13.98\,\)GHz and linewidth \(k_c/2\pi \simeq 400\,\)kHz at liquid helium temperature, measured with a critically coupled antenna. The shape of the cavity is not a regular cylinder, two symmetric sockets are carved into the cylinder to remove the angular degeneration of the normal mode, the maximum and minimum diameters are 26.7 mm and 26.1 mm, and the length is 50.0 mm. The shape of the cavity and of the mode magnetic field are shown in Fig. 3. The choice of the TM110 mode has the advantage of having a uniform maximum magnetic rf field along the cavity axis. Its volume can be increased just using a longer cavity without changing the mode resonance frequency. For this mode we calculate a form factor \(\xi =0.52\) [61]. The cavity mode is coupled to a magnetic material, thus we studied the properties of several paramagnetic samples and some ferrites. Highest values of \(n_s\) together with long relaxation times have been found for YIG (Yttrium Iron Garnet) and GaYIG (Gallium doped YIG). To avoid inhomogeneous broadening of the linewidth due to geometrical demagnetization, these garnets are shaped as highly polished spheres. Five GaYIG spheres of 1 mm diameter have been placed in the maximum magnetic field of the mode, which lies on the axis of the cavity. The spheres are housed inside a PTFE support large enough to let them rotate in all possible directions, in order to automatically align the GaYIG magnetization easy axis with the external magnetic field.
The amplitude of an external magnetic field \(B_0\) determines the Larmor frequency of the electrons. The uniformity of \(B_0\) on all the spheres must be enough to avoid inhomogeneous broadening of the ferromagnetic resonance. To achieve a magnetic field uniformity \(\le 1/Q_h\), where \(Q_h\sim 10^4\) is the quality factor of the hybrid mode, we make use of a superconducting NbTi cylindrical magnet equipped with a concentric cylindrical NbTi correction magnet. With \(B_0=0.5\,\)T we have \(f_L\simeq f_c\) and thus the hybridization of the cavity and Kittel modes, as discussed in Sect. 2. The power supply of the main magnet is a high-precision, high-stability current generator, injecting 15.416 A into the magnet with a precision better than 1 mA, while a stable current generator provides 26.0 A for the correction magnet. A simplified scheme of the cavity, material and magnet setup is represented in the left part of Fig. 4.
In the strong coupling regime, the hybrid mode frequencies are \(f_+=14.061\,\)GHz and \(f_-=13.903\,\)GHz, yielding a splitting \(g=158\,\)MHz. The coupling g scales exactly with \(\sqrt{n_sV_s}\), in fact \(g=\sqrt{5}\delta \), where \(\delta \simeq 71\,\)MHz is the measured splitting due to a single sphere. This means that all the spins are coherently participating to the material-cavity mode, and ensures that all the spheres magnetization easy axes are aligned along \(B_0\). We use g to calculate the effective number of spins in the sample using the relation described by Eq. (5), we obtain \(n_s=2.13\times 10^{28}\,\mathrm {m}^{-3}\). The weakly coupled linewidth is \(0.7\,\)MHz, yielding a critically coupled one of \(k_+/2\pi =1.4\,\)MHz, corresponding to the hybrid modes relaxation times \(\tau _-\simeq \tau _+=0.11\,\upmu \)s.
The detection electronics consists in an amplification chain which has two inputs, called Input Channel 1 and 2, (IC-1 and IC-2, respectively). Channel 1 measures the signal power, while Channel 2 has calibration and characterization purposes. A cryogenic switch is used to select the desired channel:
- IC-1:
-
the rf power inside the cavity is collected with a dipole antenna whose coupling to the cavity can be changed using an external micro-manipulator, allowing us to switch continuously from sub-critical to over-critical coupling. For optimal measurement conditions, we tune the antenna to critical coupling by doubling the sub-critical linewidth of the selected mode;
- IC-2:
-
a 50 \(\varOmega \) termination \(R_J\), enclosed in a copper block together with a heater resistance, is used as Johnson noise source. The emitted power can be used to calibrate the noise temperature of the system and the total gain, detailed in Sect. 3.1.
The detection electronics, as shown in Fig. 4, is divided into a liquid helium temperature part (LTE) and a room temperature part (HTE). The collected power is amplified by a HEMT cryogenic low-noise amplifier (A1) with gain \(G_\mathrm {A1}\simeq 38\,\)dB. To avoid the back-action noise of the amplifier, a cryogenic isolator with 18 dB of isolation is inserted in the chain. The HTE consists of a room temperature FET amplifier (A2), with \(G_\mathrm {A2}\simeq 34\,\)dB, followed by an IQ mixer used to down-convert the signal with a local oscillator (LO).
The energy distribution of DM axions is highly peaked (\(Q_a \sim 10^6\)) around the actual axion mass, so the corresponding frequency distribution can significantly overlap with one (e. g. \(f_+\)) of the two hybrid modes. As the axion mass is unknown, both modes could be monitored independently in order to double the frequency scan rate of the detector. In our simplified scheme we choose to work only with \(f_+\), thus setting the LO frequency to \(f_\mathrm {LO}=f_+-0.5\,\)MHz and its amplitude to 12 dBm. The antenna output at the hybrid mode frequency is down-converted in the 0 - 1 MHz band, allowing us to efficiently digitize the signal. The phase and quadrature outputs are fed to two low frequency amplifiers (A3\(_{I,Q}\)), with a gain of \(G_3\simeq 50\,\)dB each, and are acquired by a 16 bit ADC sampling at 2 MHz.
A weakly coupled dipole antenna is used to inject low power signals and make transmission measurements of the system using a source oscillator, SO. All the apparatus devices are referenced to a GPS disciplined, oven controlled, local oscillator. The cryogenic part of the apparatus is enclosed in a vacuum vessel immersed in liquid helium, as shown schematically in Fig. 4. Measurements are performed at temperatures \(T_c\sim T_a\simeq 5.0\,\)K and \(T_r\simeq 5.5\,\)K, as read by the cavity, amplifier and \(R_J\) thermometers, respectively.
3.1 Calibration and measurements
For the calibration of the system, the load \(R_J\) is heated to a temperature \(T_r\), as described in Sect. 3. Using IC-2 it is possible to measure the Johnson noise of \(R_J\) in the temperature range \(5 \div 25\,\)K without significantly heating other parts of the apparatus. The rf power from IC-2 is
where \(k_B\) is the Boltzmann constant, \(T_n\) is the noise temperature of the system and \(\varDelta f\) is the bandwidth. By gradually increasing \(T_r\) we linearly change the measured power level of Eq. (6) to determine the noise temperature and gain of the detection electronics, similarly to what is usually obtained with the Y-factor method [62]. In Fig. 5, four collected points are fitted with \(q(T)=aT+b\) to obtain the noise temperature \(T_n=-b/a\) and the total gain \(G_\mathrm {tot}=a\). The error of the estimated parameters is less than 1%. Using this method we verified that the noise of the system changes linearly with the temperature, and that the measured cavity output power at the temperature \(T_c\) is given by Eq. (6) with \(T_r=T_c\), assuming that IC-1 and 2 have the same losses, which is true within 0.2 dB. Typical measured values are \(T_n=9\div 11\,\)K and \(G_\mathrm {tot}=106\div 108\,\)dB, at different frequencies around 14 GHz. This procedure ensures the accuracy of the measurement and then, using IC-1, we perform measurements on the hybrid system with the calibrated electronic chain.
Multiple measurements of the effective axion field have been performed as follows. The vacuum vessel containing the system is cooled down to liquid helium temperature and when a proper thermalization is achieved the detection electronics parameters \(G_\mathrm {tot}\) and \(T_n\) are measured through IC-2. Then we switch to IC-1, set the magnetic field \(B_0\) to \(0.5\,\)T to hybridize the cavity and Kittel mode at \(f_c\simeq f_L\simeq 14\,\)GHz, and critically couple the antenna with the \(f_+\) hybrid mode using the manipulator. A dedicated DAQ software is used to control the oscillators and the ADC, and verifies the correct positioning of the LO with an automated measurement of the hybrid mode transmission spectrum. The ADC digitizes the time-amplitude down-converted signal coming from A3\(_I\) and A3\(_Q\) and the DAQ software stores collected data binary files of 5 s each. The software also provides a simple online diagnostic, extracting 1 ms of data every 5 s, and showing its 512 bin FFT together with the moving average of all FFTs.
As seen in Sect. 2, the axion wind releases a faint power in a band of \(\sim 7\) kHz around \(f_a\). This signal can be seen only if \(f_a\) falls into the detection bandwidth, which corresponds to the linewidth of the hybrid peak. The expected noise power is given by
calculated from Eq. (6) using the data collected from IC-2. Considering the losses of the system and the gain of the amplifiers,we will show that the mean of the measured power is indeed compatible with the expected noise.
3.2 Analysis and results
The signal is down-converted in its in-phase and quadrature components \(\{\phi _n\}\) and \(\{q_n\}\), with respect to the local oscillator, that are sampled separately. We applied a complex FFT to \(\{s_n\}=\{\phi _n\}+i\{q_n\}\) to get its power spectrum \(s^2_\omega \) with positive frequencies for \(f>f_\mathrm {LO}\) and negative frequencies for \(f<f_\mathrm {LO}\). In our experimental settings, the axion signal is mapped almost completely onto the positive frequencies since the hybrid mode linewidth is of order 1 MHz.
Figure 6 reports the analysis of RUN31, which we describe hereafter in some details. The \(\sim \)2.3 h of the measurement consist in 2048000 FFTs of 8192 bins each (frequency resolution of 244 Hz), which were square averaged and rebinned to the bandwidth \(\varDelta f=7.8\,\)kHz (256 bins), close to \(\varDelta f_a\). As explained, we consider only the positive part of the spectrum, consisting of 128 bins, and then calibrate \(s^2_\omega \) using Eq. (7). Some frequency intervals of the power spectrum were affected by disturbances at the ADC output, and has been ignored in the analysis procedure. A polynomial of degree 5 is fitted to the averaged spectrum and the residuals estimated. The averaged spectrum is reported in Fig. 6 together with the fitting function.
In Fig. 6 a plot of the residuals and their histogram is also given. The average value of the residuals is \(-4.6\times 10^{-23}\,\)W with standard deviation \(\sigma _P=2.2\times 10^{-22}\,\)W. The result is compatible with Dicke radiometer equation
where t is the total integration time and \(T_D=T_c+T_n\). This means that the standard deviation of the noise decreases as \(1/\sqrt{t}\) trend at least within the RUN31 time span.
The stability of \(f_+\) is monitored by injecting with SO an rf probe signal at \(f_+-0.9\,\mathrm {MHz}=f_\mathrm {LO}+0.1\,\mathrm {MHz}\). The transmitted amplitude of the probe peak is a monitor of the hybrid peak frequency since it changes if \(f_+\) drifts. Such amplitude is registered during the whole measurement and is plotted in Fig. 7, the frequency stability of this run was around 3.5%, which was enough for the purposes of this measurement. We presume that this variation is mostly due to drifts of the external static magnetic field, since with a \(B_0=0\) run the corresponding variation was much smaller.
To increase the confidence and the consistency of our estimators, additional offline tests have been performed on the acquired data. Firstly, to search for \(P_\mathrm {in}\) when distributed into two adjacent bins, the analysis procedure was repeated using a binning shifted of \(\varDelta f/2\). This test confirmed the reported result. Secondly, we calculate the residuals of the averaged spectra for each 5 s data segment verifying that there are no outliers.
To infer the axion sensitivity of our measurement, two corrections have to be introduced: (i) a loss of \(0.98\,\mathrm {dB}\) (a factor 0.8) at the cavity antenna due to imperfect matching between cavity and axion field [63]; (ii) a factor 1/2 to account for the binning search procedure. In fact, the collected power in a single bin results in \(P_\mathrm {in}/2\) because our resolution bandwidth \(\varDelta f\) overlaps partially with the axion distribution. The correct power standard deviation results \(\sigma _P^\prime =2\sigma _P/0.8\).
The measured rf power is compatible with the modeled noise for every bin and no statistically significant signal consistent with axions was found. The upper limit at the 95% C.L. is \(2\sigma ^\prime _P=1.1\times 10^{-21}\,\)W. This value can be converted to equivalent axion field with the help of Eq. (4), obtaining
where all the reported parameters have been explicitly measured. The limit holds for the central frequency of the hybrid mode, while for other frequencies the sensitivity have to be normalized: the correct sensitivity is obtained dividing \(B_m\) by the normalized amplitude of the hybrid mode Lorentzian. Several measurements have been performed for different cool downs of the setup. Probably due to mechanical instabilities and to the low resolution of the correction magnet power supply, the resulting working frequency \(f_+\) slightly changed between the runs, allowing us to perform also a limited frequency scan over a \(\sim \)3 MHz range. The maximum integration time for a 1 MHz band was 6 h, and no deviations from the \(1/\sqrt{t}\) scaling of \(\sigma _P\) were found.
4 Discussions
Our results represent also a limit on the axion–electron coupling constant. Since \(B_m\) depends on \(g_{aee}\), the explicit form of the effective magnetic field given in Eq. (3) can be recast with the help of Eq. (4), to
at 95% confidence level. The results of this preliminary measurements are far from the sensitivity requirements for a cosmological axion search [see Eq. (3)], however they can be used to detect DM axion-like particles (ALPs), which can account for the whole dark matter density [64]. During the measurement time the DM-wind amplitude was on the maximum of the daily modulation, allowing us to use the collected data to obtain an upper limit on the ALP-electron coupling at the maximum sensitivity. Through Eq. (10) we are able to exclude values of the ALP-electron coupling constant for ALP masses given by \(f_+\) through Eq. (3).
By repeating the analysis procedure described in Sect. 3.2 for seven measurement runs and averaging together overlapping bandwidths, we produce the plot in Fig. 8. The minimum measured value of \(g_{aee}\) is \(4.9\times 10^{-10}\), corresponding to an equivalent axion field limit of \(1.6\times 10^{-17}\,\)T.
4.1 Improvements and discovery potential
To push the present sensitivity towards smaller values of the coupling constant \(g_{aee}\), several improvements should be implemented. In fact, using Eq. (4), the power released by a DFSZ-axion wind in the five GaYIG spheres of our prototype is
corresponding to a rate \(r_a\sim 10^{-10}\,\mathrm {Hz}\) of 14 GHz photons, which is clearly not detectable. To have a statistically significant signal within a reasonable integration time it is mandatory to increase the signal rate, for example in the mHz range, that will give tens of counts per day. The present sensitivity to the power deposited in the system by the axion wind maintains an excess photon rate of order \(100\,\)photons/s.
Short term improvements that will be installed in the prototype include a larger volume of narrow-linewidth magnetic material, namely 10 YIG spheres of 2 mm diameter, a lower cavity temperature with dilution refrigeration and the use of a Josephson Parametric Amplifier (JPA). The increased volume will enhance the axionic signal of a factor 16. As for the noise reduction, a working temperature of 100 mK will reduce the thermal fluctuations and there are hints suggesting that it can also reduce the YIG linewidth. Ultra cryogenic temperatures allow us to use JPAs as first-stage amplifier to drastically increase the sensitivity, since its noise temperature can be of the order of 100 mK. The upgraded prototype should be capable of setting a limit on the effective magnetic field \(B_m\) two orders of magnitude better than the present one.
To achieve the QUAX goal [46], the detector requires an improvement of three to four more orders of magnitude in sensitivity, which can be obtained increasing the signal and reducing the noise. Using a \(V_s\simeq 0.1\,\)liters and \(\tau _\mathrm {min}\simeq 1\,\upmu \)s, the axionic power deposited in the system is \(\sim 10^{-27}\,\)W. This power is smaller than the quantum noise, meaning that a quantum counter, immune to such noise, must be exploited to push the sensitivity to the axion level [77, 78].
To scan different axion masses we must vary the working frequency of the haloscope. A large tuning can be achieved by changing both the cavity mode resonance frequency and the Larmor frequency (i. e. the static magnetic field \(B_0\)). A small frequency tuning is possible by varying only \(B_0\): in this case, a scanning of several MHz is possible without a significant reduction of the sensitivity.
In the favored case of a signal detection, its nature can be systematically studied by QUAX. Since the axion signal is persistent, it will be possible to infer DM properties by using the directionality of the apparatus. Moreover, this setup is able to test different axion models, measuring separately the axion-to-photon and axion-to-electron couplings. In fact the apparatus has also the features of a Sikivie haloscope [24], and can be sensitive to the axion–photon coupling by using a suitable cavity mode.
5 Conclusions
We described the operation of a prototype of the QUAX experiment, a ferromagnetic haloscope sensitive to DM axion through their interaction with electron spin. Our findings indicate the possibility of performing electron spin resonance measurements of a sizable quantity of material inside a cavity cooled down to cryogenic temperatures. By using low noise electronics we search for extra power injected in the system that could be due to DM axions. We reach a power sensitivity of \(10^{-22}\,\)W that can be translated to an upper limit on the the coupling constant \(g_{aee}<4.9\times 10^{-10}\) for an axion mass of \(58\,\upmu \)eV, which, to our knowledge, is the first measurement of the coupling between cosmological axions and electrons. The sensitivity of our apparatus is presently limited only by the noise temperature of the system and thermodynamic fluctuations, as it reaches the limit of Dicke radiometer equation. The overall behavior of the apparatus is as expected, and thus we are confident that the planned upgrades will be effective.
Change history
08 October 2018
The original version of this article unfortunately contained a mistake.
Notes
The numerical factors account for the differences between a and \(\nabla a\), see [46] for further details.
References
F. Zwicky, Helv. Phys. Acta 6, 110 (1933)
V.C. Rubin, N. Thonnard, W.K. Ford Jr., Astrophys. J. 225, L107 (1978). https://doi.org/10.1086/182804
V.C. Rubin, W.K. Ford Jr., N. Thonnard, Astrophys. J. 238, 471 (1980). https://doi.org/10.1086/158003
C. Patrignani et al., Chin. Phys. C 40(10), 100001 (2016). https://doi.org/10.1088/1674-1137/40/10/100001
R.D. Peccei, H.R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). https://doi.org/10.1103/PhysRevLett.38.1440
S. Weinberg, Phys. Rev. Lett. 40(4), 223 (1978)
M. Dine, W. Fischler, M. Srednicki, Phys. Lett. B 104(3), 199 (1981). https://doi.org/10.1016/0370-2693(81)90590-6. http://www.sciencedirect.com/science/article/pii/0370269381905906
M. Shifman, A. Vainshtein, V. Zakharov, Nucl. Phys. B 166(3), 493 (1980). https://doi.org/10.1016/0550-3213(80)90209-6. http://www.sciencedirect.com/science/article/pii/0550321380902096
J.E. Kim, Phys. Rev. Lett. 43, 103 (1979). https://doi.org/10.1103/PhysRevLett.43.103
M. Dine, W. Fischler, Physics Letters B 120(1), 137 (1983). https://doi.org/10.1016/0370-2693(83)90639-1. http://www.sciencedirect.com/science/article/pii/0370269383906391
F. Wilczek, Phys. Rev. Lett. 40(5), 279 (1978)
J. Preskill, M.B. Wise, F. Wilczek, Phys. Lett. B 120(1), 127 (1983). https://doi.org/10.1016/0370-2693(83)90637-8. http://www.sciencedirect.com/science/article/pii/0370269383906378
G.G. Raffelt, Stars as laboratories for fundamental physics: the astrophysics of neutrinos, axions, and other weakly interacting particles (University of Chicago Press, Chicago, 1996) (1996)
M.S. Turner, Phys. Rep. 197(2), 67 (1990)
L.F. Abbott, P. Sikivie, Phys. Lett. B 120(1–3), 133 (1983)
P. Fox, A. Pierce, S. Thomas, arXiv:hep-th/0409059 (2004)
D. Spergel, Astrophys. J. Suppl. 170, 377 (2007)
F. Burger, E.M. Ilgenfritz, M.P. Lombardo, M. Müller-Preussker, A. Trunin, Nucl. Phys. A 967, 880 (2017). https://doi.org/10.1016/j.nuclphysa.2017.07.006. http://www.sciencedirect.com/science/article/pii/S037594741730324X (the 26th international conference on ultra-relativistic nucleus-nucleus collisions: quark matter 2017)
E. Berkowitz, M.I. Buchoff, E. Rinaldi, Phys. Rev. D 92(3), 034507 (2015)
S. Borsányi, Z. Fodor, J. Guenther et al., Nature 539(7627), 69 (2016)
G.G. di Cortona, E. Hardy, J.P. Vega, G. Villadoro, J. High Energy Phys. 2016(1), 34 (2016). https://doi.org/10.1007/JHEP01(2016)034
P. Petreczky, H.P. Schadler, S. Sharma, Phys. Lett. B 762, 498 (2016)
C. Bonati, M. DElia, M. Mariti et al., EPJ Web Conf. 137, 08004 (2017). https://doi.org/10.1051/epjconf/201713708004
P. Sikivie, Phys. Rev. Lett. 51, 1415 (1983). https://doi.org/10.1103/PhysRevLett.51.1415
J. Jaeckel, A. Ringwald, Annu. Rev. Nucl. Particle Sci. 60(1), 405 (2010). https://doi.org/10.1146/annurev.nucl.012809.104433
P.W. Graham, I.G. Irastorza, S.K. Lamoreaux, A. Lindner, K.A. van Bibber, Annu. Rev. Nucl. Particle Sci. 65(1), 485 (2015). https://doi.org/10.1146/annurev-nucl-102014-022120
I.G. Irastorza, J. Redondo (2018). arXiv:1801.08127
J.E. Kim, AIP Conf. Proc. 1200(1), 83 (2010). https://doi.org/10.1063/1.3327743. https://aip.scitation.org/doi/abs/10.1063/1.3327743
M. Arik, S. Aune, K. Barth et al., Phys. Rev. Lett. 112, 091302 (2014). https://doi.org/10.1103/PhysRevLett.112.091302.
S. Andriamonje, S. Aune, D. Autiero et al., J. Cosmol. Astroparticle Phys. 2007(04), 010 (2007). http://stacks.iop.org/1475-7516/2007/i=04/a=010
B.M. Brubaker, L. Zhong, Y.V. Gurev et al., Phys. Rev. Lett. 118, 061302 (2017). https://doi.org/10.1103/PhysRevLett.118.061302
S.J. Asztalos, G. Carosi, C. Hagmann et al., Phys. Rev. Lett. 104, 041301 (2010). https://doi.org/10.1103/PhysRevLett.104.041301.
L.D. Duffy, P. Sikivie, D.B. Tanner et al., Phys. Rev. D 74, 012006 (2006). https://doi.org/10.1103/PhysRevD.74.012006
B.T. McAllister, G. Flower, E.N. Ivanov et al., Phys. Dark Univ. 18, 67 (2017). https://doi.org/10.1016/j.dark.2017.09.010. http://www.sciencedirect.com/science/article/pii/S2212686417300602
Eleni Petrakou, for CAPP/IBS, EPJ Web Conf. 164, 01012 (2017). https://doi.org/10.1051/epjconf/201716401012. https://doi.org/10.1051/epjconf/201716401012
A. Caldwell, G. Dvali, B. Majorovits, Phys. Rev. Lett. 118(9), 091801 (2017)
N. Du, N. Force, R. Khatiwada et al., Phys. Rev. Lett. 120, 151301 (2018). https://doi.org/10.1103/PhysRevLett.120.151301
G. Ballesteros, J. Redondo, A. Ringwald, C. Tamarit, Phys. Rev. Lett. 118, 071802 (2017). https://doi.org/10.1103/PhysRevLett.118.071802
G. Ballesteros, J. Redondo, A. Ringwald, C. Tamarit, J. Cosmol. Astroparticle Phys. 2017(08), 001 (2017). http://stacks.iop.org/1475-7516/2017/i=08/a=001
A. Ernst, A. Ringwald, C. Tamarit, J. High Energy Phys. 2018(2), 103 (2018). https://doi.org/10.1007/JHEP02(2018)103
Y. Ema, K. Hamaguchi, T. Moroi, K. Nakayama, J. High Energy Phys. 2017(1), 96 (2017). https://doi.org/10.1007/JHEP01(2017)096
L. Calibbi, F. Goertz, D. Redigolo, R. Ziegler, J. Zupan, Phys. Rev. D 95, 095009 (2017). https://doi.org/10.1103/PhysRevD.95.095009
A. Garcon, D. Aybas, J.W. Blanchard et al., Quantum Sci. Technol. 3(1), 014008 (2018). http://stacks.iop.org/2058-9565/3/i=1/a=014008
D. Budker, P.W. Graham, M. Ledbetter, S. Rajendran, A.O. Sushkov, Phys. Rev. X 4, 021030 (2014). https://doi.org/10.1103/PhysRevX.4.021030
G. Ruoso, A. Lombardi, A. Ortolan et al., J. Phys. Conf. Ser. 718(4), 042051 (2016). http://stacks.iop.org/1742-6596/718/i=4/a=042051
R. Barbieri, C. Braggio, G. Carugno et al., Phys. Dark Univ. 15, 135 (2017). https://doi.org/10.1016/j.dark.2017.01.003. http://www.sciencedirect.com/science/article/pii/S2212686417300031
P. Vorobyov, A. Kirpotin, M. Rovkin, A. Boldyrev, arXiv:hep-ph/9506371 (1995)
R. Barbieri, M. Cerdonio, G. Fiorentini, S. Vitale, Phys. Lett. B 226(3), 357 (1989). https://doi.org/10.1016/0370-2693(89)91209-4. http://www.sciencedirect.com/science/article/pii/0370269389912094
A. Kakhidze, I. Kolokolov, Zh. Eksp. Teor. Fiz 99, 1077 (1991)
F. Caspers, Y. Semertzidis, in Cosmic axions. Proceedings, Workshop, Upton, 13–14 April 1989 (1989), pp. 0173–0183
M.S. Turner, F. Wilczek, Phys. Rev. D 42, 1001 (1990). https://doi.org/10.1103/PhysRevD.42.1001
P.J. McMillan, J.J. Binney, Mon. Not. R. Astron. Soc. 402, 934 (2010). https://doi.org/10.1111/j.1365-2966.2009.15932.x
F.J. Kerr, D. Lynden-Bell, Mon. Not. R. Astron. Soc. 221(4), 1023 (1986). https://doi.org/10.1093/mnras/221.4.1023
M. Augustine, Prog. Nucl. Magnet. Resonan. Spectrosc. 40(2), 111 (2002). https://doi.org/10.1016/S0079-6565(01)00037-1. http://www.sciencedirect.com/science/article/pii/S0079656501000371
S. Bloom, J. Appl. Phys. 28(7), 800 (1957). https://doi.org/10.1063/1.1722859
N. Bloembergen, R.V. Pound, Phys. Rev. 95, 8 (1954). https://doi.org/10.1103/PhysRev.95.8
C. Kittel, Phys. Rev. 71, 270 (1947). https://doi.org/10.1103/PhysRev.71.270.2
M. Sparks, C. Kittel, Phys. Rev. Lett. 4, 232 (1960). https://doi.org/10.1103/PhysRevLett.4.232
Y. Tabuchi, S. Ishino, A. Noguchi et al., Science 349(6246), 405 (2015). https://doi.org/10.1126/science.aaa3693. http://science.sciencemag.org/content/349/6246/405
Y. Tabuchi, S. Ishino, T. Ishikawa et al., Phys. Rev. Lett. 113, 083603 (2014). https://doi.org/10.1103/PhysRevLett.113.083603
I. Stern, A.A. Chisholm, J. Hoskins, Rev. Sci. Instrum. 86(12), 123305 (2015). https://doi.org/10.1063/1.4938164
C. Braggio, G. Bressi, G. Carugno et al., Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 603(3), 451 (2009). https://doi.org/10.1016/j.nima.2009.02.021. http://www.sciencedirect.com/science/article/pii/S0168900209003921
N. Crescini et al. Understanding photon-magnon hybridization through an analytical model (in preparation)
P. Arias, D. Cadamuro, M. Goodsell et al., J. Cosmol. Astroparticle Phys. 2012(06), 013 (2012). http://stacks.iop.org/1475-7516/2012/i=06/a=013
B. Beltran, M. Kuster, G. Raffelt, Axions. Lecture Notes in Physics (Springer, Berlin, 2008). https://doi.org/10.1007/978-3-540-73518-2_3
T. Battich, A. Córsico, L. Gabriel Althaus, M. Miguel Miller Bertolami, Sci. Rep. 2016, 062 (2016)
A.H. Córsico, A.D. Romero, L.G. Althaus et al., JCAP 1607(07), 036 (2016). https://doi.org/10.1088/1475-7516/2016/07/036
D.S. Akerib, S. Alsum, C. Aquino, C. Aquino et al., Phys. Rev. Lett. 118, 261301 (2017). https://doi.org/10.1103/PhysRevLett.118.261301
D.S. Akerib, S. Alsum, C. Aquino et al., Phys. Lett. B 724(1), 46 (2013). https://doi.org/10.1016/j.physletb.2013.05.060. http://www.sciencedirect.com/science/article/pii/S0370269313004425
E. Aprile, F. Agostini, M. Alfonsi et al., Phys. Rev. D 90, 062009 (2014). https://doi.org/10.1103/PhysRevD.90.062009
E. Aprile, F. Agostini, M. Alfonsi et al., Phys. Rev. D 95, 029904 (2017). https://doi.org/10.1103/PhysRevD.95.029904
E. Armengaud, Q. Arnaud, C. Augier et al., J. Cosmol. Astroparticle Phys. 2013(11), 067 (2013). http://stacks.iop.org/1475-7516/2013/i=11/a=067
Y.S. Yoon, H.K. Park, H. Bhang, J. High Energy Phys. 2016(6), 11 (2016). https://doi.org/10.1007/JHEP06(2016)011
A.V. Derbin, I.S. Drachnev, A.S. Kayunov, V.N. Muratova, JETP Lett. 95(7), 339 (2012). https://doi.org/10.1134/S002136401207003X
M. Giannotti, I.G. Irastorza, J. Redondo, A. Ringwald, K. Saikawa, J. Cosmol. Astroparticle Phys. 2017(10), 010 (2017). http://stacks.iop.org/1475-7516/2017/i=10/a=010
M. Cicoli, M.D. Goodsell, A. Ringwald, J. High Energy Phys. 2012(10), 146 (2012). https://doi.org/10.1007/JHEP10(2012)146
S.K. Lamoreaux, K.A. van Bibber, K.W. Lehnert, G. Carosi, Phys. Rev. D 88, 035020 (2013). https://doi.org/10.1103/PhysRevD.88.035020
L.S. Kuzmin, A.S. Sobolev, C. Gatti, IEEE Trans. Appl. Supercond. 28(7), 1 (2018). https://doi.org/10.1109/TASC.2018.2850019
Acknowledgements
The authors wish to thank Fulvio Calaon, Mario Tessaro, Mario Zago, Massimo Rebeschini, Andrea Benato and Enrico Berto for the help with cryogenics and for the mechanical and electronic work on the experimental setup. We acknowledge the support of Giampaolo Galet and Lorenzo Castellani for the building of the high-precision power supply and Nicola Toniolo, Michele Gulmini and Stefano Marchini for the work on the DAQ system. We also acknowledge Paolo Falferi and Renato Mezzena for the useful experimental suggestions, and finally we thank Riccardo Barbieri and Andreas Ringwald for the stimulating theoretical discussion and advice.
Author information
Authors and Affiliations
Corresponding author
Additional information
The original version of this article was revised: Due to a typesetting mistake, the intensity scale in Fig. 1 was missing. The original article has been corrected.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
Crescini, N., Alesini, D., Braggio, C. et al. Operation of a ferromagnetic axion haloscope at \(m_a=58\,\upmu \mathrm {eV}\). Eur. Phys. J. C 78, 703 (2018). https://doi.org/10.1140/epjc/s10052-018-6163-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-018-6163-8