Neutral-current neutrino-nucleus scattering off Xe isotopes

Large liquid xenon detectors aiming for dark matter direct detection will soon become viable tools also for investigating neutrino physics. Information on the effects of nuclear structure in neutrino-nucleus scattering can be important in distinguishing neutrino backgrounds in such detectors. We perform calculations for differential and total cross sections of neutral-current neutrino scattering off the most abundant xenon isotopes. The nuclear structure calculations are made in the nuclear shell model for elastic scattering, and also in the quasiparticle random-phase approximation (QRPA) and microscopic quasiparticle phonon model (MQPM) for both elastic and inelastic scattering. Using suitable neutrino energy distributions, we compute estimates of total averaged cross sections for 8B solar neutrinos and supernova neutrinos.


Introduction
When the idea of neutrinos was first suggested by Pauli in 1930, it was thought that they would never be observed experimentally. Only two decades later interaction of neutrinos with matter was detected in the famous Cowan-Reines experiment [1]. More recently, detection and research of neutrinos has become more and more of an everyday commodity, and various more versatile ways to examine interactions of the little neutral one have emerged and are being tested in laboratories all over the world.
Coherent elastic neutrino-nucleus scattering (CEνNS) is a process where the neutrino interacts with the target nucleus as a whole instead of a single nucleon. Although CEνNS has been predicted since the 1970s [2], it was discovered only very recently by the COHERENT collaboration [3]. Due to the coherent enhancement, this experiment had the remarkable feature of detecting neutrinos with a compact 14.6 kg detector instead of a massive detector volume which is used in conventional neutrino experiments. Coherent neutrino nucleus scattering is on one hand an important potential source of information for beyond-standard-model physics [4,5,6,7,8,9,10,11], but on the other hand it may also hinder new discoveries as it will start disturbing dark matter detectors in the near future.
A great experimental effort has been put into directly detecting dark matter in the past few decades (see Ref. [12] for a review). The next-generation detectors are expected to be sensitive enough to probe cross sections low enough to start observing CEνNS as an irreducible background [13,14]. Solar neutrinos, atmospheric neutrinos, and diffuse supernova background neutrinos provide a natural source of background neutrinos, which for obvious reasons cannot be shielded against. As there are uncertainties in the fluxes of each of the aforementioned neutrino types, the sensitivity of WIMP detection is basically limited to the magnitude of this uncertainty. To make matters worse, it has been shown that for some specific WIMP masses and cross sections the recoil spectra of CEνNS very closely mimics that of scattering WIMPs [14].
It is therefore of utmost importance to device a way to go through this neutrino floor. One potential way of achieving this is having directional sensitivity in the detector [15,16]. As solar and atmospheric neutrinos have a distinct source within the solar system, it is expected that their recoil direction would be different of that of WIMPs, which are typically assumed to be gravitationally bound in a halo spanning the galaxy. Also arising from the different origin of neutrinos and WIMPs is the idea of using timing information to discriminate between neutrino and WIMP induced events in a detector [17]. Due to the motion of the Earth around the Sun, it is expected that the solar neutrino flux peaks around January, but the WIMP flux peaks in June when the velocities of the Sun and Earth are the most in phase. The recoil spectra of WIMPs and neutrinos could also be distinquished if the WIMP-nucleus interaction happens via a nonstandard operator emerging in the effective field theory framework [18,19].
Some of the leading dark matter experiments use a liquid xenon target [20,21,22,23,24], which allows for easy scalability to larger detector volumes. It is expected that the xenon detectors are the first to hit the neutrino floor. In this article we compute cross sections for elastic and inelastic neutrino-nucleus scattering for the most abundant xenon isotopes. For the coherent scattering we use the quasiparticle random-phase approximation (QRPA) framework and the nuclear shell model to model the nuclear structure and we compare the results between the two models. The wave functions of the states of odd-mass xenon isotopes are obtained by using the microscopic quasiparticle-phonon model (MQPM) on top of a QRPA calculation. Inelastic scattering is computed in the QRPA/MQPM formalism. In our calculations we consider 8 B solar neutrinos and supernova neutrinos.
A similar QRPA calculation has been made in Ref. [25] for 136 Xe, where both charged-current and neutral-current inelastic scattering was examined. Similar computations of neutral-current neutrino-nucleus scattering cross sections have been made before for the stable cadmium isotopes in Ref. [26], and for molybdenum isotopes in Ref. [27]. Both calculations used the QRPA/MQPM approach. To our knowledge this article presents the first calculation of neutral-current neutrinonucleus scattering within a complete microscopic nuclear framework for Xe isotopes other than 136 Xe.
This article is organized as follows. In section 2 we outline the formalism used to compute neutral-current neutrino-nucleus scattering. In section 3 we summarize the nuclear-structure calculations made for the target xenon isotopes. In section 4 we discuss the results of our crosssection calculations and in section 5 conclusions are drawn.

Neutral-current neutrino-nucleus scattering
In this section we summarise the formalism used to compute neutral-current neutrino-nucleus scattering processes. We examine standard-model reactions mediated by the neutral Z 0 boson, namely the processes i.e. the elastic and inelastic scattering of neutrinos off a nucleus (with A nucleons and Z protons), respectively. In the elastic process the initial and final states of the target nucleus are the same, while in the inelastic process excitation of the target nucleus takes place. The kinematics of the scattering process is illustrated in Fig. 1. We label the four-momenta of the incoming and ν ν ′ outgoing neutrino as k µ and k ′ µ . The momenta of the target nucleus before and after interacting with the neutrino are p µ and p ′ µ . The momentum transfer to the nucleus is referred to as The neutrino kinetic energy before and after scattering is E k and E k ′ The neutral-current neutrino-nucleus scattering differential cross section to an excited state of energy E ex can be written as [28] which comprises of the Coulomb-longitudinal (σ J CL ) and transverse (σ J T ) parts. They are defined as and q is the magnitude of the three-momentum transfer. The formalism and various different operators involved are discussed in detail in Refs. [28,29].
To compute the averaged cross section σ , we need to fold the computed cross sections with the energy distribution of the incoming neutrinos. We take the supernova neutrino spectrum to be of a two-parameter Fermi-Dirac character where α ν is the so-called pinching parameter, and T ν the neutrino temperature. The normalization factor F 2 (α ν ) is defined by the formula and the temperature and mean energy of neutrinos are related by We also examine solar neutrinos from 8 B beta decay. We use an 8 B neutrino energy spectrum from Ref. [30].

Nuclear structure of the target nuclei
In this section we outline the nuclear-structure calculations performed for the investigated nuclei 128,129,130,131,132,134,136 Xe. We have performed computations in the quasiparticle random-phase approximation (QRPA), microscopic quasiparticle-phonon model (MQPM), and the nuclear shell model.

QRPA/MQPM calculations
The nuclear structure of even-even Xe isotopes was computed by using the charge-conserving QRPA framework. The QRPA is based on a BCS calculation, where quasiparticle creation and annihilation operators are defined via the Bogoliubov-Valatin transformation as with the regular particle creation and annihilation operators c † α andc defined in [31]. Here α contains the quantum numbers (a, m α ) with a = (n a , l a , j a ). The excited states with respect to the QRPA vacuum are created with the phonon creation operator for an excited state ω = (J ω , M ω , π ω , k ω ), where k ω is a number labeling the excited states of given J π . In the above equation and X ω ab and Y ω ab are amplitudes describing the wave function that are solved from the QRPA equation where the matrix A is the basic Tamm-Dankoff matrix and B is the so called correlation matrix, both defined in detail in [31]. We perform the QRPA calculations using large model spaces consisting of the entire 0s-0d, 1p-0f -0g, 2s-1d-0h, and 1f -2p major shells, adding also the 0i 13/2 and 0i 11/2 orbitals. The single-particle bases are constructed by solving the Schrödinger equation for a Coulomb-corrected Woods-Saxon potential. We use the Woods-Saxon parameters given in Ref. [32]. We make an exception for 136 Xe, adopting the set of adjusted values of single-particle energies from Ref. [25]. Due to the neutron-magic nature of 136 Xe, adjusted single-particle energies are necessary to get agreement with experimental energy levels. The Bonn one-boson exchange potential [33] was used to estimate the residual two-body interaction.
The QRPA formalism involves several parameters that have to be fixed by fitting observables to experimental data. In the BCS calculation we fit the proton and neutron pairing strengths A p pair and A n pair so that the lowest quasiparticle energy matches the empirical pairing gap given by the three-point formula [34]: It should be noted, that for the neutron-magic 136 Xe this procedure cannot be done for the neutron pairing strength. We have instead used a bare value of A n pair = 1.0 for 136 Xe. The particle-particle and particle-hole terms of the two-body matrix elements are scaled by strength parameters G pp and G ph , respectively. The energies of the computed QRPA states are quite sensitive to these model parameters. We fit the lowest excited states of each J π separately to experimental values from Ref. [35] by altering the values of G pp and G ph . The values used for the model parameters are given in Table 1. Table 1: Model parameters used in the BCS and QRPA calculations. For each nucleus (column 1) the values of G pp and G ph (column 2) are given for the important J π phonons in columns 3 to 9. The QRPA process is known to produce states that are spurious, namely the first excited 0 + state and the first 1 − state. The first 0 + state has been deemed spurious in Refs. [26,36]. The first 1 − state is due to spurious center-of-mass motion as described in Ref. [31]. We have fitted the energies of these states to zero, if possible, by using the model parameters G pp and G ph , and subsequently the states have been omitted from calculations for the even-mass isotopes and also from the MQPM calculations for the odd-mass isotopes. The contributions of these spurious states to the total neutrino-nucleus scattering cross section would be tiny in any case.
Odd-mass xenon isotopes 129,131 Xe were then computed by using the MQPM formalism, in which we use a combination of one-and three-quasiparticle states by coupling a quasiparticle with a QRPA phonon to form the three-quasiparticle configurations. The MQPM basic excitation can be written in terms of quasiparticle and QRPA-phonon creation operators as [37] The amplitudes C and D are computed by solving the MQPM equations of motion. The detailed description of the process can be found in Ref. [37]. No additional model parameters are required for the MQPM calculation aside for the parameters fitted for the BCS/QRPA calculation described above. We do the MQPM calculations of 129 Xe and 131 Xe using 130 Xe and 132 Xe as reference nuclei, respectively. We select all QRPA phonons of J ≤ 6 with an energy less than 10 MeV to be used in the calculation.

Shell-model calculations
We perform shell model calculations for Xe isotopes using the shell-model code NuShellX@MSU [38]. We use the 0g 7/2 , 1d 5/2 , 1d 3/2 , 2s 1/2 , and 0h 11/2 valence space and the SN100PN interaction [39].  Table 2. For the isotopes 128−131 Xe we assume a completely filled 0g 7/2 orbital and for 128,129,131 Xe we also assume the 1d 5/2 orbital to be full. These should be reasonable approximations when aiming to describe the ground state and low-lying excited states in the xenon nuclei. The orbitals 0g 7/2 and 1d 5/2 have the lowest single-particle energies and the excitations are likely to take place from higher orbitals when the neutron number of the nuclei is quite large. Table 2: The valence-space truncations made in the shell-model calculations. The first column labels the Xe isotope, the following five columns give the minimum/maximum number of neutrons on the single-particle orbitals 0g 7/2 , 1d 5/2 , 1d 3/2 , 2s 1/2 , and 1h 11/2 , respectively. The computed energy levels of the even-mass xenon isotopes are given in Fig. 2 and the oddmass isotopes in Fig. 3. For the even-mass isotopes the experimental energy spectra are very well reproduced by the shell-model calculations. The accuracy is somewhat diminished when moving to lower masses from the closed neutron major shell of 136 Xe, but a decent correspondence between experimental and theoretical levels can be found. For the odd-mass isotopes the situation is more complex, but the positive-parity states are well reproduced by the calculations.
However, the negative-parity states 11/2 − and 9/2 − are computed to be much lower than in the experimental spectrum. This effect has been observed in earlier calculations using the SN100PN interaction in this mass region [40]. The experimental data for the xenon isotopes was obtained from [35]. To give a further measure of accuracy of our calculation, we computed the ground state magnetic moments for 129 Xe and 131 Xe. For 129 Xe the experimental magnetic moment of the 1/2 + ground state is µ exp = −0.7779763(84) µ N while the shell-model calculated value is µ sm = −1.360 µ N . For 131 Xe 3/2 + ground state the numbers are µ exp = +0.691862(4) µ N and µ sm = +1.059 µ N for experiment and shell model respectively. The sign of the magnetic moment in both cases is correct, but the magnitude of both of our calculated values is somewhat larger than that of the experimental ones.

Neutrino scattering results
In this section we present the results of our calculations for neutrino-nucleus scattering cross sections by methods described in Section 2. We have computed total cross sections for coherent and inelastic neutrino-nucleus scattering as a function of the neutrino energy, and also averaged total cross sections for solar 8 B neutrinos and supernova neutrinos scattering off the most abundant xenon isotopes. In the following calculations of averaged supernova-neutrino cross sections we have used two different neutrino temperatures corresponding to different neutrino flavors. We follow the choices of Refs. [26,36] and have the electron neutrinos described by parameters α = 3.0, E ν = 11.5 MeV, and T ν = 2.88 MeV, and the muon and tau neutrinos by α = 3.0, E ν = 16.3MeV, and T ν = 4.08MeV. Whenever we refer to supernova neutrinos in the following text these parameter values are used in the calculations.

Coherent elastic scattering
In Table 3 we present the total cross section for coherent neutrino-nucleus scattering off the target xenon isotopes as a function of neutrino energy. In Table 3 we only show calculations in the nuclear shell model, but the values for the QRPA/MQPM formalism are very similar, which is reflected on the total averaged cross sections shown later. The cross sections rise rapidly for small neutrino energies and start to saturate when approaching 100 MeV. The cross sections are larger for the higher-A isotopes, following the N 2 coherent enhancement. We present the total averaged cross section for supernova neutrinos as well as solar 8 B neutrinos in Table 4. Results for coherent scattering are shown for the shell model and QRPA/MQPM calculations. The results between the shell model and quasiparticle approaches are very similar. Some small differences can be observed in the results for the odd-mass isotopes, but those are still not very significant. The cross sections for the supernova neutrinos are larger than for 8 B neutrinos by roughly a factor of 3 or 5 depending on the neutrino flavor. This is due to the average energy of the supernova neutrinos being larger at 11.5 MeV or 16.3 MeV, while the 8 B spectrum peaks at around 7 MeV.

Inelastic scattering
Due to the limitations of the shell model in describing high-lying excited states, we compute inelastic scattering properties using only the QRPA/MQPM formalism, which is known to well depict the collective properties of excited nuclear states. The total cross section as a function of neutrino energy is given in Table 5 for each xenon isotope. For smaller neutrino energies, 0 to 30 MeV, the cross sections are slightly larger for the odd-mass isotopes than for their neighboring isotopes. The energies of solar neutrinos fit completely into this range, which leads to the averaged cross sections for solar neutrinos to be larger for the odd-mass isotopes.
The total averaged inelastic cross sections are listed in Table 4. The inelastic scattering cross sections are some orders of magnitude smaller than the coherent cross sections, as expected.
Here the cross sections of the supernova neutrinos are an order of magnitude or two larger than of 8 B solar neutrinos, again due to the supernova neutrinos having on average a higher energy. The effect of neutrino energy appears more pronounced in inelastic scattering than in coherent  scattering, however. The cross sections of the odd-mass isotopes are again slightly larger than those of the neighboring isotopes. We can compare our inelastic scattering results with those calculated in Ref. [26] for Cd isotopes using the same supernova neutrino parameters. The results for Cd isotopes in Ref. [26] in the case of electron neutrino range from 4.38 × 10 −42 cm 2 for 106 Cd to 4.96 × 10 −42 cm 2 for 111 Cd, with a general decreasing trend with increasing mass number for even-mass nuclei. Our results for Xe isotopes in Table 4 are very similar in magnitude, but the trend is rather rising than decreasing with increasing mass number. This could be a shell effect, as adding neutrons to Cd isotopes takes the nucleus further away from a closed major shell, but for the xenon nuclei it gets closer to a shell closure. Same conclusions can be made for the other neutrino flavors.    We show the contributions from different multipole channels to the total averaged cross sections in Fig. 4 for supernova electron neutrinos and Fig. 5 for solar neutrinos. It is evident that the most dominant contribution comes from an axial-vector 1 + multipole transition in all cases but one. Smaller, yet still important contributions arise from the axial-vector 1 − and 2 − channels for higher neutrino energies. This is characteristic behavior for neutral-current scattering, which has been observed in Ref. [26] for Cd isotopes, and in Ref. [27] for Mo isotopes. The contributions get more evenly distributed among the different multipoles with increasing  For the odd-mass nuclei our calculations also show a significant contribution from a vector 0 + channel, and for solar neutrinos scattering off 129 Xe this channel in fact becomes the strongest. For the even-mass isotopes this channel is more suppressed, but it becomes more significant for the lower energy solar neutrinos. Similar large 0 + contributions were observed in [26] for Cd isotopes. This is problematic as, in principle, the 0 + contribution is expected to be small because it vanishes at the limit q → 0. The particle-number violation of the quasiparticle framework can be an explanation for the large computed 0 + contribution. A detailed examination on the origins of the 0 + anomaly will be conducted in a later study. At this time one should regard the 0 + contributions with caution as they are probably at least partially spurious.
In Figs. 6 and 7 we show the dominating contributions to the inelastic scattering cross section from various final states of 128 Xe and 131 Xe, respectively. We notice that the major contributions are very similar for the solar and supernova electron neutrinos for the even-mass 128 Xe, where the leading contributions come from 1 + states at 8.4 MeV, 5.0 MeV, and 6.7 MeV. For solar neutrinos there is also a notable contribution from a 0 + state at 2.4 MeV. The situation is very much different for the odd-mass 131 Xe, where for supernova neutrinos there is a pile of 5/2 + , 3/2 + , and 1/2 + states at roughly 8 MeV giving large contributions to the total cross section in addition to the large contributions from lower-lying 5/2 + and 3/2 + states. However, for solar 8 B neutrinos this peak at 8 MeV is much smaller, and the leading contributions are more localized to the 5/2 + state at 1.8 MeV and the 3/2 + state at 2.9 MeV. It is interesting that a relatively small change in the average neutrino energy can lead to the higher-lying states to give much larger contributions to the total cross section.
Following the discussion on the anomalously large 0 + multipole contribution in 129 Xe we show the dominant final states for neutrinos scattering off 129 Xe in Fig. 8. As expected from the large 0 + multipole, the largest contributions here come from 1/2 + states at energies of roughly 2 − 3 MeV. Something in the nuclear-structure calculation seems to favor the 0 + multipole transition to 1/2 + final states over the 1 + multipole transition to 3/2 + states. Otherwise similar conclusions can be made for 129 Xe as for 131 Xe above about the location of the peaks in energy and differences between solar and supernova neutrinos.

Conclusions
We have computed various properties of cross sections of neutral-current neutrino-nucleus scattering off the most abundant Xe isotopes. The nuclear structure of our target Xe nuclei was computed in the nuclear shell model for elastic scattering, and in the QRPA framework for both elastic and inelastic scattering. For the odd-mass nuclei 129 Xe and 131 Xe an MQPM calculation was performed based on the QRPA calculation for 130 Xe and 132 Xe, respectively. We used realistic neutrino energy distributions for solar neutrinos from 8 B beta decay and supernova neutrinos to compute the averaged cross sections for each neutrino scenario. The total averaged cross section for supernova neutrinos are dependent of the values of the parameters α ν and E ν . We have shown results of only one set of parameters for electron neutrinos and one for muon/tau neutrinos. The dependence of the cross sections on the parameter α ν is typically quite mild, unless the change is large [25,26]. The values α ν = 3.0, E νe = 11.5MeV, and E νx = 16.3MeV used in this work are reasonable estimates and allow comparison of results with the 8 B solar neutrinos, for which the energy distribution is better known. A mapping of cross sections for various supernova neutrino parameters is out of scope of this work. However, we have tabulated total cross section as a function of neutrino energy, which can be used to obtain estimates of total averaged cross sections for any neutrino energy profile.
The scattering process in even-even nuclei is dominated by transitions to high-lying 1 + states, and for odd-mass nuclei typically by states differing from the initial state by one unit of angular momentum. We found that in even-mass nuclei the leading contributions from various final states are quite similar between solar neutrinos and supernova neutrinos. In odd-mass nuclei, however, the smaller energy of the solar neutrinos does not allow large contributions to the total cross section to arise from high-lying states. We also noted that the smaller energies of solar neutrinos lead into an enhancement in the vector 0 + multipole channel in comparison to the otherwise dominating 1 + axial-vector channel, especially in the odd-mass Xe nuclei. However, the large contribution from the 0 + multipole can be mostly spurious, possibly due to the particle number violation of the quasiparticle framework. This matter will be investigated further and subsequently reported elsewhere.