Axion dark matter detection by laser induced fluorescence in rare-earth doped materials.

We present a detection scheme to search for QCD axion dark matter, that is based on a direct interaction between axions and electrons explicitly predicted by DFSZ axion models. The local axion dark matter field shall drive transitions between Zeeman-split atomic levels separated by the axion rest mass energy m a c 2. Axion-related excitations are then detected with an upconversion scheme involving a pump laser that converts the absorbed axion energy (~hundreds of μeV) to visible or infrared photons, where single photon detection is an established technique. The proposed scheme involves rare-earth ions doped into solid-state crystalline materials, and the optical transitions take place between energy levels of 4f N electron configuration. Beyond discussing theoretical aspects and requirements to achieve a cosmologically relevant sensitivity, especially in terms of spectroscopic material properties, we experimentally investigate backgrounds due to the pump laser at temperatures in the range 1.9 - 4.2 K. Our results rule out excitation of the upper Zeeman component of the ground state by laser-related heating effects, and are of some help in optimizing activated material parameters to suppress the multiphonon-assisted Stokes fluorescence.

Introduction

The nature of dark matter (DM) is the most long standing question in Big Bang cosmology, and direct searches may shed light on this intriguing mystery. The non-detection of particle DM in the heavy mass range (10 GeV to 10 TeV)[1-3] has motivated the scientific community to focus theoretical and experimental efforts on much lower mass particles[4,5]. Among them, a well motivated light particle is the QCD axion[6,7], introduced by Peccei-Quinn to solve the strong CP problem[8]. Axion physical properties are described by several models that can be grouped into the KSVZ[9,10] and DFSZ classes[11-13], depending on zero or full axion coupling strength to leptons, respectively. Even so, an almost model-independent statement holds for the axion masswhere f is the Peccei-Queen symmetry-breaking energy scale, inversely proportional to the coupling strenghts with standard model particles[6,7]. A light and stable axion emerges as an ideal DM candidate if large f are considered. Due to the resulting huge occupation number, galactic halo axions can be described as a classical oscillating field a, with oscillation frequency [14]. The eV axion mass range has since long been favoured by astrophysical and cosmological bounds[15], while very recent high-temperature lattice QCD calculations suggest that eV[16].

The axion is intensively searched in haloscope experiments[17], mostly based on resonant axion-photon conversion in a static magnetic field via Primakoff effect. The Axion Dark Matter eXperiment (ADMX) is the most sensitive haloscope detector based on high quality factor microwave resonators at cryogenic temperature. ADMX searches have excluded the mass range 1.9 < m  < 3.69 μeV[18,19]. The experiment HAYSTAC (formally ADMX-High Frequency)[20], designed specifically to search for axions in the 20–100 μeV range (5–25 GHz), has recently reached cosmologically relevant sensitivity at 24 μeV (~5.8 GHz).

The axion-electron coupling, explicitly predicted by DFSZ models[11-13], can be considered to envisage another class of haloscopes, thereby providing the opportunity to discriminate among axion models in case of detection. Complementary approaches may prove crucial to determine the fractional amount of axions as DM constituent. For instance, inhomogeneous filled cavities, in which the effective axion field is converted to magnetization oscillations of a ferrimagnet, are under study[21]. In this case, single photon detection is required, and it can be realized by e.g. superconducting circuit devices acting as quantum bits properly coupled to the cavity photons[22,23], but as yet their dark count rate still exceeds the axion interaction rate.

Approaches described so far are affected by an extremely poor sensitivity for axion masses above 0.2 meV (~50 GHz), where the effective detector volume is a critical issue. Extension to the mass range up to 1 meV (250 GHz) may be rather accomplished in suitable condensed matter experiments, in which the space parameters hardly accessible to cavity technology could be tackled with the upconversion scheme investigated in this work, whereby cosmological axions cause transitions between Zeeman split levels of suitable atomic species.

As target atoms we consider rare-earth (RE) elements inserted as dopants in crystalline matrices, where they exist as trivalent ions, substitutional for one of the atoms of the host with the same valence state and similar ionic radius. Among RE ions, those with an odd number of 4f electrons are called Kramers ions[24], and have electronic doublet levels with magnetic moments of the order of 1–10 Bohr magnetons μ . Therefore, using Kramers doublets, axion-induced spin transitions can take place in the GHz range with application of moderate magnetic fields. For instance, in Er3+, the calculated splitting spans from 20 to 120 GHz with applied magnetic fields in the interval 0.4 to 2.5 T[25], which translates to a large tunability in the favoured cosmological axion mass window.

In the direct axion-electron coupling[26,27] the interaction energy is , where the term plays the role of an effective oscillating field, is the electron magnetic moment with electric charge e and is the coupling constant[17]. Resonant condition is met when the Zeeman splitting energy is m c 2. As schematized in Fig. 1(a), the axion excitation is upconverted by a pump laser to photons in visible or infrared ranges, where single photon detection with ultra-low dark count rate has been already demonstrated[28-30].

Figure 1

(a) Detection scheme: axion induced transitions take place between the Zeeman split ground state levels, then a laser pumps the excited atoms to a fluorescent level. (b) Laser-induced fluorescence experimental setup. During the tests the crystal is immersed LHe and superfluid He. Fluorescence is collected orthogonally to the laser propagation direction by means of a mirror (M) that couples light to a 10 cm-long quartz guide. Optical filters (F) are set in front of the InGaAs photodiode (PD) or photomultiplier tube (PMT) to remove stray light. (c) Portion of the energy level diagram of YLiF4:Er3+ and transitions that are relevant for the present work[58]. Downward arrows indicate fluorescence transitions when ground state absorption takes place at about 809 nm laser pump wavelength.

The proposed detection scheme is based on electronic transitions between states within a configuration of the trivalent RE, with positions of the discrete energy levels minimally perturbed by the crystal-field due to the screening action of the 5s and 5p orbitals[24]. It is immediately evident that a first requirement for the feasibility of such a scheme is related to the the linewidth of the transition driven by the laser, which must be narrower than the energy difference between the atomic levels and .

Detectability of axions in this scheme can be at first discussed by considering only the thermal excitation of the atomic level as fundamental noise limit. Backgrounds of different nature are left for experimental investigations in the second part of the work.

We consider one mole of target atoms in the ground state and, using Eq. 8 of ref.[31], we establish the transition rate to the level by axion absorption on resonance:where is the transition rate of a single target atom, N is the Avogadro number, is the axion energy, is the coupling strength to the target atom and is of the order of one[31], and is the mean square of the axion velocity. The value eV (80 GHz) is a midpoint of the Zeeman splitting frequency interval reported for Er3+ in ref.[25]. As in the considered galactic halo model axions are galaxy dominant component of dark matter, we take for its energy density the value 0.4 GeV/cm3 obtained from galaxy rotational curves.

The experiment coherence time is set by , where t is the measurement integration time (inverse of the resolution bandwidth), and is the axion gradient coherence time at the resonant frequency of the experiment that can be calculated from the axion coherence time [32]. The latter is related to the width of the axion kinetic energy distribution in the laboratory frame. If we assume a Maxwellian velocity distribution in the Galactic rest frame and we take as the local dark matter virial velocity, we get , and finally[21] where the merit factor qualifies the axion-microwave linewidth in haloscope experiments.

The lifetime of the Zeeman excited state τ is typically much longer than , and in the rare-earth doped materials considered in this work is strongly dependent on temperature, intensity of the static magnetic field, dopant concentration[33-36]. The magnetic field, beyond splitting degenerate levels and thus opening a channel for resonant axion detection, may also inhibit spin flips and thus increase the lifetime τ of the intermediate level. Lifetimes much longer than ms have been measured in several rare-earth activated optical materials at liquid helium temperature with magnetic fields comparable to those used in this work ( T) up to about 3 T[34]. Incidentally, for a given pump laser intensity, the efficiency of the mentioned upconversion process is greater for longer τ, thus allowing for mitigation of the laser power requirements when large detecting volumes are devised[37].

As one might expect, the experiment must be operated in a ultra-cryogenic environment to minimize thermal population of the Zeeman excited level. To establish the working temperature of the apparatus, we treat the pumped crystal as if it were a single photon detector with overall efficiency (including the efficiency of upconversion, the fluorescence collection efficiency and self absorption) and consider thermal excitation of the target ions. For a given temperature T of the doped crystal, the corresponding excitation rate is related to the lifetime of the Zeeman excited levelwith average number of excited ions in the energy level , and k the Boltzmann constant. It is worth noticing that the contribution of adjacent Stark sublevels (due to interaction with the crystalline field) is not considered when their energy is much higher than , as is the case analyzed in this work. To ensure a signal to noise ratio (SNR) of 3 with a statistically significant number of counts within one hour observation time[38], the thermal excitation rate settles to Hz. If a level lifetime of τ = 1 ms is taken, axions with mass greater than 80 GHz can be searched, provided the active crystal is cooled down to at least 57 mK.

Results

We demonstrate here that we are able to resolve the involved transitions in a rare-earth doped material, that is a prerequisite for the feasibility of the proposed detection scheme. We then focus our experimental investigations on possible backgrounds induced by the pump laser at cryogenic temperatures and sub-Tesla magnetic field. The samples utilized in these measurements have 1% and 0.01% Er3+ -dopant nominal concentration, which correspond to and ions/cm3, with three electrons for each ion available as axion targets. Crystals with dopant concentrations well below 0.1% have been the subject of much scientific investigation for photon-echo-based optical data storage and data processing, owing to their narrow transition linewidths and long optical coherence times (see ref.[34] and references therein). In this work we are interested in the behavior of higher concentration samples to maximize the axion interaction rate given by Eq. (2) for a given laser-pumped, active detector volume. Moreover, the 1% concentration samples allow for higher sensitivity to laser related backgrounds in the measurements described in this section.

Zeeman splitting of the 4I15/2,5/2

We consider the electronic ground state of the erbium Kramers ion, the . The interaction with the crystal-field splits the magnetic sublevels into eight () Kramers doublets, labeled by , namely the absolute value of the J projection on the crystal optical axis c, . Through application of a magnetic field parallel to the c-axis, each Kramers doublet splits into two magnetic components: and . In YLiF4:Er3+, the lowest Kramers doublet is the . To determine its splitting under application of a 370 mT magnetic field, we measure the wavelength of transitions coupling the Zeeman components of the ground and the excited level by laser induced fluorescence (LIF) measurements.

In Fig. 2 we report the results obtained with the 0.01% concentration sample, where different incident laser polarization orientations (orthogonal and parallel to the crystallographic axis c) allow for a better detectability of possible transitions between Zeeman-split ground and excited states. The fluorescence spectrum displays sharp, well-separated lines. Among the observed transitions, only ( nm) and ( nm) are unambiguously identified. With the displayed data it is then not possible to determine whether the energy levels differences   μeV, eV and eV,   μeV represent the ground state or the excited level splittings. We accomplish this task by laser excitation of the thermal population in the Zeeman-split first excited Stark level of the ground state. Independently of this limitation, the plots in Fig. 2 demonstrate that we are able to resolve the Zeeman splitting and therefore that it is possible to monitor the population of the upper Zeeman component of the ground state. Clearly, at K thermal excitation of the level still prevents us from assigning a detection sensitivity to the present apparatus, but before we get to cool the sample to hundreds of mK temperatures, a thorough investigation of the pump laser related noise is accomplished as described in the following sections.

Figure 2

The ground level () splitting with a magnetic field of 370 mT is measured via laser induced fluorescence. During the measurements the 0.01% concentration sample is immersed in liquid He. We report the registered fluorescence with no magnetic field (a) and with 370 mT field (b) for laser polarization vector orthogonal or parallel to the crystallographic axis c. The black line in these plots is only a guide to the eye. Transitions between the Zeeman split levels of the ground state and the excited are identified as shown in part (c).

The LIF measurements in Fig. 2 have been repeated with the 1% concentration sample. In this case the Zeeman transition is hardly resolved due to increased transition linewidths, ascribable to spin cross relaxation processes due to direct interactions among Er3+ ions[36,39]. However, such a limitation might be overcome in the high magnetic field and low temperature regime, required to achieve ultimate sensitivity in the proposed axion detection scheme. For instance, in a 0.1% concentration sample of YLiF4:Er3+, authors have investigated the four transitions connecting the Zeeman sublevels of the ground and lowest excited state and demonstrated that their linewidth can be as narrow as MHz[40]. The applied magnetic field was about 3 T and measurements were conducted below 4 K by Zeeman-switched optical-free-induction decay technique. These results, together with our findings, foster the development of a few liters detector with intermediate concentration active materials, matching the axion-induced transition rate in Eq. 2 to dark count rates in available single photon counters[28-30]. As a final additional remark, we note that an intermediate concentration sample would allow for increasing the axion-electron interactions of six orders of magnitude compared to a gaseous target prepared by buffer cooling techniques[41].

Laser-induced thermal noise

To assess heating effects in the active detector volume, we focus on the population of the first excited Stark (crystal-field) sublevel , that has a strong thermal coupling with the ground energy level. To enhance the sensitivity of our tests, we use the 1% concentration sample. The crystal-field splittings of Er3+ ions in YLiF4 have been calculated and measured by previous authors[42,43] and for the sublevel the separation from the ground state is cm−1 (2.11 meV). If E Sj are the Stark sublevels energies of the level, the occupation probability of the first excited Stark level is proportional to , where the sum is well approximated with the first two terms. As the fluorescence intensity F is proportional to the laser intensity I times the level occupation number , we model a possible heating effect with the term in expression , where α, β are empirical parameters determined from a fit to the data.

From the data shown in the inset of Fig. 3, we infer that the β parameter is compatible with zero within one standard deviation, which allows us to limit the temperature increase to less than 4.5 mK/[W/cm2]. We stress that such limit is obtained in an unfavorable upconversion scheme, where the de-excitation takes place also through non radiative channels as shown in Fig. 1(c). Therefore we can assign the temperature of the thermal bath to the entire crystal and calculate the ratio of the populations of the same Stark level at two different temperatures T 1 and T 2. Such ratio is then compared to the LIF peak areas.

Figure 3

Probing the first excited Stark sublevel at 1.93 and 2.16 K in the 1% concentration sample. Inset shows a plot of the fluorescence amplitude for several pump laser intensity values measured at T = 4.2 K, with wavelength set at λ = 810.108 nm in resonance with the same Stark sublevel. The linear fit rules out temperature increments greater than 4.5 mK/[W/cm2] at 4 K, and allows in first approximation to assign the bath temperature to the laser-pumped crystal volume. The ratio between the two resonant peak areas is compared to the ratio of the Boltzmann factors calculated at 1.93 and 2.16 K.

As shown in Fig. 3, with the pump laser tuned to the transition , we obtain peaks that differ only in their area for K and K. The Er(1%):YLF crystal is immersed in superfluid He, and these points are obtained under λ-point operation at which bubbling disturbances are eliminated. A small satellite line is present on the right side of the main peak at both temperatures, which hinders an accurate fitting of the data. Therefore we compare the areas of the main peaks at K and K by summing the amplitudes of the data recorded at four wavelengths around resonance. The ratio of is in agreement with the expected value, confirming the assumption made in the introduction to calculate the rate of excited atoms via thermal bath temperature.

In addition, we confirm experimentally (Fig. 4) that also the population of the Zeeman-split levels follows Boltzmann statistics. In this case, the pump laser is set to probe the populations of the Zeeman-split levels of the first excited Stark sublevel ( and ). As a larger splitting is expected for this level as compared to the ground state[43], by pumping the transitions to the Zeeman levels of the , the previously measured energy differences (Fig. 2) can also be precisely identified. In fact, from the wavelengths reported in the first column of Table 1, we obtain the splitting of the first excited Stark level eV and eV, where the indices are assigned as described in the inset of Fig. 4. The latter value is consistent with the average of eV,   μeV obtained from data in Fig. 2. Consequently, we take the average eV as the searched splitting of the ground state. To further confirm proper identification of the Zeeman split levels, we can use the ratios of reported g factors in the same material oriented with its c-axis parallel to the magnetic field[43]. We obtain , in agreement with . The ground state splitting value we measured at 0.37 mT is also in fair agreement with theoretical values reported in ref.[25].

Figure 4

Zeeman transitions from the first excited Stark level of the ground state to the lowest level of the excited . The applied magnetic field is 370 mT. Data sets corresponding to lines 1 and 4 at T = 4.2 K have been acquired with a different linear amplifier gain and lock-in amplifier sensitivity, and have to be divided by a factor 30 for direct comparison with lines 2 and 3. The physical origin of the measured background level ( mV) evidenced by the horizontal lines in the plot at T = 1.93 K is clarified with the results reported in Fig. 5.

Table 1

Lorentzian fit of the data reported in Fig. 4.

	λ c	T = 1.93 K		T = 4.2 K
		A	w	A	w
1	810.171	0.25±0.03	2.7±0.2	0.671±0.03	2.9±0.2
2	810.131	16±0.2	2.3±0.2	1.6±0.1	2.4±0.2
3	810.084	32.8±2	2.1±0.2	2.37±0.23	2.5±0.2
4	810.044	0.76±0.05	3.8±0.4	1.05±0.03	3.0±0.2

The parameters (center, area and width respectively) are expressed in nm, in (nm · mV) and pm, respectively. Errors on the peak areas are assigned by considering the error on the measured background at T = 1.93 K.

Figure 5

Observed fluorescence in the 10300–12360 cm−1 interval (shadowed band in the inset). The LIF amplitude measured for pump laser resonant with transition is shown for comparison (circled data). The out-of-resonance fluorescence is attributed to multiphonon-assisted anti-Stokes (AS) and Stokes (S) emission. While the AS component is evidently suppressed at K, the Stokes fluorescence is represented in the plot by the exponentially increasing data.

That is as far as our LIF measurement of the ground level splitting is concerned. As for the investigations of possible laser-induced deviations from Boltzmann statistics, we consider the peak areas of LIF measured for the levels and . The data displayed in Fig. 4. are fitted to a Lorentzian curve in the form and the results are reported in Table 1. The ratio and is in fair agreement with 0.69 expected ratio at K. At K we obtain and , which averages to a value compatible with the calculated one 0.37.

Multiphonon background

In the data reported in Fig. 4. (lower side) the apparatus sensitivity allows observation of an out-of-resonance fluorescence level. To understand the physical origin of this background, we measure its intensity for values of pump laser wavelength in a wide range (corresponding to 10300–12360 cm−1), as shown in Fig. 5. The LIF measured with the pump laser tuned to transition is also plotted for comparison.

A similar exponential behaviour has been previously reported in Er:YLF and has been explained in terms of multiphonon-assisted, side-band absorption[44]. The RE manifolds can be in fact excited even by a nonresonant pump photon , when the missing/excess energy is bridged by absorption/emission of phonons via Anti-Stokes and Stokes processes, respectively. The related absorbed intensity is theoretically given by:where and are the Stokes and Anti-Stokes coefficients, described in the model[45] through expressions:

In Eq. (6), is the crystal effective phonon energy, p is the number of photons needed to bridge the energy gap, is the average occupation number and the Huang-Rhys coefficient that represents the electron-phonon coupling strength. Typical values of are smaller than 200 cm−1 in bromides, greater than 400 cm−1 in oxides[46], and in YLF the value 400 cm−1 is reported.

The strong suppression of the LIF observed in our 4.2 K data as soon as the pumping wavelength differs from the pure electronic transition wavelength (Fig. 5) is ascribable to the expected suppression of the AS component with temperature (from eq. 7) and the exponential growth for increasing wavelengths is then mainly due to the Stokes process. Fitting of data with wavelength greater than 850 nm give an absorption coefficient cm−1, in agreement with the value reported in ref.[47]. Ground state absorption measurements allow to estimate a cm2 cross section of the pure electronic transition (circled data in Fig. 5) and thus to quantify the upconversion efficiency and the multiphonon side band relative amplitude. This type of background hinders the application of the present scheme to axion detection, unless a suitable combination of rare-earth dopant, pumping pathway and matrix is chosen. In particular, relevant suppression of the background should be accomplished in low phonon energy host matrices[46] or, as suggested by Eq. (5), by exploiting pumping schemes with larger . It is worth mentioning that an ultimate laser-induced background might also originate from impurity absorption, the same process that is currently limiting the efficiency of optical refrigeration[48-50].

Discussion

We have discussed a solid-state approach for direct detection of axion dark matter, and established the most important experimental parameters necessary to reach cosmologically relevant sensitivity in DFSZ models. The effect of the continuous, coherent axion field is searched in the excitation of the Zeeman upper component of the ground state of rare earth ions in crystalline matrices, at the energy scale , corresponding to transitions in the GHz range. The population of this excited level is probed by a pump laser tuned to the transition to a fluorescent level within the same atomic configuration, so as to convert the axion excitation into photons, detectable with state-of-the-art single-photon detectors. Assuming thermal excitation of the excited Zeeman level as fundamental noise limit, the active detector volume must be cooled down to ultracryogenic temperatures. The rate of thermal excitation of the atomic level is directly related to its lifetime τ, and the temperature of less than K, at which the final experiment must be performed, has been estimated for 80 GHz axion-induced transitions and τ = 1 ms by requiring SNR . As long as ms, a high upconversion efficiency is also ensured for tens of W/cm2 pumping intensity.

In the proposed scheme it is important to address a thorough experimental study of pump laser-related backgrounds. As a first step we have probed the population of atomic levels close to the ground state via LIF measurements in the temperature range 1.9–4.2 K. Our main finding is that the pump laser does not affect the thermal population of the Zeeman excited level at least up to a few W/cm2 intensity. In addition, we have shown that it is crucial to optimize the pumping pathway and crystal properties to minimize scattering of the pump photons on crystal phonons (Stokes process).

As for the detection scheme via laser induced fluorescence, at 4.2 K and with 370 mT magnetic field, we have demonstrated that the four transitions coupling the Zeeman levels of the ground and the excited can be resolved in the lowest concentration YLiF4:Er3+ sample (0.01%). This was not possible in the 1% sample. However, the spin population dynamics in Kramers ions strongly depends on applied magnetic field, temperature, dopant concentration and species, and we argue that a tradeoff between these parameters can be found for the proposed experiment feasibility. A few liters active volume ensures ~mHz axion and thermal transition rates, corresponding to statistically relevant counts of upconverted photons in a measurement time of a few hours. Fortuitously, the dark count rate of state-of-the-art single photon detectors holds below the transition rate in the detector active volume.

We are witnessing a blooming of table-top experiments pursuing new observables for axion DM direct detection[51-57]. In such a multifaceted, dynamic scenario, our complementary proposal aims to probe the uncovered few hundred μeV axion mass region by exploiting the axion-electron interaction predicted in the DFSZ models.

Methods

The mm3 volume, Er3+:YLiF4 crystals used in this work were grown with the Czochralski method. The starting raw material were powders of LiF, YF3 and ErF3 with 5 N purity (99.999%), provided by AC Materials (Tampa, Fl., USA). They have 1% and 0.01% Er3+-dopant concentration (atomic percent substitution for Y3+). The growth was carried out in high purity (5 N) Argon controlled atmosphere. The pulling rate was 0.5 mm/h, the rotation rate 5 rpm and the temperature of the melt was computer-controlled around 880 °C to maintain a constant diameter of the boule; the seed was an undoped YLF monocrystal oriented along the c axis. The crystals resulted of good quality, free of internal cracks, microbubbles or inclusions. The structural analysis by Laue X-ray diffraction checked the mono-crystalline feature, and allowed the orientation of the boules along crystallographic axes.

To allow for Zeeman studies at LHe and superfluid He temperatures, the samples were housed in an immersion dewar located between two NdFeB magnetic discs that produced a field of 370 mT at the sample position. The c-axis of the crystal was parallel to the magnetic field direction. As shown in Fig. 1(b), the sample fluorescence is collected orthogonally to the laser pump propagation direction and coupled to the photon detector by means of a mirror M and a quartz guide. With optical filters we suppress scattered pump radiation at signal wavelengths, and at the InGaAs photodiode we detect the 1.5 μm component of the overall infrared fluorescence spectrum (see Fig. 1(c)). The employed optical source is a cw Ti:sapphire laser, which can be finely tuned by rotating intracavity ethalons. Zeeman studies are conducted around 809 nm wavelength, while laser-induced backgrounds are investigated at 810.1 nm. The laser linewidth is GHz (∼1–2 pm), comparable with the detected transitions widths. The incident light polarization angle is varied by means of a half-wave plate. A typical value of laser intensity used in our measurements is 10 W/cm2, compatible with 0.1 upconversion efficiency in trivalent ions[37]. For the laser noise studies the pump laser was chopped at 15 Hz to allow phase-sensitive detection.