Pandemonium effect, nuclear databases and TAS collaboration

a) Pandemonium effect and bias in nuclear decay databases

As explained in the section Reactor Neutrino Applicatio/ns, Residual Power and Summation Calculations, the β- decay properties of fission products are important for simulating the spectrum of antineutrinos emitted by reactors and for calculating the decay heat of shutdown reactors, but the study of these quantities is subject to the accuracy of existing databases which are affected by the so-called "Pandemonium" effect [1]. This effect is due to the fact that the most widely used technique for estimating beta decay properties uses Germanium (Ge) detectors which have an efficiency that falls exponentially with gamma energy, which can lead to an underestimation of the beta decay branches towards the higher energy excited levels in the daughter nucleus. The problem is shown in figure 2: if the Ge detector does not see the high energy gamma called γ1 the feeding that was supposed to go to the higher energy level will be attributed to a lower energy level.

Figure 2 (from A. Algora) : (left) Illustration of the Pandemonium effect on the energy levels fed into the daughter nuclei. (right) Difference in measurement of a de-excitation cascade between a germanium detector and an TAS.

By way of illustration, figure 3 shows an energy spectrum of the antineutrinos emitted during the beta decay of 105Mo (calculations performed using the summation method developed by the group).

Figure3 : Energy spectrum of the antineutrinos emitted during the beta decay of 105Mo. Comparison between the JEFF3.1 database and the data from the TAS experiment.

The dashed curve shows the spectrum calculated from the decay data present in the JEFF3.1 nuclear database from measurements made with a high-resolution Ge detector. The solid line curve shows the same calculation based on a calorimetric measurement described below performed by the group that allows correction for the Pandemonium effect. At the scale of a given isotope and therefore of a fission product, it can be seen that the mean value of the distribution is almost halved as a result of the correction and the shape of the distributions is also greatly modified. This has a direct impact on nuclear databases, summation calculations for neutrinos and decay heat and beta forces for fundamental nuclear physics.

b) Total Absorption Spectroscopy (TAS) and collaboration

The solution to this problem is the Total Absorption Spectroscopy (TAS) method, which is based on the detection of de-excitation gammas through the use of a high efficiency detector covering a solid angle of almost 4π all around the selected isotope. This method allows direct access to the energy of the excited levels populated by the beta decay in the daughter nucleus. This technique is complementary to the use of Germanium (Ge) detectors for nucleus spectroscopy since it provides high detection efficiency but with lower resolution.

An experimental project started in January 2009, with the participation of members of the SEN group and a team from Valencia (J.-L. Tain and A. Algora) and Surrey (William Gelletly) in a meeting organised by the IAEA in Vienna on the possibility of improving the calculation of the decay heat by re-measuring with the TAS (Total Absorption Spectroscopy) method some selected nuclei. This meeting was the starting point of the collaboration between the two groups to identify and measure fission products of interest for decay heat, antineutrinos and, in most cases, both.Within the framework of this project, in November 2009, part of the group participated in beta decay measurements of interesting nuclei selected by the Valencia research team and the Subatech team. These measurements were carried out at the accelerator in Jyväskylä (Finland) using the TAS method with a segmented detector composed of 12 BaF2 crystals called Rocinantes (figure 4 - left).

rocinantes dtas

Figure 4 : The Rocinantes detector (left). The DTAS detector (right).

A second experiment was carried out in February 2014 again at the Jyväskylä accelerator. During this experiment, 23 nuclei of interest for the calculation of the antineutrino spectrum and the decay heat of reactors were measured [2]. The experiment was performed with a new segmented detector, the DTAS, composed of 18 NaI crystals (Figure 4 - right) [3].

Before 2009, our Valencian colleagues had developed a TAS detector consisting of a single large NaI crystal (Lucrezia) installed at the ISOLDE accelerator at CERN. The advantage of having then designed segmented detectors opened the possibility of studying the multiplicity of detected gammas of interest in particular for the nuclear structure. The SEN group is thus interested in studying the presence of possible pygmy resonances, the signature of a neutron skin.

The analysis of the TAS data, gives us direct access to the beta feeding to each level of the daughter nucleus by solving the inverse problem d=R*f, where d represents the data, R is the response matrix of the detector and f is the beta feeding that we are trying to measure. To do this we need clean data without contamination and background noise; an accurate knowledge of our measurement system obtained through simulation and minimal knowledge of the properties of the nucleus to calculate the matrix R. For this reason, the TAS method is complementary to the already existing measurements obtained with Ge detectors, since it builds on and complements the knowledge already gained from these measurements.

Through the feeding, the TAS technique gives access to the beta intensity and, therefore, to the beta force which is a microscopic observable that theoretical models can calculate, thus opening a window on the structure of the nucleus. TAS measurements will provide new constraints useful for improving the predictability of models, which are indispensable in other fields of physics such as nuclear astrophysics in nucleosynthesis calculations.

The TAS collaboration has published many results from its two campaigns of experiments in 2009 and 2014, as an illustration, see the articles [4,6]. And the SEN team is a leader in the study of the impact of these new measurements on the calculation of the antineutrino spectrum [7,8].

Recently the SEN team has applied for a grant from the Pays de la Loire Region and made a request to ANR to improve the existing DTAS with one or more crowns of LaBr3, this is the (NA)2STARS project. These new detectors, which combine high efficiency and good energy resolution, will make it possible to measure more exotic and poorly known nuclei, inaccessible to current TAS detectors.

3 - Applications to Nuclear Structure and Astrophysics

Half of the elements with mass greater than 70 is created by the astrophysical r process, which proceeds via unstable very neutron-rich nuclei in stellar explosions or other violent astrophysical events such as the coalescence of neutron stars. The identification of the location of the r-process remains one of the major challenges of nuclear astrophysics. Recent advances in the description of the interaction of neutrinos with matter and its implementation in the modelling of supernovae explosions tend to show that supernovae explosions only contribute to the production of Z elements below 50. Compact star coalescences are currently considered to be the best candidates for determining the site of the main r-process. Observational confirmation of the occurrence of the r-process in such an astrophysical event was given at the end of 2017 with the multi-messenger detection of the merging of two neutron stars and the measurement of the emitted electromagnetic radiation. Beta decay of nuclei is among the important processes that influence r-element abundances along with neutron capture reactions, photo-dissociation, temperature and density. In particular, measurements of radioactive periods of the progenitors of stable nuclei help to determine their abundance and are thus important ingredients in nucleosynthesis calculations that attempt to reproduce them. With the exception of a few key nuclei that are relatively close to the valley of stability, half-lives must be calculated using theoretical models. In this type of calculation, the strength distribution associated with beta decay as a function of energy must be calculated for all possible end states. This involves determining what proportion of the strength resides in the energy window opened by the beta decay. The TAGS technique is the preferred technique to perform these measurements, and can also be used to study the presence of low-energy collective modes that influence the paths of the nucleosynthesis process r. The TAGS measurements proposed by the Subatech team, in international collaboration with the IFIC in Valencia, Spain, CIEMAT in Madrid, the University of Surrey in England (and other international laboratories that are coming to participate in the experiments), at state-of-the-art facilities for the production of exotic nuclei such as JYFL in Finland, or ISOLDE at CERN (among others) will provide the necessary tests for the theoretical models.

4 - Applications to Reactor Neutrinos, Decay Heat and Summation Calculations

In a PWR-type reactor, the total number of fissions that occur as a function of time comes essentially from the fission of the 4 isotopes 235U, 238U initially present in the fuel and 239Pu and 241Pu produced by neutron capture and successive beta decay of 238U. Figure 5, which shows this evolution as a function of time (the reactor burn-up), shows that 235U is the main contributor to this fission rate. The increase in the production of 239Pu over time also strongly increases the associated fission rate. The resulting fission products, which are neutron-rich nuclei, are then de-excited by beta or beta-n decay (with the emission of a beta-delayed neutron) and are therefore the source of the considerable flux of antineutrinos that comes out of nuclear reactors.

Figure 5 : Evolution of the fission rate per second in a PWR-type reactor.

The beta decay properties of neutron-rich nuclei, in particular of fission products play a major role in the estimation of observables important for nuclear reactor safety as well as neutrino physics, whether fundamental or applied. Recently, the reactor anomaly - a significant discrepancy between predictions made from recent estimates of antineutrino energy spectra and fluxes measured by neutrino experiments within 100m of a reactor core - has attracted the attention of the particle physics and nuclear physics communities [9]. This issue has led to many new neutrino physics experiments in the vicinity of research reactors around the world including the Solid experiment in which the group is involved.

A better knowledge of the spectrum of reactor antineutrinos is required, and for this, a mastery of the underlying nuclear physics is essential. Nuclear physics plays a major role in the two spectrum calculation methods currently under study and the Subatech team is an international leader in the prediction of antineutrino spectra based on the summation method from nuclear data [7,8] in a highly competitive context. For a given isotope, the summation method makes it possible to calculate the expected total energy spectrum of antineutrinos produced in a nuclear reactor. Based on the knowledge of the decay data of individual nuclei, one can first calculate the expected individual antineutrino spectra for the approximately 800 fission products of a given isotope using Fermi theory. The summation method then consists of summing all these individual energy spectra weighted by their fission yield from the corresponding nuclear databases.

Our most recent predictions, of unequalled quality thanks to Total Absorption Gamma-ray Spectroscopy (TAGS) measurements carried out over the last decade, have been published in a high-impact journal [8], a viewpoint on a PRL of the Daya Bay collaboration on the topic [10] and in an article in the CNRS/in2p3 newsletter ("Deficit of antineutrinos produced in reactors: nuclear physics provides elements of an answer", with Muriel Fallot's team). These predictions have recently been taken up again in a new publication [11] to study their global impact on neutrino physics.

The aim of our TAGS experiments and instrumental projects (e-Shape, (NA)2STARS) is to provide a definitive answer to the questions associated with the emission of antineutrinos from reactors. Finally, a common issue with the study of antineutrinos is the study of fission products at the origin of the decay heat of the reactors after their shutdown. This decay heat represents about 7%!o(MISSING)f the nominal power of a reactor and must be evacuated or else the fuel cladding will melt (cf. the Fukushima accident in 2011). To date, significant differences remain between integral measurements made in experiments within reactors and calculations using nuclear data, even for the decay heat associated with the most common fissile isotopes such as 235U and 239Pu. Our TAGS experiments are aimed to improve significantly this situation in order to reduce the actual uncertainties and obtain reliable decay heat predictions for innovative reactors.

Bibliographie

[1] J. C. Hardy, L. C. Carraz, B. Jonson, and P. G. Hansen, Phyics Letter B 71 (1977) 307.

[2] V. Guadilla et al., First experiment with the NUSTAR/FAIR Decay Total Absorption γγ-Ray Spectrometer (DTAS) at the IGISOL IV facility, Nucl. Inst. and Meth. B 376 (2016) 334.

[3] V. Guadilla et al., Characterization and performance of the DTAS detector, Nucl. Instrum. Meth. A 910 (2018) 79-89.

[4] A.-A. Zakari-Issofou et al., Total Absorption Spectroscopy Study of 92Rb Decay: A Major Contributor to Reactor Antineutrino Spectrum Shape, Phys. Rev. Lett. 115, 102503 (2015);

[5] V. Guadilla et al., Total absorption γ-ray spectroscopy of the β-delayed neutron emitters 137I and 95Rb, Phys. Rev. C 100 (4) (2019) 044305.

[6] V. Guadilla et al., Large Impact of the Decay of Niobium Isomers on the Reactor ¯νe Summation Calculations, Phys. Rev. Lett. 122 (2019) 042502.

[7] M. Fallot et al., New Antineutrino Energy Spectra Predictions from the Summation of Beta Decay Branches of the Fission Products, Phys. Rev. Lett. 109, 202504 (2012).

[8] M. Estienne et al., Updated Summation Model: An Improved Agreement with the Daya Bay Antineutrino Fluxes, Phys. Rev. Lett. 123, (2019) 022502.

[9] G. Mention et al., The Reactor Antineutrino Anomaly, Phys. Rev. D 83 (2011) 073006.

[10] M. Fallot, Viewpoint: Getting to the Bottom of an Antineutrino Anomaly, Physics 10, 66 (2017).

[11] J. Berryman and P. Huber, Reevaluating Reactor Antineutrino Anomalies with Updated Flux Predictions, Phys. Rev. D 101 (2020) 015008, arXiv:1909.09267.

1 - Effet Pandémonium et biais des bases de données nucléaires de désintégration

Comme expliqué dans la section Applications aux neutrinos des réacteurs, puissance résiduelle et calculs de sommation, les propriétés de décroissance β- des produits de fission sont importantes pour la simulation du spectre d’antineutrinos émis par les réacteurs et pour le calcul de la puissance résiduelle des réacteurs à l'arrêt, mais l’étude de ces quantités est sujet à la précision des bases des données existantes qui sont affectées par ce qui est appelé effet « Pandémonium » [1] . Cet effet est dû au fait que la technique la plus employée pour l’estimation de propriétés de décroissance bêta utilise des détecteurs au Germanium (Ge) qui ont une efficacité qui chute exponentiellement avec l’énergie des gammas, ce qui peut conduire à une sous-estimation des branches de décroissance bêta vers les niveaux excités de plus haute énergie dans le noyau fil. Le problème est représenté sur la figure 1 : si le détecteur Ge ne voit pas le gamma de haute énergie appelé γ1 l’alimentation qui était censée aller au niveau de plus haute énergie sera attribuée à un niveau d’énergie plus basse.

Figure 1 (de A. Algora) : (gauche) Illustration de l'effet Pandémonium sur *les niveaux d'énergie alimentés dans le noyaux fils. (droite) différence de *mesure d'une cascade de désexcitation entre un détecteur au germanium et un *TAS.

En guise d'illustration, la figure 2 présente un spectre en énergie des antineutrinos émis lors de la désintégration bêta du 105Mo (calculs réalisés avec la méthode de sommation développée par le groupe).

Figure 2 : Spectre en énergie des antineutrinos émis lors de la *désintégration bêta du 105Mo. Comparaison entre la base de données JEFF3.1 et *les données de l'expérience TAS.

La courbe en pointillés présente le spectre calculé à partir des données de désintégration présentent dans la base de données nucléaires JEFF3.1 provenant de mesures réalisées avec un détecteur Ge de haute résolution. La courbe en trait plein montre le même calcul effectué sur la base d'une mesure calorimétrique décrite ci-dessous réalisée dans le groupe qui permet de corriger de l'effet Pandémonium. A l'échelle d'un isotope donné et de ce fait d'un produit de fission, on constate que la valeur moyenne de la distribution est quasiment divisée par deux suite à la correction et la forme des distributions est fortement modifiée également. Ceci impacte directement les bases de données nucléaires, les calculs de sommation pour les neutrinos et la puissance résiduelle et les forces bêta pour la physique nucléaire fondamentale.

2 - Spectroscopie par Absorption Totale (TAS) et collaboration

La solution à ce problème est la méthode de spectroscopie par absorption totale (TAS, Total Absorption Spectroscopy) qui est basée sur la détection des gammas de désexcitation à travers l'utilisation d’un détecteur à haute efficacité couvrant un angle solide de presque 4π tout autour de l'isotope sélectionné. Cette méthode permet d'avoir accès directement à l’énergie des niveaux excités peuplés par la désintégration bêta dans le noyau fils. Cette technique est complémentaire à l’utilisation de détecteurs Germanium (Ge) pour la spectroscopie des noyaux puisqu'elle assure une haute efficacité de détection mais avec une moindre résolution.

Un projet experimental a démarré en Janvier 2009, avec la participation de membres du groupe SEN et une équipe de Valencia (J.-L. Tain et A. Algora) et Surrey (William Gelletly) à une réunion organisée par l'AIEA à Vienne sur la possibilité d’améliorer le calcul de la puissance résiduelle en re-mesurant avec la méthode TAS certains noyaux sélectionnés. Cette réunion a été le point départ d'une collaboration entre les deux groupes de Valencia et Nantes pour identifier et mesurer les produits de fission intéressants pour la puissance résiduelle, pour les antineutrinos et, dans la plupart des cas, pour les deux.

Dans le cadre de ce projet, en novembre 2009, une partie du groupe a participé à des mesures de décroissance bêta de noyaux intéressants sélectionnés par l'équipe de chercheurs de Valencia et l’équipe de Subatech. Ces mesures ont été effectuées auprès de l’accélérateur de Jyväskylä (Finlande) en utilisant la méthode TAS avec un détecteur segmenté composé par 12 cristaux de BaF2 appelé Rocinantes (figure 3 - gauche).

rocinantes dtas

Figure 3 : Le détecteur Rocinantes (gauche). Le détecteur DTAS (droite).

Une deuxième expérience a été effectuée en février 2014 toujours auprès de l’accélérateur de Jyväskylä. Pendant cette expérience, 23 noyaux d’intérêt pour le calcul du spectre d’antineutrino et la puissance résiduelle de réacteurs ont été mesurés [2]. Cette expérience a été réalisée avec un nouveau détecteur segmenté, le DTAS, composé par 18 cristaux de NaI (Figure 3 - droite) [3].

Avant 2009, nos collègues Valenciens avaient développé un détecteur TAS constitué d'un seul gros cristal de NaI (Lucrezia) installé auprès de l'accélérateur de ISOLDE au CERN. L’avantage d'avoir conçu ensuite des détecteurs segmentés a ouvert la possibilité d’étudier la multiplicité des gammas détectés d'intérêt notamment pour la structure nucléaire. Le groupe SEN s'intéresse ainsi à l'étude de la présence d’éventuelles résonances pygmées, signature d’une peau de neutrons.

L’analyse des données TAS, nous donne directement accès à l’alimentation bêta vers chaque niveau du noyau fils en résolvant le problème inverse d=R*f, où d représente les données, R est la matrice de réponse du détecteur et f l’alimentation (feeding) que nous cherchons à mesurer. Pour cela, il nous faut des données propres sans contamination et bruit de fond ; une connaissance précise de notre système de mesure obtenue par la simulation et des connaissances minimales des propriétés du noyau pour calculer la matrice R. Pour cette raison, la méthode TAS est complémentaire aux mesures déjà existantes obtenues avec des détecteurs Ge, puisqu’elle s‘appuie sur les connaissances déjà apportées par ces mesures tout en les complétant.

A travers l’alimentation, la technique TAS, donne accès à l’intensité bêta et, donc, à la force bêta qui est une observable microscopique que les modèles théoriques peuvent calculer, ouvrant ainsi une fenêtre sur la structure du noyau. Les mesures TAS apporteront de nouvelles contraintes utiles pour améliorer la prédictivité des modèles, indispensables dans d'autres domaines de la physique comme par ex. l'astrophysique nucléaire dans les calculs de nucléosynthèse.

La collaboration TAS a publié de nombreux résultats de ses deux campagnes d’expériences de 2009 et 2014, en guise d'illustration, on pourra consulter les articles [4,6]. Et l’équipe SEN est leader dans l’étude de l’impact de ces nouvelles mesures sur le calcul du spectre antineutrino [7,8].

Dernièrement L’équipe SEN a demandé une subvention à la Région Pays de la Loire et fait une demande d'ANR pour améliorer le DTAS existant avec une ou plusieurs couronnes de LaBr3, c'est le projet (NA2)STARS. Ces nouveaux détecteurs qui rassemblent une haute efficacité et une bonne résolution en énergie permettront d’aller mesurer des noyaux plus exotiques et mal connus, inaccessibles aux détecteurs TAS actuels.