{"id":105896,"date":"2025-07-30T22:05:16","date_gmt":"2025-07-30T22:05:16","guid":{"rendered":"https:\/\/www.europesays.com\/us\/105896\/"},"modified":"2025-07-30T22:05:16","modified_gmt":"2025-07-30T22:05:16","slug":"direct-observation-of-coherent-elastic-antineutrino-nucleus-scattering","status":"publish","type":"post","link":"https:\/\/www.europesays.com\/us\/105896\/","title":{"rendered":"Direct observation of coherent elastic antineutrino\u2013nucleus scattering"},"content":{"rendered":"<p>Data taking<\/p>\n<p>Data from the CONUS+ run 1 analysis presented in this study correspond to the period from November 2023 to July 2024. Three of the four detectors, named C2, C3 and C5, are considered for analysis with a total active mass of 2.83\u2009\u00b1\u20090.02\u2009kg (ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 2\" title=\"Bonet, H. et al. Large-size sub-keV sensitive germanium detectors for the CONUS experiment. Eur. Phys. J. C 81, 267 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR2\" id=\"ref-link-section-d46697375e1304\" target=\"_blank\" rel=\"noopener\">2<\/a>). After selection cuts to remove time periods with unstable noise conditions, deficient radon flushing and a few days with a strong contribution of microphonic events, the exposure considering the active mass is 327\u2009kg\u2009days with reactor on and 60\u2009kg\u2009days with the reactor off.<\/p>\n<p>The evolution of the main environmental parameters during reactor on and off is shown in Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig4\" target=\"_blank\" rel=\"noopener\">1<\/a> for the three detectors. The shape of the noise peak reconstructed in the lowest channels of the data acquisition system (DAQ) was found to follow a Gaussian distribution<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 1\" title=\"Ackermann, N. et al. CONUS+ experiment. Eur. Phys. J. C 84, 1265 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR1\" id=\"ref-link-section-d46697375e1314\" target=\"_blank\" rel=\"noopener\">1<\/a> and its full-width at half maximum (FWHM) was monitored over time with variations below 2\u2009eVee (electron equivalent energy). The dataset was selected to ensure that the noise rate variations are below 20%.<\/p>\n<p>In a dedicated study, we confirmed the correlation between noise rate and cryocooler power. Stability was improved compared with the CONUS setup in Brokdorf (KBR) by replacing the two-fan ventilation system with a water-cooled chiller system for the pulse tube cryocoolers<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 1\" title=\"Ackermann, N. et al. CONUS+ experiment. Eur. Phys. J. C 84, 1265 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR1\" id=\"ref-link-section-d46697375e1323\" target=\"_blank\" rel=\"noopener\">1<\/a>. Even at cryocooler power variations of up to 30\u2009W, no correlation was observed with the count rate in the region of interest above 160\u2009eVee.<\/p>\n<p>Detector response and energy scale<\/p>\n<p>The trigger efficiency was determined by injecting artificial signals produced by a pulse generator with the same rise time as the physical signals. The pulses are injected through a specific circuit implemented in the HPGe preamplifier. A detailed scan allowed us to measure the detector response as a function of the energy<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 1\" title=\"Ackermann, N. et al. CONUS+ experiment. Eur. Phys. J. C 84, 1265 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR1\" id=\"ref-link-section-d46697375e1337\" target=\"_blank\" rel=\"noopener\">1<\/a>. The trigger efficiency remains more than 90% down to 140\u2009eVee for all detectors. The evolution of the trigger efficiency curve parameters during run 1 was studied with different measurements, remaining stable with differences of less than 2% throughout the run.<\/p>\n<p>The energy was calibrated using the binding energies of the K-shell (10.37\u2009keV) and L-shell (1.30\u2009keV) from the decays of 68Ge\/71Ge inside the HPGe diodes, considering a linear behaviour in this energy range. At shallow depth, the 68Ge\/71Ge radioisotopes are continuously produced by cosmic and muon-induced neutrons<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 1\" title=\"Ackermann, N. et al. CONUS+ experiment. Eur. Phys. J. C 84, 1265 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR1\" id=\"ref-link-section-d46697375e1354\" target=\"_blank\" rel=\"noopener\">1<\/a>. Thus, it is possible to monitor the stability of the energy scale during the whole measurement of this in situ activation, observing variations below 2%. A specific 252Cf irradiation was performed at the end of run 1, collecting in 45\u2009days of the measurement more than 5,000 events from the K-shell and 700 events from the L-shell in each detector. In this way, an energy calibration uncertainty below 5\u2009eVee is achieved. The HPGe diode also produces 158\u2009eV X-rays from the binding energy of the M-shell of the 68Ge\/71Ge decays. Using the ratio of the K- and L-shells, a total of 100 events per detector was expected after irradiation with 252Cf. Although an indication of such events is seen, no conclusive signal was observed because of the proximity to the noise edge and the lack of statistics.<\/p>\n<p>The linearity of the DAQ chain in the sub-keVee region was investigated using the pulse generator signals. The results are shown in Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig5\" target=\"_blank\" rel=\"noopener\">2<\/a>. Deviations from a pure linear behaviour are observed below 250\u2009eVee, at which a few eVee variations have a strong impact on the CE\u03bdNS signal. They were attributed to two DAQ-related effects<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 34\" title=\"Bonhomme, A. et al. Direct measurement of the ionization quenching factor of nuclear recoils in germanium in the keV energy range. Eur. Phys. J. C 82, 815 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR34\" id=\"ref-link-section-d46697375e1381\" target=\"_blank\" rel=\"noopener\">34<\/a>. The energy non-linearity was corrected in the CE\u03bdNS analysis and was taken into account during the event energy reconstruction.<\/p>\n<p>The energy resolution of the detectors at low energies was also evaluated with the pulse generator signals. For the C2 detector, a resolution of (48\u2009\u00b1\u20091)\u2009eVee (FWHM) was found, whereas the C5 and C3 detectors have a resolution of (47\u2009\u00b1\u20091)\u2009eVee (FWHM). Further contribution to the total energy resolution comes from the statistical fluctuation of the number of electron\u2013hole pairs produced in the Ge crystals in the case of an event, which is<\/p>\n<p>$$\\Delta {E}_{{\\rm{S}}}=2.35\\sqrt{F{\\epsilon }E}$$<\/p>\n<p>\n                    (1)\n                <\/p>\n<p>with \u03f5\u2009=\u20092.96\u2009eV, the average energy needed to create a single electron\u2013hole pair in Ge, and the Fano factor F\u2009=\u20090.11.<\/p>\n<p>Selection cuts<\/p>\n<p>Selection cuts are applied to reduce background events while keeping the CE\u03bdNS signal. Four different selection cuts are applied to the data. First, the muon-veto system allows for efficient suppression of the impact of cosmic radiation, using a 450-\u03bcs anti-coincidence window between the veto and HPGe signals. The average rate detected in the muon veto during reactor on is (274\u2009\u00b1\u20091)\u2009Hz and decreases to (214\u2009\u00b1\u20091)\u2009Hz in reactor off periods. The corresponding average dead times in both periods are 12.3% (reactor on) and 9.6% (reactor off). Second, inhibit signals are generated when the increasing baseline has reached saturation of the dynamic range of the transistor reset preamplifier for each HPGe detector<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 1\" title=\"Ackermann, N. et al. CONUS+ experiment. Eur. Phys. J. C 84, 1265 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR1\" id=\"ref-link-section-d46697375e1460\" target=\"_blank\" rel=\"noopener\">1<\/a>. An anti-coincidence window of 1\u20132.5\u2009ms (depending on the detector) is applied to veto unwanted spurious HPGe detector signals generated shortly after the resets. This cut suppressed 30% events after the muon veto anti-coincidence at low energy, becoming negligible above 5\u2009keVee. The dead time induced by this cut is calculated, combined with the previously mentioned muon-veto dead time, to avoid the overlapping of both veto windows. An additional dead time of 0.5\u20132.1% (depending on the detector) is estimated. Third, the time difference distribution of events is studied in each channel as proposed in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 2\" title=\"Bonet, H. et al. Large-size sub-keV sensitive germanium detectors for the CONUS experiment. Eur. Phys. J. C 81, 267 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR2\" id=\"ref-link-section-d46697375e1466\" target=\"_blank\" rel=\"noopener\">2<\/a>. Finally, an anti-coincidence cut is applied between different HPGe detectors with a 5-ms time window. The probability of a neutrino interacting with different detectors is negligible, whereas for other backgrounds, such as muon-induced neutrons created in the shield, simultaneous hits in several detectors at once can be expected.<\/p>\n<p>The rejection efficiencies of these selection cuts are summarized in Extended Data Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Tab2\" target=\"_blank\" rel=\"noopener\">1<\/a> for the three detectors in different energy regions. The total dead time induced by the selection cuts is between 12.8% and 14.4% in the reactor on periods.<\/p>\n<p>Background model<\/p>\n<p>The background model used in the analysis of the CONUS+ data is based on Monte Carlo simulations using the Geant4-based framework MaGe<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 39\" title=\"Boswell, M. et al. MaGe-a Geant4-based Monte Carlo application framework for low-background germanium experiments. IEEE Trans. Nucl. Sci. 58, 1212&#x2013;1220 (2011).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR39\" id=\"ref-link-section-d46697375e1487\" target=\"_blank\" rel=\"noopener\">39<\/a>, following the approach in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 29\" title=\"Bonet, H. et al. Full background decomposition of the CONUS experiment. Eur. Phys. J. C 83, 195 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR29\" id=\"ref-link-section-d46697375e1491\" target=\"_blank\" rel=\"noopener\">29<\/a>. A complete decomposition of the background in both reactor on and off data was done for all detectors used in the analysis. In the following, the main sources of background are described.<\/p>\n<p>Cosmogenic neutrons<\/p>\n<p>As described in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 30\" title=\"CONUS Collaboration. Background characterization of the CONUS+ experimental location. Eur. Phys. J. C 85, 465 (2025).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR30\" id=\"ref-link-section-d46697375e1502\" target=\"_blank\" rel=\"noopener\">30<\/a>, the impact of cosmogenic neutrons with energies up to 100\u2009MeV was investigated by first propagating the expected neutron spectrum at the KKL location (based on refs.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 40\" title=\"Goldhagen, P., Clem, J. M. &amp; Wilson, J. W. The energy spectrum of cosmic-ray induced neutrons measured on an airplane over a wide range of altitude and latitude. Radiat. Prot. Dosimetry 110, 387&#x2013;392 (2004).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR40\" id=\"ref-link-section-d46697375e1506\" target=\"_blank\" rel=\"noopener\">40<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 41\" title=\"Gordon, M. et al. Measurement of the flux and energy spectrum of cosmic-ray induced neutrons on the ground. IEEE Trans. Nucl. Sci. 51, 3427&#x2013;3434 (2004).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR41\" id=\"ref-link-section-d46697375e1509\" target=\"_blank\" rel=\"noopener\">41<\/a>) with a flux of (1.4\u2009\u00b1\u20090.2)\u2009\u00d7\u200910\u22122\u2009neutrons\u2009s\u22121\u2009cm\u22122 through a reactor building model. Then the resulting flux is tracked inside the CONUS+ room. This flux was used as the basis for the next simulation, in which neutrons were started isotropically from a half-sphere around the CONUS+ shield, and their contribution to the CONUS+ background was measured. The simulations show a contribution of (21.6\u2009\u00b1\u20093.1)\u2009counts\u2009day\u22121\u2009kg\u22121 in each detector in the energy range between 0.4\u2009keVee and 1\u2009keVee, which corresponds to approximately half of the background counts in this region.<\/p>\n<p>Cosmogenic muons<\/p>\n<p>For the muon simulations<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 29\" title=\"Bonet, H. et al. Full background decomposition of the CONUS experiment. Eur. Phys. J. C 83, 195 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR29\" id=\"ref-link-section-d46697375e1536\" target=\"_blank\" rel=\"noopener\">29<\/a>, the expected flux and muon spectrum at an overburden of 7.4\u2009m water equivalent (w.e.) were calculated from those at the surface of Earth<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 42\" title=\"Bugaev, E. V. et al. Atmospheric muon flux at sea level, underground, and underwater. Phys. Rev. D 58, 054001 (1998).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR42\" id=\"ref-link-section-d46697375e1540\" target=\"_blank\" rel=\"noopener\">42<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 43\" title=\"Reyna, D. A simple parameterization of the cosmic-ray muon momentum spectra at the surface as a function of zenith angle. Preprint at &#010;                  arxiv.org\/abs\/hep-ph\/0604145&#010;                  &#010;                 (2006).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR43\" id=\"ref-link-section-d46697375e1543\" target=\"_blank\" rel=\"noopener\">43<\/a> and propagated through the shield. The resulting spectrum was validated by comparing it with the CONUS+ data without the applied muon-veto cut, which showed good agreement. The muon-veto cut was then applied by multiplying the simulation output by a factor of 0.01 for all energies greater than 2\u2009keVee, corresponding to a muon-veto efficiency of 99%.<\/p>\n<p>For energies below 2\u2009keVee, a different approach was taken. Here, simulations show an inefficiency in the tagging ability of the muon-veto system because of the setup of the shield. The outer muon veto is located under a layer of lead in the CONUS+ shield. As such, it is possible for muons to pass through this outermost lead layer without hitting one of the muon-veto layers. These muons can induce electromagnetic showers in the outer lead layer, which propagate through the shield and are registered in the Ge detectors. However, because no muon passes through any muon-veto layer in such an event, the energy deposition in the plastic scintillator plates will be much lower, resulting in a greatly reduced tagging efficiency of these events. Simulations show that at energies below 0.4\u2009keVee, up to 80% of all muon-induced background comes from these events. Based on this simulation output, this inefficiency was modelled with a polynomial and accounted for in the final muon-veto efficiency. The resulting efficiency drops towards lower energies, with its minimum being 97% below 0.4\u2009keVee. Using this approach, the overall background contribution of the cosmic ray muons was found to be (17.4\u2009\u00b1\u20090.3)\u2009counts\u2009day\u22121\u2009kg\u22121 in each detector in the energy range between 0.4\u2009keVee and 1\u2009keVee, which corresponds to approximately a third of the background counts in this region.<\/p>\n<p>Leakage test component<\/p>\n<p>During the final run of the CONUS experiment at KBR, an additional background component had to be introduced in the background model<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 28\" title=\"Ackermann, N. et al. Final CONUS results on coherent elastic neutrino-nucleus scattering at the Brokdorf reactor. Phys. Rev. Lett. 133, 251802 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR28\" id=\"ref-link-section-d46697375e1575\" target=\"_blank\" rel=\"noopener\">28<\/a>. This component was present after ventilation of the cryostats with argon gas to avoid mechanical deformation during a regular leakage test at KBR in July 2019 (ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 29\" title=\"Bonet, H. et al. Full background decomposition of the CONUS experiment. Eur. Phys. J. C 83, 195 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR29\" id=\"ref-link-section-d46697375e1579\" target=\"_blank\" rel=\"noopener\">29<\/a>). Of the four detectors used in the CONUS experiment, two (C2 and C3) are used for the analysis presented in this work and are affected by this additional background. The simulations and the background model show that an additional component with the same shape is still present in the background of these two detectors, but is absent in C5, which was not at KBR at the time. This additional background is constant during reactor on and off periods. Therefore, the leakage test component was again included in the background model of C2 and C3 by modelling it using a function with two parameters, as in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 28\" title=\"Ackermann, N. et al. Final CONUS results on coherent elastic neutrino-nucleus scattering at the Brokdorf reactor. Phys. Rev. Lett. 133, 251802 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR28\" id=\"ref-link-section-d46697375e1583\" target=\"_blank\" rel=\"noopener\">28<\/a>. The resulting impact is below 10% in the energy region between 0.4\u2009keVee and 1\u2009keVee.<\/p>\n<p>Other background components<\/p>\n<p>The remaining background in each detector is made up of many different components, similar to the situation in KBR<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 29\" title=\"Bonet, H. et al. Full background decomposition of the CONUS experiment. Eur. Phys. J. C 83, 195 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR29\" id=\"ref-link-section-d46697375e1600\" target=\"_blank\" rel=\"noopener\">29<\/a>. There are no hints that the detector upgrades or the movement of the setup were introducing any contamination. The background model includes cosmogenically induced isotopes in the copper parts (57Co, 60Co and 54Mn) of the cryostat and the Ge crystals (57Co, 68Ge, 68Ga, 65Zn and 3H), Radon inside the detector chamber, 210Pb inside the cryostat and shield, metastable Ge states (71mGe, 73mGe and 75mGe) and inert gases coming from the reactor (85Kr, 135Xe and 3H). The results of the simulation of these components were scaled to be in accordance either with the rates of gamma lines produced by the decay of these isotopes in the spectrum (for example, for radon) or with screening measurements performed before the installation. The listed contributions are typically very subdominant in the region of interest, with the decay of radon inside the detector chamber being the only exception. This contribution results in a background rate of (1.9\u2009\u00b1\u20090.1)\u2009counts\u2009day\u22121\u2009kg\u22121 between 0.4\u2009keVee and 1\u2009keVee for the C5 detector (C2: (2.8\u2009\u00b1\u20090.1)\u2009counts\u2009day\u22121\u2009kg\u22121; C3: (2.6\u2009\u00b1\u20090.1)\u2009counts\u2009day\u22121\u2009kg\u22121). These values correspond to approximately 5% of the background in this energy region. Radon decays have a high impact on energies above 100\u2009keVee, at which they can contribute up to 60% of the measured background. Slow pulses arising from decays on the surface of the diode and in the transition layer are included in the model<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 29\" title=\"Bonet, H. et al. Full background decomposition of the CONUS experiment. Eur. Phys. J. C 83, 195 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR29\" id=\"ref-link-section-d46697375e1656\" target=\"_blank\" rel=\"noopener\">29<\/a>. Energy depositions from these events can be stopped within the transition layer, and the released charge diffuses slowly into the active volume, resulting in long rise times and incomplete charge collection. Their impact is accounted for by registering the exact coordinates of an interaction in the Ge crystal. If the coordinates place it within the transition layer of the crystal, the energy of the event is shifted towards lower energies using a sigmoid-like function. Details on this procedure can be found in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 2\" title=\"Bonet, H. et al. Large-size sub-keV sensitive germanium detectors for the CONUS experiment. Eur. Phys. J. C 81, 267 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR2\" id=\"ref-link-section-d46697375e1660\" target=\"_blank\" rel=\"noopener\">2<\/a>.<\/p>\n<p>Model differences in reactor off data<\/p>\n<p>The background model accounts for the experimental differences during data collection with the reactor off. The first of these differences is induced by the fact that during a reactor outage, the drywell head from the containment structure surrounding the reactor core is placed directly above the CONUS+ room. This drywell head is made of 3.8\u2009cm steel and, therefore, increases the overburden of the experiment by approximately 0.3\u2009m w.e., which results in a reduction of 19% in the flux of cosmogenic neutrons and a reduction of 3% in the flux of cosmic ray muons. As a result, these two background contributions are reduced accordingly. The second difference in the background model of the reactor off period comes from a more effective removal of radon in the detector chamber. During the course of run 1, shield tightness and radon-free air flushing were improved. As a result, the radon contribution in reactor off time is reduced by approximately a factor of 4\u20136 compared with that for reactor on. The radon contributions are scaled to match the count rates in the gamma lines induced by the radon decay. Moreover, reactor-correlated background components, such as reactor neutrons and high-energy gammas, for example, from 16N, were investigated. Their impact was found to be negligible in all energy regions, including the region of interest. The resulting background model for all three detectors can be seen in Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig2\" target=\"_blank\" rel=\"noopener\">2<\/a> and Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig6\" target=\"_blank\" rel=\"noopener\">3<\/a>. The time-dependent contributions of the M-shell line of the 68Ge\/71Ge decays are included in the background model. The difference between the M-shell contribution in reactor on and off periods was estimated based on the measured K-shell count rates. This difference was found to be below 1% of the observed neutrino signal.<\/p>\n<p>In Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig7\" target=\"_blank\" rel=\"noopener\">4<\/a>, the reactor on and off datasets are compared with the background models. If the calculated difference between on and off phases is added to the reactor off data, good agreement with the measured on data is found above the signal region. For all detectors, background models and data are fully consistent above 0.4\u2009keVee. As the data in the reactor off phase are statistically limited, the background model still plays an important part in the analysis. Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig8\" target=\"_blank\" rel=\"noopener\">5<\/a> shows the fractional contribution of the main components to the total background rate as predicted in the model.<\/p>\n<p>Quenching<\/p>\n<p>The ratio of the ionization energy released by nuclear recoil in a CE\u03bdNS event and the ionization energy of electrons of the same energy is given by the quenching factor. In the CONUS+ analysis, the energy-dependent signal quenching is described by the Lindhard model<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 32\" title=\"Lindhard, J., Scharff, M. &amp; Schi&#xF8;tt, H. Range Concepts and Heavy Ion Ranges (Munksgaard, 1963).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR32\" id=\"ref-link-section-d46697375e1705\" target=\"_blank\" rel=\"noopener\">32<\/a> with a quenching parameter k\u2009=\u2009(0.162\u2009\u00b1\u20090.004) as determined in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 34\" title=\"Bonhomme, A. et al. Direct measurement of the ionization quenching factor of nuclear recoils in germanium in the keV energy range. Eur. Phys. J. C 82, 815 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR34\" id=\"ref-link-section-d46697375e1712\" target=\"_blank\" rel=\"noopener\">34<\/a>. Alternative quenching descriptions are also tested as described in the supplementary material of ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 28\" title=\"Ackermann, N. et al. Final CONUS results on coherent elastic neutrino-nucleus scattering at the Brokdorf reactor. Phys. Rev. Lett. 133, 251802 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR28\" id=\"ref-link-section-d46697375e1716\" target=\"_blank\" rel=\"noopener\">28<\/a>, including a linear and cubic functional form to describe the increased quenching factor compared with the Lindhard theory found in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 33\" title=\"Collar, J. I., Kavner, A. R. L. &amp; Lewis, C. M. Germanium response to sub-keV nuclear recoils: a multipronged experimental characterization. Phys. Rev. D 103, 122003 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR33\" id=\"ref-link-section-d46697375e1720\" target=\"_blank\" rel=\"noopener\">33<\/a>. It was shown that the Migdal effect is subdominant with respect to the CE\u03bdNS signal in our region of interest<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 36\" title=\"Atzori Corona, M., Cadeddu, M., Cargioli, N., Dordei, F. &amp; Giunti, C. On the impact of the Migdal effect in reactor CE&#x3BD;NS experiments. Phys. Lett. B 852, 138627 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR36\" id=\"ref-link-section-d46697375e1725\" target=\"_blank\" rel=\"noopener\">36<\/a>.<\/p>\n<p>The signal predictions for the different quenching descriptions are shown in Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig9\" target=\"_blank\" rel=\"noopener\">6<\/a>, together with the difference between the data in the reactor on phase and the background model. A signal rate of (2,600\u2009\u00b1\u2009300) events is expected for the linear function, whereas for the cubic function, (550\u2009\u00b1\u200950) events are predicted. Both numbers are significantly higher than the neutrino rate extracted from the CONUS+ data. The standard Lindhard model provides the best description of the reactor-correlated excess at low energy. In Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig9\" target=\"_blank\" rel=\"noopener\">6<\/a>, we have included an additional data point at 200\u2009eVee, which is obtained using the information shown in fig. 4 of ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 19\" title=\"Colaresi, J., Collar, J. I., Hossbach, T. W., Lewis, C. M. &amp; Yocum, K. M. Measurement of coherent elastic neutrino-nucleus scattering from reactor antineutrinos. Phys. Rev. Lett. 129, 211802 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR19\" id=\"ref-link-section-d46697375e1740\" target=\"_blank\" rel=\"noopener\">19<\/a>. There, a CE\u03bdNS signal was approximated with an exponential using two parameters, the amplitude at 200\u2009eVee (A0.2) and a decay constant. The favoured value for A0.2 is shown within a 2\u03c3 contour in a two-dimensional plot of the two parameters. We scaled this A0.2 value to our exposure and corrected for the difference in neutrino flux. As expected, it matches the description of Lindhard with a linear increase added at low energies, but it is in clear conflict with the CONUS+ data. The lower error bar corresponds to the smallest value of A0.2 in the 2\u03c3 contour. This value is also ruled out by the CONUS+ data points.<\/p>\n<p>Likelihood fit and systematic treatment<\/p>\n<p>A likelihood function is used to determine the CE\u03bdNS signal in the CONUS+ run 1 data.<\/p>\n<p>$$\\begin{array}{l}-2\\log {\\mathcal{L}}=-2\\log {{\\mathcal{L}}}_{{\\rm{ON}}}-2\\log {{\\mathcal{L}}}_{{\\rm{OFF}}}\\\\ \\,\\,\\,\\,\\,+\\sum _{ij}{({{\\boldsymbol{\\theta }}}_{i}-{\\bar{{\\boldsymbol{\\theta }}}}_{i})}^{{\\rm{T}}}\\,{{\\rm{Cov}}}_{ij}^{-1}\\,({{\\boldsymbol{\\theta }}}_{j}-{\\bar{{\\boldsymbol{\\theta }}}}_{j})\\\\ \\,\\,\\,\\,\\,+\\sum _{i}\\frac{{({\\theta }_{i}-{\\bar{\\theta }}_{i})}^{2}}{{\\sigma }_{{\\theta }_{i}}},\\end{array}$$<\/p>\n<p>\n                    (2)\n                <\/p>\n<p>where \\({{\\mathcal{L}}}_{{\\rm{O}}{\\rm{N}}}\\) and \\({{\\mathcal{L}}}_{{\\rm{O}}{\\rm{F}}{\\rm{F}}}\\) are the binned likelihood functions for reactor on and off periods. Gaussian pull terms for the systematic uncertainties are also included. Here, the first term represents the pull terms for correlated parameters, namely, the trigger efficiency parameters, and the second term represents pull terms for uncorrelated parameters, such as the active mass of the detectors, the reactor neutrino flux and the uncertainty on the energy scale calibration of the spectra. Parameters that were experimentally determined are pulled to their measured values. The background scaling factor b, an additional fit parameter for the overall normalization of the background model, is also included and pulled to 1. The binned likelihood functions assume a Poisson distribution and have the form<\/p>\n<p>$$-\\log {\\mathcal{L}}=\\mathop{\\sum }\\limits_{i=1}^{n}-{n}_{i}\\log ({\\mu }_{i})+{\\mu }_{i}+\\log ({n}_{i}!),$$<\/p>\n<p>\n                    (3)\n                <\/p>\n<p>where N is the number of bins in the region of interest, ni is the bin content in the histograms of the measured data and \u03bci is the bin content of the model. For a single detector, the model \u03bc for reactor on data is calculated from the sum of the background model nb (scaled with b) and the predicted CE\u03bdNS spectrum ns (scaled with the signal parameter s) by taking into account the live time of the experiment (tON and tOFF), the active volume of the detector (mact), the dead-time correction (cdt), the fission flux (\u03b82) and a multiplication factor (\u03b83), summarizing uncertainties of the detector response. In total,<\/p>\n<p>$$\\begin{array}{l}{\\mu }_{i}^{{\\rm{ON}}}=\\left(s\\times {t}_{{\\rm{ON}}}\\times {\\theta }_{2}\\times {n}_{i}^{s}+b\\times {n}_{i}^{b}\\times \\frac{{t}_{{\\rm{ON}}}}{{t}_{{\\rm{OFF}}}}\\right)\\times {c}_{{\\rm{dt}}}\\times \\frac{{\\theta }_{3}}{{m}_{{\\rm{act}}}},\\\\ {\\mu }_{i}^{{\\rm{OFF}}}=b\\times {n}_{i}^{b}\\times {c}_{{\\rm{dt}}}\\times \\frac{{\\theta }_{3}}{{m}_{{\\rm{act}}}}.\\end{array}$$<\/p>\n<p>\n                    (4)\n                <\/p>\n<p>The combined fit minimizes \\(-\\log {{\\mathcal{L}}}_{{\\rm{ON}}}\\) and \\(-\\log {{\\mathcal{L}}}_{{\\rm{OFF}}}\\) for all three detectors simultaneously, whereas the signal parameter s, which indicates the number of CE\u03bdNS counts, is shared among the detectors.<\/p>\n<p>Fits using just single detectors independently give consistent results as shown in Extended Data Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Tab3\" target=\"_blank\" rel=\"noopener\">2<\/a>. To cross-check the result, fits were performed by two independent likelihoods with different approaches concerning the uncertainty on the quenching model. Likelihood A applied a predicted CE\u03bdNS spectrum with a fixed k parameter (k\u2009=\u20090.162), while introducing a fourth-order polynomial to vary the shape of the signal spectrum, applying Gaussian pull terms on each parameter. In likelihood B, the k value of the Lindhard model is a fit parameter with a pull term. Moreover, there are some differences in the treatment of the non-linearity corrections and the minimization algorithms between the two likelihood fits. The results of both likelihood implementations agree within 2%.<\/p>\n<p>For the signal prediction, we choose the Helm parameterization of the nuclear form factor<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 44\" title=\"Engel, J. Nuclear form factors for the scattering of weakly interacting massive particles. Phys. Lett. B 264, 114&#x2013;119 (1991).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR44\" id=\"ref-link-section-d46697375e2867\" target=\"_blank\" rel=\"noopener\">44<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 45\" title=\"Silk, J. et al. Particle Dark Matter: Observations, Models and Searches (Cambridge Univ. Press, 2010).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR45\" id=\"ref-link-section-d46697375e2870\" target=\"_blank\" rel=\"noopener\">45<\/a> and a reactor antineutrino spectrum based on the data-driven approach in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 31\" title=\"An, F. P. et al. Antineutrino energy spectrum unfolding based on the Daya Bay measurement and its applications. Chin. Phys. C 45, 073001 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR31\" id=\"ref-link-section-d46697375e2874\" target=\"_blank\" rel=\"noopener\">31<\/a>. The given antineutrino spectra are adjusted to the fission fraction of the Leibstadt reactor and augmented with measured antineutrino spectra at energies above 8\u2009MeV (ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 46\" title=\"An, F. P. et al. First measurement of high-energy reactor antineutrinos at Daya Bay. Phys. Rev. Lett. 129, 041801 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR46\" id=\"ref-link-section-d46697375e2878\" target=\"_blank\" rel=\"noopener\">46<\/a>) and simulation data below 1.8\u2009MeV (the threshold of inverse beta decay)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 47\" title=\"Estienne, M. et al. Updated summation model: an improved agreement with the Daya Bay antineutrino fluxes. Phys. Rev. Lett. 123, 022502 (2019).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR47\" id=\"ref-link-section-d46697375e2882\" target=\"_blank\" rel=\"noopener\">47<\/a>. We assume no uncertainty on the nuclear form factor, but account for a 2% uncertainty on the weak mixing angle<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 48\" title=\"Cadeddu, M., Dordei, F., Giunti, C., Li, Y. F. &amp; Zhang, Y. Y. Neutrino, electroweak, and nuclear physics from COHERENT elastic neutrino-nucleus scattering with refined quenching factor. Phys. Rev. D 101, 033004 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR48\" id=\"ref-link-section-d46697375e2886\" target=\"_blank\" rel=\"noopener\">48<\/a>, leading to an overall uncertainty of 3.2% on the CE\u03bdNS cross-section. The shape of the applied antineutrino spectra, the thermal power of the reactor, fission fractions and energy releases per fission of the relevant fission isotopes contribute a combined 4.6% error to the expected event rate. Further uncertainties on quenching, detector active mass, trigger efficiency and the energy threshold are taken into account and listed in Extended Data Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Tab4\" target=\"_blank\" rel=\"noopener\">3<\/a>. All contributions lead to an overall uncertainty of 17% on the signal prediction. To be conservative, we assume prediction uncertainties are fully correlated between single detectors. Currently, the dominant uncertainty is from the energy calibration. Future extended calibrations will allow us to reduce this uncertainty from 5\u2009eVee to 3\u2009eVee.<\/p>\n<p>The likelihood fit itself gives a result of (395\u2009\u00b1\u200986) CE\u03bdNS counts in the combined fit of the three detectors. Additional systematic uncertainties not implemented as pull terms are evaluated in a second step and added in quadrature, giving the final uncertainty of \u00a0\u00b1106 signal counts, as shown in the Extended Data Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Tab4\" target=\"_blank\" rel=\"noopener\">3<\/a>. The non-linearity term was obtained by varying the correction parameters and checking the impact on the likelihood result. The uncertainties of the calibration points in Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig5\" target=\"_blank\" rel=\"noopener\">2<\/a> were used to generate the parameter variations. A Gaussian fit was performed on the distribution of the central values of the new likelihood results, where its 1\u03c3 value was taken as additional systematic uncertainty. The term related to the background model was studied in the same way by varying muon flux, neutron flux and the leakage component. The variations were 14% (taken from ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 40\" title=\"Goldhagen, P., Clem, J. M. &amp; Wilson, J. W. The energy spectrum of cosmic-ray induced neutrons measured on an airplane over a wide range of altitude and latitude. Radiat. Prot. Dosimetry 110, 387&#x2013;392 (2004).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR40\" id=\"ref-link-section-d46697375e2911\" target=\"_blank\" rel=\"noopener\">40<\/a>), 6% (based on ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 43\" title=\"Reyna, D. A simple parameterization of the cosmic-ray muon momentum spectra at the surface as a function of zenith angle. Preprint at &#010;                  arxiv.org\/abs\/hep-ph\/0604145&#010;                  &#010;                 (2006).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR43\" id=\"ref-link-section-d46697375e2915\" target=\"_blank\" rel=\"noopener\">43<\/a>) and 10%, respectively. Moreover, the systematic uncertainty of the fit method was estimated by the difference between likelihood A and B.<\/p>\n<p>The fractional contribution from antineutrinos of different energies to the signal expectation is quantified in Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#Fig10\" target=\"_blank\" rel=\"noopener\">7<\/a>. At the current detector threshold, we are sensitive to antineutrino energies above 5\u2009MeV. With lower detection thresholds, the steep rise in antineutrino flux towards low energies will result in markedly higher signal expectations<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 1\" title=\"Ackermann, N. et al. CONUS+ experiment. Eur. Phys. J. C 84, 1265 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-025-09322-2#ref-CR1\" id=\"ref-link-section-d46697375e2925\" target=\"_blank\" rel=\"noopener\">1<\/a>.<\/p>\n","protected":false},"excerpt":{"rendered":"Data taking Data from the CONUS+ run 1 analysis presented in this study correspond to the period from&hellip;\n","protected":false},"author":3,"featured_media":105897,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[25],"tags":[48885,10046,10047,62570,492,159,67,132,68],"class_list":{"0":"post-105896","1":"post","2":"type-post","3":"status-publish","4":"format-standard","5":"has-post-thumbnail","7":"category-physics","8":"tag-experimental-particle-physics","9":"tag-humanities-and-social-sciences","10":"tag-multidisciplinary","11":"tag-nuclear-physics","12":"tag-physics","13":"tag-science","14":"tag-united-states","15":"tag-unitedstates","16":"tag-us"},"share_on_mastodon":{"url":"https:\/\/pubeurope.com\/@us\/114944457244940337","error":""},"_links":{"self":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts\/105896","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/comments?post=105896"}],"version-history":[{"count":0,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts\/105896\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/media\/105897"}],"wp:attachment":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/media?parent=105896"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/categories?post=105896"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/tags?post=105896"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}