Simulated neutrino signals of low and intermediate energy neutrinos on Cd detectors

Neutrino-nucleus reactions cross sections, obtained for neutrino energies in the range $\varepsilon_{\nu}\leq 100-120$ MeV (low- and intermediate-energy range), which refer to promising neutrino detection targets of current terrestrial neutrino experiments, are presented and discussed. At first, we evaluated original cross sections for elastic scattering of neutrinos produced from various astrophysical and laboratory neutrino sources with the most abundant Cd isotopes $^{112}$Cd, $^{114}$Cd and $^{116}$Cd. These isotopes constitute the main material of the COBRA detector aiming to search for neutrinoless double beta decay events and neutrino-nucleus scattering events at the Gran Sasso laboratory (LNGS). The coherent $\nu$-nucleus reaction channel addressed with emphasis here, dominates the neutral current $\nu$-nucleus scattering, events of which have only recently been observed for a first time in the COHERENT experiment at Oak Ridge. Subsequently, simulated $\nu$-signals expected to be recorded at Cd detectors are derived through the application of modern simulation techniques and employment of reliable neutrino distributions of astrophysical $\nu$-sources (as the solar, supernova and Earth neutrinos), as well as laboratory neutrinos (like the reactor neutrinos, the neutrinos produced from pion-muon decay at rest and the $\beta$-beam neutrinos produced from the acceleration of radioactive isotopes at storage rings as e.g. at CERN).


I. INTRODUCTION
In recent interdisciplinary investigations in nuclear, particle and astro-particle physics, the interactions of neutrinos with matter play key role in understanding deeply the underlying physics. Exact measurements and reliable models of neutrino-matter interactions provide unquestionable requirements for unravelling top physics issues as neutrino properties, neutrino oscillations, supernova dynamics, dark matter detection and many others [1][2][3][4]. To enable further progress, relevant nuclear model calculations, across a wide energy range and in various nuclear isotopes, may provide significant results [5][6][7].
Recently, the neutrino-nucleon and neutrino-nucleus cross-section uncertainties have reached a limiting factor in the judgement of neutrino interaction models and in interpreting many neutrino experiments [1][2][3][4] and specifically experiments like the COHERENT where recently coherent neutrino nucleus scattering events have been measured for a first time [8][9][10]. Furthermore, the presence of important nuclear effects impact the interaction cross sections as well as the final nuclear states reached through the scattering process [11][12][13]. The nuclear effects also affect the rebuilding of the incoming ν-energy spectra of the neutrino sources that are key-role input for the resolution of neutrino detection signals. Understanding neutrino-nucleus scattering processes provides to experimentalists good information to separate the background events from the detection signal [14][15][16][17][18][19][20].
The current neutrino physics searches are categorized according to the incident neutrino energy in the scattering process. Thus, the range below about 10 MeV (low-energy, from a nuclear physics viewpoint) is connected to Geo-neutrino and solar neutrino studies [14][15][16][17][18][19], the neutrino energy range of 10 up to about 100-120 MeV (intermediate energy) covers a set of ν-physics topics in the front of nuclear structure physics and astro-particle physics such as core-collapse supernovae dynamics and dark matter detection [11-13, 21, 22], while the energy range from 0.1-0.2 GeV up to about 10 GeV is related to meson decay neutrino beams such as those employed for long-baseline (high energy) neutrino experiments [23][24][25][26].
Due to the fact that neutrinos interact very weakly, they are unique messengers from astrophysical sources (the Earth, the Sun, the supernovae and other stars) [1,2,21] allowing us to investigate deep into the astrophysical objects [21,[27][28][29]. In the near future, remarkably sensitive detectors as liquid-scintillator detectors, liquid argon time projection chambers and water-Cherenkov detectors would operate aiming to study neutrino physics issues of astrophysical neutrino sources [19,20,26] (for higher energy neutrinos, like e.g. those coming from active galactic nuclei, black hole binary stars, etc. operating detectors as IceCube, KM3Net and others are appropriate) [30,31]. Each detector type has specific advantages (e.g. for supernova neutrinos, a combination of all types may allow for a better investigation of the relevant open issues).
Our present work focuses on the interpretation of various ν signals generated in nuclear detectors of terrestrial neutrino experiments through the investigation of the nuclear response of Cd detector materials to various neutrino energy spectra [3,4,26,32]. We emphasize on signals coming from geo-neutrinos, solarneutrinos, supernova-neutrinos, reactor-neutrinos and neutrinos generated from the decay of stopped pions and muons.
(ii) Reliable descriptions of the shapes of neutrino energy distributions coming out of numerical simulations of distributions in neutrino-energies ε ν ≤ 100 − 120 MeV (for the above mentioned ν-sources). (iii) Modern computational tools [35][36][37][38][39] for the required folding (convolution) procedure in order to simulate the signal expected to be recorded on the Cd detectors CdTe or CdZnTe (the detector media of COBRA experiment) [3,4,26] from neutrino sources as the geo-, reactor-, solar-, supernovaand pion/muon decay neutrinos. We mention that, the response of the Cd isotopes in the particle-bound excitation region, which coincides with the energy range of geo-neutrinos, is rather rich and this is one of the motivations for performing the present calculations. The next generation detectors (LENA, Borexino, SNO+) [19,20], are expected to give useful answers to several questions of geological importance regarding the precise geo-ν fluxes and abundances of natural radioactive elements (K, U , T h) in the Earth's interior [40][41][42][43].
In this work we pay special attention on the coherent elastic neutrino-nucleus scattering (CEvNS) that is a process in which the target nucleus recoils coherently via a combined neutral current exchange width with neutrinos or anti-neutrinos. This process is well predicted by the standard model of the electroweak interactions and has large cross sections (10 −39 cm 2 in the neutrino-energy region (ε ν ≤50 MeV). This process has very recently been observed in the COHERENT experiment at a 6.7 σ confidence level (CL), by using a low-background CsI[Na] scintillator [8][9][10]. The detector was exposed to a ν µ neutrino beam coming from the Spallation Neutron Source (SNS) at Oak Ridge, USA [9]. This facility generates the most intense (pulsed) neutron beam in the world while simultaneously a significant yield of neutrinos is generated when pions (product of proton interactions in the target) decay at rest (prompt neutrinos). In addition, the muons produced from the charged-pion decay generate the known as delayed neutrino beam [10].
Even though many groups world-wide are now studying the difficult low-energy nuclear recoil signature, only a few sources, in specific nuclear reactors and spallation neutron sources yield the required neutrino-energy beams in adequate quantities for such measurements [44][45][46]. In our present theoretical work, we do not address the improved constraints derived from this dataset on nonstandard neutrino interactions with quarks (for a comprehensive discussion on this issue the reader is referred e.g. to Refs. [47,48] and references therein). The present article is an extension of our previous calculations performed in Ref. [11,13,32] and we used the same but slightly improved nuclear method. The extension refers to the employment of new detector isotopes and the better accuracy of the calculations [11][12][13]32].
In the rest of the paper, at first (Sections 2 and 3), the main formalism is described and original cross sections calculations are presented. Then (Sections 4 and 5), a description of the main characteristics of the low and intermediate energy neutrino sources addressed here are briefly summarized and folded cross sections as well as event rates for neutral current neutrino scattering off the 112 Cd, 114 Cd and 116 Cd isotopes are presented and discussed. Finally (Section 6), the main conclusions of the present work are extracted.

II. BRIEF DESCRIPTION OF THE FORMALISM
A. Angle differential coherent ν-nucleus cross section In the description of the ν-nucleus scattering, the angle differential cross section dσ/dΩ is a useful quantity. For the elastic-scattering of a neutrino with energy ε ν on a nucleus (A,Z) the angle differential cross section (with respect to the the scattering angle ϑ) is defined as [47,[49][50][51] (G F = 1.1664 × 10 −5 GeV −2 is the Fermi weak coupling constant). In this definition, the quantity F (q 2 ) includes the nuclear structure dependence of the cross section as [33,34] where Θ w denotes the weak mixing angle, known as Weinberg angle which takes the value sin 2 Θ w ≈ 0.2313. In Eqs. (1) and (2), Q w denotes the weak charge of the target nucleus given by The latter expression shows that, the neutron coherence of neutral currents (NC), in the case of neutron rich targets, provides large cross sections. This effect can be exploited in detecting, e.g. earth and sky neutrinos by measuring nuclear recoils. Measurements of these (NC) cross sections may also provide useful information about the neutrino source [21] and yield information about the primary neutrino fluxes, i.e. before flavour conversions in the neutrino sphere of core collapse supernovae.
The sensitivity of the coherent scattering channel to the neutron number in the target nucleus, may provide nuclear structure information through investigation of νnucleus scattering and the possibility to search for nonstandard neutrino physics by taking advantage of the flavour-blind nature of the process [48,52].
The ground-state elastic nuclear form factors, F Z (q 2 ) for protons and F N (q 2 ) for neutrons entering Eq. (2), are defined by and are normalized as F N,Z (q 2 = 0) = 1. In the latter equation, ρ n,p (r) denote the neutron (n) and proton (p) charge density distributions with j 0 (qr) = sin(qr)/(qr) being the zero-order spherical Bessel function (we neglect a small correction from the single-nucleon form factors proportional to e −(qbN ) 2 /6 with b N ≈ 0.8 fm being the nucleon harmonic oscillator size parameter [53]). The proton density ρ p (r) is often taken from experiment whenever measured charge densities are available [33,49]. Moreover, assuming that F N ≈ F Z , from Eqs. (1) and (2) (in nuclei with J π = 0 + ground state), one obtains where f p V and f n V stand for the polar-vector couplings of the weak neutral current Thus, the coherent cross section depends on the square of the ground-state nuclear form factor F (q 2 ) at momentum transfer q given by From Eq. (2), we see that, since f p V = (1 − 4 sin 2 Θ w )/2 ≈ 0.0374 is small, a neutrino scattered elastically on a spin-zero nucleus couples mostly to the neutron distribution, ρ n (r). A measurement of the cross section for this process would, at some level, provide a determination of the neutron form factor F N (q 2 ) [52,54]. Some authors consider that this would be complementary to parity violating experiments [50,52] because it would provide additional data, obtained at different energy ranges and with different nuclei that could be used to calibrate nuclear structure calculations [33,34,[49][50][51].
In earlier astrophysical estimations of the coherent scattering cross sections within the Standard Model (SM) [51,55] (also in recent beyond the SM calculations [56,57]), the approximation F N (q 2 ) ≈ F Z (q 2 ) ≈ 1 was used for the total coherent cross section σ tot (ε ν ) written as (we mention that available experimental data for neutron form factors are very limited).
From an experimental point of view, and particularly for the neutrino facilities near spallation sources [45,58], it is also interesting the expression of the coherent differential cross section as a function of the nuclear recoil energy T A . This is approximately written as [58][59][60][61] where M is the nuclear mass and F denotes the ground state elastic form factor of the target nucleus. For the sake of completeness, we note that other expressions, including higher order terms with respect to T A can be found, see e.g. Refs. [9,10,[58][59][60]. The contribution, however, of these therms is negligible and thus, higher order terms in Eq. (9) does not influence essentially the calculations. Our present coherent differential cross sections are not obtained as functions of the recoil energy but as functions of the scattering angle or the momentum transfer connected through Eq. (7). It should be noted that, the signal on the coherent neutrino-nucleus scattering experiments is significantly different compared to that of the incoherent scattering where the signal could be an outgoing particle or a deexcitation product [32]. The total coherent cross section σ tot (ε ν ) is obtained by integrating numerically Eq. (5) over the angle θ (θ min = 0 to θ max = π) or Eq. (9) over T A between 61,62]. Before closing this sections, it is worth mentioning that, in our present calculations of the neutrino-nucleus cross sections part of the cross-section uncertainties are removed by performing realistic nuclear structure calculations for both proton and neutron nuclear form factors (for a recent comprehensive discussion on this issue the reader is referred e.g. to Ref. [48] where the results coming out of different nuclear models and various approximations are presented and discussed).

III. ORIGINAL CROSS SECTION CALCULATIONS
The neutral-current scattering of low and intermediate energy neutrinos ν ℓ or anti-neutrinos ν ℓ (ℓ = e, µ, τ ) off the 112,114,116 Cd isotopes (with abundances 24.13%, 28.8% and 7.5%, respectively, the first two are the most abundant Cd isotopes) are represented by (Cd * denote excited states of Cd-isotopes). We mention that, the above reactions of the Cd-isotopes and also the charged-current (CC) reactions for ℓ = e, play significant role in astrophysical environment by affecting the electron fraction Y e of the matter and its strong effect on the matter flow [22,[63][64][65][66].
In the first step of the present calculations, we evaluate original cross sections for the coherent channel (ground state to ground state transitions) of the reactions of Eq. (10) [5,11,13,62,67,68]. As can be seen from Eq. (5), the original cross section for scattering of neutrinos, ν l or anti-neutrinos, ν l , are identical (this holds only for the coherent channel). The signal (folded cross section) on the nuclear detector, however, as we will see in Sections IV and V, could be significantly different. This is due to the flavour dependent energy distributions of the ν-beam reaching the nuclear detector, that enters in the folding procedure.
In this work, the required nuclear ground state wave functions are obtained from mean-field calculations using the successful Woods-Saxon interaction plus the monopole (pairing) interaction of the Bonn C-D potential. The ground state of the studied (even-even) 112,114,116 Cd isotopes (they have ground state spin |J πi i = |0 + gs ) is computed by solving iteratively the BCS equations [11,13,32,47,69].
In Table I, we list the values of the resulting pairing parameters (g p,n pair ) and the (theoretical) energy gaps (∆ th p,n ) for protons (p) and neutrons (n) determined at the BCS level for the above isotopes. As is well known, these parameters renormalise the pairing interaction of the Bonn C-D potential in order to fit the theoretical gaps, ∆ th p,n , to the empirical ones ∆ exp p,n . The latter are provided through the application of the three point formulas (see Appendix) by using the empirical separation energies (for protons and neutrons, S p,n ) of the neighbouring nuclear isotopes [13,32]. The values of the g p,n pair adjust reliably the empirical energy gaps (see Table I) [11,13,32,69,70].
The needed proton and neutron nuclear form factors in the context of QRPA are calculated from the expressions (V k j denotes the probability amplitude for proton or neutron occupancies of the single particle (nℓ)j-level). The summation, runs over the 15 active levels of the chosen model space (the same for proton and neutrons) as well as over the fully occupied j-levels for which V k j = 1 (they describe a 40 Ca closed core). The model space assumed consists of the major harmonic oscillator shells having quantum numbers N=3, 4, 5 (N=2n + ℓ).
In Fig. 1, the quantities needed for calculating the differential and integrated coherent cross section (see Eqs. (1) and (5)) for the neutrino reactions (10) are illustrated. Figure 1(a)-(c), shows the form factors for protons (F Z ) and neutrons (F N ) obtained with our BCS calculations (for the three isotopes 112,114,116 Cd) and Fig. 1(d) shows the momentum dependence of F (q 2 ) that enters Eqs. (1) and (5).
It should be noted that, the corrections due to the nucleon finite size (e −(qbN ) 2 /6 ) and the nuclear centerof-mass motion (e (qb) 2 /4A ), which enter as an overall qdependent factor in the F N,Z (q), for the medium heavy Cd-isotopes are negligible and have been ignored. The correction due to the nucleon finite size (the larger of the two) is very well known, but not essential. For small q the influence is close to zero while at the maximum momentum q it is about 5% [53].
As can be concluded from Fig. 1, the above ground state properties of the three Cd isotopes studied are to a large extent similar which means that their nuclear structures are not significantly different (all of them have ground state spin J π = 0 + ). The differences, are mostly due to the small ratio (∆N i /N ≈ 3% − 6%) in their neutron number. Figure 2 illustrates the total integrated coherent cross sections of ν-112,114,116 Cd scattering as a function of (i) the momentum transfer q, Fig. 2(a), and (ii) the incoming neutrino energy ε ν , Fig. 2(b). As mentioned before, these original cross sections will be used below for evaluations of flux averaged folded cross sections for various neutrino spectra.
Before closing this section, it is worth mentioning that, in calculating the nuclear form factors F (q 2 ), see Fig. 1, in the context of the QRPA method, the estimated er-    Table. ror at low momentum transfer is very small, while in the momentum range of our interest 0 ≤ q ≤ 2f m −1 , it is at maximum 10-15%. On the other hand, the experimental accuracy, for the proton form factors entering Eq. (2), usually they come from electron scattering measurements, is of the order of 1% [50]. For neutron form factors, however, the available experimental data are limited and, in general, authors discuss about differences between corresponding proton and neutron nuclear form factors (in medium heavy isotopes like 112,114,115 Cd) of the order of 4 to 8% [52,54].
In the next section, we summarize the main features of the ν-energy distributions employed in this work for obtaining folded neutrino-nucleus cross sections for each ν-source.

IV. ENERGY-SPECTRA OF LOW-ENERGY AND INTERMEDIATE ν-SOURCES
In this section, we focus on the basic characteristics of the currently interesting astrophysical (solar-, supernova-, geo-neutrino) and laboratory (reactor neutrino and pion/muon decay at rest neutrino) sources, their energy spectra of which will be used in the convolution procedure (see next section) to obtain convoluted cross sections based on our original cross sections.
In general, the ν-beams of the above mentioned neutrino sources have broad energy distributions (sometimes they consist of a mixture of neutrinos and anti-neutrinos) characteristic of the considered source. Some well known mono-energetic (monochromatic) fluxes are e.g. the one coming out of the charged-pion decay at rest (corresponding to the energy ε νµ = 29.65 MeV, see Fig. 3(d panel). For the non-mono-energetic neutrino fluxes we define the energy distributions η(ε ν ) as (N ν represents the number of neutrinos of the beam). We note that, via these energy spectra η(ε ν ) of the specific neutrino sources, the original ν-nucleus cross sections (of neutral-current reactions) computed with the QRPA method, can be connected with physical observables and signals recorded at the nuclear detectors through the use of the folding (convolution) method described below. The obtained this way folded (convoluted) cross sections represent the simulated nuclear detector response of the 112,114,116 Cd isotopes, in the energy range of ν-energy distribution of the studied neutrino source.
The main properties of the aforementioned astrophysical and laboratory neutrinos are summarised in the next subsections.

A. Geoneutrinos
As it is well known, the decay of some radioactive isotopes (mainly U , T h, K) in the interior of our planet, makes the Earth a powerful source of low-energy neutrinos in the range ε ν 10 MeV [40][41][42][43]. Accurate measurements of the flux of these neutrinos [15,18] are utilized to determine the amount of heat-producing elements in the Earth's mantle. This amount may be compared to that estimated through indirect methods, an information which is important to understand the heat transfer within the Earth. The latter is responsible for earthquakes and volcanoes. The most recent measurements from Kam-LAND and Borexino [14,16,17] are useful to put limits on the parameters of various models describing the structure and evolution of our planet.
The Earth neutrinos (mainly electron anti-neutrinos ν e ), are generated through β-decay processes of neutronrich nuclei like U , T h and others. These thermonuclear reactions are accompanied by the emission of electrons (e − ) and release of energy Q β as [40] (A, Z) → (A, Z + 1) + e − + ν e + Q β .
A and Z denote the mass and atomic (proton) number, respectively, of the initial (parent) nucleus. Part of the decay energy Q β is carried away by anti-neutrinos (Q ν ) while the remainder is available for heating (Q h ). Thus, In general, the radioactive isotopes of the Earth are classified into three groups: (i) isotopes in the 238 U decay series, (ii) isotopes in 232 T h decay series, and (iii) 40 K isotope [40,41]. Thus, these isotopes are geologically important because they heat (radiogenic heat) the Earth's interior (finally each of them reaches a stable nuclear isotope) via β-decays of all intermediate radioactive isotopes. Figure 3(a panel), shows the individual anti-neutrino spectra from 40 K, 238 U series and 232 T h series (τ 1/2 = 4.47 × 10 9 y, τ 1/2 = 14.0 × 10 9 y and τ 1/2 = 1.28 × 10 9 y, respectively). Essentially, these anti-neutrino ( ν e ) energy spectra come from 82 beta decays in the U series and 70 beta decays in the T h series [40-43, 71, 72].

B. Solar neutrinos
The solar neutrino spectra (mainly ν e neutrinos) are produced through thermonuclear reactions taking place in the interior of the Sun [73][74][75]. The shape of the energy distribution (0.1 MeV≤ ε ν ≤ 18 MeV) depends on the densities and temperatures in the Sun's environment [74] and the individual process of the reaction chain (p-p neutrinos, 7 Be neutrinos, 8 B neutrinos, hep neutrinos, CNO-cycle neutrinos, etc.). In Fig.3(c panel), we show the energy spectra of the important 8 B [73] and hep [55,74] neutrino sources predicted by the standard so-lar model [55]. The 8 B ν-spectrum, is nearly symmetric, with a peak at 6.4 MeV while the hep spectrum is peaked at 9.6 MeV [55].
The detection of the solar neutrinos (produced either via the pp-chain reactions or via the CNO-cycle processes) by terrestrial experiments (SNO+ [17,19]), constitutes excellent probes for astrophysics, nuclear physics, and particle physics searches [74,75]. Besides the huge success of the solar-neutrino experiments the last decades, there are still many unsolved questions related to the metallicity of the Sun's core, the total luminosity in neutrinos, the neutrino oscillations, etc. [14,16,17,19,20].

C. Pion-muon decay at rest neutrino energy distributions
In muon factories (at J-Park, Fermilab, PSI, etc.), from pion and muon decay at rest (DAR), in addition to the monochromatic ν-beam peaked at ε νµ = 29.65 MeV), ν µ and ν e beams (with energy of a few tens of MeV) are created. Such intermediate energy neutrino sources, are also the currently available at high-intensity proton sources, like the SNS at Oak Ridge, the neutrino beam-line produced at Fermilab Booster, the future Project-X facilities at Fermilab, etc. [14,16,17,19,20].
In the farther future, such high-intensity muon beams would offer a possible site for neutrino experiments related to supernova neutrinos and for neutrino-nucleus cross section measurements in a great number of nuclei [25,[44][45][46]. In the operating pion-muon decay at rest neutrino sources (in Fermilab, at USA, J-PARC, at Japan, PSI in Switzerland, etc.) and in the neutrino facilities at the Neutron Spallation Source (Oak Ridge, USA), ν e neutrinos, andν µ anti-neutrinos are produced from the decay of muons according to the reaction The decaying muons result from the decay of pions at rest (π + → µ + + ν µ ). Thus, these neutrino beams are not completely pure as, for example, the β-beam neutrinos [23,25]. The energy-spectra of ν e andν µ neutrinos are fitted with the normalized distributions [71,76]  Obviously, the analytic expressions of Eqs. (15) and (16), are convenient for the required integrations in the folding procedure, see below [11,13,32,69]. On the other hand, their energy range and shape roughly resembles that of SN neutrinos.

D. Reactor Neutrino spectra
The fission of very heavy nuclear isotopes 235 U , 239 P u, and 238 U in the nuclear reactors produces a great number of neutron rich nuclear isotopes. Because these products are unstable, they decay via β-decay emitting antineutrinos ( ν e ) [77,78]. Hence, nuclear reactors, operate as intense ν e sources for many experiments, giving fluxes of the order of ∼ 10 13 ν/cm 2 s, at distances ∼ 10 m from the reactor core.
The energy spectrum of these anti-neutrinos, characteristic of the β − decay spectrum, is peaked at very low energies ∼ 0.3 MeV and covers the energy region below ∼ 10 MeV. Figure 3(b panel) illustrates the reactor neutrino spectra normalized so as the sum over all data-points to be equal to unity. The adopted fuel composition is 62% 235 U , 30% 239 P u and 8% 238 U [77,79].
Currently operating reactor neutrino experiments, like the TEXONO experiment in Taiwan [80,81], the MINER experiment at the Nuclear Science Center, Texas A&M University (using neutrinos from the TRIGA reactor) [82], are excellent probes of beyond the standard model neutrino physics searches (electromagnetic ν-properties) and coherent ν-nucleus scattering studies.

E. Supernova neutrino spectra
Supernovae (SN) play key role in the development of our Universe, indicated e.g from the fact that modern simulations of galaxies formation cannot reproduce the structure of the galactic disk without considering supernova data. Today, though the physics of core-collapse supernovae is not yet well-understood, investigations of SN neutrinos supply rich information for understanding their dynamics, the mechanism of SN-neutrino emission, etc., and for interpreting the supernova neutrino burst measurements [1,5,6]. Multiple physics signatures are expected from a core-collapse explosion in the next supernova observation [21,42,43,46]. The detection of a future galactic supernova will provide invaluable information on the astrophysics of core-collapse explosion while the high statistics of a galactic SN neutrino signal may allow us to unravel the relevant scenarios.
In general, the shape of SN-neutrino energydistributions is determined by the conditions pertaining during their emission from the collapsing star causing the cooling of the proto-neutron star formed in its center [63,[83][84][85][86]. For the energy distribution of SN neutrinos, some authors used available terrestrial neutrino sources with similar energy spectra, like the Neutron Spallation Source neutrinos and the boosted radioactive neutrino beams (beta beam neutrinos), in order to test the response of some ν-detectors to SN neutrinos [44][45][46]. Recent stellar modelling use analytic expressions that include various effects through a chemical potential parameter in the well-known two-parameter Fermi-Dirac (FD) distribution [87] or through the average ν-energy in the analytically simpler two parameter Power-Law (PL) distribution (see Appendix) [87][88][89].
Both parametrizations, FD and PL, yield similar distributions characterized by the temperature T or the average ν-energy ε ν [13,69,[90][91][92][93]. These analytic normalised expressions contain two parameters to include modulation effects due to various corrections required to modify the purely thermal shape initially employed [87,88,91]. The two parameter FD distribution includes the known pinching effect through the degeneracy parameter (the chemical potential divided by the neutrino temperature T), n dg = µ/T which makes the spectrum more narrow compared to the purely thermal shape of temperature T (in MeV) [13]. The two parameter PL distribution of SN-ν energy spectrum [88,89], contains as parameters the mean neutrino energy ε ν and the parameter α which adjusts the width w of the distribution [13,87,88,91] (see Appendix).
In Fig. 4, some flavour dependent ν-energy spectra (η SN ) emitted by a core-collapse Supernova, needed for our present work, are illustrated. Both FD and PL energy distributions (labelled η F D and η P L , respectively) are shown for three different values of the width parameter w = 0.7, w = 0.8 and w = 0.9 (see Appendix) and for five equivalent parametrizations. From the FD distributions (with parameters the temperature T and the width parameter w), we see that, as the temperature grows the maximum of the distribution shifts to larger ν-energy (at the same time the corresponding peak becomes smaller). Also, as the width parameter w grows (keeping the same temperature), both the maximum of the distribution shifts to smaller ε ν and its peak becomes smaller. Furthermore, the degeneracy parameter shifts the spectrum t o higher energies [13,91]. In this figure, the PL energy distributions for the corresponding values of mean neutrino energy ε ν , are also illustrated ( ε ν reflects the depth of the stars from which the neutrinos are escaping. We see that, as the ε ν grows, the maximum of the distribution shifts to higher ν-energy ε ν [88,91]. In Table II, the corresponding values of parameters for the equivalent FD and PL ν-energy spectra of Fig.  4, that have been employed in various SN scenarios are shown (for more details see the Appendix and Ref. [13]). It is worth mentioning that, due to neutrino oscillations and other phenomena, at any distance from the source the SN-ν spectra can be different compared to those originally produced at the core of the collapsing star. It is, however, expected that ν-signals with much higher statistics from future galactic SN, may allow us to assess the great number of neutrino mixing scenarios.
It is worth mentioning that, the statistics for the SN 1987A were rather poor, just a few dozen ν e events were received within about ten seconds. For the observation of the next core-collapse SN-neutrino burst, however, detectors with huge statistics and remarkably greater flavour sensitivity are in operation or have been planned to operate in the near future [94]. Among those, are the next generation detectors HyperKamiokande, Juno,   Dune, etc., which aim at measuring, among others, the diffuse SN neutrino background [95,96].

V. SIMULATED NEUTRINO SIGNALS ON NUCLEAR DETECTORS
The features of a neutrino-flux that arrives at a neutrino detector are concealed in the nuclear response of the detector-material. In the case of the COBRA detector, the semi-conductor materials CdTe or CdZnTe contain large portion of Cd isotopes [3,4,26]. Our aim in this section is to simulate some of these features by calculating convoluted cross sections as discussed in Refs. [13,32]. .
The convolution (folding) is carried out with (i) the original cross sections obtained in Sect. III, and (ii) the low and intermediate energy neutrino spectra of Section IV in order to compute, first, flux averaged total cross sections, σ tot and, then, corresponding supernova neutrino event rates and fluxes.

A. Flux averaged cross sections for Cd detector materials
For the coherent channel, which is possible only in neutral current neutrino-nucleus reactions studied in this work, the flux averaged cross section σ coh is obtained through the folding [6,11,13] For a CdTe or CdZnTe detector material, the flux averaged cross sections, computed by inserting in Eq. (17) the σ coh (ε ν ) from Fig. 2(b) and the η(ε ν ) from Figs. 3 and 4, for the isotopes 112,114,116 Cd, are listed in Tables III and IV as described below.
In Table III we list the flux averaged cross sections evaluated by adopting the neutrino distributions of the geo-neutrinos (see Fig. 3(a)), the reactor neutrinos (see Fig. 3(b)) and the solar neutrinos (see Fig. 3(c) for the 8 B and the hep solar neutrinos).
In the last two columns of Table IV we tabulate the σ coh calculated for the distributions of Eqs. (15) and (16), i.e. the ν-spectra produced by pion/muon decay at rest (DAR). In the first three columns of this Table, the flux averaged cross sections refer to various supernova neutrino scenarios described by the equivalent FD and PL distributions of Fig. 4. The corresponding parameters are listed in Table II. In supernova neutrino scenarios, usually average νenergies between 10 ≤ ε ν ≤ 12 MeV are employed for the description of ν e neutrinos, average energies between 15 ≤ ε ν ≤ 18 MeV for ν e anti-neutrinos, and average energies between 22 ≤ ε ν ≤ 26 MeV for ν x and ν x , with x = µ, τ [87][88][89]91].
Due to the dominance of the coherent channel throughout the region of the incoming neutrino energy ε ν of our present calculations, the flux averaged coherent cross section σ coh may be even two or three orders of magnitude larger than the total incoherent cross section σ incoh tot [11,13,69].

B. Number of events in ν-detectors
The present theoretical results may be connected with current neutrino experiments relying on Cd isotopes as detection materials, and specifically the COBRA experiment at Gran Sasso [3,4,26], as follows. By using the flux averaged cross sections σ(ε ν ) of Table IV, for instance those referred to the SN neutrinos of the 112,114,116 Cd isotopes, we estimate (potentially detectable) neutrino fluxes Φ ν that should arrive at each detector to create some typical scattering event rates N ev in the COBRA detector.
In general, the event rate N ev is related to the flux Φ ν reaching the nuclear detector with the expression [11,32,70]   We note that, experimentalists use the definition which takes into account the detection efficiency ǫ (usually equal to ǫ ≈ 80 -90%) of the specific detector. Here, we assume a COBRA detector of mass m det =100 kg and two cases of detector materials, i.e the semiconductors (a) CdZnTe and (b) CdTe [3,4,26].
Then, based on Eq. (18), we perform calculations assuming a total mass 100 kg of CdZnTe as COBRA detector which translates, for example, to approximately m Cd = 10.6 kg mass of 114 Cd isotope or equivalently a number of 114 Cd atoms (nuclei) equal to N Cd ≡ N114 Cd = 94.17N Avog .
In Eq. (18), as total neutrino scattering cross sections, σ tot (ε ν ) we employ the values of flux averaged cross sections σ coh of Table IV obtained through PL distribution for SN neutrino spectra (they refer to the three mean energies chosen above).
Finally, we choose four typical detection rates N ev as: In a similar way, assuming that the COBRA detector contains 100 kg of the material CdTe, we find 13.5 kg 114 Cd or about N Cd ≡ N114 Cd = 120.11N Avog atoms (nuclei) are contained in the second semiconductor material of COBRA detector. By performing similar calculations for the same SN scenarios and the same, as before, set of detection rates N ev , we find the corresponding fluxes Φ ν reaching the COBRA CdTe detector.
By performing the steps we followed for 114 Cd, for the other two Cd-isotopes, 112 Cd and 116 Cd, the resulting neutrino fluxes, for the chosen SN neutrino scenarios are listed in Table V (last four columns). Such results are useful for future use of the Cd materials in astrophysical neutrino detection. It should be stressed that, next generation experiments may be effective in the detection of much weaker signals (higher sensitivity, larger detector mass, etc.).
The above neutrino fluxes are of the same order with those of the Spallation Neutron Source (SNS) at OR-LaND, Oak Ridge [44][45][46]. We mention that the CO-HERENT experiment at Oak Ridge, with a 14.57 kg of CsI scintillator detector, by using an SNS ν µ neutrino flux (coming from π-decay at rest) as high as Φ COH ν = 1.7 × 10 11 ν µ /cm 2 s, has measured 142 CEvNS events within a period of 308.1 live days (at a distance of L = 19.3 m from the source) [8]. These results translate to event rate N COH ev = 168/y ν µ neutrinos. From the results of Table V, we may define the ratio N ev /Φ ν for the COHERENT experiment (R COH ) and for a special ν µ neutrino case of the COBRA experiment (R COB ). For a comparison of these two experiments, we choose, for example, the results referred to the 112 Cd isotope of CdTe material of the COBRA detector (sixth line from the beginning of Table V refers to ν µ neutrinos). From these two ratios we find that R = R COH /R COB = 98.95/2.94 ≈ 34, which means that, for the chosen SN neutrino scenario, the COBRA detector may observe 12 ν µ /y only if its mass is equal to m ∼ 34 times larger than the assumed above 100 kg, i.e. only if the COBRA detector has a huge total mass m det = 3.4 t CdTe material (we mention that, in the assumed scenario, the SN ν µ neutrinos correspond a Temperature T=24 MeV, see one before last column of Table V). This example indicates also the corresponding cost for detector improvement so as to be able to record neutrino signals coming from interesting astrophysical sources.
We should finally note that, in this work the detection efficiency ǫ has not been considered (equivalently we assumed ǫ = 1). Also, the neutrino mixing has not been accounted for which means that we assumed the neutrino spectra arrived at the nuclear detector are described by PL distributions (as in stars interior) of the same values of the parameters.

VI. SUMMARY AND CONCLUSIONS
In this work, we present original neutrino-nucleus cross sections obtained with realistic nuclear structure calculations (use of the QRPA method) for scattering of low and intermediate energy neutrinos off the 112,114,116 Cd isotopes. These Cd-isotopes are contents (with large abundance) of the detector materials of the COBRA detector at Gran Sasso. The neutrino energy assumed covers currently interesting laboratory (reactor, pion/muon decay at rest neutrinos) and Astrophysical (solar, super-  nova and Earth) neutrino sources. Laboratory neutrino beams are important tools for studying standard and non-standard neutrino physics while astrophysical neutrinos are key particles in investigating the structure and evolution of stars as well to deepen our knowledge on the fundamental neutrino-nucleus interactions.
By utilizing the convolution procedure, we calculated flux averaged cross sections and event rates for the above ν-sources based on specific spectral distributions describing supernova neutrino energy spectra, solar neutrinos, geo-neutrinos and laboratory neutrinos as well as reactor neutrinos and pion-muon-stopped neutrinos. The fluxaveraged total coherent cross sections, σ coh , reflect the mean neutrino signals generated in several terrestrial detectors ( 112,114,116 Cd) from such ν-sources. Important connection of our present results with current experiments may also be achieved through the evaluation of the neutrino scattering event rates on Cd detectors.
The estimated neutrino fluxes and scattering event rates for Cd-isotopes, contents of the CdTe and CdZnTe materials of the COBRA detector at LNGS, may support this experiment to reach its goal in searching for neutrino observation and detection of rare events (double beta decay, etc).

VII. ACKNOWLEDGEMENTS
The present research was financially supported (V. Tsakstara and J. Sinatkas) by the Department of Informatics Engineering of the Technological Institute of Western Macedonia. Also, Dr. Odysseas Kosmas wishes to acknowledge the support of EPSRC via grand EP/N026136/1 "Geometric Mechanics of Solids".

B.
Normalization of the ν energy distributions η(εν) adopted in this work The distributions η(ε ν ) adopted in the present work (see Sect. IV), are considered to be normalized in such a way that ∞ 0 η(ε ν )dε ν = 1.
For example, in the case of η νe (ε ν ) of Fig. 3(d), the normalization gives The Fermi-Dirac (FD) and Power-law (PL) energy distribution are commonly used in Supernova neutrino parametrizations. Both the FD and PL yield very similar distributions characterized by the temperature T or the average energy ε ν and the width w of the spectrum is defined as w = ε 2 ν − ε ν 2 /w 0 where w 0 = ε ν / √ 3 is the width of the identical FD and PL distributions [13].

Fermi-Dirac (FD) energy distribution
By introducing the degeneracy parameter n dg (equal to the ratio of the chemical potential µ divided by the neutrino temperature T , i.e. n dg = µ/T ), the Fermi-Dirac energy distribution reads In this case, the width of the spectrum is reduced compared to the corresponding thermal spectrum (pinching effect). The normalization constant F 2 (n dg ) of this distribution depends on the degeneracy parameter n dg and is given by the relation Inserting Eq. (23) into Eq. (22), we take

Power-law energy distribution
The SN-neutrino energy spectra can be fitted by using a Power-Law energy distribution of the form [88] where ε ν is the neutrino mean energy. The parameter α adjusts the width of the spectrum (see text). The normalization factor C, is calculated from the normalization condition C ∞ 0 ε ν ε ν α e −(α+1)(εν / εν ) dε ν = 1 .