The \nu-cleus experiment: A gram-scale fiducial-volume cryogenic detector for the first detection of coherent neutrino-nucleus scattering

The -cleus experiment: A gram-scale fiducial-volume cryogenic detector for the first detection of coherent neutrino-nucleus scattering

Abstract

We discuss a small-scale experiment, called -cleus, for the first detection of coherent neutrino-nucleus scattering by probing nuclear-recoil energies down to the 10 eV-regime. The detector consists of low-threshold CaWO and AlO calorimeter arrays with a total mass of about 10 g and several cryogenic veto detectors operated at millikelvin temperatures. Realizing a fiducial volume and a multi-element target, the detector enables active discrimination of , neutron and surface backgrounds. A first prototype AlO device, operated above ground in a setup without shielding, has achieved an energy threshold of  eV and further improvements are in reach. A sensitivity study for the detection of coherent neutrino scattering at nuclear power plants shows a unique discovery potential (5) within a measuring time of weeks. Furthermore, a site at a thermal research reactor and the use of a radioactive neutrino source are investigated. With this technology, real-time monitoring of nuclear power plants is feasible.

pacs:
5.55.VjNeutrino, muon, pion, and other elementary particle detectors; cosmic ray detectors

1 Introduction

The detection of coherent neutrino-nucleus scattering (CNNS) is among the most challenging tasks of modern particle and astroparticle physics. A first observation of CNNS would be an important confirmation of the Standard Model of Particles and would open the door to new physics beyond the Standard Model (BSM).

Coherent neutrino-nucleus scattering (CNNS), first proposed in 1974 (1), is an unobserved neutral-current interaction predicted by the Standard Model of Particle Physics. Neutrino-nucleus scattering via Z-exchange becomes coherent over the nuclei at low transferred momenta, for large nuclei simultaneously boosting the interaction cross-section and reducing the recoil energies. The total elastic cross section for the process can be written as (2):

(1)

where is Fermi’s coupling constant, the Weinberg angle, , and are the nucleus’ proton number, neutron number, and total mass respectively, is the neutrino energy, and the resulting nuclear-recoil energy. The nuclear form factor describes the loss of coherence as a function of transferred momentum wavenumber . It can be understood as the Fourier transform of the nuclear weak charge density, and is close to unity for small (typically at  MeV).

The process remains unobserved until now due to the small recoil energies expected which challenge detector technologies. Multiple experimental efforts for detecting CNNS are taken globally: the COHERENT experiment (3) which is currently taking data at the Spallation Neutron Source (SNS) uses a combination of conventional CsI, Ge and liquid-Xe detectors. Various other experiments are planned or being commissioned such as CONNIE (4) using CCDs (5), TEXONO (6) using ionization-based Ge detectors, and MINER (7) and RICOCHET (8) using cryogenic detector technology.

The ultra-low-threshold cryogenic calorimeters presented here (and in (9)) put a rapid detection of this process within reach technologically. The relatively large cross section compared e.g. to neutrino-electron scattering makes this experimental approach interesting in two ways: 1) CNNS is detectable with a small-scale experiment and a total target mass of 10 g within a measuring time of several weeks (see below), far less costly than traditional neutrino facilities. 2) A managable scaling of the total target mass to the still moderate range of 1-10 kg opens up a new window for precision tests of neutrino properties and interactions beyond the standard model. A recent summary of CNNS sensitivity to BSM neutrino physics is given in (10), including the following potential observations.

  • Interpreted within the standard model, a precise measurement of the CNNS cross-section allows to determine the Weinberg angle at low energy scale through Equ. 1. Transferred momenta in CNNS are on the order of few , extending the reach of other planned low momentum-transfer precision experiments (11). Together with knowledge on electroweak precision observables (e.g. from LEP), this allows to probe the running of the weak mixing angle (10) which is precisely predicted in the standard model (12). This collective measurement has sensitivity to BSM contributions well above the LHC scale.

  • The neutrino-quark sector of neutrino Non-Standard Interactions (13); (14), i.e. modified V-A quark-neutrino couplings, may measurably modify the CNNS cross-section (10).

  • Exotic Neutral Currents (10), i.e. general (pseudo-)scalar, (axial-)vector or tensor couplings can induce modifications in the CNNS cross section and energy spectrum.

  • The possibility of observing active-to-sterile neutrino oscillations using CNNS is discussed in (15).

  • For very low energy thresholds, the magnetic moment of the neutrino (causing enhanced low-energy scattering with spin exchange) can be probed beyond current limits from neutrino-electron scattering (16).

2 The detector

2.1 A fiducial-volume cryogenic detector

A detector, sensitive to CNNS, faces two main challenges: an extremely low energy threshold combined with extraordinarily small background levels. We present a new gram-scale cryogenic detector which combines the possibility of lowest nuclear-recoil thresholds (( eV)) and the advantages of a fiducial volume device. Those provide active shielding by the outermost regions against external radiation which reduces the background level in the innermost target volume (the fiducial volume). Since an exact spatial position reconstruction of events is difficult to realize in thermal detectors, so far this potential could not be exploited.

Here, a cryogenic detector is presented which realizes a fiducial volume by combining 3 individual calorimeters: 1) a target crystal (the fiducial volume) with an extremely low threshold of ( eV), 2) an inner veto as a 4 veto against surface beta and alpha decays, and 3) an massive outer veto against external gamma/neutron radiation (see Fig. 1). Additionally, the inner veto acts as an instrumented holder for the target crystal allowing to discriminate holder-related events (e.g. from stress relaxations).

Figure 1: Schematic view of the new detector which consists of 3 individual cryogenic calorimeters. The combination of, both, the outer veto against external gamma/neutron radiation, and the inner veto against surface alpha and beta decays, significantly reduces the background level in the target crystal. In this way, a fiducial-volume cryogenic detector is realized. The inner veto acts additionally as instrumented holder of the target crystal to reject possible stress-related relaxation events.

2.2 Performance model for calorimeters

In order to design the new detector, a simple model was developed to predict the practically achievable performance of calorimeters of different geometry, material and mass (17). The model is based on experimental results of cryogenic CRESST-type detectors. The main results are derived here, insofar as they drive design-choices for the fiducial-volume cryogenic detector.

The fundamental equation describing a calorimeter is that, for a system in internal thermal equilibrium, the temperature rise

(2)

where is an energy deposit and is the heat capacity of the object. Reducing yields a large increase in temperature and thus a high sensitivity to small energies.

Present cryogenic detectors of  g achieve energy thresholds of  eV(18). In this work we investigate the performance and potential of gram-scale devices.

The fundamental energy resolution of cryogenic calorimeters is given by irreducible thermal fluctuations between the absorber and the thermal bath (19):

(3)

with the absorber’s temperature , heat capacity and the Boltzmann constant . This corresponds to theoretical energy resolutions of (1 eV) at  mK for massive calorimeters with masses of  g (20). Phonon processes in cryogenic calorimeters with thin-film transition-edge-sensors (TES) as considered in this work are well described by a dedicated thermal model (21).

Equ.  2 is valid for a thermometer measuring the temperature of an absorber. In practice, the thermometers of cryogenic detectors can only measure their own temperature. Equ. 2 thus changes to

(4)

where denotes the energy absorbed in the thermometer and the heat capacity of the thermometer film. In cryogenic calorimeters at very low temperatures ( mK), the energy deposition in the thermometer film happens via the absorption of non-thermal phonons, which propagate ballistically and interact directly with the metallic film electrons. Thus they are not affected by the weak thermal coupling between thermometer phonon and electron systems at such temperatures. To achieve sufficiently low heat capacities, temperatures as low as 10mK are required for these devices. The strong electron-phonon decoupling in the thermometer film at these temperatures requires a dedicated thermal link to the heat bath. This strongly suppresses the thermal signal, which makes the non-thermal phonon component our dominant information carrier.

The thermometer’s temperature rise can therefore be written as the ratio of the time-constants of the two competing processes that reduce the non-thermal phonon population: 1) the absorption in the thermometer with a time-constant , and 2) the thermalization of non-thermal phonons at the crystal surfaces with a time-constant

(5)

It should be noted that this is only valid in the limit , which is equivalent to the statement that collection by the thermometer film does not influence the non-thermal phonon lifetime (21). All devices considered here operate in this regime. Under these conditions, the temperature signal is not influenced by the presence of the thermometer, and thermometer optimization can be considered separately from a change in absorber parameters.

For the absorber scaling law, we keep only the quantities that depend on absorber properties. The energy threshold of the device is inversely proportional to the temperature rise, so we can write

(6)

This is the basis for our scaling law which only considers varying absorber material, geometry and mass. Under these changes, scales with the average time between surface scatterings of the non-thermal phonons, which can be written

(7)

in terms of the mean phonon free path in the crystal and the mean phonon group velocity . For a fixed thermometer surface area, scales with the crystal volume and the mode-averaged absorption rate, like

(8)

is the volume spanned by the phonon modes that cross the thermometer surface per unit time and thermometer area, times the transmission probability into the thermometer. The different dimensionality, (i.e. vs. ), in the scaling laws, arises from the fact that the crystal surface area scales up with the system dimensions, whereas the thermometer area does not.

In total, the scaling law is:

(9)

The first part is purely geometric, while the second represents material parameters. The threshold of CaWO detectors is expected to be 1.72 higher than AlO of same geometry (17), while Si (1.42) and Ge (1.15) fall between these two. The scaling of two detector geometries as a function of mass are considered here. 1) For cubes of side length , and , so that which yields finally . 2) For plates of area and fixed thickness , . In the relevant range, , rises slowly from to (from MC simulation). Roughly, we can take , which also gives , but a different mass-scaling .

With values for found by Monte Carlo methods for each occurring detector geometry, the model can be used to describe the thresholds of various CRESST-type detectors. Since the model can only predict a scaling under change of absorber properties, the absolute normalization (depending e.g. on the noise level of the setup) has to be taken from the respective experiment. In the following, the noise level of the ”benchmark” CRESST setup at LNGS is considered. In Fig. 2, the model predictions for plate and cube detectors are shown as a function of detector mass, fitted to the thresholds achieved in CRESST-II CaWO detectors (green triangles) with a mass of  g (22); (18) and a sapphire cube of 262 g used in CRESST-I (blue cross) (23); (24). The model successfully predicts the energy threshold of CRESST-II light detectors studied in (17) (purple dots), which are sapphire discs with a mass of 2.2 g (diameter: 40 mm, thickness: 0.45 mm) and also the thresholds of  g CRESST-III detectors as expected from a prototype measurement (green error bar) (25). The capability to extrapolate calorimeter thresholds for different detector geometries and materials over orders of magnitude in mass can be applied to the component design for the fiducial-volume cryogenic detector. Red stars indicate the calculated performance of the calorimeters studied here.

Figure 2: Overview of the performance of different calorimeters: the nuclear-recoil energy threshold is plotted vs. the detector mass. The bands (1 ) show predictions of the performance model for CaWO (green) and AlO (blue) calorimeters (see main text). The model is fitted to data of existing CRESST-type detectors (green and blue) (17); (22); (18); (23); (24); (25). Red stars indicate the predicted performance of the detector components studied in this work. The prediction bands for the new Si calorimeters are not shown for clarity.

2.3 Design of the target calorimeter array

For the research program proposed in this work, the target calorimeter has to fulfill the following requirements:

  • A nuclear-recoil energy threshold of  eV).

  • Rates of /[kg day] are expected from CNNS at the sites studied here, as will be shown in chapter 3. Corresponding to this rate, a total target mass of  g is needed for the detection of CNNS.

  • Lowest thresholds require a sufficiently low event rate in the calorimeter. To limit the pile-up contribution to a level of (), a maximum event rate of (0.1 Hz) per detector is allowed given the typical (thermal) pulse decay times of ms (21).

  • The CNNS cross-section is proportional to the target’s neutron number squared, which highly favors heavy elements. On the contrary, the use of light nuclei facilitates a characterization of neutron backgrounds.

Considering these design features, cubic target crystals with a edge length of 5 mm equipped with a tungsten thin-film TES are ideal. A multi-target approach with a variety of elements is chosen which has great advantages for the separation of signal and background through characteristic interaction strength. Cubes of CaWO (0.76 g), AlO (0.49 g), Ge and Si crystals, which are well-known for their excellent cryogenic detector properties (18); (24), are suitable candidates. The performance model (see chapter 2.2) predicts energy thresholds of  eV for AlO and  eV for CaWO (see red stars in Fig. 2). To obtain the desired total target mass, a 3x3 detector array is foreseen (see Fig. 4). This corresponds to a total target mass of 6.84 g for the CaWO and 4.41 g for the AlO array, respectively.

For the temperature sensor, a TES is chosen similar to that which is used for the CRESST detectors (25). The TES consists of a thin W film (thickness 200 nm) with an area of 0.0061 mm and an Al phonon collector with an area of 0.15 mm attached to it (see Fig. 3). The latter increases the collection area for phonons without the penalty of increasing the heat capacity of the sensor (26) yielding an increased pulse height. The TES is weakly coupled to the heat sink via a thin Au stripe (0.01x7.0 mm, thickness: 20 nm) providing a thermal conductance of pW/K (at a temperature of 10 mK). Al and Au wire bonds with a diameter of 25 m are used to provide the electrical contacts for the TES (bonded on the phonon collectors) as well as the ohmic heater (separate bond pads), and the thermal link to the heat sink, respectively. Typically, bias currents between 100 nA and 5 A are applied to the sensor. The resistance change of the TES can be measured with a SQUID system similar to the one in the CRESST dark matter experiment (27).

Figure 3: Schematic view of the sensor design for the target calorimeter. A thin W film (thickness 200 nm) weakly coupled to the heat sink via a (0.01x7.0 mm thickness 20 nm) Au stripe. An Al phonon collector is attached to the W film to increase the collection area of the sensor. The readout current on the TES, the signals on the separate ohmic heater, and the thermal contact are connected via Al and Au wire bonds, respectively.
Figure 4: Technical drawing of the calorimeter array. 3x3 cubic calorimeters (e.g. CaWO, AlO) are installed between two dedicated Si wafers (b). The contact area to the cubic crystals is realized by pyramides (height m) which are produced by wet chemical etching. The inset shows a microscopic picture of a prototype pyramidrical structure. The outer Si wafers (a,c) act as holding structure and host the electrical wiring.

2.4 Results from a prototype calorimeter

In the framework of this project, a prototype AlO calorimeter of 0.5 g has been produced and equipped with a TES according to the design goals described in the previous section. The detector was installed in a copper holder and mounted in a detector test facility at the Max-Planck-Institut for physics in Munich. It consists of a dilution refrigerator in a surface building without dedicated shielding against ambient radioactivity. Further, no shielding against backgrounds from surfaces in the direct vicinity of the calorimeters is used. A Fe X-ray source is placed close to the detector for a calibration of the low-energy region.

In an accompanying paper (9), we present details of a 5.1 h calibration measurement performed with the 0.5 g AlO detector which achieved an energy threshold of  eV. This is independent of the type of particle interaction since it is a calorimetric device. This is the lowest nuclear-recoil energy threshold reported for massive calorimeters, beyond the fundamental nuclear-recoil reach of ionisation-based detectors (28).

The detector operates in the calorimetric mode (see section 2.2), confirmed by the pulse shape. The thermalisation times in the crystal and thermometer film are found to be  ms and  ms, respectively. This ratio fulfills the condition but leaves room for improvements (see below).

The measured threshold is higher (by a factor of 5) compared to what is predicted by the performance model for calorimeters (section 2.2). Part of the discrepancy may be explained due to worse noise level in the MPI setup (by factor 1.5-3 (25)) compared to the low-noise benchmark setup used for the calculation of the predictions (Equ. 9). The considered detector, being the first prototype of a gram-scale calorimeter, is expected to improve by further developments and adjustments of the TES sensor. The ratio of can be further decreased by reducing the thermometer area and accordingly weakening the thermal link. A corresponding reduction of the Al phonon collectors may improve the transport efficiency of quasi particles (26). Furthermore, a moderate reduction of the W-film thickness will reduce the heat capacity of the thermometer without compromising the phonon absorption.

In the calibration measurement a flat background spectrum of  counts/[kg keV day] (7-10 keV) is observed above the calibration peaks (9). This is expected due to the absence of any shielding and can be considered as an absolute upper limit for the total rate in surface experiments (here it corresponds to  Hz). It is comparable to typical total acount rates observed in (1 kg) cryogenic detectors operated in underground laboratories (22). The result clearly demonstrates that gram-scale detectors can be operated in a high-background environment – in particular at surface level – while allowing low energy thresholds and stable operating conditions.

The performance of the prototype fulfills the requirements listed in the previous section in terms energy threshold and operability at surface-level. To demonstrate the required background level – for the near future – measurements with further developed CaWO and AlO calorimeters at low-background experimental sites (e.g. a shallow laboratory) are planned. In particular, the target calorimeter(s) will be embedded in the inner and outer cryogenic shieldings which are discussed in detail in the following sections.

2.5 Low-threshold inner veto and detector holder

Background from the surfaces of the target crystals and surrounding surfaces is a big challenge for rare-event searches, and can limit the sensitivity at low energies. The inner veto provides an active discrimination against beta and alpha decays occurring on surfaces. Typical Q-values of such decays are between  keV and MeV typically shared between 2,3 or more product particles leaving the interaction point in different directions. In a configuration where the target is surrounded by a active veto, the total energy of the reaction is detected (apart from the energy transferred to neutrinos in beta decays). In this way, a high fraction of such backgrounds can be rejected by coincident events in the veto. The rejection of surface background is crucial in particular when approaching ultra-low energy thresholds, as can be seen in experimental data (see e.g.(24)). Fig. 5 shows a section view of the inner part of the detector. In the following, the functionality of the relevant detector components is briefly discussed:

Figure 5: Schematic side view of the detector array. Red: target calorimeter cubes (a) of (5x5x5) mm with a TES each. Blue: Si wafers instrumented with TESs providing a 4 surface veto (b,c,d,e,f). Two slabs have pyramides (b,c) with a height of m which are produced by wet chemical etching. These structure holds the target crystals. Slab b is flexible due to a thickness of only 200m. Purple: Sapphire balls with a diameter of 1 mm. Green: Two outer Si slabs (g,h) of 2 mm thickness press together the inner part. Slab b thereby acts as a spring. Slab h hosts the electrcal wiring which is connected to the TES with wire bonds (black).
  • Target (red): The detector consists of 9 target calorimeters (a in Fig. 5) arranged in a 3x3 detector array. Each crystal is equipped with a TES (see chapter 2.3).

  • Active components (blue): To realize a 4 veto against surface backgrounds, Si wafers read-out by a TES each are used (b-f). Two of these (b and c) are in contact with the target crystals via pyramidal Si structures on the wafers. The upper one (b) is thin enough (m) to be flexible - the wafer acts as a spring. Pressed to the target crystals, the thin wafer realizes a spring-loaded holding structure which can compensate for thermal contraction of the various components of the detector. Possible events induced by the detector holder (e.g. by thermal-stress relaxation) can be rejected since they induce also phonon signals in the TESs of b and c.

  • Passive components (green): Two Si slabs (g and h) are used as support structures for the calorimeter array. They are attached to each other by 4 posts (shown in Fig. 4) providing the necessary pressure to hold the target crystals. The lower wafer (h) is equipped with Al (Au) wiring for the electrical (thermal) connection of the target calorimeters and the inner veto devices.

The inset in Fig. 4 shows a prototype Si wafer with a pyramid structure produced at the Halbleiterlabor of the Max-Planck-Society. The structure is defined by photo-lithography techniques and the pyramid structures are then realized by wet chemical etching.

The rejection power against surface related background was estimated with a dedicated Monte Carlo (MC) study performed with the Geant4 code in version 10.2p1 (30); (31). We follow the recommendation of the Geant4 Low Energy Electromagnetic Physics Working Group (32) and implement the low energy behaviour of electromagnetic interactions via the Geant4 class G4EmStandardPhysics_option4, a selection of most accurate models. Furthermore, we enabled the atomic de-excitation via emission of fluorescence photons and Auger electrons. With one exception, we applied a production cut of throughout our geometry, i.e. for energies above this cut new secondary particles can be created in the simulation, whereas lower energies are directly deposited. The exception are fluorescence photons and Auger electrons which are produced in any case.

Exemplary for an surface contamination we simulated the -decay of by placing the lead ions at rest on the inner surface of the inner veto, facing one target calorimeter made of . The source activity is assumed to be , the maximal external -activity observed with TUM40, a module with especial low background operated in CRESST-II phase 2 (33). The black histogram in Fig. 6 shows the background spectrum seen by the target with inactive inner veto, the red histogram shows the spectrum of the remaining background in case of an active veto with a threshold of . Clearly a reduction of more than two orders of magnitude is feasible. A more detailed MC study of the complete detector array is underway and intended for future publication.

We note that the step at (Fig. 6, black histogram) is no artefact of the used production cut. Instead, it is caused by Coster-Kronig transitions as part of the atomic relaxation subsequent to the decay of the -level of to which decays in 84% of the cases (34).

Figure 6: MC simulation of the expected background from a contamination of the inner surface of the inner veto with a surface -emitter (). The histograms shows the energy deposits in the target for two cases: a passive inner veto (black) and an active inner veto (red) with a threshold of . The inset zoom to the first . Clearly a background reduction of is feasible at low energies

2.6 Outer-veto detector

Given the smallness of the calorimeter array and the inner veto system, these components can be embedded in a large cryogenic outer veto. We consider cylindrical crystals with a diameter and height of (10 cm) which are segmented into two (or more) parts with a central cavity to host the inner detector parts (see Fig. 7). Each crystal of the outer veto is instrumented with a TES. It is foreseen to use materials that are known for their excellent phonon properties, such as e.g. Ge and CaWO, and that have been demonstrated as cryogenic detectors with masses of (100 g-1 kg). Thresholds between 300 eV and 1 keV are reached with such devices, in agreement with the prediction of the performance model for calorimeters in section 2.2 (Fig. 2). CaWO is the preferred material: it has the heavy element W which provides a high cross-section for gamma radiation and the relatively light element O for an efficient moderation of neutrons. The simulations below are therefore performed using CaWO. However, when scaling up the number of detectors (see section 2.7) larger diameters of CaWO crystals are necessary which currently are not available. In this case, Ge crystals are a promising alternative, since those are readily produced in large diameters (up to 300 mm), with high radiopurity.

Figure 7: Technical drawing of the fiducial volume detector. Two calorimeter arrays are installed inside the CaWO outer veto with diameter of 10 cm. The veto is made of two parts with a height of 5 cm each which are equipped with TESs and operated as cryogenic detectors.

It is worth mentioning that the timing information of pulses in the cryogenic detectors is crucial for the efficiency of a coincidence veto. The precision with which the onset of the pulses can be determined defines the dead time in the target calorimeter. We know from neutron scattering experiments that the pulse onset of comparable cryogenic calorimeters can be determined with a uncertainty of s (35). Even an excessive rate of 100 Hz in the veto detector would introduce only a negligible deadtime of 0.1%.

Also the rejection power of the outer veto was estimated with a MC study. Here a target was placed inside the nested shields of inner and outer veto. As typical background we investigate gamma rays following the remaining spectrum at the Dortmund Low Background facility (36), a low-background site at the surface which will be discussed in section 3. FIG. 8 shows as black histogram the background spectrum observed by an unshielded target, in blue the remaining background in case of a passive outer veto, and in red the remaining background in case of an active outer veto with a threshold of . Even with only a passive veto a background supression of more than 3 orders of magnitude at low energies is feasible. Activating the outer veto increase the supression to more than 4 orders of magnitude. Importantly, the expected gamma-induced electron-recoil spectrum remains flat down to energy threshold (see inset of Fig. 8).

For a first estimate of muon-induced neutron backgrounds, a basic MC simulation was performed. Using an active CaWO outer veto, the neutron background is reduced by a factor of , independent of the recoil energy (studied in the energy range from 10 eV to 300 keV). By a clever combination of passive shielding elements like borated polyethylene, and active shielding elements like instrumented plastic or liquid scintillators, and LiF crystals, neutron background levels can be further reduced. This concerns shielding systems placed outside the cryogenic setup surrounding the cryostat at all sides. In addition we provide two technologies to further reduce and reject this potentially harmful background: 1) the outer cryogenic veto system described above and 2) the active background discrimination by the multi-target approach. The latter might be a powerful tool to reduce ultimate backgrounds, particularly neutrons. This is described in more detail in section 3. Nevertheless, we conclude that a dedicated MC simulation using measured muon spectra in combination with a calorimeter measurement at the experimental site are necessary. This is beyond the scope of this work and will be subject of a future publication.

Figure 8: MC simulation of the expected energy deposit in case of a -background similar to the remaining one in the Dortmund Low Background facility (36). The histograms show the energy deposits in the target for three cases: without any veto (black), in case of a passive outer veto (blue), and in case of an active outer veto (red) with a threshold of . The inset zoom to the first . Clearly a background reduction of at lowest energies is reasonable.

2.7 Production and scalability

A disadvantage of cryogenic detectors when compared to e.g., scintillation detectors has always been the difficulty to scale up the experiments in size. The new detector concept presented here overcomes most of these problems. In principle, the detector has been designed such, that the number of production steps of the individual detector components are independent of the number of target calorimeters involved.

The target calorimeters are produced from wafers with a thickness of 5 mm and variable diameters (CaWO up to 60 mm, AlO up to 200 mm, and Si up to 300 mm). With well-established techniques of the semiconductor industry, as e.g. photolitography, thin-film evaportation, etching or sputtering, the TES sensors are being simultaneously equipped on each target calorimeter, and the wafer is cut only afterwards into the individual (5x5x5) mm crystals. The same up-scaling is possible for the inner veto (section 2.5) which acts as a detector holder. It is entirely produced by the above-mentioned methods. The cutting of the wafers is done by means of a laser or other automated methods. The cabling for a large amount of TES sensors are implemented by photolitography in combination with sputtering on the inner veto wafers as done for the 3x3 array. Further, it has been shown (e.g. in (37)) that large amounts of SQUIDs can be realized by SQUID multiplexing.

For the first phase of the experiment, we focus on the production of 3x3 arrays with moderate requirements of size and channel numbers which is foreseen as sufficient for a discovery of CNNS (see below). In a second step, the technology mentioned above enables experiments up to the kg-scale with energy thresholds of (10 eV); an exposure allowing precision measurements of the CNNS cross-section and interesting BSM physics. Fig. 9 shows a technical drawing of a future calorimeter array of 225 crystals which correspond, using AlO, to a total mass of  g.

Figure 9: Technical drawing of a up-scaled calorimeter array using state-of-the-art wafer sizes of 150 mm diameter (e.g. AlO as target and Si as holder). In one production step, a total target mass of  g can be achieved using an array of 225 crystals.

3 Detection of Coherent Neutrino Scattering

3.1 Case 1: At a nuclear power reactor

Signal expectation

Nuclear power reactors are among the most intense (anti-) neutrino sources on earth and therefore a highly interesting site for the detection of CNNS.

A benchmark reactor with a thermal power of 4 GW, a typical value for a two-core reactor plant, yields fissions per second and an isotropic neutrino rate of  s (38). The neutrino flux can be calculated as

(10)

with the distance to the core , the fraction of the fuel component and the respective normalized neutrino-energy spectrum . Fig. 10 shows the neutrino flux for a standard fuel composition (62% of U, 30% of Pu and 8% of U (39)) from a 4 GW reactor at a distance of  m from the core.

Figure 10: Anti-neutrino flux from a benchmark pressurized-water nuclear reactor with a thermal power of 4 GW at a distance of 15 m. A standard fuel composition is used for the calculation: 62% of U, 30% of Pu and 8% of U. (39)

The differential recoil spectrum in the detector can be written as

(11)

using Equ. 1. is the number of target nuclei and the smallest neutrino energy leading to a recoil of a nucleus with the mass .

Figure 11: Count rates on CaWO (red) and AlO (green) expected from a benchmark nuclear power plant of 4 GW for the 3 experimental sites considered. The black dotted lines indicate different background levels (extrapolation to lower energies) measured in different experimental sites. From top to bottom: a) the Stanford shallow underground facility (40), b) the low-background setup at the ARC in Seibersdorf (41), c) the Dortmund low-background facility (36) and d) the Heidelberg shallow laboratory (42). The full black line (e) shows the expected (simulated) background level using the outer and inner veto of the fiducial-volume cryogenic detector. The grey band indicates the uncertainty of the background level with a lower limit at the intrinsic background level of CaWO crystals measured at LNGS (33). Reactor-correlated backgrounds are considered as negligible at the considered distances from core.

The differential recoil spectra of coherently scattered anti-neutrinos in CaWO and AlO detectors at different distances from the core of the benchmark reactor plant are shown in Fig. 11. Due to the dependency of the CNNS cross section (see Equ. 1), the heavy element W boosts significantly the rate on CaWO below 100 eV (full red line) to /[kg keV day]. The rate expected for AlO (full green line), however, stays almost constant at a value of /[kg keV day] below eV. The rates for m are about a factor of 7 lower (dashed lines). The strong material dependence of the CNNS rate is a powerful tool to discriminate the signal from irreducible backgrounds. The signal rate is significantly different for CaWO and AlO, e.g. at 10 eV the ratio is . In contrast, the background counts from external gamma radiation is comparable (within a factor of ). Further, similar neutron background spectra are expected since in both materials - for neutron induced scatters - dominantly O scatters are above energy threshold due to kinematics.

d [m] [eV] counts/[kg day] counts/[array day]
CaWO AlO CaWO AlO

15 5 790.3 112.8 5.44 0.51
10 690.2 110.1 4.75 0.49
20 547.2 105.4 3.77 0.47
40 5 111.1 15.9 0.77 0.07
10 97.1 15.5 0.67 0.07
20 77.0 14.8 0.53 0.07
100 5 17.8 2.5 0.12 0.01
10 15.5 2.5 0.11 0.01
20 12.3 2.4 0.08 0.01

Table 1: Integrated CNNS count rate from a nuclear reactor with a total thermal power of 4 GW at different distances between and 5 keV. The rates are integrated up to 5 keV.

The integrated count rates for different energy thresholds and distances are listed in Table 1. The signal is integrated up to an energy of 5 keV where the contribution to the signal is negligible. Count rates are given per day and kg as well as per day and detector array (CaWO: 6.84 g, AlO: 4.41 g). A signal rate of up to  counts per array and day is expected for CaWO target calorimeters.

Due to the relatively high rates predicted at such sites, the detection of CNNS with a small-scale detector of low threshold (10 eV) at a moderate distance from the core is clearly in reach.

Background level

We consider a shallow experimental site with a small overburden to shield against cosmogenic backgrounds at m from the reactor core. Possible candidate sites are, e.g. a room in the basement of a building outside the reactor containment, an artificial overburden outside the reactor building or even a site outside the reactor plant. At such places, the reactor-correlated gamma and neutron backgrounds are considered as negligible due to the large distance and significant shieldings. In the following we concentrate on uncorrelated backgrounds which at shallow sites are dominated by muon-induced events (36). Plenty of experimental data describing detectors operated in shallow or above-ground low-background environments exist in the literature, mostly for Ge detectors. The following total background levels are reached in selected experiments: 0.4 counts/[kg keV day] at the shallow underground lab in Heidelberg (42), 5 counts/[kg keV day] at the Dortmund low-background facility (36),  counts/[kg keV day] at the ARC in Seibersdorf (41) and 200 counts/[kg keV day] in the CDMS experiment operated at the Stanford underground facility (40) (black dotted lines). All values correspond to the rates in the lowest energy bin of the respective experiment. The grey band in Fig. 11 indicates the uncertainty in the observed background level depending on the individual site, the overburden and the shielding design. The lowest energy threshold (eV) among the listed experiments is achieved by CDMS (40).

We use the highest background level reported as a conservative upper limit for the sensitivity studies. Even more conservative, we do not consider the additional background-rejection capability of the inner and outer cryogenic veto. As shown in chapter 2.5 and 2.6 by a dedicated MC study, the cryogenic fiducial-volume detector reduces surface, gamma and neutron backgrounds by factors of and , respectively, in the target volume. In the following, the (flat) background rate of 200 counts/[kg keV day] is referred to as the benchmark.

In case of CaWO the CNNS signal is orders of magnitude above the conservative benchmark background whereas in case of AlO the signal-to-background ratio is much smaller (factor of 1-5), see Fig. 11. The multi-target approach, therefore, is a powerful tool to actively discriminate neutrino-induced signals from backgrounds. In particular, it allows to identify possible ultimate exponentially shaped, signal-like backgrounds.

Experimental site and discovery potential

An extensive likelihood study is performed to investigate the discovery potential of CNNS with the proposed small-scale experiment. We consider one CaWO (total mass: 6.84 g) and one AlO (total mass: 4.41 g) calorimeter array inside the inner and outer active cryogenic veto (see Fig. 7). The benchmark background level is assumed and, conservatively, the rejection capability of the surface veto is not used for the background estimation. Three different thresholds are studied (5, 10 and 20 eV) which, however, have only a minor impact on the discovery potential. We define three scenarios:

  • Near case: A distance of 15 m from the reactor core - a site within the reactor containment. Highest count rates are expected, but there are tough requirements for the shielding against correlated backgrounds. The access is restricted and strict safety regulations have to be considered.

  • Medium case: A distance of 40 m from the reactor core - outside the containment and the reactor building. Possibly a shallow site in a adjoining building or a dedicated site with an artificial overburden. Easier access and a better infrastructure.

  • Far case: A distance of 100 m from the reactor core - far away from the critical reactor components, possibly outside the entire power-plant area. Straightforward access and plenty of possible sites.

Figure 12: Artist view of a typical nuclear power plant with possible experimental sites (red boxes) for the 3 different scenarios (see text).
Figure 13: Discovery potential of CNNS vs. time at a 4 GW reactor core from the likelihood ratio analysis described in the text. The combination of one CaWO and one AlO calorimeter array is investigated assuming the benchmark background level of  counts/[kg keV day]. The full lines indicate the median discovery probability for a nominal energy threshold of 10 eV, the bands show the 90 % confidence intervals. Three cases for the reactor distance are considered (see text): near case (blue), medium case (green) and far case (grey). Varying the threshold to 5 and 20 eV, respectively, has only a minor impact on the discovery potential (see text).

For each case, spectra are randomly generated for a large number of varying exposures. The results of this MC simulation are studied with a likelihood ratio analysis. In every MC experiment, one spectrum each is generated for the CaWO and AlO arrays. The unbinned likelihood of a model’s parameters is calculated as a product over the individual likelihoods for each event in both spectra and the Poisson likelihood for observing this total event numer (Extended Maximum Likelihood method). The single event likelihood is proportional to the sum of the signal and background rates for the given parameter values. Two very simple models are considered: the free model has two parameters, namely the level of the flat background and the strength of the CNNS signal relative to the standard model expectation. In the null model, the CNNS signal strength is held at zero. The maximum likelihood of each model at the best fit parameter values is denoted and respectively. Since the two models are nested with one additional parameter in the free model, the likelihood ratio test statistic

follows a -distribution with one degree of freedom (by Wilks’ theorem). The square root of the test statistic therefore follows a standard normal distribution, so that the statistical significance in of the claim of a CNNS signal with nonzero cross-section in addition to the assumed flat background is directly given by for each pair of spectra.

Fig. 13 shows the resulting discovery potential of the 3 scenarios. The full lines indicate the median discovery potential as derived from the MC data, using an energy threshold of 10 eV. The bands give the 90 % confidence intervals. All three scenarios show a very promising potential for the discovery (5) of CNNS - in the near case within  day, in the medium case within  weeks and in the far case within  year of measuring time. Improving the threshold to 5 eV reduces the measuring time necessary for a 5 discovery by a factor of , in average for the three scenarios discussed. For a threshold of 20 eV, times longer measurements are required.

Systematic deviations of background and signal rates have only minor influence on the discovery potential: Repeating the simulations with 20 % higher and lower background level yield 10 % higher and lower times to discovery, respectively. A 5 % stronger signal makes discovery faster by 5 %, while a 5 % weaker signal requires 7.5 % more measuring time.

To study the impact of a non-flat background on the discovery potential we use data from a CDMS detector operated at a shallow laboratory (40). The measured spectrum was fitted with an exponential below 10 keV and extrapolated exponentially to beyond the energy threshold of 500 eV. The fitted background model corresponds to an exponential component rising above a flat background of 2.5/ keV kg d at around 2 keV and reaching 700/ keV kg d at zero energy. Using this background level in the likelihood study has only a minor impact. The measuring time required for a 5 discovery in the three scenarios increase moderately by a factor of .

Background studies including dedicated measurements on the individual sites and detailed MC simulations are required to find the most suitable site. At the medium and far sites, expected backgrounds are rather straightforward, while for the near site a proper understanding of the possibly remaining reactor-correlated backgrounds is needed. The near site, however, would - despite a rapid discovery of CNNS - allow a precision measurement (statistical error on a percent level) of the cross-section predicted by the Standard Model within a measuring time of one year. Impressively, this can be performed by a detector with a total target mass of  g, given the necessary control of systematics.

3.2 Case 2: At a thermal research reactor

To study the possibility of detecting CNNS at a thermal research reactor, both the signal and background spectrum were adapted to the altered conditions. The signal expectation was calculated for the fuel composition found at FRM2 (96% of U, 0% of Pu and 4% of U) (43) which does not change the signal shape appreciably. Our fiducial model is a 20 MW reactor at a distance of 3 m to the detector, which corresponds to a neutrino flux reduced by with respect to the medium distance case at the power reactor. The close proximity to the research reactor implies great uncertainty regarding the composition and shape of backgrounds correlated to the reactor power, and thus to the neutrino signal. We use the results of a detailed background measurement and simulation campaign (7) by the MINER collaboration for our background estimate. In the framework of the MINER experiment it is planned to deploy an array of Ge cryogenic detectors in close proximity to a 1 MW research reactor (7). Between the reactor core and the detectors, several layers of shielding made of graphite, borated high-density polyethylene, Pb and Cu are placed. These conditions are assumed to be comparable to the setup we are investigating in this context. Fig. 14 shows the expected CNNS rates at the research reactor at a distance of 3 m from the core, along with the neutron (scaled to 10%, see below) and gamma backgrounds from (7). Compared to those, cosmogenic backgrounds and environmental radioactivity appear to be sub-dominant in shallow labs (compare to section 3.1.2) and are neglected in this basic feasibility study.

Figure 14: Expected CNNS count rates at the FRM2 research reactor with a thermal power of 20 MW at a distance 3 m. The isotopic fractions of the neutrino emission is adjusted to 96% from U and 4% from U to account for the different fuel composition at FRM2. The thick red and green lines show the CNNS rates on CaWO and AlO, respectively. In black, the simulated research-reactor background spectra for neutrons (scaled to 10%, see text) and gammas from (7) are shown.

The similarity of the signal shape to the reactor correlated neutron background makes the detection of CNNS challenging in this environment even with extensive reactor ON/OFF measurements. Under these conditions, a multi target approach can be beneficial because of the material dependence of CNNS (through the cross-section ). In both CaWO and AlO, neutrons are expected to scatter predominantly off the light oxygen nuclei, leaving a comparable signature. On the other hand, the AlO array does not contribute meaningfully to the CNNS signal measurement, but yields an important measurement of the background rate as a function of energy. This information helps to break the degeneracy between the CNNS signal and the neutron background, which are very similar in shape.

Fig. 15 shows the expected counts in the AlO array (left) and in the CaWO array (right) after one year of measurement in the described conditions. The similarity of the signal and background shapes in CaWO is apparent. The signal in AlO is overwhelmed by the Poisson fluctuations, so the AlO array only contributes to the determination of the background level.

Figure 15: Expected spectra at a thermal research reactor after one year of measurement in the AlO array (left) and in the CaWO array (right). Shown are the expected background counts (blue) from the MINER background scaled to 10%, and the expected CNNS signal counts (in red). The error bars show the expected fluctuations. The background spectra are assumed to be identical in the two materials, scaled only by the respective exposure (higher in CaWO due to the higher density of the material). The CNNS signal is strongly enhanced in CaWO due to the neutron-rich W nucleus. The AlO array with a similar neutron response can serve for an in-situ background characterization.

For the likelihood study, we assume a fixed live time of one year with the CaWO and AlO arrays and show the detection significance (computed as above) as a function of the background level instead. The simulated background spectrum is a scaled version of the MINER neutron background. The shape of the background spectrum is assumed to be known for the likelihood models. To show the added benefit of the AlO array for background characterization, we plot the significance obtained by each detector material separately (green, blue) along with the combined significance (obtained as above, black) in Fig. 16. Full lines indicate the median discovery potential, the bands are 90 % confidence intervals. The background level has to be reduced significantly with respect to the level reported in (7) to allow a -detection within one year. The “background-only” information provided by the AlO array considerably relaxes the background requirements, so that a detection after one year becomes very likely with a background below of the MINER neutron background, and feasible with a background below .

Figure 16: CNNS discovery potential at a research reactor vs. background strength. The background is based on a detailed neutron simulation performed within the MINER collaboration (Fig. 14 in (7)). The detection significance is plotted separately for the AlO (green) and CaWO (blue) arrays. The combination of both (black) considerably enhances the detection significance. Full lines represent the median discovery potential, the bands constraint by thin lines are 90 % confidence intervals as derived by the MC simulation.

3.3 Case 3: Using a neutrino source

Figure 17: Neutrino flux from the Ce neutrino source with an assumed activity of 150 kCi. The low-energy spectrum up to  keV corresponds to the initial Ce decay, with the high-energy spectrum up to  MeV stems from the decay of the daughter nucleus Pr.

The detection of CNNS using a radioactive neutrino source is a scenario which poses quite different challenges compared to the other considered cases. To evaluate the new situation, we assume a neutrino source similar to the source proposed for the SOX experiment (44), specifically we show in Fig. 17 the neutrino spectrum of a Ce source with an initial activity of 150 kCi. The low-energy neutrinos ( keV) originate in the initial decay of Ce with a half life of 285 days, while the broad neutrino spectrum up to a Q value about 3 MeV originates in the fast decay of the daughter nucleus Pr. The low-energy neutrinos do not produce a detectable W recoil in the detectors considered here. With a 10 eV threshold, a CaWO detector is sensitive mostly to neutrinos above 1 MeV. Therefore, even with an optimistic shielding scenario (source distance of 1 m) to stop residual gammas and neutrons from source impurities, the recoil rates are more than an order of magnitude below the far case of the power reactor scenario, as shown in Fig. 18.

Figure 18: Count rates expected from a benchmark Ce neutrino source with an initial activity () of  Bq. The dotted lines indicate the worst-case scenario with a distance of 1.5 m and a measurement after two half-lifes of the source ( y). The black dotted line indicates the background level achieved with CaWO in the CRESST setup at LNGS, Italy (33).

The background level shown there is the measured intrinsic background of CRESST crystals, which constitutes the best reasonably achievable background level. We point out that any parasitic, radioactive contamination of the neutrino source which may produce additional neutron and gamma background is neglected in this study. With this optimistic assumption, detection of CNNS becomes feasible, but still suffers from low count rates.

Fig. 19 shows the discovery potential of CNNS at 1 m from the benchmark radioactive neutrino source as a function of the exposure collected at full activity (150 kCi), obtained with the likelihood ratio method described above. Detection comes in reach with an exposure of  kg d. Such exposures are feasible, however require a larger detector mass. The exposure has to be collected within a few half-lives of the source isotope ( d), which necessitates a larger array, e.g. cubes ( g).

Figure 19: Discovery potential of CNNS vs. exposure using a neutrino source, derived by a dedicated likelihood analysis. The full line represents the median discovery potential, the band is the 90 % confidence interval. An exposure of 10 kg d is required for a detection of CNNS.

4 Summary and Outlook

The smallness of gram-scale calorimeters offers the following significant advantages: 1) very low energy thresholds down to the 10 eV-regime and presumably below, 2)  the possibilty of an encapsulation of the small calorimeters by other cryogenic devices which act as anti-coincidence vetos and 3)  the ability to operate the detectors at the surface in a relatively high-rate environment. These advantages are demonstrated experimentally by a measurement with a prototype 0.5 g AlO calorimeter which reaches a threshold of  eV. Based on that, we propose a new detector concept here: a gram-scale fiducial-volume cryogenic detector for the detection of CNNS, called the -cleus experiment.

A basic version of the detector, consisting of two 3x3 calorimeter arrays made of CaWO and AlO crystals with a total mass of  g, has a high potential for a rapid discovery of CNNS. We study various experimental scenarios for this new technology: an installation at a nuclear power plant, at a thermal research reactor and close to a radioactive neutrino source. We conclude that all three methods allow a detection of CNNS, however the first scenario clearly shows the highest potential.

We investigate the operation of the 10 g detector at a distance of 40 m from a nuclear power reactor with a thermal power of 4 GW. This corresponds to an experimental site outside the reactor containment, and is therefore rather straight-forward in terms of background levels, infrastructure and access. As shown by a dedicated likelihood analysis, the rate is still sufficiently high to achieve a 5 discovery within a measuring time of  weeks.

The detector placed at a well-shielded site within the reactor containment, e.g. at a distance of 15 m from the core, would give the unique possibility for precision measurements of the CNNS cross-section and probe, e.g. the Weinberg angle at low momentum transfers  (10). Furthermore, since a discovery of CNNS is possible within a day, this technology can be used for real-time monitoring of nuclear reactors. Such a small-scale experimental setup could provide a comprehensive surveillance system for non-proliferation and accident prevention at the nuclear reactors world-wide.

The detector is designed such to be scaled up in a relatively simple way due to the use of production techniques of the semiconductor industry. We point out that total target masses of (1 kg) are feasible with the design principle given here. This enables new approaches in rare events searches, such as the search for MeV-scale dark matter and flavour-independent precision measurements of the solar neutrino flux. Operating such a kg-scale detector with energy thresholds in the 10 eV regime at a power reactor opens the door to a variety of new physics beyond the Standard Model of Particle physics.

References

  1. D.Z. Freedman, Phys. Rev. D 9, 1389 (1974)
  2. A. Drukier, L. Stodolsky, Phys. Rev. D 30, 2295 (1984)
  3. D. Akimov, et al., (2015). URL https://arxiv.org/abs/1509.08702
  4. A. Aguilar-Arevalo, et al., J. Phys. Conf. Ser. 761(1), 012057 (2016). DOI 10.1088/1742-6596/761/1/012057
  5. G. Fernandez Moroni, J. Estrada, E.E. Paolini, G. Cancelo, J. Tiffenberg, J. Molina, Phys. Rev. D91(7), 072001 (2015). DOI 10.1103/PhysRevD.91.072001
  6. S. Kerman, V. Sharma, M. Deniz, H.T. Wong, J.W. Chen, H.B. Li, S.T. Lin, C.P. Liu, Q. Yue, Phys. Rev. D93(11), 113006 (2016). DOI 10.1103/PhysRevD.93.113006
  7. G. Agnolet, et al., ArXiv e-prints (2016). URL http://adsabs.harvard.edu/abs/2016arXiv160902066M
  8. J. Billard, et al., (2016). URL https://arxiv.org/abs/1612.09035
  9. R. Strauss, et al., Phys. Rev. D 96, 022009 (2017). URL {https://inspirehep.net/record/1591648}
  10. M. Lindner, et al., JHEP 2017(3), 979 (2017)
  11. N. Berger, et al., ArXiv e-prints (2015). URL http://adsabs.harvard.edu/abs/2015arXiv151103934B
  12. J. Erler, M.J. Ramsey-Musolf, Phys. Rev. D 72, 073003 (2005)
  13. T. Ohlsson, Reports on Progress in Physics 76(4), 044201 (2013)
  14. A. Santamaria, C. Peña Garay, S. Davidson, N. Rius Dionis, Journal of High Energy Physics, 2003, vol. 0303, p. 011 (2003)
  15. A.J. Anderson, et al., Phys. Rev. D 86, 013004 (2012)
  16. H.T. Wong, et al., Modern Physics Letters A 20(15), 1103 (2005)
  17. J. Rothe, Master’s thesis, LMU München, https://publications.mppmu.mpg.de (MPP-2016-391) (2016)
  18. G. Angloher, et al., Eur. Phys. J. C76(1), 25 (2016)
  19. S.H. Moseley, et al., J. Appl. Phys. 56(5), 1257 (1984)
  20. M. Pyle, E. Figueroa-Feliciano, B. Sadoulet (2015). URL http://adsabs.harvard.edu/abs/2015arXiv150301200P
  21. F. Pröbst, et al., J. Low Temp. Phys. 100(1-2), 69 (1995)
  22. G. Angloher, et al., Eur. Phys. J. C74(12), 3184 (2014)
  23. O. Meier, et al., Nucl. Instrum. Meth. A444, 350 (2000)
  24. G. Angloher, et al., Astropart. Phys. 18(1), 43 (2002)
  25. R. Strauss, et al., Nucl. Instrum. Meth. A845, 414 (2017)
  26. G. Angloher, et al., J. Low Temp. Phys. 184(1), 323 (2016)
  27. G. Angloher, et al., Eur. Phys. J. C72, 1971 (2012)
  28. R. Agnese, et al., Phys. Rev. D 95, 082002 (2017)
  29. R. Strauss, et al., Eur.Phys.J. C74(7), 2957 (2014)
  30. S. Agostinelli, et al., Nucl. Instrum. Meth. A 506(3), 250 (2003). DOI 10.1016/S0168-9002(03)01368-8
  31. J. Allison, et al., IEEE Transactions on Nuclear Science 53(1), 270 (2006). DOI 10.1109/TNS.2006.869826
  32. Geant4 Low Energy Electromagnetic Physics Working Group. Low energy physics lists. https://twiki.cern.ch/twiki/bin/view/Geant4/LowePhysicsLists (2016)
  33. R. Strauss, et al., JCAP 1506(06), 030 (2015)
  34. R.B. Firestone, Table of Isotopes, 8th edn. (Wiley, 1996)
  35. R. Strauss, et al., Eur. Phys. J. C 74(7), 2957 (2014). DOI 10.1140/epjc/s10052-014-2957-5
  36. H. Gastrich, et al., Appl. Radiat. Isot. 112, 165 (2016)
  37. M.A. Dobbs, et al., Review of Scientific Instruments 83(7), 073113 (2012)
  38. V. Kopeikin, L. Mikaelyan, V. Sinev, Phys. Atom. Nucl. 67, 1892 (2004). [Yad. Fiz.67,1916(2004)]
  39. O. Tengblad, et al., Nuclear Physics A 503, 136 (1989)
  40. D.S. Akerib, et al., Phys. Rev. D82, 122004 (2010)
  41. M. Schwaiger, F. Steger, T. Schroettner, C. Schmitzer, Applied Radiation and Isotopes 56(1–2), 375 (2002)
  42. G. Heusser, et al., Eur. Phys. J. C75(11), 531 (2015)
  43. J. Wuttke, et al., Review of Scientific Instruments 83(7), 075109 (2012). DOI 10.1063/1.4732806
  44. M. Vivier, et al., Journal of Physics: Conference Series 718(6), 062066 (2016)
103625
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
Edit
-  
Unpublish
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel
Comments 0
Request answer
""
The feedback must be of minumum 40 characters
Add comment
Cancel
Loading ...