New Constraints on the Axion–Electron Coupling Constant for Solar Axions

The resonant excitation of the 83Kr first excited nuclear level (E = 9.4 keV) by solar axions whose fluxes depend on the axion–electron coupling constant gAe is sought. The γ- and X-ray photons and the conversion and Auger electrons from the excited-level relaxation are detected with a gas proportional counter of a low-background detector in the underground Baksan Neutrino Observatory (Institute for Nuclear Research, Russian Academy of Sciences). As a result, a new constraint \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{|}}{{g}_{{Ae}}}(g_{{AN}}^{3} - g_{{AN}}^{0}){\text{|}}$$\end{document} ≤ 1.50 × 10–17 (90% C.L.) has been obtained for the axion–electron and axion–nucleon coupling constants, which corresponds to new constraints on the axion mass mA ≤ 320 eV and mA ≤ 4.6 eV in the KSVZ and DFSZ axion models, respectively.


INTRODUCTION
Light pseudoscalar particles, axions, were introduced in the theory to solve the CP problem of strong interactions [1][2][3]. Despite negative results of intensive experimental searches for axions, they are still well-justified candidates for the dark-matter constituents [4]. The anomalous transparency of the Universe to high-energy γ rays [5] and the overly fast cooling of some star systems compared to the predictions of theoretical models [6] can be treated as promising astrophysical indications of the existence of axions.
The interaction of axions with matter is specified by the quantity f A at which the Peccei-Quinn symmetry [1] is broken and is determined by the effective axion-photon (g Aγ ), axion-electron (g Ae ), and axionnucleon (g AN ) coupling constants. The axion mass m A and the parameter f A are related to the respective characteristics of the meson as . The corresponding expression for m A in terms of f A has the form [7,8] (1) The original PQWW model of the "standard" axion [1][2][3] implies that PQ symmetry is broken at the electroweak scale GeV. The constant f A in two classes of new models KSVZ [9,10] and DFSZ [11,12] of the "invisible" axion can be arbitrarily large up to the Planck mass m P ~ 10 19 GeV, thus reducing both the expected axion mass and the interaction of the axion with matter. Stars should be intense sources of axions. Intense fluxes of axions can be formed in the Sun in a number of processes whose probabilities depend on the axion coupling constants g Aγ , g Ae , and g AN . The constant g Aγ specifies the probability of conversion of photons to axions in the electromagnetic field of the solar plasma (Primakoff axions). The constant g AN determines the emission of axions in nuclear magnetic transitions that are thermally excited at high temperatures in the center of the Sun or appear in nuclear reactions of the рр chain and CNO cycle. The constant g Ae specifies the axion fluxes from bremsstrahlung and Compton process , as well as in discharge and recombination processes in atoms , . The spectra and fluxes of axions appearing in the above processes were calculated in [13][14][15][16][17] and are shown in Fig. 1.
Here, we report new results for the axion-electron coupling constant g Ae determined from the complete set of data obtained in experiments on the search for the resonant absorption of solar axions by 83 Kr nuclei with a large krypton proportional counter [18,19]. The 10 GeV = 5.69 (5) eV.
A A m f cross section for the resonant excitation of a nuclear level depends on the constant g AN ; consequently, the rate of absorptions of axions in an experiment depends on the product of the constants g Ae and g AN . The dimensionless coupling constant g Ae in the DFSZ model is expressed in terms of the parameter f A , which determines the axion mass m A , and the free parameter cos 2 β: where is the mass of the electron and β is an arbitrary angle. For the maximum value cos 2 β = 1 and taking into account Eq. (1), we obtain g Ae = 2.99 × , where m A is in electronvolts. In the KSVZ model, the axion does not interact with the electron; the effective axion-electron coupling constant calculated for the one-loop correction is [20,21] (3) where is the fine-structure constant, = 0.56 and = 0.029 are the ratios of the masses of the u, d, and s quarks, Λ ≈ 1 GeV is the scale cutoff in QCD, and E/N is the modeldependent parameter about unity, which is equal to 8/3 and 0 in the DFSZ axion model and the original KSVZ axion model, respectively. The interaction of the hadron axion with the electron is suppressed by a factor of at least . According to Eqs. (2) and (3), the coupling constant g Ae is proportional to m A , with the coefficient of proportionality including unknown parameters cos 2 β and E/N for the DFSZ and KSVZ axions, respectively. The axion coupling constants g Aγ , g Ae , and g AN determine both the probabilities of production of axions in various processes and cross sections for the reactions used to detect them.
The best known experiments are aimed at the search for solar axions produced through thermalphoton conversion in the field of the solar plasma. Under the assumption of the axion-photon coupling,  researchers attempt to detect axions through the inverse axion-photon conversion in a laboratory magnetic field [22,23] or in a crystal field [24,25]. Photon count rates expected in these experiments are proportional to . In this work, to detect solar axions appearing in reactions involving the electron, the resonant absorption of axions by 83 Kr nuclei is employed. In this case, the expected count rate is proportional to the dimensionless product . In our previous works, using a krypton proportional counter, we attempted to detect monochromatic solar axions emitted in the relaxation of the 83 Kr first excited nuclear level, which is excited because of the high temperature of the Sun [18], and Primakoff axions, which resonantly excite the 83 Kr nucleus [19]. Theoretical and experimental studies of the axion problem are reviewed in [8].

RESONANT EXCITATION OF THE 83 Kr
FIRST EXCITED NUCLEAR LEVEL BY SOLAR AXIONS Figure 1 shows the energy spectrum of axions produced in the processes caused by the axion-electron coupling that was calculated with the axion-electron coupling constant g Ae = 10 -11 [17] and is used in further calculations. The average energy of axions is 1.6 keV and their flux almost vanishes at energies above 15 keV. The flux of axions at an energy of 9.4 keV, which corresponds to the 83 Kr first excited nuclear level, is 1.32 × 10 11 cm -2 s -1 keV -1 [17], which is a factor of 60 lower than the maximum flux reached at an energy of 0.7 keV.
The 9.405-keV 83 Kr first excited nuclear level has the spin and parity 7/2 + . The transition to the 9/2 + ground state is an М1 magnetic transition (a small fraction of the Е2 transition is ) and can be accompanied by the emission of a γ-ray photon and by the emission or absorption of a pseudoscalar particle, an axion. The electron conversion coefficient important for target-detector experiments, where conversion and Auger electrons are absorbed in the target, is e/γ = 17.1 for the transition under study [26].
The ratio of the probabilities of the axion and electromagnetic transitions ω A /ω γ was calculated in the long-wavelength approximation in [27,28] in the form (4) where and are the momenta of the photon and axion, respectively; δ is the ratio of probabilities of the Е2 and М1 transitions; and are the isoscalar and isovector nuclear magnetic moments, respectively; and β* and η are the parameters determined by particular nuclear matrix elements.
The parameters β* and η for the nucleus with an odd number of nucleons and an unpaired neutron are estimated in the single-particle approximation as and [29].
The axion-nucleon coupling constant g AN is the sum of the isoscalar and isovector parts. In the KSVZ axion model, the constants and can be expressed in terms of the mass of the axion as [19] (5) To calculate Eqs. (5), we used the particular axial vector baryon coupling constants and and the polarization structure function of the proton [19], as well as the commonly accepted ratios = 0.56 and = 0.029 of the masses of the u, d, and s quarks (rather than more modern values and [8]) to correctly compare to the previous results.
We note that the detection of the resonant absorption of axions on the М1 transition in the 83 Kr nucleus is complicated by a methodological problem of a negative value of the parameter in Eq. (4). This negative value, together with the existing wide interval of possible S and z values, leads to a large uncertainty of the ratio ω A /ω γ given by Eq. (4). The effect of uncertainties in the parameters S and z on the final result through a factor of in Eq. (4) was discussed in [18,19,30], where it was shown in particular that the value at the relation can be more than an order of magnitude smaller than the value calculated with S = 0.5 and z = 0.56.
In the DFSZ axion model, the axion-nucleon coupling constants and depend on an additional unknown parameter cos 2 β, but they are of the same order of magnitude [20,31]. At cos 2 β = 1, the parameter appearing in Eq. (4) for the axion emission probability ω A /ω γ is a factor of 2.05 larger than this parameter for the KSVZ axion.
The cross section for axion absorption at the energy E A is given by an expression similar to the cross section for the resonant absorption of γ-ray photons with the correction to the ratio ω A /ω γ . The total cross section for axion absorption can be obtained by integrating over the solar-axion spectrum [30]. As a result, the rate of absorption of solar axions R A by the 83 Kr nucleus is given by the expression where cm 2 is the maximum γ-ray absorption cross section and keV is the width of the 83 Kr first excited nuclear level.
The flux of these axions is proportional to , and the ratio ω A /ω γ depends on the parameter . The resulting expression for the rate of axion absorption R A by the 83 Kr nucleus in units of inverse seconds per atom in terms of the coupling constant has the model-independent form (7) Substituting Eqs. (5) for the constants and in terms of the axion mass m A obtained in the KSVZ model, the axion absorption rate can be expressed in terms of g Ae and m A , the latter being measured in electronvolts, The total number of detected axions depends on the number of 83 Kr nuclei in the target, measurement time, and detector efficiency. The probability of observing the 9.4-keV peak depends on the background level of the experimental setup.

EXPERIMENTAL SETUP
The experimental setup is described in detail in our previous works [18,19]. Here, we present only the basic characteristics. A detector with a gas proportional counter is located in the low-background underground Baksan Neutrino Observatory (Institute for Nuclear Research, Russian Academy of Sciences) at a depth of 4900 mwe, where the muon flux is (2.6 ± 0.09) m -2 day -1 , which is lower than that on the ground by a factor of 5 × 10 6 [32].
The cylindrical copper counter has a total volume of 10.8 L. A gold-plated tungsten wire running along the counter axis serves as an anode. To exclude the influence of edge effects on charge collection, the anode diameter is increased, which limits the working volume of the chamber to 8.8 L. The counter is filled with krypton enriched to 99.9% in the 83 Kr isotope at a pressure of 1.8 bar. The mass of the 83 Kr isotope in the working volume of the counter is 58 g. The passive shield of the counter consists of sequential copper, lead, and polyethylene layers. A 12.5-MHz digitizer is used to measure the amplitude of the pulse, the duration of its rising edge, and the secondary photoemission pulse. The procedure of analysis of the pulse shape is described in [33,34]. The detection efficiency for the γ-and X-ray photons, as well as Auger and conversion electrons, arising from the relaxation of the 83 Kr 9.4-MeV excited level as estimated through a Monte Carlo simulation with the Geant4 package is [19]. 4. RESULTS The measurements were carried out. The energy spectrum of the proportional gas counter signals detected over a live time of 776.6 day is shown in Fig. 2. The most intense peak observed in the spectrum is due to the X-ray K Cu lines ( keV, keV, and keV) from the copper frame of the counter.
The second peak observed at an energy of 13.5 keV consists of several peaks with close energies. The long-lived 81 Kr isotope ( yr) is formed from the stable 82 Kr and 80 Kr isotopes under the action of neutrons and decays through electron capture into the 81 Br ground state with a probability of 99.7%. The absorption of the characteristic X rays and Auger electrons from bromide in the sensitive volume of the detector is responsible for the 13.47-keV peak corresponding to the binding energy of the electron in the K shell of the Br atom. X-ray photons emitted by krypton (K α12 = 12.65 keV) and bromine (K α12 = 11.92 keV) beyond the sensitive volume of the gas chamber make an additional contribution to the broadened 13.5-keV peak.
The measured spectrum is approximated in the range of 4-20 keV by the sum of a continuous background and four Gaussian peaks. The function describing the continuous background is the sum of a constant component and an exponential function of the energy E: (9) where a, b, and c are varying parameters. Three Gaussians describe the known 8.04-keV K α12 Cu and 8.905-keV K β Cu peaks and the wide 13.5-keV peak. The fourth Gaussian describes the 9.405-keV axion peak whose position and width are set to the parameters of the E 1 peak K α12 Cu.
The approximation of the energy spectrum in the range of 4-20 keV corresponding to χ 2 = 156.3/147 and P = 0.28 is shown by the solid line in Fig. 2. The 9.4-keV "axion" peak is not manifested statistically.
To establish an upper bound on the number of counts in the 9.4-keV peak, we used the standard method to determine the profile and the probability function . The upper bound on the number of peak events thus determined is for 90% C.L.
The determined upper bound on the number of events in the 9.4-keV peak provides bounds on the axion-electron coupling constant g Ae , -, and the axion mass m A according to Eqs. (7) and (8). The expected number of detected axions is (10) α = . 1 8 048 K α = . According to Eqs. (7) and (10), under the condition , which is valid for the masses of the axion m A < 2 keV, we obtain the upper bound (11) for 90% C.L. Bound (11) is a model-independent constraint on the coupling constants of the axion or any pseudoscalar axion-like particle to the electron and nucleon.
Using Eqs. (8) and (11), one can establish the following upper bound on the product of the constant g Ae and mass m A for the KSVZ axion: The upper bound (12) for the DFSZ axion at cos 2 β = 1 is almost two times lower eV. The upper bound (12) on the allowed values makes it possible to compare the results with the results of other experiments on the search for solar axions, in particular, with the results of the search for the axio-electric effect in atoms [36,37] (Fig. 3).
The upper bound (12) excludes a new region of the coupling constant g Ae at relatively large masses of the axion m A and is a factor of almost 4 lower than the experimental result on the detection of resonant excitation of the 169 Tm first excited nuclear level [16,35]. Since the 8.4-keV М1 transition in the 169 Tm nucleus is primarily a proton transition ( ), the ratio ω A /ω γ does not include uncertainty. Figure 3 shows the band of possible g Ae and m A values in the KSVZ and DFSZ axion models. Using Eqs. (2) and (3)  keV are the most stringent laboratory constraints and are close to astrophysical constraints. It is interesting that the analysis of the luminosity of white dwarfs indicates a nonzero |g Ae | value in the interval of (0.7-2.2) × 10 -13 [8,38], although the corresponding upper bound is above the bound obtained for red giants in some globular clusters [39] (Fig. 3).

CONCLUSIONS
The resonant absorption of 9.4-keV solar axions by 83 Kr nuclei with the excitation of the first excited nuclear level has been sought. Gamma-and X-ray photons, as well as conversion and Auger electrons, have been detected with a large gas proportional counter that is filled with the 83 Kr isotope and is placed in a low-background detector in the underground Baksan Neutrino Observatory (Institute for Nuclear Research, Russian Academy of Sciences). As a result, a new upper bound 1.50 × 10 -17 (90% C.L.) has been obtained for the axionelectron and axion-nucleon coupling constants. This bound corresponds to constraints on the axionelectron coupling constant g Ae and the axion mass m A eV and 10 -10 eV and on the axion mass m A ≤ 320 eV and m A ≤ 4.6 eV in the KSVZ and DFSZ (cos 2 β = 1) axion models, respectively.