Search for axioelectric effect of 5.5 MeV solar axions using BGO detectors

A search for axioelectric absorption of solar axions produced in the $ p + d \rightarrow {^3\rm{He}}+\gamma (5.5 \rm{MeV})$ reactions has been performed with a BGO detector placed in a low-background setup. A model-independent limit on an axion-nucleon and axion-electron coupling constant has been obtained: $| g_{Ae}\times g_{AN}^3|<2.9\times 10^{-9}$ for 90% confidence level. The constrains of the axion-electron coupling have been obtained for hadronic axion with masses in (0.1 - 1) MeV range: $|g_{Ae}| \leq (1.4 - 9.7)\times 10^{-7}$.


INTRODUCTION
The axion is a pseudo-Goldstone boson arising from breaking the global Peccei-Quinn (PQ) symmetry [1][2][3]. The original so-called Weinberg-Wilczek-Peccei-Quinn axion model contained some certain direct predictions for the coupling constants between axions and photons (g Aγ ), electrons (g Ae ) and nucleons (g AN ) as well as for the axion mass (m A ). This model was quickly disproved by the reactor and accelerator experiments [4].
The particle has been retained by two new theoretical models of the "invisible" axion in the form required to solve the CP problem of strong interactions; at the same time its interaction with matter is suppressed. These are models of the "hadronic"(or Kim-Shifman-Vainstein-Zakharov (KSVZ)) [5,6] and the "GUT" (or Dine-Fischler-Srednicki-Zhitnitskii (DFSZ)) [7,8] axion. The scale of Peccei-Quinn symmetry violation is arbitrary in both models and can be extended up to the Planck mass ≈ 10 19 GeV. The axion mass in these models is determined by the axion decay constant f A : ; m A (eV) ≈ 6.0 × 10 6 where m π and f π are, respectively, the mass and decay constant of neutral π-meson and z = m u /m d is quark-mass ratio.
Since the axion coupling constants g Aγ , g Ae and g AN are proportional to the axion mass, the interaction of an axion with matter is suppressed.
The main experimental efforts are focused on searches for an axion with a mass in the range of 10 −6 to 10 −2 eV. As this range is free of astrophysical and cosmological constraints, relic axions with such a mass could be considered to be the most likely candidates for dark matter particles. Experimental bounds on the mass of the axion follow from the constraints on the g Aγ , g Ae , and g AN coupling constants, which significantly depend on the theoretical model used.
New possibilities for solving a strong CP problem are based on the concept of the existence of the mirror particles [9] and supersymmetry [10]. These models suppose the existence of axions with the mass of about 1 MeV, and this existence is forbidden by neither laboratory experiments nor astrophysical data.
This work is devoted to the search for solar axions with an energy of 5.5 MeV, produced in the p + d → 3 He + A reaction. The range of axion masses under study is expanded up to 5 MeV. The axion flux is proportional to the pp-neutrino flux, which is estimated with a high degree of accuracy [11]. The axion interaction exploited in this study is the axioelectric effect A + e + Z → e + Z. As the cross section of the reaction depends on the charge of a nucleus as Z 5 , the BGO detector containing bismuth (stable nucleus with the largest charge Z = 83) is the most suitable from this point of view.
Recently, the high energy solar axions and axions from a nuclear reactor have been sought by the Borexino [12,13], the CAST [14] and the Texono [15] collaborations. We have previously published a search for 5.5 MeV axions with two small BGO detectors [16].

AXION PRODUCTION IN NUCLEAR MAGNETIC TRANSITIONS AND THE AXIOELECTRIC EFFECT
Since the temperature in the center of the Sun is about 1.3 keV, the Sun should be an intense source of low energy ∼(1-10) keV axions. The reactions of the main solar chain and CNO cycle can produce axions with higher energies. The most intensive flux is expected as a result of the reaction: p + d → 3 He + γ when 5.5 MeV axion is emitted instead γ−quantum.
According to the standard solar model (SSM), 99.7% of all deuterium is produced as a result of the two protons fusion, p + p → d + e + + ν e , while the remaining 0.3% is due to the p + p + e − → d + ν e reaction. The expected solar axion flux can thus be expressed in terms of the pp-neutrino flux, which is 6.0 × 10 10 cm −2 s −1 [11]. The proportionality factor between the axion and neutrino fluxes is determined by the axion-nucleon coupling constant g AN , which consists of isoscalar g 0 AN and isovector g 3 AN components. In the p(d, 3 He)γ reaction, the M1-type transition corresponds to the capture of a proton with a zero orbital momentum. The probability of proton capture from the S state at proton energies below 80 keV has been measured in [17]; at the proton energy of ∼ 1 keV, M1 fraction of the total p(d, 3 He)γ cross section is χ = 0.55. The proton capture from the S state corresponds to the isovector transition, and the ratio of the probability of a nuclear transition with the axion production (ω A ) to the probability of a magnetic transition (ω γ ) depends only on g 3 AN [18][19][20][21][22]: where p γ and p A are, respectively, the photon and axion momenta; α ≈ 1/137 is the fine-structure constant; and µ 3 = µ p − µ n ≈ 4.71 is isovector nuclear magnetic momenta.
Within the hadronic axion model, the constant g 3 AN can be written in terms of the axion mass [23,24]: The value of g 3 AN in the DFSZ model depends on an additional unknown parameter cos 2 β , but have the same order of magnitude. The numerical value lies in the range of (0.3 − 1.5) from the the value of g 3 AN for the hadronic axion [23].
The calculated values of the ω A /ω γ ratio as a function of the axion mass are shown in Fig.1. The axion flux on the Earth's surface is where Φ ν pp = 6.0 × 10 10 cm −2 s −1 is the pp-neutrino flux. To detect 5.5 MeV axions, we chose the reaction of axioelectric effect A + Z + e → Z + e which is caused by the axion-electron interaction. The cross section of the axioelectric effect depends on the nuclear charge according to the law, and therefore, it is reasonable to search for this process using detectors with a large Z. The cross section for bismuth atoms exceeds the cross section of the Compton conversion of an axion (A + e − → γ + e − ) by almost two orders of magnitude. The detection efficiency for the produced electron is almost 100% and the background level at 5.5 MeV is much lower than in the range of natural radioactivity. As a result, the sensitivity to constants g Ae and g AN can be high even in an experiment using a relatively small target mass. In the axioelectric effect (an analog of the photoelectric effect), an axion disappears and an electron with the energy of E e = E A − E b , where E b is the electron binding energy, is emitted from the atom. The axioelectric effect cross section for K-shell electrons has been calculated (on the assumption that E A ≫ E b and Z ≪ 137) in [25]: The dependence of the cross section on the axion mass for the coupling constant g Ae = 1 is shown in Fig. 1. The K-shell electrons make the main contribution to the cross section. The contribution from the other electrons was incorporated by introducing a factor of 5/4, by analogy with the photoelectric effect.

INTERACTION OF AXIONS WITH SOLAR MATTER AND AXION DECAYS
The flux of 5.5 MeV axions on the Earth's surface is proportional to the pp-neutrino flux only when the axion lifetime exceeds the time of a flight from the Sun and when the flux is not reduced as a result of the axion absorption by matter inside the Sun. The requirement that most axions escape the Sun and reach the Earth thus limits the axion coupling strengths accessible to terrestrial experiments [12,13,16,19].
Axions leaving the center of the Sun pass through the matter layer of ≈ 6.8 × 10 35 electrons and ≈ 5 × 10 35 protons per cm −2 in order to reach the Sun's surface. The maximum cross section of the axio-electric effect on atoms for 5.5 MeV axions is σ Ae ≈ g 2 Ae Z 5 3.8 × 10 −33 cm 2 (see Fig.1 for Bi). Axion loss due to axioelectric absorption by Fe atoms (Z=26, N Fe /N H = 2.8 × 10 −5 ) imposes an upper limit on |g Ae | < 10 −4 above which the sensitivity of terrestrial experiments using solar axions is reduced. The abundance of heavy (Z > 50) elements in the Sun is N Z /N H ∼ 10 −9 in relation to hydrogen [26]. If |g Ae | < 10 −3 , the change in the axion flux does not exceed 10%.
The other process -Compton conversion of an axion into a photon -imposes stronger limit on the sensitivity of Earthbound experiments to the constant g Ae . The cross section of this reaction for 5.5-MeV axions depends weakly on the axion mass and can be written as σ cc ≈ g 2 Ae 4 × 10 −25 cm 2 [13,18]. For |g Ae | values below 10 −6 , the axion flux is not substantially suppressed.
The axion-photon interaction, as determined by the constant g Aγ , leads to the conversion of an axion into a photon in a field of nucleus (Primakoff conversion). The cross section of the reaction is σ pc (5.5MeV) ≈ g 2 Aγ Z 2 2 × 10 −29 cm 2 . Taking into account the density of protons and 4 He nuclei, the condition that axions efficiently escape the Sun imposes the constraint |g Aγ | < 10 −4 GeV −1 . Constraint for the other elements are negligible due to their low concentration in the Sun.
The axion-nucleon interaction leads to axion absorption in a threshold reaction similar to photo-dissociation: A + Z → Z 1 + Z 2 . For axions with energy 5.5-MeV this can occur for only a few nuclei: 17 O, 13 C, 3 He and 2 H [19]. It was shown in [19] that the absorption of axions can be mainly due to the process A + 17 O → 16 O + n. The cross section of the reaction depends on the combination of iso-scalar and iso-vector coupling constants. If | − g 0 AN + g 3 AN | < 10 −2 , the change in the axion flux does not exceed 10%. Because the absorption of axions due to isovector-transition A( 3 He, d)p is insignificant, the range of values |g 3 AN | available for the study is practically unlimited.
For the axions with a mass above 2m e , the main decay mode is the decay into an electron-positron pair. The condition that 90% of all axions reach the Earth limits the sensitivity of the solar axion experiments to |g Ae | < (10 −12 − 10 −11 ) [16]. On the other hand the measured value of the interplanetary positron flux allowed us to establish a new constraint on the axion-electron coupling constant: g Ae ≤ (1 − 5) × 10 −17 for axions with masses in the range of (1.2 − 5.4) MeV [16].
If the axion mass is less than 2m e , A → e + + e − decay is forbidden, but the axion can decay into two γ quanta. The probability of the decay depends on the axion-photon coupling constant and the axion mass -τ A→γγ = 64π/g 2 Aγ m 3 A . The present-day experimental constraint on g Aγ is 10 −9 GeV −1 , which corresponds to τ cm = 10 5 s for 1-MeV axions. This means that the axion flux is not practically reduced due to the A → 2γ decay even for axions with 5 MeV mass.

EXPERIMENTAL SETUP
We used a 2.46 kg BGO crystal, manufactured from bismuth orthogermanate Bi 4 Ge 3 O 12 (1.65 kg of Bi), to search for the 5.5 MeV axions. The BGO crystal was grown at the Nikolaev Institute of Inorganic Chemistry and it was shaped as a cylinder, 76 mm in diameter and 76 mm in height. The detector signal was measured by an R2887 photoelectron multiplier, which had an optical contact with a crystal end surface.
The spectrometric channel of the BGO scintillation detector included an amplifier with the shaping time of 1 µs and a 12-digit ADC. The amplification was selected so that the ADC channel scale was 2.7 keV. The standard calibration sources ( 60 Co and 207 Bi), in combination with the natural radioactivity lines of 40 K and the uranium and thorium families, were used for the energy calibration of the detector. The energy dependence of the detector resolution σ can be presented as σ /E ≈ 3.9% × E −1/2 , where E is in MeV. The detection efficiency ε for 5.5-MeV electrons in the BGO crystal was calculated with GEANT4. The response function consists of Gaussian peak and a flat tail, the number of events in Gaussian peak is ε = 0.67.
The external γ activity was suppressed using a passive shield that consisted of successive layers of lead (90 mm) and bismuth (15 mm Bi 2 O 3 ). The total thickness of the passive shield was ≈ 110 g cm −2 .
The setup was located on the Earth's surface. In order to suppress the cosmic-ray background we used an active veto, which consisted of five 50 × 50 × 12 cm plastic scintillators. The active veto included two energy thresholds: the first was set at a level ∼ 5 MeV, which corresponds to 600 s −1 counts rate and 4% dead time for 70 µs inhibit pulse. The second stage with ∼ 0.1 MeV threshold and ± 6 µs inhibit pulse preformed a selection among the BGO impulses for energies exceeding 2.3 MeV. This energy limit was established in order to make an estimate of the second stage dead time using the 2.614 MeV natural radioactivity peak. Our estimate to the active shielding dead time is 23%.

RESULTS
The measurements were performed over 29.8 days in live time by 2-hour series. The measurements were divided into series in order to monitor the time stability of the BGO detector and the active shielding. The energy spectrum of the BGO detector in the range of (0-11) MeV is shown in Fig.  2. The spectrum of the BGO signals that were not accompanied by the active shielding signal is designated as 1.
In the spectrum, one can identify two pronounced peaks at 1.460 MeV and 2.614 MeV; these are due to the natural radioactivity of the 40 K (located in the PMT's glass housing) and of 208 Tl from the 232 Th family (Fig. 2). The positions and intensities of these peaks were used for monitoring of time stability.
Bismuth has the largest nuclear charge among the stable isotopes (Z = 83), and the cross section of (e + e − )-pair production upon the interaction of γ quanta is therefore the largest for this element. The annihilation peak at 0.511 MeV is pronounced in the spectrum. The peak at 2.10 MeV is related to the escape of one annihilation γ quantum from the detector upon the detection of 2.614 MeV γ rays. The visible kink at ≈ 7.5 MeV is due to the γ quanta produced as a result of the capture of thermal neutrons by the components of the passive shield.
The positions and dispersion of the 1.46 MeV and 2.614 MeV peaks determined during the measurements were used to find the energy scale and resolution of the BGO detector. The energy calibration of the spectrometric channel was found as a linear function: E = A × N + B, where E is the released energy and N is the channel number. For higher energies the energy calibration was checked with a 239 Am-9 Be neutron source. The reaction 9 Be(α, n) 12 C * produce γquanta with energy 4.439 MeV corresponding to the energy of the first exited states 12 C nuclei. The position of 4.439 MeV peak is restored with accuracy better than 5 keV while 1461 keV and 2614 keV calibration peaks are used. (Fig.2,  inset).
The dependence of the energy resolution of a scintillation detector vs energy can be written as σ = C × √ E. The parameter C was found to be 0.04 MeV 1/2 . The values of σ determined from the background spectrum are in good agreement with the measurements performed with 60 Co, 207 Bi and 239 Pu 9 Be standard calibration sources. The expected standard deviation of the 5.5 MeV peak due to the axion absorption is σ = 93 keV. Figure 3 shows the energy range of (4.5 − 6.5) MeV, in which the axion peak was expected. The spectrum measured in the range of (4.5 − 6.5) MeV was fitted by a sum of exponential and two Gaussian functions: Here a, b and c are parameters of the function describing the smooth background. The position and dispersion of the first Gaussian peak corresponded to the desired-peak parameters: E 1 = 5.49 MeV is the axion peak position, σ 1 = 0.093 MeV is the Gaussian peak standard deviation. Because a small unknown peak can be seen at ≈ 5.8 MeV, the second Gaussian was added to the fitting function. The position and area of the second peak were free, while the dispersion σ 2 = 0.095 MeV was fixed.
The position of first peak (E 1 ) and dispersion (σ 1 ) were fixed and six parameters were varied, three of which described the continuous background while three others described area of the two peaks (S 1 , S 2 ), and the second peak position (E 2 ). The total number of the degrees of freedom in the range of (4.5-6.5) MeV was 95.
The fit results, corresponding to the minimum χ 2 = 99.4/95, are shown in Fig. 3. The position and area of the second peak are: E 2 = 5.77 ± 0.06 MeV and S 2 = 97 ± 55 counts. We attribute this peak to the intense 5.824 MeV gamma-ray line resulting in the capture of thermal neutrons by 113 Cd [27]. The cross section of thermal neutron capture is 2 × 10 4 barns. The intensity of the 5.49 MeV peak was found to be S 1 = −18 ± 58, this corresponds to the upper limit on the number of counts in the peak, S lim = 85 at a 90% confidence level [28].
The expected number of axioelectric absorption events are: where σ Ae is the axioelectric effect cross section, given by expression (5); Φ A is the axion flux (4); N Bi = 4.76 × 10 24 is the number of Bi atoms; T = 2.57 × 10 6 s is the measurement time; and ε = 0.67 is the detection efficiency for 5.5 MeV electrons. The axion flux Φ A is proportional to the constant (g 3 AN ) 2 , and the cross section σ Ae is proportional to the constant g 2 Ae , according to expressions (4) and (5). As a result, the S abs value depends on the product of the axionelectron and axion-nucleon coupling constants: (g Ae ) 2 ×(g 3 AN ) 2 .  Fig. 4 The limits on the g Ae coupling constant obtained by 1-present work, 2 -present work for |g Ae × g 3 AN |, 3-solar [12] and reactor experiments [15,29], 4-beam dump experiments [30,31]. The allowed |g Ae | and |g Ae × g 3 AN | values lie below the corresponding curves. The relations between g Ae and m A for DFSZ-and KSVZ-models are also shown .
The experimentally found condition S abs ≤ S lim imposes some constraints on the range of possible |g Ae × g 3 AN | and m A values. The range of excluded |g Ae × g 3 AN | values is shown in Fig. 4, at m A → 0 the limit is The dependence of |g Ae × g 3 AN | on m A is related only to the kinematic factor in formulae (2) and (5). These constraints are completely model-independent and valid for any pseudoscalar particle with coupling |g Ae | less than 10 −6(4) .
Within the hadronic axion model, g 3 AN and m A quantities are related by expression (3), which can be used to obtain a constraint on the g Ae constant, depending on the axion mass (Fig. 4). For m A = 1 MeV, this constraint corresponds to |g Ae | ≤ 1.4 × 10 −7 . Figure 4 also shows the constraints on the constant |g Ae | that were obtained in the Borexino experiment for 478-keV 7 Li solar axions [12] and in the Texono reactor experiment for 2.2-MeV axions produced in the n + p → d + A reaction [15]. Recently, Borexino coll. reported new more stringent limits on g Ae coupling for 5.5 MeV solar axions [13]. Unlike our work, these limits on g Ae were obtained in assumption that the axion interacts with electron through the Compton conversion process.
The sensitivity of the experiment to the constant g Ae depends on the target mass M, specific background level B, detector resolution σ , and measurement time T : g lim Ae ∼ (σ B/MT ) 1/2 . The level g Ae ≈ 10 −10 can be achieved with longer measurements with the detector mass enlarged by three orders of magnitude and the background level reduced by two orders of magnitude, the latter can be done by placing the setup in an underground laboratory.

CONCLUSIONS
A search for the axioelectric absorption of 5.5 MeV axions produced in the p + d → 3 He + γ reaction was performed using a BGO detector with a mass of 2.5 kg, located in a low-background setup equipped with passive and active shielding. As a result, a model-independent limit on axionnucleon and axion-electron coupling constant has been obtained: |g Ae ×g 3 AN | < 2.9×10 −9 (90% c.l.). Within the hadronic axion model the constraints on the axion-electron coupling constant |g Ae | ≤ (1.4 − 9.7) × 10 −7 for axions with masses 0.1 < m A < 1 MeV were obtained for 90% c.l..