An overview of the current theoretical studies on neutrino-atom scattering processes is presented. The ionization channel of these processes, which is studied in experiments searching for neutrino magnetic moments, is brought into focus. Recent developments in the theory of atomic ionization by impact of reactor antineutrinos are discussed. It is shown that the stepping approximation is well applicable for the data analysis practically down to the ionization threshold.
1. Introduction
In particle physics, the neutrino plays a remarkable role of a “tiny” particle. The scale of neutrino mass mν is much lower than that of the charged fermions (mνf≪mf, f=e,μ,τ). Interaction of neutrinos with matter is extremely weak as compared to that in the case of other known elementary fermions. According to the Standard Model (SM), it can be mediated only via exchange of the W± and Z0 bosons. However, the recent development of our knowledge of neutrino mixing and oscillations, supported by the discovery of flavor conversions of neutrinos from different sources (see [1–4]), substantiates the assumption that neutrinos can possess electromagnetic properties and, hence, take part in electromagnetic interactions (see, e.g., the review articles [5–7]). These properties include, in particular, the electric charge, the charge radius, the anapole moment, and the dipole electric and magnetic moments. Such nontypical neutrino features are of particular interest, because they open a door to “new physics” beyond the SM (BSM). In spite of appreciable efforts in searches for neutrino electromagnetic characteristics, up to now there is no experimental evidence favoring their nonvanishing values.
Among the electromagnetic properties of neutrinos, the most studied and well understood theoretically are neutrino magnetic moments (NMM), along with electric dipole moments. For the most recent and complete review on theoretical and experimental aspects of NMM, as well as for the corresponding references, see [7]. The effective Lagrangian, which describes the coupling of NMM to the electromagnetic field Fαβ, can be written in the form
(1)Lint=12ψ¯iσαβ(μij+ϵijγ5)ψjFαβ+h.c.,
where the magnetic moments μij, in the presence of mixing between different neutrino states, are associated with the neutrino mass eigenstates νi. The interplay between the magnetic moment and neutrino mixing effects is important. Note that the electric (transition) moments ϵij do also contribute to the coupling. A Dirac neutrino may have nonzero diagonal electric moments in models where CP invariance is violated. For a Majorana neutrino the diagonal magnetic and electric moments are zero. Therefore, NMM can be used to distinguish between Dirac and Majorana neutrinos [8–10].
In the Standard Model the magnetic moment of a massless neutrino is zero. In the minimal extension of the SM, the explicit evaluation of the one-loop contributions to the Dirac NMM in the leading approximation over small parameters bi=mi2/MW2 (mi are the neutrino masses, i=1,2,3) that however exactly accounts for the parameters al=ml2/MW2 (l=e,μ,τ) yields the following result [11–14]:
(2)μijD=eGFmi82π2(1+mjmi)∑l=e,μ,τf(al)UljUli*,
where Uli is the neutrino mixing matrix, and
(3)f(al)=34[1+11-al-2al(1-al)2-2al2lnal(1-al)3].
A Majorana neutrino can also have nondiagonal (transition) magnetic moments μijM=2μijD (i≠j). The obtained value for NMM is proportional to the neutrino mass and is, in general, of the order ~10-21–10-19μB.
Much larger values for NMM can be obtained in different other SM extensions (see [6, 7] for the detailed discussion). However, there is a problem [15] for any BSM theory of how to get a large NMM value and simultaneously to avoid an unacceptable large contribution to the neutrino mass. Recently, this problem has been reconsidered for a class of BSM theories and it has been shown in a model-independent way that in principle it is possible to avoid the above mentioned contradiction in the case of Dirac [16] and Majorana [17] neutrinos. It has been shown that in this kind of theoretical models the NMM can naturally reach values of ~10-15–10-14μB. These values are at least two orders of magnitude smaller than the present laboratory experimental limits (see below). There is also a huge gap of many orders of magnitude between these values and the prediction of the minimal extension of the SM. Therefore, if any direct experimental confirmation of nonzero NMM is obtained in the laboratory experiments, it will open a window to “new physics.”
The neutrino magnetic moments are being searched in reactor [18, 19], accelerator [20, 21], and solar [22, 23] experiments on low-energy elastic neutrino-electron scattering (for more details see the review articles [6, 7] and references therein). The current best upper limit on the NMM value obtained in such direct laboratory measurements is
(4)μν≤2.9×10-11μB,
where μB=e/(2me) is a Bohr magneton. This bound, which is due to the GEMMA experiment [19] with a HPGe detector at Kalinin nuclear power station, is by an order of magnitude larger than the tightest constraint obtained in astrophysics [24]:
(5)μν≤3×10-12μB.
And it by many orders of magnitude exceeds the value derived in the minimally extended SM that includes right-handed neutrinos [12, 14]:
(6)μν≤3×10-19μB(mν1eV),
where mν is a neutrino mass. At the same time, there are different theoretical BSM scenarios that predict much higher μν values. For example, the effective NMM value in a class of extra-dimension models can be as large as about 10-10μB [25]. Future higher precision reactor experiments can therefore be used to provide new constraints on large extra-dimensions.
The paper is organized as follows. Section 2 outlines the current status of searches for NMM and the problem of atomic-ionization effects in reactor experiments. Section 3 is devoted to the theoretical background for neutrino scattering on atomic electrons. In Section 4, we discuss the case of neutrino scattering on one bound electron. Hydrogen-like states and a semiclassical limit are considered. Section 5 focuses on ionization of many-electron atoms by neutrino impact. The case of a helium atomic target and the Thomas-Fermi and ab initio approaches are discussed. Finally, Section 6 summarizes this review.
2. Searches for Neutrino Magnetic Moments of Reactor Antineutrinos
The strategy of experiments searching for NMM is as follows. One studies an inclusive cross section for elastic (anti)neutrino-electron scattering which is differential in the energy transfer T. In the ultrarelativistic limit mν→0, it is given by an incoherent sum of the SM contribution dσSM/dT, which is due to weak interaction that conserves the neutrino helicity, and the helicity-flipping contribution dσ(μ)/dT, which is due to μν,
(7)dσdT=dσSMdT+dσ(μ)dT.
The SM term is well documented and is given by [26]
(8)dσSMdT=GF2me2π×[(gV+gA)2+(gV-gA)2(1-TEν)2i+(gA2-gV2)meTEν2],
where Eν is the incident antineutrino energy, gA=1/2 and gV=(4sin2θW+1)/2 for νe, and gA=-1/2 and gV=(4sin2θW-1)/2 for νμ and ντ, with θW being the Weinberg angle. For antineutrinos one must substitute gA→-gA.
The possibility for neutrino-electron elastic scattering due to NMM was first considered in [27], and the cross section of this process was calculated in [28] (the related brief historical notes can be found in [29]). Here we would like to recall the paper by Domogatsky and Nadezhin [30], where the cross section of [28] was corrected and the antineutrino-electron cross section was considered in the context of the earlier experiments with reactor antineutrinos of Cowan and Reines [31] and Cowan et al. [32], which were aimed to reveal the NMM effects. Discussions on the derivation of the cross section and on the optimal conditions for bounding the NMM value, as well as a collection of the cross section formulas for elastic scattering of neutrinos (antineutrinos) on electrons, nucleons, and nuclei, can be found in [29, 33]. The result relevant to the μν component in (7) reads [29, 30, 33]
(9)dσ(μ)dT=4παμν2(1T-1Eν),
where α is the fine-structure constant. Thus, the two components of the cross section (7) exhibit qualitatively different dependencies on the recoil-electron kinetic energy T. Namely, at low T values the SM cross section is practically constant in T, while that due to μν behaves as 1/T. This means that the experimental sensitivity to NMM value critically depends on lowering the energy threshold of the detector employed for measurement of the recoil-electron spectrum.
The current reactor experiments with germanium detectors [18, 19] have reached threshold values of T as low as few keV and are to further improve the sensitivity to low-energy deposition in the detector [34–36]. At low energies, however, one can expect a modification of the free-electron formulas due to the binding of electrons in the germanium atoms, where, for example, the energy of the Kα line, 9.89 keV, indicates that at least some of the atomic binding energies are comparable to the already relevant to the experiment values of T. Thus a proper treatment of the atomic effects in neutrino scattering is necessary and important for the analysis of the current and, even more, of the future data with a still lower threshold. Furthermore, there is no known means of independently calibrating experimentally the response of atomic systems, such as the germanium, to the scattering due to the interactions relevant for the neutrino experiments. Therefore, one has to rely on a pure theoretical analysis in interpreting the neutrino data. For the first time this problem was addressed in [37], where a 2-3-time enhancement of the electroweak cross section in the case of ionization from a 1s state of a hydrogen-like atom with nuclear charge Z had been numerically determined at neutrino energies Eν~αZmec2. Subsequent numerical calculations within the relativistic Hartree-Fock method for ionization from inner shells of various atoms showed much lower enhancement (~5–10%) of the electroweak contribution [38–43]. It was found that in the scattering on realistic atoms, such as germanium, the so-called stepping approximation works with a very good accuracy. The stepping approach, introduced in [40] from an interpretation of numerical data, treats the process as scattering on individual independent electrons occupying atomic orbitals and suggests that the cross section follows the free-electron behavior in (8) and (9) down to T equal to the ionization threshold for the orbital and that below that energy the electron on the corresponding orbital is “deactivated” thus producing a sharp “step” in the dependence of the cross section on T.
The interest in the role of atomic effects was renewed in several more recent papers. The early claim [44] of a significant enhancement of the NMM contribution in the case of germanium due to the atomic effects has been later disproved [45, 46] and it was argued [47–50] that the modification of the free-electron formulas (8) and (9) by the atomic-binding effects is insignificant down to very low values of T. This conclusion appeared to be also in contradiction to the results of [51], where it was deduced by means of numerical calculations that the μν contribution to ionization of the helium atomic target by impact of electron antineutrinos strongly enhances relative to the free-electron approximation. However, from calculations performed in [52] it follows that the stepping approximation is well applicable practically down to the ionization threshold for helium.
3. General Theoretical Framework
As indicated in the introduction, the most sensitive and widely used method for the experimental investigation of the neutrino electromagnetic properties is provided by direct laboratory measurements of low-energy elastic scattering of neutrinos and antineutrinos with electrons in reactor, accelerator, and solar experiments. In this section, we deliver a theoretical background for such studies.
3.1. Neutrino-Electron Interactions
Let us consider the elastic-scattering process
(10)ν+e-⟶ν+e-,
where an incident neutrino with energy Eν transfers to a free electron, which is initially at rest in the laboratory frame, the energy-momentum q. There are two recoil-electron observables: the kinetic energy T, which amounts to the energy transfer, and the outgoing angle χ measured with respect to the incident neutrino direction. In the ultrarelativistic limit mν=0, these kinematical variables are related by
(11)cosχ=Eν+meEνTT+2me.
The maximal value of the kinetic electron energy is thus realized when χ=0° and is given by
(12)Tmax=2Eν22Eν+me.
Within the SM, the scattering process (10) takes place due to exchange of the weak bosons, as shown in Figure 1. The W-boson channel corresponds to the charged current interaction and is absent in the cases of the muon and tau neutrinos. If |q2|≪mW2, where mW is the W-boson mass, the scattering amplitude is given by [26]
(13)MW=GF2u¯ν2γα(1-γ5)uν1u¯e2γα(1-γ5)ue1,
where uν1 (ue1) and uν2 (ue2) are initial and final neutrino (electron) spinors. The Z0 boson mediates the neutral current interaction. The corresponding scattering amplitude in the case |q2|≪mZ2, where mZ is the Z0-boson mass, reads [26]
(14)MZ=GF2u¯ν2γα(1-γ5)uν1u¯e2γα(gV-gAγ5)ue1.
Using the matrix elements (13) and (14), one arrives, after averaging over the initial and summing over the final electron spins, at the SM single-differential cross section (8).
Elastic neutrino-electron scattering due to the weak interaction. Exchange by the W (a) and Z0 (b) bosons is shown.
Figure 2 shows the electromagnetic channel of the scattering process (10). In general, the matrix element of the neutrino electromagnetic current Jμ can be considered between different neutrino initial ψi(p) and final ψj(p′) states of different masses, p2=mi2 and p′2=mj2,
(15)〈ψj(p′)|Jμ|ψi(p)〉=u¯j(p′)Λμ(q)ui(p).
In the most general case consistent with Lorentz and electromagnetic gauge invariance, the vertex function is defined as (see [5, 6] and references therein)
(16)Λμ(q)=[fQ(q2)ij+fA(q2)ijγ5](q2γμ-γμq)+fM(q2)ijiσμνqν+fE(q2)ijiσμνqνγ5,
where fQ(q2), fA(q2), fM(q2), and fE(q2) are, respectively, the charge, anapole, dipole magnetic, and dipole electric neutrino form factors, which are matrices in the space of neutrino mass eigenstates [14].
Contribution of the neutrino electromagnetic vertex function to neutrino elastic scattering on a charged lepton [5].
Consider the diagonal case i=j. The hermiticity of the electromagnetic current and the assumption of its invariance under discrete symmetries’ transformations put certain constraints on the form factors, which are in general different for the Dirac and Majorana neutrinos. In the case of Dirac neutrinos, the assumption of CP invariance combined with the hermiticity of the electromagnetic current Jμ implies that the electric dipole form factor vanishes, fE=0. At zero momentum transfer only fQ(0) and fM(0), which are called the electric charge and the magnetic moment, respectively, contribute to the Hamiltonian Hint~JμAμ that describes the neutrino interaction with the external electromagnetic field Aμ. The hermiticity also implies that fQ, fA, and fM are real. In contrast, in the case of Majorana neutrinos, regardless of whether CP invariance is violated or not, the charge, dipole magnetic, and electric moments vanish, fQ=fM=fE=0, so that only the anapole moment can be nonvanishing among the electromagnetic moments. Note that it is possible to prove [8–10] that the existence of a nonvanishing magnetic moment for a Majorana neutrino would bring about a clear evidence for CPT violation.
In the off-diagonal case i≠j, the hermiticity by itself does not imply restrictions on the form factors of Dirac neutrinos. It is possible to show [8] that if the assumption of the CP invariance is added, the form factors fQ, fM, fE, and fA should have the same complex phase. For the Majorana neutrino, if CP invariance holds, there could be either a transition magnetic or a transition electric moment. Finally, as in the diagonal case, the anapole form factor of a Majorana neutrino can be nonzero.
The neutrino dipole magnetic and electric form factors (and the corresponding magnetic and electric dipole moments) are theoretically the most well-understood among the form factors. The value of the magnetic form factor fM(q2) at q2=0 defines the NMM, μν=fM(0). In the low-energy limit, the NMM contribution to the effective electromagnetic vertex can be expressed in the following form:
(17)Λα=μν2meσαβqβ.
Thus, the corresponding scattering amplitude is
(18)M(μ)=4πμνα2meq2u¯ν2σαβqβuν1u¯e2γαue1.
This leads to the NMM single-differential cross section given by (9).
3.2. Neutrino Scattering on Atomic Electrons
Consider the process where a neutrino with energy-momentum pν=(Eν,pν) scatters on an atom at energy-momentum transfer q=(T,q). In what follows the recoil of the atomic nucleus is neglected because of the typical of current experiments situation T≫2Eν2/MN, where MN is the nuclear mass. The atomic target is supposed to be unpolarized and in its ground state |0〉 with the corresponding energy E0. It is also supposed that T≪me and αZ≪1, where Z is the nuclear charge and α is the fine-structure constant, so that the initial and final electronic systems can be treated nonrelativistically. The neutrino states are described by the Dirac spinors assuming mν=0.
Thus, the magnetic moment interaction of the neutrino field ψ with the atomic electrons is described by the Lagrangian
(19)Lint=μν2meψ¯(pν′)σαβψ(pν)qαAβ,
where pν′ is the final neutrino four-momentum. The electromagnetic field A=(A0,A) of the atomic electrons is A0(q)=4παρ(q)/q2, A(q)=4παj(q)/q2 (hereafter we use the notation q=|q|), where ρ(q) and j(q) are the Fourier transforms of the electron number density and current density operators, respectively,
(20)ρ(q)=∑a=1Zexp(iq·ra),(21)j(q)=-i2m∑a=1Z[exp(iq·ra)∂∂ra+∂∂raexp(iq·ra)],
and the sums run over the positions ra of all the Z electrons in the atom. The double-differential cross section can be presented as
(22)d2σ(μ)dTdq2=(d2σ(μ)dTdq2)∥+(d2σ(μ)dTdq2)⊥,
where
(23)(d2σ(μ)dTdq2)∥=4παμν2q2(1-T2q2)S(T,q2),(24)(d2σ(μ)dTdq2)⊥=4παμν2q2(1-q24Eν2)R(T,q2),
where S(T,q2), also known as the dynamical structure factor [53], and R(T,q2) are
(25)S(T,q2)=∑nδ(T-En+E0)|〈n|ρ(q)|0〉|2,(26)R(T,q2)=∑nδ(T-En+E0)|〈n|j⊥(q)|0〉|2,
with j⊥ being the j component perpendicular to q and parallel to the scattering plane, which is formed by the incident and final neutrino momenta. The sums in (25) and (26) run over all the atomic states |n〉 with energies En of the electron system, with |0〉 being the initial state.
The longitudinal term (23) is associated with atomic excitations induced by the force that the neutrino magnetic moment exerts on electrons in the direction parallel to q. The transverse term (24) corresponds to the exchange of a virtual photon which is polarized as a real one, that is, perpendicular to q. It resembles a photoabsorption process when q→T and the virtual-photon four-momentum thus approaches a real-photon value. Due to selections rules, the longitudinal and transverse excitations do not interfere (see [54] for detail).
The factors S(T,q2) and R(T,q2) are related to, respectively, the density-density (or polarization) and current-current Green’s functions
(27)S(T,q2)=1πImF(T,q2),R(T,q2)=1πImL(T,q2),
where
(28)F(T,q2)=∑n|〈n|ρ(q)|0〉|2T-En+E0-iϵ=〈0|ρ(-q)1T-H+E0-iϵρ(q)|0〉,(29)L(T,q2)=∑n|〈n|j⊥(q)|0〉|2T-En+E0-iϵ=〈0|j⊥(-q)1T-H+E0-iϵj⊥(q)|0〉,H being the Hamiltonian for the system of electrons. From these relations it follows that, due to the parity selection rule, the functions S(T,q2) and R(T,q2) are even with respect to q.
For small q values, in particular, such that q~T, only the lowest-order nonzero terms of the expansion of (27) in powers of q2 are of relevance (the so-called dipole approximation). In this case, one has [45, 47]
(30)R(T,q2)=T2q2S(T,q2).
Taking into account (30), the experimentally measured singe-differential inclusive cross section is, to a good approximation, given by (see, e.g., [47, 49, 50])
(31)dσ(μ)dT=4παμν2∫T2(2Eν-T)2S(T,q2)dq2q2.
The standard electroweak contribution to the cross section can be similarly expressed in terms of the same factor S(T,q2) [45, 50] as
(32)dσSMdT=GF24π(1+4sin2θW+8sin4θW)×∫T2(2Eν-T)2S(T,q2)dq2,
where the factor S(T,q2) is integrated over q2 with a unit weight, rather than q-2 as in (31).
The kinematical limits for q2 in an actual neutrino scattering are explicitly indicated in (31) and (32). At large Eν, typical for the reactor neutrinos, the upper limit can in fact be extended to infinity, since in the discussed here nonrelativistic case the range of momenta q~Eν is indistinguishable from infinity on an atomic scale. The lower limit can be shifted to q2=0, since the contribution of the region of q2<T2 can be expressed in terms of the photoelectric cross section [45] and is negligibly small (at the level of below one percent in the considered range of T). For this reason one can discuss the momentum-transfer integrals in (31) and (32) running from q2=0 to q2=∞:
(33)I1(T)=∫0∞S(T,q2)dq2q2,I2(T)=∫0∞S(T,q2)dq2.
For a free electron, which is initially at rest, the density-density correlator is the free particle Green’s function
(34)F(FE)(T,q2)=(T-q22me-iϵ)-1,
so that the dynamical structure factor is given by
(35)S(FE)(T,q2)=δ(T-q22me),
and the discussed here integrals are in the free-electron limit as follows:
(36)I1(FE)=∫0∞S(FE)(T,q2)dq2q2=1T,I2(FE)=∫0∞S(FE)(T,q2)dq2=2me.
Clearly, these expressions, when used in the formulas (31) and (32), result in the free-electron cross sections for the case T≪Eν,
(37)dσ(μ)dT=4παμν2T,dσSMdT=GF2me2π(1+4sin2θW+8sin4θW),
correspondingly.
4. Scattering on One Bound Electron
In this section, we consider neutrino scattering on an electron bound in an atom following consideration of [47, 49, 50]. The binding effects generally deform the density-density Green’s function, so that both the integrals (33) are somewhat modified. Namely, the binding effects spread the free-electron δ-peak in the dynamical structure function (35) at q2=2meT and also shift it by the scale of characteristic electron momenta in the bound state.
4.1. Ionization from a Hydrogen-Like Orbital
Consider the situation when the initial electron occupies the discrete nl orbital in a Coulomb potential V(r)=-αZ/r. The dynamical structure factor for this hydrogen-like system is given by
(38)S(nl)(T,q2)=mek(2π)312l+1∑m=-ll∫dΩk|〈φk-|ρ(q)|φnlm〉|2,
where φnlm is the bound-state wave function, φk- is the outgoing Coulomb wave for the ejected electron with momentum k, and k=|k|=2meT-pn2, with pn=αZme/n being the electron momentum in the nth Bohr orbit. The closed-form expressions for the bound-free transition matrix elements in (38) can be found, for instance, in [55]. In principle, they allow for performing angular integrations in (38) analytically. This task, however, turns out to be formidable for large values of n. Therefore, below we restrict our consideration to the n=1,2 states only, which nevertheless is enough for demonstrating the validity of the semiclassical approach described in Section 4.2.
Using results of [56], we can present the function (38) when n=1,2 as
(39)S(nl)(T,q2)=28mepn63[1-exp(-2πη)]×q2fnl(q2)[(q2-k2+pn2)2+4pn2k2]2n+1×exp[-2ηarctan(2pnkq2-k2+pn2)],
where the branch of the arctangent function that lies between 0 and π should be used, η=αZme/k is the Sommerfeld parameter, and
(40)f1s(q2)=3q2+k2+p12,(41)f2s(q2)=8[(p22+k2)4,(p22+k2)43q10-(32p22+11k2)q8+(82p24+72p22k2+14k2)q6+(20p26-62p24k2-20p22k4-6k6)q4+(p22+k2)×(475p26-475p24k2-7p22k4-k6)q2+(4p22+k2)(p22+k2)4(p22+k2)4(p22+k2)4],(42)f2p(q2)=2p22[(443p22+4k2)36q8-48(p22+k2)q6+(152p24-48p22k2-8k4)q4+(p22+k2)×(171215p24+156815p22k2+16k4)q2+(443p22+4k2)(p22+k2)3].
Figure 3 shows the magnetic single-differential cross section (31) for ionization from the 1s orbital, which is normalized to the free-electron value (9), as a function of T/ɛb, with the electron binding energy given by ɛb=α2Z2me/2. As can be seen, the numerical results for Eν≫ɛb are close to the free-electron ones in magnitude. This can be qualitatively explained by noticing the following facts. First, in an attractive Coulomb potential there is an infinite set of bound states, with the discrete spectrum smoothly transforming into the continuum at the ionization threshold. Second, the average value of the 1s electron momentum is pe=p1 and the average change in the electron momentum after ejection, Δpe, is such that Δpe2=2meT, which is analogous to the free-electron case.
The ratio of single-differential cross sections for magnetic neutrino scattering from the 1s hydrogen-like and free-electron states, respectively, versus T/ɛb at different values of Eν. The Eν=50ɛb and Eν=100ɛb curves are practically indistinguishable.
Thus, taking into account the results in Figure 3, one might expect the atomic-binding effects to play a subsidiary role when Eν≫ɛb. The authors of [44], however, came to the contrary conclusion that the single-differential cross section dramatically enhances due to atomic ionization when T~ɛb. The enhancement mechanism proposed in [44] is based on an analogy with the photoionization process. As mentioned above, when q→T the virtual-photon momentum approaches the physical regime T2-q2=0. In this limit, we have for the transverse component of the double-differential cross section (24)
(43)(d2σ(μ)dTdq2)⊥=μν2πσγ(T)T,
where σγ(T) is the photoionization cross section at the photon energy T [57]. The limiting form (43) was used in [44] in the whole integration interval, when deriving the single-differential cross section. Such procedure is obviously incorrect, for the integrand rapidly falls down as q2 ranges from T2 up to (2Eν-T)2, especially when q2≳ra-2, where ra is a characteristic atomic size (within the Thomas-Fermi model ra-1=Z1/3αme [58]). This fact reflects a strong departure from the real-photon regime. For this reason we can classify the enhancement of the differential cross section determined in [44] as spurious.
Insertion of (39) into the integrals (33) and integration over q2, using the change of variable
(44)2pnkq2-k2+pn2=tanx
and the standard integrals involving the products of the exponential function and the powers of sine and cosine functions, yield [50]
(45)I1(1s)(T)=I2(1s)(T)2meT=T-11-exp(-2π/y1-1)×{1-exp(-πy1-1)h×exp[-2y1-1arctan(y1-22y1-1)]h×(1-4y1+163y12)},(46)I1(2s)(T)=T-11-exp(-4π/y2-1)×{1-exp(-2πy2-1)h×exp[-4y2-1arctan(y2-22y2-1)]h×(1-8y2+803y22-44815y23+179215y24)},(47)I2(2s)(T)=2me1-exp(-4π/y2-1)×{1-exp(-2πy2-1)h×exp[-4y2-1arctan(y2-22y2-1)]h×(1-8y2+803y22-44815y23+102415y24)},(48)I1(2p)(T)=T-11-exp(-4π/y2-1)×{1-exp(-2πy2-1)h×exp[-4y2-1arctan(y2-22y2-1)]h×(1-8y2+803y22-70415y23+332845y24)},(49)I2(2p)(T)=2me1-exp(-4π/y2-1)×{1-exp(-2πy2-1)h×exp[-4y2-1arctan(y2-22y2-1)]h×(1-8y2+803y22-70415y23+51215y24)},
where yn=2meT/pn2≡T/|En|. The largest deviations of these integrals from the free-electron analogs (36) occur at the ionization threshold T=|En|. The corresponding relative values in this specific case are [50]
(50)I1(1s)I1(FE)=I2(1s)I2(FE)=1-73e-4=0.9572635093,(51)I1(2s)I1(FE)=1-163915e-8=0.9633451168,I2(2s)I2(FE)=1-87115e-8=0.9805208034,I1(2p)I1(FE)=1-210145e-8=0.9843376226,I2(2p)I2(FE)=1-10315e-8=0.9976964900.
The above results indicate a clear tendency: the larger the n and l are, the closer the I1(nl) and I2(nl) are to the free-electron values. The departure from the free-electron behavior does not exceed several percent at most. These observations provide a solid base for the semiclassical approach described below.
4.2. Semiclassical Approach
In the one-electron approximation, the Hamiltonian has the form H=p2/2me+V(r), and the density-density Green’s function from (28) can be written as
(52)F(T,q2)=〈0|e-iq·r[T-H(p,r)+E0]-1eiq·r|0〉=〈0|[T-H(p+q,r)+E0]-1|0〉=〈0|[T-q22me-p·qme-H(p,r)+E0]-1|0〉,
where the infinitesimal shift T→T-iϵ is implied.
Clearly, a nontrivial behavior of the latter expression in (53) is generated by the presence of the operator (p·q) in the denominator and the fact that it does not commute with the Hamiltonian H. Thus an analytical calculation of the Green’s function as well as the dynamical structure factor is feasible in only few specific problems. In Section 4.1 such calculation has been presented for ionization from the 1s, 2s, and 2p hydrogen-like states. In particular, we have seen that the deviations of the discussed integrals (33) from their free values are very small: the largest deviation is exactly at the ionization threshold, where, for instance, each of the 1s integrals is equal to the free-electron value multiplied by the factor (1-7e-4/3)≈0.957 (see (50)). It can be also noted from (45) that both integrals are modified in exactly the same proportion, so that their ratio is not affected at any T: I2(T)/I1(T)=2meT. We find, however, that this exact proportionality is specific for the ionization from the ground state in the Coulomb potential.
The problem of calculating the integrals (33), however, can be solved in the semiclassical limit, where one can neglect the noncommutativity of the momentum p with the Hamiltonian and rather treat this operator as a number vector. Taking also into account that (H-E0)|0〉=0, one can then readily average the latter expression in (53) over the directions of q and find the formula for the dynamical structure factor:
(53)S(T,q2)=me2pq[θ(T-q22me+pqme)-θ(T-q22me-pqme)],
where θ is the standard Heaviside step function. The expression in (53) is nonzero only in the range of q satisfying the condition -pq/me<T-q2/2me<pq/me, that is, between the (positive) roots of the binomials in the arguments of the step functions: qmin2=2meT+p2-p and qmax2=2meT+p2+p. One can notice that the previously mentioned “spread and shift” of the peak in the dynamical structure function in this limit corresponds to a flat pedestal between qmin2 and qmax2. The calculation of the integrals (33) with the expression (53) is straightforward and yields the free-electron expressions (36) for the discussed here integrals in the semiclassical (WKB) limit:
(54)I1(WKB)=1T,I2(WKB)=2me.
The appearance of the free-electron expressions here is not surprising, since (53) can be also viewed as the one for scattering on an electron boosted to the momentum p. The difference from the pure free-electron case however is in the range of the energy transfer T. Namely, the expressions (54) are applicable in this case only above the ionization threshold, that is, at T≥|E0|. Below the threshold the electron becomes “inactive.”
We believe that the latter conclusion explains the so-called stepping behavior observed empirically [40] in the results of numerical calculations. Namely, the calculated cross section dσ/dT for ionization of an electron from an atomic orbital follows the free-electron dependence on T all the way down to the threshold for the corresponding orbital with a very small, at most a few percent, deviation. This observation led the authors of [40] to suggest the stepping approximation for the ratio of the atomic cross section (per target electron) to the free-electron one:
(55)f(T)≡dσ/dT(dσ/dT)FE=1Z∑iniθ(T-|Ei|),
where the sum runs over the atomic orbitals with the binding energies Ei and the filling numbers ni. Clearly, the factor f(T) simply counts the fraction of “active” atomic electrons at the energy T, that is, those for which the ionization is kinematically possible. For this reason we refer to f(T) as an atomic factor. We conclude here that the stepping approximation is indeed justified with a high accuracy in the approximation of the scattering on independent electrons, that is, if one neglects the two-electron correlations induced by the interference of terms in the operator ρ(q) in (20) corresponding to different electrons. The effects of such interference will be discussed in the next section.
5. Scattering on Many-Electron Atoms
In considering the neutrino scattering on actual many-electron atoms one needs to evaluate the dependence of the number of active electrons on T and generally also evaluate the effect of the two-electron correlations. The latter can be studied, for example, in the case of a helium atom, where the electron-electron correlations are known to play a very significant role.
5.1. Helium
Recently, the authors of [51] deduced by means of numerical calculations that the μν contribution to ionization of the He atomic target by impact of electron antineutrinos from reactor and tritium sources strongly departures from the stepping approximation, exhibiting large enhancement relative to the free-electron case. According to [51], the effect is maximal when the T value approaches the ionization threshold in helium, TI=24.5874eV, where the relative enhancement is as large as almost eight orders of magnitude. It was thus suggested that this finding might have an impact on searches for μν, provided that its value falls within the range 10-13-10-12μB. In this section, following consideration of [52], we show that (i) the result of [51] is erroneous and (ii) the stepping approximation for helium is well applicable, except the energy region T~TI where the differential cross section substantially decreases relative to the free-electron case.
We consider the process where an electron antineutrino with energy Eν scatters on a He atom at energy and spatial-momentum transfers T and q, respectively. In what follows we focus on the ionization channel of this process in the kinematical regime T≪Eν, which mimics a typical situation with reactor (Eν~1 MeV) and tritium (Eν~10 keV) antineutrinos when the case T→TI is concerned. The He target is assumed to be in its ground state |Φi〉 with the corresponding energy Ei. Since for helium one has αZ≪1, where Z=2 is the nuclear charge, the state |Φi〉 can be treated nonrelativistically. As we are interested in the energy region T~TI, the final He state |Φf〉 (with one electron in continuum) can also be treated in the nonrelativistic approximation.
Under the above assumptions, the dynamical structure factor (25) is given by
(56)S(T,q2)=∑f|〈Φf(r1,r2)|eiq·r1+eiq·r2|Φi(r1,r2)〉|2×δ(T-Ef+Ei).
Here the f sum runs over all final He states having one electron ejected in continuum, with Ef being their energies.
For evaluation of the dynamical structure factor (56) we employ the same models of the initial and final He states as in [51]. The initial state is given by a product of two 1s hydrogen-like wave functions with an effective charge Zi,
(57)Φi(r1,r2)=φ1s(Zi,r1)φ1s(Zi,r2),φ1s(Zi,r)=Zi3πa03e-Zir/a0,
where a0=1/(αme) is the Bohr radius. The final state has the form
(58)Φf(r1,r2)=12[φk-(Zf,r1)φ1s(Z,r2)+φk-(Zf,r2)φ1s(Z,r1)],
where φk-(Zf,r) is an outgoing Coulomb wave for the ejected electron with spatial momentum k. Zf is the effective charge experienced by the ejected electron in the field of the final He+ ion. Contributions to the dynamical structure factor from excited He+ states are neglected due to their very small overlap with the K-electron state in the He atom.
To avoid nonphysical effects connected with nonorthogonality of states (57) and (58), we use the Gram-Schmidt orthogonalization
(59)|Φf〉⟶|Φf〉-〈Φi∣Φf〉|Φi〉.
Substitution of (57) and (58) into (56) thus yields
(60)S(T,q2)=∫dk(2π)3|F(k,q)|2×δ(T-k22me+2α2me-Zi2α2me),
where k=2me(T+2α2me-Zi2α2me), and
(61)F(k,q)=2〈φk-(Zf,r1)φ1s(Z,r2)|eiq·r1+eiq·r2-2ρ1s(q)|φ1s(Zi,r1)φ1s(Zi,r2)〉
is the inelastic form factor, with
(62)ρ1s(q)=∫φ1s(Zi,r)eiq·rφ1s(Zi,r)dr.
It is straightforward to perform the further calculation of the dynamical structure factor analytically (see, e.g., the textbook [58]). The resulting expression is
(63)S(T,q2)=216α4me5ZfZi6(1-e2πη)(2+Zi)6×[A1(k,q)+B(k,q)A2(k,q)+B2(k,q)],
where η=-αZfme/k and (introducing pi=αZime)
(64)A1(k,q)=((exp(2ηarccospi2+q2-k2(pi2+q2+k2)2-4k2q2))h×([(pi2+q2+k2)2-4k2q2]3)-1(exp(2ηarccospi2+q2-k2(pi2+q2+k2)2-4k2q2)))×{-ηpi(pi2+q2+k2)13](pi+ηk)2(pi2+q2+k2)2h+4kq2[13kpi2-23η2kpi2-η2k3hi-ηpi(pi2+q2+k2)13]},A2(k,q)=((2exp(ηarccospi2+q2-k2(pi2+q2+k2)2-4k2q2))h×((pi2+q2+k2)2-4k2q2)-1(2exp(ηarccospi2+q2-k2(pi2+q2+k2)2-4k2q2)))×[(η2ln(k+q)2+peff2(k-q)2+pi2)picos(η2ln(k+q)2+pi2(k-q)2+pi2)h+pi2-q2+k22qh×sin(η2ln(k+q)2+pi2(k-q)2+pi2)],B(k,q)=e2ηarctg(k/pi)(Zi-Zf)αme(k2+pi2)2h×{(2+Zi)4α4me4[(2+Zi)2α2me2+q2]2-32pi4(4pi2+q2)2}.
Finally, the usual choice of the effective charges is Zi=27/16≈1.69 and Zf=1 (see, e.g., [59] and references therein). The value Zi=27/16 follows from the variational procedure that minimizes the ground-state energy Ei, while the value Zf=1 ensures the correct asymptotic behavior of the final state. However, the authors of [51] utilized in their calculations the values Zi=1.79 and Zf=1.1 derived from fitting the photoionization cross section data on helium with the present model of the He states.
The departures of the differential cross sections (31) and (32) from the free-electron approximation are characterized by the respective atomic factors
(65)fSM=dσSM/dTdσSMFE/dT,fNMM=dσ(μ)/dTdσ(μ)FE/dT,
where dσSMFE/dT and dσ(μ)FE/dT are the SM and μν contributions to the differential cross section for scattering of an electron antineutrino on two free electrons. Let us recall that following [51] one should expect the fNMM value to be of about 108 at T→TI.
Numerical results for atomic factors (65) are shown in Figure 4. They correspond to the kinematical regime T≪αme≪2Eν, which is typically realized both for reactor and for tritium antineutrinos when T<200eV. Note that in such case one can safely set the upper limit of integrals in (31) and (32) to infinity, as the dynamical structure factor S(T,q2) rapidly falls down when q≳αme and practically vanishes in the region q≫αme. It can be seen from Figure 4 that atomic factors exhibit similar behaviors for both sets of the Zi and Zf parameters discussed in the previous section. Namely, their values are minimal (~0.5) at the ionization threshold, T=TI, and tend to unity with increasing T. The latter tendency is readily explained by approaching the free-electron limit. It can be also seen that a more or less serious deviation (>10%) of the present results from the stepping approximation is observed only in the low-energy region T<100 eV. This deviation can be attributed to the effect of the electron-electron correlation in a helium atom. Indeed, if the electrons do not interact with each other, then they occupy two 1s hydrogen-like states (with opposite spins), in which case the departure of the atomic-factor values from unity is, according to the results of Section 4.1, less than 5%.
Atomic factors (65) as functions of the energy transfer [52].
Thus, the calculations presented in Figure 4 do not confirm the huge enhancement of the μν contribution with respect to the free-electron approximation. Moreover, in accord with various calculations for other atomic targets [38–41, 43, 47, 49, 50], we find that at small energy-transfer values the electron binding in helium leads to the appreciable reduction of the differential cross section relative to the free-electron case. We attribute the erroneous prediction of [51] to the incorrect dynamical model that draws an analogy between the NMM-induced ionization and photoionization. Indeed, as discussed in Section 3.2, the virtual photon in the NMM-induced ionization process can be treated as real only when q→T. However, the integration in (31) involves the q values ranging from T up to 2Eν-T. Since Eν≫T, the real-photon picture appears to be applicable only in the vicinity of the lower integration limit. When moving away from that momentum region, one encounters a strong departure from the real-photon approximation which treats the integrand as a constant in the whole integration range, assuming it to be equal to its value at q=T; that is,
(66)1q2S(T,q2)=1T2S(T,T2).
Such an approach is manifestly unjustified, and it gives rise to the spurious enhancement of the μν contribution to the differential cross section.
5.2. Thomas-Fermi Model
In the Thomas-Fermi model (see, e.g., [58]) the atomic electrons are described as a degenerate free-electron gas in a master potential ϕ(r) filling the momentum space up to the zero Fermi energy, namely, up to the momentum p0(r) such that p02/2me-eϕ=0. The electron density n(r)=p03/(3π2) then determines the potential ϕ(r) from the usual Poisson’s equation. In the discussed picture at an energy transfer T the ionization is possible only for the electrons whose energies in the potential are above -T, that is, with momenta above pT(r) with pT2/2me-eϕ=-T. The electrons with lower energy are inactive. Calculating the density of the inactive electrons as pT3/(3π2) and subtracting their total number from Z, one readily arrives at the formula for the atomic factor, that is, the effective fraction of the active electrons Zeff/Z as a function of T,
(67)f(T)=Zeff(T)Z=1-∫0x0(T)[χ(x)x-TT0]3/2x2dx,
where χ(x) is the Thomas-Fermi function, well known and tabulated, of the scaling variable x=2(4/3π)2/3meαZ1/3, the energy scale T0 is given by
(68)T0=2(43π)2/3meα2Z4/3≈30.8Z4/3eV,
and, finally, x0(T) is the point where the integrand becomes zero, namely, corresponding to the radius beyond which all the electrons are active at the given energy T. The energy scale T0 in germanium (Z=32) evaluates to T0≈3.1keV. The Thomas-Fermi atomic factor for germanium calculated from the formula (67) is shown by the dashed line in the plot of Figure 5. The discussed statistical model is known to approximate the average bulk properties of the atomic electrons with a relative accuracy O(Z-2/3) and as long as the essential distances r satisfy the condition Z-1≪meαr≪1, which condition in terms of the scaling variable x reads as Z-2/3≪x≪Z1/3. In terms of the formula (67) for the number of active electrons, the lower bound on the applicability of the model is formally broken at T~Z2/3T0, that is, at the energy scale of the inner atomic shells. However, the effect of the deactivation of the inner electrons is small, of order Z-1 in comparison with the total number Z of the electrons. On the other hand, at low T, including the most interesting region of T~T0, the integral in (67) is determined by the range of x of order one, where the model treatment is reasonably justified.
The atomic factor f for germanium in the stepping approximation with the actual energies of the orbitals (solid line) and its interpolation in the Thomas-Fermi model (dashed line) [49].
The energies of the inner K, L, and M orbitals in the germanium atom are well known (see, e.g., [61]) and provide the necessary data for a description of the neutrino scattering by the stepping formula (55) down to the values of the energy transfer T in the range of the binding of the M electrons, that is, at T>|EM|≈0.18keV. The corresponding steps in the atomic factor are shown in Figure 5. One can see that the stepping atomic factor (55) mimics upon average over the energy intervals between the electron shells in germanium the Thomas-Fermi result. Thus, it can be considered as refinement of the latter due to accounting for the quantization of the electron binding energies. It can be mentioned that if one applies formulas of Section 4.1 to the onset of the K shell step, namely, just above 10.9 keV, the difference from the shown in the plot step function would be practically invisible in the scale of Figure 5.
5.3. Ab Initio Approaches
While the treatments based on a generic model of many-electron atomic targets allow determining characteristic features and behaviors of the differential cross sections (31) and (32), to obtain accurate numerical results one needs to resort to ab initio calculations. Such calculations can be realized using the Hartree-Fock (HF) method (see, e.g., [58]) and its modifications. In the HF approximation, atomic electrons occupy one-electron states in a spherically-symmetric mean-field potential which is derived self-consistently from the solution of the HF equations. Accordingly, each one-electron state independently contributes to the atomic-ionization process. For the first time this approach was formulated in [38, 39], where it was illustrated with numerical calculations of neutrino-impact ionization of the F and Mo atomic targets. The wave functions and energies of atomic bound states were calculated within the relativistic HF method [62, 63] with local exchange-correlation potential [64]. The wave functions of outgoing electrons were obtained by a numerical solution of the Dirac equation in the same mean-field potential as for the wave functions of discrete states. Performed in [40] numerical calculations for ionization of the iodine atoms by impact of reactor antineutrinos led the authors to suggest the stepping approximation (55).
In a very recent theoretical study [60], the authors adopted the multiconfiguration relativistic random-phase approximation (MCRRPA) [65, 66] to evaluate the germanium atomic factors. This particular method is based on the time-dependent HF approximation [67]; however, several important features make it a better tool beyond the usual HF approximation to describe transitions of open-shell atoms of high atomic number Z. First, for open-shell atoms, typically there is more than one configuration which has the desired ground-state properties. Therefore, a proper HF reference state should be formed by a linear combination of these allowed configurations, that is, a multiconfiguration reference state. Second, for atoms of high Z, the relativistic corrections can no longer be ignored. By using a Dirac equation, instead of a Schrödinger one, the leading relativistic terms in the atomic Hamiltonian are treated nonperturbatively from the onset. Third, two-body correlations in addition to the HF approximation are generally important for excited states and transition matrix elements. The random-phase approximation (RPA) is devised to account for part of the additional two-body correlations (particles can be in the valence or core states) not only for the excited but also for the reference state, and in a lot of cases, it gives good agreement with experiment [68]. Furthermore, it has been shown that RPA equations preserve gauge invariance [69]; this provides a measure of stability of their solutions.
The MCRRPA has been applied successfully to photoexcitation and photoionization of divalent atoms such as Be, Mg, Zn, and Sr (some of the results are summarized in [70]). Following similar treatments, the authors of [60] treated the electronic configuration of germanium as a core filled up to the 4s orbits, with two valence electrons in the 4p orbits. As the Ge ground state is a 3P0 state, it is a linear combination of two configurations, namely, [Zn] 4p1/22 and [Zn] 4p3/22. The wave function was calculated using the multiconfiguration Dirac-Fock package [71]. The atomic excitations due to weak and magnetic scattering were solved by the MCRRPA equation, and consequently transition matrix elements were yielded. In that calculation, all the current operators were expanded by spherical multipoles, and the resulting final scattering states were represented in the spherical wave basis and subject to the incoming-wave boundary condition.
Compared with the previous works on the same subject [38–41, 43] which are also in the similar spirit of the relativistic HF method, the MCRRPA approach differs in several respects. First, due to the near degeneracy of the NII(4p3/2) and NIII(4p1/2) levels, using a multiconfiguration reference state is necessary. Second, the nonlocal Fock term is treated exactly, without resorting to the local exchange potentials. Third, the excited states are calculated with two-body correlation built in by MCRRPA, not simply by solving a Coulomb wave function with a static one-hole mean field.
Figure 6 shows numerical results from [60] for ionization of germanium by impact of an electron antineutrino. As can be seen, in the energy region T≳1 keV the results are very well explained by the stepping-approximation formula (55). At the same time, in the sub-keV region, that is, where the electrons from the K and L shells in germanium stay “inactive,” both the SM and NMM contributions appear to be significantly suppressed as compared to the stepping approximation. The latter finding seems to disagree with the semiclassical approach discussed in Section 4.2, according to which the ionization involving more loosely-bound electron states, such as those belonging to the M and N shells, is expected to follow more closely the free-electron scenario. Notice that a similar suppression of the atomic-factor values close to the ionization threshold was observed in the case of helium (see Figure 4), and it was attributed to the two-electron correlation effect. Thus, we can suggest that the correlation effects beyond the approximation of independent electrons lead to the suppression of atomic factors in the low-energy region. This feature will be important for the next-generation experiments with Ge detectors having energy thresholds in the sub-keV region [34–36].
The SM (weak) and NMM contributions to the differential cross section of ν¯e-Ge ionization at Eν=1 MeV [60] in comparison with the corresponding stepping-approximation results (FEA). The NMM value is set to be the current upper limit μν=2.9×10-11μB [19].
6. Summary and Perspectives
In this review, we have considered the neutrino-atom ionizing collisions with focus on the most important theoretical issues related to the problem. The main results discussed in the paper can be summed up as follows.
The differential over the energy transfer cross section given by the free-electron formulas (8) and (9) and the stepping behavior of the atomic factor given by (55) provide a reasonable description of the neutrino-impact ionization of a complex atom, such as germanium, down to quite low-energy transfer. The deviation from this approximation due to the onset of the ionization near the threshold is less than 5% (of the height of the step) for the K electrons, if one applies the analytical behavior of this onset that one finds for the ground state of a hydrogen-like ion. It is also found that the free-electron expressions for the differential cross section are not affected by the atomic binding effects in the semiclassical limit and for independent electrons. These analytical results can support the numerically determined behaviors of the electroweak and magnetic contributions to the neutrino-impact ionization of various atomic targets within the mean-field model [38–41]. At the same time, very recent numerical calculations of the ν¯e-impact ionization processes of helium [52] and germanium [60] exhibit suppression of the SM and NMM differential cross sections relative to the stepping approximation with lowering the energy-transfer value. This suppression can be assigned to the electron-electron correlation effects.
A theoretical analysis [52] of ionization of helium by electron-antineutrino impact shows no evidence of the recently predicted enhancement [51] of the electromagnetic contribution as compared with the free-electron case. In contrast, in line with previous studies on other atomic targets, it is found that the magnitudes of the differential cross sections decrease relative to the free-electron approximation when the energy transfer is close to the ionization threshold. Thus, no sensitivity enhancement can be expected when using the He atomic target in searches for NMM. And the stepping approximation appears to be valid, within a few-percent accuracy, down to the energy-transfer values as low as almost 100 eV. We thus conclude that for practical applications, that is, for the analysis of data of the searches for NMM, one can safely apply the free-electron formulas and the stepping approximation at the energy transfer down to this range.
When analyzing the low-T data of the current high-sensitivity experiments searching for neutrino electromagnetic properties, one must go beyond the free-electron approximation for the elastic neutrino-electron scattering and take into account the atomic-ionization effects, at least, in the case of K electrons. At the present time, the experiment GEMMA-II with reactor antineutrinos is in preparation [19]. Its sensitivity to the NMM value is expected to be at the level of 1×10-11μB. To achieve such a sensitivity level, which is the region of astrophysical interest [24], it is planned to reduce the effective energy threshold of a Ge detector from 2.8 to 1.5 keV. This threshold value will be very close to the binding energies of the L electrons in germanium (1.2–1.4 keV [61]). Recently, a p-type point-contact Ge detector [34–36] has been implemented in the TEXONO experiment with reactor antineutrinos. The energy threshold of this detector is about 0.3 keV, whose value is comparable to the binding energies of the MI–III electrons in germanium (0.12–0.18 keV [61]). This means that an accurate analysis of the corresponding data will require numerical calculations based on the ab initio methods.
With lowering the T value down to T=2Eν2/(Eν+MN), an additional collision channel apart from ionization opens up, namely, the coherent elastic neutrino-nucleus scattering [72], which has not been observed experimentally so far. The early treatments of the atomic effects in the coherent elastic neutrino-nucleus scattering within SM can be found in [37, 73, 74]. It should be noted that any deviation of the measured cross section of the coherent elastic neutrino-nucleus scattering from the very precisely known SM value [75] will provide a signature of the BSM physics (see [76–79]). In this context, the accurate calculations of the NMM-induced contribution to the cross section of the coherent elastic neutrino-nucleus scattering appear to be of particular importance. The NMM-induced coherent neutrino scattering by single atoms as well as by crystals was discussed theoretically only in [80]. However, further studies are necessary for the correct interpretation of future measurements at low T values [81].
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was partially supported by the Russian Foundation for Basic Research (Grants no. 14-01-00420-a and no. 14-22-03043-ofi-m).
GiuntiC.KimC. W.2007Oxford, UKOxford University PressGonzalez-GarciaM. C.MaltoniM.Phenomenology with massive neutrinos20084601–3112910.1016/j.physrep.2007.12.0042-s2.0-42049088297BilenkyS.2010817New York, NY, USASpringerLecture Notes in PhysicsXingZ. Z.ZhouS.2011Zhejiang University PressGiuntiC.StudenikinA.Neutrino electromagnetic properties200972122089212510.1134/S10637788091201262-s2.0-74549116431BrogginiC.GiuntiC.StudenikinA.Electromagnetic properties of neutrinos201220124745952610.1155/2012/4595262-s2.0-84867357885GiuntiC.StudenikinA.Neutrino electromagnetic interactions: a window to new physicshttp://arxiv.org/abs/1403.6344NievesJ. F.Electromagnetic properties of Majorana neutrinos198226113152315810.1103/PhysRevD.26.31522-s2.0-0000244034KayserB.Majorana neutrinos and their electromagnetic properties19822671662167010.1103/PhysRevD.26.16622-s2.0-0001587472KayserB.CPT, CP, and C phases, and their effects, in Majorana-particle processes19843051023103310.1103/PhysRevD.30.10232-s2.0-26144467947PetcovS. T.The processes μ→e+γ,μ→e+e+e-,ν′→ν+γ in the Weinberg-Salam model with neutrino mixing197725340FujikawaK.ShrockR. E.Magnetic moment of a massive neutrino and neutrino-spin rotation1980451296396610.1103/PhysRevLett.45.9632-s2.0-0001049687PalP. B.WolfensteinL.Radiative decays of massive neutrinos198225376677310.1103/PhysRevD.25.7662-s2.0-0000093858ShrockR. E.Electromagnetic properties and decays of Dirac and Majorana neutrinos in a general class of gauge theories1982206335937910.1016/0550-3213(82)90273-52-s2.0-49049142567BeacomJ. F.VogelP.Neutrino magnetic moments, flavor mixing, and the super-kamiokande solar data199983255222522510.1103/PhysRevLett.83.52222-s2.0-0000403122BellN. F.CiriglianoV.Ramsey-MusolfM. J.VogelP.WiseM. B.How magnetic is the dirac neutrino?20059515415180210.1103/PhysRevLett.95.1518022-s2.0-28844508330BellN. F.GorchteinM.Ramsey-MusolfM. J.VogelP.WangP.Model independent bounds on magnetic moments of Majorana neutrinos2006642437738310.1016/j.physletb.2006.09.0552-s2.0-33750453031WongH. T.LiH. B.LinS. T.LeeF. S.SinghV.WuS. C.ChangC. Y.ChangH. M.ChenC. P.ChouM. H.DenizM.FangJ. M.HuC. H.HuangH. X.JonG. C.KuoW. S.LaiW. P.LeeS. C.LiJ.LiaoH. Y.LinF. K.LinS. K.LuJ. Q.ShengH. Y.SuR. F.TongW. S.XinB.YehT. R.YueQ.ZhouZ. Y.ZhuangB. A.Search of neutrino magnetic moments with a high-purity germanium detector at the Kuo-Sheng nuclear power station20077511601200110.1103/PhysRevD.75.0120012-s2.0-33846454005BedaA. G.BrudaninV. B.EgorovV. G.MedvedevD. V.PogosovV. S.ShirchenkoM. V.StarostinA. S.The results of search for the neutrino magnetic moment in GEMMA experiment201220121235015010.1155/2012/3501502-s2.0-84864957223AuerbachL. B.BurmanR. L.CaldwellD. O.ChurchE. D.DonahueJ. B.FazelyA.GarveyG. T.GunasinghaR. M.ImlayR.LouisW. C.MajkicR.MalikA.MetcalfW.MillsG. B.SandbergV.SmithD.StancuI.SungM.TayloeR.VandalenG. J.VemonW.WadiaN.WhiteD. H.YellinS.Measurement of electron-neutrino electron elastic scattering2001631111200110.1103/PhysRevD.63.1120012-s2.0-0034884641SchwienhorstR.CiampaD.EricksonC.GrahamM.HellerK.RusackR.SielaffJ.TrammellJ.WilcoxJ.KodamaK.UshidaN.AndreopoulosC.SaoulidouN.TzanakosG.YagerP.BallerB.BoehnleinD.FreemanW.LundbergB.MorfinJ.RameikaR.YunJ. C.SongJ. S.YoonC. S.ChungS. H.BerghausP.KubantsevM.ReayN. W.SidwellR.StantonN.YoshidaS.AokiS.HaraT.RheeJ. T.HoshinoK.JikoH.MiyanishiM.KomatsuM.NakamuraM.NakanoT.NiwaK.NonakaN.OkadaK.SatoO.AkdoganT.PaoloneV.RosenfeldC.KulikA.KafkaT.OliverW.PatzakT.SchnepsJ.A new upper limit for the tau-neutrino magnetic moment20015131-2232910.1016/S0370-2693(01)00746-82-s2.0-0035954639LiuD. W.AshieY.FukudaS.Limits on the neutrino magnetic moment using 1496 days of Super-Kamiokande-I solar neutrino data2004932502180210.1103/PhysRevLett.93.021802ArpesellaC.BackH. O.BalataM.Direct measurement of the 7Be solar neutrino flux with 192 days of Borexino data200810196091302RaffeltG. G.New bound on neutrino dipole moments from globular-cluster stars199064242856285810.1103/PhysRevLett.64.28562-s2.0-0000825649MohapatraR. N.PalP. B.2004723rdRiver Edge, NJ, USAWorld ScientificLecture Notes in PhysicsComminsE. D.BucksbaumP. H.19833rdCambridge, UKCambridge University PressCarlsonJ. F.OppenheimerJ. R.The impacts of fast electrons and magnetic neutrons193241676379210.1103/PhysRev.41.7632-s2.0-33645322687BetheH.Ionization power of a neutrino with magnetic moment1935311108115KyuldjievA. V.Searching for effects of neutrino magnetic moments at reactors and accelerators1984243338739710.1016/0550-3213(84)90482-62-s2.0-0000698314DomogatskyG.NadezhinD.Modern theory of star evolution and experiments of F. Reines on anti ν e-scattering detection197112678CowanC. L.ReinesF.Neutrino magnetic moment upper limit1957107252853010.1103/PhysRev.107.5282-s2.0-0040046566CowanC. L.Jr.ReinesF.HarrisonF. B.Upper limit on the neutrino magnetic moment1954965, article 129410.1103/PhysRev.96.12942-s2.0-36149024779VogelP.EngelJ.Neutrino electromagnetic form factors198939113378338310.1103/PhysRevD.39.33782-s2.0-30244462612WongH. T.Low energy neutrino and dark matter physics with sub-kev germanium detectors20112081463147010.1142/S02182718110196452-s2.0-79961150381LiH. B.LiaoH. Y.LinS. T.LiuS. K.SinghL.SinghM. K.SomaA. K.WongH. T.WuY. C.ZhaoW.AsryanG.ChuangY. C.DenizM.FangJ. M.HsuC. L.HuangT. R.Kiran KumarG.LeeS. C.LiJ.LiJ. M.LiY. J.LiY. L.LinC. W.LinF. K.LiuY. F.MaH.RuanX. C.ShenY. T.SinghV.TangC. J.TsengC. H.XuY.YangS. W.YuC. X.YueQ.ZengZ.ZeyrekM.ZhouZ. Y.Limits on spin-independent couplings of wimp dark matter with a p-type point-contact germanium detector201311026526130110.1103/PhysRevLett.110.2613012-s2.0-84879483944LiH. B.SinghL.SinghM. K.Differentiation of bulk and surface events in p-type point-contact germanium detectors for light WIMP searches20145618GaponovY. V.DobryninY. L.TikhonovV. N.Neutrino and antineutrino scattering on hydrogen-like atom1975222328334DobretsovV. Y.DobrotsvetovA. B.FayansS. A.Inelastic neutrino scattering by atomic electrons1992551180FayansS. A.DobretsovV. Y.DobrotsvetovA. B.Effect of atomic binding on inelastic ve scattering19922911-2162-s2.0-0002229128KopeikinV. I.MikaelyanL. A.SinevV. V.FayansS. A.Scattering of reactor antineutrinos by electrons19976011185918642-s2.0-0031316991FayansS. A.MikaelyanL. A.SinevV. V.Weak and magnetic inelastic scattering of antineutrinos on atomic electrons20016481475148010.1134/1.13989402-s2.0-0035537976GounarisG. J.PaschosE. A.PorfyriadisP. I.The ionization of H, He or Ne atoms using neutrinos or antineutrinos at keV energies20025251-2637010.1016/S0370-2693(01)01415-02-s2.0-0037122819KopeikinV. I.MikaelyanL. A.SinevV. V.Inelastic scattering of tritium-source antineutrinos on electrons of germanium atoms200366470771110.1134/1.15755702-s2.0-0038106144WongH. T.LiH.LinS.Enhanced sensitivities for the searches of neutrino magnetic moments through atomic ionization2010105606180110.1103/PhysRevLett.105.0618012-s2.0-77955314974VoloshinM.Neutrino scattering on atomic electrons in searches for the neutrino magnetic moment2010105204201801WongH. T.LiH.-B.LinS.-T.Enhanced sensitivities for the searches of neutrino magnetic moments through atomic ionization207410506180110.1103/PhysRevLett.105.061801KouzakovK. A.StudenikinA. I.Magnetic neutrino scattering on atomic electrons revisited2011696325225610.1016/j.physletb.2010.12.0432-s2.0-78651070755KouzakovK. A.StudenikinA. I.Electromagnetic neutrino-atom collisions: the role of electron binding2011217135335610.1016/j.nuclphysbps.2011.04.1372-s2.0-80051663068KouzakovK. A.StudenikinA. I.VoloshinM. B.Testing neutrino magnetic moment in ionization of atoms by neutrino impact2011931162362710.1134/S00213640111100752-s2.0-79961076101KouzakovK. A.StudenikinA. I.VoloshinM. B.Neutrino-impact ionization of atoms in searches for neutrino magnetic moment201183111111300110.1103/PhysRevD.83.1130012-s2.0-79960777340MartemyanovV. P.TsinoevV. G.Ionization of helium atoms under the effect of the antineutrino magnetic moment201174121671167510.1134/S10637788110701182-s2.0-84873813021KouzakovK. A.StudenikinA. I.On sensitivity of neutrino-helium ionizing collisions to neutrino magnetic moments2014114458461van HoveL.Correlations in space and time and born approximation scattering in systems of interacting particles19549524926210.1103/PhysRev.95.249MR00639642-s2.0-0007300482FanoU.Penetration of protons, alpha particles, and mesons196313166BelkicD.Bound-free non-relativistic transition form factors in atomic hydrogen19811412, article 0051907191410.1088/0022-3700/14/12/0052-s2.0-0001130552HoltA. R.Matrix elements for bound-free transitions in atomic hydrogen1969211, article 3111209121310.1088/0022-3700/2/11/3112-s2.0-2742593856AkhiezerA. I.BerestetskiiV. B.19652ndNew YorkWileyMR0191374LandauL. D.LifshitsE. M.19773rdOxford, UKPergamonKouzakovK. A.ZaytsevS. A.PopovY. V.TakahashiM.Singly ionizing 100-MeV/amu C6+He collisions with small momentum transfer2012863803271010.1103/PhysRevA.86.0327102-s2.0-84866705310ChenJ.-W.ChiH.-C.HuangK.-N.Atomic ionization of germanium by neutrinos from an ab initio approach2014731159162HenkeB. L.GulliksonE. M.DavisJ. C.X-ray interactions: photoabsorption, scattering, transmission, and reflection at E=50–30,000 eV, Z=1–92199354218134210.1006/adnd.1993.10132-s2.0-0004932883GrantI. P.Relativistic self-consistent fields1961262555576MR0129818GrantI. P.Relativistic self-consistent fields19658652352710.1088/0370-1328/86/3/311MR01877902-s2.0-11744294474MoruzziV. L.JanakJ. F.WilliamsA. R.1978Oxford, UKPergamonHuangK.-N.JohnsonW. R.Multiconfiguration relativistic random-phase approximation: theory198225263464910.1103/PhysRevA.25.6342-s2.0-5344232294HuangK.-N.Relativistic many-body theory of atomic transitions. The relativistic equation-of-motion approach198226273473910.1103/PhysRevA.26.7342-s2.0-35949017791JorgensenP.Molecular and atomic applications of time-dependent Hartree-Fock theory19752635938010.1146/annurev.pc.26.100175.0020432-s2.0-0016655864AmusiaM. Ya.CherepkovN. A.Many-body electron correlations in scattering processes1975547179LinD. L.Gauge properties of the Hartree-Fock and random-phase approximations197716260060510.1103/PhysRevA.16.6002-s2.0-25044476944HuangK.-N.ChiH.-C.ChouH.-S.The MCRRPA theory and its applications to photoexcitation and photoionization1995336565644DesclauxJ. P.A multiconfiguration relativistic DIRAC-FOCK program197591314510.1016/0010-4655(75)90054-52-s2.0-0002010637FreedmanD. Z.Coherent effects of a weak neutral current1974951389139210.1103/PhysRevD.9.13892-s2.0-35949032073GaponovYu. V.TikhonovV. N.Elastic scattering of low-energy neutrinos by atomic systems197726314SehgalL. M.WanningerM.Atomic effects in coherent neutrino scattering1986171110711210.1016/0370-2693(86)91008-72-s2.0-33645513496DrukierA.StodolskyL.Principles and applications of a neutral-current detector for neutrino physics and astronomy198430112295230910.1103/PhysRevD.30.22952-s2.0-33750966650ScholbergK.Prospects for measuring coherent neutrino-nucleus elastic scattering at a stopped-pion neutrino source2006733903300510.1103/PhysRevD.73.0330052-s2.0-33244476534BarrancoJ.MirandaO. G.RashbaT. I.Probing new physics with coherent neutrino scattering off nuclei20052005, article 0211410.1088/1126-6708/2005/12/0212-s2.0-33645528368BarrancoJ.MirandaO. G.RashbaT. I.Sensitivity of low energy neutrino experiments to physics beyond the standard model200776907300810.1103/PhysRevD.76.0730082-s2.0-35448975989DavidsonS.Peña-GarayC.RiusN.SantamariaA.Present and future bounds on non-standard neutrino interactions200333, article 0112-s2.0-23044473363AugustinJ.MüllerB.GreinerW.Magnetic neutrino scattering by crystals19904151683168610.1103/PhysRevD.41.16832-s2.0-35949008574FormaggioJ. A.Figueroa-FelicianoE.AndersonA. J.Sterile neutrinos, coherent scattering, and oscillometry measurements with low-temperature bolometers201285101300910.1103/PhysRevD.85.0130092-s2.0-84856655850