Constraining the axion-photon coupling using radio data of the Bullet cluster.

Axion is one of the most popular candidates of the cosmological dark matter. Recent studies considering the misalignment production of axions suggest some benchmark axion mass ranges near [Formula: see text] μeV. For such axion mass, the spontaneous decay of axions can give photons in radio band frequency [Formula: see text] GHz, which can be detected by radio telescopes. In this article, we show that using radio data of galaxy clusters would be excellent to constrain axion dark matter. Specifically, by using radio data of the Bullet cluster (1E 0657-55.8), we find that the upper limit of the axion-photon coupling constant can be constrained to [Formula: see text] GeV[Formula: see text] for [Formula: see text] μeV, which is tighter than the limit obtained by the CERN Axion Solar Telescope (CAST).

Introduction

Observational data reveal that some unknown dark matter particles exist in our universe. Some proposals have suggested a class of hypothetical particles called weakly interacting massive particles (WIMPs) which can account for the dark matter[1]. These particles are fermions which may interact with ordinary matter or self-annihilate to give high-energy particles like photons. However, recent direct-detection experiments[2] and large hadron collider experiments[3] show null result of these particles. Also, a large parameter space of these particles has been ruled out by indirect detections, such as gamma-ray observations[4-6], cosmic-ray observations[7,8] and radio observations[9-12].

On the other hand, some early studies suggested that a class of light scalar or pseudo-scalar particles might exist[13,14]. In particular, a boson candidate called axion can help solve the CP-violation problem of the strong interaction[13]. Axions are massive and very stable particles which can probably account for the dark matter in our universe[15]. Axions can couple with photons which lead to the conversion between axions and photons () via the Primakoff effect in the presence of an external electric or magnetic field[16]. They can also decay into photons (), though the decay rate is very low. Examining the coupling between axions and photons provides a possible way for detecting or constraining axions. Some experiments like the International Axion Observatory (IAXO)[17] and Any Light Particle Search (ALPS)[18] are going to search for the signals of axions.

Some benchmark cosmological axion models suggest that the axion mass is of the order μeV[15,19,20]. For such axion mass, the photons emitted via the spontaneous decay process would be in radio band ( GHz). Therefore, using radio data to search for axion is possible. In this article, we show that using radio data of galaxy clusters would be excellent to constrain axion dark matter. We find that the axion–photon coupling constant can be constrained down to GeV using radio data of the Bullet cluster (1E 0657-55.8). Our limits are tighter than the CERN Axion Solar Telescope (CAST) limit[21] by nearly an order of magnitude within μeV.

Cosmological decaying axion model

There are many proposals which suggest axions to account for all cosmological dark matter. The possible mass of axions can range from to 10 eV[15]. For example, the simplest standard thermal freeze-out arguments suggest that eV[15,22]. Some other models considering symmetry arguments propose μeV[15,23]. Specifically, some recent studies follow the misalignment mechanism and suggest μeV[19,20]. A more recent calculation of the misalignment production of axions gives[15,24]:where is the number of relativistic degrees of freedom during the realignment process, is the initial angle of misalignment, is the Hubble parameter defined in km s Mpc. Following the standard assumption of the mean alignment angle and assuming axions being all cold dark matter (), we get μeV[15]. This proposal with such a narrow range of has become one of the most important benchmark cosmological axion models for further investigation. However, the assumption of the mean alignment angle does not involve the uncertainties from the QCD parameters that connect the axion mass and decay constant. Also, the presence of topological defects might lead to a significantly higher estimate of the dark matter axion mass[25]. Therefore, the narrow range of suggested by[15] might not be completely justified. Nevertheless, the range of μeV being axion dark matter is still a popular range which is worth to have further investigation. In the followings, we will focus on this particular range μeV and constrain the relevant parameters.

The frequency of the emitted photons in spontaneous decay of axions () is given by . Therefore, for μeV, the emission frequency is GHz. On the other hand, the decay time can be predicted theoretically as[15]Current general tightest upper limit for the axion–photon coupling constant is obtained by the CAST solar axion experiment: GeV for eV[21]. Therefore, for μeV, would be larger than s, which suggests that axions are very stable compared with the age of our universe ( s). It also reveals that detecting decaying signal of axions might be very difficult.

Fortunately, the spontaneous decay of axions could be greatly enhanced by the stimulated emission mechanism[15,26,27]. If the background contains a large amount of photons same with the emission frequency, the decay rate would be increased by a large amount. Such enhancement is characterized by the photon occupation number . The total radio flux density of the spontaneous decay in a structure with total mass M is given by[15]where D is the luminosity distance to the structure, z is the redshift of the structure and is the frequency width of the decay. There are a few components which can contribute to the photon background to enhance the decay. For example, the cosmic microwave background (CMB) contains a large amount of photons at GHz. The photon occupation number for the CMB photons iswhere K. This gives . The other components such as synchrotron radio background or hot gas thermal Bremsstrahlung radiation background may be able to give a larger enhancement.

Note that the stimulated decay of axions considered above is a little bit different from the resonant decay of the axion field. As stated in[28], the resonant decay assumes that photons produced from the decays of axions contribute to the photon occupation number responsible for stimulating the decay. Therefore, it will give an exponential growth of photons and the axions would almost completely decay in a short time. For the stimulated decay considered here, the photon occupation number is assumed only to arise from the background photons (e.g. CMB or hot gas). The decaying photons would not further stimulate the decay to trigger exponential growth of photons. Moreover, as discussed in[29], the resonant decay of axions requires that the axion momentum spread is not too large and the gravitational potential well in which axions are bound is not too strong. These constraints might prevent axions from decaying too fast so that most of the axions would remain in the universe[30]. However, these conditions do not affect the stimulated decay due to the background photons.

Data analysis

After an extensive search of archival radio data, we find that the radio data of the Bullet cluster (1E 0657-55.8) reported in[31] are the best for constraining the benchmark cosmological axion model. It is because it is a large galaxy cluster (the reason for focusing galaxy clusters will be discussed below) and its radio data involve a continuum frequency range covering GHz with relatively small uncertainties. Some previous studies have also investigated the effect of decaying axion-like dark matter using observational data of the Bullet cluster[32]. The redshift and the luminosity distance of the Bullet cluster are and Mpc (assuming Hubble parameter and following the standard CDM model with and ) respectively[31]. The observing radio frequency is 1.1–3.1 GHz, including wide band and narrow band (290 MHz sub-bands) methods[31]. Due to the cosmological redshift, the frequency of the decaying photons would decrease by a factor of when the photons arrive the radio telescope, which becomes GHz for μeV. Therefore, the observing frequency can cover most of the possible frequency range of the decaying photons.

As mentioned above, the spontaneous decay would be enhanced by the background photons via stimulated emission. In the Bullet cluster, a large amount of hot gas particles exist which are emitting thermal Bremsstrahlung radiation. The amount of radiation can be characterized by the hot gas temperature T. For thermal Bremsstrahlung emission, the hot gas is in local thermodynamic equilibrium such that the emissivity is the product of the absorptivity and the radiation spectral intensity (i.e. ). In view of this, the hot gas particles and photons are interacting via absorption and emission processes. The radiation spectral intensity represents the background photon distribution in the hot gas and it follows the Planck spectrum: . Note that the Bremsstrahlung emission (the photons escaped from the hot gas to reach us) does not follow the Planck spectrum and it mainly consists of X-ray photons (represented by the emissivity ) because the optical depth is very small for X-ray photons. However, our focus is the photons inside the hot gas, not the photons escaped from the hot gas. The background photon distribution inside the hot gas volume still maintains the Planck spectrum via absorption and emission equilibrium (represented by the spectral intensity ). Therefore, the photon occupation number in the hot gas is alsoFor the Bullet cluster, we have keV[33,34] and . Therefore, the stimulated emission can be greatly enhanced in the thermal Bremsstrahlung hot gas. For the synchrotron radio background, the corresponding photon occupation number is just , which is not significant. A simple physical intuition may think that the peak energy of thermal Bremsstrahlung emission is of the order keV, which is much higher than the photon energy required for stimulated emission ( eV) so that the hot gas component may have less effect on stimulated emission. However, the background photon distribution of the hot gas (not the escaped X-ray photons) also follows the Planck spectral distribution, which gives a very large photon occupation number at low energy. This is because the absorption of low-energy photons ( eV) in the keV-temperature hot gas is still important. For the synchrotron radio background component, the power-law distribution suppresses the photon occupation number so that this component does not give a significant enhancement in stimulated emission.

The frequency width can be written in terms of the velocity dispersion of axion dark matter : [26]. The average velocity dispersion of axion dark matter can be found by the Virial theorem:where is the molecular weight and is the proton mass. Putting keV, we get km/s, which gives a very narrow frequency width GHz. Therefore, a very sharp radio line would be observed if the spontaneous decay is strong enough. Due to the narrow frequency width, the radio data observed or considered must be in wide band or continuous band. Otherwise, the possible radio signal of the decay would be ignored if the observing frequencies do not match the specific decaying photon frequency.

From the whole radio spectrum ( GHz) reported in[31], we can analyse the axion decay for μeV. This range of would be divided into 6 sub-bands because there are 6 sub-bands of frequencies in the radio observations. For the narrow mass range μeV suggested in[15], there are two continuous sub-bands (1.6–1.9 GHz and 1.9–2.2 GHz) which can cover the observing frequency range of the axion decay (1.8–2.1 GHz).

Firstly, we can get the most conservative limit of the axion–photon coupling constant if we assume that all radio flux density detected in these sub-bands originates from axion decay. Putting , Mpc[31] and the virial mass of the Bullet cluster [34] into Eq. (3), we getNote that the value of the luminosity distance D is calculated based on the standard CDM model. The uncertainty of D would not be very significant because the uncertainty of the observed redshift z is small. However, the uncertainty in M and T might have some impact on the calculation of as is directly proportional to M and . Here, the percentage uncertainty in T is about 14%[33], which would only give about 7% uncertainty in . Nevertheless, the value of M is somewhat model-dependent, which might give a significant systematic uncertainty in . The value of M adopted in this study is calculated based on the parameters following the Navarro-Frenk-White (NFW) dark matter density profile[34,35], which is one of the most robust profiles to describe the dark matter density profiles in galaxy clusters[36].

The radio flux density upper limits for the frequency range GHz have been obtained from observations[31]. We can therefore obtain the conservative upper limit of for different using Eq. (7) (see Fig. 1). The conservative upper limit is GeV, which is tighter than the limit obtained by CAST ( GeV)[21].

Figure 1

The blue dashed line and blue solid line represent the conservative upper limit of and the two-component model upper limit of (with the mean value of in each bin) respectively. The cyan dot represents the overall best-fit across all bins. The region inside the cyan dotted contour is the range from the overall best-fit parameters. The region between the black dotted lines indicate the narrow mass range μeV suggested in[15]. The red dashed line, orange dotted line and the pink dotted line indicate the upper limits form the CAST experiment[21], projected upper limit of ALPS-IIc[18] and projected upper limit of IAXO[20] respectively. The region bounded by the green lines represent the HAYSTAC ruled out region[37].

Nevertheless, the synchrotron radiation of cosmic rays in a galaxy cluster usually dominates the radio emission. Thus, assuming that all radio flux density originates from axion decay is unrealistic. Moreover, considering the cosmic-ray contribution can further constrain the upper limit of . Using the radio flux density of other sub-bands can predict the background radio emission due to cosmic rays. Based on the radio spectrum of the Bullet cluster obtained in[31], the entire spectrum between 1.1 and 3.1 GHz can be best fitted by a power-law spectrum with a constant spectral index . This would become our null hypothesis for comparison (no axion decay scenario).

We now examine a two-component model: cosmic-ray contribution plus axion decay contribution. The total observed radio flux density should be a sum of two components . Here, we assume that is a sharp Gaussian function centred at the central frequency of each sub-band with frequency width . Here, we consider that the axion decay component is non-zero for GHz only. Moreover, the spectral index and the normalization constant are set as free parameters because adding the axion decay component would slightly alter these values obtained from the null hypothesis. Therefore, there are three free parameters (, and ) in the spectral fits. The goodness of fits can be examined by the value defined aswhere and are the observed radio flux density and their uncertainties respectively. As mentioned in[38], the Chernoff’s theorem states that the test statistic is asymptotically distributed according to when the null hypothesis is true[39,40]. Therefore, based on the test statistic value, we can determine the corresponding statistical significance.

After considering the two-component model, we can get the best-fit scenarios for different sub-bands (see Table 1). In particular, we find that there exist strong possible signatures of axion decay in the sub-bands μeV and μeV. The statistical significance is and respectively compared with the null hypothesis. The overall best-fit axion–photon coupling constant is GeV at μeV (with ). Furthermore, we can also release as a free parameter to get a contour of (, ) (i.e. with for 2 degrees of freedom). The region inside the contour indicate the range of (, ) from the best-fit parameters (see Fig. 1). Note that we have neglected the look-elsewhere effect here as we have constrained our range of in the analysis. We also show the spectral fits of the two-component model in Fig. 2 for the best-fit parameters. The peaks indicate the contribution of axion decay for the best-fit scenarios for different sub-bands. Besides, we can also determine the upper limits of ruled out at by comparing the two-component model with the null hypothesis (no axion decay). We assume the central value of for each mass bin (i.e. frequency bin) to get the upper limits of and corresponding parameters (see Table 2). By including the cosmic-ray emission, the upper limits of can be further constrained down to GeV (see Fig. 1).

Table 1

The fitting parameters of the null hypothesis ( model) and the best-fit two-component model ( model).

Model	m_a (μeV)	S_0 (mJy)	α	g_aγγ (GeV^-1)	χ^2
M_1		82	1.48		50.9
M_2	14–17	60	1.13	2.4×10^-11	16.4
	17–20	82	1.48	0	50.9
	20–23	82	1.48	0	50.9
	23–26	82	1.48	0	50.9
	26–29	85	1.54	6.3×10^-12	48.7
	29–32	90	1.65	8.9×10^-12	28.8

Here, the value of is the best-fit value for each sub-band.

Figure 2

The data points with error bars of the Bullet cluster radio spectrum are extracted from[31]. The coloured lines show the spectra of the two-component model with the best-fit parameters for different sub-bands of .

Table 2

The fitting parameters of the two-component model ( model) just ruled out at compared with the null hypothesis (with the mean value of in each bin).

m_a (μeV)	S_0 (mJy)	α	g_aγγ (GeV^-1)
14–17	44	0.79	3.4×10^-11
17–20	83	1.49	9.0×10^-12
20–23	82	1.48	3.8×10^-12
23–26	84	1.51	4.6×10^-12
26–29	88	1.62	9.6×10^-12
29–32	96	1.77	1.3×10^-11

Here, the values of are the upper limits.

Discussion

In this article, we have used the radio spectral data of the Bullet cluster to investigate the potential axion decay signal and constrain the axion–photon coupling constant. We have found the overall best-fit scenario (the largest statistical significance compared with the null hypothesis) which is for the sub-band μeV. The radio excess in this sub-band may indicate a potential signal of axion decay. Besides, we further constrain the axion–photon coupling constant down to GeV for the popular axion mass range μeV. The upper limits obtained are tighter than the CAST limit, especially for the mass range μeV suggested in[15].

In fact, some experiments are going to search for the signal of axions. For example, the IAXO and ALPS are performing experiments to search for axions and constrain the parameters of axions. The projected limits of constrained by these experiments are GeV (IAXO)[20] and GeV (ALPS-IIc)[18] respectively. Our upper limits for μeV are somewhat tighter than the ALPS-IIc projected limit and close to the IAXO limit (see Fig. 1). Therefore, using radio data of galaxy clusters is excellent to constrain axion decay as the limits can be comparable to the current axion searching experiments. It can be used as a complementary measurement for the axion search.

Some specific models suggest GeV for μeV[41]. Nevertheless, our upper limits obtained are still far from this predicted order of magnitude. Although some detections such as the ADMX haloscope[42,43] and HAYSTAC microwave cavity axion experiment[37] can constrain it down to GeV. However, the axion mass ranges constrained for these detections are μeV (ADMX) and μeV (HAYSTAC) respectively, which only cover a certain narrow mass ranges. Future plan of HAYSTAC experiment might be able to cover the entire range ( μeV)[44]. Besides, based on this current study, future radio observations of galaxy clusters with a very good sensitivity at frequency GHz (e.g. using the Square Kilometer Array) could also be helpful in detecting axion decay signals or constraining axion–photon coupling constant down to GeV.