Next Article in Journal
Study of a Minimally Deformed Anisotropic Solution for Compact Objects with Massive Scalar Field in Brans–Dicke Gravity Admitting the Karmarkar Condition
Next Article in Special Issue
Weak Neutrino (Antineutrino) Charged-Current Responses and Scaling for Nuclear Matter in the Relativistic Mean Field
Previous Article in Journal
Probing the Neutron Skin of Unstable Nuclei with Heavy-Ion Collisions
Previous Article in Special Issue
High-Energy Lepton Scattering and Nuclear Structure Issues
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cross Sections for Coherent Elastic and Inelastic Neutrino-Nucleus Scattering

1
Department of Physics and Astronomy, Ghent University, Proeftuinstraat 86, B-9000 Gent, Belgium
2
Fermi National Accelerator Laboratory, Batavia, IL 60510, USA
3
Department of Physics, University of Florida, Gainesville, FL 32611, USA
*
Author to whom correspondence should be addressed.
Universe 2023, 9(5), 207; https://doi.org/10.3390/universe9050207
Submission received: 27 March 2023 / Revised: 19 April 2023 / Accepted: 19 April 2023 / Published: 25 April 2023
(This article belongs to the Special Issue Many Body Theory)

Abstract

:
The prospects of extracting new physics signals in coherent elastic neutrino–nucleus scattering (CE ν NS) processes are limited by the precision with which the underlying nuclear structure physics, embedded in the weak nuclear form factor, is known. We present calculations of charge and weak nuclear form factors and CE ν NS cross sections on 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb nuclei. We obtain the proton and neutron densities, and charge and weak form factors by solving Hartree–Fock (HF) equations with a Skyrme (SkE2) nuclear potential. We validate our approach by comparing 208 Pb and 40 Ar charge form factor predictions with available elastic electron scattering data. Since CE ν NS experiments at stopped-pion sources are also well suited to measure inelastic charged–current and neutral–current neutrino–nucleus cross sections, we also present calculations for these processes, incorporating a continuum Random Phase Approximation (CRPA) description on top of the HF–SkE2 picture of the nucleus. Providing both coherent as well as inelastic cross sections in a consistent framework, we aim at obtaining a reliable and detailed comparison of the strength of these processes in the energy region below 100 MeV. Furthermore, we attempt to gauge the level of theoretical uncertainty pertaining to the description of the 40 Ar form factor and CE ν NS cross sections by comparing relative differences between recent microscopic nuclear theory and widely-used phenomenological form factor predictions. Future precision measurements of CE ν NS will potentially help in constraining these nuclear structure details that will in turn improve prospects of extracting new physics.

1. Introduction

Coherent elastic neutrino–nucleus scattering (CE ν NS), where the only detectable reaction product is a low momentum recoiling nucleus, was suggested soon after the experimental discovery of a weak neutral current in neutrino interactions [1]. Even though for neutrino energies of some tens of MeV the CE ν NS cross section is a few orders of magnitude larger than competing inelastic processes, the difficulty in detecting the ∼keV scale recoil of a nucleus has hindered experimental detection of this process for decades. In 2017, the COHERENT collaboration detected the first CE ν NS signal using a stopped-pion beam in the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory with a CsI detector [2,3], followed up by another recent measurements in a liquid argon (LAr) detector [4,5] and in CsI detector [6].
The detection of CE ν NS has opened up a slew of opportunities in high–energy physics, astrophysics and in nuclear physics, inspiring new probes into beyond–Standard–Model (BSM) physics and new experimental methods. Several extensions of the SM that can be explored at low energy such as non–standard interactions (NSI) [7,8,9,10], sterile neutrinos [11,12], CP–violation [13], as well as exploration of nuclear effects, are being studied [14,15,16,17]. Several experimental programs have been or are being set up to detect CE ν NS and BSM signals in the near future using stopped-pion neutrino sources in COHERENT at the SNS [2], Coherent CAPTAIN–Mills (CCM) at Los Alamos National Laboratory (LANL) [18] and at the proposed European Spallation Source (ESS) facility [19], as well as reactor-produced neutrinos in CONNIE [20], MINER [21], ν GEN [22], NUCLEUS [23], RICOCHET [24], TEXONO [25], NEON [26] and vIOLETA [27].
The main source of uncertainty in the evaluation of the CE ν NS cross section is the accuracy with which the underlying nuclear structure and nucleon dynamics that determine the distributions of the nucleon density in the nuclear ground state, embedded in the form factor, are known in the target nucleus. The ground state proton (charge) density distributions are relatively well constrained through elastic electron scattering experiments pioneered by Hofstadter and collaborators at the Stanford Linear Accelerator [28], followed by other measurements in the following decades [29,30,31]. CE ν NS is however primarily sensitive to the neutron density distributions of the nucleus, which are only poorly constrained. Hadronic probes have been used to extract neutron distributions, these measurements are however plagued by ill-controlled model-dependent uncertainties associated with the strong interaction [32]. More (experimentally) challenging electroweak probes such as parity-violating electron scattering (PVES) [32,33] and CE ν NS provide relatively model-independent ways of determining neutron distributions. In recent years, one such PVES experiment, PREX at Jefferson lab, has measured the weak charge of 208 Pb at a single value of momentum transfer [34,35], while a follow up PREX–II experiment is ongoing to improve the precision of that measurement. Another PVES experiment, CREX at Jefferson lab, is underway to measure the weak form factor of 48 Ca [36]. Future ton and multi-ton CE ν NS detectors will enable more precise measurements and will potentially offer a powerful avenue to constrain neutron density distributions and weak form factors of nuclei at low momentum transfers where the process remains coherent [14,15,37].
As long as no precision measurements of neutron density distributions of nuclei are available, the weak nuclear form factor has to be modeled in order to evaluate the CE ν NS cross section and event rates. The accuracy of such an assumption is vital to the CE ν NS program since any experimentally measured deviation from the expected CE ν NS event rate can point to new physics or to unconstrained nuclear physics. It is therefore crucial to treat the underlying nuclear structure physics that is embedded in nuclear form factors with utmost care. Phenomenological approaches, such as the Klein–Nystrand form factor [38] adapted by the COHERENT collaboration, or the Helm form factor [39] where density distributions are represented by analytical expressions, are widely used in the CE ν NS community. Empirical values of the proton rms radius, measured in elastic electron scattering, are often used to evaluate the proton form factor and often similar parameterizations are assumed for the neutron form factor. Microscopic nuclear physics approaches which provide a more accurate description of the nuclear ground state and density distributions such as density functional theory [37], coupled–cluster theory from first principles [40], relativistic mean–field model [41], Hartree–Fock plus Bardeen–Cooper–Schrieffer model [42] as well as effective field theory approaches [43,44] have also been reported in recent years.
In this work we will present a microscopic many-body nuclear theory model where the nuclear ground state is described in a Hartree–Fock (HF) approach with a Skyrme (SkE2) nuclear potential. We calculate proton and neutron density distributions, charge and weak form factors, and CE ν NS cross sections on 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb, and confront our predictions with the available experimental data. In view of the worldwide interest in liquid-argon-based neutrino and dark matter experiments, we pay special attention to the 40 Ar nucleus. We attempt to gauge the level of theoretical uncertainty pertaining to the description of the 40 Ar form factor and CE ν NS cross section by comparing relative differences between recent nuclear theory and widely-used phenomenological form factor predictions.
CE ν NS experiments at stopped-pion sources are also well-suited to measure inelastic neutrino–nucleus cross sections. These measurements, in particular on 40 Ar, will provide powerful constraints on supernova detection capabilities of future kiloton neutrino experiments. To this end, we also present inelastic charged–current (CC) and neutral–current (NC) cross section calculations on 40 Ar, incorporating a continuum Random Phase Approximation (CRPA) description on top of the initial HF–SkE2 picture of the nucleus.
The remainder of this manuscript is organized as follows. In Section 2, we lay out the general formalism of calculating the CE ν NS and inelastic neutrino–nucleus scattering cross section. In Section 3, we present results of proton and neutron densities, charge and weak form factors, and CE ν NS cross sections on 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb obtained within our HF–SkE2 approach. We focus on 40 Ar in Section 3, and compare our predictions with experimental data and other theoretical calculations. We also present inelastic cross sections on 40 Ar in Section 3. We present conclusions of this study in Section 4.

2. Formalism

In this section, we lay out the general formalism for calculating cross sections of the coherent elastic and inelastic neutrino-nucleus scattering process.

2.1. CE ν NS Cross Section

A neutrino with four momentum k i = ( E i , k i ) scatters off the nucleus, which is initially at rest in the lab frame with p A = ( M A , 0 ) , exchanging a Z 0 boson. The neutrino scatters off, carrying away four momentum k f = ( E f , k f ) while the nucleus remains in its ground state and receives a small recoil energy T, so that p A = ( M A + T , p A ) with | p A | = ( M A + T ) 2 M A 2 and T = q 2 / 2 M A . Here, M A is the rest mass of the nucleus, q = | q | is the absolute value of the three–momentum transfer which is of the order of keV for neutrino energies of tens of MeV, Q 2 q 2 = | k f k i | 2 , and the velocity dependent factor in the denominator refers to the relative velocity of the interacting particles. The process is schematically shown in Figure 1.
The initial elementary expression for the cross section reads
d 6 σ = 1 v i v A m i E i m f E f d 3 k f ( 2 π ) 3 M A M A + T d 3 p A ( 2 π ) 3 × ( 2 π ) 4 ¯ f i M 2 δ ( 4 ) ( k i + p A k f p A ) .
This expression can be integrated to yield the expression for the cross section differential in neutrino scattering angle θ f :
d σ d cos θ f = m i E i m f E f M A M A + T E f 2 2 π f r e c 1 ¯ f i M 2 .
The recoil factor reads
f r e c = E i E f M A M A + T .
Working out the Feynman amplitude one gets
¯ f i M 2 = G F 2 2 L μ ν W μ ν ,
with the nuclear tensor W μ ν reading
W μ ν = ¯ f i ( J n u c μ ) J n u c ν .
The summation symbols in these expressions denote summing and averaging over initial and final polarizations respectively. The nuclear tensor depends on the nuclear current transition amplitudes:
J n u c μ = Φ 0 | J ^ μ ( q ) | Φ 0 .
Under the assumption that the nuclei of interest are spherically symmetric with J π = 0 + and taking the z-axis to be along the direction of q , one only needs to take into account the zeroth and third component of the nuclear current’s vector part, which are furthermore connected through vector current conservation (CVC):
q μ J ^ μ ( q ) = 0 .
Through performing the necessary algebra, one arrives at the final expression
d σ d cos θ f = G F 2 2 π E f 3 E i Q 4 q 4 ( 1 + cos θ f ) | J 0 V | 2
where J 0 V is the transition amplitude induced by the nuclear current. One can then safely approximate Q 4 q 4 1 and express the differential cross section as a function of the neutrino scattering angle θ f as:
d σ d cos θ f = G F 2 2 π E f 3 E i ( 1 + cos θ f ) Q W 2 4 F W 2 ( Q 2 )
where G F is the Fermi coupling constant, and Q W the weak nuclear charge:
Q W 2 = [ g p V Z + g n V N ] 2 = [ ( 1 4 sin 2 θ W ) Z N ] 2
with coupling constants g n V = 1 and g p V = ( 1 4 sin 2 θ W ) . N and Z are the nucleus’ neutron and proton number, and θ W is the weak mixing angle. The value is such that sin 2 θ W = 0.23857 , which is valid at low momentum transfers [45].
Here we have introduced the elastic form factor, F W 2 ( Q 2 ) , which we will discuss later in this subsection. In elastic scattering the entire nuclear dynamics is encoded in this form factor. Equivalently one can express the differential cross section as a function of the nuclear recoil T, which reads:
d σ d T = G F 2 π M A 1 T E i M A T 2 E i 2 Q W 2 4 F W 2 ( Q 2 ) ,
In Equations (9) and (11), we have expressed the CE ν NS kinematic distribution both in neutrino scattering angle, θ f , and in nuclear recoil energy T. In most experiments the only signal of a CE ν NS event is a nuclear recoil energy deposition. In principle, future experiments with more advanced detector technologies may be able to detect both nuclear recoil and angular distribution simultaneously. Such capabilities are already being explored in some dark-matter experiments and will greatly enhance the physics capabilities of future CE ν NS experiments [46].
The scattering process’ cross section is proportional to the squared magnitude of the transition amplitude induced by the nuclear current. Since the relevant ground state to ground state transition for spherically symmetrical nuclei is 0 + 0 + , only the vector part of the current will contribute. The amplitude can be expressed as
J 0 V = Φ 0 | J ^ 0 V ( q ) | Φ 0 = e i q · r Φ 0 | J ^ 0 V ( r ) | Φ 0 = 1 2 1 4 sin 2 θ W f p ( q ) F p ( Q 2 ) f n ( q ) F n ( Q 2 ) ,
where we have inserted the impulse approximation (IA) expression for the nuclear current, as a sum of single-body operators:
J ^ 0 V ( r ) = i F Z ( Q 2 , i ) δ ( 3 ) ( r r i ) ,
with
F Z ( Q 2 , i ) = 1 2 sin 2 θ W ( F p F n ) τ 3 ( i ) sin 2 θ W ( F p + F n ) ,
where we used the convention τ 3 ( i ) = + 1 for proton, −1 for neutrons. Furthermore, f p ( q ) and f n ( q ) are the Fourier transforms of the proton and neutron densities, respectively. F p and F n are proton and neutron form factors, for which we adopt the standard Galster parametrization. Note that using a more sophisticated parametrization of the form factor, other than Galster, will not affect the results at the energies relevant to this work. The overall structure of the transition amplitude consists of products of the weak charge with two factors: the nuclear form factor, determined by the spatial distribution of the nucleons in the nucleus, as well as the nucleon form factor. We arrive at the expression:
F W ( Q 2 ) = 1 Q W 1 4 sin 2 θ W f p ( q ) F p ( Q 2 ) f n ( q ) F n ( Q 2 ) = 2 Q W J 0 V ,
such that the form factor becomes 1 in the static limit. Note that in writing down the functional dependence we can make use of the non-relativistic approximation Q | q | , valid in the energy regime considered.
We employ a microscopic many–body nuclear theory model where the nuclear ground state is described in a Hartree–Fock (HF) approach with a Skyrme (SkE2) nuclear potential, which we will refer to as HF–SkE2. We solve the HF equations to obtain single-nucleon wave functions for the bound nucleons in the nuclear ground state. We evaluate proton ( ρ p ( r ) ) and neutron ( ρ n ( r ) ) density distributions from those wave functions. The proton density is utilized to calculate charge form factor (which can also be referred to as electromagnetic form factor), F c h ( Q 2 ) , while both proton and neutron densities are utilized to compute weak, F W ( Q 2 ) , nuclear form factor, as shown in Equation (15). This approach involves more realistic nuclear structure calculations of proton and neutron density distributions making it more reliable compared to the phenomenological approaches that rely on the approximation ρ n ( r ) ρ p ( r ) , utilizing empirical values of ρ p ( r ) extracted from electron scattering experiments.

2.2. Inelastic Cross Sections

CE ν NS experiments at stopped-pion sources are also sensitive to inelastic neutrino–nucleus interactions. In several astrophysical environments elastic and inelastic processes come in competition. To this end, we also present calculations of inelastic charged–current (CC) and neutral–current (NC) cross sections, calculated within the same framework. These results are obtained by including effects of long-range correlations through a continuum Random Phase Approximation (CRPA) description on top of the HF–SkE2 initial picture of the nucleus.
The inelastic neutrino–nucleus scattering process is schematically shown in Figure 1. A neutrino with four momentum k i = ( E i , k i ) scatters off the nucleus, which is initially at rest in the lab frame, exchanging a W + (CC) or a Z 0 (NC) boson. The nucleus receives four momentum Q = ( ω , q ) , where ω = E i E f and q = k i k f , while the scattered lepton carries away four momentum k f = ( E f , k f ) . Since we concern ourselves with inclusive calculations, the hadronic part of the final states are integrated out. The inelastic neutrino–nucleus differential cross section of this process can be written as
d 3 σ d ω d Ω = σ W E f k f ζ 2 ( Z , E f ) × v C C R C C + v C L R C L + v L L R L L + v T R T + h v T R T ,
with the Mott-like cross section prefactor σ W defined as
σ W C C = G F cos θ c 2 π 2 , σ W N C = G F 2 π 2 ,
where G F is the Fermi constant and cos θ c the Cabibbo angle. The factor ζ 2 ( Z , E f ) is introduced in order to take into account the distortion of the scattered lepton wave function in the Coulomb field of the final nucleus with Z protons, in the case of CC interaction [47]. In the NC case ζ 2 ( Z , E f ) equals 1. The influence of the lepton helicity on the cross section is encoded in h which is + for neutrinos and − for antineutrinos.
The v-factors are leptonic functions that are entirely determined by lepton kinematics. The R-factors are the nuclear response functions that depend on the energy and momentum transfer ( ω , q) and contain all the nuclear information involved in this process. The indices L and T correspond to longitudinal and transverse contributions, relative to the direction of the momentum transfer. The nuclear responses are function of the transition amplitude, J μ n u c l ( ω , q ) , between the initial | Φ 0 and final | Φ f state:
J μ n u c l ( ω , q ) = Φ f | J ^ μ ( q ) | Φ 0 ,
where the nuclear current, J ^ μ ( q ) , is the Fourier transform of the nuclear current operator in coordinate space:
J ^ μ ( q ) = d x e i x · q J ^ μ ( x ) .
These are computed within a HF–CRPA framework. For a detailed discussion of the nuclear response we refer the reader to our previous work in refs. [47,48,49,50,51,52,53,54,55,56,57,58,59,60,61]. Here we briefly describe the essence of our approach. The CRPA description goes beyond a pure spectator approach, incorporating long-range correlations in the cross section calculations. Within many–body theory, the random phase approximation achieves this by modeling excitations as superpositions of particle–hole ( p h 1 ) and hole-particle ( h p 1 ) states out of a correlated ground state:
| Ψ R P A C = C X c , c | p h 1 Y c , c | h p 1 ,
where the summation index C denotes a set of quantum numbers defining an excitation channel unambiguously:
C = n h , l h , j h , m j h , ε h ; l p , j p , m j p , τ z .
The indices p and h represent the quantum numbers related to the particle or the hole state, ε h denotes the binding-energy of the hole state and τ z defines the isospin character of the particle-hole pair. Since the RPA approach describes nuclear excitations as the coherent superposition of individual particle-hole states out of a correlated ground state, it allows the description of collective effects in the nucleus.
Besides Equation (19), the RPA approach can also be formulated in a propagator description of many–body theory, where the central object containing the information on the excited states of the many–body system is the polarization propagator. In the Lehmann representation, the CRPA approach involves solving the RPA equation for the local polarization propagator Π R P A ( x 1 , x 2 , E e x c ) in coordinate space:
Π R P A ( x 1 , x 2 , E e x c ) = Π ( 0 ) ( x 1 , x 2 , E e x c ) + 1 d x d x Π ( 0 ) ( x 1 , x , E e x c ) × V ˜ ( x , x ) Π R P A ( x , x 2 , E e x c ) ,
where E e x c is the excitation energy of the target nucleus and x is the shorthand notation for the combination of the spatial, spin, and isospin coordinates. In this equation, the antisymmetrized residual interaction V ˜ ( x , x ) , is the same SkE2 Skyrme interaction we have utilized to calculate the single particle wave functions (and therefore, nuclear densities) of the CE ν NS cross sections, keeping the scheme self-consistent. Π ( 0 ) ( x 1 , x 2 , E e x c ) denotes the zeroth-order contribution to the polarization propagator which is equivalent to the HF contribution. The (local) polarization propagator Π R P A ( x 1 , x 2 , E e x c ) , which describes the propagation of particle–hole pairs, is obtained by adding the iteration of first-order contributions to the bare local polarization propagator Π ( 0 ) ( x 1 , x 2 , E e x c ) . By solving this equation, one obtains the CRPA transition amplitudes needed to calculate the inelastic neutrino–nucleus cross sections.
It is worth mentioning that the effect of long-range correlations included through the CRPA approach, vital for inelastic calculations at low energies, are found to be negligible in evaluating ground state densities of nuclei [62] and are therefore not included in the elastic scattering calculations discussed in Section 2.1.
The HF–CRPA framework offers an elegant formalism that accounts for collective excitations in the continuous spectrum as well as describes quasielastic neutrino–nucleus scattering in the low and medium energy regime. Our model has been developed over decades and has been utilized extensively to calculate various electron- and neutrino–nucleus cross sections suited for astrophysical processes as well as accelerator-based neutrino oscillation experiments [47,48,49,50,51,52,53,54,55,56,57,58,59,60,61].

3. Results and Discussion

Since the weak charge of the proton is strongly suppressed by the weak mixing angle (Equation (15)) the nuclear weak charge is predominately carried by the neutrons. The weak form factor F W ( Q 2 ) , and hence the CE ν NS cross section, are both dominated by the distribution of neutrons within the nucleus. As proton densities are well-constrained by experimental elastic electron scattering data [31] while little reliable neutron density data is available, phenomenological approaches approximate ρ n ( r ) ρ p ( r ) and thus assume F n ( Q 2 ) F p ( Q 2 ) , making the nuclear form factor more of a global factor [63]. Within the HF–SkE2 approach we treat proton and neutron densities and their corresponding form factors separately and do not have to rely on such assumptions. The densities are defined in terms of the reduced radial single particle wave functions as
ρ q ( r ) = 1 4 π r 2 a v a , q 2 ( 2 j a + 1 ) | ϕ a , q ( r ) | 2 ,
with v a , q 2 being the occupation probability of orbital a of nature q (i.e., proton p or neutron n). In Table 1, we show shells and single-particle energy levels, in the case of 40 Ar nucleus, as yielded by a HF calculation using SkE2 potential.
In Figure 2, we present proton (panel (a)) and neutron (panel (b)) density distributions of 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb obtained using our HF–SkE2 approach. Naturally the heavier the nucleus, the more broadly the densities are distributed. Panel (c) and (d) show the charge and weak form factors for all the nuclei. In both the charge and weak form factor cases, the heavier the nuclei the faster the form factor encounters its first minimum at rising q values. Lighter nuclei have their minima spread over a larger q range. 12 C has its first minimum at q 1.8 fm 1 while 208 Pb has its first minimum around q 0.65 fm 1 . Although the charge and weak form factors have a similar overall structure, the minima and maxima of both occur at slightly different values of the momentum transfer, with larger differences in heavier nuclei. To further illustrate this, in Figure 3, we show the “weak–skin" form factor [32] for all these nuclei, defined as the difference between the charge and weak form factors:
F W , skin ( q ) = F ch ( q ) F W ( q ) ,
which, near the origin, is proportional to the experimentally observable weak skin [32]. The figure illustrates that the charge and weak form factors significantly differ from each other.
In the left panel of Figure 4, we show our predictions for the charge form factor of 208 Pb. The predictions are compared with the experimental charge form factor obtained from a Fourier–Bessel fit to the elastic electron scattering data of Ref. [29]. Our predictions describe the experimental data remarkably well. Our predictions almost overlap with data for q 1.8 fm 1 . We also performed a comparison with the relativistic mean–field (RMF) predictions of Yang et al. [41]. There are no visible differences between both models up to q 1.8 fm 1 . The right panel shows our predictions for the weak form factor, again compared with the RMF predictions of [41]. We also show the single data point measured at a momentum transfer of q = 0.475 fm 1 by the PREX collaboration [34,35]. This remains the only measurement of the weak form factor obtained with an electroweak probe. The error bars on the data point are too large to discriminate between theoretical predictions. The follow-up PREX–II measurement at Jefferson lab aims to reduce the error bars by at least a factor of three.
The total CE ν NS cross section as a function of neutrino energy for 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb is shown in Figure 5. All nuclei show a similar behavior: there is a rapid rise of the cross section for incoming neutrino energies up to about ∼30 MeV, then the steep increase slows down and flattens out on the log scale thereafter. The cross section increases with the atomic number, with nearly two to three orders of magnitude difference between 12 C and 208 Pb, reflecting the N 2 scaling behavior shown in Equation (10).
To demonstrate the dominance of the CE ν NS strength over the quasielastic one for a neutrino energy of a few tens of MeV, in Figure 6 we compare CE ν NS cross sections to ν e –nucleus charged–current quasielastic (CCQE) and neutral–current quasielastic (NCQE) cross sections. For the energies relevant for pion decay–at–rest neutrinos, E 52 MeV, the CE ν NS cross section is roughly two orders of magnitude larger than inelastic cross sections.

Constraining 40 Ar

In view of the worldwide interest in liquid argon (LAr)-based detectors in neutrino and dark matter experiments, in this section we will focus on 40 Ar. In the COHERENT collaboration’s expanding series of detectors at SNS, the collaboration has recently presented new measurements from a 24 kg, single-phase, LAr CENNS–10 detector [5] while a ton-scale LAr experiment is underway. A 10 ton LAr scintillation detector, Coherent CAPTAIN-Mills (CCM), was recently built at LANL to study CE ν NS on 40 Ar and to search for low-mass dark matter that coherently scatters off 40 Ar nuclei [18]. Several other neutrino [64,65] and dark matter experiments [66,67,68,69] employ LAr detectors, making it vital to study ground state properties of the 40 Ar nucleus.
In Figure 7 (left) we compare our argon charge form factor ( F c h ( q ) ) predictions with the elastic electron scattering data of Ref. [70]. Our predictions describe experimental data remarkably well for q 2 fm 1 , validating our approach. We also compare with the predictions of Payne et al. [40] where form factors are calculated within a coupled–cluster approach, using a chiral NNLO sat interaction. At higher q, q 2 fm 1 , both predictions diverge from experimental data. Note that for neutrino energies relevant for pion decay-at-rest the region above q 0.5 fm 1 does not contribute to CE ν NS cross sections. We also show a comparison with two phenomenological form factors which are widely used in the CE ν NS community: the Klein–Nystrand [38] form factor that is adapted by the COHERENT collaboration and the Helm form factor [39]. Note that we also show an adapted version of the Klein–Nystrand form factor that will be described in more detail in our discussion of form factor predictions, later in this section.
After validating our approach, we make predictions for the weak form factor of 40 Ar in Figure 7 (right). There is no data available for the weak form factor on argon yet. We compare our results with the prediction of Payne et al. [40], Yang et al. [41] and Hoferichter et al. [43] as well as with the Helm form factor [39], the Klein–Nystrand [38] and an adapted version of the Klein–Nystrand form factor. Overall, the shape and structure of the weak form factor is similar to the charged one, but the positions of minima and maxima are somewhat different. In our HF–SKE2 approach the first minimum of F c h ( q ) is at q 1.23 fm 1 while for F W ( q ) it lies at q 1.19 fm 1 , pointing to the fact that the neutron distribution extends further out compared to the proton one. To quantify differences between the charge and weak form factor, in Figure 3 we show the “weak-skin” form factor of 40 Ar using Equation (23).
In order to appreciate which values of momentum transfer q are involved at different neutrino energies, as well as to see at which q values the differences in the nuclear modeling start causing discrepancies in reaction strength predictions, we plot cumulative cross sections for 40 Ar at several neutrino energies and for different models in Figure 8. This is defined as the total cross section strength, integrated up to a cutoff value q c u t o f f in the momentum transfer:
σ ( q c u t o f f ) = 0 q c u t o f f d σ ( q ) d q d q
The model differences become stronger for increasingly high energies with discrepancies originating from the higher-qq regions of the elastic form factor. The range of cutoff values also coincides with all kinematically available momentum transfers. At 100 MeV e.g., 40 Ar is only probed up to q 1 fm 1 .
In Figure 9, we show differential cross sections on 40 Ar as a function of recoil energy T, and scattering angle θ f , for different incoming neutrino energies according to Equations (9) and (11). For comparison, we have also plot the case with no nuclear structure effects i.e., F ( Q 2 ) = 1 . The effects of nuclear structure physics are more prominent as the neutrino energy increases. Most of the cross section strength lies at the lower-end of the recoil energy spectrum and for forward scattering as the cross section falls off rapidly at higher T (top panels) and higher θ f values (bottom panels). Most CE ν NS detectors are sensitive only to the recoil energy deposited in the detector but, in principle, in the future more advanced detector technologies might enable measurement of both nuclear recoil and angular distribution simultaneously. Utilizing such additional information can be valuable in disentangling new physics signals in CE ν NS experiments [46].
In Figure 7 (right), we come back to the differences between various predictions. Different form factor approaches are based on different representations of the nuclear densities, with no experimental data to constrain neutron distributions. Identifying the size of the differences between various theoretical predictions is crucial as experiments have to assign any deviation from expected event rates either to new physics or to unconstrained nuclear physics. We compare six predictions. These include four nuclear theory approaches: the HF–SkE2 calculation of this work, the predictions of Payne et al. [40], and the RMF calculations of Yang et al. [41] where form factors predictions are informed by properties of finite nuclei and neutron star matter, and the predictions of Hoferichter et al. [43] where form factors are calculated using a large-scale nuclear shell model. They also contain two phenomenological approaches: the Helm [39] and Klein–Nystrand [38] form factors where density distributions are represented by analytical expressions.
In the Helm approach [39] the density distribution is described as a convolution of a uniform nucleonic density with a given radius and a Gaussian profile characterized by the folding width s, accounting for the nuclear skin thickness. The resulting form factor is expressed as:
F Helm ( q 2 ) = 3 j 1 ( q R 0 ) q R 0 e q 2 s 2 / 2 ,
where j 1 ( x ) = sin ( x ) / x 2 cos ( x ) / x is a spherical Bessel function of the first kind. R 0 is an effective nuclear radius given as: R 0 2 = ( 1.23 A 1 / 3 0.6 ) 2 + 7 3 π 2 r 0 2 5 s 2 with r 0 = 0.52 fm and s = 0.9 fm, fitted [71,72] to muon spectroscopy and electron scattering data compiled in [30]. The Klein–Nystrand (KN) form factor, adapted by the COHERENT Collaboration, is obtained from the convolution of a short–range Yukawa potential with range a k = 0.7 fm over a Woods–Saxon distribution approximated as a hard sphere with radius R A = 1.23 A 1 / 3 fm [38]. The resulting form factor is expressed as:
F KN ( q 2 ) = 3 j 1 ( q R A ) q R A 1 1 + q 2 a k 2 .
An adapted version of the KN form factor is often used, where R A is defined as R A = 5 3 r 0 2 10 a k 2 utilizing measured proton rms radii r 0 of the nucleus [16,63]. We show both the standard and the adapted (ad.) KN form factor. For the adapted one we use r 0 = 3.427 fm, the measured proton rms radii of 40 Ar [31].
We attempt to quantify differences between different form factors and the CE ν NS cross section due to different underlying nuclear structure details. We consider quantities that emphasize the relative differences between the results of different calculations, arbitrarily using HF–SkE2 as a reference calculation, as follows:
| Δ F W i ( q ) | = | F W i ( q ) F W HF ( q ) | | F W HF ( q ) | ,
and
Δ σ W i ( E ) = | σ W i ( E ) σ W HF ( E ) | σ W HF ( E ) ,
where i refers to calculations from different approaches as discussed above.
The relative differences are shown in Figure 10. We show only the low-momentum part of the weak form factor to a maximum value of q = 0.5 fm 1 (∼100 MeV) that corresponds to a maximum incoming neutrino energy of E ∼ 50 MeV, as shown in Figure 8. The relative differences are shown on a linear scale. At smaller energies the momentum transfer is low and hence the differences between form factors are also small. For higher energies the available momentum transfer increases and therefore the differences between the form factors become more prevalent. The differences in model predictions amount to < 7.5 % over the entire momentum transfer range. The differences rise rapidly at the higher end of the q range. This translates into relative differences in CE ν NS cross sections, Δ σ ( E ) , of < 5 % over the whole energy range, where E 55 MeV, relevant for neutrinos from pion decay-at-rest. Note that most of the strength in the cross section lies at the lower T end (and therefore at the lower q end), as we have seen in Figure 9.
The CE ν NS cross section on 40 Ar as a function of the neutrino energy is shown in Figure 11 (left). We also show recent flux-averaged measurements performed by the COHERENT collaboration [5]. Measurements from two analyses are included, with the horizontal bars indicating the minimum value set by the nuclear recoil threshold energy for each analysis. The flux–averaged measured cross section is 2.2 ± 0.7 × 10 39 cm 2 (average of both analyses), while the HF–SkE2 predicted flux–averaged cross section is 1.82 × 10 39 cm 2 . The total experimental error is dominated by statistics, amounting to ∼ 30%. Future measurements by ton–scale LAr detector at SNS and 10-ton LAr detector CCM at LANL will be able to provide more precise measurements of the CEνNS cross section on 40 Ar. In Figure 11 (right), we also show flux–folded cross sections as a function of neutron number for all five nuclei— 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb—considered in this paper. As expected, the deviation of F ( Q 2 ) = 1 from the full HF–SkE2 calculation becomes more prominent as the number of neutrons, and hence the influence of nuclear structure effects, increases. Also included is the 40 Ar data measured by COHERENT [5].
CE ν NS liquid argon detectors at stopped-pion sources are well suited to measure inelastic cross sections as well. Inelastic cross section measurements on 40 Ar will provide powerful constraints on supernova detection capabilities of future kiloton liquid argon experiments such as DUNE [65]. In view of this, in Figure 12 we present CC (left) inelastic ( ν e , 40 Ar) and NC (right) inelastic ( ν , 40 Ar) cross sections for energies relevant to pion decay–at–rest neutrinos. These cross sections are calculated by incorporating the CRPA approach on top of the initial HF–SkE2 nuclear picture.
The top panels in Figure 12 show total cross section as a function of incoming neutrino energy along with separate contributions coming from the dominating individual multipoles. In both CC and NC case, most strength arises from 1 , 1 + and 2 multipoles. The 0 + and 0 transitions contribute only minimally to the total reaction strength for excitations into the continuum and are not shown here. Still, it is clear that a considerable part of the strength stems from forbidden transitions. The middle panels show the differential cross sections as a function of excitation energy ω for two incoming neutrino energies E ν = 30 MeV and 50 MeV. As the energy increases, more resonance peaks show up as an increasing number of excitations becomes accessible. Differential cross sections are folded with a Lorentzian of width 3 MeV in order to account for the finite width of the resonances [56]. The bottom panels show the differential cross sections as a function of the direction of the outgoing lepton scattering angle cos θ f for two incoming neutrino energies E ν = 30 MeV and 50 MeV. The differential cross sections in scattering angles favor backward scattering.

4. Conclusions

The experimental observation of coherent elastic neutrino–nucleus scattering processes by the COHERENT collaboration has inspired physicists across many fields. The power of CE ν NS as a probe of BSM physics and its potential for determining neutron density distributions is becoming more and more apparent. The main uncertainty in the evaluation of the CE ν NS cross sections is driven by the weak form factor that encodes the entire nuclear structure contribution to the CE ν NS cross section.
We presented microscopic nuclear physics calculations of charge and weak nuclear form factors and the CE ν NS cross section on 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb nuclei. We obtain neutron (proton) densities and weak (charge) form factors by solving the Hartree–Fock equations with a Skyrme (SkE2) nuclear potential. Our predictions for 208 Pb and 40 Ar charge form factors describe elastic electron scattering data remarkably well.
After validating 40 Ar charge form factor calculations, we make predictions for the 40 Ar weak form factor. Thereby, we calculate differential cross section as a function of recoil energy and neutrino scattering angle. We attempt to gauge the level of theoretical uncertainty pertaining to the description of 40 Ar form factor and CE ν NS cross section by comparing relative differences between recent nuclear theory and widely-used phenomenological form factor predictions. We compare our 40 Ar prediction with recent measurements of the COHERENT collaboration. Future precise measurements of CE ν NS with ton and multi–ton detectors will aid in constraining neutron densities and weak nuclear form factor that will in turn improve prospects of extracting new physics through CE ν NS.
Furthermore, we calculate inelastic charged–current and neutral–current cross section on 40 Ar within the same formalism, and comparing the strength of coherent and inelastic processes. We present total and differential cross sections as a function of excitation energy and lepton scattering angle for neutrino energy relevant for pion decay–at–rest neutrinos. CE ν NS experiments at stopped-pion sources are well-suited to measure these inelastic cross sections and can provide powerful constraints on supernova detection capabilities of future kiloton liquid argon experiments.

Author Contributions

Conceptualization, N.V.D., V.P., H.R. and N.J.; Methodology, N.V.D., V.P., H.R. and N.J.; Formal analysis, N.V.D., V.P., H.R. and N.J.; Writing—original draft, N.V.D., V.P., H.R. and N.J.; Writing—review & editing, N.V.D., V.P., H.R. and N.J. All authors have read and agreed to the published version of the manuscript.

Funding

N.V.D. and N.J. are supported by the Research Foundation Flanders (FWO–Flanders). V.P. and H.R. acknowledge the support from US DOE under grant DE-SC0009824. This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics.

Acknowledgments

We thank S. Bacca, J. Yang and M. Hoferichter for providing their calculations for comparison. We thank Richard Van de Water for fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Freedman, D.Z. Coherent Neutrino Nucleus Scattering as a Probe of the Weak Neutral Current. Phys. Rev. D 1974, 9, 1389. [Google Scholar] [CrossRef]
  2. Akimov, D. et al. [COHERENT Collaboration] Observation of Coherent Elastic Neutrino-Nucleus Scattering. Science 2017, 357, 1123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Akimov, D. et al. [COHERENT Collaboration] COHERENT Collaboration data release from the first observation of coherent elastic neutrino-nucleus scattering. arXiv 2018, arXiv:1804.09459. [Google Scholar]
  4. Akimov, D. et al. [COHERENT Collaboration] First Constraint on Coherent Elastic Neutrino-Nucleus Scattering in Argon. Phys. Rev. D 2019, 100, 115020. [Google Scholar] [CrossRef] [Green Version]
  5. Akimov, D. et al. [COHERENT Collaboration] First Measurement of Coherent Elastic Neutrino-Nucleus Scattering on Argon. Phys. Rev. Lett. 2021, 126, 012002. [Google Scholar] [CrossRef] [PubMed]
  6. Akimov, D. et al. [COHERENT Collaboration] Measurement of the Coherent Elastic Neutrino-Nucleus Scattering Cross Section on CsI by COHERENT. Phys. Rev. Lett. 2022, 129, 081801. [Google Scholar] [CrossRef]
  7. Liao, J.; Marfatia, D. Repulsive baryonic interactions and lattice QCD observables at imaginary chemical potential. Phys. Lett. B 2017, 775, 47–54. [Google Scholar] [CrossRef]
  8. Dent, J.B.; Dutta, B.; Liao, S.; Newstead, J.L.; Strigari, L.E.; Walker, J.W. Accelerator and reactor complementarity in coherent neutrino-nucleus scattering. Phys. Rev. D 2018, 97, 035009. [Google Scholar] [CrossRef] [Green Version]
  9. Aristizabal Sierra, D.; Rojas, N.; Tytgat, M. Neutrino non-standard interactions and dark matter searches with multi-ton scale detectors. J. High Energy Phys. 2018, 2018, 197. [Google Scholar] [CrossRef] [Green Version]
  10. Denton, P.B.; Farzan, Y.; Shoemaker, I.M. Testing large non-standard neutrino interactions with arbitrary mediator mass after COHERENT data. J. High Energy Phys. 2018, 2018, 37. [Google Scholar] [CrossRef] [Green Version]
  11. Kosmas, T.; Papoulias, D.; Tortola, M.; Valle, J. Probing light sterile neutrino signatures at reactor and Spallation Neutron Source neutrino experiments. Phys. Rev. D 2017, 96, 063013. [Google Scholar] [CrossRef] [Green Version]
  12. Blanco, C.; Hooper, D.; Machado, P. Constraining Sterile Neutrino Interpretations of the LSND and MiniBooNE Anomalies with Coherent Neutrino Scattering Experiments. arXiv 2019, arXiv:1901.08094. [Google Scholar] [CrossRef]
  13. Aristizabal Sierra, D.; De Romeri, V.; Rojas, N. CP violating effects in coherent elastic neutrino-nucleus scattering processes. J. High Energy Phys. 2019, 9, 69. [Google Scholar] [CrossRef] [Green Version]
  14. Cadeddu, M.; Giunti, C.; Li, Y.; Zhang, Y. Average CsI Neutron Density Distribution from COHERENT Data. Phys. Rev. Lett. 2018, 120, 072501. [Google Scholar] [CrossRef] [Green Version]
  15. Ciuffoli, E.; Evslin, J.; Fu, Q.; Tang, J. P-wave contributions to Bψππ decays in the perturbative QCD approach. Phys. Rev. D 2018, 97, 113003. [Google Scholar] [CrossRef] [Green Version]
  16. Aristizabal Sierra, D.; Liao, J.; Marfatia, D. Impact of form factor uncertainties on interpretations of coherent elastic neutrino-nucleus scattering data. J. High Energy Phys. 2019, 2019, 141. [Google Scholar] [CrossRef] [Green Version]
  17. Papoulias, D.; Kosmas, T.; Sahu, R.; Kota, V.; Hota, M. Constraining nuclear physics parameters with current and future COHERENT data. Phys. Lett. B 2020, 800, 135133. [Google Scholar] [CrossRef]
  18. Aguilar-Arevalo, A.A. et al. [CCM Collaboration] First dark matter search results from Coherent CAPTAIN-Mills. Phys. Rev. D 2022, 106, 012001. [Google Scholar] [CrossRef]
  19. Coloma, P.; Esteban, I.; Gonzalez-Garcia, M.C.; Menendez, J. Determining the nuclear neutron distribution from Coherent Elastic neutrino-Nucleus Scattering: Current results and future prospects. J. High Energy Phys. 2020, 2020, 30. [Google Scholar] [CrossRef]
  20. Aguilar-Arevalo, A.; Bertou, X.; Bonifazi, C.; Butner, M.; Cancelo, G.; Vázquez, A.C.; Vergara, B.C.; Chavez, C.R.; Da Motta, H.; D’Olivo, J.C.; et al. Results of the engineering run of the Coherent Neutrino Nucleus Interaction Experiment (CONNIE). JINST 2016, 11, P07024. [Google Scholar] [CrossRef] [Green Version]
  21. Agnolet, G. et al. [MINER Collaboration] Background studies for the MINER Coherent Neutrino Scattering reactor experiment. Nucl. Instrum. Meth. A 2017, 853, 53. [Google Scholar] [CrossRef] [Green Version]
  22. Belov, V.; Brudanin, V.; Egorov, V.; Filosofov, D.; Fomina, M.; Gurov, Y.; Korotkova, L.; Lubashevskiy, A.; Medvedev, D.; Pritula, R.; et al. The νGeN experiment at the Kalinin Nuclear Power Plant. J. Instrum. 2015, 10, P12011. [Google Scholar] [CrossRef]
  23. Strauss, R.; Rothe, J.; Angloher, G.; Bento, A.; Gütlein, A.; Hauff, D.; Kluck, H.; Mancuso, M.; Oberauer, L.; Petricca, F.; et al. The ν-cleus experiment: A gram-scale fiducial-volume cryogenic detector for the first detection of coherent neutrino-nucleus scattering. Eur. Phys. J. C 2017, 77, 506. [Google Scholar] [CrossRef] [Green Version]
  24. Billard, J.; Carr, R.; Dawson, J.; Figueroa-Feliciano, E.; Formaggio, J.A.; Gascon, J.; Heine, S.T.; De Jesus, M.; Johnston, J.; Lasserre, T.; et al. Coherent neutrino scattering with low temperature bolometers at Chooz reactor complex. J. Phys. G 2017, 44, 105101. [Google Scholar] [CrossRef] [Green Version]
  25. Wong, H.T. Neutrino-nucleus coherent scattering and dark matter searches with sub-keV germanium detector. Nucl. Phys. A 2010, 844, 229c–233c. [Google Scholar] [CrossRef]
  26. Choi, J.J.; Jeon, E.J.; Kim, J.Y.; Kim, K.W.; Kim, S.H.; Kim, S.K.; Kim, Y.D.; Ko, Y.J.; Koh, B.C.; Ha, C.; et al. Exploring coherent elastic neutrino-nucleus scattering using reactor electron antineutrinos in the NEON experiment. Eur. Phys. J. C 2023, 83, 226. [Google Scholar] [CrossRef]
  27. Akindele, O.A.; Berryman, J.M.; Bowden, N.S.; Carr, R.; Conant, A.J.; Huber, P.; Langford, T.J.; Link, J.M.; Littlejohn, B.R.; Fernandez-Moroni, G.; et al. High Energy Physics Opportunities Using Reactor Antineutrinos. arXiv 2022, arXiv:2203.07214. https://arxiv.org/abs/2203.07214. [Google Scholar]
  28. Hofstadter, R. Electron Scattering and Nuclear Structure. Rev. Mod. Phys. 1956, 28, 214. [Google Scholar] [CrossRef]
  29. De Vries, H.; De Jager, C.W.; De Vries, C. Nuclear charge-density-distribution parameters from elastic electron scattering. Atom. Data Nucl. Data Tabl. 1987, 36, 495. [Google Scholar] [CrossRef]
  30. Fricke, G.; Bernhardt, C.; Heilig, K.; Schaller, L.A.; Schellenberg, L.; Shera, E.B.; de Jager, C.W. Nuclear Ground State Charge Radii from Electromagnetic Interactions. Atom. Data Nucl. Data Tabl. 1995, 60, 177. [Google Scholar] [CrossRef] [Green Version]
  31. Angeli, I.; Marinova, K. Table of experimental nuclear ground state charge radii: An update. Atom. Data Nucl. Data Tabl. 2013, 99, 69. [Google Scholar] [CrossRef]
  32. Thiel, M.; Sfienti, C.; Piekarewicz, J.; Horowitz, C.; Vanderhaeghen, M. Neutron skins of atomic nuclei: Per aspera ad astra. J. Phys. G 2019, 46, 093003. [Google Scholar] [CrossRef] [Green Version]
  33. Donnelly, T.; Dubach, J.; Sick, I. Isospin dependences in parity-violating electron scattering. Nucl. Phys. A 1989, 503, 589. [Google Scholar] [CrossRef]
  34. Abrahamyan, S. et al. [PREX Collaboration] New Measurements of the Transverse Beam Asymmetry for Elastic Electron Scattering from Selected Nuclei. Phys. Rev. Lett. 2012, 108, 112502. [Google Scholar]
  35. Horowitz, C.J.; Ahmed, Z.; Jen, C.M.; Rakhman, A.; Souder, P.A.; Dalton, M.M.; Liyanage, N.; Paschke, K.D.; Saenboonruang, K.; Silwal, R.; et al. Weak charge form factor and radius of 208Pb through parity violation in electron scattering. Phys. Rev. C 2012, 85, 032501. [Google Scholar] [CrossRef] [Green Version]
  36. Kumar, K.S. Electroweak probe of neutron skins of nuclei. Ann. Phys. 2020, 412, 168012. [Google Scholar] [CrossRef]
  37. Patton, K.; Engel, J.; McLaughlin, G.C.; Schunck, N. Neutrino-nucleus coherent scattering as a probe of neutron density distributions. Phys. Rev. C 2012, 86, 024612. [Google Scholar] [CrossRef] [Green Version]
  38. Klein, S.; Nystrand, J. Exclusive vector meson production in relativistic heavy ion collisions. Phys. Rev. C 1999, 60, 014903. [Google Scholar] [CrossRef] [Green Version]
  39. Helm, R.H. Inelastic and Elastic Scattering of 187-Mev Electrons from Selected Even-Even Nuclei. Phys. Rev. 1956, 104, 1466. [Google Scholar] [CrossRef]
  40. Payne, C.G.; Bacca, S.; Hagen, G.; Jiang, W.; Papenbrock, T. Coherent elastic neutrino-nucleus scattering on from first principles. Phys. Rev. C 2019, 100, 061304. [Google Scholar] [CrossRef] [Green Version]
  41. Yang, J.; Hernandez, J.A.; Piekarewicz, J. Electroweak probes of ground state densities. Phys. Rev. C 2019, 100, 054301. [Google Scholar] [CrossRef] [Green Version]
  42. Co’, G.; Anguiano, M.; Lallena, A. Nuclear structure uncertainties in coherent elastic neutrino-nucleus scattering. J. Cosmol. Astropart. Phys. 2020, 4, 44. [Google Scholar]
  43. Hoferichter, M.; Menéndez, J.; Schwenk, A. Coherent elastic neutrino-nucleus scattering: EFT analysis and nuclear responses. Phys. Rev. D 2020, 102, 074018. [Google Scholar] [CrossRef]
  44. Tomalak, O.; Machado, P.; Pandey, V.; Plestid, R. Flavor-dependent radiative corrections in coherent elastic neutrino-nucleus scattering. J. High Energy Phys. 2021, 2021, 097. [Google Scholar] [CrossRef]
  45. Tanabashi, M. et al. [Particle Data Group] Review of Particle Physics. Phys. Rev. D 2018, 98, 030001. [Google Scholar] [CrossRef] [Green Version]
  46. Abdullah, M.; Aristizabal Sierra, D.; Dutta, B.; Strigari, L.E. Coherent Elastic Neutrino-Nucleus Scattering with directional detectors. arXiv 2020, arXiv:2003.11510. [Google Scholar] [CrossRef]
  47. Van Dessel, N.; Jachowicz, N.; Nikolakopoulos, A. Forbidden transitions in neutral- and charged-current interactions between low-energy neutrinos and argon. Phys. Rev. C 2019, 100, 055503. [Google Scholar] [CrossRef] [Green Version]
  48. Ryckebusch, J.; Waroquier, M.; Heyde, K.; Moreau, J.; Ryckbosch, D. An RPA model for the description of one-nucleon emission processes and application to 16O(γ, N) reactions. Nucl. Phys. A 1988, 476, 237. [Google Scholar] [CrossRef]
  49. Ryckebusch, J.; Heyde, K.; Van Neck, D.; Waroquier, M. Aspects of the final-state interaction and long-range correlations in quasi-elastic (e, e’p) and (e, e’n) reactions. Nucl. Phys. A 1989, 503, 694. [Google Scholar] [CrossRef]
  50. Jachowicz, N.; Rombouts, S.; Heyde, K.; Ryckebusch, J. Cross sections for neutral-current neutrino-nucleus interactions: Applications for 12C and 16O. Phys. Rev. C 1999, 59, 3246. [Google Scholar] [CrossRef] [Green Version]
  51. Jachowicz, N.; Heyde, K.; Ryckebusch, J.; Rombouts, S. Continuum random phase approximation approach to charged-current neutrino-nucleus scattering. Phys. Rev. C 2002, 65, 025501. [Google Scholar] [CrossRef]
  52. Jachowicz, N.; Heyde, K.; Ryckebusch, J. Cross sections for neutral-current neutrino scattering on 208 Pb. Phys. Rev. C 2002, 66, 055501. [Google Scholar] [CrossRef]
  53. Jachowicz, N.; Vantournhout, K.; Ryckebusch, J.; Heyde, K. Identifying Neutrinos and Antineutrinos in Neutral-Current ScatteringReactions. Phys. Rev. Lett. 2004, 93, 082501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Jachowicz, N.; McLaughlin, G. Reconstructing supernova-neutrino spectra using low-energy beta beams. Phys. Rev. Lett. 2006, 96, 172301. [Google Scholar] [CrossRef] [Green Version]
  55. Pandey, V.; Jachowicz, N.; Ryckebusch, J.; Van Cuyck, T.; Cosyn, W. Quasielastic contribution to antineutrino-nucleus scattering. Phys. Rev. C 2014, 89, 024601. [Google Scholar] [CrossRef] [Green Version]
  56. Pandey, V.; Jachowicz, N.; Van Cuyck, T.; Ryckebusch, J.; Martini, M. Low-energy excitations and quasielastic contribution to electron-nucleus and neutrino-nucleus scattering in the continuum random-phase approximation. Phys. Rev. C 2015, 92, 024606. [Google Scholar] [CrossRef] [Green Version]
  57. Pandey, V.; Jachowicz, N.; Martini, M.; González-Jiménez, R.; Ryckebusch, J.; Van Cuyck, T.; Van Dessel, N. Impact of low-energy nuclear excitations on neutrino-nucleus scattering at MiniBooNE and T2K kinematics. Phys. Rev. C 2016, 94, 054609. [Google Scholar] [CrossRef] [Green Version]
  58. Van Dessel, N.; Jachowicz, N.; González-Jiménez, R.; Pandey, V.; Van Cuyck, T. A dependence of quasielastic charged-current neutrino-nucleus cross sections. Phys. Rev. C 2018, 97, 044616. [Google Scholar] [CrossRef] [Green Version]
  59. Nikolakopoulos, A.; Jachowicz, N.; Van Dessel, N.; Niewczas, K.; González-Jiménez, R.; Udías, J.M.; Pandey, V. Electron versus muon neutrino induced cross sections in charged current quasielastic processes. Phys. Rev. Lett. 2019, 123, 052501. [Google Scholar] [CrossRef]
  60. Van Dessel, N.; Nikolakopoulos, A.; Jachowicz, N. Lepton kinematics in low-energy neutrino-argon interactions. Phys. Rev. C 2020, 101, 045502. [Google Scholar] [CrossRef] [Green Version]
  61. Nikolakopoulos, A.; Pandey, V.; Spitz, J.; Jachowicz, N. Modeling quasielastic interactions of monoenergetic kaon decay-at-rest neutrinos. arXiv 2020, arXiv:2010.05794. [Google Scholar] [CrossRef]
  62. Tohyama, M. Application of extended random-phase approximation with ground-state correlations to collective excitations of 16O. J. Phys. Conf. Ser. 2014, 529, 012026. [Google Scholar] [CrossRef]
  63. Papoulias, D.; Kosmas, T.; Kuno, Y. Recent Probes of Standard and Non-standard Neutrino Physics With Nuclei. Front. Phys. 2019, 7, 191. [Google Scholar] [CrossRef]
  64. Antonello, M. et al. [MicroBooNE, LAr1-ND and ICARUS-WA104 Collaboration] A Proposal for a Three Detector Short-Baseline Neutrino Oscillation Program in the Fermilab Booster Neutrino Beam. arXiv 2015, arXiv:1503.01520. [Google Scholar]
  65. Abi, B. et al. [DUNE Collaboration] Deep Underground Neutrino Experiment (DUNE), Far Detector Technical Design Report, Volume II: DUNE Physics. arXiv 2020, arXiv:2002.03005. [Google Scholar]
  66. Amaudruz, P. et al. [DEAP-3600 Collaboration] First Results from the DEAP-3600 Dark Matter Search with Argon at SNOLAB. Phys. Rev. Lett. 2018, 121, 071801. [Google Scholar] [CrossRef] [Green Version]
  67. Agnes, P. et al. [DarkSide Collaboration] Low-Mass Dark Matter Search with the DarkSide-50 Experiment. Phys. Rev. Lett. 2018, 121, 081307. [Google Scholar] [CrossRef] [Green Version]
  68. Calvo, J. et al. [ArDM Collaboration] Backgrounds and pulse shape discrimination in the ArDM liquid argon TPC. J. Cosmol. Astropart. Phys. 2017, 3, 3. [Google Scholar]
  69. Hime, A. [MiniCLEAN Collaboration]. The MiniCLEAN Dark Matter Experiment. arXiv 2011, arXiv:1110.1005. [Google Scholar]
  70. Ottermann, C.R.; Schmitt, C.H.; Simon, G.G.; Borkowski, F.; Walther, V.H. Elastic electron scattering from 40Ar. Nucl. Phys. A 1982, 379, 396. [Google Scholar] [CrossRef]
  71. Duda, G.; Kemper, A.; Gondolo, P. Model Independent Form Factors for Spin Independent Neutralino-Nucleon Scattering from Elastic Electron Scattering Data. J. Cosmol. Astropart. Phys. 2007, 4, 12. [Google Scholar] [CrossRef]
  72. Lewin, J.; Smith, P. Review of mathematics, numerical factors, and corrections for dark matter experiments based on elastic nuclear recoil. Astropart. Phys. 1996, 6, 87–112. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (Left) Diagrammatic representation of the CE ν NS process where a single Z 0 boson is exchanged between neutrino and target nucleus. (Right) Diagrammatic representation of the inelastic neutrino-nucleus scattering where a single W + (CC) or Z 0 (NC) boson is exchanged between neutrino and target nucleus.
Figure 1. (Left) Diagrammatic representation of the CE ν NS process where a single Z 0 boson is exchanged between neutrino and target nucleus. (Right) Diagrammatic representation of the inelastic neutrino-nucleus scattering where a single W + (CC) or Z 0 (NC) boson is exchanged between neutrino and target nucleus.
Universe 09 00207 g001
Figure 2. Panels (a,b) represent proton and neutron densities of different nuclei obtained using the HF–SkE2 approach. Panels (c) through (d) represent charge and weak form factors for the different nuclei.
Figure 2. Panels (a,b) represent proton and neutron densities of different nuclei obtained using the HF–SkE2 approach. Panels (c) through (d) represent charge and weak form factors for the different nuclei.
Universe 09 00207 g002
Figure 3. The “weak-skin” form factor depicts the difference between the charge and weak form factors.
Figure 3. The “weak-skin” form factor depicts the difference between the charge and weak form factors.
Universe 09 00207 g003
Figure 4. (Left): the charge form factor of 208 Pb compared with elastic electron scattering data of Ref. [29]. (Right): the weak form factor of 208 Pb along with the single point measured by the PREX collaboration at the momentum transfer of q = 0.475 fm 1 [34,35]. Both form factors are compared with relativistic mean–field predictions of Yang et al. [41].
Figure 4. (Left): the charge form factor of 208 Pb compared with elastic electron scattering data of Ref. [29]. (Right): the weak form factor of 208 Pb along with the single point measured by the PREX collaboration at the momentum transfer of q = 0.475 fm 1 [34,35]. Both form factors are compared with relativistic mean–field predictions of Yang et al. [41].
Universe 09 00207 g004
Figure 5. Total CE ν NS cross sections for a set of nuclear targets obtained within the HF–SkE2 approach.
Figure 5. Total CE ν NS cross sections for a set of nuclear targets obtained within the HF–SkE2 approach.
Universe 09 00207 g005
Figure 6. CE ν NS cross section strength compared to CCQE and NCQE scattering cross sections for several nuclei, above particle emission threshold.
Figure 6. CE ν NS cross section strength compared to CCQE and NCQE scattering cross sections for several nuclei, above particle emission threshold.
Universe 09 00207 g006
Figure 7. (Left) The 40 Ar charge form factor predictions compared to elastic electron scattering data taken from Ref. [70], a comparison is also performed with the coupled–cluster theory predictions of Payne et al. [40] as well as with Klein–Nystrand [38] (standard and adapted) and Helm [39] form factors. (Right) The 40 Ar weak form factor predictions compared with calculations of Payne et al. [40], Yang et al. [41], Hoferichter et al. [43] and with the predictions of Klein–Nystrand [38] (standard and adapted) and Helm [39] form factors.
Figure 7. (Left) The 40 Ar charge form factor predictions compared to elastic electron scattering data taken from Ref. [70], a comparison is also performed with the coupled–cluster theory predictions of Payne et al. [40] as well as with Klein–Nystrand [38] (standard and adapted) and Helm [39] form factors. (Right) The 40 Ar weak form factor predictions compared with calculations of Payne et al. [40], Yang et al. [41], Hoferichter et al. [43] and with the predictions of Klein–Nystrand [38] (standard and adapted) and Helm [39] form factors.
Universe 09 00207 g007
Figure 8. The 40 Ar cumulative cross section as a function of q c u t o f f compared with calculations done using Payne et al. [40], Yang et al. [41], as well the Klein–Nystrand [38] (standard and adapted) and Helm [39] form factors.
Figure 8. The 40 Ar cumulative cross section as a function of q c u t o f f compared with calculations done using Payne et al. [40], Yang et al. [41], as well the Klein–Nystrand [38] (standard and adapted) and Helm [39] form factors.
Universe 09 00207 g008
Figure 9. Differential cross section on argon as a function of recoil energy and scattering angle.
Figure 9. Differential cross section on argon as a function of recoil energy and scattering angle.
Universe 09 00207 g009
Figure 10. Relative differences in the weak form factor and CE ν NS cross section predictions of Payne et al. [40], Yang et al. [41], Hoferichter et al. [43], Helm [39], Klein–Nystrand [38] and the adapted Klein–Nystrand [16,63], all with respect to HF–SkE2.
Figure 10. Relative differences in the weak form factor and CE ν NS cross section predictions of Payne et al. [40], Yang et al. [41], Hoferichter et al. [43], Helm [39], Klein–Nystrand [38] and the adapted Klein–Nystrand [16,63], all with respect to HF–SkE2.
Universe 09 00207 g010
Figure 11. (Left) The CE ν NS cross section on 40 Ar as a function of neutrino energy, recent flux–folded measurement by the COHERENT collaboration [5] is shown along with the flux-folded HF–SkE2 prediction. (Right) Flux–averaged CE ν NS cross sections as a function of neutron number for the 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb nuclei. We also show 40 Ar data measured by COHERENT [5].
Figure 11. (Left) The CE ν NS cross section on 40 Ar as a function of neutrino energy, recent flux–folded measurement by the COHERENT collaboration [5] is shown along with the flux-folded HF–SkE2 prediction. (Right) Flux–averaged CE ν NS cross sections as a function of neutron number for the 12 C, 16 O, 40 Ar, 56 Fe and 208 Pb nuclei. We also show 40 Ar data measured by COHERENT [5].
Universe 09 00207 g011
Figure 12. Charged-current (left) and neutral-current (right) inelastic cross section: total as a function of neutrino energy shown along with contributions from different multipoles (top panel), differential as a function of excitation energy (middle panel) and as a function of lepton scattering angle (bottom panel) for fixed neutrino energies, E ν = 30 and 50 MeV.
Figure 12. Charged-current (left) and neutral-current (right) inelastic cross section: total as a function of neutrino energy shown along with contributions from different multipoles (top panel), differential as a function of excitation energy (middle panel) and as a function of lepton scattering angle (bottom panel) for fixed neutrino energies, E ν = 30 and 50 MeV.
Universe 09 00207 g012
Table 1. Single –particle energies in 40 Ar, as provided by a HF calculation using the SkE2 interaction.
Table 1. Single –particle energies in 40 Ar, as provided by a HF calculation using the SkE2 interaction.
p / n i n i , l i , j i ε i (MeV) v i 2 # N
p1 1 s 1 / 2 −43.70291.002
p2 1 p 3 / 2 −31.44961.004
p3 1 p 1 / 2 −27.39211.002
p4 1 d 5 / 2 −17.70271.006
p5 2 s 1 / 2 −12.08221.002
p6 1 d 3 / 2 −10.92430.502
n1 1 s 1 / 2 −48.30471.002
n2 1 p 3 / 2 −35.20201.004
n3 1 p 1 / 2 −31.02471.002
n4 1 d 5 / 2 −21.10351.006
n5 2 s 1 / 2 −16.11161.002
n6 1 d 3 / 2 −14.02661.004
n7 1 f 7 / 2 −7.21080.252
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Van Dessel, N.; Pandey, V.; Ray, H.; Jachowicz, N. Cross Sections for Coherent Elastic and Inelastic Neutrino-Nucleus Scattering. Universe 2023, 9, 207. https://doi.org/10.3390/universe9050207

AMA Style

Van Dessel N, Pandey V, Ray H, Jachowicz N. Cross Sections for Coherent Elastic and Inelastic Neutrino-Nucleus Scattering. Universe. 2023; 9(5):207. https://doi.org/10.3390/universe9050207

Chicago/Turabian Style

Van Dessel, Nils, Vishvas Pandey, Heather Ray, and Natalie Jachowicz. 2023. "Cross Sections for Coherent Elastic and Inelastic Neutrino-Nucleus Scattering" Universe 9, no. 5: 207. https://doi.org/10.3390/universe9050207

APA Style

Van Dessel, N., Pandey, V., Ray, H., & Jachowicz, N. (2023). Cross Sections for Coherent Elastic and Inelastic Neutrino-Nucleus Scattering. Universe, 9(5), 207. https://doi.org/10.3390/universe9050207

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop