Beta-decay studies for applied and basic nuclear physics

In this review we will present the results of recent beta-decay studies using the total absorption technique that cover topics of interest for applications, nuclear structure and astrophysics. The decays studied were selected primarily because they have a large impact on the prediction of a) the decay heat in reactors, important for the safety of present and future reactors and b) the reactor electron antineutrino spectrum, of interest for particle/nuclear physics and reactor monitoring. For these studies the total absorption technique was chosen, since it is the only method that allows one to obtain beta decay probabilities free from a systematic error called the Pandemonium effect. The measurements presented and discussed here were performed mainly at the IGISOL facility of the University of Jyvaskyla (Finland) using isotopically pure beams provided by the JYFLTRAP Penning trap. Examples are presented to show that the results of our measurements on selected nuclei have had a large impact on predictions of both the decay heat and the anti-neutrino spectrum from reactors. Some of the cases involve beta delayed neutron emission thus one can study the competition between gamma- and neutron-emission from states above the neutron separation energy. The gamma-to-neutron emission ratios can be used to constrain neutron capture (n,gamma)cross sections for unstable nuclei of interest in astrophysics. The information obtained from the measurements can also be used to test nuclear model predictions of half-lives and Pn values for decays of interest in astrophysical network calculations. These comparisons also provide insights into aspects of nuclear structure in particular regions of the nuclear chart.


Introduction
Our knowledge of the properties of atomic nuclei is derived almost entirely from studies of nuclear reactions and radioactive decays. The ground and excited states of nuclei exhibit many forms of decay but the most common are alpha, beta and gamma-ray emission. Our focus here is on beta decay in its various manifestations. A glance at the Segre Chart reveals that it is the most common way for the ground states of nuclei to decay and it is frequently the observation of such beta decays that brings us our first knowledge of a particular nuclear species and its properties.
The study of beta decay is intrinsically much more difficult than the study of either alpha or gamma decay. The reason for this is straightforward. Alpha particles and gamma rays are emitted with discrete energies determined by the differences in energy between the initial and final states involved. Thus characteristic alpha and gamma ray spectra exhibit a series of discrete lines. It requires sophisticated detection and analysis techniques to determine the excitation energies of the states involved, their lifetimes and the transition rates between states. Beta decay carries the same information, but the difficulties of measurement and interpretation are compounded because the spectrum is continuous, not discrete. In 1930 this was explained by Pauli's hypothesis [1] of the existence of a neutral, zero mass particle called in his letter the neutron that is emitted with the beta particle. The sharing of momentum and energy then explains the continuous spectrum. Shortly afterwards Fermi [2] was able to formulate a theory of beta decay based on this idea and coined the name neutrino (little neutral one) for the particle.
A knowledge of beta decay transition probabilities is of particular importance for application to a) tests of nuclear model calculations, b) the radioactive decay heat in Fig. 1. Schematic picture of how the beta feeding is determined in a beta decay experiment employing Ge detectors. The beta feeding (I β (i)) to level i is determined from the difference of the total intensity feeding the level and those deexciting it. The sum over (k) represents all transitions feeding or de-exciting the level. Iγ k stands for the gamma intensity of transition k and ICE k represents the conversion electron intensity.
reactors, c) the reactor electron anti-neutrino spectrum and d) reaction network calculations for nucleosynthesis in explosive stellar events. In this article we will provide examples of our recent studies of beta decays that involve the use of total absorption gamma spectroscopy (TAGS) to tackle the topics listed above. The TAGS method was adopted in our measurements because it overcomes the difficulties inherent in the conventional use of Ge detector arrays for this purpose. Such arrays are an important and essential tool for constructing nuclear decay schemes since they are very well suited to the study of gamma-gamma coincidences, the main basis for building such schemes. The normal practice is then to derive beta decay transition probabilities for each level populated from the difference in the total intensity of all the gamma rays feeding the level and the sum of the intensities of all those deexciting it, corrected by the effect of internal conversion (see Figure 1). In principle this allows us to obtain the beta branching to every level, assuming that we are able to determine by some other means the number of decays that go directly to the daughter ground state, which are not accompanied by gamma emission.
Unfortunately this "simple" procedure does not necessarily give us the correct answers. States at high excitation energies in the daughter nucleus can be populated if the Q β value of the decay is large. In this case both the number of levels that can be directly populated by the beta decay is large and the number of levels available to which they can gamma decay is also large. As a result, in general, individual gamma rays (emitted by levels at high excitation energy) have low intensity. Ge detectors, indeed even gamma-ray arrays, have limited detection efficiencies particularly at higher energies and thus weak transitions are often not detected in experiments. It is clear that this means that we have a problem that has become known as the Pandemonium effect [3] (see Figure 2 for a simplified picture).
We can overcome this problem using the total absorption gamma spectroscopy technique, where we take a different approach. The method involves a large 4π scintillation detector and is based on the detection of the full deexcitation gamma cascade for each populated level, rather than the individual gamma rays. The power of TAGS to find the missing beta intensity has been demonstrated in a number of papers [4,5,6,7,8,9,10,11,12,13]. The use of the TAGS method began at ISOLDE [14]. Its development and history are described in [15,16].
Looking at a wider picture we see that many entries in the international databases, that rely on measurements with Ge detectors alone, will have systematic errors. As we shall see in the sections that follow this means that the results cannot be relied on for certain applications. The answer to the resulting difficulties lies in the use of TAGS. In the remainder of this article we will describe the TAGS method in more detail and then use our results to illustrate how it can be applied.
The structure of this article is the following: in Section 2 details of the experimental method and the analysis of the spectra are described. Sections 3, 4, 5 and 6 deal with beta decay studies related to a) radioactive decay heat (DH), b) reactor antineutrino spectra c) nuclear models, and d) astrophysical applications respectively. Finally, in Section 7, a summary will be presented.

TAGS measurements
In Section 1 it was already explained why we need TAGS measurements. Figure 3 shows how the simple beta decay presented in Figure 2 is detected by typical detectors Fig. 3. Schematic picture of how the simple beta decay depicted in Figure 2 is seen ideally by different detectors used in beta decay studies. Left panel, representation of a total absorption detector, rigth panel, ideally detected spectra with a beta dectector (a silicon detector), a Ge detector and a total absorption detector after the simple decay represented in Figure 2. used in beta decay experiments. Because a TAGS detector acts like a calorimeter, in an ideal TAGS experiment the detected spectrum will be proportional to the beta intensity distribution. This spectrum is obtained in ideal conditions, where there is no penetration of the beta particles, or the radiation generated by them, into the detector and the TAGS detector is 100% efficient up to the full energy of the gamma rays that follow the beta decay. That means that only the full absorption peak corresponding to the sum energy of the gamma cascade is detected in the case of a β − decay.
A real experiment does not quite match this ideal. In order to achieve very high detection efficiencies, large, close to 4π, detector volumes are needed. Thus inorganic scintillation material has been the natural choice. Because of its good average properties NaI(Tl) has been used in all except one (see later) of the existing spectrometers. Nevertheless we need some opening to take the sources to the centre of the spectrometer, either in the form of a radioactive beam or deposited onto a tape transport system. The latter may also be needed to remove the sources after some measuring time. We may also need ancillary detectors for detecting coincidences and selecting the events in which we are interested. In addition TAGS detectors require, in general, some form of encapsulation. All these requirements mean that we have dead material and holes in our detector system. Accordingly the gamma detection efficiency of our system will not be 100%. The consequence is that to obtain the beta intensity distribution we need to solve the inverse problem represented by the following equation: where d i is the content of bin i in the measured TAGS spectrum, R ij is the response matrix of the TAGS setup and represents the probability that a decay that feeds level j in the level scheme of the daughter nucleus gives a count in bin i of the TAGS spectrum, f j is the beta feeding to the level j (our goal) and C i is the contribution of the contaminants to bin i of the TAGS spectrum. The index j in the sum runs over the levels populated in the daugther nucleus in the beta decay. The response matrix R ij depends on the TAGS setup and on the assumed level scheme of the daughter nucleus. The dependence on the level scheme of the daughter nucleus is introduced through the branching ratio matrix B. This matrix contains the information of how the different levels in the assumed level scheme decay to the lower lying levels. To calculate the response matrix R ij (B) the branching ratio matrix B has to be determined first. There are different ways to extract the feeding distribution from equation 1 or, in other words, to solve the TAGS inverse problem. One can assume the existence of "pseudo" levels that are added manually (with their decaying branches) to the known level scheme, calculate their response and see their effect in the calculated spectrum (see for example [17,18]). In our analysis until now we have followed an alternative way for which the level scheme of the daughter nucleus is divided into two regions, a low excitation part and a high excitation part. Conventionally the levels of the low excitation part and their gamma decay branchings are taken from high resolution measurements available in the literature, since it is assumed that the gamma branching ratios of these levels are well determined. Above a certain energy, the cut-off energy, a continuum of possible levels divided into 40 keV bins is assumed. From this energy up to the decay Q β value, the statistical model is used to generate a branching ratio matrix for the high excitation part of the level scheme. The statistical model is based on a level density function and gamma strength functions of E1, M1, and E2 character. Once the branching ratio matrix (B) is defined, the response of the setup R ij to that branching matrix B (or level scheme) is calculated using previously validated Monte Carlo simulations of the relevant electromagnetic interactions in the experimental setup. The validation of the Monte Carlo simulations is performed by reproducing measurements of well known radioactive sources, that are made under the same experimental conditions as the real experiment. The Monte Carlo simulations require a careful implementation of all the details of the geometry of the setup, a proper knowledge of the materials employed in the construction of the setup and testing to find the best Monte Carlo tracking options and physics models that reproduce the measured sources. It should be noted that from high resolution measurements we use only the branching ratios of the levels, and not the information on the feeding of these levels.
Once the response function is determined we can solve Equation 1 using appropriate algorithms to determine the feeding (or beta intensity) distribution. In our analyses we follow the procedure developed by the Valencia group. In [19] several algorithms were explored. From those that are possible, the expectation maximization (EM) algorithm is conventionally used, since it provides only positive solutions for the feeding distributions and no additional reg-ularization parameters (or assumptions) are required to solve the TAGS inverse problem.
Clearly, the first level scheme (or defined branching ratio matrix) considered is not necessarily the one that will provide a nice description of the measured TAGS spectrum. For that reason, as part of the analysis the cut-off energy and the parameters that define the branching ratio matrix can be varied until the best description of the experimental data is obtained. Also assumptions on the spin and parity of the ground state of the parent nucleus and on the spins and parities of the populated levels can be changed when they are not known unambigously, since to connect the levels in the continuum to levels in the known part of the level scheme we need information about their spin and parity. All these changes provide different branching ratio matrices (or daughter level schemes) that are considered during the analysis and for all of them the corresponding response matrixes are calculated and Equation 1 is solved. The final analysis is then based on the level scheme (or branching ratio matrix) that is consistent with the available information from high resolution measurements and at the same time provides the best description of the experimental data. So in practical terms the following steps are followed until the best description of the data is obtained: a) define a branching ratio matrix B, b) calculate the corresponding response matrix R ij (B), and c) solve the corresponding Equation 1 using an appropriate algorithm d) compare the generated spectrum after the analysis (R(B)f + C) with the experimental spectrum d.
We have only mentioned briefly how the response function R ij (B) is calculated. More specifically the response for each level can be determined recursively starting from the lowest level in the following way [20]: where R j is the response to level j, g jk is the response of the gamma transition from level j to level k which is calculated using Monte Carlo simulations, b jk is the branching ratio for the gamma transition connecting level j to level k, and R k is the response to level k. Here the index k runs for all the levels below the level j. For simplicity we have not included here in the formula the convolution with the response of the beta particles and only the gamma part of the response is presented. Note that in this last notation R j is a vector that contains as elements the R ij matrix elements mentioned above for all possible i-s (or channels) of the TAGS spectrum and the branching ratio matrix enters in the formula of the response matrix through the decay branches b jk -s. In the real calculation of the responses the internal conversion process is also taken into account.
Prior to the analysis, the contaminants in the TAGS spectrum (C i ) have to be isolated and their individual contributions evaluated. The nucleus to be studied is produced by nuclear reactions together with a number of additional nuclei. Two alternative separation methods are normally used to isolate the nucleus of interest. On-line mass separators are used with low energy radioactive beams to reduce the contamination represented by mass isobars. In-flight separation is used at high-energy fragmentation facilities to reduce the number of nuclear species in the "cocktail beam" to suitable levels. Even if we can isolate the nucleus of interest, daughter (and other descendants) activity can contaminate the measured spectrum depending on the half-life of the studied decay. This contamination can be determined through dedicated measurements on the decay of the contaminant nuclei under the same conditions as the one of interest. Another source of contamination of the spectrum is the pile-up of signals. The pile-up can distort the full TAGS spectrum and can generate counts in regions of the spectra where there should not be counts, as for example in the region beyond the Q β value of the decay. Also it can distort the spectrum in regions where we expect reduced statistics as for example close to the Q β value of the decay. This is the reason why estimating this contribution is of importance. Algorithms have been developed to evaluate this contribution [21,22]. Its determination is based on the random superposition of true detector pulses, measured during the experiment, within the time interval defined by the acquisition gate of the data acquisition system.
Another possible contamination appears when the decay is accompanied by beta delayed particle emission, since this process can lead promptly to the emission of gamma rays from the final nucleus populated by the beta delayed particle emission. The case of the emission of beta delayed neutrons is even more complex. Neutrons interact easily with the detector material and release their energy through inelastic and capture processes. The proper evaluation of this contamination is of great relevance in the study of beta decays far from stability on the neutron-rich side of the Segre chart and requires careful Monte Carlo simulations of the neutron-detector interactions [22,23]. The reproduction of this contamination is complicated because it has two components: one, which is prompt with the beta decay, is composed of gamma rays emitted in the final nucleus after the beta-delayed neutron emission when an excited state is populated, the other component due to neutron interactions in the detector is delayed, since the speed of neutrons is much lower than that of gamma-rays. To simulate these effects properly an event generator [24], that takes into account relative contribution of the two components is required. It is also necessary to know the energy spectrum of the emitted beta-delayed neutrons. In addition, the Monte Carlo simulation code should include an adequate physics model of the neutron interactions. As an example, in Figure 4 [25] the contribution of the calculated beta delayed neutron contamination to the TAGS decay spectrum of 95 Rb is presented. Two available neutron energy spectra were used in the simulations [26,27], and clearly only one reproduces the experimental TAGS data at high excitation energies. This figure shows the relevance of the neutron spectrum used in the simulations (for more details see [25]). Due to these complications we have built Rocinante [24,28] Fig. 4. Impact of the neutron energy spectrum (In) in the simulations of the contamination associated with the beta delayed neutrons in the TAGS spectrum (for more details see [25]). Only the spectrum measured by Kratz et al. [26] reproduces the TAGS spectrum at high excitation energies. The Monte Carlo (MC) spectra are normalized to the experimental spectrum around the neutron capture peak indicated with an arrow. The prompt 836.9 keV γ-ray peak from the first excited state in the final nucleus after the beta-delayed neutron emission 94 Sr is highlighted. Reprinted figure with permission from [25], Copyright (2019) by the American Physical Society.
material, aimed at the measurement of beta-delayed neutron emitters. BaF 2 has a neutron capture cross-section one order-of-magnitude smaller than the NaI(Tl), that is conventionally used. This spectrometer was also the first of a new generation of segmented devices designed to exploit the cascade multiplicity information to improve the TAGS analysis, as will be mentioned later.
It is important to first identify the different distortions or contaminations, but it is also important to determine properly their corresponding weight in the measured spectrum. Depending on the distortion, different strategies have been followed. For example, the contribution from contaminant decays can be evaluated if there is a clear peak identified in the spectrum that comes from this contamination that can be used for normalization. Another option is the assessment of this contribution from the solution of the Bateman equations, using the information on half-lives and measurement conditions (collection and measuring cycle times). In the case of high-energy fragmentation experiments where the contamination is due to beta-gamma events uncorrelated with the implanted ion it can be evaluated from correlations backward in time. The pileup distortion can be evaluated based on the number of counts in the TAGS spectra which lie beyond the highest Q β value in the decay chain and which are clearly above the contribution of the background, since we can assume that those counts can only come from this contribution. When this option is not possible because of inadequate statistics, a procedure is given in [21] for the normalisation of this contribution based on the counting rate and the length of the analogue to digital converter (ADC) gate. And finally if there is a contamination arising from beta-delayed neutrons, this contribution can be normalized to the broad high-energy structure generated by neutron captures in the detector material when possible, otherwise it should be normalized to the Pn value of the decay.
In Figure 5, we present as an example a total absorption spectrum measured during our first experiment in Jyväskylä of the decay of 104 Tc [7,8] which is relevant for the decay heat application (see Section 3). In the upper panel of this figure we show the spectrum of this decay compared with the reproduction of the spectrum after the analysis and the contribution of the contaminants (back-ground+daughter activity+pileup).
In this measurement a TAGS detector that consisted of two NaI(Tl) cylindrical crystals with dimensions: = 200 mm × l = 200 mm, and = 200 mm × l = 100 mm was used (courtesy of Dr. L. Batist). The longer crystal has a longitudinal hole of = 43 mm for the positioning of the sources in the approximate geometrical centre of the spectrometer using a tape system. In the experiment the crystals were separated by 5 mm. This separation and the ideal position of the sources inside the spectrometer was studied previous to our experiment using Monte Carlo simulations in order to maximize the gamma efficiency of the setup [29]. This TAGS had a 57% peak and 92% total efficiency for the 662 keV gamma transition emitted in 137 Cs decay and 27% peak and 70% total efficiency for a 5 MeV gamma transition. This last value was obtained from previously validated Monte Carlo simulations. The efficiency of this setup is modest compared with recently developed total absorption spectrometers such as DTAS [30], or MTAS [31]. This detector was designed at the Nuclear Institute of St. Petersburg (Russia) [32].
This measurement was analyzed in singles, since the precision and reproducibility of the tape positioning system was considered not sufficiently good to allow coincidence counting. The positioning of the sources is critical in the determination of the efficiency of the Si detector used as an ancilliary detector for coincidences with the beta particles emitted in the decay. The efficiency of the beta detector as a function of end-point energy has a direct impact on the normalization of the combined betagamma cascade response of the spectrometer (R ij ). Using singles has the advantage of providing much higher statistics in the analysis compared with gated spectra. However, the use of gated spectra is preferred, eventually unavoidable, in order to reduce contamination from ambient background and the selection of events of particular interest.
The lower panel of Figure 5 shows the feeding distribution deduced for the 104 Tc decay obtained from the TAGS measurement compared with the distribution obtained from high resolution measurements. From the figure it is clear that the feeding distribution obtained with the TAGS is shifted to higher energies in the daughter nucleus, which is typical of a case suffering from the Pandemonium effect [7,8,33]. Similar measurements will be discussed in more detail in Section 3.
As mentioned earlier, in this measurement a detector was used that was composed of two crystals. The new generation of available detectors such as Rocinante [24], SUN [34], MTAS [31] and DTAS [30] exploit segmentation to a greater degree to extract additional information from the decay under study. Using the segmentation it is possible to measure the detector fold (number of detectors fired in an event) which is related to gamma-cascade multiplicity as a function of excitation energy and ultimately to the de-excitation branching ratio matrix B. The lack of knowledge of the matrix B is the largest source of uncertainty in TAGS analysis and this can be greatly improved with segmented detectors. Our current approach to the iterative procedure for updating B described earlier in this Section, is to include in step d) the comparison to fold-gated TAGS spectra and single module spectra. Reconstructed fold-gated spectra are obtained by MC simulation using the appropriate event generator since it is not possible to define a fold-gated response in a manner similar to Equation 2. This also prevents us from including them as part of the inverse problem (Equation 1). A different approach has been taken by the ORNL group to analyze MTAS data [35,36]. They use the coincidence between one module and the sum of all the modules to define total energy gated single detector spectra that are fitted by the sum of a number of de-excitation cascades, usually taken from high resolution spectroscopy and supplemented when necessary with "pseudo" levels with guess branching ratios and modified iteratively until the best reproduction is achieved. Yet another approach is used by the NSCL group to extract B from SUN data [37]. They start from the same total energy gated single detector spectra but apply the so called Oslo-method [38] to obtain the branching ratio matrix for a subset of levels. Because of this, the TAGS analysis is not performed with this B but uses the "pseudo" level approach including in the fit the total absorption spectrum and the spectrum of detector multiplicities [39]. It should be noticed that the traditional Oslo method is not strictly applicable to TAGS data because the assumed equivalence of total deposited energy with excitation energy does not hold in general, due to the non-ideal detector response. Currently we are working in a method to solve the full non-linear inverse problem represented by Equation 1, to obtain feedings and branching ratios from the complete data set provided by a segmented spectrometer: sum energy spectrum gated by detector fold, sum energy spectrum versus single crystal spectrum and crystal-crystal correlations.
In Figure 6 we show the spectrum of the beta decay of 100 Tc measured in a recent campaign of measurements at the IGISOL IV facility of the Univ. of Jyväskylä [40,41] with the segmented DTAS detector. This single decay is part of the A=100 system of relevance for double beta decay studies ( 100 Ru -100 Tc -100 Mo). Previous to this study, only a high resolution measurement existed for this single decay and there were doubts whether feeding at high excitation energy is not detected in the high resolution measurements. Single decays, like this one can be of relevance for fixing model parameters used in theoretical calculations for neutrino and neutrinoless double beta decay studies. Our TAGS results show only a modest improvement with relation to the earlier high resolution  104 Tc with the spectrum generated after the analysis (reconstructed spectrum). This last spectrum is obtained by multiplying the response function of the decay with the determined feeding distribution (R(B)f f inal ). The lower panel shows the beta-decay feeding distribution obtained compared with that previously known from high resolution measurements [7,8,33]. Reprinted figure with permission from [7], Copyright (2010) by the American Physical Society.
results, revealing that this decay did not suffer seriously from the Pandemonium systematic error (see Figure 7). This decay is not only important in the framework of double beta decay studies, it has also recently attracted attention in another neutrino related topic [42]. The decay is a relevant contributor in a newly identified flux-dependent correction to the antineutrino spectrum produced in nuclear reactors that takes into account the contribution of the decay of nuclides that are produced by neutron capture of long lived fission products. In this particular case 99 Tc is produced as a fission product, which after neutron capture becomes 100 Tc that beta decays. The effect has a nonlinear dependence on the neutron flux, because first a fission is required and later a neutron capture. Effects like this one are considered in order to explain features of the predicted antineutrino spectrum for reactors not yet fully understood (see Section 4).
The study of this decay was the first time that the DTAS detector was used at a radioactive beam facility. Prior to the analysis of this case a full characterization of the detector was performed [22,40]. This included a check on the ability to reproduce with MC simulations the spectrum of decays obtained with different detector multiplicity (fold) conditions. As an example we present in Figure 8 the reproduction of the multiplicities for the 22 Na source used in the characterization of the detector. The DTAS is constructed in a modular way that adds extra versatility to the setup [30]. Depending on the installation, it can be used in an 18 detector configuration for ISOL type facilities or in an 16 detector configuration for fragmentation facilities, where the positioning of the implantation detectors normally requires more space.  Fig. 7. Comparison of the obtained TAGS feeding distribution in the decay of 100 Tc with the data available from high resolution measurement. This is an example of a case that did not suffer from the Pandemonium effect. Reprinted figure with permission from [41], Copyright (2017) by the American Physical Society.

Decay heat
Nuclear reactor applications require beta decay data. The relevance of beta decay is shown by the fact that each fission is followed by approximately six beta decays. The energy balance released in fission is presented in Table 1 for 235 U and 239 Pu fissile isotopes [43]. In the case of 235 U, for example, 7.4 % of the energy released comes from the beta decay of the fission products (FP) (gamma and beta energy). Depending on the composition of the fuel in the reactor this percentage can change, but it is of the order of 7% of the total released energy for a working reactor. Once the reactor is shut-down, the decay energy becomes dominant and the related heat has to be removed. If for some reason this is not possible, it can produce accidents like the one caused originally by the tsunami that followed the Great East Japan Earthquake (2011) in the Fukushima Daiichi power plant. Clearly one needs to estimate this source of energy for the safety of any nuclear installation  [22,40]. How well the different multiplicity spectrum is reproduced, is a stringent test of the quality of the branching ratio matrix used in the analysis.
(design of a reactor, storage of the nuclear waste, loss of coolant accident (LOCA), etc.). Decay heat is defined as the amount of energy released by the decay of fission products not taking into account the energy taken by the neutrinos. The first method to estimate the decay heat was introduced by Way and Wigner [44], which was based on statistical considerations of the fission process. Their results provide a good estimate of the heat released, but the precision reached is not sufficient for present-day safety standards. Nowadays the most extended way to estimate the decay heat is to perform summation calculations, which relies on the increased amount of available nuclear data. In this method, the power function of the decay heat f (t) is obtained as the sum of the activities of the fission products times the energy released per decay: where E i is the mean decay energy of the ith nuclide (β or charged-particle and γ or electromagnetic components), λ i is the decay constant of the ith nuclide, and N i (t) is the number of nuclides of type i at the cooling time t (for simplicity the α-decays of minor actinides are not included here). These calculations require extensive libraries of cross sections, fission yields and decay data, since the method first requires the solution of a system of coupled differential equations to determine the inventory of nuclei N i (t) produced in the working reactor and after shut-down.
Several ingredients of this method depend on decay data. The determination of the activities of the fission products (λ i N i (t)) requires a knowledge of the half-lives of the decaying isotopes. The other important quantities are the mean energies released per decay (E β,i , E γ,i ). The mean energies released per decay i can be obtained by direct measurements as in the systematic studies by Rud- stam et al. [45] and Tengblad et al. [46]. These integral measurements (energy per decay) require specific setups that are only sensitive to the energy of interest and a careful treatment of all possible systematic errors. Alternatively the mean energies can be deduced from available decay data in nuclear databases such as the Evaluated Nuclear Structure Data File (ENSDF) [47] if the decay properties are properly known. The term "properly known" beta decay implies a knowledge of the Q β value of the decay, the half-life, the beta distribution probability to the levels in the daughter nucleus and the decay branching ratios of the populated levels. If all this information is Table 3. List of parent nuclides identified in [50] that should be measured using the total absorption technique to improve the predictions of the decay heat in reactors based on 233 U/ 232 Th fuel. The list does not contain several relevant cases already measured by [17] and already included in Table 2 (marked with †), for more details see [50]. Rel. (relevance) stands for the priority of the measurement. Isotopes marked with asterisks show the measurements performed by our collaboration. For more details in the notation see Table 2. available, then it is possible to deduce the mean energies released by the decay using the following relations: where E j is the energy of the level j in the daughter nucleus, I j is the probability of a beta transition to level j, and < E β > j is the mean energy of the beta continuum populating level j. As can be seen from the formula the mean gamma energy is approximated by the sum of the energy levels populated in the decay weighted by the beta transition probability. This approximation assumes that each populated level decays by gamma deexcitation and ignores conversion electrons which are taken into account in the complete treatment of the mean energy calculations. The mean beta energy, because of the continuum character of the beta distribution emitted in the population of each level, requires the determination of the released mean energy < E β > j for each end-point energy of the beta transition (Q β − E j ). Then the mean beta energy (E β ) is obtained as the weighted sum of the mean beta energies populating each level by the beta transition probability. For the determination of < E β > j for each level one needs to make assumptions about the type of the beta transition (allowed, first forbidden, etc.) and the knowledge of the Q β value of the decay is needed to determine the beta transition end-points.
Pandemonium can have an impact in the determination of the mean energies from data available in databases. If the beta decay data suffers from the Pandemonium effect the beta decay probability distribution is distorted. This distortion, which implies increased beta probability to lower lying levels in the daugther nucleus, causes an underestimation of the mean gamma energy and an over-estimation of the mean beta energy. This is why TAGS measurements are relevant to this application.
In fission more than 1000 fission products can be produced. But not all of them are equally important. When addressing a particular problem, like the decay heat, it is of interest to identify which are the most relevant contributors among the large number of fission products. A group of experts working for the International Atomic Energy Agency (IAEA) [48] identified high priority lists of nuclei that are important contributors to the decay heat in reactors and that should be measured using the TAGS technique. These lists included nuclides that are produced with high yields in fission and for which the decay data was suspected of suffering from the Pandemonium effect. One argument used for this last selection was if the decay data shows no levels fed in the daughter nucleus in the upper 1/3 excitation energy window of the Q β value. It is worth noting that this can be considered only as an indication of Pandemonium and not a rigorous rule. Another way of looking for questionable (or odd) data was to look for cases that show different mean energies in the different international databases. The final lists were published in two separate reports, first for U/Pu fuels [49] and then later for the possible future Th/U fuel [50].
In 2004 we started a research programme aimed at studying the beta decay of nuclei making important contributions to the decay heat in reactors. For the planning of any nuclear physics experiment the first step is to decide the best facility to perform it in terms of the availability of the beams, their cleanness and their intensity. Since some of the most important contributors to the priority list for 235 U and 239 Pu fission were refractory elements like Tc, Mo and Nb, the options were very limited. In a classical ISOL facility like ISOLDE, the development of a particular beam can take some time if the beam of a particular element is not available. It is a lengthy and complex task to find the optimum chemical and physical conditions in the ion source for the extraction of a particular element. We decided that the best option concerning the availability of the beams was to perform the measurements at the IGISOL facility in Jyväskylä (Finland) [51]. The reason for that was the development of the ion-guide technique. The ion-guide technique, and more specifically the fission ion guide, allows the extraction of fission products independently of the element. In this technique, fission is produced by bombarding a thin target of natural U with a proton beam. The fission products that fly out of the target are stopped in a gas and transported through a differential pumping system into the first accelerator stage of the mass separator. The dimensions of the ion guide and the pressure conditions are optimized in such a way that the process is fast enough for the ions to survive as singly charged ions. As a result the system is chemically insensitive and very fast (sub-ms) [52] allowing the extraction of any element including those that are refractory.
Another important advantage of performing experiments at IGISOL is the availability of the JYFLTRAP Penning trap [53] developed for high precision mass measurements at this facility. JYFLTRAP can also be used as a high resolution mass separator for trap assisted spectroscopy measurements, providing a mass resolving power ( M ∆M ) of the order of 100 000 to be compared with the resolving power of approximately 500 of the IGISOL separator magnet. The purity of the beams is particularly important for calorimetric measurements like those with TAGS since it reduces systematic errors that can be associated with contamination of the primary radioactive beams. This advantage has also been used in other types of calorimetric measurements at IGISOL such as the measurements of beta delayed-neutrons using 3 He counters embedded in a polyethylene matrix [54].
Three experimental campaigns have been performed at the IGISOL facility to study the beta decay of important contributors to the decay heat and to the antineutrino spectrum in reactors using the TAGS technique [55,56,57,58]. One of the total absortion setups used in the experiments is presented in Figure 9 (Rocinante TAGS). In a typical experiment, the radioactive beam extracted from IGISOL is first mass separated using the separator magnet and then further separated using the JYFLTRAP Penning trap. Then the beam is transported to the measuring position, at the centre of the total absorption spectrometer where it is implanted in a tape from a tape transport system. The tape is moved in cycles, which are optimized depending on the half-life of the decay of interest. As mentioned earlier, the reason for using a tape transport system is to reduce the effect of undesired daughter, grand-daughter, etc., decay contaminants in the measured spectrum. If necessary, these contaminants have to be substracted from the measured TAGS spectrum and require dedicated measurements. In this kind of measurement the TAGS detector is usually combined with a beta detector as shown in the inset to the Figure 9. The beta detector is used to select coincidences of the beta particles with the TAGS spectrum, which essentially eliminates the effect of the ambient background.
In Figure 10 we show an example of the recently measured 86 Br decay, which was considered priority one for the decay heat problem. The spectrum shows the TAGS spectrum obtained with the beta-gate condition, and the contribution of the different contaminants. The analysis was performed as described earlier. Known levels up to the excitation value of 3560 keV were taken from the compiled high resolution data from [59]. From 3560 keV to the Q β =7633(3) keV value, the statistical model is used to generate a branching ratio matrix using a level density function resulting from a fit to levels from [60] and corrected to reproduce the level density at low excitation energy, and E1, M1, and E2 gamma strength functions taken from [61]. For more details see [62,63]. In Figure 10 we also show the comparison of the beta gated TAGS spectrum with the results from the analysis. The reconstructed spectra are obtained by multiplying the response function of the detector with the final feeding distribution obtained from the analysis. In this particular case two results are presented. Response A corresponds to the conventionally calculated branching ratio matrix that fits better the experimental spectrum. Response B corre- sponds to a modified branching ratio matrix to reproduce the measured gamma-ray intensities de-exciting each level as measured in high resolution experiments. In Figure 11 the feeding distributions obtained are compared with the available high resolution measurements. This comparison shows that 86 Br decay was suffering from the Pandemonium effect.
In Table 4 we show a summary of the mean energies deduced from TAS analyses obtained in our measurements performed at Jyväskylä. It shows that most of the cases addressed from the high priority list were suffering from the Pandemonium effect. Two cases, that originally were suspected to be Pandemonium cases ( 102 Tc, 101 Nb), were not. Those cases also show the necessity of the TAGS measurements to confirm the suspicion of Pandemonium. Clearly the non existence of feeding in the last Q β /3 excitation is a good indicator to select cases, but it is not always conclusive. The 102 Tc, 101 Nb cases have strong ground state feedings, which reduces the impact of the undetected gamma branches at high excitation. This is clearly reflected in the differences of the deduced mean energies, when they are compared with the mean energies deduced from high resolution data. The values for one-third of the Q β -value (Q β /3) are also given for comparison with the mean energies. The Q β /3 value is an approach sometimes used by databases, when there is lack of experimental data, which in practical terms divides the available decay energy in equal parts between the mean gamma, beta and antineutrino energies. In the table the mean energies deduced from the high resolution data (ENDSF database) are also given for comparison [64].
Our results can also be compared where possible with the results of Greenwood et al. [17] and Rudstam et al. [45]. Greenwood and co-workers performed a systematic study of fission products at the Idaho National Engineering Laboratory (INEL), Idaho Falls, USA, using a 252 Cf source and the He-jet technique. They employed a total absorption spectrometer built of NaI(Tl) with the following dimensions, 25.4 cm ×30.5 cm length with a 5.1 cm × Table 4. Mean gamma and beta energies deduced from our analyses of beta decays studied at Jyväskylä in comparison with the values deduced from high resolution measurements (ENSDF database). The highest level identified in the decay studies using high resolution and the decay Q β values are also given for completeness ( 20.3 cm long axial well. The analysis technique employed in those studies was based on the creation of a level scheme using the information from high resolution measurements and complementing it with the addition of "pseudo" levels by hand when necessary. The response of the detector was then calculated for the assumed level scheme using the Monte Carlo code CYLTRAN (for more details see [65]). The systematic study of Greenwood et al. provided TAGS data for 48 decays including the decay of three isomeric states. Since a different analysis technique was used, it is interesting to compare the results of Greenwood with the more recent results and look for possible systematic deviations for cases where the comparison is possible. As mentioned earlier Rudstam et al. performed systematic measurements of gamma and beta spectra and deduced mean energies [45,46] of interest for the prediction of decay heat. Beta spectra of interest for the prediction of the antineutrino spectrum from reactors were also measured [46]. The direct measurements were performed at ISOLDE and at the OSIRIS separator using setups optimized for the detection of the gamma-and beta-rays emitted by the fission products. In the case of the mean gamma energies, the setup required an absolute gamma efficiency calibration with the assumption that the decay used in the calibration did not suffer from the Pandemonium effect [45]. In their measurements the beta decay of 91 Rb was used as a calibration. This decay with a Q β of 5907(9) keV shows a complex decay scheme from high resolution measurements that populates levels up to 4700 keV in the daughter nucleus. So it was assumed naturally by the authors of [45] that this decay was not a Pandemonium case (in the publication [45] it is mentioned that the gamma spectrum extends up to 4500 keV).
In a contribution to the Working Party on International Evaluation Co-operation of the NEA Nuclear Science Committee (WPEC 25) group activities [48], the late O. Bersillon performed a comparison between the Greenwood and Rudstam mean gamma and beta energies for Differences between the mean gamma energies obtained with TAGS measurements (see Table 5 and the text for more details) and the direct measurements of Rudstam et al. [45] after renormalization by a factor of 1.14, which was deduced from the comparison of our newly determined mean energy for the decay of 91 Rb [62] with that employed by Rudstam et al. [45]. Blue points represent Greenwood TAGS data and red points represent results from our collaboration. A systematic difference with a mean value of -180 keV remains.
Reprinted figure with permission from [62], Copyright (2017) by the American Physical Society.
those decays that were available in both data sets. This comparison was revisited in [66]. A clear systematic difference was shown, pointing to possible systematic errors in one or in both data sets [24,62,66]. For that reason it was decided to measure the 91 Rb decay using the TAGS technique, to see if the decay was suffering from the Pandemonium effect or not and accordingly check if this decay was adequate as an absolute calibration point to obtain mean gamma energies in [45]. The 2669(95) keV mean gamma value obtained from our measurements [62] can be compared with the value used by Rudstam et al. (2335(33) keV) showing that this decay suffered from Pandemonium and also showing the necessity of renormalizing the Rudstam data by a factor of 1.14. With this renormalization, the mean value of the differences between the two data sets (TAGS vs Rudstam) reduces from -360 keV to -180 keV, but still the discrepancy remains [62]. This is shown in Figure 12 and  [7,8,24,62,67] is presented in comparison with the decay heat measurements reported by Tobias [68] and Dickens [69] for the electromagnetic component (E γ ) of the decay heat. Similarly in Figs. 15-16 the impact of the same decays is compared for the charged particle component (E β ) of the decay heat. To show the impact of the total absorption measurements the data base JEFF3.1.1 is used, in which no total absorption data is included. Calculations are performed using the bare JEFF 3.1.1 and a modified version of the JEFF3.1.1 database with the inclusion of the total absorption data for the mean energies.  The results presented were provided by Dr. L. Giot [70]. In the figures, it can be noted the large impact of the mentioned decays and the relevance of the total absorption measurements for a proper assesment of the decay heat based on summation calculations.  Fig. 15. Impact of the inclusion of the total absorption measurements performed for 13 decays in the beta component of the decay heat calculations for 239 Pu (see Figure 13 for more details).

Neutrino applications
Nuclear reactors constitute an intense source of electron antineutrinos, with typically 10 20 antineutrinos per second emitted by a 1GWe reactor. The reactor at Savannah River was the site of the discovery of the neutrino in 1956 by Reines and Cowan [71], thus confirming Pauli's predictions of twenty-six years earlier [1]. Just like the decay heat described above, antineutrinos arise from the beta decays of the fission products in-core. Their energy spectrum and flux depend on the distribution of the fission products which reflects the fuel content of a nuclear reactor. This property combined with the fact that neutrinos are sensitive only to the weak interaction could make an- tineutrino detection a new reactor monitoring tool [72]. Both particle and applied physics are the motivations of their study at power or research reactors nowadays with detectors of various sizes and designs placed at short or long distances. In the last decade, three large neutrino experiments with near and far detectors, were installed at Pressurized Water Reactors [73,74,75], to try to pin down the value of the θ 13 mixing angle parameter governing neutrino oscillations. These experiments have sought the disappearance of antineutrinos by comparing the flux and spectra measured at the two sites, both distances being carefully chosen to maximise the oscillation probability at the far site. The three experiments [73,74,75] have now achieved a precise measurement of the θ 13 mixing angle, paving the way for future experiments at reactors looking at the neutrino mass hierarchy or for experiments at accelerators for the determination of the delta phase, that governs the violation of the CP symmetry in the leptonic sector, thus shedding light on why there is an abundance of matter rather than antimatter in the Universe.
Though they have used one or several near detectors in order to measure the initial flux and energy spectrum of the emitted antineutrinos, the prediction of the latter quantities still enters in the systematic uncertainties of their measurements, because their detectors are usually not placed on the isoflux lines of the several reactors of the plant at which they are installed [76]. In addition, the Double Chooz experiment started to take data with the far detector alone, implying the need to compare the first data with a prediction of the antineutrino emission by the two reactors of the Chooz plant. Two methods employed to calculate reactor antineutrino energy spectra were revisited at that time, i.e. the conversion and the summation methods.
The conversion method consists in converting the integral electron spectra measured at the research reactor at the Institute Laue-Langevin (ILL) in Grenoble (France) by Schreckenbach and co-workers [77,78,79,80] with 235 U, 239 Pu and 241 Pu thin targets under a thermal neutron flux. These spectra exhibit rather small uncertainties and remain a reference as no other comparable measurement has been performed since. Being integral measures, no information is available on the individual beta decay branches of the fission products. This prevents the use of the conservation of energy to convert the beta into antineutrino spectra. Schreckenbach et al. developed a conversion model, in which they used 30 fictitious beta branches spread over the beta energy spectrum to convert their measurements into antineutrinos. In 2011 Mueller et al. [81] revisited the conversion method and improved it through the use of more realistic end-points and Z distributions of the fission products, available thanks to the wealth of nuclear data accumulated over 30 years, and through the application of the corrections to the Fermi theory at branch level in the calculation of the beta and antineutrino spectra. After these revisions, the prediction of detected antineutrino flux at reactors compared with the measurements made at existing short baseline neutrino experiments revealed a deficit of 3%. The result was confirmed immediately by Huber [82] who carried out a similar calculation though he did not explicitly use beta branches from nuclear data. This antineutrino deficit was even increased by the revision of the neutron lifetime and the influence of the long-lived fission products recalculated at the time, to finally amount to 6% [83]. A new neutrino anomaly was born: the reactor anomaly. Several research leads were followed since to explain this deficit. An exciting possibility is the oscillation of reactor antineutrinos into sterile neutrinos [83], which has triggered several new experimental projects worldwide [84,85,86]. In 2015, the mystery deepened when Daya Bay in China, Double Chooz in France and RENO in Korea, reported the detection of a distortion (colloquially called bump) in the measured antineutrino energy spectrum with respect to the converted spectrum, which could not be explained by any neutrino oscillation. The three experiments rely on the same detection technique and similar detector designs, which make it possible that they would all suffer from the same detection bias as suggested in [87]. But the three collaborations have thoroughly investigated this hypothesis without success. In the face of the observed discrepancies between the converted spectra and the measured reactor antineutrino spectra, it is worth considering more closely the existing methods used to compute them and other possible explanations like the case of nonlinear corrections discussed in Section 2 in relation to the 100 Tc decay [42].
Converted spectra rely on the unique measurements performed at the high flux ILL research reactor with the high resolution magnetic spectrometer BILL [88], using thin actinide target foils exposed to a thermal neutron flux that was well under control. This device was exceptional as it allowed the measurement of electron spectra ranging from 2 to 8 MeV in 50 keV bins (smoothed over 250 keV in the original publications) with an uncertainty dominated by the absolute normalization uncertainty of 3% at 90% C.L. except for the highest energy bins with poor statistics [77,78,79,80]. The calibration of the spectrometer was performed with conversion electron sources or (n,e − ) reactions on targets of 207 Pb, 197 Au, 113 Cd and using the beta decay of 116 In providing calibration points up to 7.37 MeV. The irradiation duration ranged from 12 hours to 2 days. Two measurements of the 235 U electron spectrum were performed, the first one lasting 1.5 days and the second one 12 hours. The normalisation of the two spectra disagree because they were normalized using two different (n,e − ) reactions on 197 Au and 207 Pb respectively in chronological order. The measurement retained by the neutrino community is the second one. The conversion procedure consists in successive fits of the electron total spectrum with beta branches starting with the largest end-points. The total electron spectrum is fitted iteratively bin by bin starting with the highest energy bins, and the contributions to the fitted bin are subtracted from the total spectrum. The reformulation of the finite size corrections, as well as a more realistic charge distribution of the fission products and a much larger set of beta branches have been the key for the newly obtained converted antineutrino spectra of [81] and [82]. But the possibility remains that the electron and/or converted spectra suffer from unforeseen additional uncertainties. Indeed the normalisation of the electron spectra relies on the (n,e − ) reactions quoted above and on internal conversion coefficient values that may have both been re-evaluated since [89]. In addition, the exact position of the irradiation experiment in the reactor is not well known and may have an impact on the results as well [89]. Another concern is associated with the conversion model itself, where uncertainties may not take into account missing underlying nuclear physics. In Mueller's conversion model, forbidden non-unique transitions are replaced by forbidden unique transitions (when the spins and parities are known!). The shapes of the associated beta and antineutrino spectra are not well known and the forbidden transitions dominate the flux and the spectrum above 4 MeV. Several theoretical works have attempted to estimate the uncertainties introduced by this lack of knowledge [90,91,92]. The latest study [92] reports a potential effect compatible with the observed shape and flux anomalies. Another source of uncertainties comes from the weak magnetism correction entering in the spectral calculation [82,93] that is not well constrained experimentally in the mass region of the fission products. These two extra uncertainties affect converted spectra and are not included in the published uncertainties. Eventually the conversion process itself could be discussed, as the iterative fitting procedure is not the only possible conversion method and it is suspected of inducing additional uncertainties [94].
In order to identify what could be at the origin of these anomalies, the understanding of the underlying nuclear physics ingredients is mandatory. Indeed, only the decomposition of the reactor antineutrino spectra into their individual contributions and the study of the missing underlying nuclear physics will allow us to understand fully the problem and provide reliable predictions. The best tool to address these questions is to use the summation method. This method is based on the use of nuclear data combined in a sum of all the individual contributions of the beta branches of the fission products, weighted by the amounts of the fissioning nuclei. Two types of datasets are thus involved in the calculation: fission product decay data, and fission yields. This method was originally developed by [95] followed by [96] and then by [46,97]. The β/ν spectrum per fission of a fissile isotope S k (E) can be broken down into the sum of all fission product β/ν spectra weighted by their activity λ i N i (t) similarly to what is done for decay heat calculations: Eventually, the β/ν spectrum of one fission product (S i ) is the sum over the β branches (or beta transition probabilities) of all β decay spectra (or associatedν spectra), S b i (in equation 6), of the parent nucleus to the daughter nucleus weighted by their respective beta branching ratios according to: where f b i represents the beta transition probability of the b branch, Z i and A i the atomic number and the mass number of the daughter nucleus respectively and E b 0i is the endpoint of the beta transition b. In 1989 the measurement of 111 beta spectra from fission products by Tengblad et al. [46] was used for a new calculation of the antineutrino energy spectra through the summation method. But the overall agreement with the integral beta spectra measured by Hahn et al. [80] was at the level of 15-20% showing that a large amount of data were missing at that time. Lately, the summation calculations were re-investigated using updated nuclear databases. Indeed the summation method is the only one able to predict antineutrino spectra for which no integral beta measurement has been performed. The existing aggregate beta spectra needed to apply the conversion method are relatively few and were measured under irradiation conditions that are not exactly the same as those existing in power reactors. Among the discrepancies, the energy distributions of the neutrons generating the fissions in the ILL experiments are different from those in actual power reactors, and even more from the ones in innovative reactor designs such as fast breeder reactors. The aggregate beta spectra were measured for finite irradiation times much shorter than the typical times encountered in power reactors. These few spectra and the specific conditions are not usable for innovative reactor fuels or require corrections for longer irradiation times (called off-equilibrium corrections) or more complex neutron energy distributions incore. Until the recent measurement of the 238 U beta spectrum at Garching by [98], the conversion method could not be applied to obtain a prediction of the 238 U fast fission antineutrino spectrum. This was one of the motivations for the first re-evaluation of the summation spectra that was performed in Mueller et al., the second being to provide off-equilibrium corrections [81] to the converted spectra. In this work, several important conclusions were already listed regarding summation calculations for antineutrinos.
The evaluated nuclear databases do not contain enough decay data to supply detailed beta decay properties for all the fission products stored in the fission yields databases. The evaluated databases have thus to be supplemented by other data or by model calculations for the most exotic nuclei. The relative ratio of the aggregate beta spectra with the obtained summation spectra from databases exhibited a shape typical of the Pandemonium effect, with an overestimate of the high energy part of the spectra in the nuclear data. The maximum amount of data free of the Pandemonium effect should thus be included in the summation calculations. The difficulty comes from the fact that these Pandemonium-free data are usually not included in the evaluated databases. One has thus to gather the existing decay data and compute the associated antineutrino spectra. The Pandemonium-free data are mostly existing TAGS measurements [17] and the electron spectra directly measured by Tengblad et al. [46]. They were included in an updated summation calculation performed in [99], in which the seven isotopes measured with the TAGS technique that had so much impact on the 239 Pu electromagnetic decay heat, i.e. 105 Mo, 102,104−107 Tc, and 101 Nb [7], were taken into account. The calculation revealed that these TAGS results had a very large impact on the calculated antineutrino energy spectra, reaching 8% in the Pu isotopes at 6 MeV. But it appeared that summation calculations still overestimate the beta spectra at high energy, indicating that there were large contributions from nuclei where the data suffer from the Pandemonium effect in the decay databases. The situation is thus similar to that already encountered in the decay heat summation calculations. These conclusions reinforced the necessity for new experimental TAGS campaigns and spread the message worldwide. New summation calculations were developed and other experimental campaigns were launched using the TAGS technique [100,101,102]. In [100] a careful study of the existing evaluated fission yield databases was performed. It appeared that the choice of the fission yield database had a large impact on the summation spectra obtained, because of mistakes identified in the ENDF/B-VII.1 fission yields for which corrections were proposed. Once corrected, the ENDF/B-VII.1 fission yields provide spectral shapes in close agreement with the JEFF3.1 fission yields. In 2012, the agreement obtained was at the level of 10% with respect to the integral beta spectra measured at ILL and the number of nuclei requiring new TAGS measurements was considered as achievable. Lists of priority for new TAGS measurements were established first by the Nantes group [67] (which triggered our first experimental campaign devoted to reactor antineutrinos in 2009), then by the BNL team [100] and eventually a table based on the Nantes summation method was published in the frame of TAGS consultant meetings organized by the Nuclear Data Section of the IAEA [103]. A portion of the table from [103] is shown in Table 6, with the mea- Table 6. List of nuclides identified by the IAEA TAGS Consultants that should be measured using the total absorption technique to improve the predictions of the reactor antineutrino spectra. These nuclides are of relevance for conventional reactors based on 235 U and 239 Pu nuclear fuels. The list contains 34 nuclides [103]. Relevance (Rel.) stands for the priority of the measurement. Isotopes marked with asterisks show the measurements performed by our collaboration, m stands for metastable or isomeric state.

Isotope
Rel  [24,62] and Figure 9) placed after the JYFLTRAP Penning trap of the IGISOL facility [51]. Only 92,93 Rb were in the top ten of the nuclei contributing significantly to the reactor antineutrino spectrum. 92 Rb itself contributes 16% of the antineutrino spectrum emitted by a pressurized water reactor (PWR) between 5 and 8 MeV. Its contributions to the 235 U and 239 Pu antineutrino spectra are 32% and 25.7% in the 6 to 7 MeV bin and 34 and 33% in the 7 to 8 MeV bin. In 2009, the ground state (GS) to GS beta intensity of this decay was set to 56% in ENSDF. This value was revised in 2014 to 95.2%. A maximum of 87%±2.5% was deduced from our TAGS data, having a large impact on the antineutrino spectra [67]. This value was confirmed by the Oak Ridge measurements [101]. The 92 Rb case is worth noting because it is not a case suffering from Pandemonium, but its GS to GS beta branch was underestimated in former evaluations. In the analysis of this nucleus, the sensitivity of the reconstructed spectrum (and thus of the χ 2 obtained ) to the value of the GS-GS branch was very high, because of the large penetration of the electrons in the TAGS. The quoted uncertainties were obtained by varying the input parameters entering into the analysis, such as the calibration parameters, the thickness of the beta detector, the level density, the normalisation of the backgrounds, etc.. The beta decay data for 92 Rb that were used in the previous summation calculations [99] were those from Teng-blad et al. The impact of replacing these data with the new TAGS results amounts to 4.5% for 235 U, 3.5% for 239 Pu, 2% for 241 Pu, and 1.5% for 238 U. A similar impact was found on the summation model developed by Sonzogni et al. [100] but a much larger impact (more than 25% in 235 U) was found on another model in which no Pandemonium-free data were included [104].
Though not motivated by neutrino physics, 86−88 Br and 91,94 Rb were measured in the same experiment. 86−88 Br and 91 Rb were not in the priority list of [103]. The main motivation for those cases was the study of moderate beta delayed neutron emitters using complementary techniques and the study of the decay used as normalization in the measurements by Rudstam et al [45] already mentioned in Section 3. These TAGS measurements confirmed the Pandemonium problem in the existing data. 86,87 Br and 91 Rb did not show a large impact on the antineutrino spectra [12,24,62]. On the other hand, 94 Rb, ranked as priority 2 in the table from [103], and 88 Br exhibited a quite large impact on the spectra. This was verified with two different summation calculations [24]. In one of the models, the new TAGS data replaced high resolution spectroscopy data, and thus the observed impact was typical of a correction of the Pandemonium effect i.e. a decrease of the high energy part of the aggregate antineutrino energy spectrum. The impact reached 4% in 235 U and 239 Pu in the case of 94 Rb and was more modest as regards 88 Br with a 2-3% decrease in the latter actinides between 8 and 9 MeV. The latter range explains the reason why 88 Br did not belong to the priority list of [103] that was established for a contribution to the PWR antineutrino spectrum larger than 1% between 3 and 8 MeV. In the second model the TAGS data replaced data measured by Tengblad et al. that were considered in the model because they were assumed to be Pandemonium-free. The replacement of 87 Br had little impact, that of 94 Rb led to a 3% decrease at 8 MeV but that of 88 Br brought a 7% increase between 8 and 9 MeV, with a cancellation of the last two effects below 8 MeV.
The cumulative impact of the TAGS beta intensities measured with the Rocinante detector at Jyväskylä on the antineutrino energy spectra generated after the thermal fissions of 235 U, 239 Pu and 241 Pu, and fast fission of 238 U are presented in Figure 17 with respect to that built with the most recent evaluation decay databases JEFF3.3 [105] and ENDF/B-VIII.0 [106] for the same nuclei and containing only the TAGS data from [7,17]. The decrease of the two plutonium spectra above 1.5 MeV is remarkable, reaching 8%. The impact on the two uranium isotopes amounts to about 2% and 3.8% in the 3 to 4 MeV range in 235 U and 238 U respectively. These results were provided by Dr. M. Estienne [107].
In our 2014 experimental campaign, we were almost exclusively focussed on nuclei of importance for the prediction of the reactor antineutrino spectrum and for decay heat calculations using the DTAS detector [30]. Twentythree isotopes were measured, among them many isomers which require the separation power of the Jyväskylä Penning trap. An illustration of the experimental challenge is given by the case of the Niobium isomers 100,100m,102,102m Nb.  [103]. 100 Nb is among the main contributors to the antineutrino flux in the region of the shape distortion, along with 92 Rb, 96 Y and 142 Cs. The results showed that the high resolution measurements for 100,100m Nb and 102gs Nb were affected by the Pandemonium effect, while the beta-intensity distribution for 102m Nb was determined for the first time [13]. The impact of these measurements on the summation calculations was evaluated (see Figure 18) and resulted in a large impact between 3 and 7 MeV, with a strong decrease of the spectrum peaked at 4.5 MeV and a strong increase peaked at 6.5 MeV, in the region of the shape distortion. In the calculation, the TAGS data replaced high resolution spectroscopy data extracted from JEFF3.3 [105] and ENDF/B-VIII.0 [106]. As a result, the discrepancy between the summation antineutrino spectra including these data and the experimental reactor antineutrino spectra is diminished in the region of the shape distortion, though the distortion has not vanished completely [13]. The results presented in this last figure were also provided by Dr. M. Estienne [107].
In parallel to the TAGS campaigns, the reactor antineutrino experiments have published their near detector measurement of the emitted antineutrino flux and spectrum from PWRs. In 2017, the Daya Bay experiment could measure the reactor antineutrino flux associated with various fuel compositions [108], and found a flux coming from 239 Pu fission in agreement with the prediction of the Huber-Mueller model, while the flux associated with 235 U fission exhibited a deficit of 7% thus nearly explaining by itself the reactor anomaly. This new result does not favour the idea of oscillation into sterile neutrinos, as it would affect equally the antineutrinos arising from both fuels. It would rather confirm the hypothesis of an additional systematic uncertainty associated with the 235 U energy spectrum. These recent findings reinforced the necessity of an alternative approach to the converted spectra which could be brought by the use of nuclear data. It was thus timely to perform a comparison of the summation method spectra with the Daya Bay results. The first comparison was performed in [109] showing a discrepancy with the measured antineutrino flux of only 3.5%, nearly twice as small as that with the Huber-Mueller model. We have performed an update of our summation model in [110] using the above-mentioned Pandemoniumfree datasets improved by the TAGS campaigns of the last decade, the most recent evaluated databases (JEFF3.3, ENDF/B-VIII.0) and updated gross theory spectra [111] for the unknown beta decay properties.
After folding with the Inverse Beta Decay (IBD) cross section [112] the summation spectrum built with the actinides spectra weighted with the fission fractions published by Daya Bay, the resulting detected spectrum was compared with that of Daya Bay [113] and that built using the Huber-Mueller model. In Figure 19, the top panel Energy (MeV) Fig. 19. Comparison of the summation antineutrino spectrum obtained using the fission fractions published in [113] and all the TAGS data quoted in this section, with the experimental spectrum from reference [113]. Gr ee nw oo d Fig. 20. Comparison of the Inverse Beta Decay (IBD) yield computed with the summation antineutrino spectrum obtained using the fission fractions published in [108] for 239 Pu and using all the TAGS data quoted in this section (included successively), with the experimental IBD yield from [108]. Greenwood represents the result of the summation model [110] when only the TAGS results of Greenwood et al. [17] are included (for more details of the model see [110] [24,62], and SM-2018 contains the impact of 100,100m,102,102 Nb decays [13] (always considered in addition to the earlier version of the summation model).
Reprinted figure with permission from [110], Copyright (2019) by the American Physical Society.
shows the ratio of the Daya Bay antineutrino spectrum over that computed with the Huber-Mueller model (open diamonds with error bars) superposed with the ratio of Daya Bay over the summation method spectrum including the TAGS results from our first campaign (dashed line) and over the summation method spectrum including the TAGS results from both campaigns (plain line). The normalisation of the summation method spectrum is clearly in better agreement with the experimental data than that of the Huber-Mueller spectrum (closser to 1), which is considered nowadays the model reference. The inclusion of the TAGS measurements of the Niobium isomers [13] has further improved the shape agreement especially in the energy region of the shape distortion. The bottom panel of Fig. 19 shows the ratio of the summation method spectrum with that of Huber-Mueller. Here again the latest TAGS measurements have flattened the ratio which shows a rather good shape agreement, though located below one at about 95-96%. Still the summation method spectrum does not reproduce the shape distortion seen by the reactor antineutrino experiments at PWRs. Figure 20 summarizes the detected antineutrino flux (called IBD yield) as a function of the fission fraction of 239 Pu obtained with the summation method spectra depending on the TAGS results included in the calculation. The explicit labels of the lines describe the TAGS results introduced one after the other. It is noticable that the inclusion of more TAGS data systematically decrease the detected antineutrino flux to end with an 1.9% discrepancy with the Daya Bay measured IBD yield. This is a consequence of the correction of the Pandemonium effect in nuclear databases and emphasizes the importance of the TAGS method and measurements. More details are given in [110] in which the individual IBD yields associated with 235 U, 239 Pu, 241 Pu, and 238 U obtained with the summation model are also compared with the Daya Bay results. The agreement is good in general for all four isotopes. This is at variance with the Huber-Mueller model for which a large discrepancy is observed in the case of 235 U while the three other cases are in very good agreement with the experiment.
In our 2014 TAGS campaign at Jyväskylä devoted to reactor antineutrinos and decay heat, 96,96m Y, 140,142 Cs, 138 I, 137 I, 95 Rb, 95 Sr, 103 Mo and 103 Tc were measured as well. These future TAGS results may complete the picture that starts to be drawn of the reactor antineutrino energy spectra. In parallel to the nuclear physics effort, reactor antineutrino experiments at short baseline from research reactors start to release their first results. Up to now, only the NEOS and Neutrino-4 collaborations have released a combined result which signals the presence of an oscillation [114]. Neither the STEREO [85] or the PROSPECT [86] experimental results confirm this oscillation signal. Furthermore, the PROSPECT experiment has released their first spectral measurement of the antineutrino energy spectrum from Highly Enriched Fuel (HEU) which is equivalent to a pure spectrum from 235 U. It is remarkable that their shape-only result does not show such a pronounced shape distortion as the large experiments at PWRs. It thus excludes the idea that the shape anomaly arises solely from 235 U. It is worth mentioning also that PROSPECT is the first detector using 6 Li instead of Gd to capture the neutron formed in the IBD process since the Bugey experiment [115], which did not see a shape anomaly either. Lately Double Chooz has also released their fourth measurement of the θ 13 mixing angle obtained by cumulating neutron captures on Gd and the H and C contained in the target and gamma catcher volumes [116]. They observe a shape distortion which could be fitted either with a single or a double Gaussian with a slope. One of their conclusions is that the one sigma envelope for today's prediction appears insufficient to accommodate the mismatch between data and model for both rate and shape. A better understanding of the origin of model deviations remains critical and the role of nuclear data is definitely crucial at the time at which experiments are being set up to measure the mass hierarchy of neutrinos. In a recent publication [117], the global reactor antineutrino data set was re-analyzed using three reactor antineutrino flux predictions, the Huber-Mueller model, the summation method of [110] and the model of [92] which includes a theoretical calculation of the form factors for the first forbidden transitions. Relative to the traditional Huber-Mueller predictions, the two new calculations result in diverging evidence for a sterile neutrino when total IBD rate measurements are considered. The summation calculation of [110] decreases the significance from 2.3 to 0.95 σ, while that of [92] increases the significance to 2.8 σ. However, the spectral anomaly is robust with any of the flux models. The accurate determination of the reactor antineutrino spectra is also mandatory to monitor future reactors with antineutrino detection. Predictions for innovative reactor designs and fuels can only be obtained through the use of nuclear data and the summation method. The fine structures present in the antineutrino energy spectra induced by the end-points of the individual beta branches from the fission products [104,118] could provide a benchmark for nuclear data and an insight of what is going on inside a reactor. But they also degrade the sensitivity of detectors such as JUNO [119] by mimicking a periodic oscillation pattern. These fine structures may be directly observed by the JUNO-TAO one-ton detector that will be located a few metres away from a PWR core [120]. In parallel, an experimental confirmation of the observed first hint of coherent elastic neutrino-nucleus scattering [121] would definitely open new possibilities for neutrino applications.

Nuclear structure
In an article about the application of TAGS, it would be remiss of us to neglect its application to the study of nuclear structure. It is not our intention here to tell the reader all about beta decay. That can be found in text books (see for example [122,123,124,125]). Instead what we want to do is provide a few examples that show how useful TAGS can be in testing nuclear models that contribute to our understanding of the underlying structure of the atomic nucleus, and in particular present some cases recently studied in the framework of reactor applications.
In Section 2 we already mentioned the 100 Tc case, of relevance for double beta decay studies. Here we will focus on two nuclear structure applications of TAGS that are of significance, namely the study of of the quenching of Gamow-Teller transitions and the study of nuclear shapes.
One essential concept in beta decay, important in the examples that follow, is the beta strength function (see [14]), a quantity that can be deduced from experiment. It is defined as: where I β (E) is the beta intensity to the level at excitation E in the daughter nucleus, f is the statistical rate Fermi integral which depends on the energy available in the decay (Q β − E) and T 1/2 is the half-life of the decay. S β (E) is in practical terms the reciprocal of the f t values given conventionally in the literature. We will concentrate on our contribution to the determination of the strength by measuring the beta feeding in a reliable way, which is the main subject of this article. The other two quantities, namely Q β − E or the T 1/2 are obtained from measurements dedicated to this purpose. In the following we will focus on allowed Gamow Teller transitions, since allowed Fermi transitions are normally concentrated in a single state, and not affected by P andemonium very much. The experimental beta strength is related to the experimental B(GT) in the case of Gamow Teller transitions through the following equation: where g A and g V are the axial-vector and vector coupling constants. The B(GT ) as defined above can be related to the transition probability calculated theoretically between the parent state and the states populated in the daughter defined as follows: where σ and t represent the spin and isospin operators acting on the individual nucleons and Ψ i and Ψ f the initial and final nuclear states.
In consequence, a comparison of the B(GT ) theo should reproduce the B(GT ) determined in the experiment. Actually, the quality of the comparison reflects the goodness of the nuclear model in describing the involved nuclear states. In addition there is a model independent rule, called the Ikeda sum rule, that tells us how much strength we should observe. Curiously enough, the strength obtained experimentally seems to be systematically lower than theory. This is called the Gamow Teller quenching problem and it has been discussed for more than four decades and is not yet fully resolved (see for instance [126,127] and more recently [128] and references therein). The discussion of this mismatch between theory and experiment involves theoretical as well as experimental arguments. One main difficulty is that the full strength, in the case of GT transitions, is normally concentrated in a resonance at relatively high excitation energy, for example in the range of 8-15 MeV for A∼100, that is difficult to access in beta decay. In consequence, charge-exchange reactions with hadronic probes, such as (p,n), ( 3 He,t) or (t, 3 He), have been used to measure the full strength. Extracting the GT strength from these probes is more complicated than extracting it from beta decay. Among other reasons, this is because it relates to the penetrability of the hadronic probe in the nucleus and because of the difficulty of selecting the GT process in a clean way. Beta decay, however, has its own difficulties. The principal one, as mentioned above, is to know how much of the strength lies within the Q-value window, and the other difficulty is to be sure that we measure all the strength inside the Q β window, because of the Pandemonium effect. It is in overcoming this second difficulty that the TAGS technique has had an impact in tackling this problem.
In order to avoid the first problem, one can choose cases where most of the strength is expected to be located at relatively low energy, inside the Q β window, and this happens, in principle, in the beta decay of nuclei south east of 100 Sn on the Nuclear Chart and in the rare-earth nuclei above the spherical nucleus 146 Gd. These cases, even although they do not have a direct relation with the cases of relevance for reactor applications presented in this article, can provide information on the necessary corrections that are required for a proper theoretical description of the beta decay process. The reason, in both the 100 Sn and above 146 Gd regions, is that there is only one main component in the GT strength on the β + side, namely πg 9/2 →νg 7/2 in the 100 Sn region and πh 11/2 →νh 9/2 in the rare-earth case. All other proton occupied orbitals have no empty neutron orbital partner. Unfortunately, the expected B(GT) cannot be directly compared with the Ikeda sum rule in these cases because this rule involves the B(GT) values for both the β + and the β − decays and here only the former can be measured. So, it has to be compared with theory. A relatively simple but realistic calculation of the expected beta strength on the β + side was carried out by Towner [129] for decays in both of these regions of the Segre Chart. In this work [129] a hindrance factor h is defined as the ratio between the summed GT strength from theory and experiment. Initially he adopted the extreme single particle approach (s.p), namely considering only the two pairs of orbitals, πg 9/2 -νg 7/2 and πh 11/2 -νh 9/2 . He then made a series of corrections to this approach taking into account pairing, core polarization and higher-order effects and then looked at how hindered the corrected theoretical strength would be in comparison with the extreme s.p picture. This result defines a theoretical hindrance factor that can be compared later with the hindrance obtained from the ratio of the extreme single particle approximation and experiment. The theoretical hindrance was calculated for the range of cases from n=1 to n=10 active protons in the g 9/2 orbital in the 100 Sn region and n=1 to n=12 in the h 11/2 orbital in the rare earths.
A series of experiments were carried out at GSI (Germany) with heavy ion beams from the UNILAC at energies slightly above the Coulomb barrier on the appropriate targets to study relevant beta decays in these regions of interest. At these energies, the reaction was dominated by the fusion evaporation channels which are fewer than in the fission examples described in previous sections. Consequently, the separation achieved with the relatively simple Mass SEParator (MSEP) [130], provides clean enough samples to perform the experiments. The GSI TAS [131] was built and briefly used at the Berkeley SuperHILAC and after the accelerator was closed, it was installed at the GSI MSEP. This spectrometer enjoyed two advantages over the INEL Idaho TAGS [65]. Firstly, the crystal was more than twice the size, secondly it included a small cooled high purity Ge detector for X-Ray detection. The first improvement was important for a better response of the spectrometer to the absorption of the gamma cascades, and the second to clean the EC (Electron Capture) component of the decays further and hence to obtain a very good Z separation. The results can be seen in references [132] ( 100 In), [133] ( 97 Ag), [134] ( 148 Dy) [135,136] ( 150 Ho) and [137] ( 148 Tb(2 − and 9 + isomers), and 152 Tm(2 − and 9 + isomers)). For the very special case of 100 Sn, where the expectations are that all the strength is concentrated in a a single 1 + state in the daughter, we will consider the results of Hinke et al. [138] and the more recent work by Lubos et al. [139] measured with Ge detectors. To study this case further, a TAGS experiment has been proposed and partially carried out at RIKEN [140] and is currently under analysis [141]. In general, one observes that the calculations [129] reproduce the observed hindrance factor fairly well for the 100 Sn region. One should mention here that this case is simple enough that it has been calculated from first principles reproducing the experimental value also fairly well [128]. At the time of this last work [128] only the [138] results were available and included in the comparison, but the more recent results of [139] show even better agreement with these calculations, similar in quality to the agreement with the systematic extrapolation of the strength from [142]. In contrast, in the rare-earth region higher hindrance factors have been observed experimentally compared to Towner's calculations. One should note that, somewhat surprisingly, the GT resonance is clearly observed in this case. Moreover, the tail below the resonance that has been observed experimentally and discussed at length in the Charge Exchange reaction experiments, is clearly seen in these beta decay experiments for the first time. As an example, we show in Figure 21 the case of the decay of the 150 Ho 2 − isomer to 150 Dy [135], with the TAGS measurement represented by the spectrum under the black line. A similar result was obtained in the 152 Tm case reported in [137]. In this article, it was suggested that the missing strength, or in other words the explanation of the disagreement with Towner, is probably located in that part of the tail of the resonance which is cut off by the Q β window.
The same Figure 21 can also be used to illustrate the importance of the TAGS experiments and the limitations of the Ge detectors in terms of observing beta feeding at high excitation energy that was discussed in the introduction. We see the B(GT) distribution measured with the Cluster cube [135] and with the GSI TAS [136]. The Cluster cube was an array of six Euroball Ge cluster detectors in compact geometry. It was equivalent to forty two individual Ge detectors and had an efficiency of 10.2(5) % at a gamma-ray energy of 1332 keV. The figure shows clearly the importance of both types of measurement. In the Ge measurements 1064 gamma rays were identified and the coincidences between detectors allowed the construction of a decay scheme with 295 levels in 150 Dy [135]. Figure 21 shows in blue the B(GT) strength to each of these levels as deduced from the beta feedings in the decay scheme. Inspection of the total absorption spectrum reveals that the Ge array loses sensitivity as a function of excitation energy in the daughter nucleus when compared with the TAGS and our ability to determine the feeding or beta intensity distribution diminishes. Once converted into B(GT) strength we conclude that we lose 50% of the total strength compared with the TAGS measurement. Moreover, the tail of the resonance, one of the foci of interest in the experiments discussed above cannot be seen with the Ge detector array.
In summary, returning to the discussion of the missing B(GT) strength in beta decay, even though the experiments explained above are probably the best cases to study the Gamow Teller quenching, they rely very much on comparison with theoretical calculations which, in general, cannot locate with sufficient precision whether the calculated strength is inside the accessible beta window or not. However, these measurements have demonstrated that the tail below the Gamow Teller resonance exists, and this observation is free of background ambiguities. Moreover, in the future, when either the calculations are accurate enough to tell us how much of the strength should be located within the accessible beta-window, or, alternatively when we are able to perform charge exchange reactions using radioactive beams, we have here very reliable measurements of that part of the spectrum that lies below the Q β energy. This can be used for normalisation purposes as well as for control of the reaction mechanism.
Another application of TAGS measurements, first pioneered at CERN-ISOLDE with the Lucrecia TAGS, relates to the shapes of nuclear ground states. The concept of nuclear shape is deceptively simple. In practice it is difficult to measure. The measurements with TAGS are based on a theoretical idea put forward by Hamamoto and Zhang [143], that was developed further by Sarriguren et al. [144] and Petrovici et al. [145]. They showed that the beta strength distribution for transitions to excited states in the daughter nucleus depends on the shape assumed for the ground state of the decaying nucleus. Intuitively one can see why this might be so if we look at the ordering of deformed single particle orbits on a Nilsson diagram. The levels on the prolate and oblate sides are in different order  Fig. 21. Comparison of the beta strength deduced from the high resolution measurement of the beta decay of the 150 Ho 2 − isomer using the cluster cube setup [135] (shown in blue) with the strength obtained from a total absorption measurement [136] (black) . The measurements were performed at the Mass Separator at GSI. Reprinted figure with permission from [135], Copyright (2003) by the American Physical Society. and thus filling them up to the Fermi level to determine the ground state configuration of a particular nucleus involves different single particle contributions. Their beta decay strength distributions will also be different since they are dictated by angular momentum and parity selection rules as well as the overlap of the wavefunctions of the states involved. The calculations by Hamamoto and Zhang, Sarriguren et al. and Petrovici et al., are of course rather more sophisticated than this simple picture which is just used for understanding the underlying physics phenomena.
In particular regions of the nuclide chart, nuclei can have several minima in the potential energy surface with different shapes for the ground state. The calculated B(GT) distributions for each of these states with a defined deformation are quite different in some but not all cases (see for example [144]). Where they are different, the experimental B(GT) distribution measured with TAGS can then be compared with the theoretical distributions and the ground state shape inferred. A number of studies of this kind ( [5], [6], [9], [10], [11]), have been carried out for nuclei with A∼ 80 and A∼190. A summary of these activities at ISOLDE can be found in [16].
This method has also been applied for some of the cases studied at IGISOL. See for example 100,102 Zr and 103 Mo [40,146]. In Figure 22 [40,147] we present a comparison of the deduced strength for the decay of 103 Mo with Quasiparticle Random Phase Approximation (QRPA) calculations performed by P. Sarriguren [148]. From this comparison a preference for an oblate shape in the ground state of 103 Mo can be inferred. In the calculations a quenching factor of ( g A g V ) ef f = 0.77( g A g V ) has been applied, which is equivalent to a hindrance factor of 1.69 with respect to the QRPA calculations used in the comparison.
Another example of the importance of TAGS measurements in testing models is provided by 105 Mo. The decay of 105 Mo was calculated using the FRDM-QRPA model [149,150]. The best theoretical description of this decay  [40,147] in comparison with QRPA calculations assuming prolate or oblate deformations in the ground state of 103 Mo [148]. In the model a quenching factor of ( was obtained assuming a ground state deformation of 2 =-0.31 for 105 Mo. The experimental half-life of this decay is 35.6 s, and this value can be better reproduced if first forbidden transitions are included in the model calculation (T theo 1/2 =30.3 s), but in that case, the experimental beta distribution is not reproduced so well. This can be seen in Figure 23 where the experimental feeding distribution is compared with the theoretically deduced distributions with and without first forbidden transitions. A better reproduction of the beta distribution by theory is obtained if no first forbidden component is included in the model. But in that case the experimental half-life is not so nicely reproduced (T theo 1/2 =150 s). This clearly shows a limitation of the performance of this model in a region which is dominated by shape effects and where triaxiality can play a role. QRPA calculations assume that both the parent and the daughter have the same deformation, which might not always be applicable in regions where shape transitions are common. This example shows the relevance of having in addition to the experimental half-life the possibility of comparing the theoretical strength (or the deduced theoretical feeding) with reliable experimental data, like that provided by TAGS measurements. Based on the description of the half-life only, we might conclude it is necesary to introduce the first forbidden component for the description of this decay, which does not reproduce well the experimental beta feeding. The relevance of this kind of model validation will be further discussed in the next Section in relation to astrophysical applications.

Astrophysical applications
TAGS measurements are important also in the context of nuclear astrophysics. We select here some examples related to the astrophysical r process. As mentioned in the previous Section, some of those examples will show a strong interrelation with nuclear structure studies.
The r process is driven by a huge instantaneous flux of neutrons that creates by successive neutron captures  very neutron-rich nuclei, up to the heaviest ones, that then beta decay towards stability. About half of the observed abundance of elements heavier than Fe in the Universe is synthesized in this way. The identification of the astrophysical site where the process occurs is the subject of very active investigations. Core Collapse Supernovae were the classical favoured site in spite of persistent difficulties met when trying to reproduce observations with calculations. On the other hand Neutron Star Mergers have recently become [151] a confirmed site for heavy element formation after the first observation of the gravitational waves generated and the analysis of the subsequent electromagnetic radiation. Much remains still to be done in order to understand the role of both scenarios combining astrophysical observations and calculations that require nuclear physics input (see [152] for a recent review). Some of the key input parameters in r process calculations are the decay properties of neutron-rich nuclei, more specifically half lives (T 1/2 ) and beta-delayed neutron emission probabilities (P n ) that control the nucleosynthesis flow. For such exotic nuclei the neutron sep-aration energy (S n ) becomes smaller than the decay Q β value and neutron emission from populated neutron unbound states occurs. In spite of current efforts at the most advanced radioactive beam facilities to determine this information experimentally [153], most of the nuclei involved cannot be accessed in the laboratory and need theoretical estimates. The key point here is that both quantities are derived from the beta strength distribution S β (E) (see also Equation 7) 1 (11) Equation 11 above includes the competition between neutron (Γ n ) and gamma (Γ γ ) emission. It should be noted that models often assume that neutron emission prevails always and the competition is ignored.
The quality of beta strength calculations is usually asserted by global comparisons with measured half-lives and to a lesser extent with measured neutron emission probabilities. However it is found that different theoretical models with comparable quality predict quite different T 1/2 and P n (see for example [154]). This clearly indicates that the quality assessment based on these integral quantities (Equations 10 and 11) is not good enough. This comes as no surprise since several (theoretical) strength distributions can lead to the same half-life for a particular nucleus. But the underlying nuclear structure can then predict very different numbers for neighboring nuclei. Comparing TAGS measurements of the beta strength with different models then becomes the only reliable validation method (see the 105 Mo case discussed in the previous Section). Moreover it can give us hints on how to improve the nuclear structure calculations.
Another example is related to the determination of neutron capture (n, γ) cross-sections for very exotic neutronrich nuclei, that also controls the nucleosynthesis flow in the r process. These are even more difficult to determine experimentally because of the need to prepare suitable targets. Direct measurements will require very imaginative techniques thus current efforts concentrate on indirect methods [155]. Theoretical estimates are based on the statistical Hauser-Feshbach model [156] that uses average quantities: nuclear level densities, photon strength functions and neutron transmission coefficients. These are parameterized using data measured mostly close to stability and consequently there is considerable uncertainty on the values needed in r process calculations.
We have proposed a way to obtain experimental constraints on the quantities that intervene in the Hauser-Feshbach estimate for very exotic neutron rich nuclei [12,24]. It is based on the analogy between radiative neutron capture reactions and the process of beta-delayed neutron emission. The former depends mostly on Γ γ and weakly on Γ n and the latter can provide the ratio Γ γ /Γ n provided that we are able to measure the (expected weak) gamma emission from neutron unbound states. This is where the sensitivity of the TAGS technique comes into play. The main advantage of the method is that the measurements can be extended into regions quite far from stability.
In [12,24] the gamma-neutron competition was studied for the 87,88 Br and 94 Rb decays and more recently in [25] for 95 Rb and 137 I . The results are summarized in Table  7, which shows P γ , the gamma emission probability above S n defined by analogy with P n (also shown). Observation of Table 7 reveals that in most of the cases P γ is large, even larger than P n . The large P γ for 137 I is confirmed by the TAGS measurement of [35]. The reason for this surprising result is to be found in the nuclear structure of the nuclei in the decay chain. A large mismatch between spin and parity of unbound states in the daughter nucleus and the available states in the final nucleus means that neutron emission is hindered by the centrifugal barrier. Other measurements have also found large P γ values in the decay of 70 Co [157] and 83 Ga [158] and different nuclear structure effects were invoked to explain it. This notable result warns us about the neglect of gamma-neutron competition in theoretical estimates of P n (see also [159]), but does not tell us about the statistical parameters of the Hauser-Feshbach model. The most interesting case from that point of view is the decay of 94 Rb where level densities in the daughter nucleus are large and neutron emission is not hindered by angular momentum mismatch, making it a good test case for the statistical model. In this case gamma emission above S n represents only 5% of neutron emission but even so it is more than one order-of-magnitude larger than Hauser-Feshbach calculations using standard statistical parameters. This is a challenging outcome. However in order to translate a constraint in Γ γ /Γ n into a constraint on (n, γ) cross-section we need additional information. If we follow the assumption that extrapolating far from stability nucleon optical parameters (that determine Γ n ) is more reliable than extrapolating photon strength functions (that determine Γ γ ) then this result would indicate one order-of-magnitude increase in the calculated capture cross-section. Clearly more investigations are required and new TAGS measurements on suitable isotopes are planned.
A different method to obtain constraints on (n, γ) crosssections for unstable nuclei using TAGS measurements has been proposed by the NSCL group [37]. It was already mentioned in Section 2 in connection with the extraction of the branching ratio matrix, which is the first step of the Oslo method [38]. The goal of the Oslo method is to obtain the shape of the nuclear level density and photon strength function from the branching ratio matrix (in their terminology: primary gamma ray intensities) and it was originally applied to nuclear reaction experiments. Going from a relative quantity (branching ratios) to absolute quantities (photon strength functions and nuclear level densities) requires the use of normalization parameters coming from external sources. In the case of beta decay TAGS measurements away from the stability these are systematics, extrapolations or theory. Thus the impact of this method is mainly that of the shape of the photon strength function.
A closely related topic and also very interconnected with nuclear structure is the potential provided by beta decay in relation to the study of collective phenomena. Beta decay could constitute a new means to investigate the presence and maybe some of the properties of lowlying collective modes, such as pygmy dipole modes predicted to appear at lower energies as nuclei become more neutron rich. Collective modes are of crucial importance in nuclear structure as they reflect the ability of the nucleons to move coherently and provide insights into the properties of the nuclear force. The study of collective modes puts constraints on theoretical models as well. They are also the only observables that we can study on earth providing access to the intrinsic properties of nuclear matter, entering into the modelling of astrophysics phenomena like supernovae or neutron stars. Pygmy dipole resonances (PDR) could be the consequence of the appearance of neutron skins in medium to heavy neutron-rich nuclei. The PDR might deliver information on neutron-star properties [160]. Important information on the equation of state (EOS) of neutron-rich matter via strength-neutron-skin thickness correlation could be obtained [161].
The presence of low-lying PDR could influence processes of nucleosynthesis, especially (n,γ), (γ,n) reactions playing an important role in the r-process [162] as mentioned earlier. Several questions remain unanswered about the collective modes when nuclei become more exotic. One limitation up to now has been the low intensity of the accessible exotic beams which limits the possible studies using standard nuclear or electromagnetic probes. In this context beta decay constitutes a new probe for lowlying collective modes. Further away from stability, as the energy window opened by beta decay increases, the energy of the pygmy modes decreases, allowing their excitation through the Gamow-Teller operator when the spin and parity conservation conditions are fulfilled. Beta decay then offers new possibilities to study systematically the presence of low-lying collective modes with the existing exotic beam intensities. Our collaboration was first to propose an experiment on this topic [163]. Later on, the theoretical demonstration was provided by two models [158,164]. The quasi-particle model of [164] predicts that other components of the collective mode are excited through beta decay than those excited by the usual nuclear and electromagnetic probes. In particular, beta decay would feed preferentially two-particle two-hole components of the collective mode, being thus complementary to nuclear reactions. In the experimental results of [164] and [158], high resolution setups with a relatively small detection efficiency were used and the data may suffer from the Pandemonium effect. The TAGS technique, using modern segmented spectrometers, seems to be very well adapted to tackle this problem, especially to evidence high energy gamma-rays feeding the daughter ground state or the first excited state. In parallel, ways to obtain experimental evidence of the collectivity of the states fed by beta decay, Table 7. Pγ obtained from our measurements [24,25] in comparison with the Pn values of the decays. Pγ is defined as the gamma emission probability above the Sn value (in analogy to Pn). The values are given in % (see the text for more details). which would not rely on theoretical predictions, should also be investigated.

Summary, future and conclusions
In this article we have presented a review of the impact of our total absorption studies of beta decays that are relevant for reactor applications. The measurements presented have been performed at the IGISOL facility of the University of Jyväskylä employing the high isotopic purity beams provided by the JYFL Penning Trap. These measurements are not only relevant for the decay heat predictions and for the predictions of the reactor neutrino from reactors, but also provide results of interest for nuclear structure and astrophysics. In particular they offer the possibility of testing nuclear models in a more stringent way and can provide additional information for the estimation of (n,γ) cross sections of astrophysical interest for cases not directly accessible using reactions. Considerable progress has been made, but the ultimate goal of the work presented in this article has not yet been reached. From the comparisons of the measured decay heat with the predictions of summation calculations, it is clear that there is still work to be done, in particular for the 235 U fuel. The situation is similar in relation to the prediction of the antineutrino spectrum in reactors, where the remaining discrepancies still require to measurements of a number of decays. Our collaboration is still working on these subjects and has approved proposals to continue our studies at the IGISOL IV facility in order to measure new decays that are important in the next relevant order. In this publication we have concentrated mainly on the discussion of results obtained by our collaboration, but other groups are also involved in similar research programmes at other facilities that provide experimental results relevant to the topics discussed here (see for example [35,101,102,157]). The upgrade and advent of a new generation of radioactive beam facilities like FRIB (Michigan, USA), RIBF (RIKEN, Japan), FAIR (Germany), Spiral2 (France), etc. extends the possibilities of TAGS measurements to more exotic domains than those offered by the present and future ISOL facilities. These measurements represent new challenges concerning the purity of the beams and require the development of detectors adapted to the experimental conditions of such a facilities. From those facilities new and exciting results will appear in the near future.