Exploring CEvNS with NUCLEUS at the Chooz Nuclear Power Plant

Coherent elastic neutrino-nucleus scattering (CE$\nu$NS) offers a unique way to study neutrino properties and to search for new physics beyond the Standard Model. Nuclear reactors are promising sources to explore this process at low energies since they deliver large fluxes of (anti-)neutrinos with typical energies of a few MeV. In this paper, a new-generation experiment to study CE$\nu$NS is described. The NUCLEUS experiment will use cryogenic detectors which feature an unprecedentedly low energy threshold and a time response fast enough to be operated in above-ground conditions. Both sensitivity to low-energy nuclear recoils and a high event rate tolerance are stringent requirements to measure CE$\nu$NS of reactor antineutrinos. A new experimental site, denoted the Very-Near-Site (VNS) at the Chooz nuclear power plant in France is described. The VNS is located between the two 4.25 GW$_{\mathrm{th}}$ reactor cores and matches the requirements of NUCLEUS. First results of on-site measurements of neutron and muon backgrounds, the expected dominant background contributions, are given. In this paper a preliminary experimental setup with dedicated active and passive background reduction techniques is presented. Furthermore, the feasibility to operate the NUCLEUS detectors in coincidence with an active muon-veto at shallow overburden is studied. The paper concludes with a sensitivity study pointing out the promising physics potential of NUCLEUS at the Chooz nuclear power plant.


Introduction
The existence of neutral-current neutrino interactions implies the existence of elastic neutrinonucleus scattering [1]. The cross-section of this process [2] is given in the Standard Model of Particle Physics (SM) by where G F is the Fermi constant, Q W is the nuclear weak charge, m N is the total mass of the nucleus, E R is the nuclear-recoil energy, and F(q 2 ) is the nuclear form factor as a function of the momentum transfer q. The weak charge is given by Q W = N − Z · (1 − 4 · sin 2 θ W ), with N and Z being the number of neutrons and protons of the target nucleus, and θ W the Weinberg angle. The maximum nuclear recoil energy E R is given by E max R = 2 E 2 ν /(m N + 2 E ν ), where E ν is the incident neutrino energy. For momentum transfers smaller than the inverse of the nuclear radius (typically for E ν 30 MeV) the form factor is close to unity and the scattering is coherent over all nucleons in a nucleus [1]: this process is called coherent elastic neutrino nucleus scattering (CEνNS). For larger momentum transfers the full coherence is no longer given and the form factor drops below unity.
The coherent cross-section is boosted by a factor of N 2 , and can exceed the cross-section of the standard neutrino detection methods, such as inverse beta-decay (IBD), by more than two orders of magnitude [5]. Contrary to the enhancement of the cross-section, the experimental signature, i.e. the nuclear recoil, is suppressed in energy by the nucleus mass, m N . Figure 1 shows the differential count rate of coherently scattered antineutrinos in different targets. For E R ≤ 100 eV, the CEνNS count rate in a CaWO 4 or Li 2 WO 4 crystal exceeds the rate on lighter nuclei such as present in Al 2 O 3 by more than an order of magnitude. However, as the CEνNS count rate decreases rapidly with increasing energy, an energy threshold well below 100 eV is necessary to explore the CEνNS signal with heavy targets.
While for IBD the kinetic energy of the antineutrino needs to be at least 1.8 MeV, CEνNS does not have an energy threshold. Thus, CEνNS provides a unique probe of the Standard Model (SM) at low energies: e.g. measuring the Weinberg angle, θ W , at low momentum transfer [6] or exploring fundamental neutrino properties such as the existence of a neutrino magnetic dipole moment [7]. Moreover, as a neutral current interaction, CEνNS is flavor insensitive, and, thus, a new probe to search for sterile neutrinos [8]. Any deviation from the SM prediction may reveal new physics beyond the SM, such as non-standard neutrino interactions, i.e. modified V-A quarkneutrino couplings, or new exotic neutral currents, see e.g. [2,9]. CEνNS offers a broad spectrum of possible applications in nuclear physics [10,11] and supernovae detection [6,12]. The enhancement of the CEνNS cross-section allows for a miniaturization of neutrino detectors, from the typical Very Near Site tonne-size to kilogram-or even gram-scale in the case of NUCLEUS, and thus, a possible practical application in nuclear reactor monitoring [13]. Furthermore, CEνNS of solar and atmospheric neutrinos will become an irreducible background for future dark matter experiments searching for weakly interacting massive particles, which thus profit from an independent measurement of the cross-section. The first observation of CEνNS was reported by the COHERENT collaboration at a 6.7σ confidence level [5]. The COHERENT detector [14], operating at the Spallation Neutron Source (SNS) in a neutrino flux of 4.3·10 7 ν/(s · cm 2 ), uses a 14.6 kg CsI[Na] scintillating crystal target. In contrast to reactor neutrinos, the neutrino beam produced at the SNS features a well defined energy spectrum reaching up to 50 MeV, i.e. partially above the coherence regime of elastic neutrinonucleus scattering.
This paper describes the NUCLEUS experiment at the Chooz Very-Near-Site (VNS), designed to study CEνNS using reactor antineutrinos. The VNS is presented in Section 2, together with on-site muon and neutron attenuation measurements. Section 3 focuses on the concept of the NUCLEUS experiment at the VNS, whereas the NUCLEUS target detectors are described in details in Reference [13]. Section 4 demonstrates that the NUCLEUS detectors can be operated at the VNS with a moderate dead time induced by the muon-veto. The potential to measure the CEνNS process with the NUCLEUS experiment at the VNS is presented in Section 5.

The Very-Near-Site at the Chooz Nuclear Power Plant
The VNS is a new experimental site between the two power reactors of the Chooz-B plant, planned to host the NUCLEUS experiment. The nuclear power plant, shown in Figure 2, is operated by the French company Electricité de France (EDF). The two N4-type pressurized water reactors, hereafter identified as B-1 and B-2, are separated by 160 m, and their respective cores are located about 7 m above the Chooz ground level. Each reactor runs at a nominal thermal power of 4.25 GW th and features a high duty cycle. The two cores are switched off for refuelling approximately one month per year during alternating periods.
Nuclear reactors are one of the strongest artificial continuous neutrino sources. Reactor-ν e are produced via the beta-decays of the fission products of 235 U, 238 U, 239 Pu and 241 Pu with energies up to 10 MeV. The meanν e energy is around 1.5 MeV [15]. Thus, reactor-ν e are expected to be in the fully coherent region for typical nuclear target elements, maximizing the CEνNS cross-section. Assuming an average number of sixν e per fission with an average energy release of 200 MeV, about 8 · 10 20ν e /s are produced by a 4.25 GW th power reactor [16]. Until 1998, the Chooz power plant hosted the long baseline (∼ 1000 m) neutrino oscillation experiment CHOOZ [17]. The main result was an upper limit on neutrino oscillation in theν e disappearance mode. Since 2008 the Chooz laboratory site is used by the neutrino oscillation experiment Double Chooz to measure the neutrino mixing angle θ 13 [18].
Thanks to these previous neutrino experiments, the Chooz reactor site has an existing infrastructure for scientists, including office, meeting and storage rooms, which can be made available for future projects. The acquired expertise and relationship developed in past collaborations with EDF greatly helps to establish the new site for a future neutrino experiment.

Description of the VNS
The VNS is a 24 m 2 room situated in the basement of a five-story office building at the Chooz power plant (see inset of Figure 2). The small size of the room requires that the full setup does not exceed a volume of several cubic meters and a weight of the order of O(10 t). To support loads on the floor of O(5 t/m 2 ), as expected for the NUCLEUS shielding, a weight support platform will be installed. Furthermore, minor modifications of the VNS are foreseen to meet the safety regulations.
Mechanical vibrations are a known cause of disturbance in cryogenic detectors. Vibrations can arise within the cryostat necessary to operate the detectors. In addition, vibrations may be introduced to the setup at the VNS by the turbines and generators of the nuclear power plant in the close vicinity. Further external sources of vibrations can be of seismic origin, as well as introduced by the elevator system, or movements in the office building. A vibration measurement campaign is planned for 2019. Based on the campaign results, a system to attenuate the vibrations will be designed.
The baseline of the VNS to the two reactor cores is 72 m to B-1, and 102 m to B-2 respectively. The expected neutrino flux at the VNS is about 10 12ν e /(s · cm 2 ). At this distance, no neutron background from the reactor cores is expected (see Section 3.1). However, for an expected overburden of less than 10 meters of water equivalent (m.w.e.), see section below, special care of the cosmic-ray induced background needs to be taken. Therefore, a campaign to characterize the neutron-and muon-background at the VNS was performed from October 2017 until May 2018. The campaign results are presented below and used to estimate the upper limit on the trigger rate of the NUCLEUS muon-veto in Section 4. Furthermore, they will be used to optimize the design of the compact shielding and to estimate the expected background count rate in the CEνNS target detectors in the future.

Measurement of the Muon Flux Attenuation
The muon rate at the location of the VNS was measured with the cosmic wheel developed by the Centre de physique des particules de Marseille (CPPM) for the "Scienceà l'école" outreach program [19]. The device consists of three (26 × 14) cm 2 plastic scintillator plates separated by 10 cm. Each plastic scintillator is read out by a photo-multiplier tube (PMT). Time-coincident signals in all three scintillators are interpreted as a cosmic muon event. Therefore, the number of background events from external gamma-rays or random coincidences are expected to be negligible. The mechanical structure allows to measure different zenith positions with an opening angle of 70 • . Figure 3 shows the measured muon rate above ground (open red markers) and four different azimuthal orientations at the VNS (filled markers). The angular zenithal distribution of the muon rate above ground follows the expected cos 2 θ law for cosmic muons at the surface [20], shown by the red line. Based on the attenuation a µ of the muon rate at the VNS, R VNS , with respect to the surface, R surface , the overburden m 0 is approximated with the following relation: where d = 1 + m 0 /10, and m 0 is given in m.w.e. [21]. Figure 4 shows a projection of the overburden over the VNS as obtained from the muonattenuation measurements. Since the VNS is not centered with respect to the building, the angular    Table 1). The averaged attenuation factor is given by a µ = 1/4 · i R VNS,i / R surface = 0.71 ± 0.01, where the error is given by the standard deviation. This attenuation corresponds to an overburden of m 0 = (2.9 ± 0.1) m.w.e.

Measurement of the Fast Neutron Flux Reduction
Neutrons are a potentially harmful background (discussed in Section 3.1). To evaluate the reduction of the neutron flux at the VNS, an array of neutron counters from Technische Universität München was deployed at Chooz. Each detector consists of a hexagonal cell with a diameter of 91 mm and a height of ∼ 50 mm. It is filled with approximately 300 g of liquid scintillator (EJ-301 from Eljen Technology [22]) which has a high discrimination capability between gamma-ray induced electronic and neutron-induced proton recoil. In order to collect the scintillation light the detectors are equipped with Philips XP-3461-B PMTs which are optically coupled to a perspex window. Cell and PMT are enclosed in a plastic tube. A schematic drawing of the detectors can be seen in Figure 5.
The setup is not only sensitive to neutrons but also to atmospheric protons. However, the flux of protons is reduced by a factor of ∼ 38 compared to the neutron flux whereby the protons are neglected in this measurement [23]. An array of seven such neutron detectors was operated at the Chooz power plant, at the surface and at the VNS location. The surface measurements were performed in a hut close to the Double Chooz near site (see Figure 2) and at negligible overburden. During the measurement campaign, the detectors were operated for a cumulative time of 15 h at the surface and 21 h at the VNS location.
A reduced data set of three detectors, selected for discrimination threshold and performance stability, was used for the final analysis. The detectors were calibrated by identifying the Compton edges of various gamma sources ( 22 Na, 137 Cs), wherefore the energy is given in units of electron equivalent, keV ee , in the following. The discrimination thresholds of the used detectors are 100 keV ee , 400 keV ee and 650 keV ee , respectively. Neutron-induced nuclear recoils in the liquid scintillator are unambiguously identified by a generic pulse-shape parameter based on the difference in the pulse decay time with respect to electron recoils. The cut is defined using an AmBe neutron calibration data set and then applied to the background data.
The data show a statistically compatible shape of the neutron recoil spectra at the VNS and surface. Thus, an energy-independent neutron reduction factor for the range 100 keV ee -2 MeV ee is found: at the VNS the neutron flux is reduced by a factor of 8.1 ± 0.4 compared to the surface (see Figure 6). A detailed study of the data is in preparation. Figure 6: Neutron-induced recoil rate in the liquid scintillator detectors up to an energy of 2 MeV ee observed at the VNS (green histogram, solid line) and a nearby surface location (blue histogram, dashed line). For the calculation the best three detectors in terms of discrimination threshold and stability are used. Each bin displays the exposure-weighted average of the rates from all detectors with a sufficiently low threshold, along with error bars following Poissonian statistics. The histogram uses logarithmic binning with small bins at low energy to display the shape of the spectrum and large bins at high energy to account for the low statistics.

The NUCLEUS Experiment at the VNS
A high signal-to-background ratio is one of the key requirements for the success of a CEνNS experiment. The primary CEνNS target material in NUCLEUS will be CaWO 4 . The enhanced cross-section for heavy nuclei such as tungsten, together with an envisioned detector threshold of O(10 eV), will give a strong CEνNS signal (see Figure 1 and Equation 1). One of the main experimental challenges is to achieve a sufficiently low background level in the CEνNS target detector. To reach the desired background count rate of the order of 10 2 counts/(keV·kg·day), the NUCLEUS experiment will apply passive and active background suppression methods.
In the following section, the potential background sources for the NUCLEUS experiment at the VNS are discussed. A first concept of the NUCLEUS experiment is presented in Section 3.2 with focus on active and passive background suppression. The BASKET R&D project, presented in Section 3.3, gives the opportunity to complement the CEνNS detection concept of NUCLEUS with a unique idea to characterize the neutron background, one of the most challenging backgrounds for CEνNS experiments.

Sources of Background
Experimental sites close to nuclear reactors are usually located at shallow depth, typically corresponding to an overburden of O(10 m. w. e.). The background sources of interest for the NUCLEUS experiment are here categorized into external, internal, cosmic-ray induced and reactor-correlated backgrounds. The dominant sources are briefly summarized below, for a more detailed description the reader is referred to [24,25].
The external background is dominated by gamma-rays from α-and β-decaying nuclides in the 232 Th, 238 U and 235 U decay chains, as well as 40 K which is present in material surrounding the experiment. Gamma-rays originating in the concrete of the laboratory building can be shielded by high-Z material, such as lead. Materials used for components close to the target detectors, such as mounting and support structures or electronic components, need to be carefully selected, e.g. by exploiting γ-ray spectroscopy, and undergo cleaning and purification processes [25].
Nuclides from the thorium and uranium decay chains may also be present in the detector target itself. In the case of NUCLEUS, this internal contamination can be reduced e.g. by special measures during crystal production [26]. In the uranium and thorium decay chains, radon is produced which escapes by diffusion or nuclear recoil accompanying the ejection of an α-particle. Radon is the strongest source of airborne radioactivity [25]. To prevent radon deposition on the detector surface, special cleaning procedures are needed during detector assembly and installation.
The above-mentioned decays typically produce γ-rays with energies up to 2.6 MeV, which is the highest naturally occurring γ-line from the decay of 208 Tl (a progeny of 232 Th). γ-rays of MeV energies predominantly interact via Compton scattering, thus the energy deposited in the target detector ranges from zero up to the energy of the Compton edge. Since γ-rays likely interact via multiple Compton scattering, a fiducialization of the NUCLEUS detector is a promising tool to actively suppress the external γ-background (see Section 3.2 for further details).
The energies of α-and β-particles typically reach up to several MeV. In β-decays, the available energy is split between the emitted ν and the e +/− . While a β-particle is typically absorbed within a few mm in a solid, the neutrino remains undetected. If the β-decay occurs on the surface of the detector or is followed by the emission of γ-rays, the fiducialization of the detector volume allows to discriminate this background from CEνNS events. The range of an α-particle does not exceed a few tens of µm in a solid. If released in the bulk of the target detector, α-particles are fully absorbed, and the event exhibits an energy well above the (sub-)keV signal of a CEνNS event. However, if the α-decay happens on or close by the detector surface a significant fraction of the energy remains undetected and such an energy deposit may mimic a (sub-)keV event. By instrumenting the surrounding material surface events can be discriminated from CEνNS events (see Section 3.2).
Most of the secondary cosmic rays from air showers like electrons, γ-rays, protons and pions are absorbed by the building structure or a dedicated lead shielding. The remaining muons and neutrons are the main components of cosmic-ray-induced background at shallow overburden [24]. Besides a source of background via ionization, pair production and Bremsstrahlung, muons also produce neutrons in nuclear reactions, especially in high-Z materials like lead. The flux of neutrons caused by air showers decreases rapidly with increasing overburden. Already with an overburden of 3 m.w.e. the neutron flux starts being dominated by neutrons produced by muons in a massive lead shielding [24]. The experimental signature of a neutron is a nuclear recoil, and thus, neutrons are a particularly dangerous source of background for any CEνNS experiments. Furthermore, they can be highly penetrating and difficult to shield. Neutrons may also produce radio-nuclides by inelastic scattering and neutron-capture in the detector or surrounding materials.
To reduce the muon-induced background, an active muon veto, typically based on plastic scintillators is used: all events in coincidence with a muon event are discarded. In this way, muon-induced events, which are e.g. induced by ionization or Bremsstrahlung, or by neutrons produced in the shielding material, are vetoed. Neutron-induced radio-isotopes with lifetimes much shorter than the coincidence window are rejected as well. In order to decrease the time for thermalization and capture of neutrons in the shielding, moderators, e.g. polyethylene (PE), are used. To attenuate neutrons produced inside the shielding, it is most efficient to place a neutron absorber inside the high-Z shielding material [24].
Nuclear reactors emit a large number of γ-rays and neutrons. γ-rays are produced in the fission process as well as the subsequent β-decay, neutron capture or de-excitation of the fission products. High energy γ-rays may produce neutrons in photo-nuclear reactions e.g. in the concrete of the reactor building [27]. The VNS is located at a distance of more than 70 m from the cores, with O(10 m) of rock and concrete in the line of sight. Furthermore, the reactor vessel is shielded by a thick layer of steel. Therefore, γ-ray induced neutron background is negligible.
Fission induced neutrons are emitted with a mean energy of 2 MeV. Most of these primary neutrons are thermalized in the reactor. The expected reactor-correlated fast neutron background at the VNS is estimated referring to the experience obtained with the NUCIFER [28] experiment, a project designed for online monitoring of nuclear reactors. The NUCIFER detector was based on Gd-loaded liquid scintillator and was operated 7.2 m from the core of the Osiris research reactor (70 MW) in Saclay, France. The estimated neutron elastic scattering rate on hydrogen in the whole NUCIFER target volume was 4 · 10 −5 events/(day) for energies above 2 MeV [28]. Scaling this result to the relevant parameters of NUCLEUS, i.e. reactor power of 2·4.25 GWth, distance of 72 m between the detector and the core, and assuming the same scattering cross-section on each nuclide of CaWO 4 as on hydrogen, indicates that the neutron background from the reactor core is negligible.  (1), where the target detectors are installed in its lower part, the experimental volume (2). The latter is surrounded by a passive shielding (3) consisting of alternating layers of borated PE and lead. The top part will be completed by a cold shield inside the cryostat. The outermost layer is the active muon-veto (4) made from plastic scintillator panels. The full setup with a footprint of approximately 1 m 2 is placed on a weight support platform (5). Further support structures, e.g. for the cryostat, are omitted for clarity.

The NUCLEUS Experiment
A promising technology for the study of CEνNS are cryogenic calorimeters. Cryogenic detectors measure the temperature rise ∆T = ∆E/C following an energy deposition ∆E in a target with a heat capacity C. For crystalline material kept at a temperature O(mK), C is small enough to achieve a measurable ∆T even for small energy depositions. The detector technology of NUCLEUS is based on CRESST cryogenic detectors [29] which have world-leading sensitivity in the field of direct low mass dark matter search and are driven by similar requirements as the study of CEνNS: a sub-keV energy threshold as well as a low intrinsic background. The CRESST experiment has developed CaWO 4 and Al 2 O 3 crystals as targets with masses between 24 g and 300 g.
The NUCLEUS concept expands this technology to an ultra-low-threshold (O(10 eV)) gramscale cryogenic calorimeter [30] as a target detector for CEνNS combined with cryogenic veto detectors for active background discrimination. The experiment will proceed in two phases: in a first step, NUCLEUS-10g, a detector with a 10 g CEνNS target will be deployed. A later stage, NUCLEUS-1kg, foresees an upgrade to a total target mass of 1 kg. Figure 7 shows a first design of the proposed CEνNS experiment at the VNS. The NUCLEUS target detectors are installed in a cryostat. The experimental volume is surrounded by a shielding for active and passive background reduction.

The NUCLEUS Target Detectors
The NUCLEUS detector concept is a fiducial-volume cryogenic detector (see Figure 8) based on three individual calorimeter systems: (1) an array of gram-scale cryogenic calorimeters with ultralow energy threshold as target detector, (2) a low-threshold inner veto holding and encapsulating the target, and (3) a surrounding kg-scale cryogenic detector used as outer veto.
CaWO 4 , Al 2 O 3 , Si and Ge are well known materials for cryogenic detectors and suitable candidates for the target volume of the experiment. Detectors with a total mass of O(1 g) enable ultra-low thresholds [30] and imply a total event rate per crystal of 1 Hz which is crucial to keep the deadtime at a value of O(1%). In the CRESST experiment [29], transition edge sensors (TESs) have proven to be a highly sensitive thermometer able to measure the low recoil energies produced in CEνNS. With a 0.5 g prototype detector made from a (5 mm) 3 Al 2 O 3 cubic crystal, an unprecedented ultra-low threshold of E th = (19.7 ± 0.9) eV has been reached [13], one order of magnitude lower than previous devices. Such gram-scale crystals exploit the full potential of the CEνNS signal (see Figure 1). Arrays of these gram-scale cryogenic calorimeters allow deploying a larger total target mass. Figure 8 shows a schematic view of NUCLEUS-10g: a 3×3 array of CaWO 4 (6 g) and a 3×3 array of Al 2 O 3 (4 g). This multi-target approach will provide an in-situ background characterization: while the CEνNS scattering rate (∼N 2 ) is strongly enhanced for CaWO 4 , fast neutrons are expected to induce comparable signatures in CaWO 4 and Al 2 O 3 due to scattering on O nuclei [13]. Furthermore, coincident background events in several target detectors can efficiently be discriminated from CEνNS events, which are single scatters.
The target detectors will be embedded within cryogenic veto detectors which allow for an efficient active background reduction. The NUCLEUS detector support structures are made of Si-wafers read out by TESs which enable a 4π low-threshold inner veto. Background events originating from surface α and β contamination of the detector components can be rejected by anticoincidence [13]. Both target and inner veto will be enclosed in a kg-scale outer veto which is operated in anti-coincidence mode to reject neutron and γ-ray backgrounds. Ge and CaWO 4 are established materials for cryogenic detectors with masses up to 1 kg and typically reach energy thresholds of a few keV. Ge single crystals are commercially available in extremely high purity and large sizes. CaWO 4 is especially interesting due to the presence of tungsten with a high cross section for γ-radiation, and oxygen which enables efficient neutron moderation. An outer veto from CaWO 4 would require R&D into large-diameter crystal production. Studies based on Monte Carlo simulations show that with the outer and inner vetos a background suppression of O(10 3 ) can be achieved [13].
The NUCLEUS detector concept foresees to scale the total target mass from gram to kilogram. The production of the cryogenic detectors is based on well-established techniques of the semiconductor industry. In particular, the capability to produce multiple target elements simultaneously before cutting and polishing has been demonstrated, which is crucial for the later stages of NUCLEUS.  (1), an inner veto (2) and an outer veto (3) with a diameter of 10 cm. The assembly is held mechanically by a non-instrumented support structure (4). The target is operated in anti-coincidence with the inner and the outer veto. See text for details.

The Cryostat
The base temperature of O(10 mK) required to operate cryogenic detectors is typically provided by "wet" dilution refrigerators relying on a continuously evaporating bath of liquid helium at 4 K. At the VNS, the handling of cryogenic liquids can be avoided using a "dry" dilution refrigerator which replaces the helium bath with a pulse-tube cryocooler integrated into the cryostat. Though much simpler in operation and infrastructure requirements, this feature makes a "dry cryostat" a challenging environment for cryogenic detectors due to vibrations induced by the pulse-tube, which will degrade the detector performance if unaddressed. This challenge can be met by dedicated work on vibration decoupling [31], which is mandatory for the NUCLEUS setup. The CEνNS target detectors are installed in the experimental volume inside the cryostat. The cryogenic infrastructure is chosen such as to host both NUCLEUS phases. The experimental cryogenic volume is surrounded by a compact passive and active shielding.

Active and Passive Shielding
Many examples in literature demonstrate that a low background count rate can be reached with a shallow overburden of O(10 m.w.e.). Two promising examples are GIOVE, a highly sensitive HPGe spectrometer in the shallow depth laboratory (15 m.w.e.) at the Max-Planck-Institut für Kernphysik [32], and the Dortmund low-background facility with an overburden of 10 m.w.e [23]. Both set-ups feature an active muon-veto and passive shields of lead, copper, and borated PE. With these graded shielding structures background rates of 0.13 counts/(keV·kg·d) for GIOVE and 1.4 counts/(keV·kg·d) for the Dortmund low-background facility are reached in the 40-2700 keV energy range.
The goal of the NUCLEUS experiment is to achieve a background count rate of 10 2 counts/(keV· kg · d) in the sub-keV region. At present, this energy range has never been probed in a shallow depth laboratory. Based on the experience of the aforementioned experiments operating in the keV-region, the passive shielding of NUCLEUS will consist of several alternating layers of borated PE and lead. To avoid a hole directly above the detectors, a part of the shielding has to be mounted inside the cryostat and kept cold during operation. Neutrons produced in nuclear reactions are moderated by the PE and captured on boron. External gamma-rays and such produced in radiative neutron captures (for boron the most intense line is at 478 keV [33]) are attenuated by the lead. As copper can be produced with very low intrinsic radioactive contamination compared to lead (see e.g. [34]), copper may be used as innermost shielding layer. The cubic passive shielding will have an approximated edge length of 1 m. Its final dimension and layout, i.e. the sequence and thickness of the PE and Pb layers, as well as the boron doping of PE, will be based on dedicated Monte Carlo studies and is the subject of future publications.
As discussed in Section 3.1, cosmic-ray induced events are an unavoidable source of background. Cosmic muons cannot be sufficiently attenuated by a passive shielding, thus, muon-induced background events need to be identified and removed from the data. An active muon veto consisting of 5 cm thick plastic scintillator plates is planned. With an energy loss of ∼ 2 MeV/cm a 4 GeV muon (mean energy of the cosmic muons at the surface) deposits about 10 MeV in a veto panel [20]. The most energetic gamma-ray from natural radioactivity is at 2.6 MeV, see Section 3.1. Thus, cosmic muons can efficiently be discriminated by a simple energy cut. An event in the muon-veto that exceeds the energy threshold, causes a muon-trigger. Events in the target detector in coincidence with a muon trigger are disregarded.

The BASKET Detectors
BASKET (Bolometers At Sub-KeV Energy Threshold) is an R&D program started in 2017 at CEA in collaboration with CSNSM (Centre de Sciences Nucléaires et Sciences de la Matière) to develop innovative detectors for neutrinoless double beta decay and CEνNS [35]. For the latter, the project focuses on the development of Li 2 WO 4 crystals as a new absorber material together with new thermal sensors to optimize the time response and energy resolution. First tests on an 11-g crystal, which was read out with a Neutron Transmutation Doped Ge sensor and a Ge Neganov-Luke light detector, showed that Li 2 WO 4 exhibits good bolometric and scintillation properties. Furthermore, it could be demonstrated that α-particles can be discriminated from β-particles and γ-rays with a discrimination power of more than five sigma [36]. Further tests with 1-g crystals coupled to metallic magnetic calorimeters [37] and NbSi TES [38] are currently performed. The goal is to reach an energy threshold of O(10 eV) and a time response of O(100 µs). A possible synergy of BASKET with the NUCLEUS experiment is the deployment of gram-size Li 2 WO 4 crystals as target detectors. The CEνNS rate is similar to the one expected for CaWO 4 . Moreover, an additional target compound to CaWO 4 and Al 2 O 3 may yield supplementary information on the background yet to be explored at a few tens of eV.
Furthermore, Li 2 WO 4 is a very interesting material to monitor neutron backgrounds and could be used for the outer veto of NUCLEUS. The heavy element tungsten has a high attenuation power for γ-rays. Neutrons can be tagged using the neutron capture on 6 Li which has a natural abundance of about 7.4%, allowing for an in-situ characterization of the neutron background [39]. The reaction 6 Li + n → 3 H + 4 He, with an energy release of E = 4.78 MeV + E n , where E n is the kinetic energy of the neutron, produces only heavy charged particles. The latter are well separated from electronic recoils (β, γ, µ) by the reduced light signal. Thus, the identified 6 Li neutron capture Table 2: Number of simulated muon triggers, N µ , for the NUCIFER experiment for different muon energies E. In total 10 4 muons were simulated for each energy. The weights are given within the logarithmic energy bin centered at the energy E with respect to the total flux. The total number of triggers is given by the weighted sum of the mono-energetic simulations.

Muon Induced Dead Time Considerations
To reach the required background level for NUCLEUS at a shallow experimental site, a compact passive shielding in combination with an efficient muon veto is compulsory (see Section 3.2.3). At the Earth's surface the count rate of cosmic muons is ∼100 Hz/m 2 [20]. Hence, to avoid significant detector dead time the timing of the coincident events has to be sufficiently fast. While typical muon-veto panels feature pulse rise times < 1 µs, pulses of cryogenic detectors as the one used for NUCLEUS are typically orders of magnitude slower which makes the detector operation at shallow depth challenging.

Monte Carlo Simulations of Cosmic Muons at the VNS
A Monte Carlo (MC) simulation tool based on GEANT4 [40] was used to simulate cosmic muons at the VNS and to estimate an upper limit on the trigger rate of the muon-veto. Muons are randomly generated on a plane tangent to a half-sphere centered at the detector whereas the distribution of the point of tangency follows the cos 2 θ-distribution of muons on the surface. The direction of the generated muons is parallel to the normal vector of the plane, pointing towards the detector.
To validate the MC simulations, a shielding configuration similar to the one used in the NU-CIFER experiment [28] was simulated and the amount of simulated hits in the shield was compared to the published experimental data. The differential flux of incident atmospheric muons dΦ/dEdΩ is approximated according to an adapted Gaisser parametrization [41,42]. For simplicity, only mono-energetic muons between 4.5 MeV and 450 GeV in logarithmic binning were simulated. The results were weighted according to the flux within the given energy range with respect to the total flux. The trigger rate of the muon-veto is calculated by where N MC is the number of generated muons in the simulation, N µ is the number of energy depositions in the muon-veto (i.e. a muon trigger) which pass the 10 MeV threshold [28] of the NUCIFER experiment, Φ µ = 70/(sr · s · m 2 ) is the total muon flux at the surface [20], A = (3.5 m) 2 is the surface of the plane where the muons are generated, and a µ is the muon-attenuation at the experimental site. In total 10 4 muons were generated for each muon energy bin. Table 2 summarizes the simulated muon trigger rates in the NUCIFER setup. The muon simulations with the simplified NUCIFER geometry yield a total muon trigger rate of (535 ± 9) Hz, by taking into account a veto-efficiency of 97% the trigger rate is reduced to (519 ± 9) Hz. The MC result is about 50% higher than the measured muon trigger rate of 350 Hz [28]. This deviation may originate in the simplifications of the detector geometry, and in further unaccounted inefficiencies. It is believed that the applied approximations to the muon energy spectrum give the largest systematic uncertainty. The simulation assumes a uniform overburden of 12 m.w.e, and does not account for any changes in the muon energy spectrum with respect to the surface. Furthermore, the rough binning may be over-simplified.
As the goal is to estimate an upper limit on the trigger rate for the actual NUCLEUS muonveto, a verified overestimation of 50% in the simulated trigger rate is not an obstacle. To estimate the muon rate, a simplified cubic geometry with six plates, i.e. at the top, bottom and four side faces, of equal size was implemented in the simulation tool. Each plate of the muon-veto is made from standard PE (C 2 H 4 ) with a density of 0.96 g/cm 3 . For simplicity the generated muons have an energy of 4 GeV which corresponds to the mean energy of cosmic muons at the surface. For the muon attenuation, a uniform factor of a µ = 0.71 is assumed, see Section 2.2. No trigger threshold is applied, which is legitimate as only an upper limit on the trigger rate is estimated. The muon rate is calculated using Equation 3. As seen in Figure 9, the rate is proportional to the surface of the muon-veto. For a shielding covering a volume of 1 m 3 , as envisioned by NUCLEUS, the simulation yields an upper limit on the muon trigger rate of (487 ± 38) Hz.

Estimation of Muon-Induced Dead Time
The overall detector dead time is governed by the timing of the cryogenic pulses, since the timing of the muon-veto panels is orders of magnitude faster and, therefore, can be neglected. In the following, the pulse timing, i.e. the uncertainty σ τ of the determination of the pulse's onset is assumed to be proportional to the rise time. In cryogenic detectors, the rise time is given by the thermalization of the signal phonons at the surface of the absorber crystal, and is thus expected to scale roughly with the linear dimension of the target. The pulse rise time of the 0.5 g NUCLEUS detector prototype is τ r ≈ 300 µs [30], significantly faster than CRESST-II detectors with a mass of 300 g and a typical rise time of several ms [43], and with other cryogenic detector technology.
For the estimation of the detector dead-time, a time window of ±5σ τ is defined around every event in the muon-veto and removed from the total exposure time. For example, in case of a cubic muon-veto of 1 m side length with a count rate of R = (487 ± 38) Hz the pulse's onset has to be determined with σ τ 20 µs to restrict the dead time to 10%. The precision of the pulse onset of a fast cryogenic detector with a rise time of 100 µs has been measured using a pulsed neutron beam from an accelerator in [44]. The uncertainty of the onset determination is σ τ = (4.8 ± 0.4) µs corresponding to a time window of 48 µs around every event in the muon-veto to be removed from the exposure time.
Using the results of the dedicated MC simulation (see section 4.1) for the muon count rate of different shielding dimensions, the resulting dead time can be calculated with respect to the side length of the cubic-shaped muon-veto. Figure 10 shows the resulting dead time for a number of detector technologies with different pulse rise times. For a cryogenic bolometer of CRESST-II with a mass of 300 g, its slow signal implies a dead time above 10% even for small setups (red line in Figure 10). For the existing NUCLEUS prototype detector, the dead time is predicted to be about 7% for a cubic muon-veto of 1 m side length.
The rise time in NUCLEUS detectors can be influenced by changing the phonon collection area of the TES. A larger area leads to a quicker collection of signal phonons and thus a faster pulse rise time. For the NUCLEUS experiment, a rise time of 100 µs is envisioned. Under these conditions, the dead time, e.g. for a cubic-shaped muon-veto of 1 m side length, is about 2% and stays below 10% even for a dimension of (2 m) 3 .

Sensitivity of NUCLEUS to CEνNS at the VNS
The NUCLEUS experiment will proceed in a staged approach, as presented in Section 3.2: a first phase called NUCLEUS-10g, and an upgrade called NUCLEUS-1kg. At the VNS, with distances of 72 m and 102 m to the two 4.25 GW th cores of the CHOOZ nuclear power plant, NUCLEUS-10g will observe 0.33 (0.03) CEνNS-induced nuclear recoil events per day above 10 eV in the CaWO 4 (Al 2 O 3 ) array. The high CEνNS rate in the CaWO 4 array, induced by the large tungsten nuclei, allows for a fast discovery of CEνNS at the VNS. The strongly suppressed CEνNS rate in the Al 2 O 3 array provides an effective "neutrino-off" reference spectrum, which is of great advantage in case of an unknown background shape, as discussed in Section 5.3.

Sensitivity of NUCLEUS-10g for Different Flat Backgrounds
The results of a sensitivity study for NUCLEUS-10g, similar to the sensitivity study of Reference [13], are shown in Figure 11 for different background indices. For each background index and exposure time, 600 random spectra with a flat background and the Standard Model CEνNS expectation at the VNS are generated for the two target materials CaWO 4 and Al 2 O 3 in the energy range from 10 eV to 2 keV. To each pair of spectra, a flat-background plus CEνNS signal strength model ("free model") and a flat-background-only model ("null model") are fit. The "null model" has only one free parameter, the flat background rate, which is common to both target materials. The "free model" has one additional parameter, the signal normalization, for which a value of 1 corresponds to the Standard Model CEνNS expectation in each target material. The ratio of the maximum likelihood for these nested models (L free and L null respectively) is converted to a significance of CEνNS observation using Z = 2 · log (L free /L null ). The median, 5 th and 95 th percentile of the resulting Z distribution are used to draw the color bands (corresponding to a 90% probability interval).
For the most pessimistic background index of 1000 counts/(keV·kg·day), the signal rate in CaWO 4 exceeds the background rate only close to the assumed energy threshold of 10 eV, so that long exposures are necessary to distinguish the contributions by spectral shape. Below around 10 counts/(keV·kg·day), the experiment is background-free over most of the region of interest, so that the necessary live time becomes limited by signal statistics.
For a benchmark background rate of 100 counts/(keV·kg·day), a 5σ observation of CEνNS at VNS can be achieved in less than 40 days of measurement time.

Achievable Precision of NUCLEUS-10g and -1kg
As an absolute-rate experiment with the goal of studying the CEνNS cross-section, NUCLEUS has to rely on reactor neutrino flux predictions which will contribute a significant systematic uncertainty to the measurement. At energies above the IBD threshold, the reactor antineutrino spectrum is modeled with an uncertainty of 2-3% [45,46]. As nuclear recoils on tungsten induced by neutrinos below the IBD threshold extend up to 38 eV in nuclear recoil energy, NUCLEUS has the potential to become sensitive to this as yet unobserved low-energy part of the spectrum. An effort to predict the fluxes of reactor antineutrinos below the IBD threshold is thus necessary to interpret the NUCLEUS data at these lowest recoil energies.
Assuming a flat background of 100 counts/(keV·kg·day) and an energy threshold of 10 eV, the statistical precision achievable with NUCLEUS-10g after one year is 11%. Therefore, NUCLEUS-10g is expected to be limited by statistics over its complete measurement time. On the contrary, for NUCLEUS-1kg (with an assumed improved background index of 1 counts/(keV·kg·day) and an energy threshold of 4 eV) the achievable statistical precision with one year of data is 0.64%, so that this stage of the experiment will be limited by systematics early on in the measurement campaign. Figure 12 shows the precision (one standard deviation) on the CEνNS cross-section achievable by NUCLEUS, compared to the best current precision achieved by COHERENT [5] (black dashed line). For NUCLEUS-10g (blue), the solid line shows statistical precision only, while the dashed line adds a 10% systematical uncertainty in quadrature. NUCLEUS-10g can approach the expected systematic limit within few years of live time. The continuous green line in Figure 12 shows the statistical precision achievable with a 1 kg Ge target at the VNS. Adding an optimistic 1% systematical uncertainty (dashed green line) shows that NUCLEUS-1kg can reach a percent-level measurement of the CEνNS cross-section within few years of measurement time. The final precision on the CEνNS cross-section that can be reached by NUCLEUS-1kg depends strongly on the control of systematics, such as neutrino flux normalization, precise knowledge of the energy scale and modeling of backgrounds.

NUCLEUS-10g in the Case of a Non-flat Background
Since NUCLEUS is aiming to measure a signal in a previously unexplored energy range in a radiogenically challenging environment, assuming a flat background down to threshold may well be overly optimistic. Accepting a background model with a free shape parameter drastically reduces the detection significance in a single target material, as the signal can then also be fit reasonably well with the background-only model.
The multi-target approach of NUCLEUS-10g (i.e. deploying CaWO 4 as well as Al 2 O 3 target arrays) has the potential to mitigate the impact of a background of unknown shape. With a CEνNS signal more than an order of magnitude below the one in CaWO 4 (see Fig. 1), Al 2 O 3 can be conservatively regarded as measuring background under identical experimental conditions to the CaWO 4 array. Precise knowledge about the rate and shape of the spectra associated with the various background sources in the two materials is needed to apply the knowledge from this "sideband region" (Al 2 O 3 spectrum) to constraining the fit in the signal region (CaWO 4 spectrum).  This requires careful simulations and calibration campaigns using both target materials in different environments.
Assuming for simplicity an identical appearance of backgrounds in Al 2 O 3 and CaWO 4 allows to study the power of this multi-target approach in NUCLEUS-10g. To demonstrate the advantage of combining materials with different expected CEνNS signal strength, a scenario is considered here in which the background is non-flat and, more importantly, of unknown shape.
As a simple non-flat background model an exponential + flat model with the functional form R(E) = C + A · e −E/B is considered, where C is the flat background rate, A is the amplitude and B is the slope of the exponential. Changing the slope of the exponential background while keeping the total integral fixed and the flat background level constant at 100 counts/(keV·kg·day), shows the strongest impact on the sensitivity of NUCLEUS-10g after one year to occur for a slope of around 45 eV. Therefore, to be conservative, a background of this slope is chosen for the non-flat scenario. An amplitude of 1500 counts/(keV·kg·day) is set for the exponential background (see Figure 13) so that the background level is higher than the signal rate at all observed energies. In this scenario, there are 120 (636) CEνNS (background) interactions in the 6 g CaWO 4 array expected between 10 eV and 2 keV within one year. In the 4 g Al 2 O 3 array, these numbers change to 12 (415) CEνNS (background) interactions.
The CaWO 4 spectrum alone is fit nearly equally well by both the (three parameter) backgroundonly model and the (four parameter) model including the signal. Thus no sensitivity to the CEνNS cross-section can be derived from CaWO 4 alone. The simultaneous fit of both spectra uses the same background model in both materials (three parameters) and the SM CEνNS expectation for each material scaled by a common strength parameter. In this way, the Al 2 O 3 array helps constraining the background rate and shape, so that an expected preference for the background + CEνNS model of more than 4 σ can be extracted after one year.
This result shows the crucial advantage of the multi-target approach: the sensitivity is significantly improved even without any prior knowledge of the background shape (apart from the flat + exponential -parameterization). Even in case of a background with a signal-like shape and strength, CEνNS becomes observable through its characteristic dependence on the nuclear composition of the target. This advantage is complementary to the traditional concept of particle identification, i.e. distinction of electron recoil and nuclear recoil events. While the multi-target approach is statistical only, it can be used to constrain neutron background which features an identical experimental signature as CEvNS, i.e. a nuclear recoil, but a distinct dependence on the nuclear composition of different materials.

Conclusion
The VNS at the Chooz nuclear plant is a promising new experimental site for the planned CEνNS experiment NUCLEUS. The site is located in close distance to the two 4.25 GWth reactor cores of the Chooz nuclear power plant in France, providing an antineutrino flux of ∼ 10 12 ν/(s · cm 2 ). First muon attenuation measurements on-site indicate a shallow overburden of about 3 m.w.e. at the VNS.
The NUCLEUS detector concept provides a suitable technology for a new-generation CEνNS experiment at the VNS. Thanks to its unprecedentedly low energy threshold of ≤ 20 eV, a strong CEνNS signal is expected which allows a significant miniaturization of the detector size. Dedicated MC simulations show that the muon-induced dead time is expected to stay well below 10%, given the fast rise-time of the NUCLEUS detectors. The small size of the target detectors and the envisioned compact setup, which consists of passive and active shielding material, greatly reduce the overall size of the experiment matching the space constraints of the VNS. Using active and passive background reduction techniques, the NUCLEUS experiment aims to reach a background count rate of ≤ 100 counts/(keV·kg·day).
In the first phase of the experiment, NUCLEUS-10g aims at a first observation of CEνNS at a nuclear reactor with a total target mass of 10 g. Due to the demonstrated threshold in the 10 eV regime, NUCLEUS will probe the reactor antineutrino spectrum for the first time below 1.8 MeV -below the threshold of inverse beta-decay. The reach of NUCLEUS strongly depends on the achieved background level at low energy. NUCLEUS-10g will explore this background for the first time at energies below 100 eV and at shallow overburden. A statistical analysis shows that CEνNS can be observed within two weeks of live time assuming our benchmark background model. Even in a scenario in which shape and rate of the background are similar to that of the CEνNS signal, NUCLEUS-10g reaches a sensitivity of more than 4σ for CEνNS after one year of measuring time. NUCLEUS-10g, which is expected to take data from 2021 on, has the potential to achieve a final precision of ∼10% on the CEνNS cross-section.
The second phase, NUCLEUS-1kg, is planned for commissioning in 2023 and aims for precision measurements of the CEνNS cross-section at a level of 1%. Harnessing the full potential of NUCLEUS-1kg requires a significant reduction of systematic uncertainties. Particularly, to take advantage of the achievable statistical precision, an improved prediction of the reactor antineutrino spectrum is needed. This opens the door for the study of physics beyond the Standard Model of Particle Physics such as the search for non-standard neutrino interactions, a neutrino magnetic dipole moment or exotic neutral currents, as well as a test of the Reactor Antineutrino Anomaly [47].
NUCLEUS will demonstrate a new detector technology with ultra-low thresholds and allow to study fundamental properties of the neutrino, pushing the low-energy frontier in neutrino observations.