Impact Factor 2.638 | CiteScore 2.3
More on impact ›

Original Research ARTICLE

Front. Phys., 12 April 2019 |

Simulated Neutrino Signals of Low and Intermediate Energy Neutrinos on Cd Detectors

  • 1Department of Informatics Engineering, Technological Institute of Western Macedonia, Kastoria, Greece
  • 2Electrical Engineering Department, Technological Institute of Western Macedonia, School of Applied Science, Kozani, Greece
  • 3Division of Theoretical Physics, University of Ioannina, Ioannina, Greece
  • 4Modelling and Simulation Centre, MACE, University of Manchester, Manchester, United Kingdom

Neutrino-nucleus reactions cross sections, obtained for neutrino energies in the range εν ≤ 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 112Cd, 114Cd, and 116Cd. 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 ν-nucleus reaction channel addressed with emphasis here, dominates the neutral current ν-nucleus scattering, events of which have only recently been observed for a first time in the COHERENT experiment at Oak Ridge. Subsequently, simulated ν-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 ν-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 β-beam neutrinos produced from the acceleration of radioactive isotopes at storage rings as e.g., at CERN).

1. 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 unraveling top physics issues as neutrino properties, neutrino oscillations, supernova dynamics, dark matter detection and many others [14]. To enable further progress, relevant nuclear model calculations, across a wide energy range and in various nuclear isotopes, may provide significant results [57].

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 [14] and specifically experiments like the COHERENT where recently coherent neutrino nucleus scattering events have been measured for a first time [810]. Furthermore, the presence of important nuclear effects impact the interaction cross sections as well as the final nuclear states reached through the scattering process [1113]. 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 [1420].

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 [1419], 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 [1113, 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 [2326].

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, 2729]. 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, solar-neutrinos, supernova-neutrinos, reactor-neutrinos, and neutrinos generated from the decay of stopped pions and muons.

The main ingredients to this aim are: (i) The original differential and integrated cross sections of the neutral-current reactions of neutrinos, 112, 114, 116Cd(ν,ν′)112, 114, 116Cd*, and anti-neutrinos, 112, 114, 116Cd(ν~,ν~)112,114,116Cd*, computed for the coherent channel by using a refinement of the quasi particle random phase approximation (QRPA) [11, 13, 3234]. (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 [3539] 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-, supernova-, and 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, Th) in the Earth's interior [4043].

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−39cm2 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 [810]. 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 [4446]. In our present theoretical work, we do not address the improved constraints derived from this dataset on non-standard neutrino interactions with quarks (for a comprehensive discussion on this issue the reader is referred e.g., to Papoulias and Kosmas [47, 48] and references therein). The present article is an extension of our previous calculations performed in Tsakstara and Kosmas [11, 13] and Tsakstara [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 [1113, 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 112Cd, 114Cd, and 116Cd isotopes are presented and discussed. Finally (section 6), the main conclusions of the present work are extracted.

2. Brief Description of the Formalism

2.1. Angle Differential Coherent ν-Nucleus Cross Section

In the description of the ν-nucleus scattering, the angle differential cross section / 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 scattering angle ϑ) is defined as Papoulias and Kosmas [47], Engel [49], De Vries et al. [50], and Drukier and Stodolsky [51]

dσdΩ=GF24π2εν2(1+cosϑ)Qw24F(q2)2    (1)

(GF=1.1664×10-5 GeV−2 is the Fermi weak coupling constant). In this definition, the quantity F(q2) includes the nuclear structure dependence of the cross section as Kosmas et al. [33] and Kosmas [34]

F(q2)=1Qw[(14sin2Θw)ZFZ(q2)NFN(q2)]    (2)

where Θw denotes the weak mixing angle, known as Weinberg angle which takes the value sin2Θw0.2313. In Equations (1) and (2), Qw denotes the weak charge of the target nucleus given by

Qw=(1-4sin2Θw)Z-N.    (3)

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 flavor 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 non-standard neutrino physics by taking advantage of the flavor-blind nature of the process [48, 52].

The ground-state elastic nuclear form factors, FZ(q2) for protons and FN(q2) for neutrons entering Equation (2), are defined by

Fk(q2)=k4πj0(qr)ρn,p(r)r2dr,k=N,Z    (4)

and are normalized as FN,Z(q2=0)=1. In the latter equation, ρn,p(r) denote the neutron (n) and proton (p) charge density distributions with j0(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 bN ≈ 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 FNFZ, from Equations (1) and (2) (in nuclei with Jπ = 0+ ground state), one obtains

dσ(εν,ϑ)dcosϑ=GF22π(1+cosϑ)εν2[fVpZ+fVnN]2FZ2(q2).    (5)

where fVp and fVn stand for the polar-vector couplings of the weak neutral current

fVp=12-2sin2ΘW,  fVn=-12.    (6)

Thus, the coherent cross section depends on the square of the ground-state nuclear form factor F(q2) at momentum transfer q given by

q=2ενsin(ϑ/2),    (7)

From Equation (2), we see that, since fVp=(1-4sin2Θw)/20.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 FN(q2) [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, 4951].

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 FN(q2)FZ(q2)1 was used for the total coherent cross section σtotν) written as

σtot(εν)=GF28π[(14sin2ΘW)ZN]2εν2.    (8)

(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 TA. This is approximately written as Vergados and Giomataris [58], Vergados et al. [59], Giomataris and Vergados [60], and Vogel and Engel [61]

dσ(εν,TA)dTA=GF24πQW2M(1-MTA2εν2)F(2MTA2),    (9)

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 TA can be found, see e.g., [9, 10, 5860]. The contribution, however, of these therms is negligible and thus, higher order terms in Equation (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 Equation (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 de-excitation product [32].

The total coherent cross section σtotν) is obtained by integrating numerically Equation (5) over the angle θ (θmin = 0 to θmax = π) or Equation (9) over TA between


to TAmax= [47, 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 Papoulias and Kosmas [48] where the results coming out of different nuclear models and various approximations are presented and discussed).

3. Original Cross Section Calculations

The neutral-current scattering of low and intermediate energy neutrinos ν or anti-neutrinos ν~ (ℓ = e, μ, τ) off the 112, 114, 116Cd isotopes (with abundances 24.13, 28.8, and 7.5%, respectively, the first two are the most abundant Cd isotopes) are represented by

νl(ν~l)+112,114,116Cd112,114,116Cd*+νl(ν~l),    (10)

(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 Ye of the matter and its strong effect on the matter flow [22, 6366].

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 Equation (10) [5, 11, 13, 62, 67, 68]. As can be seen from Equation (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 4 and 5, could be significantly different. This is due to the flavor 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, 116Cd isotopes (they have ground state spin |Jiπi=|0gs+) is computed by solving iteratively the BCS equations [11, 13, 32, 47, 69].

In Table 1, we list the values of the resulting pairing parameters (gpairp,n) and the (theoretical) energy gaps (Δp,nth) for protons (p) and neutrons (n) determined at the BCS level for the above isotopes. As is well-known, these parameters renormalize the pairing interaction of the Bonn C-D potential in order to fit the theoretical gaps, Δp,nth, to the empirical ones Δp,nexp. The latter are provided through the application of the three point formulas (see Appendix) by using the empirical separation energies (for protons and neutrons, Sp, n) of the neighboring nuclear isotopes [13, 32]. The values of the gpairp,n adjust reliably the empirical energy gaps (see Table 1) [11, 13, 32, 69, 70].


Table 1. Pairing parameters gpairp (for protons), and gpairn (for neutrons) determining the monopole pairing interactions for each of the studied isotopes.

The needed proton and neutron nuclear form factors in the context of QRPA are calculated from the expressions

Fk(q2)=1kjĵ(n)j|j0(qr)|(n)j(Vjk)2,k=N,Z    (11)

(Vjk 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 Vjk=1 (they describe a 40Ca closed core). The model space assumed consists of the major harmonic oscillator shells having quantum numbers N = 3, 4, 5 (N = 2n + ℓ).

In Figure 1, the quantities needed for calculating the differential and integrated coherent cross section (see Equations 1, 5) for the neutrino reactions (10) are illustrated. Figures 1A–C, shows the form factors for protons (FZ) and neutrons (FN) obtained with our BCS calculations (for the three isotopes 112, 114, 116Cd) and Figure 1D shows the momentum dependence of F(q2) that enters Equations (1) and (5).


Figure 1. Neutron and proton nuclear form factors FN,Z(q2) (A) for 114Cd, (B) for 112Cd, and (C) for 116Cd isotopes. (D) The ground-state elastic nuclear form factor F(q2) for 112, 114, 116Cd isotopes.

It should be noted that, the corrections due to the nucleon finite size (e-(qbN)2/6) and the nuclear center-of-mass motion (e(qb)2/4A), which enter as an overall q-dependent factor in the FN, 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 Figure 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 (ΔNi/N ≈ 3 − 6%) in their neutron number.

Figure 2 illustrates the total integrated coherent cross sections of ν-112, 114, 116Cd scattering as a function of (i) the momentum transfer q, Figure 2A, and (ii) the incoming neutrino energy εν, Figure 2B. As mentioned before, these original cross sections will be used below for evaluations of flux averaged folded cross sections for various neutrino spectra.


Figure 2. Total cross sections of coherent (ground state to ground state g.s. → g.s.) transitions for the neutral current reactions 112, 114, 116Cd(Cd(νl,νl)112,114,116Cd*, l = e, μ, τ, as a function of (A) the momentum transfer q, and (B) the incoming neutrino energy εν.

Before closing this section, it is worth mentioning that, in calculating the nuclear form factors F(q2), see Figure 1, in the context of the QRPA method, the estimated error at low momentum transfer is very small, while in the momentum range of our interest 0 ≤ q ≤ 2fm−1, it is at maximum 10–15%. On the other hand, the experimental accuracy, for the proton form factors entering Equation (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, 115Cd) 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.

4. 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 Figure 3D. For the non-mono-energetic neutrino fluxes we define the energy distributions η(εν) as

dNν(εν)dενη(εν)    (12)

(Nν represents the number of neutrinos of the beam).


Figure 3. (A) Spectra of the U-Series, Th-Series, and 40K Geo-Neutrinos. Neutrinos from 40K electron capture are also shown in this figure. (B) Normalized reactor neutrino spectra. (C) Normalized energy spectrum of 8B and hep νe solar neutrinos. (D) Energy-spectra of νe and ν~μ neutrino beams, generated from the muon-decay at rest (see e.g., [71, 72]).

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, 116Cd isotopes, in the energy range of ν-energy distribution of the studied neutrino source.

The main properties of the aforementioned astrophysical and laboratory neutrinos are summarized in the next subsections.

4.1. Geoneutrinos

As it is well-known, the decay of some radioactive isotopes (mainly U, Th, K) in the interior of our planet, makes the Earth a powerful source of low-energy neutrinos in the range εν ≼ 10 MeV [4043]. 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 KamLAND 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 neutron-rich nuclei like U, Th, 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β.    (13)

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 (Qh). Thus, Qβ = Qν + Qh.

In general, the radioactive isotopes of the Earth are classified into three groups: (i) isotopes in the 238U decay series, (ii) isotopes in 232Th decay series, and (iii) 40K 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 3A, shows the individual anti-neutrino spectra from 40K, 238U series, and 232Th series (τ1/2=4.47×109 y, τ1/2=14.0×109 y and τ1/2=1.28×109 y, respectively). Essentially, these anti-neutrino (ν~e) energy spectra come from 82 beta decays in the U series and 70 beta decays in the Th series [4043, 71, 73].

4.2. Solar Neutrinos

The solar neutrino spectra (mainly νe neutrinos) are produced through thermonuclear reactions taking place in the interior of the Sun [7476]. The shape of the energy distribution (0.1 MeV ≤ εν ≤ 18 MeV) depends on the densities and temperatures in the Sun's environment [75] and the individual process of the reaction chain (p-p neutrinos, 7Be neutrinos, 8B neutrinos, hep neutrinos, CNO-cycle neutrinos, etc.). In Figure 3C, we show the energy spectra of the important 8B [74] and hep [55, 75] neutrino sources predicted by the standard solar model [55]. The 8B ν-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 [75, 76]. 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].

4.3. 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, 4446]. 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

μ+e++νe+ν~μ.    (14)

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, 72]

ηνe(εν)=96εν2Mμ-4(Mμ-2εν),    (15)
ην~μ(εν)=16εν2Mμ-4(3Mμ-4εν),    (16)

see Figure 3D, where Mμ = 105.6 MeV, is the muon rest mass. The ν~μ spectrum is peaked at ενmax=52.8 MeV = Mμ/2 while that of νe is peaked at ενmax=35.2 MeV = Mμ/3 [7, 72].

Obviously, the analytic expressions of Equations (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.

4.4. Reactor Neutrino Spectra

The fission of very heavy nuclear isotopes 235U, 239Pu, and 238U in the nuclear reactors produces a great number of neutron rich nuclear isotopes. Because these products are unstable, they decay via β-decay emitting anti-neutrinos (ν~e) [77, 78]. Hence, nuclear reactors, operate as intense ν~e sources for many experiments, giving fluxes of the order of ~1013ν~/cm2 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 3B 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% 235U, 30% 239Pu, and 8% 238U [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.

4.5. 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 energy-distributions 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, 8386]. 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 [4446]. Recent stellar modeling 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) [8789].

Both parametrizations, FD and PL, yield similar distributions characterized by the temperature T or the average ν-energy 〈εν〉 [13, 69, 9093]. These analytic normalized 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), ndg = μ/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 Figure 4, some flavor dependent ν-energy spectra (ηSN) emitted by a core-collapse Supernova, needed for our present work, are illustrated. Both FD and PL energy distributions (labeled ηFD and ηPL, 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].


Figure 4. Supernova neutrino energy spectra (ηSN) coming out of the analytic expressions of: (i) the two-parameter Fermi-Dirac distribution (FD) and (ii) the two-parameter Power-law (PL) distribution (see Appendix). The five sets of values of their parameters refer to equivalent distributions (for details see the text).

In Table 2, the corresponding values of parameters for the equivalent FD and PL ν-energy spectra of Figure 4, that have been employed in various SN scenarios are shown (for more details see the Appendix and [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.


Table 2. Corresponding values of parameters for equivalent Fermi-Dirac (FD) and Power-Law (PL) distributions (SN neutrino energy spectra) of Figure 4.

It is worth mentioning that, the statistics for the SN 1987A were rather poor, just a few dozen ν~e events were received within about 10 s. For the observation of the next core-collapse SN-neutrino burst, however, detectors with huge statistics and remarkably greater flavor 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].

5. 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 Tsakstara and Kosmas [13] and Tsakstara [32].

The convolution (folding) is carried out with (i) the original cross sections obtained in section 3, and (ii) the low and intermediate energy neutrino spectra of section 4 in order to compute, first, flux averaged total cross sections, 〈σtot〉 and, then, corresponding supernova neutrino event rates and fluxes.

5.1. 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]

σcoh=0σcoh(εν)η(εν)dεν.    (17)

For a CdTe or CdZnTe detector material, the flux averaged cross sections, computed by inserting in Equation (17) the σcohν) from Figure 2B and the η(εν) from Figures 3, 4, for the isotopes 112, 114, 116Cd, are listed in Tables 3, 4 as described below.


Table 3. Calculated values for the flux averaged coherent cross sections 〈σcoh〉 (in units 10−39 cm2) for 112Cd, 114Cd, and 116Cd isotopes.


Table 4. Flux averaged coherent cross sections 〈σcoh〉, as in Table 3 but now referred to: (i) three different Supernova neutrino spectra determined from the parameters of: (a) Fermi Dirac parametrizations and (b) Power-Law parametrizations, and (ii) the energy spectra of Pion/muon decay at rest (DAR) neutrinos.

In Table 3 we list the flux averaged cross sections evaluated by adopting the neutrino distributions of the geo-neutrinos (see Figure 3A), the reactor neutrinos (see Figure 3B) and the solar neutrinos (see Figure 3C for the 8B and the hep solar neutrinos).

In the last two columns of Table 4 we tabulate the 〈σcoh〉 calculated for the distributions of Equations (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 Figure 4. The corresponding parameters are listed in Table 2.

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 = μ, τ [8789, 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 σtotincoh [11, 13, 69].

5.2. 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 4, for instance those referred to the SN neutrinos of the 112, 114, 116Cd isotopes, we estimate (potentially detectable) neutrino fluxes Φν that should arrive at each detector to create some typical scattering event rates Nev in the COBRA detector.

In general, the event rate Nev is related to the flux Φν reaching the nuclear detector with the expression [11, 32, 70]

dNνdtNev=NCdσtot(εν)Φν(εν).    (18)

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 mdet=100 kg and two cases of detector materials, i.e., the semiconductors (a) CdZnTe and (b) CdTe [3, 4, 26].

In the first step, we choose three SN neutrino scenarios in which the mean energies are: (i) 〈εν〉 = 12 MeV (corresponding to SN electron neutrinos νe), (ii) 〈εν〉 = 16 MeV (corresponding to SN electron anti-neutrinos ν~e), and (iii) 〈εν〉 = 24 MeV (corresponding to SN νx, ν~x, with x = μ, τ (anti)neutrinos of heavy leptons).

Then, based on Equation (18), we perform calculations assuming a total mass 100 kg of CdZnTe as COBRA detector which translates, for example, to approximately mCd = 10.6 kg mass of 114Cd isotope or equivalently a number of 114Cd atoms (nuclei) equal to NCdN114Cd=94.17NAvog.

In Equation (18), as total neutrino scattering cross sections, σtotν) we employ the values of flux averaged cross sections 〈σcoh〉 of Table 4 obtained through PL distribution for SN neutrino spectra (they refer to the three mean energies chosen above).

Finally, we choose four typical detection rates Nev as: (a) Nev=1 event s−1=3.15 × 107 events y−1, (b) Nev=1 event hr−1=8.76 × 103 events y−1, (c) Nev=1 event d−1=3.65 × 102 events y−1, and (d) Nev= 12 events y−1 and from Equation (18) we compute the corresponding SN ν fluxes Φν.

In a similar way, assuming that the COBRA detector contains 100 kg of the material CdTe, we find 13.5 kg 114Cd or about NCdN114Cd=120.11NAvog 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 Nev, we find the corresponding fluxes Φν reaching the COBRA CdTe detector.

By performing the steps we followed for 114Cd, for the other two Cd-isotopes, 112Cd and 116Cd, the resulting neutrino fluxes, for the chosen SN neutrino scenarios are listed in Table 5 (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.).


Table 5. Neutrino fluxes Φνν) and corresponding event rates Nev estimated to be recorded on 112, 114, 116Cd isotopes of two detector materials (CdTe and CdZnTe) of the COBRA experiment [3, 4, 26].

The above neutrino fluxes are of the same order with those of the Spallation Neutron Source (SNS) at ORLaND, Oak Ridge [4446]. We mention that the COHERENT 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×1011νμ/cm2 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 NevCOH = 168/y νμ neutrinos.

From the results of Table 5, we may define the ratio Nevν for the COHERENT experiment (RCOH) and for a special νμ neutrino case of the COBRA experiment (RCOB). For a comparison of these two experiments, we choose, for example, the results referred to the 112Cd isotope of CdTe material of the COBRA detector (sixth line from the beginning of Table 5 refers to νμ neutrinos). From these two ratios we find that R = RCOH/RCOB = 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 mdet = 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 5). 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.

6. 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, 116Cd 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, supernova 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 flux-averaged total coherent cross sections, 〈σcoh〉, reflect the mean neutrino signals generated in several terrestrial detectors (112, 114, 116Cd) 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.).

Author Contributions

JS and VT: nuclear structure calculations; OK: computation of folding (convolution) procedure.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


The present research was financially supported VT and JS by the Department of Informatics Engineering of the Technological Institute of Western Macedonia. Also, OK wishes to acknowledge the support of EPSRC via grand EP/N026136/1 Geometric Mechanics of Solids.


1. Ejiri H. Nuclear spin isospin responces for low-energy neutrinos. Phys Rep. (2000) 338:265–351. doi: 10.1016/S0370-1573(00)00044-2

CrossRef Full Text | Google Scholar

2. Ejiri H, Engel J, Kudomi N. Supernova-neutrino studies with 100Mo. Phys Lett B. (2002) 530:27–32. doi: 10.1016/S0370-2693(02)01349-7

CrossRef Full Text | Google Scholar

3. Zuber K. COBRA-double beta decay searches using CdTe detectors. Phys Lett B. (2001) 519:1–7. doi: 10.1016/S0370-2693(01)01056-5

CrossRef Full Text | Google Scholar

4. Zuber K. Spectroscopy of low energy solar neutrinos using CdTe detectors. Phys Lett B. (2003) 571:148–54. doi: 10.1016/j.physletb.2003.07.070

CrossRef Full Text | Google Scholar

5. Donnelly TW, Peccei RD. Neutral current effects in nuclei. Phys Rep. (1979) 50:185. doi: 10.1016/0370-1573(79)90010-3

CrossRef Full Text | Google Scholar

6. Kosmas TS, Oset E. Charged current neutrino-nucleus reaction cross sections at intermediate energies. Phys Rev C. (1996) 53:1409–15. doi: 10.1103/PhysRevC.53.1409

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Kolbe E, Kosmas TS. Recent highlights on neutrino-nucleus interactions. Springer Trac Mod. Phys. (2000) 163199–225. doi: 10.1007/BFb0109652

CrossRef Full Text | Google Scholar

8. Akimov D, Albert JB, An P, Awe C, Barbeau PS, Becker B, et al. Observation of coherent elastic neutrino-nucleus scattering. Science. (2017) 357:1123–6. doi: 10.1126/science.aao0990

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Scholberg K. Prospects for measuring coherent neutrino-nucleus elastic scattering at a stopped-pion neutrino source. Phys Rev D. (2006) 73:033005. doi: 10.1103/PhysRevD.73.033005

CrossRef Full Text | Google Scholar

10. Scholberg K. Coherent elastic neutrino-nucleus scattering. J Phys Conf Ser. (2015) 606:012010. doi: 10.1088/1742-6596/606/1/012010

CrossRef Full Text | Google Scholar

11. Tsakstara V, Kosmas TS. Low-energy neutral-current neutrino scattering on 128, 130Te isotopes. Phys Rev C. (2011) 83:054612. doi: 10.1103/PhysRevC.83.054612

CrossRef Full Text | Google Scholar

12. Tsakstara V, Kosmas TS, Wambach J. Studying low-energy astrophysical neutrinos with neutrino nucleus cross-section calculations and beta beam neutrino spectra. Prog Part Nucl Phys. (2011) 66:424–9. doi: 10.1016/j.ppnp.2011.01.045

CrossRef Full Text | Google Scholar

13. Tsakstara V, Kosmas TS. Analyzing astrophysical neutrinos through realistic nuclear structure calculations and the convolution procedure. Phys Rev C. (2011) 84:064620. doi: 10.1103/PhysRevC.84.064620

CrossRef Full Text | Google Scholar

14. Abe S, Ebihara T, Enomoto S, Furuno K, Gando Y, Ichimura K, et al. Precision Measurement of Neutrino Oscillation Parameters with KamLAND. Phys Rev Lett. (2008) 100:221803. doi: 10.1103/PhysRevLett.100.221803

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Gando A, Gando Y, Ichimura K, Ikeda H, Inoue K, Kibe Y, et al. Partial radiogenic heat model for Earth revealed by geoneutrino measurements. Nat Geosci. (2011) 4:647–51. doi: 10.1038/ngeo1205

CrossRef Full Text | Google Scholar

16. Bellini G, Benziger J, Bonetti S, Buizza Avanzini M, Caccianiga B, Cadonati L, et al. Observation of geo-neutrinos. Phys Lett B. (2010) 687:299–304. doi: 10.1016/j.physletb.2010.03.051

CrossRef Full Text | Google Scholar

17. Bellini G, Benziger J, Bick D, Bonetti S, Bonfini G, Buizza Avanzini M, et al. Precision measurement of the 7Be solar neutrino interaction rate in Borexino. Phys Rev Lett. (2011) 107:141302. doi: 10.1103/PhysRevLett.107.141302

CrossRef Full Text | Google Scholar

18. Bellini G, Benziger J, Bick D, Bonfini G, Bravo D, Buizza Avanzini M, et al. Measurement of geo-neutrinos from 1353 days of Borexino. Phys Lett B. (2013) 722:295–300. doi: 10.1016/j.physletb.2013.04.030

CrossRef Full Text | Google Scholar

19. Zuber K. Status of the double beta experiment COBRA. Prog Part Nucl Phys. (2006) 57:235–40. doi: 10.1016/j.ppnp.2005.12.010

CrossRef Full Text | Google Scholar

20. Wurm M, Beacom JF, Bezrukov LB, Bick D, Blümer J, Choubey S, et al. The next-generation liquid-scintillator neutrino observatory LENA. Astropart Phys. (2012) 35:685–732. doi: 10.1016/j.astropartphys.2012.02.011

CrossRef Full Text | Google Scholar

21. Kolbe E, Langanke K, Martinez-Pinedo G, Vogel P. Neutrino-nucleus reactions and nuclear structure. J Phys G. (2003) 29:2569–96. doi: 10.1088/0954-3899/29/11/010

CrossRef Full Text | Google Scholar

22. Sajjad Athar M, Singh SK. νe(ν̄e)-40Ar absorption cross sections for supernova neutrinos. Phys Lett B. (2004) 591:69–75. doi: 10.1016/j.physletb.2004.04.025

CrossRef Full Text

23. Zucchelli P. A novel concept for a ν̄e/νe neutrino factory: the beta-beam. Phys Lett B. (2002) 532:166–72. doi: 10.1016/S0370-2693(02)01576-9

CrossRef Full Text

24. Volpe C. What about a beta beam facility for low-energy neutrinos?. J Phys G. (2004) 30:L1–6. doi: 10.1088/0954-3899/30/7/L01

CrossRef Full Text | Google Scholar

25. Volpe C. Topical review on beta-beams. J Phys G. (2007) 34:R1–44. doi: 10.1088/0954-3899/34/1/R01

CrossRef Full Text

26. Zuber K. The status of the COBRA double-beta-decay experiment. Prog Part Nucl Phys. (2010) 64:267–9. doi: 10.1088/1742-6596/203/1/012070

CrossRef Full Text | Google Scholar

27. Smponias T, Kosmas OT. High energy neutrino emission from astrophysical jets in the galaxy. Adv High Energy Phys. (2015) 2015:921757. doi: 10.1155/2015/921757

CrossRef Full Text | Google Scholar

28. Smponias T, Kosmas OT. Neutrino emission from magnetized microquasar jets. Adv High Energy Phys. (2017) 2017:4962741. doi: 10.1155/2017/4962741

CrossRef Full Text | Google Scholar

29. Kosmas OT. Smponias T. Simulations of gamma-ray emission from magnetized microquasar jets. (2018) Adv High Energy Phys. 2018:9602960. doi: 10.1155/2018/9602960

CrossRef Full Text | Google Scholar

30. Aartsen MG, Ackermann M, Adams J, Aguilar JA, Ahlers M, Ahrens M, et al. Search for prompt neutrino emission from gamma-ray bursts with IceCube. Astrophys J Lett. (2015) 805:L5. doi: 10.1088/2041-8205/805/1/L5

CrossRef Full Text | Google Scholar

31. Adrian-Martinez S, Ageron M, Aharonian F, Aiello S, Albert A, Ameli F, et al. Letter of intent for KM3NeT 2.0. J Phys G Nucl Part Phys. (2016) 43:084001. doi: 10.1088/0954-3899/43/8/084001

CrossRef Full Text | Google Scholar

32. Tsakstara V. Convoluted-signals on 114Cd isotope from astrophysical and laboratory neutrino sources. Adv High Energy Phys. (2015) 2015:632131. doi: 10.1155/2015/632131

CrossRef Full Text | Google Scholar

33. Kosmas TS, Vergados JD, Civitarese O, Faessler A. Study of the muon number violating (μ, e conversion in a nucleus by using quasi-particle RPA. Nucl Phys A. (1994) 570:637–56. doi: 10.1016/0375-9474(94)90077-9

CrossRef Full Text | Google Scholar

34. Kosmas TS. Exotic μe conversion in nuclei: energy moments of the transition strength and average energy of the outgoing e. Nucl Phys A. (2001) 683:443–62. doi: 10.1016/S0375-9474(00)00471-1

CrossRef Full Text | Google Scholar

35. Tsoulos IG, Kosmas OT, Stavrou VN. DiracSolver: a tool for solving the Dirac equation. Comput Phys Commun. (2019) 236:237–43. doi: 10.1016/j.cpc.2018.10.010

CrossRef Full Text | Google Scholar

36. Kosmas O, Leyendecker S. Analysis of higher order phase fitted variational integrators. Adv Comput Math. (2016) 42:605–619. doi: 10.1007/s10444-015-9436-x

CrossRef Full Text | Google Scholar

37. Kosmas O, Leyendecker S. Variational integrators for orbital problems using frequency estimation. Adv Comput Math. (2019) 45:1–219. doi: 10.1007/s10444-018-9603-y

CrossRef Full Text | Google Scholar

38. Kosmas O, Papadopoulos D. Multisymplectic structure of numerical methods derived using non standard finite difference schemes. J Phys Confer Ser. 490:012205. doi: 10.1088/1742-6596/490/1/012205

CrossRef Full Text

39. Kosmas OT. Charged particle in an electromagnetic field using variational integrators. AIP Confer Proc. (2011) 1389:1927. doi: 10.1063/1.3636989

CrossRef Full Text | Google Scholar

40. Vogel P, Beacom JF. Angular distribution of neutron inverse beta decay ν̄e+pe++n. Phys Rev D. (1999) 60:053003. doi: 10.1103/PhysRevD.60.053003

CrossRef Full Text

41. Dye S. Geoneutrinos and the radioactive power of the Earth. Rev Geophys. (2012) 50:3007. doi: 10.1029/2012RG000400

CrossRef Full Text | Google Scholar

42. Fiorentini G, Mantovani E, Ricci B. Neutrinos and energetics of the Earth. Phys Lett B. (2003) 557:139–46. doi: 10.1016/S0370-2693(03)00193-X

CrossRef Full Text | Google Scholar

43. Fiorentini G, Ianni A, Korga G, Lissia M, Mantovani F, Miramonti L, et al. Nuclear physics for geo-neutrino studies. Phys Rev C. (2010) 81:034602. doi: 10.1103/PhysRevC.81.034602

CrossRef Full Text | Google Scholar

44. Avignone FT, Efremenko YV. ORLaND - a neutrino facility at the Spallation Neutron Source. Nucl Phys B. (2000) 87:304–8. doi: 10.1016/S0920-5632(00)00687-3

CrossRef Full Text | Google Scholar

45. Avignone FT, Efremenko YV. Neutrino-nucleus cross-section measurements at intense. pulsed spallation sources. J Phys G. (2003) 29:2615–28. doi: 10.1088/0954-3899/29/11/012

CrossRef Full Text | Google Scholar

46. Burman RL, Louis WC. Neutrino physics at meson factories and spallation neutron sources. J Phys G. (2003) 29:2499–512. doi: 10.1088/0954-3899/29/11/006

CrossRef Full Text | Google Scholar

47. Papoulias DK, Kosmas TS. Nuclear aspects of neutral current non-standard ν-nucleus reactions and the role of the exotic μe transitions experimental limits. Phys Lett B. (2014) 728:482–8. doi: 10.1016/j.physletb.2013.12.028

CrossRef Full Text | Google Scholar

48. Papoulias DK, Kosmas TS. Standard and nonstandard neutrino-nucleus reactions cross sections and event rates to neutrino detection experiments. Adv High Energy Phys. (2015) 215:763648. doi: 10.1155/2015/763648

CrossRef Full Text | Google Scholar

49. Engel J. Nuclear form-factors for the scattering of weakly interacting massive particles. Phys Lett B. (1991) 264:114–9. doi: 10.1016/0370-2693(91)90712-Y

CrossRef Full Text | Google Scholar

50. De Vries H, De Jager CW, De Vries C. Nuclear charge-density-distribution parameters from elastic electron scattering. Atomic Data Nucl Data Tables. (1987) 36:495–536. doi: 10.1016/0092-640X(87)90013-1

CrossRef Full Text | Google Scholar

51. Drukier A, Stodolsky L. Principles and applications of a neutral-current detector for neutrino physics and astronomy. Phys Rev D. (1984) 30:2295–309. doi: 10.1103/PhysRevD.30.2295

CrossRef Full Text | Google Scholar

52. Horowitz CJ, Ahmed Z, Jen C-M, Rakhman A, Souder PA, Dalton MM, et al. Weak charge form factor and radius of 208Pb through parity violation in electron scattering. Phys Rev C. (2012) 85:032501. doi: 10.1103/PhysRevC.85.032501

CrossRef Full Text | Google Scholar

53. Kosmas TS Vergados JD. Nuclear densities with fractional occupation probabilities of the states. Nucl Phys A. (1992) 536:72–86. doi: 10.1016/0375-9474(92)90246-G

CrossRef Full Text | Google Scholar

54. Chiang HC, Oset E, Kosmas TS, Faessler A, Vergados JD. Coherent and incoherent (mu−, e−) conversion in nuclei. Nucl Phys A. (1993) 559:526–42. doi: 10.1016/0375-9474(93)90259-Z

CrossRef Full Text | Google Scholar

55. Bahcall JN. Neutrino Astrophysics. Melbourne, VIC: Cambridge University Press (1989).

Google Scholar

56. Barranco J, Miranda OG, Rashba TI. Probing new physics with coherent neutrino scattering off nuclei. J High Energy Phys. (2005) 05:021. doi: 10.1088/1126-6708/2005/12/021

CrossRef Full Text | Google Scholar

57. Barranco J, Bolanos A, Miranda OG, Rashba TI. Tensorial NSI and Unparticle physics in neutrino scattering. Int J Mod Phys A. (2012) 27:1250147. doi: 10.1142/S0217751X12501473

CrossRef Full Text | Google Scholar

58. Vergados JD, Avignone FT, Giomataris I. Coherent neutral current neutrino-nucleus scattering at a spallation source: a valuable experimental probe. Phys Rev D. (2009) 79:113001. doi: 10.1103/PhysRevD.79.113001

CrossRef Full Text | Google Scholar

59. Vergados JD, Giomataris Y. Dedicated supernova detection by a network of neutral current spherical TPC detectors. Phys Atom Nucl. (2007) 70:140–9. doi: 10.1134/S1063778807010164

CrossRef Full Text | Google Scholar

60. Giomataris Y, Vergados JD. A network of neutral current spherical TPCs for dedicated supernova detection. Phys Lett B. (2006) 634:23–9. doi: 10.1016/j.physletb.2006.01.040

CrossRef Full Text | Google Scholar

61. Vogel P, Engel J. Neutrino electromagnetic form factors. Phys Rev D. (1989) 39:3378–83. doi: 10.1103/PhysRevD.39.3378

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Kolbe E. Differential cross sections for neutrino scattering on 12C. Phys Rev C. (1996) 54:1741–8. doi: 10.1103/PhysRevC.54.1741

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Langanke K, Martinez-Pinedo G. Nuclear weak-interaction processes in stars. Rev Mod Phys. (2003) 75:819–62. doi: 10.1103/RevModPhys.75.819

CrossRef Full Text | Google Scholar

64. Juodagalvis A, Langanke K, Martinez-Pinedo G, Hix WR, Dean DJ, Sampaio JM., Neutral-current neutrino-nucleus cross sections for A~50−65 nuclei. Nucl Phys A. (2005) 747:87–108. doi: 10.1016/j.nuclphysa.2004.09.005

CrossRef Full Text | Google Scholar

65. Langanke K. Weak interaction, nuclear physics and supernovae. Acta Phys Polon B. (2008) 39:265–82.

Google Scholar

66. Haxton WC. Nuclear response of water Cherenkov detectors to supernova and solar neutrinos. Phys Rev D. (1987) 36:2283–92. doi: 10.1103/PhysRevD.36.2283

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Donnelly TW, Walecka JD. Semi-leptonic weak and electromagnetic interactions in nuclei with application to 16O. Phys Lett B. (1972) 41:275–80. doi: 10.1016/0370-2693(72)90577-1

CrossRef Full Text | Google Scholar

68. Donnelly TW, Walecka JD. Elastic magnetic electron scattering and nuclear moments. Nucl Phys A. (1973) 201:81–106. doi: 10.1016/0375-9474(73)90689-1

CrossRef Full Text | Google Scholar

69. Tsakstara V, Kosmas TS. Nuclear responses of 64, 66Zn isotopes to supernova neutrinos. Phys Rev C. (2012) 86:044618. doi: 10.1103/PhysRevC.86.044618

CrossRef Full Text | Google Scholar

70. Balasi KG, Ydrefors E, Kosmas TS. Theoretical study of neutrino scattering off the stable even Mo isotopes at low and intermediate energies. Nucl Phys A. (2011) 868–9:82–98. doi: 10.1016/j.nuclphysa.2011.08.003

CrossRef Full Text | Google Scholar

71. Anderson AJ, Conrad JM, Figueroa-Feliciano E, Ignarra C, Karagiorgi G, Scholberg K, et al. Measuring active-to-sterile neutrino oscillations with neutral current coherent neutrino-nucleus scattering. Phys Rev D. (2012) 86:013004. doi: 10.1103/PhysRevD.86.013004

CrossRef Full Text | Google Scholar

72. Louis WC. Searches for muon-to-electron (anti) neutrino flavor change. Prog Part Nucl Phys. (2009) 63:51–73. doi: 10.1016/j.ppnp.2009.01.002

CrossRef Full Text | Google Scholar

73. de Gouvea A, et al. Working Group Report: Neutrinos (Intensity Frontier Neutrino Working Group). arXiv:1310.4340 [hep-ex] (2013).

74. Bahcall JN, Holstein BR. Solar neutrinos from the decay of 8B. Phys Rev C. (1986) 33:2121–7. doi: 10.1103/PhysRevC.33.2121

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Bahcall JN, Urlich RK. Solar models, neutrino experiments, and helioseismology. Rev Mod Phys. (1988) 60:297–372. doi: 10.1103/RevModPhys.60.297

CrossRef Full Text | Google Scholar

76. Bahcall JN, Basu S, Pinsonneault M, Serenelli AM. Helioseismological implications of recent solar abundance determinations. Astrophys J. (2005) 618:1049–56. doi: 10.1086/426070

CrossRef Full Text | Google Scholar

77. Davis BR, Vogel P, Mann FM, Schenter RE. Reactor antineutrino spectra and their application to antineutrino-induced reactions. Phys Rev C. (1979) 19:2259–66. doi: 10.1103/PhysRevC.19.2259

CrossRef Full Text | Google Scholar

78. Declais Y, de Kerret H, Lefiévre B, Obolensky M, Etenko A, Kozlov Yu, et al. Study of reactor antineutrino interaction with proton at Bugey nuclear power plant. Phys Lett B. (1994) 338:383–9. doi: 10.1016/0370-2693(94)91394-3

CrossRef Full Text | Google Scholar

79. Tengblad O, Aleklett K, Von Dincklage R, Lund E, Nyman G, Rudstam G. Integral ν̄-spectra derived from experimental β-spectra of individual fission products. Nucl Phys A. (1989) 503:136–60. doi: 10.1016/0375-9474(89)90258-3

CrossRef Full Text

80. Kerman S, Sharma V, Deniz M, Wong HT, Chen J-W, Li HB, et al. Coherency in neutrino-nucleus elastic scattering. Phys Rev D. (2016) 93:113006. doi: 10.1103/PhysRevD.93.113006

CrossRef Full Text | Google Scholar

81. Deniz M, Sevda B, Kerman S, Ajjaq A, Singh L, Wong HT, et al. Constraints on scalar-pseudoscalar and tensorial non-standard interaction and tensorial unparticle couplings from neutrino-electron scattering. Phys Rev D. (2017) 95:033008. doi: 10.1103/PhysRevD.95.033008

CrossRef Full Text | Google Scholar

82. Agnolet G, Baker G, Barker D, Beck R, Carroll TJ, Cesar J, et al. Background Studies for the MINER coherent neutrino scattering reactor experiment . Nucl Instrum Methods Phys Res A. (2017) 853:53–60. doi: 10.1016/j.nima.2017.02.024

CrossRef Full Text | Google Scholar

83. Bethe HA. Supernova mechanisms. Rev Mod Phys. (1990) 62:801–66. doi: 10.1103/RevModPhys.62.801

CrossRef Full Text | Google Scholar

84. Fröhlich C, Martínez-Pinedo G, Liebendörfer M, Thielemann FK, Bravo E, Hix WR, et al. Neutrino-induced nucleosynthesis of A>64 nuclei: the vp process. Phys Rev Lett. (2006) 96:142502. doi: 10.1103/PhysRevLett.96.142502

CrossRef Full Text | Google Scholar

85. Janka H-Th, Mueller B. Neutrino-driven type-II supernova explosions and the role of convection. Phys Rep. (1995) 256:135–56.

Google Scholar

86. Janka H-Th, Langanke K, Mareka A, Martinez-Pinedo G, Mullera B. Theory of core-collapse supernovae. Phys Rep. (2007) 442:38–74. doi: 10.1016/j.physrep.2007.02.002

CrossRef Full Text | Google Scholar

87. Janka HT, Hillebrand W. Neutrino emission from type II supernovae – an analysis of the spectra. Astron Astrophys. (1989) 224:49–56.

Google Scholar

88. Keil MT, Raffelt GG, Janka H-T. Monte Carlo study of supernova neutrino spectra formation. Astrophys J. (2003) 590:971–91. doi: 10.1086/375130

CrossRef Full Text | Google Scholar

89. Raffelt GG, Keil MT, Buras R, Janka H-T, Rampp M. Supernova neutrinos: flavor-dependent fluxes and spectra. Proc Neutrino Oscill Origin. (2004) 03:380–7. doi: 10.1142/9789812703101_0044

CrossRef Full Text | Google Scholar

90. Duan H, Fuller GM, Qian Y-Z. Stepwise spectral swapping with three neutrino flavors. Phys Rev D. (2008) 77:085016. doi: 10.1103/PhysRevD.77.085016

CrossRef Full Text | Google Scholar

91. Jachowicz N, McLaughlin GC. Reconstructing supernova-neutrino spectra using low-energy beta beams. Phys Rev Lett. (2006) 96:172301. doi: 10.1103/PhysRevLett.96.172301

PubMed Abstract | CrossRef Full Text | Google Scholar

92. Tsakstara V, Kosmas TS, Divari PC, Sinatkas J. The interpretation of SN?v signals in terrestrial experiments through the folding procedure. AIP Conf Proc. (2009) 1180:140–4. doi: 10.1063/1.3266092

CrossRef Full Text | Google Scholar

93. Tsakstara V, Kosmas TS. Neutrino-nucleus reactions in terrestrial experiments and astrophysics. Prog Part Nucl Phys. (2010) 64:407–10. doi: 10.1016/j.ppnp.2009.12.060

CrossRef Full Text | Google Scholar

94. Scholberg K. Supernova signatures of neutrino mass ordering. J Phys G.(2018) 45:014002. doi: 10.1088/1361-6471/aa97be

CrossRef Full Text | Google Scholar

95. Vissani F, Pagliaroli G. The diffuse supernova neutrino background: expectations and uncertainties derived from SN1987A. Astron Astrophys. (2011) 528:L1. doi: 10.1051/0004-6361/201016109

CrossRef Full Text | Google Scholar

96. Moeller K, Suliga AM, Tamborra I, Denton PB. Measuring the supernova unknowns at the next-generation neutrino telescopes through the diffuse neutrino background. J Cosmol Astropart Phys. (2018) 1805:066. doi: 10.1088/1475-7516/2018/05/066

CrossRef Full Text

97. Bohr A, Mottelson BR. Nuclear Structure. New York, NY: Benjamin (1969).

A. Appendix

A.1. Three-Point Formulas for Empirical Energy Gaps Δn,pexp of Neutrons and Protons

The empirical energy gaps for neutrons, Δnexp, and protons, Δnexp, needed at the BCS level to construct the ground state wave function of the detector nucleus (A, Z), are computed through the respective separation energies for neutrons, Sn or protons, Sp of the isotope (A, Z) and also those of the neighboring nuclear isotopes with N ± 1 neutrons or Z ± 1 protons, respectively, by employing the expressions

Δnexp=-14[Sn(N-1,Z)-2Sn(N,Z)+Sn(N+1,Z)]    (A1)
Δpexp=-14[Sp(N,Z-1)-2Sp(N,Z)+Sp(N,Z+1)]    (A2)

The above equations are known as the three-point formulas (see, e.g., [69, 97]).

A.2. Normalization of the ν Energy Distributions η(εν) Adopted in This Work

The distributions η(εν) adopted in the present work (see section 4), are considered to be normalized in such a way that

0η(εν)dεν=1.    (A3)

For example, in the case of ηνeν) of Figure 3D, the normalization gives


where we have used ενmin=0 and ενmax=Mμ/2

A.3. Parametrization of Supernova Neutrino Energy Spectra

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


where w0=εν/3 is the width of the identical FD and PL distributions [13].

A.3.1. Fermi-Dirac (FD) Energy Distribution

By introducing the degeneracy parameter ndg (equal to the ratio of the chemical potential μ divided by the neutrino temperature T, i.e., ndg = μ/T), the Fermi-Dirac energy distribution reads

ηFD[x,T,ndg]=F(ndg)1Tx21+e(x-ndg),  x=ενT.    (A4)

In this case, the width of the spectrum is reduced compared to the corresponding thermal spectrum (pinching effect). The normalization constant F2(ndg) of this distribution depends on the degeneracy parameter ndg and is given by the relation

1F(ndg)0x2ex-ndg+1dx.    (A5)

Inserting Equation (A5) into Equation (A4), we take

ηFD[εν,T,ndg]=[0x2ex-ndg+1dx]-1(εν2/T3)1+e(εν/T-ndg).    (A6)

A.3.2. Power-Law Energy Distribution

The SN-neutrino energy spectra can be fitted by using a Power-Law energy distribution of the form [88]

ηPL[εν,α]=C(ενεν)αe-(α+1)(εν/εν),    (A7)

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

C0(ενεν)αe-(α+1)(εν/εν)dεν=1.    (A8)

From the later equation we find

C=(α+1)α+1Γ(α+1)εν,    (A9)

therefore, Equation (A7) becomes

ηPL[εν,α]=(α+1)α+1Γ(α+1)εναενα+1e-(α+1)(εν/εν).    (A10)

Keywords: nuclear detector responses, neutrino nucleus cross sections, supernova neutrino detection, neutral-current neutrino-nucleus processes, quasi-particle random phase approximation

PACS numbers: 26.50.+x, 25.30.Pt, 97.60.Bw, 25.30.-c, 23.40.Bw, 21.60.Jz

Citation: Sinatkas J, Tsakstara V and Kosmas O (2019) Simulated Neutrino Signals of Low and Intermediate Energy Neutrinos on Cd Detectors. Front. Phys. 7:42. doi: 10.3389/fphy.2019.00042

Received: 20 September 2018; Accepted: 06 March 2019;
Published: 12 April 2019.

Edited by:

Yoshitaka Kuno, Osaka University, Japan

Reviewed by:

Luca Fiorini, University of Valencia, Spain
Andreas Redelbach, Frankfurt Institute for Advanced Studies, Germany

Copyright © 2019 Sinatkas, Tsakstara and Kosmas. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: J. Sinatkas,