Abstract

Event detection rates for WIMP-nucleus interactions are calculated for Ga, Ge, As, and I (direct dark matter detectors). The nuclear structure form factors, which are rather independent of the underlying beyond the Standard Model particle physics scenario assumed, are evaluated within the context of the deformed nuclear shell model (DSM) based on Hartree-Fock nuclear states. Along with the previously published DSM results for Ge, the neutrino-floor due to coherent elastic neutrino-nucleus scattering (CENS), an important source of background to dark matter searches, is extensively calculated. The impact of new contributions to CENS due to neutrino magnetic moments and mediators at direct dark matter detection experiments is also examined and discussed. The results show that the neutrino-floor constitutes a crucial source of background events for multi-ton scale detectors with sub-keV capabilities.

1. Introduction

In the last few decades, the measurements of the cosmic microwave background (CMB) radiation offered a remarkably powerful way of modelling the origin of cosmic-ray anisotropies and constraining the geometry, the evolution, and the matter content of our universe. Such observations have in general indicated the consistency of the standard cosmological model [1] and the fact that our universe hardly contains % luminous matter, whereas the remainder consists of nonluminous dark matter (%) and dark energy (%) [2]. After the discovery of the CMB fluctuations by the Cosmic Background Explorer (COBE) satellite [3], the extremely high precision of the WMAP satellite and especially of the Planck third-generation space mission has helped us to produce maps for the CMB anisotropies and other cosmological parameters (see [4] for details). We also mention that high-resolution ground-based CMB data, like those of the Atacama Cosmology Telescope (ACT) [5] and the South Pole Telescope (SPT) [6] have recovered the underlying CMB spectra observed by the space missions.

Focusing on the topic we address in this work, it is worth noting that the CMB data, the Supernova Cosmology project [7], and so on suggest that most of the dark matter of the universe is cold. Furthermore, the baryonic cold dark matter (CDM) component can be considered to consist of either massive compact halo objects (MACHOs) like neutron stars, white dwarfs, Jupiters, etc. or Weakly Interacting Massive Particles (WIMPs) that constantly bombard Earth’s atmosphere. Several results of experimental searches suggest that the MACHO fraction should not exceed a portion of about 20% [1]. On the theoretical side, within the framework of new physics beyond the Standard Model (SM), supersymmetric (SUSY) theories provide promising nonbaryonic candidates for dark matter [8] (for a review see [9]). In the simple picture, the dark matter in the galactic halo is assumed to be Weakly Interacting Massive Particles (WIMPs). The most appealing WIMP candidate for nonbaryonic CDM is the lightest supersymmetric particle (LSP) which is expected to be a neutral Majorana fermion traveling with nonrelativistic velocities [10].

In recent years, there have been considerable theoretical and experimental efforts towards WIMP detection through several nuclear probes [1113]. Popular target nuclei include among others the Ga, Ge, As, I, Xe, and Pb isotopes [14, 15]. Towards the first ever dark matter detection, a great number of experimental efforts take place aiming at measuring the energy deposited after the galactic halo WIMPs scatter off the nuclear isotopes of the detection material. Because of the low count rates, due to the fact that the WIMP-nucleus interaction is remarkably weak, the choice of the detector plays very important role and for this reason spin-dependent interactions require the use of targets with nonzero spin. The Cryogenic Dark Matter Search (CDMS) experimental facility [16] has been designed to directly detect the dark matter by employing a Ge as the target nucleus, setting the most sensitive limits on the interaction of WIMP with terrestrial materials to date. The development of upgrades is under way and will be located at SNOLAB. Another prominent dark matter experiment is the EDELWEISS facility in France [17] which uses high purity germanium cryogenic bolometers at milli-Kelvin temperatures. There are also other experimental attempts using detectors like I, Xe, Cs, etc. (see [1820]).

Inevitably, direct detection experiments are exposed to various neutrino emissions, such as those originating from astrophysical sources (e.g., Solar [21], Atmospheric [22, 23], and diffuse Supernova [24] neutrinos), Earth neutrinos (Geoneutrinos [25]), and in other cases even artificial terrestrial neutrinos (e.g., neutrinos from nearby reactors [26]). The subsequent neutrino interactions with the material of dark matter detectors, namely, the neutrino-floor, may perfectly mimic possible WIMP signals [27]. Thus, the impacts of the neutrino-floor on the relevant experiments looking for CDM as well as on the detector responses to neutrino interactions need be comprehensively investigated. Since Geoneutrino fluxes are relatively low, astrophysical neutrinos are recognised as the most significant background source that remains practically irreducible [28]. The recent advances of direct detection dark matter experiments, mainly due to the development of low threshold technology and high detection efficiency, are expected to reach the sensitivity frontiers in which astrophysical neutrino-induced backgrounds are expected to limit the observation potential of the WIMP signal [29].

In this work, we explore the impact of the most important neutrino background source on the relevant direct dark matter detection experiments by concentrating on the dominant neutrino-matter interaction channel, e.g., the coherent elastic neutrino-nucleus scattering (CENS) [30, 31]. It is worthwhile to mention that events of this process were recently measured for the first time by the COHERENT experiment at the Spallation Neutrino Source [32], completing the SM picture of electroweak interactions at low energies. Such a profound discovery motivated our present work and we will make an effort to shed light on the nuclear physics aspects. Neutrino nonstandard interactions (NSIs) [33] may constitute an important source of neutrino background and have been investigated recently in [34, 35]. Thus, apart from addressing the SM contributions to CENS [36], we also explore the impact of new physics contributions that arise in the context of electromagnetic (EM) neutrino properties [37, 38] as well as of those emerging in the framework of gauge interactions [39] due to the presence of new light mediators [40, 41]. The aforementioned interaction channels may lead to a novel neutrino-floor as demonstrated by [42]. The latter could be detectable in view of the constantly increasing sensitivity of the upcoming direct detection experiments with multi-ton mass scale and sub-keV capabilities [43].

Direct detection dark matter experiments are currently entering a precision era, and nuclear structure effects are expected to become rather important and should be incorporated in astroparticle physics applications [45]. For this reason, our nuclear model is at first tested in its capabilities to adequately describe the nuclear properties before being applied to problems like dark matter detection. This work considers the deformed shell model (DSM), on the basis of Hartree-Fock (HF) deformed intrinsic states with angular momentum projection and band mixing [46], all with a realistic effective interaction and a set of single-particle states and single-particle energies, which is established to be rather successful in describing the properties of nuclei in the mass range = 60–90 (see [47] for details regarding DSM and its applications). In particular, the DSM is employed for calculating the required nuclear structure factors entering the dark matter and neutrino-floor expected event rates by focusing on four interesting nuclei regarding dark matter investigations such as Ga, Ge, As, and I. Let us add that details of nuclear structure and dark matter event rates for Ge obtained using DSM have been reported recently [48].

The paper has been organised as follows. Section 2 gives the main ingredients of WIMP-nucleus scattering, while Section 3 provides the formulation for neutrino-nucleus scattering (neutrino-floor) within and beyond the SM. Then in Section 4 we describe briefly the methodology of the DSM, and the main results of the present work are presented and discussed in Section 5. Finally, the concluding remarks are drawn in Section 6.

2. Searching WIMP Dark Matter

The Earth is exposed to a huge number of WIMPs originating from the galactic halo. Their direct detection through nuclear recoil measurements after scattering off the target nuclei at the relevant dark matter experiments is of fundamental interest in modern physics and is expected to have a direct impact on astroparticle physics and cosmology. In this section we discuss the mathematical formulation of WIMP-nucleus scattering. The formalism introduces an appropriate separation of the SUSY and nuclear parts entering the event rates of WIMP-nucleus interactions in our effort to emphasise the important role played by the nuclear physics aspects. In particular we perform reliable nuclear structure calculations within the context of DSM based on Hartree-Fock states.

2.1. WIMP-Nucleus Scattering

For direct detection dark matter experiments, the differential event rate of a WIMP with mass scattering off a nucleus with respect to the momentum transfer can be cast in the form [1]where denotes the number of target nuclei per unit mass, stands for the mass number of the target nucleus, and is the proton mass. In the above expression the WIMP flux is , with being the local WIMP density. The distribution of WIMP velocity relative to the detector (or Earth) and also the motion of the Sun and Earth, , is taken into account and assumed to resemble a Maxwell-Boltzmann distribution to ensure consistency with the LSP velocity distribution. Note that, by neglecting the rotation of Earth in its own axis, accounts for the relative velocity of WIMP with respect to the detector. For later convenience a dimensionless variable is introduced with denoting the oscillator length parameter, and the corresponding WIMP-nucleus differential cross section in the laboratory frame reads [10, 15, 4850]withThe first three terms account for the spin contribution due to the axial current, while the fourth term accounts for the coherent contribution arising from the scalar interaction. The coherent contribution is expressed in terms of the nuclear form factors given asThe coherent part in the approximation of nearly equal proton and neutron nuclear form factors is given asThe respective values of the nucleonic-current parameters , for the isoscalar and isovector parts of the vector current (not shown here), , for the isoscalar and isovector parts of the axial-vector current, and , for the isoscalar and isovector parts of the scalar current depend on the specific SUSY model employed [51]. The spin structure functions with for the isoscalar and isovector contributions, respectively, take the formwithHere, and with for protons and for neutrons, while represents the solid angle for the position vector of the -th nucleon and stands for the well-known spherical Bessel function. The quantities are the static spin matrix elements (see, e.g., [8]). In this context, the WIMP-nucleus event rate per unit mass of the detector is conveniently written asThe functions enter the definition of the WIMP-nucleus event rate through the three-dimensional integrals, given bywithIn the latter expression, , , and account for the spin-dependent parts of (3), while is associated with the coherent contribution.

In this work, the nuclear wave functions and entering (7) are calculated within the nuclear DSM of [46, 47]. For a comprehensive discussion on the explicit form of the function , the integration limits of (9) and the various parameters entering into these, the reader is referred to [48].

3. Neutrino-Nucleus Scattering

The neutrino-floor stands out as an important source of irreducible background to WIMP searches at a direct detection experiment. In this work we explore the neutrino-floor due to neutrino-nucleus scattering since the corresponding floor coming from neutrino-electron scattering is relatively low [52]. Motivated by the novel neutrino interaction searches using reactor neutrinos of [42], here we consider various astrophysical neutrino sources in our calculations that involve the conventional and beyond the SM interactions channels (see below).

3.1. Differential Event Rate at Dark Matter Detectors

For a given interaction channel , the differential event rate of CENS processes at a dark matter detector is obtained through the convolution of the normalised neutrino energy distribution of the background neutrino source in question (i.e., Solar, Atmospheric and Diffuse Supernova Neutrinos, as seen below) with the CENS cross section, as follows [53]:where is the maximum neutrino energy of the source in question (for the case of Solar neutrinos see, e.g., Table 1) and is the minimum neutrino energy that is required to yield a nuclear recoil with energy . In the latter expression with being the exposure time, is the number of target nuclei and is the assumed neutrino flux.

3.1.1. Standard Model Interactions

Assuming SM interactions only, at low and intermediate neutrino energies , the weak neutral-current CENS process is adequately described by the four-fermion effective interaction Lagrangian [33, 36]where denote the chiral projectors, represents the neutrino flavour, and is a first generation quark. By including the radiative corrections of [54], the -handed couplings of the quarks to the -boson are expressed aswith , , , , , and .

In this work we restrict our study only to low momentum transfer in order to satisfy the coherent condition , where is the nuclear size and is the magnitude of the three-momentum transfer [31]. Focusing on the dominant CENS channel, the relevant SM differential cross section with respect to the nuclear recoil energy takes the form [41]with denoting the neutrino energy and the mass of the target nucleus. The relevant vector () and axial-vector () weak charges entering the CENS cross section are given by the relations [55]Here, () stands for the number of protons (neutrons) with spin up (+) and spin down (−), respectively, while () represent the axial-vector couplings of protons (neutrons) to the boson. At the nuclear level, the relevant vector (axial-vector) couplings of protons () and neutrons () take the formThe axial-vector nucleon form factor takes into account the spin structure of the nucleon and is defined as [56]where is the free axial-vector coupling constant and the axial mass is taken to be GeV, while strange quark effects have been neglected.

We note that for spin-zero nuclei the axial-vector contribution vanishes, while for the odd- nuclei considered in the present study it is negligible and of the order of . The weak charges in (15) encode crucial information regarding the finite nuclear size through the proton and neutron nuclear form factors, which in our work are obtained within the context of the DSM (see below), as functions of the momentum transfer . Contrary to similar studies assuming the conventional Helm-type form factors, the present work also takes into account the nuclear effects due to the nonspherical symmetric nuclei employed in dark matter searches.

3.1.2. Electromagnetic Neutrino Contributions

Turning our attention to new physics phenomena we now address potential contributions to CENS in the framework of nontrivial neutrino EM interactions that may lead to a new neutrino-floor at low detector thresholds. In this framework, the presence of an effective neutrino magnetic moment leads to an EM contribution of the differential cross section that has been written as [41]Neglecting axial effects, the EM contribution to CENS at a direct detection dark matter is encoded in the factorwhere the relevant EM charge is written in terms of the electron mass , the fine-structure constant , and the effective neutrino magnetic moment as [57]In contrast to the dependence of the SM case, (19) and (20) imply the existence of a coherence along with a characteristic enhancement of the total cross section. This implies a potential distortion of the expected recoil spectrum at very low recoil energies that may be detectable at future direct dark matter detection with sub-keV operation thresholds.

For the sake of completeness, we stress that the effective neutrino magnetic moment is expressed through neutrino amplitudes of positive and negative helicity states, e.g., the 3-vectors and and the neutrino transition magnetic moment matrix, , in flavour basis, as [37, 58]Then, the effective neutrino magnetic moment is written in mass basis through a proper rotation; for a detailed description of this formalism see [59].

3.1.3. Novel Mediator Contribution

We now explore novel mediator fields that could be accommodated in the context of simplified scenarios [60, 61] predicting the existence of a new vector mediator with mass [62]. Such beyond the SM interactions may constitute a new neutrino-floor at direct detection dark matter experiments [39].

The presence of a mediator gives rise to subleading contributions to the SM CENS rate, described by the Lagrangian [63]where only left-handed neutrinos are assumed (right-handed neutrinos in the theory would lead to vector-axial-vector cancellations). The resulting cross section reads [41]with the factor being written in terms of the neutrino-vector coupling , asThe relevant charge in this case is expressed through the vector quark couplings to the boson, as [39] Let us mention that emerging degeneracies can be either reduced through multidetector measurements [61] or broken in the framework of NSIs [64]. For completeness we note that, despite being not present for the low energies considered here, these couplings could be changed by currently unknown in-medium effects (see, e.g., [23] and references therein).

3.2. Neutrino Sources
3.2.1. Solar Neutrinos

In terrestrial searches for dark matter candidates at low energies, the Solar neutrinos emanating from the interior of the Sun generated through various fusion reactions produce a dominant background for direct CDM detection experiments. Assuming WIMP masses less than 10 GeV, an estimated total Solar neutrino flux of about [65] hitting the Earth is expected to appreciably limit the sensitivity of such experiments [27]. On the other hand, the theoretical uncertainties of Solar neutrinos are presently quite large and depend strongly on the assumed Solar neutrino model. To maintain consistency with existing Solar data, in this work we consider the high metallicity Standard Solar Model (SSM) [21]. We note however that the dominant Solar neutrino component coming from the primary proton-proton channel ( neutrinos) that accounts for about 86% the Solar neutrinos flux has been recently measured by the Borexino experiment with an uncertainty of 1% [66]. Through CENS, the direct detection dark matter experiments are mainly sensitive to two sources of Solar neutrinos, namely, the B and the neutrinos which cover the highest energy range of the Solar neutrino spectrum. Since B neutrinos are generated from the decay , while neutrinos from , both sources occur in the aftermath of the chain. Following previous similar studies [28], in this work, we explore the neutrino-floor extending our analysis to the lowest neutrino energies, by considering the neutrino line which belongs to the chain and the -capture reaction on Be that leads to two monochromatic beams at 384.3 and 861.3 keV as well as the well-known CNO cycle. The latter neutrinos appear as three continuous spectra (N, O, and F) with end point energies close to the neutrinos.

3.2.2. Atmospheric Neutrinos

Atmospheric neutrinos are decay products of the particles (mostly pions and kaons) produced as a result of cosmic-ray scattering in the Earth’s atmosphere. The generated secondary particles decay to , , , and constituting a significant background to dark matter searches especially for WIMP masses above 100 GeV. In particular, the effect is crucial on the discovery potential of WIMPs with spin-independent cross section of the order of . The direct detection dark matter experiments, due to the lack of directional sensitivity, are in principle sensitive to the lowest energy (less than 100 MeV) atmospheric neutrinos. For this reason, in our present work atmospheric neutrinos are considered by employing the low-energy flux coming out of the FLUKA code simulations [22].

3.2.3. Diffuse Supernova Neutrinos

The weak glow of MeV neutrinos emitted from the total number of core-collapse supernovae, known as the Diffuse Supernova Neutrino Background (DSNB), creates an important source of neutrino background specifically for the WIMPs mass range 10–30 GeV [24]. Despite the appreciably lower flux compared to Solar neutrinos, DSNB neutrino energies are higher than those of the Solar neutrino spectrum. In our simulations, the adopted DSNB distributions (usually of Fermi-Dirac or power-law type) correspond to temperatures 3 MeV for , 5 MeV for , and 8 MeV for the other neutrino flavours denoted as or , .

Figure 1 shows the unoscillate neutrino flux considered in the present study, illustrating the Solar neutrino spectra of the dominant neutrino sources assuming the high metallicity Standard Solar Model (SSM) as defined in [21]. Also shown is the low-energy atmospheric neutrino flux as obtained from the FLUKA simulation [22] as well as the DSNB spectrum [67]. The corresponding neutrino types, maximum energies, and fluxes are listed in Table 1.

4. Deformed Shell Model

In the formalism of the WIMP-nucleus or neutrino-nucleus event rates of Sections 2 and 3, both for the case of elastic or inelastic interaction channels, the nuclear physics and particle physics (SUSY model) parts appear almost completely separated. In the present work our main focus drops on the nuclear physics aspects which are contained in the nuclear structure factors discussed in Section 2. Special attention is paid on the factors of (9) that depend on the spin structure functions and the nuclear form factors. These quantities have been calculated using the DSM method [48] (for a comprehensive discussion of DSM see [47]) given the kinematics and the assumptions describing the WIMP particles.

The construction of the many-body wave functions for the initial and final nuclear states in the framework of DSM involves performance of the following steps. (i) At first, one chooses a model space consisting of a given set of spherical single-particle (sp) orbits, sp energies, and the appropriate two-body effective interaction matrix elements. For Ga and As, the spherical sp orbits are , , , and with energies 0.0, 2.20, 2.28, and 5.40 MeV and 0.0, 0.78, 1.08, and 3.20 MeV, respectively, while the assumed effective interaction is the modified Kuo interaction [68]. Similarly for I, the sp orbits, their energies, and the effective interaction are taken from a recent paper [69]. (ii) Assuming axial symmetry and solving the HF single-particle equations self-consistently, the lowest-energy prolate (or oblate) intrinsic state for the nucleus in question is obtained. An example is shown in Figure 2 for Ga. (iii) The various excited intrinsic states then are obtained by making particle-hole (-) excitations over the lowest-energy intrinsic state (lowest configuration). (iv) Then, because the HF intrinsic nuclear states ( is azimuthal quantum number and distinguishes states with the same ) do not have definite angular momentum, angular momentum projected states are constructed asIn the previous expression, , , ) represents the Euler angles, denotes the known general rotation operator, and the Wigner -matrices are defined as . Here, is the normalisation constant which by assuming axial symmetry is defined aswhere the functions are the diagonal elements of the matrix . Finally, the good angular momentum states are orthonormalised by band mixing calculations and then, in terms of the index , it is possible to distinguish between different states having the same angular momentum ,Within the DSM method, for the evaluation of the reduced nuclear matrix element entering (6) and (7), we first calculate the single-particle matrix elements of the relevant operators , aswhere is the 9- symbol. For more details, the reader is referred to [7072]. It should be noted that in the DSM method one considers an adequate number of intrinsic states in the band mixing calculations.

DSM calculations are performed in the same spirit as in spherical shell model where one takes a model space and a suitable effective interaction (single-particle orbitals, single-particle energies, and a two-body effective interaction). This procedure has been found to be quite successful in describing the spectroscopic properties and electromagnetic properties of many nuclei in the mass region A = 60–90 and has also been applied to double beta decay nuclear transition matrix elements [47]. In addition, this model has been used recently in calculating the event rates for dark matter detection [48]. With the proper choice of effective interaction, one will not be considering core excitations. This is a standard prescription in shell model as well as in DSM. To go beyond this, one has to use no-core shell model or DSM with much larger set of single-particle orbitals (inclusion of core orbitals), such refinements are planned to be employed in future calculations.

We note that the many-body nuclear calculations performed take into account in the usual way the inert core orbits (completely filled by the protons and neutrons) and the extra-core nucleons moving in the assumed model space under the influence of an effective interaction. The explicit values of the occupied nucleon single-particle deformed orbits of the HF intrinsic states considered in our calculations are listed in Table 2.

5. Results and Discussion

5.1. Nuclear Physics Aspects

To maximise the significance of our WIMP-nucleus and neutrino-floor calculations, the reliability of the obtained nuclear wave functions is tested by comparing the extracted energy level spectrum and magnetic moments with available experimental data. The consistency of this method obtained for Ge has been already presented in [48]. Furthermore, in the DSM calculations for Ga and As, we restrict ourselves to prolate solutions only, since the oblate solution does not reproduce the energy spectra and electromagnetic properties of these nuclei. It also does not mix with the prolate solution. Hence, we neglect the oblate solutions in the calculations. For each of these nuclei, we consider only four intrinsic prolate states which should be sufficient to explain the systematics of the ground state and close lying excited state. Due to size restrictions, in Figure 3 we illustrate only the calculated spectrum for Ga.

For As, the ground state is and there are also two and levels around 0.12 MeV. In addition, there is a collective band consisting of , , and levels at 0.279, 1.095, 2.150, and 3.091 MeV, respectively. All these levels are well reproduced by the DSM method. Turning to the I spectrum, there are four observed collective bands with band heads , , and . There are evidences suggesting that low-lying states in I have oblate deformation [73]. Hence, for this nucleus, we consider only oblate configurations and take the six lowest oblate intrinsic states in the band mixing calculation. These intrinsic states are found to provide adequate description of the energy spectrum and electromagnetic properties for this nucleus. The calculations for this nucleus utilise a new effective interaction developed by an Italian group very recently [69]. The new effective interaction is seen to reproduce well the I spectrum; details will be presented elsewhere.

We thus conclude that concerning the evaluation of the WIMP-nucleus and CENS event rates we are interested in this work, the required ground state wave functions obtained through the DSM method are reliable and the intrinsic states used in the subsequent analysis are considered sufficient. From the perspective of nuclear physics, spin contributions constitute significant ingredients in the evaluation of WIMP-nucleus event rates. For this reason, the first stage of our work involves the calculation of the magnetic moment, which is decomposed into an orbital and spin part. The relevant results for the proton and neutron contributions to the orbital and spin parts concerning the ground states of the four nuclear isotopes studied in this paper are given in Table 3. A comparison between the obtained magnetic moments and the respective experimental data is also provided. Despite the fact that these calculations adopt bare values of -factors neglecting quenching effects, the obtained DSM results of the ground state magnetic moments are consistent with the experimental values.

Having successfully reproduced the energy spectrum and the magnetic moments within the context of the DSM wave functions, we evaluate important nuclear physics inputs entering the WIMP-nucleus and CENS cross sections. Figures 4 and 5 present a comparison between the DSM nuclear form factors and the effective Helm-type ones employed in various similar studies, where, as can be seen, the DSM results differ from the Helm-type ones. The behaviour of the proton form factor for Ga is found to be different from those of the other nuclei and this may be due to the nearby proton shell closure and the neutron subshell closure. Calculations with several different effective interactions are under way to rule out the possibility of any deficiency of the effective two-body interaction used. We furthermore illustrate the spin structure functions of WIMP-Ga elastic scattering calculated using (6) and (7). The variation of , , and with respect to the parameter is shown in Figure 6, while similar results are obtained for As and I (for the Ge case see [48]).

The consistency of our nuclear physics DSM calculations has been extensively explored in this work and compared with existing experimental data (see Figures 2 and 3 and Table 3) making the considered form factors reliable. Specifically we have tested the reliability of this model to describe nuclear structure properties such as excitation spectra and nuclear magnetic moments. We mention that DSM has been tested in the past in many nuclei in the = 60–90 region [47] (see above).

5.2. WIMP-Nucleus Rates and the Neutrino-Floor

The WIMP-nucleus event rates and the neutrino-floor due to neutrino-nucleus scattering are calculated for a set of interesting nuclear targets such as Ga, Ge, As, and I (see Table 3). In evaluating the neutrino-induced backgrounds, we consider only the dominant CENS channel, since neutrino-electron events are expected to produce less events by about one order of magnitude [27]. For the case of a Ga target, in Figure 7 we provide the coefficients associated with the spin dependent and coherent interactions given in (9) as functions of the WIMP mass by assuming three typical values of the detector threshold energy = 0, 5, 10 keV. For the special case of = 0, all plots peak at 35 GeV, while for higher threshold energies are shifted towards higher values of the WIMP mass. The calculations take also into account the annual modulation which is represented by the curve thickness. As can be seen from the figure, the modulation signal varies with respect to the WIMP mass, being larger for GeV, while its magnitude is slightly different for the spin dependent and coherent channels.

Proceeding further, in Figure 8 we evaluate the expected event rates for the four target nuclei assuming elastic WIMP scattering for WIMP candidates with mass GeV, by adopting the nucleonic-current parameters , , , and . As in the previous discussion, the thickness of the graph accounts for the annual modulation. We find that there is a strong dependence of the event rate on the studied nuclear isotope. Again the modulation is found to decrease for heavier mass. Among the four studied nuclei, we come out with a larger event rate for the case of a Ga nuclear detector, since , , and are all positive and have similar values. For Ge, is negative and its magnitude is comparable to and , while for As, is positive but small, and finally for I, and are relatively smaller and is large. The coherent contribution has more or less similar values for all nuclei considered.

For each component of the Solar, Atmospheric, and DSNB neutrino distributions we calculate the expected neutrino-floor due to CENS, by considering the target nuclei presented in Table 3. In our calculations, we neglect possible recoil events arising from Geoneutrinos as they are expected to be at least one order of magnitude less that the aforementioned neutrino sources (see, e.g., [25, 26]). In order to make a quantitative estimate of the neutrino-floor, here we do not consider neutrino oscillations and we assume that CENS is a flavour blind process in the SM. The differential event rate due to CENS, for the various dark matter detectors considered in the present study, is presented in Figure 9. It can be noticed that the neutrino background is dominated by Solar neutrinos at very low recoil energies. We stress that, for the typical keV-recoil thresholds of the current direct detection dark matter experiments only the and B sources constitute a possibly detectable background. From our results we conclude that, for recoil energies above about 10 keV, Atmospheric neutrinos dominate the neutrino background event rates, having a tiny contribution coming from the DSNB spectrum.

The number of expected background events due to CENS for each component of the Solar, Atmospheric, and DSNB neutrino fluxes is illustrated in Figure 10. Similar to the differential event case, at low energies the neutrino background is dominated by the Solar neutrino spectrum with the dominant components being the and B neutrino sources. The results imply that future multi-ton scale detectors with sub-keV sensitivities may be also sensitive to Be and neutrinos. We comment however that such sensitivities will be further limited due to the quenching effect of the nuclear recoil spectrum which is not taken into account here. Moreover, it is worth mentioning that neutrino-induced and WIMP-nucleus scattering processes provide similar recoil spectra; e.g., the recoil spectrum of B neutrinos may mimic that of a WIMP with mass 6 GeV (100 MeV) [28].

At this point, we consider additional interactions in the context of new physics beyond the SM that may enhance the CENS rate at a direct detection dark matter experiment. Specifically we study the impact of neutrino EM properties as well as the impact of new interactions due to a mediator, on the neutrino floor. In our calculations we assume the existence of a neutrino magnetic moment , extracted from CENS data in [41] as well as the corresponding limit from scattering data of the GEMMA experiment, e.g., [74]. Regarding the interaction we consider typical values such as , , and , [75]. Following [64], by assuming universal couplings, our calculations involve the product of neutrino and quark couplings defined as (for a comprehensive study involving the flavour dependence of the couplings the reader is refereed to [76])The corresponding results are presented in Figure 11, indicating that such new physics phenomena may constitute a crucial source of background even for multi-ton scale detectors with sub-keV capabilities. We stress, however, that the latter conclusion depends largely on the assumed parameters, which currently are unknown.

Before closing, we estimate the difference in the calculated number of neutrino-floor events between the conventional Helm-type and DSM predictions by defining the ratioFor each nuclear system, the corresponding results are presented in Figure 12 indicating that the differences can become significant, especially in the high energy tail of the detected recoil spectrum.

6. Conclusions

In this work, we studied comprehensively the expected event rates in WIMP-nucleus and neutrino-floor processes by performing reliable calculations for a set of prominent nuclear materials of direct dark matter detection experiments. The detailed calculations involve crucial nuclear physics inputs in the framework of the deformed shell model based on Hartree-Fock nuclear states. This way, the nuclear deformation and the spin structure effects of odd- isotopes that play significant role in searching for dark matter candidates are incorporated. The chosen nuclear detectors involve popular nuclear isotopes in dark matter investigations such as the Ga, Ge, As, and I isotopes. The DSM results indicate that Ga needs further investigation by employing another effective two-body interaction than the one used in the chosen set of nuclear isotopes.

The deformed shell model (DSM) employed for the nuclear structure calculations in this work is very well tested in many examples in the past [47] for nuclei with =60–90. Therefore, in our study we have chosen the dark matter candidates Ga, Ge, and As. In addition, to extend DSM to heavier nuclei of interest in dark matter detection, we have considered I and the results, reported in the present paper, are quite encouraging. In the near future we will consider Xe isotopes that are also of current interest. For lighter candidate nuclei, such as Na, Si, and Ar, clearly shell model will be better choice and DSM may also be tested for these isotopes.

More importantly, by exploiting the expected neutrino-floor due to Solar, Atmospheric, and DSNB neutrinos, which constitute an important source of background to dark matter searches, the impacts of new physics CENS contributions based on novel electromagnetic neutrino properties and mediator bosons have been estimated and discussed. Our results also indicate that the addressed novel contributions may lead to a distortion of the expected recoil spectrum that could limit the sensitivity of upcoming WIMP searches. Such aspects could also provide key information concerning existing anomalies in -meson decay at the LHCb experiment [77] and offer new insights into the LMA-Dark solution [78, 79].

Finally, the present results indicate that the addressed nuclear effects may become significant, leading to alterations especially in the high energy tail of the expected neutrino-floor as described by effective nuclear calculations, thus motivating further studies in the context of advanced nuclear physics methods such as the deformed shell model or the Quasiparticle Random Phase approximations and others. Such a comprehensive study using available data of the COHERENT experiment is under way and will be presented elsewhere.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

R. Sahu is thankful to SERB of Department of Science and Technology (Government of India) for financial support. D. K. Papoulias is grateful to Professor Naumov for stimulating discussions.