Implications of the first detection of coherent elastic neutrino-nucleus scattering (CEvNS) with Liquid Argon

The CENNS-10 experiment of the COHERENT collaboration has recently reported the first detection of coherent-elastic neutrino-nucleus scattering (CEvNS) in liquid Argon with more than $3 \sigma$ significance. In this work, we exploit the new data in order to probe various interesting parameters which are of key importance to CEvNS within and beyond the Standard Model. A dedicated statistical analysis of these data shows that the current constraints are significantly improved in most cases. We derive a first measurement of the neutron rms charge radius of Argon, and also an improved determination of the weak mixing angle in the low energy regime. We also update the constraints on neutrino non-standard interactions and electromagnetic properties with respect to those derived from the first COHERENT-CsI data.


INTRODUCTION
Over the past few decades, neutrino physics has entered the precision era and is steadily approaching a full description of the three neutrino oscillation picture [1]. Although the next generation of long-baseline experiments such as DUNE and Hyper-Kamiokande will play a key role in settling the open issues, some degeneracies are likely to remain. Moreover, beyond the challenges associated to the three-neutrino parameters, such as the atmospheric octant, the presence of leptonic CP violation, the neutrino mass ordering, and the absolute scale of the neutrino mass, there are a number of less standard, but important, neutrino properties to pin down.
Here we show how the recent confirmation of the CEvNS process by the COHERENT collaboration using a 24 kg liquid Argon (LAr) detector at the SNS [38], after collecting 6.12 GWh data substantially improves the sensitivity of a number of key weak interaction measurements. The reported results have shown a greater than 3σ preference in favor of CEvNS.
This paper is organized as follows. In Sec. 2 we set up the notation and summarize briefly the CEvNS formalism. Then, in Sec. 3 we present the sensitivities on various parameters of the electroweak interaction in the SM and beyond, obtained using the recent CENNS-10 results. Finally, our main conclusions are summarized in Sec. 4.

BASICS
Proposed more than forty years ago by Freedman [39], this neutral current process is characterized by a cross section that increases as N 2 , where N being the number of neutrons in the nucleus: where G F denotes the Fermi constant, T A is the nucleus kinetic energy, E ν the neutrino energy and Q V W the vector weak charge written in the form Here, Z and N are the number of protons and neutrons in the nucleus while the neutral current vector couplings, g p,n V , are given by with the weak mixing angle taken in the MS scheme, i.e. sin 2 θ W ≡ŝ 2 Z = 0.2312. Finally, F p,n (Q) 2 stands for the nuclear form factors for protons and neutrons respectively, for which we employ the well-known Helm parametrization where the magnitude of the three-momentum transfer is Q = √ 2m A T A , j 1 denotes the spherical Bessel function of order one and R 2 0 = 5 3 (R 2 p,n − 3s 2 ) with R n = 3.36 fm (R p = 3.14 fm) denoting the neutron (proton) rms radius and s = 0.9 fm.
To analyze the recent results of the liquid Argon detector reported by the COHERENT collaboration, we will consider in this work a neutrino flux arriving to the detector from the SNS at ORNL, given by three different components from the π + decay at rest. These are a "prompt" monoenergetic muon neutrino signal given by and the "delayed" neutrino flux composed of muon antineutrinos and electron neutrinos with the normalization factor η given by η = rN POT /4πL 2 . The number of protons on target (POT) corresponding to the 6.12 GWh exposure is N POT = 1.37 × 10 23 , r = 0.08 denotes the produced neutrinos per POT and L = 27.5 m is the CENNS-10 baseline.
In the present work, prompted by the importance of the liquid Argon detector of the COHERENT collaboration, we study relevant implications of the measurement for different standard and new physics scenarios.

THE ANALYSIS
As a first step, we simulate the CEvNS signal at the LAr detector using the efficiency function 1 A(T A ), corresponding to the "Analysis A" of Ref. [38]. We then convert the nuclear recoil spectrum into electron recoil space through the reported quenching factor QF(T A ). Following Ref. [38], we evaluate the number of events in the region of interest, i.e. for recoil energies below 120 keV nr ≈ 30 keV ee . Our simulated energy spectrum is shown in Fig. 1 as a function of the electron recoil energy, T er . Given the reliability of our simulated CENNS-10 signal, we perform a statistical analysis of the recent liquid Argon CEvNS measurement with the ultimate goal of probing important observables such as the neutron mean radius, the electroweak mixing angle sin 2 θ W , as well as to constrain new physics parameters. Before displaying our results in the next sections, here we discuss the general procedure we followed to investigate the different scenarios.
To test the sensitivity of the experiment to observables under study, we have performed a χ 2 analysis, minimizing the function: with N meas = 159 denoting the number of the measured events from the fit of "Analysis A" in Ref. [38] and N theor (X) being the theoretical prediction. Here, the argument X represents the set of parameters to be tested, such as the weak mixing angle, the neutron rms radius or NSI parameters. Here we note that we have successfuly calibrated our procedure with the one discussed in the recent CENNS-10 result. In Eq.(8), the statistical uncertainty is given by σ stat = √ N meas + N BRN , where N BRN = 563 represents the number of background events due to beam related neutrons (BRN). The parameter α quantifies the normalization, that has a systematic error σ α = 8.5% (for details, see Ref. [38]).

Standard electroweak and nuclear physics
One of the most important parameters of the SM is the weak mixing angle, that is measured with great accuracy at the Z peak. At low energies, however, the existing measurements are less precise but still very relevant given the prediction of an increase of about 3% in its value due to radiative corrections. Moreover, any deviation from the SM prediction for this value would be a signature of new physics. With this motivation in mind, we have performed a χ 2 analysis for this SM parameter using the liquid Argon data. Our corresponding results are shown in the left panel of Fig. 2, where one can see a notable improvement with respect to the previous determination of sin 2 θ W from the COHERENT-CsI data (see Ref. [21]). The new measurement of the weak mixing angle, derived from the CENNS-10 data at 90% C.L. read Another very useful standard information that can be obtained from the CEvNS interaction is the neutron mean radius R n for the Argon isotope. Although there are theoretical predictions for this value, its direct determination can facilitate a better understanding of the CEvNS background at dark matter oriented experiments [40]. We have performed the corresponding analysis for this observable and the obtained result is illustrated in the right panel of Fig. 2. Although not directly comparable, the average neutron rms radius of CsI is also shown for the reader's convenience. On the other hand, it is worthwhile to notice that R n is more severely constrained for the case of Argon. The obtained 90% C.L. limit for the neutron rms radius for Argon reads R n < 4.33 fm .
We should note that this provides the first experimental determination of the neutron radius in Argon. In addition, comparing with the previous result on CsI [5] (see also Refs. [6,41]), one sees that the level of precision seems somewhat improved.

Non-standard Interactions
Most models that try to explain the neutrino mass pattern predict a modification of the V − A couplings predicted by the SM. In many cases, the corrections are expected to be negligible, like in the most simple type I seesaw model while, in other cases, there could be a relatively larger signal, as for the linear and inverse seesaw cases. A large family of new physics models can be phenomenologically described using the formalism of NSI, that modify the neutral current SM Lagrangian through the contribution [42][43][44] where f corresponds to an elementary fermion that, in the case of CEvNS, reduces to the quarks of the first family, f = {u, d}. α and β denote the neutrino flavors {e, µ, τ }, P X the left and right chirality projectors P L,R , and ε f P αβ are the couplings that quantify the relative strength of the NSI. Due to the presence of these new interactions, the weak charge of the CEvNS reaction is modified according to the substitution Q V W → Q V NSI in Eq. (1), with the NSI charge given by As it has already been noticed [45], CEvNS reaction is sensitive to the NSI parameters and, therefore, it can provide important information to probe the so-called LMA-Dark so-  lution [46]. The results of our χ 2 analysis, for one NSI parameter at a time, are shown in Fig. 3, both for the flavor changing (left panel) and the non-universal case (right panel). We can see from this figure that the sensitivities on the flavor changing NSI parameters are only marginally improved with respect to the previous CsI case. This is due to the detection of a larger number of events with respect to the SM prediction. Although the excess is below one standard deviation, still the preferred value for the flavor changing parameters are non-zero, as can be seen in Fig. 3. In the case of non-universal NSI, one can see the improvement in their restriction in comparison with the first CEvNS detection.
We can go one step further in the analysis and study more general restrictions on NSI parameters [45,47]. For example, the new CENNS-10 measurement can be used to constrain pairs of NSI parameters, allowing us to seek for possible correlations between them, as shown in Fig. 4, where the constraints on the two-dimensional parameter space of non-universal and flavor-changing NSI couplings with d quark (ε dV ee , ε dV τ e ) are given (left panel), as well as those for the case of the non-universal NSI couplings with d and u quarks (ε dV ee , ε uV ee ). In both cases one can appreciate the improvement in the determination of the parameters. The good consistency of our analysis with the available results of Ref. [38] is evident from the right panel of this figure.

Electromagnetic properties
The discovery of neutrino oscillations with solar, atmospheric, reactor and accelerator neutrinos [1] has set a new milestone in particle physics pointing to the existence of massive neutrinos, hence constituting the clearest signature of new physics beyond the SM. The nonvanishing neutrino mass stands out as the best motivation for non-trivial electromagnetic (EM) neutrino properties. Indeed, the expansion of the EM neutrino vertex (for details see Ref. [48]) yields two main phenomenological parameters, namely the neutrino magnetic moment and the neutrino charge radius [49].
The neutrino magnetic moment µ να with α = e, µ, τ is a flavor-dependent quantity which, for the case of scattering experiments, is usually expressed in the mass basis [50,51]. Due to the helicity-violating nature of the EM cross section, there is no interference with the SM one given in Eq.(1) and yields an additive contribution of the form where the flavor index has been dropped. For sufficient low detection threshold, the latter cross section leads to an enhancement of the recoil spectrum, i.e. a feature that is not expected for the other types of new physics considered in this work. Note that, for the case of Majorana neutrinos, only transition magnetic moments are expected, and the corresponding sensitivities from the analysis of neutrino-electron scattering [52] and CEvNS [23] experiments have already been given in the literature. Here, for simplicity, we only consider flavor-dependent effective neutrino magnetic moments. The resulted sensitivity profiles relevant to µ νe , µ νµ and µν µ from the analysis of the recent CENNS-10 data are shown in the We notice that the LAr sensitivity on the effective neutrino magnetic moment is not drastically improved in comparison to the one reached with the first COHERENT-CsI data [21]. This is due to the fact that the CENNS-10 result of N meas = 159 events is by about 30 events larger than the 130 events expected in the Standard Model. Therefore, even though the systematic uncertainties of CENNS-10 experiment are better compared to the first CsI measurement, these extra 30 events in the χ 2 translate into a finite neutrino magnetic moment contribution, thus weakening the limits.
We now turn our attention to the neutrino charge radius r 2 να with α = e, µ, τ being the flavor index. Likewise the neutrino magnetic moment, r 2 να is also required to be expressed in the mass basis through a rotation with the lepton mixing matrix [22]. Being a helicitypreserving quantity, its impact to the SM cross section is simply taken as a shift on the weak mixing angle according to Note that only the proton coupling g p V , proportional to the number of protons, interacts with the charge radius 2 . On the other hand, when antineutrinos are involved, both the g p V and r 2 να change sign and, therefore, Eq.(15) holds for both neutrinos and antineutrinos, a result that is consistent with Ref. [21]. From our analysis of the recent CENNS-10 data, we obtain the sensitivity to the neutrino charge radii corresponding to r 2 νe , r 2 νµ and r 2 νµ , as shown in the right panel of Fig. 5. The extracted constraints at 90% C.L. read in units of 10 −32 cm 2 . Note that, comparing with the corresponding results extracted from the CsI case, the precision of the new determinations derived from the new LAr measurement are significantly improved since now the allowed regions are narrower, and appear as two separate intervals. At this point, and as we did for the NSI case, we have performed a combined analysis allowing for several parameters being non-zero at a time. In particular, we have chosen to probe the parameter space of neutrino magnetic moments (µ να , µ ν β ) and neutrino charge radii ( r 2 να , r 2 ν β ), by considering the simultaneous presence of two of them, i.e. α = β. The corresponding results are presented in Fig. 6, where we also compare with those derived from the 2017 data release of the COHERENT-CsI measurement. The improvement obtained with the most recent data is more than evident.

CONCLUSIONS
We have analyzed the recent results of the CENNS-10 detector subsystem of the COHER-ENT collaboration that led to the first detection of CEvNS in LAr. Through a dedicated statistical analysis, taking into account the available information from Ref. [38], we have shown that this new measurement typically leads to improved sensitivities with respect to the first COHERENT-CsI measurement in 2017. Specifically, we have presented an improved determination of the weak mixing angle, as well as the first ever determination of the 40 Ar neutron rms radius. Turning to new physics, we have derived the constraints on non-universal as well as flavor-changing NSI imposed by this new data release. Moreover, concerning neutrino electromagnetic properties, we have found only minor improvement of the sensitivity to neutrino magnetic moments. In contrast, we have now an indication for finite neutrino-charge radii.