Total Absorption Spectroscopy Study of (92)Rb Decay: A Major Contributor to Reactor Antineutrino Spectrum Shape

Origin of the Reactor Antineutrino Anomalies in Light of a New Summation Model with Parametrized β^{-} Transitions

We investigate the possible origins of the reactor antineutrino anomalies in norm and shape within the framework of a summation model where β^{-} transitions are simulated by a phenomenological model of Gamow-Teller decay strength. The general trends of divergence from the Huber-Mueller model on the antineutrino side can be reproduced in both norm and shape. From the exact electron-antineutrino correspondence of the summation model, we predict similar distortions in the electron spectra, suggesting that biases on the reference spectra of fission electrons could be the cause of the anomalies.

First Measurement of High-Energy Reactor Antineutrinos at Daya Bay

This Letter reports the first measurement of high-energy reactor antineutrinos at Daya Bay, with nearly 9000 inverse beta decay candidates in the prompt energy region of 8-12 MeV observed over 1958 days of data collection. A multivariate analysis is used to separate 2500 signal events from background statistically. The hypothesis of no reactor antineutrinos with neutrino energy above 10 MeV is rejected with a significance of 6.2 standard deviations. A 29% antineutrino flux deficit in the prompt energy region of 8-11 MeV is observed compared to a recent model prediction. We provide the unfolded antineutrino spectrum above 7 MeV as a data-based reference for other experiments. This result provides the first direct observation of the production of antineutrinos from several high-Q_{β} isotopes in commercial reactors.

Updated Summation Model: An Improved Agreement with the Daya Bay Antineutrino Fluxes

A new summation method model of the reactor antineutrino energy spectrum is presented. It is updated with the most recent evaluated decay databases and with our total absorption gamma-ray spectroscopy measurements performed during the last decade. For the first time, the spectral measurements from the Daya Bay experiment are compared with the antineutrino energy spectrum computed with the updated summation method without any renormalization. The results exhibit a better agreement than is obtained with the Huber-Mueller model in the 2-5 MeV range, the region that dominates the detected flux. A systematic trend is found in which the antineutrino flux computed with the summation model decreases with the inclusion of more pandemonium-free data. The calculated flux obtained now lies only 1.9% above that detected in the Daya Bay experiment, a value that may be reduced with forthcoming new pandemonium-free data, leaving less room for a reactor anomaly. Eventually, the new predictions of individual antineutrino spectra for the ^{235}U, ^{239}Pu, ^{241}Pu, and ^{238}U are used to compute the dependence of the reactor antineutrino spectral shape on the fission fractions.

Fuel-Composition Dependent Reactor Antineutrino Yield at RENO

We report a fuel-dependent reactor electron antineutrino (ν[over ¯]_{e}) yield using six 2.8 GW_{th} reactors in the Hanbit nuclear power plant complex, Yonggwang, Korea. The analysis uses 850 666 ν[over ¯]_{e} candidate events with a background fraction of 2.0% acquired through inverse beta decay (IBD) interactions in the near detector for 1807.9 live days from August 2011 to February 2018. Based on multiple fuel cycles, we observe a fuel ^{235}U dependent variation of measured IBD yields with a slope of (1.51±0.23)×10^{-43}  cm^{2}/fission and measure a total average IBD yield of (5.84±0.13)×10^{-43}  cm^{2}/fission. The hypothesis of no fuel-dependent IBD yield is ruled out at 6.6σ. The observed IBD yield variation over ^{235}U isotope fraction does not show significant deviation from the Huber-Mueller (HM) prediction at 1.3 σ. The measured fuel-dependent variation determines IBD yields of (6.15±0.19)×10^{-43} and (4.18±0.26)×10^{-43}  cm^{2}/fission for two dominant fuel isotopes ^{235}U and ^{239}Pu, respectively. The measured IBD yield per ^{235}U fission shows the largest deficit relative to the HM prediction. Reevaluation of the ^{235}U IBD yield per fission may mostly solve the reactor antineutrino anomaly (RAA) while ^{239}Pu is not completely ruled out as a possible contributor to the anomaly. We also report a 2.9 σ correlation between the fractional change of the 5 MeV excess and the reactor fuel isotope fraction of ^{235}U.

Large Impact of the Decay of Niobium Isomers on the Reactor ν[over ¯]_{e} Summation Calculations

Even mass neutron-rich niobium isotopes are among the principal contributors to the reactor antineutrino energy spectrum. They are also among the most challenging to measure due to the refractory nature of niobium, and because they exhibit isomeric states lying very close in energy. The β-intensity distributions of ^{100gs,100m}Nb and ^{102gs,102m}Nb β decays have been determined using the total absorption γ-ray spectroscopy technique. The measurements were performed at the upgraded Ion Guide Isotope Separator On-Line facility at the University of Jyväskylä. Here, the double Penning trap system JYFLTRAP was employed to disentangle the β decay of the isomeric states. The new data obtained in this challenging measurement have a large impact in antineutrino summation calculations. For the first time the discrepancy between the summation model and the reactor antineutrino measurements in the region of the shape distortion has been reduced.

Summation Calculations for Reactor Antineutrino Spectra, Decay Heat and Delayed Neutron Fractions Involving New TAGS Data and Evaluated Databases

Three observables of interest for present and future reactors depend on the β decay properties of the fission products: antineutrinos from reactors, the reactor decay heat and delayed neutron emission. In these proceedings, we present new results from summation calculations of the first two quantities quoted above, performed with evolved independent yields coupled with fission product decay data, from various nuclear data bases or models. New TAGS results from the latest experiment of the TAGS collaboration at the JYFL facility of Jyväskylä will be displayed as well as their impact on the antineutrino spectra and the decay heat associated to fission pulses of the main actinides.

Analysis of the Daya Bay Reactor Antineutrino Flux Changes with Fuel Burnup

We investigate the recent Daya Bay results on the changes in the antineutrino flux and spectrum with the burnup of the reactor fuel. We find that the discrepancy between current model predictions and the Daya Bay results can be traced to the original measured ^{235}U/^{239}Pu ratio of the fission β spectra that were used as a base for the expected antineutrino fluxes. An analysis of the antineutrino spectra that is based on a summation over all fission fragment β decays, using nuclear database input, explains all of the features seen in the Daya Bay evolution data. However, this summation method still allows for an anomaly. We conclude that there is currently not enough information to use the antineutrino flux changes to rule out the possible existence of sterile neutrinos.

Recent advances in beta decay measurements

Three observables of interest for present and future reactors depend on the β decay data of the fission products: the reactor decay heat, antineutrinos from reactors and delayed neutron emission. Concerning the decay heat, significant discrepancies still exist between summation calculations in − their two main ingredients: the decay data and the fission yields − performed using the most recent evaluated databases available. It has been recently shown that the associated uncertainties are dominated by the ones on the decay data. But the results subtantially differ taking into account or not the correlations between the fission products. So far the uncertainty propagation does not include as well systematic effects on nuclear data such as the Pandemonium effect which impacts a large number of nuclei contributing to the decay heat. The list of nuclei deserving new TAGS measurements has been updated recently in the frame of IAEA working groups. The issues listed above impact in the same way the predicted energy spectra of the antineutrinos from reactors computed with the summation method, the interest of which has been recently reinforced by the Daya Bay latest publication. Nuclear data should definitely contribute to refine and better control these calculations. Lastly, a lot of nuclear data related to delayed neutrons are missing in nuclear databases. Despite the progresses already done these last years with new measurements now requiring to be included in evaluated databases, the experimental efforts which still need to be done are significant. These different issues will be addressed here before to comment on recent experimental results and on their impacts on the quoted observables. Some perspectives will also be presented. Solving the issues listed above will require to bring together experimental, simulation, evaluation and theoretical activities.

Dissecting Reactor Antineutrino Flux Calculations

Current predictions for the antineutrino yield and spectra from a nuclear reactor rely on the experimental electron spectra from ^{235}U, ^{239}Pu, ^{241}Pu and a numerical method to convert these aggregate electron spectra into their corresponding antineutrino ones. In the present work we investigate quantitatively some of the basic assumptions and approximations used in the conversion method, studying first the compatibility between two recent approaches for calculating electron and antineutrino spectra. We then explore different possibilities for the disagreement between the measured Daya Bay and the Huber-Mueller antineutrino spectra, including the ^{238}U contribution as well as the effective charge and the allowed shape assumption used in the conversion method. We observe that including a shape correction of about +6%  MeV^{-1} in conversion calculations can better describe the Daya Bay spectrum. Because of a lack of experimental data, this correction cannot be ruled out, concluding that in order to confirm the existence of the reactor neutrino anomaly, or even quantify it, precisely measured electron spectra for about 50 relevant fission products are needed. With the advent of new rare ion facilities, the measurement of shape factors for these nuclides, for many of which precise beta intensity data from TAGS experiments already exist, would be highly desirable.

Impact of Modular Total Absorption Spectrometer measurements of β decay of fission products on the decay heat and reactor ν[over ¯]_{e} flux calculation

We report the results of a β-decay study of fission products ^{86}Br, ^{89}Kr, ^{89}Rb, ^{90gs}Rb, ^{90m}Rb, ^{90}Kr, ^{92}Rb, ^{139}Xe, and ^{142}Cs performed with the Modular Total Absorption Spectrometer (MTAS) and on-line mass-separated ion beams. These radioactivities were assessed by the Nuclear Energy Agency as having high priority for decay heat analysis during a nuclear fuel cycle. We observe a substantial increase in β feeding to high excited states in all daughter isotopes in comparison to earlier data. This increases the average γ-ray energy emitted by the decay of fission fragments during the first 10 000 s after fission of ^{235}U and ^{239}Pu by approximately 2% and 1%, respectively, improving agreement between results of calculations and direct observations. New MTAS results reduce the reference reactor ν[over ¯]_{e} flux used to analyze reactor ν[over ¯]_{e} interaction with detector matter. The reduction determined by the ab initio method for the four nuclear fuel components, ^{235}U, ^{238}U, ^{239}Pu, and ^{241}Pu, amounts to 0.976, 0.986, 0.983, and 0.984, respectively.

NEOS Data and the Origin of the 5 MeV Bump in the Reactor Antineutrino Spectrum

We perform a combined analysis of recent NEOS and Daya Bay data on the reactor antineutrino spectrum. This analysis includes approximately 1.5 million antineutrino events, which is the largest neutrino event sample analyzed to date. We use a double ratio which cancels flux model dependence and related uncertainties as well as the effects of the detector response model. We find at 3-4 standard deviation significance level, that plutonium-239 and plutonium-241 are disfavored as the single source for the so-called 5 MeV bump. This analysis method has general applicability and, in particular, with higher statistics data sets, will be able to shed significant light on the issue of the bump. With some caveats, this should also allow us to improve the sensitivity for sterile neutrino searches in NEOS.

Decays of the Three Top Contributors to the Reactor ν[over ¯]_{e} High-Energy Spectrum, ^{92}Rb, ^{96gs}Y, and ^{142}Cs, Studied with Total Absorption Spectroscopy

We report total absorption spectroscopy measurements of ^{92}Rb, ^{96gs}Y, and ^{142}Cs β decays, which are the most important contributors to the high energy ν[over ¯]_{e} spectral shape in nuclear reactors. These three β decays contribute 43% of the ν[over ¯]_{e} flux near 5.5 MeV emitted by nuclear reactors. This ν[over ¯]_{e} energy is particularly interesting due to spectral features recently observed in several experiments including the Daya Bay, Double Chooz, and RENO Collaborations. Measurements were conducted at Oak Ridge National Laboratory by means of proton-induced fission of ^{238}U with on-line mass separation of fission fragments and the Modular Total Absorption Spectrometer. We observe a β-decay pattern that is similar to recent measurements of ^{92}Rb, with a ground-state to ground-state β feeding of 91(3)%. We verify the ^{96gs}Y ground-state to ground-state β feeding of 95.5(20)%. Our measurements substantially modify the β-decay feedings of ^{142}Cs, reducing the β feeding to ^{142}Ba states below 2 MeV by 32% when compared with the latest evaluations. Our results increase the discrepancy between the observed and the expected reactor ν[over ¯]_{e} flux between 5 and 7 MeV, the maximum excess increases from ∼10% to ∼12%.

Effects of Fission Yield Data in the Calculation of Antineutrino Spectra for ^{235}U(n,fission) at Thermal and Fast Neutron Energies

Fission yields form an integral part of the prediction of antineutrino spectra generated by nuclear reactors, but little attention has been paid to the quality and reliability of the data used in current calculations. Following a critical review of the thermal and fast ENDF/B-VII.1 ^{235}U fission yields, deficiencies are identified and improved yields are obtained, based on corrections of erroneous yields, consistency between decay and fission yield data, and updated isomeric ratios. These corrected yields are used to calculate antineutrino spectra using the summation method. An anomalous value for the thermal fission yield of ^{86}Ge generates an excess of antineutrinos at 5-7 MeV, a feature which is no longer present when the corrected yields are used. Thermal spectra calculated with two distinct fission yield libraries (corrected ENDF/B and JEFF) differ by up to 6% in the 0-7 MeV energy window, allowing for a basic estimate of the uncertainty involved in the fission yield component of summation calculations. Finally, the fast neutron antineutrino spectrum is calculated, which at the moment can only be obtained with the summation method and may be relevant for short baseline reactor experiments using highly enriched uranium fuel.

Measurement of the Reactor Antineutrino Flux and Spectrum at Daya Bay

This Letter reports a measurement of the flux and energy spectrum of electron antineutrinos from six 2.9 GWth nuclear reactors with six detectors deployed in two near (effective baselines 512 and 561 m) and one far (1579 m) underground experimental halls in the Daya Bay experiment. Using 217 days of data, 296 721 and 41 589 inverse β decay (IBD) candidates were detected in the near and far halls, respectively. The measured IBD yield is (1.55±0.04) ×10(-18)  cm(2) GW(-1) day(-1) or (5.92±0.14) ×10(-43)  cm(2) fission(-1). This flux measurement is consistent with previous short-baseline reactor antineutrino experiments and is 0.946±0.022 (0.991±0.023) relative to the flux predicted with the Huber-Mueller (ILL-Vogel) fissile antineutrino model. The measured IBD positron energy spectrum deviates from both spectral predictions by more than 2σ over the full energy range with a local significance of up to ∼4σ between 4-6 MeV. A reactor antineutrino spectrum of IBD reactions is extracted from the measured positron energy spectrum for model-independent predictions.