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

TL;DR

This paper evaluates neutrino-nucleus reaction cross sections for cadmium isotopes relevant to current neutrino detection experiments and simulates expected neutrino signals from various astrophysical and laboratory sources.

Contribution

It provides new cross section calculations for Cd isotopes and simulates neutrino signals for multiple neutrino sources using modern techniques.

Findings

Cross sections for neutrino interactions with Cd isotopes are calculated.

Simulated neutrino signals for astrophysical and laboratory sources are presented.

Coherent neutrino-nucleus reactions dominate the neutral current scattering in Cd detectors.

Abstract

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).

Tables

Table 1. Table 1: Pairing parameters g p ​ a ​ i ​ r p superscript subscript 𝑔 𝑝 𝑎 𝑖 𝑟 𝑝 g_{pair}^{p} (for protons), and g p ​ a ​ i ​ r n superscript subscript 𝑔 𝑝 𝑎 𝑖 𝑟 𝑛 g_{pair}^{n} (for neutrons) determining the monopole pairing interactions for each of the studied isotopes. The obtained theoretical values of the energy gaps (in units of MeV), Δ p t ​ h subscript superscript Δ 𝑡 ℎ 𝑝 \Delta^{th}_{p} (for protons) and Δ n t ​ h subscript superscript Δ 𝑡 ℎ 𝑛 \Delta^{th}_{n} (for neutrons), are also shown for comparison with the empirical ones. As can be seen, the corresponding empirical energy gaps, Δ p , n e ​ x ​ p subscript superscript Δ 𝑒 𝑥 𝑝 𝑝 𝑛 \Delta^{exp}_{p,n} are well reproduced. Values of the harmonic oscillator size parameter, b 𝑏 b , for each of the isotopes 112 , 114 , 116 112 114 116 {}^{112,114,116}\! Cd are also given in this Table.

Isotope	Z, N	Abundance ($\%$)	b (fm)	$g_{pair}^n$	$g_{pair}^p$	${\Delta }_p^{exp}$	${\Delta }_p^{th}$	${\Delta }_n^{exp}$	${\Delta }_n^{th}$
${}_{}{}^{112}Cd$	48, 64	24.13	2.208	1.001	1.064	1.516	1.512	1.320	1.322
${}_{}{}^{114}Cd$	48, 66	28.73	2.214	0.956	0.975	1.441	1.441	1.351	1.351
${}_{}{}^{116}Cd$	48, 68	7.50	2.219	1.069	1.043	1.432	1.432	1.371	1.372

Table 2. Table 2: Corresponding values of parameters for equivalent Fermi-Dirac (FD) and Power-Law (PL) distributions (SN neutrino energy spectra) of Fig. 4 . The selected flavour dependent mean ν 𝜈 \nu -energy values (describing the PL distribution), ⟨ ε ν ⟩ delimited-⟨⟩ subscript 𝜀 𝜈 \langle\varepsilon_{\nu}\rangle in MeV, have been chosen as model values for ν e subscript 𝜈 𝑒 \nu_{e} (10-12) MeV, ν ~ e subscript ~ 𝜈 𝑒 \widetilde{\nu}_{e} (15-18) MeV and ν x subscript 𝜈 𝑥 \nu_{x} , where x = ν μ subscript 𝜈 𝜇 \nu_{\mu} , ν τ subscript 𝜈 𝜏 \nu_{\tau} , ν ~ μ subscript ~ 𝜈 𝜇 \widetilde{\nu}_{\mu} , ν ~ τ subscript ~ 𝜈 𝜏 \widetilde{\nu}_{\tau} (22-26) MeV.

Equivalent Fermi-Dirac and Power-Law Supernova Neutrino Spectra
Parameter			Temperature (in MeV)
Width ($w$)	Pinching ($\alpha$)	Degeneracy ($n_{dg}$)	$⟨{\epsilon }_{\nu }⟩=$10	$⟨{\epsilon }_{\nu }⟩=$12	$⟨{\epsilon }_{\nu }⟩=$16	$⟨{\epsilon }_{\nu }⟩=$20	$⟨{\epsilon }_{\nu }⟩=$24 (MeV)
0.7	5.1	4.4	2.14	2.57	3.42	4.28	5.13
0.8	3.7	2.7	2.58	3.10	4.14	5.17	6.20
0.9	2.7	1.1	2.98	3.57	4.77	5.96	7.15

Table 3. Table 3: Calculated values for the flux averaged coherent cross sections ⟨ σ c ​ o ​ h ⟩ delimited-⟨⟩ subscript 𝜎 𝑐 𝑜 ℎ \langle\sigma_{coh}\rangle (in units 10 − 39 superscript 10 39 10^{-39} cm ) 2 {}^{2}) for 112 Cd, 114 Cd and 116 Cd isotopes. The neutrino sources distributions of neutrino beams coming from: (i) Geo-neutrinos, (ii) Reactor neutrinos and (iii) Solar neutrinos have been used in the folding procedure.

Flux Averaged Cross Sections $⟨{\sigma }_{coh}⟩$ $({10}^{-39}c{m}^2$)
	Geo-Neutrinos			Reactor Neutrinos			Solar Neutrinos
Isotope	40K	238U	232Th	235U	238U	239Pu	8B	hep
112Cd	0.142	1.410	0.911	0.180	0.477	9.001	7.970	9.333
114Cd	0.151	1.504	0.973	0.192	0.509	9.604	8.504	9.957
116Cd	0.161	1.602	1.036	0.205	0.542	1.023	9.055	10.598

Table 4. Table 4: Flux averaged coherent cross sections ⟨ σ c ​ o ​ h ⟩ delimited-⟨⟩ subscript 𝜎 𝑐 𝑜 ℎ \langle\sigma_{coh}\rangle , 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.

Flux Averaged Cross Sections $⟨{\sigma }_{coh}⟩$ $({10}^{-39}c{m}^2$)
Isotope	Supernova Neutrinos								Pion-muon DAR Neutrinos
	Fermi-Dirac (FD)				Power Law (PL)						${\nu }_e$    Spect.	${\tilde{\nu }}_{\mu }$ Spect.
	$T=$	3.10	4.14	6.20	$⟨{\epsilon }_{\nu }⟩=$	12	16	24
112Cd		2.484	4.184	8.132		2.489	4.180	8.142			12.338	14.960
114Cd		2.648	4.458	8.648		2.653	4.453	8.658			13.110	15.881
116Cd		2.817	4.739	9.178		2.823	4.734	9.189			13.801	16.824

Table 5. Table 5: Neutrino fluxes Φ ν ​ ( ε ν ) subscript Φ 𝜈 subscript 𝜀 𝜈 \Phi_{\nu}(\varepsilon_{\nu}) and corresponding event rates N e ​ v subscript 𝑁 𝑒 𝑣 N_{ev} estimated to be recorded on 112,114,116 Cd isotopes of two detector materials (CdTe and CdZnTe) of the COBRA experiment Zuber-plb-1 ; Zuber-plb-2 ; Zuber-PPNP-10 . They refer to the case of supernova neutrinos with mean energies ⟨ ε ν ⟩ delimited-⟨⟩ subscript 𝜀 𝜈 \langle\varepsilon_{\nu}\rangle = 12, 16 and 24 MeV. N A ​ v ​ o ​ g subscript 𝑁 𝐴 𝑣 𝑜 𝑔 N_{Avog} is the Avogadro’s number. In the last column, N e ​ v C ​ O ​ H superscript subscript 𝑁 𝑒 𝑣 𝐶 𝑂 𝐻 N_{ev}^{COH} and Φ ν C ​ O ​ H superscript subscript Φ 𝜈 𝐶 𝑂 𝐻 \Phi_{\nu}^{COH} describe COHERENT experiment values (see the text).

Supernova neutrino coherent fluxes ${\Phi }_{\nu }$ (s${}_{}{}^{-1}\text{c}{m}^{-2})$ and event rates $N_{ev}$
					$N_{ev}$ = 1/s	$N_{ev}$ = 1/hr	$N_{ev}$ = 1/d	$N_{ev}$ = 12/y	$N_{ev}^{COH}$ = 168/y
Isotope	Detector	Atoms ($N_{Avog}$)	m (kg)	$⟨{\epsilon }_{\nu }⟩$ (MeV)	${\Phi }_{\nu }(\times {10}^5)$	${\Phi }_{\nu }(\times {10}^9)$	${\Phi }_{\nu }(\times {10}^{10})$	${\Phi }_{\nu }(\times {10}^{11})$	${\Phi }_{\nu }^{COH}(\times {10}^{11})$
				12	5.14	1.85	4.44	13.33
112Cd	CdTe	100.08	11.25	16	3.06	1.10	2.65	7.94
				24	1.57	0.57	1.36	4.08	1.70
				12	5.15	1.85	4.45	13.34
112Cd	CdZnTe	100.05	11.25	16	3.06	1.10	2.64	7.94
				24	1.57	0.57	1.36	4.08	1.70
				12	3.98	1.45	3.44	10.32
114Cd	CdTe	121.29	13.50	16	2.37	0.86	2.05	6.15
				24	1.22	0.45	1.05	3.16	1.70
				12	4.26	1.85	3.67	11.00
114Cd	CdZnTe	121.26	10.60	16	2.53	1.10	2.18	6.55
				24	12.98	0.57	1.12	3.36	1.70
				12	14.10	5.08	12.19	36.56
116Cd	CdTe	32.18	3.62	16	8.41	3.03	7.27	21.80
				24	4.33	1.56	3.74	11.23	1.70
				12	16.00	5.76	13.83	41.48
116Cd	CdZnTe	32.18	3.62	16	9.53	3.43	8.23	24.70
				24	4.89	1.76	4.23	12.66	1.70

Equations

$$\frac{d\sigma }{d\Omega }=\frac{{\mathrm{G}}_{\mathrm{F}}^2}{4{\pi }^2}{\epsilon }_{\nu }^2(1+\cos \vartheta )\frac{Q_w^2}{4}\mathcal{F}(q^2{)}^2$$

$${\cal F}(q^{2})=\frac{1}{Q_{w}}\Bigg{[}(1-4\sin^{2}\Theta_{w})ZF_{Z}(q^{2})-NF_{N}(q^{2})\Bigg{]}\,$$

$$Q_w=(1-4{\sin }^2{\Theta }_w)Z-N.$$

$$F_k(q^2)=\frac{k}{4\pi }\int j_0(qr){\rho }_{n,p}(r)r^{2}dr,k=N,Z$$

$$\frac{d\sigma(\varepsilon_{\nu},\vartheta)}{dcos\vartheta}=\frac{G^{2}_{F}}{2\pi}(1+\cos\vartheta)\varepsilon^{2}_{\nu}\Bigg{[}f_{V}^{p}Z+f_{V}^{n}N\Bigg{]}^{2}F_{Z}^{2}(q^{2})\,.$$

$$f_V^p=\frac{1}{2}-2{\sin }^2{\Theta }_W,f_V^n=-\frac{1}{2}.$$

$$q=2{\epsilon }_{\nu }\sin (\vartheta /2),$$

$$\sigma_{tot}(\varepsilon_{\nu})=\frac{G^{2}_{F}}{8\pi}\Bigg{[}(1-4\sin^{2}\Theta_{W})Z-N\Bigg{]}^{2}\varepsilon^{2}_{\nu}\,.$$

$$\frac{d\sigma ({\epsilon }_{\nu },T_A)}{d{T}_A}=\frac{G_F^2}{4\pi }{Q}_W^{2}M\left(1-\frac{M{T}_A}{2{\epsilon }_{\nu }^2}\right)F(2M{T}_A^2),$$

$$T_A^{min}=\frac{T_A}{2}+\sqrt{\frac{T_A}{2}(M_A+\frac{T_A}{2})},$$

$${\nu }_l({\tilde{\nu }}_l)+{}^{112,114,116}Cd\to {}^{112,114,116}C{d}^{\ast }+{\nu }_l^{\prime }({\tilde{\nu }}_l^{\prime }),$$

$$F_k(q^2)=\frac{1}{k}\sum _j\hat{j}⟨(n\ell )j|j_0(qr)|(n\ell )j⟩(V_j^k{)}^2,k=N,Z$$

$$\frac{d{N}_{\nu }({\epsilon }_{\nu })}{d{\epsilon }_{\nu }}\equiv \eta ({\epsilon }_{\nu })$$

$$(A,Z)\to (A,Z+1)+e^{-}+{\tilde{\nu }}_e+Q_{\beta }.$$

$${\mu }^{+}\to e^{+}+{\nu }_e+{\tilde{\nu }}_{\mu }.$$

$${\eta }_{{\nu }_e}({\epsilon }_{\nu })=96{\epsilon }_{\nu }^2{M}_{\mu }^{-4}\left(M_{\mu }-2{\epsilon }_{\nu }\right),$$

$${\eta }_{{\tilde{\nu }}_{\mu }}({\epsilon }_{\nu })=16{\epsilon }_{\nu }^2{M}_{\mu }^{-4}\left(3M_{\mu }-4{\epsilon }_{\nu }\right),$$

$$⟨{\sigma }_{coh}⟩={\int }_0^{\infty }{\sigma }_{coh}({\epsilon }_{\nu })\eta ({\epsilon }_{\nu })d{\epsilon }_{\nu }.$$

$$\frac{d{N}_{\nu }}{dt}\equiv N_{ev}=N_{Cd}{\sigma }_{tot}({\epsilon }_{\nu }){\Phi }_{\nu }({\epsilon }_{\nu }).$$

$$N_{ev}=ϵN_{Cd}{\sigma }_{tot}({\epsilon }_{\nu }){\Phi }_{\nu }({\epsilon }_{\nu }),$$

$${\Delta }_n^{exp}=-\frac{1}{4}\left[S_n(N-1,Z)-2S_n(N,Z)+S_n(N+1,Z)\right]$$

$${\Delta }_p^{exp}=-\frac{1}{4}\left[S_p(N,Z-1)-2S_p(N,Z)+S_p(N,Z+1)\right]$$

$${\int }_0^{\infty }\eta ({\epsilon }_{\nu })d{\epsilon }_{\nu }=1.$$

$$(96/M_{\mu }^4)\left[M_{\mu }{\int }_0^{M_{\mu }/2}{\epsilon }_{\nu }^{2}d{\epsilon }_{\nu }-2{\int }_0^{M_{\mu }/2}{\epsilon }_{\nu }^{3}d{\epsilon }_{\nu }\right]=1$$

$$w=\sqrt{⟨{\epsilon }_{\nu }^2⟩-⟨{\epsilon }_{\nu }{⟩}^2}/w_0$$

$${\eta }_{FD}[x,T,n_{dg}]=F(n_{dg})\frac{1}{T}\frac{x^2}{1+e^{(x-n_{dg})}},x=\frac{{\epsilon }_{\nu }}{T}.$$

$$\frac{1}{F(n_{dg})}\equiv {\int }_0^{\infty }\frac{x^2}{e^{x-n_{dg}}+1}dx.$$

$${\eta }_{FD}[{\epsilon }_{\nu },T,n_{dg}]={\left[{\int }_0^{\infty }\frac{x^2}{e^{x-n_{dg}}+1}dx\right]}^{-1}\frac{({\epsilon }_{\nu }^2/T^3)}{1+e^{({\epsilon }_{\nu }/T-n_{dg})}}.$$

$$\eta_{PL}[\langle\varepsilon_{\nu}\rangle,\alpha]=C\left(\frac{\varepsilon_{\nu}}{\langle\varepsilon_{\nu}\rangle}\right)^{\alpha}e^{-(\alpha+1)(\varepsilon_{\nu}/\langle\varepsilon_{\nu}\rangle)}\,,\$$

$$C\int_{0}^{\infty}\left(\frac{\varepsilon_{\nu}}{\langle\varepsilon_{\nu}\rangle}\right)^{\alpha}e^{-(\alpha+1)(\varepsilon_{\nu}/\langle\varepsilon_{\nu}\rangle)}d\varepsilon_{\nu}\,=1\,.\$$

$$C=\frac{(\alpha +1{)}^{\alpha +1}}{\Gamma (\alpha +1)⟨{\epsilon }_{\nu }⟩},$$

Full text

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

††preprint: APS/IOA-165

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 Ejiri-PR00 ; Ejiri-plb-02 ; Zuber-plb-1 ; Zuber-plb-2 . To enable further progress, relevant nuclear model calculations, across a wide energy range and in various nuclear isotopes, may provide significant results DonPe ; Kos-Ose ; kolbe-kos .

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 Ejiri-PR00 ; Ejiri-plb-02 ; Zuber-plb-1 ; Zuber-plb-2 and specifically experiments like the COHERENT where recently coherent neutrino nucleus scattering events have been measured for a first time CEnN_neut ; Scholberg ; Scholb_IoP-2010 . Furthermore, the presence of important nuclear effects impact the interaction cross sections as well as the final nuclear states reached through the scattering process TSAK-prc1 ; kos-eric-11 ; tsak_prc2 . The nuclear effects also affect the rebuilding of the incoming $\nu$-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 KamLAND ; KamLAND-11 ; Borexino-Col ; Borexino ; Borexino-13 ; Zub-SNO+ ; LENA .

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 KamLAND ; KamLAND-11 ; Borexino-Col ; Borexino ; Borexino-13 ; Zub-SNO+ , the neutrino energy range of 10 up to about 100-120 MeV (intermediate energy) covers a set of $\nu$-physics topics in the front of nuclear structure physics and astro-particle physics such as core-collapse supernovae dynamics and dark matter detection TSAK-prc1 ; kos-eric-11 ; tsak_prc2 ; Kolbe-Lag-Pin-03 ; Athar-Sing , 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 bet-beam-Zuc ; bet-beam-Vol ; Volpe07 ; Zuber-PPNP-10 .

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) Ejiri-PR00 ; Ejiri-plb-02 ; Kolbe-Lag-Pin-03 allowing us to investigate deep into the astrophysical objects Kolbe-Lag-Pin-03 ; Smpon-Ody-2015 ; Smpon_Ody-2017 ; Smpon-Ody-2018 . 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 Zub-SNO+ ; LENA ; Zuber-PPNP-10 (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) IceCube_2015 ; KM3Net_2016 . 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 $\nu$ 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 Zuber-plb-1 ; Zuber-plb-2 ; Zuber-PPNP-10 ; Tsak_AHEP . 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$(\nu,\nu^{\prime})^{112,114,116}$Cd∗, and anti-neutrinos, 112,114,116Cd$(\widetilde{\nu},\widetilde{\nu}^{\prime})^{112,114,116}$Cd∗, computed for the coherent channel by using a refinement of the quasi particle random phase approximation (QRPA) Tsak_AHEP ; TSAK-prc1 ; tsak_prc2 ; Kosmas-94 ; Kosmas01 . (ii) Reliable descriptions of the shapes of neutrino energy distributions coming out of numerical simulations of distributions in neutrino-energies $\varepsilon_{\nu}\leq 100-120$ MeV (for the above mentioned $\nu$-sources). (iii) Modern computational tools ody05 ; ody01 ; ody02 ; ody03 ; ody04 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) Zuber-plb-1 ; Zuber-plb-2 ; Zuber-PPNP-10 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+) Zub-SNO+ ; LENA , are expected to give useful answers to several questions of geological importance regarding the precise geo-$\nu$ fluxes and abundances of natural radioactive elements ($K$, $U$, $Th$) in the Earth’s interior Vogel-Beacom ; Dye-2011 ; Fiorentini-2003 ; Fiorentini-2010 .

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}\,{\textrm{c}m}^{2}$ in the neutrino-energy region ($\varepsilon_{\nu}\leq$50 MeV). This process has very recently been observed in the COHERENT experiment at a 6.7 $\sigma$ confidence level (CL), by using a low-background CsI[Na] scintillator CEnN_neut ; Scholberg ; Scholb_IoP-2010 . The detector was exposed to a $\nu_{\mu}$ neutrino beam coming from the Spallation Neutron Source (SNS) at Oak Ridge, USA Scholberg . 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 Scholb_IoP-2010 .

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 ORLaND ; Avignone03 ; Burman-03 . 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 Refs. Papoul-Kosm ; Papoul_tsk_AHEP-2015 and references therein). The present article is an extension of our previous calculations performed in Ref. TSAK-prc1 ; tsak_prc2 ; Tsak_AHEP 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 TSAK-prc1 ; kos-eric-11 ; tsak_prc2 ; Tsak_AHEP .

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.

II Brief Description of the Formalism

II.1 Angle differential coherent $\nu$-nucleus cross section

In the description of the $\nu$-nucleus scattering, the angle differential cross section $d\sigma/d\Omega$ is a useful quantity. For the elastic-scattering of a neutrino with energy $\varepsilon_{\nu}$ on a nucleus (A,Z) the angle differential cross section (with respect to the the scattering angle $\vartheta$) is defined as Engel-91 ; De-Vries ; Drukier ; Papoul-Kosm

[TABLE]

($G_{F}=1.1664\times 10^{-5}$ GeV*-2* is the Fermi weak coupling constant). In this definition, the quantity ${\cal F}(q^{2})$ includes the nuclear structure dependence of the cross section as Kosmas01 ; Kosmas-94

[TABLE]

where ${\Theta_{w}}$ denotes the weak mixing angle, known as Weinberg angle which takes the value sin${}^{2}\Theta_{w}\approx$ 0.2313. In Eqs. (1) and (2), $Q_{w}$ denotes the weak charge of the target nucleus given by

[TABLE]

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 Kolbe-Lag-Pin-03 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 $\nu$-nucleus scattering and the possibility to search for non-standard neutrino physics by taking advantage of the flavour-blind nature of the process Papoul_tsk_AHEP-2015 ; Horowitz-PRC-2012 .

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

[TABLE]

and are normalized as $F_{N,Z}(q^{2}=0)=1$. In the latter equation, $\rho_{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^{-(qb_{N})^{2}/6}$ with $b_{N}\approx 0.8$ fm being the nucleon harmonic oscillator size parameter tsk-NPA-92 ). The proton density $\rho_{p}(r)$ is often taken from experiment whenever measured charge densities are available Kosmas-94 ; Engel-91 .

Moreover, assuming that $F_{N}\approx F_{Z}$, from Eqs. (1) and (2) (in nuclei with $J^{\pi}=0^{+}$ ground state), one obtains

[TABLE]

where $f_{V}^{p}$ and $f_{V}^{n}$ stand for the polar-vector couplings of the weak neutral current

[TABLE]

Thus, the coherent cross section depends on the square of the ground-state nuclear form factor ${\cal F}(q^{2})$ at momentum transfer $q$ given by

[TABLE]

From Eq. (2), we see that, since $f_{V}^{p}=(1-4\sin^{2}\Theta_{w})/2\approx 0.0374$ is small, a neutrino scattered elastically on a spin-zero nucleus couples mostly to the neutron distribution, $\rho_{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})$ (Horowitz-PRC-2012, ; Chiang-Oset-Kosmas, ). Some authors consider that this would be complementary to parity violating experiments (Horowitz-PRC-2012, ; De-Vries, ) because it would provide additional data, obtained at different energy ranges and with different nuclei that could be used to calibrate nuclear structure calculations Kosmas01 ; Kosmas-94 ; Engel-91 ; De-Vries ; Drukier .

In earlier astrophysical estimations of the coherent scattering cross sections within the Standard Model (SM) Drukier ; Bachall (also in recent beyond the SM calculations Omar-1 ; Omar-2 ), the approximation $F_{N}(q^{2})\approx F_{Z}(q^{2})\approx 1$ was used for the total coherent cross section $\sigma_{tot}(\varepsilon_{\nu})$ written as

[TABLE]

(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 Avignone03 ; JDV-Avig , 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 Giom-JDV ; JDV-Giomat ; JDV-Avig ; Vogel-Engel

[TABLE]

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. Scholberg ; Scholb_IoP-2010 ; Giom-JDV ; JDV-Giomat ; JDV-Avig . 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 de-excitation product Tsak_AHEP .

The total coherent cross section $\sigma_{tot}(\varepsilon_{\nu})$ is obtained by integrating numerically Eq. (5) over the angle $\theta$ ($\theta_{min}=0$ to $\theta_{max}=\pi$) or Eq. (9) over $T_{A}$ between

[TABLE]

to $T_{A}^{max}=\infty$ Papoul-Kosm ; Vogel-Engel ; Kolbe-96 .

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. Papoul_tsk_AHEP-2015 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 $\nu_{\ell}$ or anti-neutrinos $\widetilde{\nu}_{\ell}$ ($\ell=e,\mu,\tau$) 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

[TABLE]

(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 $\ell=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 Lang-Pin-03 ; Athar-Sing ; Juoda-Lang ; Lang-08 ; Hax-87 .

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) DonPe ; TSAK-prc1 ; tsak_prc2 ; Don-Wal72 ; Don-Wal73 ; Kolbe-96 . As can be seen from Eq. (5), the original cross section for scattering of neutrinos, $\nu_{l}$ or anti-neutrinos, $\widetilde{\nu}_{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 $\nu$-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 $|J_{i}^{\pi_{i}}\rangle=|0^{+}_{gs}\rangle$) is computed by solving iteratively the BCS equations TSAK-prc1 ; tsak_prc2 ; prc3 ; Papoul-Kosm ; Tsak_AHEP .

In Table 1, we list the values of the resulting pairing parameters ($g_{pair}^{p,n}$) and the (theoretical) energy gaps ($\Delta^{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, $\Delta^{th}_{p,n}$, to the empirical ones $\Delta^{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 tsak_prc2 ; Tsak_AHEP . The values of the $g_{pair}^{p,n}$ adjust reliably the empirical energy gaps (see Table 1) TSAK-prc1 ; tsak_prc2 ; prc3 ; Bal-Ydr-Kos-11 ; Tsak_AHEP .

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

[TABLE]

($V_{j}^{k}$ denotes the probability amplitude for proton or neutron occupancies of the single particle $(n\ell)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_{j}^{k}=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+\ell$).

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,116Cd) and Fig. 1(d) shows the momentum dependence of ${\cal F}(q^{2})$ that enters Eqs. (1) and (5).

It should be noted that, the corrections due to the nucleon finite size ($e^{-(qb_{N})^{2}/6}$) and the nuclear center-of-mass motion ($e^{(qb)^{2}/4A}$), which enter as an overall q-dependent 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% tsk-NPA-92 .

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^{\pi}=0^{+}$). The differences, are mostly due to the small ratio ($\Delta N_{i}/N\approx 3\%-6\%$) in their neutron number.

Figure 2 illustrates the total integrated coherent cross sections of $\nu$-112,114,116Cd scattering as a function of (i) the momentum transfer $q$, Fig. 2(a), and (ii) the incoming neutrino energy $\varepsilon_{\nu}$, 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 ${\cal F}(q^{2})$, see Fig. 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\leq q\leq 2fm^{-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% De-Vries . 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% Horowitz-PRC-2012 ; Chiang-Oset-Kosmas .

In the next section, we summarize the main features of the $\nu$-energy distributions employed in this work for obtaining folded neutrino-nucleus cross sections for each $\nu$-source.

IV Energy-spectra of low-energy and Intermediate $\nu$-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 $\nu$-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 $\varepsilon_{\nu_{\mu}}=29.65$ MeV, see Fig. 3(d panel). For the non-mono-energetic neutrino fluxes we define the energy distributions $\eta(\varepsilon_{\nu})$ as

[TABLE]

($N_{\nu}$ represents the number of neutrinos of the beam).

We note that, via these energy spectra $\eta(\varepsilon_{\nu})$ of the specific neutrino sources, the original $\nu$-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 $\nu$-energy distribution of the studied neutrino source.

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

IV.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 $\varepsilon_{\nu}\preceq$ 10 MeV Vogel-Beacom ; Dye-2011 ; Fiorentini-2003 ; Fiorentini-2010 . Accurate measurements of the flux of these neutrinos KamLAND-11 ; Borexino-13 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 KamLAND ; Borexino-Col ; Borexino 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 $\widetilde{\nu}_{e}$), are generated through $\beta$-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_{\beta}$ as Vogel-Beacom

[TABLE]

A and Z denote the mass and atomic (proton) number, respectively, of the initial (parent) nucleus. Part of the decay energy $Q_{\beta}$ is carried away by anti-neutrinos ($Q_{\nu}$) while the remainder is available for heating ($Q_{h}$). Thus, $Q_{\beta}=Q_{\nu}+Q_{h}$.

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}Th$ decay series, and (iii) ${}^{40}K$ isotope Vogel-Beacom ; Dye-2011 . 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 $\beta$-decays of all intermediate radioactive isotopes.

Figure 3(a panel), shows the individual anti-neutrino spectra from ${}^{40}K$, ${}^{238}U$ series and ${}^{232}Th$ series ($\tau_{1/2}=4.47\times 10^{9}$ y, $\tau_{1/2}=14.0\times 10^{9}$ y and $\tau_{1/2}=1.28\times 10^{9}$ y, respectively). Essentially, these anti-neutrino ($\widetilde{\nu}_{e}$) energy spectra come from 82 beta decays in the U series and 70 beta decays in the $Th$ series Vogel-Beacom ; Dye-2011 ; Fiorentini-2003 ; Fiorentini-2010 ; Anderson-2012 ; Gouvea-neutrinos .

IV.2 Solar neutrinos

The solar neutrino spectra (mainly $\nu_{e}$ neutrinos) are produced through thermonuclear reactions taking place in the interior of the Sun Bachall-Holstein ; Bachall-Urlich ; Bachall-05a . The shape of the energy distribution (0.1 MeV$\leq\varepsilon_{\nu}\leq 18$ MeV) depends on the densities and temperatures in the Sun’s environment Bachall-Urlich and the individual process of the reaction chain (p-p neutrinos, 7Be 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$ Bachall-Holstein and $hep$ Bachall-Urlich ; Bachall neutrino sources predicted by the standard solar model Bachall . The ${}^{8}B$ $\nu$-spectrum, is nearly symmetric, with a peak at 6.4 MeV while the $hep$ spectrum is peaked at 9.6 MeV Bachall .

The detection of the solar neutrinos (produced either via the pp-chain reactions or via the CNO-cycle processes) by terrestrial experiments (SNO+ Zub-SNO+ ; Borexino ), constitutes excellent probes for astrophysics, nuclear physics, and particle physics searches Bachall-Urlich ; Bachall-05a . 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. KamLAND ; Borexino-Col ; Borexino ; Zub-SNO+ ; LENA .

IV.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 $\nu$-beam peaked at $\varepsilon_{\nu_{\mu}}$= 29.65 MeV), $\widetilde{\nu}_{\mu}$ and $\nu_{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. KamLAND ; Borexino-Col ; Borexino ; Zub-SNO+ ; LENA .

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 ORLaND ; Avignone03 ; Burman-03 ; Volpe07 . 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), $\nu_{e}$ neutrinos, and ${\tilde{\nu}}_{\mu}$ anti-neutrinos are produced from the decay of muons according to the reaction

[TABLE]

The decaying muons result from the decay of pions at rest ($\pi^{+}\to\mu^{+}+{\nu}_{\mu}$). Thus, these neutrino beams are not completely pure as, for example, the $\beta$-beam neutrinos bet-beam-Zuc ; Volpe07 . The energy-spectra of $\nu_{e}$ and ${\tilde{\nu}}_{\mu}$ neutrinos are fitted with the normalized distributions BOONE ; Anderson-2012

[TABLE]

[TABLE]

see Fig. 3(d panel), where $M_{\mu}=105.6$ MeV, is the muon rest mass. The ${{\tilde{\nu}}_{\mu}}$ spectrum is peaked at $\varepsilon_{\nu}^{max}=52.8$ MeV $=M_{\mu}/2$ while that of ${\nu_{e}}$ is peaked at $\varepsilon_{\nu}^{max}=35.2$ MeV $=M_{\mu}/3$ BOONE ; kolbe-kos .

Obviously, the analytic expressions of Eqs. (15) and (16), are convenient for the required integrations in the folding procedure, see below TSAK-prc1 ; tsak_prc2 ; prc3 ; Tsak_AHEP . On the other hand, their energy range and shape roughly resembles that of SN neutrinos.

IV.4 Reactor Neutrino spectra

The fission of very heavy nuclear isotopes ${}^{235}U$, ${}^{239}Pu$, and ${}^{238}U$ in the nuclear reactors produces a great number of neutron rich nuclear isotopes. Because these products are unstable, they decay via $\beta$-decay emitting anti-neutrinos ($\widetilde{\nu}_{e}$) Davis-Vogel-79 ; Bugey_experiment . Hence, nuclear reactors, operate as intense $\widetilde{\nu}_{e}$ sources for many experiments, giving fluxes of the order of $\sim 10^{13}\,\widetilde{\nu}/$cm2 s, at distances $\sim$ 10 m from the reactor core.

The energy spectrum of these anti-neutrinos, characteristic of the $\beta^{-}$ decay spectrum, is peaked at very low energies $\sim$ 0.3 MeV and covers the energy region below $\sim$ 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}Pu$ and 8$\%$ ${}^{238}U$ Davis-Vogel-79 ; Tengblad-89 .

Currently operating reactor neutrino experiments, like the TEXONO experiment in Taiwan TEXONO-Kerman ; TEXONO-Sevda , the MINER experiment at the Nuclear Science Center, Texas A&M University (using neutrinos from the TRIGA reactor) MINER-Agnolet , are excellent probes of beyond the standard model neutrino physics searches (electromagnetic $\nu$-properties) and coherent $\nu$-nucleus scattering studies.

IV.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 DonPe ; Kos-Ose ; Ejiri-PR00 . Multiple physics signatures are expected from a core-collapse explosion in the next supernova observation Burman-03 ; Kolbe-Lag-Pin-03 ; Fiorentini-2003 ; Fiorentini-2010 . 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 Lang-Pin-03 ; Bethe90 ; Lank-prl ; Janka95 ; Janka07 . 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 $\nu$-detectors to SN neutrinos ORLaND ; Avignone03 ; Burman-03 . 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 Janka89 or through the average $\nu$-energy in the analytically simpler two parameter Power-Law (PL) distribution (see Appendix) Janka89 ; Raffelt03 ; Raffelt-proc03 .

Both parametrizations, FD and PL, yield similar distributions characterized by the temperature $T$ or the average $\nu$-energy $\langle\varepsilon_{\nu}\rangle$ tsak_prc2 ; prc3 ; Duan08 ; Jachowicz06 ; Vtsak-Zn-fold ; TSAK10-erice2-2010 . These analytic normalised expressions contain two parameters to include modulation effects due to various corrections required to modify the purely thermal shape initially employed Janka89 ; Raffelt03 ; Jachowicz06 . 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}=\mu/T$ which makes the spectrum more narrow compared to the purely thermal shape of temperature $T$ (in MeV) tsak_prc2 . The two parameter PL distribution of SN-$\nu$ energy spectrum Raffelt03 ; Raffelt-proc03 , contains as parameters the mean neutrino energy $\langle\varepsilon_{\nu}\rangle$ and the parameter $\alpha$ which adjusts the width $w$ of the distribution tsak_prc2 ; Janka89 ; Raffelt03 ; Jachowicz06 (see Appendix).

In Fig. 4, some flavour dependent $\nu$-energy spectra ($\eta_{SN}$) emitted by a core-collapse Supernova, needed for our present work, are illustrated. Both FD and PL energy distributions (labelled $\eta_{FD}$ and $\eta_{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 $\nu$-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 $\varepsilon_{\nu}$ and its peak becomes smaller. Furthermore, the degeneracy parameter shifts the spectrum t o higher energies tsak_prc2 ; Jachowicz06 . In this figure, the PL energy distributions for the corresponding values of mean neutrino energy $\langle\varepsilon_{\nu}\rangle$, are also illustrated ($\langle\varepsilon_{\nu}\rangle$ reflects the depth of the stars from which the neutrinos are escaping. We see that, as the $\langle\varepsilon_{\nu}\rangle$ grows, the maximum of the distribution shifts to higher $\nu$-energy $\varepsilon_{\nu}$ Raffelt03 ; Jachowicz06 .

In Table 2, the corresponding values of parameters for the equivalent FD and PL $\nu$-energy spectra of Fig. 4, that have been employed in various SN scenarios are shown (for more details see the Appendix and Ref. tsak_prc2 ). It is worth mentioning that, due to neutrino oscillations and other phenomena, at any distance from the source the SN-$\nu$ spectra can be different compared to those originally produced at the core of the collapsing star. It is, however, expected that $\nu$-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 $\widetilde{\nu}_{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 Scholberg-Exp-Statis . Among those, are the next generation detectors HyperKamiokande, Juno, Dune, etc., which aim at measuring, among others, the diffuse SN neutrino background DSNB-Vissani ; DSNB-Moller .

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 Zuber-plb-1 ; Zuber-plb-2 ; Zuber-PPNP-10 . Our aim in this section is to simulate some of these features by calculating convoluted cross sections as discussed in Refs. tsak_prc2 ; Tsak_AHEP . .

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, $\langle\sigma_{tot}\rangle$ and, then, corresponding supernova neutrino event rates and fluxes.

V.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 $\langle\sigma_{coh}\rangle$ is obtained through the folding Kos-Ose ; TSAK-prc1 ; tsak_prc2

[TABLE]

For a CdTe or CdZnTe detector material, the flux averaged cross sections, computed by inserting in Eq. (17) the $\sigma_{coh}(\varepsilon_{\nu})$ from Fig. 2(b) and the $\eta(\varepsilon_{\nu})$ from Figs. 3 and 4, for the isotopes 112,114,116Cd, are listed in Tables 3 and 4 as described below.

In Table 3 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 8B and the $hep$ solar neutrinos).

In the last two columns of Table 4 we tabulate the $\langle\sigma_{coh}\rangle$ calculated for the distributions of Eqs. (15) and (16), i.e. the $\nu$-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 2.

In supernova neutrino scenarios, usually average $\nu$-energies between 10 $\leq\langle\varepsilon_{\nu}\rangle\leq$ 12 MeV are employed for the description of $\nu_{e}$ neutrinos, average energies between $15\leq\langle\varepsilon_{\nu}\rangle\leq 18$ MeV for $\widetilde{\nu}_{e}$ anti-neutrinos, and average energies between 22 $\leq\langle\varepsilon_{\nu}\rangle\leq$ 26 MeV for $\nu_{x}$ and $\widetilde{\nu}_{x}$, with x = $\mu$, $\tau$ Janka89 ; Raffelt03 ; Raffelt-proc03 ; Jachowicz06 .

Due to the dominance of the coherent channel throughout the region of the incoming neutrino energy $\varepsilon_{\nu}$ of our present calculations, the flux averaged coherent cross section $\langle\sigma_{coh}\rangle$ may be even two or three orders of magnitude larger than the total incoherent cross section $\langle\sigma_{tot}^{incoh}\rangle$ TSAK-prc1 ; tsak_prc2 ; prc3 .

V.2 Number of events in $\nu$-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 Zuber-plb-1 ; Zuber-plb-2 ; Zuber-PPNP-10 , as follows. By using the flux averaged cross sections $\langle\sigma(\varepsilon_{\nu})\rangle$ of Table 4, for instance those referred to the SN neutrinos of the 112,114,116Cd isotopes, we estimate (potentially detectable) neutrino fluxes $\Phi_{\nu}$ 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 $\Phi_{\nu}$ reaching the nuclear detector with the expression TSAK-prc1 ; Bal-Ydr-Kos-11 ; Tsak_AHEP

[TABLE]

We note that, experimentalists use the definition

[TABLE]

which takes into account the detection efficiency $\epsilon$ (usually equal to $\epsilon\approx$ 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 Zuber-plb-1 ; Zuber-plb-2 ; Zuber-PPNP-10 .

In the first step, we choose three SN neutrino scenarios in which the mean energies are: (i) $\langle\varepsilon_{\nu}\rangle=12$ MeV (corresponding to SN electron neutrinos ${\nu}_{e}$), (ii) $\langle\varepsilon_{\nu}\rangle=16$ MeV (corresponding to SN electron anti-neutrinos $\widetilde{\nu}_{e}$), and (iii) $\langle\varepsilon_{\nu}\rangle=24$ MeV (corresponding to SN $\nu_{x}$, $\widetilde{\nu}_{x}$, with $x=\mu,\tau$ (anti)neutrinos of heavy leptons).

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 114Cd isotope or equivalently a number of 114Cd atoms (nuclei) equal to $N_{Cd}\equiv N_{{}^{114}Cd}=94.17N_{Avog}$.

In Eq. (18), as total neutrino scattering cross sections, $\sigma_{tot}(\varepsilon_{\nu})$ we employ the values of flux averaged cross sections $\langle\sigma_{coh}\rangle$ 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 $N_{ev}$ as: (a) $N_{ev}$=1 event s*-1*=3.15$\times 10^{7}$ events y*-1*, (b) $N_{ev}$=1 event hr*-1*=8.76$\times 10^{3}$ events y*-1*, (c) $N_{ev}$=1 event d*-1*=3.65$\times 10^{2}$ events y*-1*, and (d) $N_{ev}$= 12 events y*-1* and (from Eq. (18) we compute the corresponding SN $\nu$ fluxes $\Phi_{\nu}$.

In a similar way, assuming that the COBRA detector contains 100 kg of the material CdTe, we find 13.5 kg 114Cd or about $N_{Cd}\equiv N_{{}^{114}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 $\Phi_{\nu}$ 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.).

The above neutrino fluxes are of the same order with those of the Spallation Neutron Source (SNS) at ORLaND, Oak Ridge ORLaND ; Avignone03 ; Burman-03 . We mention that the COHERENT experiment at Oak Ridge, with a 14.57 kg of CsI scintillator detector, by using an SNS $\nu_{\mu}$ neutrino flux (coming from $\pi$-decay at rest) as high as $\Phi_{\nu}^{COH}=1.7\times 10^{11}\nu_{\mu}/{\textrm{c}m}^{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) CEnN_neut . These results translate to event rate $N_{ev}^{COH}$ = 168/y $\nu_{\mu}$ neutrinos.

From the results of Table 5, we may define the ratio $N_{ev}/\Phi_{\nu}$ for the COHERENT experiment ($R^{COH}$) and for a special $\nu_{\mu}$ neutrino case of the COBRA experiment ($R^{COB}$). 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 $\nu_{\mu}$ neutrinos). From these two ratios we find that $R=R^{COH}/R^{COB}=98.95/2.94\approx 34$, which means that, for the chosen SN neutrino scenario, the COBRA detector may observe 12 $\nu_{\mu}/y$ only if its mass is equal to m $\sim$ 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 $\nu_{\mu}$ 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 $\epsilon$ has not been considered (equivalently we assumed $\epsilon=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,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 $\nu$-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, $\langle\sigma_{coh}\rangle$, reflect the mean neutrino signals generated in several terrestrial detectors (112,114,116Cd) from such $\nu$-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”.

Appendix

VII.1 Three-point formulas for empirical energy gaps $\Delta^{exp}_{n,p}$ of neutrons and protons

The empirical energy gaps for neutrons, $\Delta^{exp}_{n}$, and protons, $\Delta^{exp}_{n}$, 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, $S_{n}$ or protons, $S_{p}$ of the isotope $(A,Z)$ and also those of the neighbouring nuclear isotopes with $N\pm 1$ neutrons or $Z\pm 1$ protons, respectively, by employing the expressions

[TABLE]

The above equations are known as the three-point formulas (see, e.g. Ref. Bohr-Motel ; prc3 ).

VII.2 Normalization of the $\nu$ energy distributions $\eta(\varepsilon_{\nu})$ adopted in this work

The distributions $\eta(\varepsilon_{\nu})$ adopted in the present work (see Sect. IV), are considered to be normalized in such a way that

[TABLE]

For example, in the case of $\eta_{\nu_{e}}(\varepsilon_{\nu})$ of Fig. 3(d), the normalization gives

[TABLE]

where we have used $\varepsilon_{\nu}^{min}=0$ and $\varepsilon_{\nu}^{max}={M_{\mu}/2}$

VII.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 $\langle\varepsilon_{\nu}\rangle$ and the width $w$ of the spectrum is defined as

[TABLE]

where $w_{0}=\varepsilon_{\nu}/\sqrt{3}$ is the width of the identical FD and PL distributions tsak_prc2 .

VII.3.1 Fermi-Dirac (FD) energy distribution

By introducing the degeneracy parameter $n_{dg}$ (equal to the ratio of the chemical potential $\mu$ divided by the neutrino temperature $T$, i.e. $n_{dg}=\mu/T$), the Fermi-Dirac energy distribution reads

[TABLE]

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

[TABLE]

Inserting Eq. (23) into Eq. (22), we take

[TABLE]

VII.3.2 Power-law energy distribution

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

[TABLE]

where $\langle\varepsilon_{\nu}\rangle$ is the neutrino mean energy. The parameter $\alpha$ adjusts the width of the spectrum (see text). The normalization factor $C$, is calculated from the normalization condition

[TABLE]

From the later equation we find

[TABLE]

therefore, Eq. (25) becomes

[TABLE]