INT-PUB-20-026 Coherent elastic neutrino{nucleus scattering: EFT analysis and nuclear responses Martin Hoferichter,1, 2, ∗ Javier Men´endez,3, 4, y and Achim Schwenk5, 6, 7, z 1Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, 3012 Bern, Switzerland 2Institute for Nuclear Theory, University of Washington, Seattle, WA 98195-1550, USA 3Department of Quantum Physics and Astrophysics and Institute of Cosmos Sciences, University of Barcelona, Spain 4Center for Nuclear Study, The University of Tokyo, 113-0033 Tokyo, Japan 5Institut f¨urKernphysik, Technische Universit¨atDarmstadt, 64289 Darmstadt, Germany 6ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum f¨urSchwerionenforschung GmbH, 64291 Darmstadt, Germany 7Max-Planck-Institut f¨urKernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany The cross section for coherent elastic neutrino{nucleus scattering (CEνNS) depends on the re- sponse of the target nucleus to the external current, in the Standard Model (SM) mediated by the exchange of a Z boson. This is typically subsumed into an object called the weak form factor of the nucleus. Here, we provide results for this form factor calculated using the large-scale nuclear shell model for a wide range of nuclei of relevance for current CEνNS experiments, including ce- sium, iodine, argon, fluorine, sodium, germanium, and xenon. In addition, we provide the responses needed to capture the axial-vector part of the cross section, which does not scale coherently with the number of neutrons, but may become relevant for the SM prediction of CEνNS on target nuclei with nonzero spin. We then generalize the formalism allowing for contributions beyond the SM. In particular, we stress that in this case, even for vector and axial-vector operators, the standard weak form factor does not apply anymore, but needs to be replaced by the appropriate combination of the underlying nuclear structure factors. We provide the corresponding expressions for vector, axial-vector, but also (pseudo-)scalar, tensor, and dipole effective operators, including two-body- current effects as predicted from chiral effective field theory. Finally, we update the spin-dependent structure factors for dark matter scattering off nuclei according to our improved treatment of the axial-vector responses. I. INTRODUCTION from relativistic mean-field methods (RMF) [22, 23], even though results based on nonrelativistic energy-density CEνNS, suggested as a probe of the weak current as functionals are also available [24{26]. For argon, there early as 1974 [1], was finally observed by the COHER- is a recent first-principles calculation based on coupled- ENT collaboration in 2017 [2]. After the initial detec- cluster theory [27]. tion in CsI, also the scattering off argon has just been Here, we provide nuclear structure results for CEνNS, observed [3, 4]. With future advances in COHERENT extending large-scale nuclear shell model calculations de- and other experiments [5{13], the CEνNS process will veloped in the context of nuclear structure factors in soon develop into another sensitive low-energy probe of direct-detection searches for dark matter [28{36]. First, physics beyond the Standard Model (BSM) [14]. the level of agreement between the shell-model, the RMF, and coupled-cluster results suggests that the form factor A crucial input in interpreting the measured cross sec- uncertainties are not as severe as claimed in Ref. [37], tion is the response of the nucleus. If BSM constraints but in addition the shell-model approach also allows us are to be extracted, the nuclear structure has to be pro- to address the spin-dependent (SD) responses, which are vided from elsewhere. In fact, since the weak charge similar, but somewhat different to the ones in SD dark of the proton is small, the SM CEνNS process mainly matter searches. probes the nuclear neutron distribution, which is sig- To this end, we first derive the decomposition of the arXiv:2007.08529v2 [hep-ph] 23 Oct 2020 nificantly less constrained experimentally than the elec- tromagnetic charge distribution. Apart from CEνNS, cross section into Wilson coefficients of effective opera- the only direct information of the neutron distribution tors, hadronic matrix elements, and nuclear structure fac- comes from parity-violating electron scattering (PVES) tors. In the SM the effective operators just parameterize off lead [15, 16]. Accordingly, assuming the absence of a the Z-boson exchange, but this approach can be conve- significant BSM signal, the measured CEνNS cross sec- niently extended to include BSM effects. The hadronic tion could be used to constrain the neutron distribution matrix elements determine the hadronization of these op- instead [17{21]. Currently, the nuclear input used in the erators at the single-nucleon level, and finally the nuclear interpretation of CEνNS experiments is mainly derived structure factors take into account the many-body nu- clear matrix element of these single-nucleon currents. As a first step, we demonstrate how the standard weak form factor emerges when combining all these ingredients into a single object. However, this analysis shows that even ∗E-mail: [email protected] yE-mail: [email protected] for the coherent part of the nuclear response four different zE-mail: [email protected] underlying structure factors contribute to the cross sec- 2 tion. Therefore, in principle the weak form factor needs where we adopted the following conventions: neutrino to be modified as well when allowing for BSM effects in indices are suppressed throughout, indicating that oscil- the Wilson coefficients, since their contribution does not lation effects are usually negligible at the scale of CEνNS factorize. experiments, so that incoming and outgoing flavors are While for the dominant vector operators corrections understood to be identical. In comparison to the case beyond the single-nucleon currents are small [38, 39], of a dark-matter spin-1=2 particle [36] the number of since the magnetic-moment form factors happen to be operators is reduced by a factor of 2 when assuming kinematically suppressed, this is no longer true for the that the neutrino is left-handed. This is implemented axial-vector [29] or for scalar currents, see [32{34, 36, 40{ in Eq. (1) in terms of projectors PL = (1 γ5)=2, and 45] for the analogous effects in the case of dark matter given that observing chirality-violating effects− would re- scattering off nuclei. As long as the SM dominates, such quire right-handed neutrino beams (suppressed by tiny effects will only become relevant for CEνNS once ex- neutrino masses), we will not consider the opposite chi- periments become sensitive to SD responses. Otherwise, rality in the following. mainly the limits on scalar operators would be affected, With these conventions the set of operators is then sim- but in CEνNS such contributions are suppressed due to ilar to the dark-matter case: at dimension-5 level there the need for right-handed neutrinos and lack of interfer- is a single (dipole) operator involving the photon field ence with the SM. strength tensor Fµν . At dimension 6 we have the vector The paper is organized as follows: in Sec. II we first and axial-vector operators already present in the SM, as review the necessary formalism, regarding effective oper- well as the tensor operator. Introducing quark masses for ators, hadronic matrix elements, and nuclear responses, renormalization-group invariance, the scalar and pseu- with some details of the multipole expansion and nuclear- doscalar operators would be counted as dimension-7. The µ structure calculations summarized in the appendices. In operator involving the QCD trace anomaly θµ would also Sec. III we first introduce the charge and weak form be of dimension 7. We have already rewritten the gluon a µν factors in the context of electron scattering, before dis- term Gµν Ga in terms of this operator (we will not con- cussing the application to CEνNS. In particular, we sider the operators involving the dual field strength ten- ~a present an improved treatment of the axial-vector re- sor Gµν or the photon field strength). In particular, sponses both for CEνNS and dark matter. In Sec. IV we already integrated out the heavy quarks [47] and ab- we discuss how the nuclear responses need to be adapted sorbed their effect into when considering SM extensions, before concluding with 1 X C0S = CS CS ; (2) a summary in Sec. V. g g − 12π Q Q=c;b;t S where Cg is the original coefficient of the II. FORMALISM a µν νP¯ Lν αsGµν Ga gluon operator and we used the relation A. Effective Lagrangians µ X 9 a µν 2 θ = mqqq¯ αsG G + α (3) µ − 8π µν a O s As a first step, we review the operator basis for q 1 CEνNS [36, 46] in the transition. Elsewhere, the sum over q in princi- ple refers to all quark species, but in practice we will (5) µν = CF νσ¯ PLνFµν ; L restrict the analysis to the light quarks q = u; d; s. Fi- X (6) = CV νγ¯ µP ν qγ¯ q + CAνγ¯ µP ν qγ¯ γ q nally, there are operators with derivatives acting on the q L µ q L µ 5 neutrino fields (in analogy to the spin-2 operator for L q dark matter), but we will concentrate on the more fre- T µν + Cq νσ¯ PLν qσ¯ µν q ; quently considered operators in Eq. (1). We note that the dimensional counting is not unambiguous regarding (7) X S 8π 0S = C + C νP¯ Lν mqqq¯ the quark mass mq, e.g., sometimes the tensor operator L q 9 g q is introduced at dimension-7 by adding a factor mq in X P 8π 0S µ this operator as well [46] (the notation in Eq. (1) follows + C νP¯ Lν mqqiγ¯ 5q C νP¯ Lν θ ; (1) q − 9 g µ Ref.