coherent scattering of neutrinos on complex nuclei was proposed long ago @xcite as a prominent probe to study neutral - current ( nc ) @xmath0-nucleus processes , but up to now no events have been experimentally measured .
neutrino detection , constitutes an excellent probe to search for a plethora of conventional neutrino physics applications and new - physics open issues @xcite . in principle
, low - energy astrophysical and laboratory neutrino searches provide crucial information towards understanding the underling physics of the fundamental electroweak interactions within and beyond the sm @xcite .
well - known neutrino sources include ( i ) supernova neutrinos ( with energies up to 60 - 100 mev ) and ( ii ) laboratory neutrinos ( with energies up to 52.8 mev ) emerging from stopped - pion and muon decays at muon factories ( fermilab , psi , jparc , etc . ) and at the spallation neutron source ( sns ) at oak ridge national lab @xcite .
recently , it became feasible @xcite to detect neutrinos by exploiting the neutral current interactions and measuring the nuclear recoil signals through the use of very low threshold - energy detectors @xcite .
to this purpose , great experimental effort has been put and new experiments have been proposed to be performed at facilities with stopped - pion neutrino beams , based on promising nuclear detectors like those of the coherent experiment @xcite and others @xcite at the sns , or alternative setups at the booster neutrino beam ( bnb ) at fermilab @xcite .
the nuclear @xmath0-detectors adopted by the relevant experiments include liquid noble gases , such as @xmath2ne , @xmath3ar , @xmath4xe as well as , @xmath5ge and csi[na ] detection materials @xcite .    on the theoretical side ,
the @xmath0-signals of low - energy neutrinos , expected to be recorded in sensitive nuclear detectors @xcite , could be simulated through nuclear calculations of @xmath0-nucleus scattering cross sections .
such results may provide useful information relevant for the evolution of distant stars , the core collapse supernovae , explosive nucleosynthesis and other phenomena @xcite .
in fact , coherent neutral current @xmath0-nucleus scattering events are expected to be observed by using the high intensity stopped - pion neutrino beams @xcite and nuclear targets for which recoil energies are of the order of a few to tens of kev , and therefore appropriate for detection of wimps @xcite , candidates of cold dark - matter @xcite .
such detectors are e.g. the supercdms @xcite , gerda @xcite and other multi - purpose detectors @xcite . for low - energies ,
the dominant vector components of nc interactions lead to a coherent contribution of all nucleons ( actually all neutrons ) in the target nucleus @xcite .
it is worth mentioning that , after the discovery @xcite of neutrino oscillations in propagation , the challenge of neutral and charged lepton flavour violation ( lfv ) is further investigated by extremely sensitive experiments @xcite searching for physics beyond the current standard model ( sm ) @xcite . to this end
, neutrino - nucleus coherent scattering experiments may probe new physics beyond the sm involved in exotic neutrino - nucleus interactions @xcite , an undoubtable signature of non - standard physics .
therefore , new data and insights will be provided to the physics of flavour changing neutral current ( fcnc ) processes , in the leptonic sector , in non - standard neutrino oscillation effects @xcite , in neutrino transition magnetic moments @xcite , in sterile neutrino search @xcite and others @xcite .
furthermore , such experimental sensitivity may also inspire advantageous probes to shed light on various open issues in nuclear astrophysics @xcite .    in recent works @xcite ,
neutral - current ( nc ) non - standard interactions ( nsi ) involving ( anti)neutrino scattering processes on leptons , nucleons and nuclei have been investigated .
the reactions of this type that take place in nuclei are represented by @xmath6 ( @xmath7 with @xmath8 ) .
it has been suggested @xcite that , theoretically the latter processes , can be studied with the same nuclear methods as the exotic clfv process of @xmath9 conversion in nuclei @xcite .
the corresponding lagrangians may be derived within the context of various extensions of the sm @xcite , like the four fermion contact interaction , seesaw model @xcite , left - right symmetric models @xcite , gluonic operator model @xcite , etc .
it is well - known that neutrino nsi may have rather significant impact in many areas of modern physics research and thus , motivate a great number of similar studies @xcite . especially in astrophysical applications , constraints coming out of some supernova explosion scenarios @xcite , may be affected and eventually lead to the necessity of further investigation of nsi in both lfv and clfv processes that may occur in solar and supernova environment @xcite .
such open issues motivated our present work too .
one of our main purposes in this paper , which is an extension of our previous study @xcite , is to comprehensively study the above issues by performing nuclear structure calculations for a set of experimentally interesting nuclei .
we estimate reliably the nuclear matrix elements describing both interaction channels , the exotic and the standard model ones , but we mainly focus on the sm component of the neutrino - nucleus processes , i.e. we consider @xmath10 in the reactions of eq .
( [ neutrin - nsi ] ) .
exotic neutrino - nucleus events are also computed . by exploiting our accurate original cross sections ,
we obtain simulated @xmath0-signals and flux averaged cross sections which are experimentally interesting quantities for both supernova and sns neutrinos .
the total number of events expected to be recorded over the energy - threshold for the studied nuclear targets are also presented for both cases .
we stress that , we have devoted special effort to obtain results of high accuracy by constructing the nuclear ground state within the context of the quasi - particle random phase approximation ( qrpa ) , i.e. by solving iteratively the bcs equations for realistic pairing interactions ( the bonn c - d potential ) @xcite , and achieving high reproducibility of the available experimental data @xcite .
in addition , we made comparisons with the results of other methods evaluating the nuclear form factors that enter the coherent rate @xcite as the one which employs fractional occupation probabilities ( fop ) of the states ( on the basis of analytic expressions ) @xcite , and other well - known methods @xcite .
in this section we present briefly the necessary formalism for describing all channels of the nsi processes of the reactions ( [ neutrin - nsi ] ) , derived by starting from the corresponding nuclear - level feynman diagrams .    -nucleus reactions , ( _ b _ ) non - standard z - exchange @xmath0-nucleus reactions , and ( _ c _ ) z - exchange and photon - exchange @xmath11 in the presence of a nucleus ( muon - to - electron conversion ) . the non - standard ( clfv or lfv ) physics enters in the complicated vertex denoted by the bullet @xmath12 @xcite.,scaledwidth=65.0% ]    in fig .
[ fig.1 ] the exchange of a @xmath13-boson between a lepton and a nucleon is represented , for the sm @xmath0-nucleus scattering , fig .
[ fig.1](a ) , and for the exotic @xmath0-nucleus scattering , fig .
[ fig.1](b ) .
as already mentioned in the introduction , the non - standard @xmath0-nucleus processes @xcite and the exotic clfv @xmath9 conversion in nuclei @xcite , can be predicted within the context of the same new - physics models @xcite .
for this reason , in fig .
[ fig.1](c ) we also show the exchange of a @xmath13-boson or a virtual @xmath14-photon leading to the nuclear @xmath9 conversion , @xcite .
thus , the leptonic vertex in the cases of fig . [ fig.1](b),(c ) is a complicated one .
a general effective lagrangian that involves sm interactions ( @xmath15 ) and nsi ( @xmath16 ) with a non - standard flavour preserving ( fp ) term , a non - universal ( nu ) term and a flavour changing ( fc ) term reads @xmath17 each of the components @xmath15 and @xmath16 , the individual terms @xmath18 and @xmath19 as well as the nuclear matrix elements that arise from each part , are discussed below .
the quark - level lagrangian for neutral current non - standard neutrino interactions @xmath16 , at the four fermion approximation , ( energies @xmath20 ) is parametrized as @xcite @xmath21\left[\bar{f}\gamma^\rho p f\right]\ , , \label{nsi_lagr}\ ] ] where @xmath22 denotes a first generation sm quark , @xmath23 are three light neutrinos with majorana masses and @xmath24 are the chiral projectors . in the latter lagrangian ( [ nsi_lagr ] ) , two classes of non - standard terms are considered ( i ) flavour preserving non - sm terms that are proportional to @xmath25 ( known as non - universal , nu interactions ) and ( ii ) flavour - changing ( fc ) terms proportional to @xmath26 , @xmath27 .
these couplings are defined with respect to the strength of the fermi coupling constant @xmath28 @xcite . in the present work
we examine spin - zero nuclei thus , the polar - vector couplings defined as @xmath29 are mainly of interest . for the axial - vector couplings it holds @xmath30 .    following ref .
@xcite , the nuclear physics aspects of the neutrino - matter nsi can be explored by transforming the quark - level lagrangian ( [ nsi_lagr ] ) , eventually to the nuclear level where the hadronic current is written in terms of nc nucleon form factors that are functions of the four momentum transfer . generally , for inelastic @xmath0-nucleus scattering , the magnitude of the three momentum transfer , @xmath31 , is a function of the scattering angle of the outgoing neutrino @xmath32 ( in laboratory frame ) , the initial , @xmath33 , and final , @xmath34 , nuclear energies , as well as the excitation energy of the target nucleus , @xmath35 , and takes the form @xmath36 @xcite .
our analysis in the present paper , concentrates on the dominant coherent ( elastic ) channel where only @xmath37 transitions occur ( @xmath38 , @xmath39 ) and the momentum transfer in terms of the incoming neutrino energy , @xmath40 , becomes @xmath41 or equivalently @xmath42 .
the nsi coherent differential cross section of neutrinos scattering off a spin - zero nucleus , with respect to the scattering angle @xmath32 reads @xcite @xmath43 @xmath44 , denotes the flavour of incident neutrinos and @xmath45 represents the nuclear ground state ( for even - even nuclei assumed here , @xmath46 ) .
the nuclear matrix element , that enters the cross section of eq .
( [ nsi_dcostheta ] ) , is written as @xcite @xmath47 ^ 2 \\ & \ , \ , \ , \ , + \sum_{\beta \neq \alpha } \left[\left ( 2 \epsilon_{\alpha \beta}^{uv}+ \epsilon_{\alpha \beta}^{dv }   \right ) z f_z ( q^2)+ \left(\epsilon_{\alpha \beta}^{uv}+ 2 \epsilon_{\alpha\beta}^{dv } \right)n f_n ( q^2)\right]^2 , \end{aligned } \label{gv}\ ] ] ( @xmath48 ) where @xmath49 denote the nuclear ( electromagnetic ) form factors for protons ( neutrons ) .
we stress that , in the adopted nsi model , the coherent nc @xmath0-nucleus cross section is not flavour blind as in the sm case . obviously , by incorporating the nuclear structure details , in eqs .
( [ nsi_dcostheta ] ) and ( [ gv ] ) , the cross sections become more realistic and accurate @xcite .
the structure of the lagrangian ( [ tot_lagr ] ) , implies that in the r.h.s .
( [ gv ] ) the first term is the nu matrix element , @xmath50 , and the summation is the fc matrix element , @xmath50 , hence we write @xmath51    from experimental physics perspectives , it is rather crucial to express the differential cross section with respect to the recoil energy of the nuclear target , @xmath52 . in recent years it became feasible for terrestrial neutrino detectors to detect neutrino events by measuring nuclear recoil @xcite .
therefore , it is important to compute also the differential cross sections @xmath53 . in the coherent process ,
the nucleus recoils ( intrinsically it remains unchanged ) with energy which , in the approximation @xmath54 takes the maximum value @xmath55 , with @xmath56 denoting the nuclear mass @xcite .
then , to a good approximation , the square of the three momentum transfer , is equal to @xmath57 , and the coherent nsi differential cross section with respect to @xmath58 can be cast in the form @xmath59 we note that , compared to previous studies @xcite , we have also taken into consideration the interaction @xmath60 quark [ see eq .
( [ gv ] ) ] , in addition to the momentum dependence of the nuclear form factors @xcite . both eqs .
( [ nsi_dcostheta ] ) and ( [ nsi_dt ] ) are useful for studying the nuclear physics of nsi of neutrinos with matter
.    furthermore , by performing numerical integrations to eq .
( [ nsi_dcostheta ] ) over the scattering angle @xmath32 or to eq .
( [ nsi_dt ] ) over the recoil energy @xmath52 , one can obtain integrated ( total ) coherent nsi cross sections , @xmath61 .
following eq .
( [ nu - fc - terms ] ) , the individual cross sections @xmath62 and @xmath63 may be evaluated accordingly .      the effective ( quark - level ) sm @xmath0-nucleus interaction lagrangian ,
@xmath15 at low and intermediate neutrino energies , is written as @xmath64 \left [ \bar{f } \gamma^{\rho } p f \right ] \ , , \label{sm_lagr}\ ] ] where @xmath65 and @xmath66 are the left- and right - handed couplings of the @xmath67-quark to the @xmath13-boson and @xmath68 and @xmath69 are the corresponding couplings of the @xmath70-quark ( @xmath71 is the weinberg mixing angle ) .    for coherent @xmath0-nucleus scattering ,
the sm angle - differential cross section reads @xmath72 the operator @xmath73 in the nuclear matrix element of the latter equation is the coulomb operator which is equal to the product of the zero - order spherical bessel function times the zero - order spherical harmonic @xcite .
this matrix element can be cast in the form @xcite
@xmath74 ^ 2 \ , , \label{sm - me}\ ] ] where , the polar - vector couplings of protons @xmath75 and neutrons @xmath76 with the @xmath13 boson ( see fig . [ fig.1](a ) ) , are written as @xmath77 and @xmath78 , respectively .
as can be easily seen , the vector contribution of all protons is very small ( @xmath79 ) , hence the coherence in eq .
( [ sm - me ] ) essentially refers to all neutrons only of the studied nucleus .
after some straightforward elaboration the differential cross section with respect to the nuclear recoil energy , @xmath52 , takes the form @xmath80 the lagrangian @xmath81 of eq .
( [ tot_lagr ] ) , contains the flavour preserving ( fp ) part , equal to @xmath82 , which can be evaluated through the coulomb matrix element @xmath83 subsequently , the total coherent cross section may be computed on the basis of the matrix element @xmath84    in a previous work @xcite we evaluated original differential cross sections @xmath85 and @xmath86 , as well as individual angle - integrated cross sections of the form @xmath87 , with @xmath44 , and @xmath88 ( fc stands for the six flavour changing processes @xmath89 ) .    in this work ,
we perform standard model cross sections calculations ( for convenience from now on we drop the index @xmath90 and always consider @xmath91 ) for a set of nuclei throughout the periodic table up to @xmath92pb .
we adopt various nuclear models ( see section [ nucl - methods ] ) to compute the nuclear form factors .
then , for a great part of the cross section results ( except differential cross sections ) we evaluate folded cross sections , and event rates .
at first , we study the nuclear structure details of the matrix elements entering eq .
( [ sm - me ] ) , such results reflect the dependence of the coherent cross section on the incident - neutrino energy @xmath93 and the scattering angle @xmath32 ( or the recoil energy @xmath58 ) .
we mention that for the even - even nuclei this study involves realistic qrpa calculations for the differential cross sections @xmath94 and @xmath95 , performed after constructing the nuclear ground state @xmath96 by solving iteratively the bardeen cooper schriffer ( bcs ) equations .
the solution of these equations provides the probability amplitudes @xmath97 and @xmath98 of the @xmath99-th single nucleon level to be occupied or unoccupied , respectively .
moreover , the latter equations provide the single quasi - particle energies , based on the single particle energies of the nuclear field ( a coulomb corrected woods - saxon potential in our case ) as well as the pairing part of the residual two - body interaction ( bonn c - d potential in our case ) .
then , the nuclear form factors for protons ( neutrons ) are obtained as @xcite @xmath100\ , \langle j\vert j_0(qr)\vert j\rangle\left(\upsilon^j_{n_n}\right)^2 \ , , \ ] ] with @xmath101=\sqrt{2 j + 1}$ ] , @xmath102 .
for each nuclear system studied , the chosen active model space , the harmonic oscillator ( h.o . ) parameter @xmath103 and the values of the two parameters @xmath104 for proton ( neutron ) pairs that renormalise the monopole ( pairing ) residual interaction ( obtained from the bonn c - d two - body potential describing the strong two - nucleon forces ) , are presented in table [ table1 ] .
the adjustment of @xmath104 is achieved through the reproducibility of the pairing gaps @xmath105 ( see e.g. @xcite ) .
c|cccccc nucleus & model - space & @xmath103 & @xmath106 & @xmath107 & @xmath108 & @xmath109 + @xmath110 & 8 ( no core ) & 1.522 & 4.68536 & 4.84431 & 1.12890 & 1.19648 + @xmath111 & 8 ( no core ) & 1.675 & 3.36181 & 3.49040 & 1.06981 & 1.13636 + @xmath112 & 10 ( no core ) & 1.727 & 3.81516 & 3.83313 & 1.15397 & 1.27600 + @xmath113 & 10 ( no core ) & 1.809 & 3.03777 & 3.14277 & 1.15568 & 1.23135 + @xmath114 & 15 ( no core ) & 1.902 & 1.75518 & 1.76002 & 0.94388 & 1.01348 + @xmath115 & 15 ( no core ) & 1.952 & 1.91109 & 1.55733 & 1.05640 & 0.99890 + @xmath116 & 15 ( no core ) & 2.086 & 1.52130 & 1.56935 & 0.95166 & 1.17774 + @xmath117 & 15 ( core @xmath118 ) & 2.262 & 1.19766 & 1.20823 & 0.98207 & 1.13370 +      the nuclear form factor , which is the fourier transform of the nuclear charge density distribution @xmath119 , is defined as @xmath120 with @xmath121 being the zero order spherical bessel function . due to the significance of the nuclear form factors in our calculations and for the benefit of the reader
we devote a separate discussion to summarise some useful possibilities of obtaining these observables .
\i ) _ use of available experimental data _
+ for many nuclei and especially for odd - a isotopes , the proton nuclear form factors @xmath122 , are computed by means of a model independent analysis ( using a fourier - bessel expansion model or others ) of the electron scattering data for the proton charge density @xmath119 @xcite wherever , possible .
the absence of similar data for neutron densities , restricts us to assume that @xmath123 . in this work
, we consider this method only for the case of the very heavy doubly closed @xmath124 nucleus .
ar and @xmath125ti nuclei .
the introduction of fractional occupation probabilities ( fop ) of the states provides higher reproducibility of the experimental data , compared to the simple shell - model and that of eq . (
[ f - bessel ] ) .
the bcs nuclear neutron form factor @xmath126 is also presented and compared . ]
\ii ) _ fractional occupation probabilities in a simple shell - model _
@xcite the form factor @xmath127 , for h.o .
wavefunctions has been written as @xcite @xmath128 the radial nuclear charge density distribution @xmath119 , entering the definition of eq .
( [ definition - ff ] ) , is written in the following compact form @xcite @xmath129 where @xmath130 , with @xmath103 denoting the h.o .
size parameter .
@xmath131 stands for the number of quanta of the highest occupied proton ( neutron ) level .
the coefficients @xmath132 are expressed as @xmath133 where @xmath134 is the gamma function .
for the coefficients @xmath135 , @xmath136 and further information see ref .
@xcite .
up to this point , that the proton occupation probabilities entering eq .
( [ definition - ff ] ) eq .
( [ forfac ] ) have been considered equal to unity for the states below the fermi surface and zero for those above the fermi surface . in ref .
@xcite , the authors introduced depletion and occupation numbers , to parametrise the partially occupied levels of the states .
these parameters satisfy the relation @xmath137 within this context , the `` active '' surface nucleons ( above or below the fermi level ) have non - zero occupation probability @xmath138 , smaller than unity , while the `` core '' levels have occupation probability @xmath139 . in this paper
we extend the work of ref .
@xcite where three parameters @xmath140 are used to describe the partial occupation probabilities of the surface orbits .
we improve the formalism by introducing more parameters , increasing this way the number of `` active '' nucleons in the studied nuclear system and come out with higher reproducibility of the experimental data @xcite . to this aim
, we introduce four parameters @xmath141 in eq .
( [ depletion ] ) .
then , the assumed `` active '' single - particle levels are five and eq .
( 16 ) of ref .
@xcite becomes @xmath142 \\ &
+ \pi(\chi , z_{c } ) \left[\frac{z^{\prime}-z}{z^{\prime}-z_{c}}-\frac{\alpha_{2}}{z_{c}-z_{1}}-\frac{\alpha_{3}}{z^{\prime}-z_{c } } \right ] \\ & + \pi(\chi , z^{\prime } ) \left[\frac{z - z_{c}}{z^{\prime}-z_{c } } + \frac{\alpha_{3}}{z^{\prime}-z_{c}}- \frac{\alpha_{4}}{z^{\prime \prime}-z^{\prime } } \right ] \\ & + \pi(\chi , z^{\prime \prime } ) \left[\frac{\alpha_{4}}{z^{\prime \prime}-z^{\prime } } - \frac{\lambda}{z^{\prime \prime \prime}-z^{\prime \prime } } \right ] + \pi(\chi , z^{\prime \prime \prime } ) \frac{\lambda}{z^{\prime \prime \prime}-z^{\prime \prime}}\ , , \end{aligned } \label{fractional}\ ] ] with @xmath143 . by substituting the polynomial @xmath144 of eq .
( [ rho(r ) ] ) with that of the latter expression and using the experimental data @xcite , we we fit the parameters @xmath145 ( and similarly for the form factor of eq .
( [ forfac ] ) ) . as an example , for the @xmath3ar isotope we have , @xmath146 .
the resulting fractional occupation probabilities that fit the experimental charge density distribution are @xmath147 , @xmath148 , @xmath149 , @xmath150 .
similarly for the @xmath125ti nucleus , we have @xmath151 and the fitting parameters are @xmath152 , @xmath153 , @xmath154 , @xmath155 . in fig .
[ fig.2 ] the prediction of the method is compared with that of the simple shell - model and the experimental data .
we note that in the momentum transfer range of our interest ( i.e. @xmath156 ) the form factor has excellent behaviour .
we however mention that even though the fop method presents very high reproducibility of the experimental data , it is not always applicable , e.g. for deformed nuclei ( where bcs appears to be still successful ) .
\iii ) _ use of effective expressions for the nuclear form factors _ +
we finally discuss one of the most accurate effective methods for calculating the nuclear form factor by ref .
@xcite @xmath157\ , , \label{f - bessel}\ ] ] where @xmath158 is the known first - order spherical - bessel function and @xmath159 , with @xmath160 and @xmath161 being the radius and surface thickness parameters of the nucleus respectively .
the radius parameter is usually given from the semi - empirical form @xmath162 while @xmath161 is of the order of @xmath163 ( see ref .
@xcite ) .
it is worth noting that , by inserting the form factors @xmath49 obtained as described above in eq .
( [ sm - me ] ) , the resulting cross sections have a rather high confidence level . in the next part of the paper the results show that the momentum dependence of the nuclear form factors becomes crucial , especially for intermediate and high energies . in some cases , differences of even an order of magnitude may occur as compared to the calculations neglecting the momentum dependence of the nuclear form factors .
the next phase of our calculational procedure is related to the total coherent @xmath0-nucleus cross sections , obtained through numerical integration of eq .
( [ sm_dcostheta ] ) over angles [ or eq .
( [ sm_dt ] ) over @xmath52 ] as @xmath165 the results for the standard model cross sections , for a set of different promising targets throughout the periodic table , are presented in fig .
[ fig.3 ] . as can be seen
, the present nuclear structure calculations indicate that between light and heavy nuclear systems , the cross sections may differ by even two orders of magnitude ( or more ) as a consequence of the dependence on the nuclear parameters ( i.e. mass , form factors , etc . ) .
we also see that for heavier nuclei the cross sections flatten more quickly ( at lower neutrino energies ) compared to that of lighter nuclear isotopes .
the latter conclusion originates mainly from the fact that , for heavy nuclei the suppression of the cross sections due to the nuclear form factors becomes more significant .
thus , for heavy material the nuclear effects become important even at low energies .
such original cross section results are helpful for the simulations of the standard and non - standard model signals of @xmath0-detection experiments ( see below ) .     in units
@xmath166 for a set of nuclei as a function of the incoming neutrino energy @xmath40 , for the sm neutrino processes @xmath167.,scaledwidth=50.0% ]      as discussed previously , our present calculations may also be useful for ongoing and future neutrino experiments related to supernova ( sn ) neutrino detection , since as it is known , the neutrinos emitted in sn explosions transfer the maximum part of the total the energy released .
then , the total neutrino flux , @xmath168 , arriving at a terrestrial detector as a function of the sn neutrino energy @xmath40 , the number of emitted ( anti)neutrinos @xmath169 at a distance @xmath70 from the source ( here we consider @xmath170 ) , reads @xcite @xmath171 ( @xmath172 ) where @xmath173 denotes the energy distribution of the ( anti)neutrino flavour @xmath174 .
the emitted sn - neutrino energy spectra @xmath175 may be parametrised by maxwell - boltzmann distributions that depend only on the temperature @xmath176 of the ( anti)neutrino flavour @xmath177 or @xmath178 ( the chemical potential is ignored ) we have @xmath179 ( @xmath180 , @xmath181 , @xmath182 , @xmath183 @xcite ) . for each flavour
, the total number of emitted neutrinos @xmath169 is obtained from the mean neutrino energy @xcite @xmath184 and the total energy released from a sn explosion , @xmath185 @xcite .
the spallation neutron source ( sns ) at oak ridge national lab @xcite produces neutrons by firing a pulsed proton beam at a liquid mercury target @xcite .
the main aim of the coherent proposal @xcite ( or of other similar conceptual @xcite ) concerns with possible detection of neutrino - nucleus coherent scattering events at the sns .
our simulations here , are mainly motivated by previous studies @xcite and the hope to provide our accurate nuclear structure calculations .    in stopped pion - muon sources ,
neutrinos are produced by the pion decay chain .
pion decay at rest @xmath186 ( @xmath187 ) produces monochromatic muon neutrinos @xmath188 at 29.9 mev , followed by electron neutrinos @xmath189 and muon antineutrinos @xmath190 that are produced by the muon - decay @xmath191 ( @xmath192 ) @xcite . for pulsed beams in time - scales narrower than @xmath193 , @xmath194 s and @xmath190 s
will be delayed with the beam while @xmath188 s will be prompt with the beam @xcite .
the emitted @xmath194 and @xmath195 neutrino spectra are described by the high precision normalized distributions , known as the michel spectrum @xcite @xmath196 ( @xmath197 is the muon rest mass ) .
the maximum neutrino energy in the latter distributions is @xmath198 ( see e.g. @xcite ) .
the spallation neutron source ( sns ) at oak ridge national lab is currently the most powerful facility to detect for a first time neutrino - nucleus coherent scattering events , since it provides exceptionally intense fluxes @xmath199 at 20 m and @xmath200 at 40 m from the source @xcite .
the simulated laboratory neutrino signals @xmath201 coming out of our calculations for the adopted nuclear targets are discussed below .      by weighting the integrated cross section @xmath202 with the neutrino distributions of eq .
( [ eta - sn ] ) , for sn neutrinos , or eq . ( [ labor - nu ] ) , for laboratory neutrinos , the total signal produced on a terrestrial detector is described by @xcite @xmath203 the resulting signals , @xmath204 , obtained by inserting in eq .
( [ eq.signal ] ) the cross sections @xmath205 of fig .
[ fig.3 ] , are plotted in fig .
[ fig.4 ] .
-detector , ( _ left _ ) for supernova neutrinos ( @xmath206 ) , evaluated with maxwell - boltzmann distributions at @xmath170 , and ( _ right _ ) for sns neutrinos ( @xmath207 ) , at 20 m from the source . for the case of sns
neutrinos the figure takes into account only the delayed beam , evaluated with the generic flux of @xmath208 .
different nuclear detectors have been studied . ]    in our previous work @xcite it was shown that the simulated cross sections reflect the characteristics of the incident neutrino spectrum of the specific neutrino flavour @xmath174 and therefore , such a simulated signal is characterised by its own position of the maximum peak and width of the distribution @xmath209 .
we , however , remind that within the framework of the sm , coherent neutrino scattering is a flavour blind and a particle - antiparticle blind process . for this particular case our results
are shown in fig .
[ fig.4 ] for supernova and laboratory ( sns ) neutrinos .    in neutrino simulations ,
another useful quantity is the flux averaged cross section @xcite which in our notation is written as @xmath210    c|cccccccccccc nucleus & @xmath110 & @xmath111 & @xmath112 & @xmath113 & @xmath114 & @xmath115 & @xmath211 & @xmath116 & @xmath117 & @xmath124 + @xmath212 & 1.46 & 2.51 & 3.91 & 7.52 & 18.59 & 25.43 & 33.29 & 70.63 & 207.56 & 514.93 + @xmath213 & 3.07 & 5.33 & 8.13 & 15.52 & 37.91 & 51.50 & 67.02 & 139.83 & 395.59 & 949.50 +    [ table2 ]    the results for @xmath214 , obtained by using the angle - integrated cross sections of fig .
[ fig.3 ] are listed in table [ table2 ] for both neutrino sources .      from experimental physics perspectives , predictions for the differential event rate , @xmath215 , of a @xmath0-detector
are crucial @xcite .
the usual expression for computing the yield in events is based on the neutrino flux , and is defined as @xcite @xmath216 where @xmath217 accounts for the total number of nuclei ( atoms ) in the detector material @xmath218 times the total time of exposure @xmath219 . using the latter equation ,
one concludes that , the lower the energy recoil , the larger the potentially detected number of events ( see fig .
[ fig.5 ] and fig .
[ fig.6 ] ) . in principle , in order to maximize the potential detection of a rare event process like the @xmath0-nucleus scattering , detector materials with very low energy - recoil threshold and low - background are required .
( _ right _ ) , for supernova neutrinos at @xmath170 . here ,
1 ton of perfectly efficient @xmath2ne and @xmath3ar detectors have been considered and also possible neutrino oscillation in propagation effects are neglected . for heavier nuclear targets
the differences become rather significant . in this figure
, @xmath220 stands for eq .
( [ f - bessel ] ) and fop for the method of fractional occupation probabilities of the states . for more details
see the text . ]
but for @xmath5ge and @xmath4xe .
]    in the last stage of our study we make predictions for the total number of coherent scattering events , the most important quantity , both from theoretical and experimental perspectives . to this purpose
, we evaluate the number of expected counts , for the studied detector materials , by performing numerical integration of eq .
( [ diff.rate ] ) over the nuclear recoil threshold @xmath221 ( see table [ table3 ] ) .
as has been discussed previously @xcite , sn neutrino detection might become possible by the massive dark matter detectors @xcite which have very good energy resolution and low threshold capabilities @xcite .
these experiments are designed ( or planned ) to search for wimps @xcite and/or other rare events such as the neutrinoless double beta decay . the latter ,
use heavy nuclei as nuclear detectors , e.g. ge ( gerda @xcite and supercdms @xcite experiments ) .
in addition we report that sn neutrino events can be potentially detected by experiments using noble gases like ne ( clean detector @xcite ) , ar ( warp programme @xcite ) and xe ( xenon 100 collaboration @xcite ) .
c|cccccc nucleus & @xmath58 & @xmath222 kev & @xmath223 kev & @xmath224 kev & @xmath225 kev + @xmath110 & 2.52 & 2.25 & 2.05 & 1.60 & 1.14 + @xmath111 & 3.29 & 2.84 & 2.51 & 1.83 & 1.19 + @xmath112 & 4.03 & 3.35 & 2.87 & 1.96 & 1.16 + @xmath114 & 9.46 & 6.63 & 5.01 & 2.53 & 1.00 + @xmath115 & 10.73 & 7.04 & 5.06 & 2.27 & 0.76 + @xmath211 & 12.00 & 7.36 & 5.04 & 2.01 & 0.57 + @xmath116 & 18.58 & 9.61 & 5.82 & 1.70 & 0.30 + @xmath117 & 30.68 & 9.84 & 4.16 & 0.46 & 0.01 + @xmath124 & 46.93 & 7.86 & 1.95 & 0.03 & @xmath226 +    [ table3 ]    as mentioned in section [ nucl - methods ] , in order to test our nuclear calculations we have also employed other nuclear methods . to this purpose , we have compared our original results evaluated with the bcs method with those obtained as discussed in subsection [ other - methods ] and concluded that for the case of the coherent channel all available nuclear methods are in good agreement , but their results differ significantly from those obtained assuming @xmath227 ( see fig . [ fig.5 ] and
[ fig.6 ] ) .
we stress however , that , since the cross section is mostly sensitive to the neutron distribution of the target nucleus , the most accurate method ( at low and intermediate energies ) is the bcs method which provides realistic proton as well as neutron form factors .
all other methods employed here consider only the proton distribution and assume @xmath228 , which especially for heavy nuclei , is a rather crude approximation .
we remark , however , that the aforementioned nuclear methods offer reliable results on the differential and total event rates for low energies ( see fig .
[ fig.5 ] and fig .
[ fig.6 ] ) , but in order to correctly estimate the neutron form factor , methods like the bsc are probably more appropriate . our present nuclear structure calculations for laboratory ( sns ) neutrinos @xcite ( see fig .
[ fig.7 ] ) , are in good agreement with previous results @xcite .
they imply that a comparably large number of coherent neutrino scattering events is expected to be measured by using lne , lar , lxe , ge and csi[na ] materials adopted by the coherent collaboration @xcite .
the predictions of the bcs method for these nuclei are illustrated in fig .
[ fig.7 ] and compared with those of other promising nuclear targets . because the neutrino flux produced at the sns is very high , ( of the order of @xmath229 per flavour at 20 m from the source @xcite ) , even kg - scale experiments expect to measure neutrino - nucleus coherent scattering events at significantly higher rates than those of supernova neutrinos .    it is worth noting that , the choice of the target nucleus plays also a crucial role , since a light nuclear target may yield almost constant number of events throughout the energy range , but small number of counts . on the other hand ,
a heavy nuclear target provides more counts , but yields low - energy recoil making the detection more difficult .
this leads to the conclusion that the most appropriate choice for a nuclear detector might be a combination of light and heavy nuclear isotopes , like the scintillation detectors discussed in ref .
@xcite .    ) .
the upper ( lower ) panel assumes the delayed ( prompt ) flux of laboratory stopped - pion neutrino sources .
this figure assumes a perfectly efficient detector and negligible neutrino oscillation effects.,scaledwidth=50.0% ]      the multi - target approach of the coherent experiment @xcite aiming on neutrino detection can also explore non - standard physics issues such as nsi @xcite , neutrino magnetic moment @xcite and sterile neutrino @xcite . in this subsection
we find it interesting to evaluate the non - standard neutrino - nucleus events that could be potentially detected by this experiment in each of the proposed nuclear targets .
the high intensity sns neutrino beams @xcite and the two promising @xmath0-detectors , liquid @xmath230 ( 391 kg ) and liquid @xmath114 ( 456 kg ) @xcite , firstly proposed by the clear @xcite and clean @xcite designs ( located at distance 20 m from the source ) , constitute excellent probes to search for the exotic @xmath0-reactions . other possibilities @xcite include medium and heavy weight targets like @xmath231 ( 100 kg ) inspired by the dark matter supercdms @xcite detector ( located at 20 m ) and @xmath232 ( 100 kg located at 40 m ) .     and ( _ right _ ) 456 kg of liquid @xmath114 , both located at a distance of 20 m ( @xmath233 ) from the source .
a perfectly efficient detector and negligible neutrino oscillation effects are assumed.,scaledwidth=100.0% ]    , but for 100 kg of @xmath5ge at 20 m ( @xmath234 ) and 100 kg of liquid @xmath4xe at 40 m ( @xmath235 ) from the source.,scaledwidth=100.0% ]    in fig .
[ fig.8 ] and fig .
[ fig.9 ] the resulting number of exotic events are illustrated and compared with the sm predictions .
we note , however , that especially for the case of the flavour changing ( fc ) channel @xmath236 , by using the extremely high sensitivity of the ongoing @xmath1 conversion experiments ( comet @xcite and mu2e @xcite ) , very robust bounds have been set on the vector parameters @xmath237 @xcite . to this end
, we conclude that , if the mu2e and comet experiments will not detect muon to electron conversion events , then the new @xmath237 parameters extracted in @xcite will lead to undetectable coherent rates at the sns facility for this channel .    for our present calculations we use the current bounds @xcite set by the sensitivity of the psi experiment @xcite and found countable number of events for the near detectors in the case of the corresponding @xmath236 reaction .
the other exotic parameters , i.e. @xmath238 with @xmath239 and @xmath240 have been taken from ref . @xcite . as discussed in ref .
@xcite , we do not take into account the @xmath241 contribution , since the corresponding limits are poorly constrained and eventually predict unacceptably high rates .    before closing
, it is worth noting that , the present calculations indicate significant possibility of detecting exotic neutrino - nucleus events through coherent scattering in the aforementioned experiments .
since neutrino - physics enters a precision era @xcite , a difference from the standard model predictions leads to an undoubtable evidence of non - standard neutrino - nucleus interactions ( nsi ) .
we recall that , in order to experimentally constrain simultaneously all the exotic parameters at high precision , the detector material should consist of maximally different ratio @xmath242 @xcite .
our future plans include estimation of the incoherent channel which may provide a significant part of the total cross section , especially for energies higher than @xmath243 mev ( depending on the nuclear target @xcite and the particle model predicting the exotic process ) .
initially , in this paper the evaluation of all required nuclear matrix elements , related to standard model and exotic neutral - current @xmath0-nucleus processes is formulated , and realistic nuclear structure calculations of @xmath0-nucleus cross sections for a set of interesting nuclear targets are performed .
the first stage involves cross sections calculations for the dominant coherent channel in the range of incoming neutrino - energies @xmath244 mev ( it includes @xmath0-energies of stopped pion - muon neutrino decay sources , supernova neutrinos , etc ) .
additionally , new results for the total number of events expected to be observed in one ton of various @xmath0-detector materials are provided and the potentiality of detecting supernova as well as laboratory neutrino - nucleus events is in detail explored .
the calculations are concentrated on interesting nuclei , like @xmath112 and @xmath114 , @xmath5ge and @xmath4xe which are important detector materials for several rare event experiments , like the coherent at oak ridge national laboratory , and also experiments searching for dark matter events as the gerda , supercdms , xenon 100 , clean , etc . by comparing our results with those of other methods
, we see that the nuclear physics aspects ( reflecting the accuracy of the required @xmath0-nucleus cross sections ) , appreciably affect the coherent @xmath37 transition rate , a result especially useful for supernova @xmath0-detection probes .    in the present work ,
the qrpa method that considers realistic nuclear forces has been adopted in evaluating the nuclear form factors , for both categories of @xmath0-nucleus processes , the conventional and the exotic ones .
also , a comparison with other simpler methods as ( i ) effective methods and ( ii ) the method of fractional occupation probabilities , which improves over the simple shell - model and gives higher reproducibility of the available experimental data , is presented and discussed .
we conclude that among all the adopted methods the agreement is quite good , especially for light and medium nuclear isotopes .
however , since coherent neutrino - nucleus scattering can probe the neutron nuclear form factors , methods like the bcs provide more reliable results .    in view of the operation of extremely intensive neutrino fluxes ( at the sns , psi , j - parc , fermilab , etc .
) , the sensitivity to search for new physics will be largely increased , and therefore , through coherent neutrino - nucleus scattering cross section measurements , several open questions ( involving non - standard neutrino interactions , neutrino magnetic moment , sterile neutrino searches and others ) may be answered . towards this purpose ,
we have comprehensively studied the non - standard neutrino - nucleus processes and provided results for interesting nuclear detectors .
our predictions for the total number of events indicate that , within the current limits of the respective flavour violating parameters , the coherent experiment may come out with promising results on nsi
. moreover , this experiment in conjunction with the designed sensitive muon - to - electron conversion experiments ( mu2e , comet ) , may offer significant contribution for understanding the fundamental nature of electroweak interactions in the leptonic sector and for constraining the parameters of beyond the sm lagrangians .
one of us , ( dkp ) wishes to thank dr .
kosmas for technical assistance .
00 d.z . freedman , phys .
d 9 ( 1974 ) 1389 .
a. drukier and l. stodolsky , phys . rev .
d 30 ( 1984 ) 2295 .
k. langanke , acta phys . pol .
b 39 ( 2008 ) 265 .
h. ejiri , phys .
338 ( 2000 ) 265 .
kosmas and e. oset , phys .
c 53 ( 1996 ) 1409 . j. schechter and j.w.f .
valle , phys .
d 22 ( 1980 ) 2227 ; + j. schechter and j.w.f .
valle , phys .
d 25 ( 1982 ) 774 .
k. scholberg , phys .
d 73 ( 2006 ) 033005 . a.a .
aguilar - arevalo et al . , phys . rev .
d 79 ( 2009 ) 072002 ; arxiv0806.1449 .
louis , prog .
63 ( 2009 ) 51 - 73 .
a. bolozdynya et al . , arxiv:1211.5199 .
d. akimov et al .
( csi collaboration ) , arxiv:1310.0125 .
k. scholberg et al . , arxiv:0910.1989 . j.
measuring coherent - nc@xmath0as at fermilab . " ( 2011 ) .
url : http://if-neutrino.fnal.gov/neutrino1-pagers.pdf .
s.j . brice et .
al . , phys .
d 89 ( 2014 ) 072004 .
collar et al.,arxiv:1407.7524 , submitted to nucl .
k. hirata , et al .
58 ( 1987 ) 1490 ; + r.m .
bionta , et al .
58 ( 1987 ) 1494 .
keil , g.g .
raffelt and h .-
janka , astrophys .
j. 590 ( 2003 ) 971 .
haxton , phys .
lett . 60 ( 1988 ) 768 .
giannaka , t.s .
kosmas , j. phys .
( 2013 ) 012124 .
avignone and y.v .
efremenko , j. phys .
g29 , ( 2003 ) 2615 .
y. efremenko and w. hix , j. phys .
( 2009 ) 012006 .
horowitz , k.j .
coakley , and d.n .
mckinsey , phys .
d 68 ( 2003 ) 023005 . a.j .
anderson , j.m .
conrad , e. figueroa - feliciano , k. scholberg and j. spitz , phys . rev .
d 84 ( 2011 ) 013008 .
kosmas and j.d .
vergados , phys .
d 55 ( 1997 ) 1752 .
vergados , t. s. kosmas , phys . atom .
61 ( 1998 ) 1066 .
p. brink et al . , arxiv : astro - ph/0503583v1 .
gerda collaboration , acta phys.polon , b41 ( 2010 ) 1469 .
mckinsey , k.j .
coakley astropart .
( 2005 ) 355 .
r. brunetti et al .
, new astron .
rev.49 ( 2005 ) 265 .
xenon 100 collaboration , arxiv:1107.2155v1 .
y. giomataris , j.d .
vergados , phys .
b 634 ( 2006 ) 23 .
j. monroe and p. fisher , phys .
d 76 ( 2007 )
033007 . m. biassoni and c. martinez , astropart .
36 ( 2012 ) 151 .
super - kamiokande collaboration , y. fukuda , et al . ,
86 ( 2001 ) 5651 .
sno collaboration , q.r .
ahmad , et al .
( 2002 ) 011302 .
kamland collaboration , k. eguchi , et al .
( 2003 ) 021802 , hep - ph/0212021 .
forero , m. tortola , j.w.f .
valle phys .
d 86 ( 2012 ) 073012 ; arxiv:1405.7540 [ hep - ph ] . j.w.f .
valle , nucl .
229 - 232 ( 2012 ) 23 - 29 and references therein .
cui et al .
( comet collaboration ) , kek report 2009 - 10 ; + comet collaboration , the comet proposal to jparc , jparc proposal , 2007 ; + comet collaboration , a. kurup , nucl .
suppl . 218
( 2011 ) 38 .
abrams et al .
( mu2e collaboration and project ) , mu2e conceptual design report , arxiv:1211.7019 ; + mu2e collaboration , f. cervelli , j. phys .
( 2011 ) 012073 .
bernstein , p.s .
cooper , phys . rept .
532 ( 2013 ) 27 .
barlow , nucl .
suppl . 218
( 2011 ) 44 .
y. kuno , nucl .
225 - 227 ( 2012 ) 228 - 231 .
y. kuno and y. okada , rev .
73 ( 2001 ) 151 .
s. davidson , c. pena - garay , n. rius , and a. santamaria , jhep 03 ( 2003 ) 011 .
j. barranco , o.g .
miranda and t.i .
rashba , jhep 0512 ( 2005 ) 021 .
papoulias and t.s .
kosmas , phys .
b 728 ( 2014 ) 482 .
a. friedland , c. lunardini , and c. pena - garay , phys .
b 594 ( 2004 ) 347 .
a. friedland , c. lunardini and m. maltoni , phys .
d 70 ( 2004 ) 111301 .
a. friedland , c. lunardini , phys.rev .
d 72 ( 2005 ) 053009 . k.j .
healey , a.a .
petrov , d. zhuridov , phys.rev .
d87 ( 2013 ) 11 , 117301 .
anderson , j.m .
conrad , e. figueroa - feliciano , c. ignarra , g. karagiorgi , k. scholberg , m.h .
shaevitz , j. spitz , phys .
d 86 ( 2012 ) 013004 . p.s .
amanik and g.c .
mclaughlin , j.phys .
g 36 ( 2009 ) 015105 .
amanik , g.m .
fuller and b. grinstein , astropart .
( 2005 ) 160 .
a. esteban - pretel , r. tomas , and j.w.f .
valle , phys .
d 81 ( 2010 ) 063003 . d.k .
papoulias and t.s .
kosmas , j. phys .
410 ( 2013 ) 012123 .
chiang , e. oset , t.s .
kosmas , a. faessler and j.d .
vergados , nucl .
a 559 ( 1993 ) 526 .
kosmas and j.d .
vergados , phys .
( 1996 ) 251 .
kosmas , nucl .
phys . a 683 ( 2001 ) 443 ; + t.s .
kosmas , prog . part .
48 ( 2002 ) 307 .
kosmas , g.k .
leontaris and j.d .
vergados , prog . part .
nucl . phys .
33 ( 1994 ) 397 .
f. deppisch , t.s .
kosmas and j.w.f .
valle , nucl .
b 752 ( 2006 ) 80 .
forero , s. morisi , m. tortola , j.w.f .
valle , jhep 1109 ( 2011 ) 142 .
s.p das , f.f .
deppisch , o. kittel , j.w.f .
valle , phys .
d 86 ( 2012 ) 055006 .
petrov , d.v .
zhuridov , phys .
d 89 ( 2014 ) 033005 .
amanik and g.m .
fuller , phys .
d 75 ( 2007 ) 083008 .
j. barranco , o.g .
miranda , c.a .
moura and j.w.f .
valle , phys . rev .
d 73 ( 2006 ) 113001 . o.g .
miranda , m.a .
tortola and j.w.f .
valle , jhep 0610 ( 2006 ) 008 .
j. barranco , o.g .
miranda and t. i. rashba , phys .
d 76 ( 2007 ) 073008 . t.s .
kosmas and j.d .
vergados , phys .
b 215 ( 1988 ) 460 ; + t.s . kosmas and j.d .
vergados , nucl .
phys . a 510 ( 1990 ) 641 .
kosmas , j.d .
vergados , o. civitarese and a. faessler , nucl .
phys . a 570 ( 1994 ) 637 .
kosmas , s. kovalenko and i. schmidt , phys .
b 511 ( 2001 ) 203 ; + t.s .
kosmas , s. kovalenko and i. schmidt , phys .
b 519 ( 2001 ) 78 .
chasioti and t.s .
kosmas , nucl .
a 829 ( 2009 ) 234 .
v. tsakstara and t.s .
kosmas , phys .
c 83 ( 2011 ) 054612 ; + v. tsakstara and t.s .
kosmas , phys .
c 84 ( 2011 ) 064620 .
h. de vries , c.w . de jager and c. de vries , at . data and nucl .
data tables 36 ( 1987 ) 495536 .
donnelly and j.d.walecka , nucl .
phys . a 274 ( 1976 ) 368 .
alberico , et al .
c 79 ( 2009 ) 065204 .
kosmas and j.d .
vergados , nucl .
phys . a 536 ( 1992 ) 72 .
j. engel , phys .
b 264 ( 1991 ) 114 . c. dohmen et al .
( sindrum - ii collaboration ) , phys .
b 317 ( 1993 ) 631 .
p. wintz , in proceedings of the first international symposium on lepton and baryon number violation ( 1998 ) , edited by h.v .
klapdor - kleingrothaus and i.v .
krivosheina ( institute of physics publishing , bristol and philadelphia ) , p.534 . unpublished .