Probing the axion-nucleon coupling with the next generation of axion helioscopes.

A finite axion-nucleon coupling, nearly unavoidable for QCD axions, leads to the production of axions via the thermal excitation and subsequent de-excitation of 57 Fe isotopes in the sun. We revise the solar bound on this flux adopting the up to date emission rate, and investigate the sensitivity of the proposed International Axion Observatory IAXO and its intermediate stage BabyIAXO to detect these axions. We compare different realistic experimental options and discuss the model dependence of the signal. Already BabyIAXO has sensitivity far beyond previous solar axion searches via the nucleon coupling and IAXO can improve on this by more than an order of magnitude.

Introduction

Axions [1, 2] are the direct and testable prediction following from the Peccei–Quinn (PQ) mechanism [3, 4] proposed as a possible solution to the strong CP problem of the Standard Model (SM) of particle physics. This problem, still one of the most puzzling in modern particle physics, concerns the apparent absence of CP-violating effects in Quantum Chromodynamics (QCD). To emphasize their role in explaining the observed QCD behavior, such particles are also dubbed QCD axions. Other kinds of pseudoscalar particles, with properties very similar to the QCD axion but no relation to the strong CP problem, emerge naturally in several extensions of the SM, in particular in compactified string theories [5-12]. To distinguish them from the QCD axion, these are often called axion-like particles (ALPs). The considerations discussed in this paper apply, in large part, to the QCD axion as well as to ALPs, except when specific QCD axion models are discussed. In the interest of brevity, we will refer to both QCD axions and ALPs as “axions” throughout this paper.

The axion phenomenology and experimental landscape are discussed in several recent reviews [13-17]. In summary, axions are expected to couple to the SM fields with model dependent couplings,where F is the electromagnetic field tensor and f are the SM fermions (for the present work, it is sufficient to consider couplings to protons, neutrons and electrons, so we can assume ).

Particularly appealing is the axion coupling to photons, since it allows very promising experimental strategies for the axion search. However, the couplings to electrons and nucleons are also employed in experimental axion searches.

In this work, we examine the axion coupling to nucleons (for a recent calculation see [18]) and explore the potential to probe it with axion helioscopes [19] which, as suggested by the name, search for axions produced in the sun.

From a theoretical point of view the nucleon coupling is of particular interest because it receives, similar to the photon coupling, an unavoidable contribution from the defining coupling of QCD axions to gluons. While it is possible to conceive models where this is (partially) cancelled, thereby having QCD axions with a suppressed coupling to nucleons [20], a large nuclear coupling is expected in most models which solve the strong CP problem and in certain cases can be considerably enhanced (see, e.g., Ref. [21]). A sizable nucleon coupling is therefore an expected feature of a QCD axion. In this sense, measuring the axion–nucleon coupling would be a good indication of the QCD axion nature. Furthermore, such a detection by a helioscope would likely be accompanied by a signal from Primakoff (and possibly Compton/Bremsstrahlung) axions. All these spectra have well defined shapes and so, given enough statistics and energy resolution, the data analysis might help understand specific properties of the detected particle (see discussion in Ref. [22] for a similar analysis in the case of the axion–photon and axion–electron couplings).

The axion coupling to nucleons can be probed indirectly in astrophysics, through the effects on the cooling of neutron stars (NS) [23-28] and from the analysis of the observed neutrino signal from SN 1987A [29-35]. A direct detection is also possible, in principle, through experiments such as CASPEr-gradient [36] or ARIADNE [37].1

Axions coupled to nucleons could also be produced in the sun, for example through the decay of excited nuclear states. As pointed out already in Weinberg’s seminal paper [1], being pseudoscalars, axions can be emitted in magnetic nuclear transitions. The best known example, and the one adopted here, is the decay of the first excited state of Fe with the emission of a 14.4 keV axion. Other nuclear transitions turn out to generate a substantially smaller axion flux (see also Appendix 1).

The search for these Fe axions has a long history. The currently most powerful helioscope, the CERN Axion Solar Telescope (CAST) [39], as well as CUORE [40] and, more recently, XENON1T [41], have searched for axions produced in this transition and provided constraints on the axion–nucleon coupling.

In this paper, we assess the potential of the next generation of axion helioscopes, BabyIAXO [42, 43] and IAXO [44-46], to detect such axions. An important motivating factor is that with BabyIAXO the transition from concept to real experiment is imminent, as its construction is expected to start in 2022. We therefore believe it is of great importance to assess its potential for an axion detection in all possible channels. Furthermore, the BabyIAXO and IAXO potential to probe this coupling are expected to be considerably superior to that of CAST [39], allowing to probe large regions of additional axion parameter space. We provide a guide for the best setups required to maximize the efficiency to detect axions from Fe. We also use this opportunity to include new theoretical developments. In particular, the matrix elements for the relevant transition have recently been revised [47], and show a  30% increase in the axion emission rate. In addition, we use a more recent solar model than previous publications. Finally, we discuss the model dependence of the Fe signal and identify a class of QCD axion models which yield an enhanced signal compared to the standard KSVZ [48, 49] and DFSZ [50, 51] axion models.

The paper is organized as follows. In Sect. 2, we revisit the problem of the solar axion production in nuclear transitions, provide an updated expression for the flux from Fe and discuss its dependence on the specific axion model; in Sect. 3, we present a brief discussion of the current stellar bounds; Sect. 4 is dedicated to the sensitivity estimates for BabyIAXO and IAXO, assuming different experimental setups; finally, in Sect. 5 we provide our conclusions. A discussion of other nuclei that could contribute spectral lines to the solar axion flux can be found in Appendix 1.

Axions from nuclear transitions in the sun

Axion helioscopes [19] such as the planned International Axion Observatory (IAXO) [42-46] are searching for very light (sub-eV) axions produced in the sun. Axions can be produced in the solar plasma by various processes. Helioscopes search primarily for axions originating from the Primakoff effect (cf, e.g. [52-58]) or from axion–electron interactions [22, 59–61]. This makes them sensitive to both the axion–photon coupling and to the product of photon and electron couplings [22, 62].

Axions can also be emitted in solar nuclear processes. There are two noteworthy mechanisms for this: nuclear fusion and decays as well as thermal excitation and subsequent de-excitation of the nuclei of stable isotopes.

Axions produced in nuclear fusion and nuclear decay processes typically carry an energy of the order of  MeV. The axion flux from the  (5.5 MeV) reaction, which provides one of the most intense axion fluxes from nuclear reactions, has been experimentally searched for by Borexino [63]. More recently, the SNO data [64] has also been scrutinized for traces of these 5.5 MeV axions. Helioscopes can, in principle, also be equipped with a -ray detector to look for such axions. For example, the CAST helioscope installed a -ray calorimeter for some time to gain sensitivity to these high-energy axions [65]. However, in general all these analyses provided somewhat weak bounds on the axion–nucleon coupling, since the solar flux from nuclear reactions is not expected to be very large. Indeed, Ref. [65] estimated an axion flux of the order of for the reaction mentioned above. This is more than 10 orders of magnitude smaller than the flux we will find below, Eq. (2.16).

A perhaps more promising direction is to look for low lying nuclear excitations of stable isotopes with a significant abundance inside the sun that can be thermally excited. Two candidates have been proposed for this in the past, Fe [66, 67] and Kr [68]. The former has a first nuclear excitation energy of 14.4 keV and the latter of 9.4 keV. In the solar core, at temperatures  keV, these excited states have a small but non-vanishing occupation number that can be calculated from the Boltzmann distribution. The amount of axions produced is proportional to the occupation number, the isotope abundance and the inverse lifetime of the excited state. By combining a list of possible elements and their nuclear transitions [69] with the solar abundances in [70], it becomes clear that for IAXO it is the Fe transition that would produce the strongest signal (see Appendix 1 for more details).

Effective axion coupling in the Fe nuclear transition

To describe the axion interactions with nuclei it is convenient to rewrite the relevant terms in the Lagrangian (1.1) asHere, is the nucleon doublet, , are the iso-scalar and iso-vector couplings, respectively, and is the Pauli matrix.

The axion-to-photon branching ratio for the decay rates of the first excited state of Fe can then be expressed as [71, 72]where are the axion and photon momenta, and are the isoscalar and isovector nuclear magnetic moments (expressed in nuclear magnetons), is the E2/M1 mixing ratio for the Fe nuclear transition, while and are constants dependent on the nuclear structure [72]. The constant can be safely neglected to the level of precision required in this work. The other constants were recently reevaluated in Ref. [47]. Assuming , which applies to ultrarelativistic axions, and adopting the most recent values in Ref. [47] for the relevant nuclear constants, we findwhere the second line is expressed in terms of the more common axion couplings to neutrons,and to protons,From the form of Eq. (2.4), it is convenient to define the effective nucleon coupling as2which is the coupling combination that controls the axion emission rate in this transition. Notice that the updated branching ratio is 27% larger than the one found in Ref. [71] and used in the previous experimental analyses, including CAST [39], CUORE [40] and more recently XENON1T [41]. Furthermore, the relative importance of the coupling to protons has strengthened in this new analysis even though, as evident from Eq. (2.4), it is still considerably less relevant than the coupling to neutrons.

Axion model dependence of the Fe transition rate

Before proceeding with the experimental sensitivity, it is interesting to look at the axion emission rate for some benchmark axion models.

Expressing the axion–nucleon couplings in terms of the dimensionless axion-quark coefficients [73] (defined via the Lagrangian term , with the axion decay constant), Eq. (2.4) can be cast as3where is the axion mass. In the KSVZ model [48, 49] one has , yieldingwhile in the DFSZ model [50, 51], and corresponding towith defined in the perturbative domain [15].

Note that in KSVZ models the axion emission rate is accidentally suppressed as in these models the neutron coupling is very small. In DSFZ models it can get enhanced by up to a factor of  4 with respect to the KSVZ.

Equation (2.9) suggests that a strong enhancement of the axion emission rate can be achieved for “down-philic” axions, with . More generally this holds if a sizable cancellation with the model-independent factor normalized to 1 in Eq. (2.9) is avoided. This last possibility naturally happens in a class of non-universal DFSZ models with that have been analyzed in Refs. [20, 74, 75] (for a summary of axion couplings in those models see Table 3 in [76]). For instance, in the non-universal M1 model of Ref. [20] one has , , , , with (also, and ), leading towhich at small yields an enhancement of the axion emission rate with respect to the KSVZ model.4 Other non-universal DFSZ models, among those mentioned above, feature a similar enhancement of the coupling, but they have different values for the axion–photon coupling (that is important in detection for the IAXO setup). In particular, the non-universal model of Ref. [74] features the largest axion coupling to photons among the general class of non-universal DFSZ models with two Higgs doublets (see Table 5 in [76]).

Solar Fe axion flux and limits

To calculate the resulting flux from thermal excitation and subsequent de-excitation in detail, we follow the derivation by Moriyama [66], also used by the CAST collaboration in their search for this source [39], but adopt updated nuclear matrix elements derived in Ref. [47]. The axion emission rate per unit mass of solar matter is given bywhere is the Fe number density per solar mass, the occupation number of the first excited state, the lifetime of the excited state, the internal conversion coefficient and the branching ratio of axion to photon emission, given in Eq. (2.4). A detailed description of these parameters is given in Appendix 1.

The axion signal from this nuclear transition is expected to be very narrow. The natural line width is negligible compared to the Doppler broadening ofcorresponding to a full width at half maximum (FWHM) of eV. Consequently, the spectral axion flux at a distance of one astronomical unit from the sun is an integral of a Gaussian peak over the solar radius ,5Although, depending on the optics adopted, IAXO’s field of view may cover only part of the sun, we integrate over the entire solar radius since – as illustrated in Fig. 1 – more than 99% of the total flux originates from a circle around the solar centre with radius 0.15 .

Fig. 1

Radial dependence of the solar axion flux from Fe transitions and the Primakoff effect. We show the fraction of the flux inside the field of view which is a circle around the solar center with radius r

The solar model (in our case B16-AGSS09 [77]) fixes temperature T, density and iron abundance at every radius.

It is noteworthy that this axion source is particularly sensitive to the temperature. This is because of the thermal occupation number of the excited nuclear state. For instance, the latest high and low metallicity solar models (B16-GS98 and B16-AGSS09 [77]) only differ by about 1% in their respective core temperature but this alone results in a difference of 12% in the total flux.

We do not expect the narrow, Doppler-broadened peak to be resolved.6 Hence, only the total flux is of interest and the integral over can be performed to get the total solar axion flux from the Fe nuclear transition [39],This numerical result is larger than the one previously derived by the CAST collaboration [39] because we worked with the updated nuclear matrix elements. In addition, we integrated over a more recent solar model, namely B16-AGSS09 [77], whose core temperature and iron abundance is smaller compared to the values adopted previously. This explains why the overall flux is not equally enhanced as the updated nuclear matrix elements. The following calculations are all done with this updated axion flux and the cited bounds derived by the CAST collaboration in [39] were rescaled accordingly.

Confronting stellar bounds

Let us start by quickly updating the energy loss constraint. Using the revised axion rate, Eq (2.4), and the solar model B16-AGSS09 we find a total axion luminosity (energy loss rate),via the Fe transition. The recent study in Ref. [79] constrains an exotic energy loss to a maximum of 3% of the standard solar luminosity . This leads to an improved bound on the axion effective coupling with nucleons,Notice that the result in Eq. (3.2), known as the solar bound on , is about a factor of two more stringent than the previous constraint, [39]. Besides the enhanced emission rate [47] and the updated solar model [77], this is also due to Ref. [39] having excluded only .

The axion–nuclear coupling can also be constrained from other stellar observations. In particular, strong bounds were derived from X-ray observations of various NS [23-28]. These bounds are, however, subject to several uncertainties and do not always agree with each other, not even when referring to the same star [76]. In any case, all these analyses suggest a limit of on some combination of axion–nucleon couplings.7

A similar bound can be deduced from the analysis of the neutrino signal observed in coincidence with the SN 1987A event [29, 31–33, 35, 80–84]. Here, we refer specifically to the most recent analysis [33], which derived the boundThis is a very strong constraint and, as we will see, only advanced setups may allow the exploration of smaller couplings. It is worth noticing, however, that the coupling expected in the Fe transition differs from the coupling in Eq. (3.3) and that, for some specific models, the coupling relevant for the solar axion searches may be enhanced or suppressed with respect to what is constrained by the supernova (SN) argument (see Fig. 2). For now, we will not investigate this argument further but we stress that, due to the peculiar effective coupling appearing in the Fe transition, the comparison with the astrophysical bounds is model dependent.

Fig. 2

Effective axion coupling entering in the Fe transition as a function of the angle , which defines the model dependent couplings (Cf. Sec. 2.2). In the figure, the colored lines indicate the value of the coupling calculated assuming that the specific axion model saturates the SN bound, given in Eq. (3.3)

That said, Eq. (3.3) can be translated into a limit on by choosing the ratio between the proton and neutron couplings such that the left hand side of this equation is minimal, while keeping constant. This yields,The only model in Fig. 2 which saturates this bound is the M1 model.

Sensitivity estimates

In this section, we present a detailed discussion of the helioscopes potential to detect the Fe axion flux. A particular emphasis is given to BabyIAXO, which is expected to start construction soon at DESY in Hamburg. Estimates for the more advanced IAXO and IAXO+ configurations [45] will also be presented.

Before going into specific experimental aspects let us make a few general considerations that can provide some help with the experimental design. To detect the Fe line we want to maximize the signal to noise ratio in the single relevant energy bin which contains all of the signal events. There are two main contributions to the background for the measurement of the nucleon line. One is the usual background rate of the detector itself (e.g., cosmic rays, environmental gammas, intrinsic detector radioactivity) that usually grows linearly with the area of the detector and it is typically proportional to the spectral size of the signal bin. This means that, in the case of the expected narrow signals, it can be reduced by making use of good energy resolution detectors. Second, there is the physics background due to Primakoff production. This grows with the photon coupling. Importantly, as the Primakoff spectrum is continuous, it also grows linearly with worsening energy resolution. Combining these two effects leads us to the following figure of merit,Here, are the optics and detector efficiencies, is the energy resolution of the detector,8b is the (spectral) background rate per area, and a is the signal spot area on the detector. quantifies the Primakoff flux in the Fe signal bin and is implicitly defined in Eq. (4.6). From this we can identify the parameters to be optimized. While good energy resolution is critical this should not be offset by too large a background. Similarly, with focusing optics we can reduce the detector area and therefore the background contribution. However, again there is a balance because such X-ray optics may have an efficiency significantly smaller than 1. We will see all this more concretely when considering explicit setups below.

BabyIAXO

Conceived as an intermediate stage towards the full IAXO experiment, BabyIAXO [42, 43] is nevertheless expected to substantially advance the exploration of the axion parameter space. With its sensitivity, the experiment will be able to study QCD axion models, and to investigate stellar cooling hints and other well motivated sections of the parameter space [15, 45].

The BabyIAXO experiment is mainly designed to measure Primakoff axions in the energy range from 1 keV to 10 keV with the peak of the solar axion flux spectrum at 3 keV. BabyIAXO will consist of two magnetic bores of 10 m length and 70 cm diameter, each with an average magnetic field strength of about 2 T. Together with newly developed X-ray optics and detector systems providing higher energy resolution and lower background, BabyIAXO will be the first helioscope to exceed the sensitivity of the CAST experiment. The magnet bores of BabyIAXO are of similar diameter as those of IAXO (60 cm) and IAXO+ (80 cm), thus the experience from the optics and detector development for BabyIAXO can later be applied to IAXO [42]. The first detection system for BabyIAXO is chosen to be microbulk micromegas technology. This detector technology has proven background levels as low as   [85] and a high detection efficiency for Primakoff photons <10 keV. A variety of other detector types, like silicon drift detectors (SDD), metallic magnetic calorimeters (MMC) and transition edge sensors (TES), are also studied for BabyIAXO aiming to optimize the energy resolution for precision measurements of the axion spectrum  [42].

BabyIAXO configurations for Fe detection

For the measurement of Fe axions at 14.4 keV, the detection efficiency in the baseline BabyIAXO configuration is not optimal and an enhancement of the detection system is necessary. Indeed, GEANT4 simulations of the Micromegas detector to be used in the baseline BabyIAXO configuration show that the ionization probability of 14.4 keV photons at a conversion length of 3 cm at the given Argon gas mixture and pressure is only around 15%. Additionally, the current designs of the BabyIAXO telescopes are not optimized for energies above 10 keV [86]. The expected optics efficiency at 14.4 keV is a mere (cf. “baseline”/BabyIAXO configuration in Table 1).

Table 1

List of experimental parameters adopted for all helioscope configurations which are considered in Figs. 3 and 4. As usual B is the magnetic field of the helioscope, L its length and A the area. t is the time that the helioscope is pointed at the sun. For these parameters we use values based on [45]. As already mentioned below Eq. (4.1) are the efficiencies of the optics and detector, b is the spectral background rate per detector area, and is the relative spectral resolution of the detector. Setup BabyIAXO is the baseline BabyIAXO, BabyIAXO is a version without optics, BabyIAXO assume optics optimized for the 14.4 keV line, with BabyIAXO including also a good energy resolution. In addition we show parameters from more advanced setups of IAXO and IAXO+

Label	BabyIAXO				IAXO		IAXO+
	Baseline	No	Optimized	High energy	Low	High energy	Low	High energy
		optics	optics	resolution	background	resolution	background	resolution
	BabyIAXO_0	BabyIAXO_1	BabyIAXO_2	BabyIAXO_3	IAXO_b	IAXO_r	IAXOb+	IAXOr+
B [T]	2	2	2	2	2.5	2.5	3.5	3.5
L [m]	10	10	10	10	20	20	22	22
A [m^2]	0.77	0.38	0.38	0.38	2.3	2.3	3.9	3.9
t [year]	0.75	0.75	0.75	0.75	1.5	1.5	2.5	2.5
b [1keV^2s]	10^-7	10^-6	10^-7	10^-5	10^-8	10^-6	10^-9	10^-6
ϵ_d	0.15	0.9	0.5	0.99	0.99	0.99	0.99	0.99
ϵ_0	0.013	1	0.3	0.3	0.3	0.3	0.3	0.3
a [cm^2]	0.6	3800	0.3	0.3	1.2	1.2	1.2	1.2
r_ω=ΔE_d14.4keV	0.12	0.12	0.12	0.02	0.02	514400	0.02	514400

In the following we discuss different adaptations to the detector system of BabyIAXO for the measurement of 14.4 keV photons with enhanced efficiency, energy resolution and lower background.

Let us start with a relatively minimal modification. The sensitivity of the Micromegas detector for 14.4 keV photons can be enhanced by changing the gas mixture and adjusting the pressure. Gas mixtures with an inert gas of higher Z-value, like Xenon, show a higher conversion efficiency for 14.4 keV photons, compared to Argon. A suitable high ionization probability for 14.4 keV photons of >90% should be reachable with a 10 cm photon conversion length and a Xenon based gas mixture at atmospheric pressure. Currently, microbulk Micromegas [87] are developed to have a square active area of 25 cm  25 cm [88] and the whole bore opening can be covered with a few detector tiles. As BabyIAXO features two bores, one could be used for operation without X-ray optics and this detector could be operated in parallel on the second magnetic bore. Removing the optics and covering of the whole area improves the detection efficiency, as this avoids the losses from inefficient X-ray optics at this energy. However, the detector will have a higher background due to the larger conversion volume than the smaller micromegas. This “no optics” configuration is denoted as “no optics”/BabyIAXO in Table 1 and in the figures.

The second detection concept for 14.4 keV photons at BabyIAXO is based on silicon drift detectors (SDD). This detector consists of a thick negative doped layer, that is fully depleted by a negative bias voltage, and positive doped contacts and strips on both sides of the layer. The incoming X-ray radiation generates electrons in the depleted zone that are then drifted towards the anode at the end of the layer. An SDD detector of 300 m thickness may reach 50% sensitivity for 15 keV photons [89]. This detector type is also considered as an additional detector for the baseline measurements at BabyIAXO [42]. Efficient X-ray optics focusing the photons on the small detector structure are mandatory, as SDDs are only produced in small pixels with a detector size in the range of millimeters. Preliminary studies showed that the optimization of the X-ray optical system for high energy X-rays is possible, using multilayer coating techniques [90]. The Nuclear Spectroscopic Telescope Array (NuSTAR) [91] has already used X-ray optics sensitive to photons in the range of 5-80 keV. A dedicated instrument with non-imaging optics can further boost the throughput at 14.4 keV. A realistic figure for the optics efficiency is , which is the value adopted in our analysis. In our study here, we assume that the optimized optics will be used in just one of the two BabyIAXO bores. This assumption will be lifted in our analysis of IAXO and IAXO+. With this included we obtain the “optimized optics”/BabyIAXO configuration (cf. Table 1).

Finally, a third detection concept suitable for the detection of 14.4 keV photons at BabyIAXO is based on Cadmium–Zinc–Telluride (CZT) semiconductors. The detection principle is similar to silicon-based ionization detectors. With a higher Z-value compared to silicon, CZT provides a higher ionization probability for the photons and CZT of only 300 m thickness have an ionization probability for 14.4 keV of >99%. Indeed CZT detectors are already used in experiments focusing on the detection of hard X-rays like NuSTAR [91]. These detectors can reach an energy resolution down to 2% in the relevant energy range [92], which makes them optimal for discrimination of the Primakoff background. However, the good energy resolution comes with the price of a somewhat increased background. Furthermore, currently it is only possible to produce detectors with an active area of up to 2 cm  2 cm [93], which again makes an efficient X-ray optic necessary for the use of CZT detectors at BabyIAXO. Coupled with the same NuSTAR-like optics adopted in BabyIAXO, this is our “energy resolution”/BabyIAXO configuration. Notice that, just like in the case of BabyIAXO, the detector and optics discussed above will be implemented in only one of the two BabyIAXO bores (see Table 1).

IAXO and IAXO+

The full scale helioscope IAXO is expected to begin construction during the BabyIAXO data-taking period. It adopts realistic components which will allow a considerably better performance than BabyIAXO in the entire mass range. IAXO+ is a more aggressive setup which will allow to increase the sensitivity and a deep exploration of physically motivated axion parameters (see. e.g., Refs. [45, 76]).

Here, we consider two possible setups for IAXO/IAXO+, summarized in Table 1. The first, indicated as IAXO, uses the benchmark configuration parameters discussed in the most recent IAXO publication [45]. In this case, we assume an energy resolution of 2% in the relevant energy range, which is a very realistic figure for CZT detectors, as discussed in Sect. 4.1. In addition we assume, somewhat optimistically, that significant improvements in the backgrounds of these detectors can be achieved.

Alternatively, we consider a further configuration, indicated as IAXO, with optimized energy resolution. Though generally realistic, we nevertheless note that the energy resolution in this configuration is non-trivial and may require new technologies such as microcalorimeters operated at mK temperatures [78, 94]. Magnetic microcalorimeters (MMCs) are also studied in the IAXO collaboration for precision measurements of the axion spectrum in the energy range <10 keV. First measurements have shown an energy resolution of 6.1 eV (FWHM) at the 5.9 keV Fe-peak [78]. To achieve the sufficient energy resolution, advanced cryogenic setups operating in the milli-Kelvin regime will need to be implemented and a sufficiently low background still needs to be established. We therefore have allowed for a somewhat larger background rate.

For IAXO+ we anticipate also additional improvements of the detection system combined with the enhanced magnetic field, area etc. outlined in [45].

Results (massless axions)

For our sensitivity study in the coupling space spanned by axion–photon and effective axion–nucleon coupling , we are going to assume that the axion is very light (20 meV) or massless. We will briefly comment on the effects of the mass later, in Sect. 4.5. For the moment, let us nevertheless note that in case of a discovery of an axion with a non-vanishing mass one can employ a suitable amount of buffer gas so that further measurements such as the one of the nucleon coupling can essentially be done as in the case of vanishing mass.

In the massless case, the conversion probability of an axion to a photon is energy independent,where B is the magnetic field and L the length of the conversion volume. Other experimental parameters that we used in the calculation include the cross-section of the magnet bores A, the total observation time t, optics and detector efficiencies, and , the size of the focal spot a, the background level b and the relative energy resolution at 14.4 keV .

We have considered several different combinations of parameters which are listed in Table 1. The benchmark values for the BabyIAXO, IAXO and IAXO+ magnets were taken from the most recent IAXO publication [45].

Because the FWHM of the Doppler-broadened iron peak is only  eV, it can be assumed that the whole signal is always in one energy bin.9 We can therefore calculate the expected number of signal events directly from the total flux given in Eq. (2.16)As already discussed at the beginning of the section, to find the Fe detection sensitivity, we also have to take all possible background sources into account. First, the detector background, which is quantified by the background level b and which can be measured accurately at times when the magnet bores are not pointed at the sun. Second, the tail of the Primakoff spectrum, which may act as an additional background. The expected background events are therefore given byThe Primakoff flux as well as the efficiencies and are in general functions of . In case of a sufficiently small energy resolution , we can average over the energy and describe the Primakoff background using the constant , as done in Eq. (4.6). If the Primakoff background to the Fe-peak is detectable at 14.4 keV, it is clear that there will be a much stronger Primakoff signal at smaller energies and we will have measured very precisely. Therefore, we either know the expected contribution from Primakoff axions to the number of background events or it is negligible compared to the intrinsic detector background.

The p-value of the Poisson-distributed observed number of counts k in the signal bin is given bywhere . In order to find the sensitivity in parameter space, we have to calculate the expectation value of p assuming a Poisson distribution for k with mean for each possible value of the two couplings. Note that is the sum of the usual detector background and Primakoff background. We regard the experiment as sensitive to a set of axion couplings if the expected p value, , is smaller than 0.05.10 The resulting sensitivity curves are plotted in Fig. 3.

Fig. 3

Model independent prediction for the sensitivity to the axion couplings for light axions ( meV). The different regions refer to the setups presented in Table 1. The dark red region is the solar bound discussed in the text (cf. Eq. (3.2)). The dark blue region represents the latest CAST exclusion regions from searches for the Primakoff flux [95] and the Fe peak [39] , which we have rescaled to match the updated axion flux from Fe transitions. The dashed horizontal green line indicates the expected sensitivity to the pure Primakoff flux. The supernova limit, Eq. (3.4), is not shown

Fig. 4

Model independent prediction of the IAXO sensitivity to the Fe peak, assuming that axions are produced only through the axion coupling to nucleons. In contrast to Fig. 3, we directly show the sensitivity of the various setups to the coupling combination under the assumption that the Primakoff background is negligible. The oscillations at higher masses are due to the form factor in the conversion probability. For comparison we show the effect of decoherence for a Primakoff spectrum as a dashed black line. The dark blue region represents the rescaled CAST result [39]. The dark yellow region indicates the parameter space expected for the DFSZ model. In brighter shades of yellow, we show the flavor non-universal DFSZ models M1 and . The dashed yellow line, on the other hand, shows the expected coupling for a nucleophilic QCD axion model of the kind presented in Ref. [21], with . All models with would be already accessible to BabyIAXO. Note that the experimental sensitivity estimates here do not assume the use of a buffer gas, which would extend the sensitivity to higher masses

Since we have assumed massless axions for these sensitivity estimates, we do not show the parameter space of DFSZ and M1 models (cf. Sect. 2.1). These models require large masses – meV – at those couplings, and helioscopes quickly lose sensitivity above  20 meV (cf. Fig. 4). This problem could be eased by filling the helioscope bore with a buffer gas [96].

Our analysis shows the potential of BabyIAXO to study areas of the parameter space, well beyond the solar bound and the region probed by CAST. The different green shaded areas show the experimental potential for different configurations, summarized in Table 1. The most efficient setups for BabyIAXO are the ones with optimized optics (labeled BabyIAXO in Table 1). As evident from the table, these setups allows for an enormous reduction of the total background by limiting the focal spot area a by  4 orders of magnitude with respect to the no-optics solution. Adopting this setting, BabyIAXO would be able to extend its detection potential also to regions of the parameter space below the Primakoff sensitivity, for . It is, therefore, possible, at least in principle, that BabyIAXO could discover axions through the Fe channel, before the Primakoff flux can be detected. If, on the other hand, axions have couplings in the green shaded area above the Primakoff sensitivity line (dashed green), one might have the opportunity to extract both couplings and derive information about the underlying axion model.

The more optimistic IAXO and IAXO+ configurations can explore an even larger area of parameter space. A noteworthy feature of these exclusion curves is their behaviour at different values of . At values of the detector background dominates over the Primakoff background. Therefore the figure of merit in Eq. (4.1) becomes . With the parameters given in Table 1, the configurations with minimized background slightly outperform the ones with optimized energy resolution in this regime (cf. Fig. 3). However, at the Primakoff background starts to play a role and eventually dominates. In this regime the figure of merit is given by . The detector background becomes negligible and the configurations with optimized energy resolution are significantly more sensitive than the ones with minimized background. Therefore, the ideal detector for the Fe line crucially depends on the value of . If BabyIAXO detects Primakoff axions, a detector with good energy resolution may be required to supress the Primakoff background to the Fe line. If, on the other hand, BabyIAXO only puts a stronger bound on , the energy resolution becomes less important and the low background detectors may be advantageous.

Effects of a finite axion mass

On the production side, the axion mass only becomes relevant at scales of the solar temperature  keV. Therefore, we can safely regard the solar axion flux as independent of the mass. However, on the detection side, a finite axion mass can cause decoherence between the photon and axion wave functions inside the magnet bores and lead to a signal suppression. The full expression for the conversion probability of axions into photons in the helioscope reads [96, 97]where q is the transferred momentum given by and is the energy of the axion. The suppression becomes relevant for masses above  20 meV, the exact value depending on the length of the magnet.

The decoherence effect can be compensated at the cost of some absorption by feeding a buffer gas into the bores [96, 98]. This is why the sensitivity study to effectively massless (i.e. ) axions in the previous section serves as a good benchmark. Nevertheless, we also want to explicitly investigate the sensitivity to massive axions without a buffer gas. To do this, we assume that the background from Primakoff axions is negligible and that instead the detector background dominates. In this case, the background is independent of any axion properties and the signal depends on the product of the two couplings, , as well as the mass. The statistical analysis is equivalent to the one in the massless case. From the expected signal and background we calculate the expected p-value and draw the sensitivity curves where p is expected to be smaller than 0.05. The results for a selection of viable setups are plotted in Fig. 4. The regions shaded in yellow indicate the coupling relations for the DFSZ, M1 and models (cf. Sect. 2.1).

The sensitivity curves are very similar to typical helioscope exclusion plots in the coupling vs. mass plane with two noteworthy exceptions. Because the Fe line is highly energetic at 14.4 keV, the transferred momentum q is smaller than for axions of the same mass from other solar processes. As a result the decoherence effect only becomes relevant at slightly higher masses in comparison to – for instance – Primakoff axions. To illustrate this effect, we have plotted the expected sensitivity of the IAXO setup with a decoherence factor from a Primakoff spectrum as a dashed black line in Fig. 4. Furthermore, the oscillations of the form factor for large qL are clearly visible in the Fe exclusion lines while they are washed out in the case of the broadband Primakoff spectrum.

Discussion and conclusion

In this work, we have presented the first dedicated investigation of the BabyIAXO and IAXO potential to detect 14.4 keV axions from Fe transitions in the sun. The analysis is based on a recent reevaluation of the matrix elements for the nuclear transition [47] and an updated solar model [77]. We have carefully considered different realistic setups for the detectors and optics to evaluate which combination gives the best performance for 14.4 keV axions. The configuration parameters are listed in Table 1.

Our results, summarized in Figs. 3 and 4, show that already BabyIAXO will be able to study a large section of interesting axion parameter space, well beyond the region accessed by CAST, particularly in the configurations with an optimized X-ray optics. The potential will be greatly improved with IAXO and IAXO+.

The figures also show representative QCD axion models. The yellow band in Fig. 4 represents the parameter area spanned by generalized DFSZ models [20, 76] (QCD axion models with two Higgs doublets plus a singlet scalar field, allowing for flavor non-universality). There, we present the specific examples of the classical DFSZ axion model, and of the M1 and models discussed in Sect. 2.1. We are not showing the well known KSVZ axion model since an accidental cancellation reduces its effective coupling to nucleons relevant in the Fe transition (see Sect. 2.1).

Although BabyIAXO is expected to have enough sensitivity to explore large sections of the parameter space for these models [76], a sizable axion flux from Fe transitions requires large axion masses, where BabyIAXO loses sensitivity. This problem can be eased with the use of a buffer gas [96], a technique already tested in CAST [95]. The sensitivities shown in the figures do not account for this option in BabyIAXO nor in its scaled up versions. A dedicated study may show if a buffer gas may allow to probe the M1 or other DFSZ-like models through the Fe line in the near future.

Less minimal models for QCD axions may present larger couplings to nucleons and be better accessible through the Fe channel. For example, the nucleophilic QCD axion models presented in Ref. [21] (and shown in Fig. 4), have exponentially large couplings to nucleons11 and are efficiently produced in Fe transitions even at lower axion mass. As shown in Fig. 4, practically, the entire class of these models will be accessible already to BabyIAXO, even without the need for a buffer gas.

One should nevertheless keep in mind that most of the region shown in the figures is in tension with astrophysical considerations, in particular, SN1987A (cf. Eq. (3.4)). As these are affected by their own uncertainties (e.g. relying on a single supernova event) it would nevertheless be comforting to have independent confirmation in more controlled setups. In addition, the IAXO+ setup shown in Fig. 3 approaches a level of sensitivity comparable to Eq. (3.4). This shows a pathway to pushing beyond the astrophysical limits.

In conclusion, helioscopes of the next generation may offer a unique chance to probe an interesting range of the - parameter space. While this potential is expected to be greatly improved with its scaled up versions, IAXO and IAXO+, already BabyIAXO will have enough sensitivity to detect axions with couplings to nucleons over an order of magnitude below the solar bound (see, Fig. 3). Furthermore, our analysis shows that already BabyIAXO, especially if equipped with optimized optics, has the potential to detect through the Fe channel axions too weakly coupled to photons to give a sizable Primakoff flux (region below the dashed green line in Fig. 3). In a more likely scenario, a detection of axions through Fe will be accompanied by a (larger) signal from Primakoff axions, allowing to extract important information about its couplings to both photons and nucleons.

Ending on an optimistic outlook, we note that discovery of an axion and its nucleon-coupling induced lines could perhaps also shed light on properties of the sun.12 For example, the strong temperature dependence may make this a good way to measure the sun’s core temperature.

Table 2

Isotopes with a nuclear M1 transition and  keV. The element abundances are taken from Ref. [70]. All other values are tabled in the appendix of Ref. [69]. The values in the last row were calculated by evaluating Eqs. (A.1) and (A.2) with the solar core temperature keV

	^57Fe	^83Kr	^169Tm	^187Os	^201Hg
E^∗ [keV]	14.4	9.4	8.4	9.7	1.6
J_0	1/2	9/2	1/2	1/2	3/2
J_1	3/2	7/2	3/2	3/2	1/2
τ_0 [ns]	141	212	5.9	3.4	144
α	8.56	17.09	285	264	47,000
ϵ	10^-4.5	10^-8.75	10^-11.9	10^-10.6	10^-10.83
a [%]	2.14	11.55	100	1.6	13.2
N_a(r=0)	1	1.8 ×10^-3	1.3 ×10^-4	3.0 ×10^-5	1.9 ×10^-6
[relative to ^57Fe]