Development of a Reference Database for Beta-Delayed Neutron Emission
P. Dimitriou,1, ∗ I. Dillmann,2, 3 B. Singh,4 V. Piksaikin,5 K.P. Rykaczewski,6 J.L. Tain,7 A. Algora,7 K. Banerjee,8 I.N. Borzov,9, 10 D. Cano-Ott,11 S. Chiba,12 M. Fallot,13 D. Foligno,14 R. Grzywacz,15, 6 X. Huang,16 T. Marketin,17 F. Minato,18 G. Mukherjee,8 B.C. Rasco,19, 6, 15, 20 A. Sonzogni,21 M. Verpelli,1 A. Egorov,5 M. Estienne,13 L. Giot,13 D. Gremyachkin,5 M. Madurga,15 E.A. McCutchan,21 E. Mendoza,11 K.V. Mitrofanov,5 M. Narbonne,13 P. Romojaro,11 A. Sanchez-Caballero,11 and N.D. Scielzo22 1NAPC-Nuclear Data Section, International Atomic Energy Agency, A-1400 Vienna, Austria 2TRIUMF, Vancouver, British Columbia V6T 2A3, Canada 3Department of Physics and Astronomy, University of Victoria, Victoria, British Columbia V8P 5C2, Canada 4Department of Physics and Astronomy, McMaster University, Hamilton, Ontario L8S 4M1, Canada 5State Scientific Center of Russian Federation, Institute of Physics and Power Engineering, 249033 Obninsk, Russian Federation 6Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 7IFIC, CSIC-Universitat de Valencia, 46071 Valencia, Spain 8Variable Energy Cyclotron Centre, Kolkata 700064, India 9National Research Centre Kurchatov Institute, Department of Nuclear Astrophysics, 1 Akademika Kurchatova pl., Moscow 123182, Russian Federation 10Joint Institute for Nuclear Research, Bogolubov Laboratory of Theoretical Physics, 6 Joliot-Curie, 141980 Dubna, Moscow region, Russian Federation 11Centro de Investigaciones Energ´eticas, Medioambientales y Tecnol´ogicas (CIEMAT), Avenida Complutense 40, Madrid 28040, Spain 12Tokyo Institute of Technology, 2-12-1-N1-9, Ookayama, Meguru-ku, Japan 13Subatech, CNRS/in2p3, Univ. of Nantes, IMTA, 44307 Nantes, France 14CEA, DEN, DER/SPRC/LEPh Cadarache, F-13108 Saint Paul-Lez-Durance, France 15Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, USA 16China Nuclear Data Center, China Institute of Atomic Energy, Beijing 102413, China 17Department of Physics, Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia 18Nuclear Data Center, Nuclear Science and Engineering Center, Japan Atomic Energy Agency, Tokai-mura Ibaraki-ken, Japan 19JINPA, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 20Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA 21National Nuclear Data Center, Brookhaven National Laboratory, Brookhaven, Upton NY, USA 22Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA (Dated: February 3, 2021) Beta-delayed neutron emission is important for nuclear structure and astrophysics as well as for reactor applications. Significant advances in nuclear experimental techniques in the past two decades have led to a wealth of new measurements that remain to be incorporated in the databases. We report on a coordinated effort to compile and evaluate all the available β-delayed neutron emis- sion data. The different measurement techniques have been assessed and the data have been com- pared with semi-microscopic and microscopic-macroscopic models. The new microscopic database has been tested against aggregate total delayed neutron yields, time-dependent group parameters in 6-and 8-group re-presentation, and aggregate delayed neutron spectra. New recommendations of macroscopic delayed-neutron data for fissile materials of interest to applications are also presented. The new Reference Database for Beta-Delayed Neutron Emission Data is available online at: http://www-nds.iaea.org/beta-delayed-neutron/database.html. arXiv:2102.01165v1 [nucl-ex] 1 Feb 2021 CONTENTS 4. Neutron counting per fission (“fiss,n”)8 5. Neutron and γ counting (“n,γ”)8 I. INTRODUCTION2 6. Double γ counting (“γ,γ”)8 7. Ion and neutron counting (“ion,n”)9 II. MICROSCOPIC DATA: Methods and 8. Double ion counting (“ion,ion”)9 measurements4 B. Delayed neutron spectra9 A. Delayed neutron emission probabilities4 1. 3He spectrometers9 1. Beta-neutron coincidences (“β–n”)6 2. Gaseous proton recoil spectrometers 10 2. Neutron and β counting (“n,β”)7 3. Neutron energy spectroscopy with A 3. Relative neutron counting (“P1n Z”)8 time-of-flight spectrometers 10 C. New methods (for Pn values and energy spectra) 11 1. Total absorption γ-ray spectroscopy 11 ∗ Corresponding author: [email protected] Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
2. Ion-recoil methods 12 neutron data 45 3. Optical Time Projection Chamber 12 C. Estimation of energy dependence of total 4. Measurement in a storage ring 12 delayed neutron yields 45 5. Ion counting using a MR-TOF 13 D. Outlook for future measurements 13 VI. SUMMATION CALCULATIONS 50 A. A numerical (T1/2,Pn) data file for practical III. COMPILATION AND EVALUATION 14 applications 50 A. Compilation and evaluation methodology 14 B. Total delayed neutron yields 51 1. Basic summation calculations 51 B. New evaluated (T1/2,Pn) data 14 1. Z = 2 – 28 region 14 2. Time-dependent calculations 54 2. Z = 29 – 57 region 17 C. Time-dependent delayed-neutron integral 3. Z = 58 – 87 region 17 spectra 57 D. Time-dependent delayed-neutron C. P1n Standards 17 parameters 60 D. Systematics for P1n values 17 1. Kratz-Herrmann formula 18 2. McCutchan systematics 18 VII. INTEGRAL CALCULATIONS 67 3. Improved systematics along Z chains 18 A. Comparison with integral experiments 67 4. Systematics based on the Effective B. Impact on reactor calculations 69 Density Model 18 VIII. SYSTEMATICS OF MACROSCOPIC DATA: E. Delayed neutron spectra for individual Time-dependent parameters 71 precursors 19 A. New approach for estimation of temporary 1. Digitization of delayed neutron (DN) parameters for unmeasured nuclides 71 energy spectra 19 B. The (A /Z) · 92 systematics of the average 2. Assessment of spectra as a reference for c half-life 72 calibration purposes 20 C. Systematics (Ac − 2 · Z) of relative abundances in the 8-group model 72 IV. THEORY OF MICROSCOPIC DATA 20 A. Self-consistent QRPA-based models 22 IX. RECOMMENDED DATA IN 6- and 8-GROUP 1. DF+CQRPA 22 MODELS 72 2. Relativistic Density Functional+QRPA 22 B. Comparison with new evaluated (T1/2,Pn) X. REFERENCE DATABASE 79 data for important β-delayed neutron A. Microscopic data 80 emitters 23 B. Macroscopic data 80 1. Light fission products 23 2. Heavy fission products 24 XI. SUMMARY AND CONCLUSIONS 80 C. Comparison for isotopes relevant to nuclear A. Microscopic experiments 84 astrophysics 24 B. Compilation and evaluation of microscopic 1. β-decay properties in the 78Ni region 27 data 84 2. β-decay properties in the 132Sn region 28 C. Microscopic theory 85 D. Global comparisons of theoretical results 30 D. Summation and time-dependent E. Delayed neutron spectra calculations 35 calculations 86 E. Integral experiments and reactor V. MACROSCOPIC DATA: Methods and calculations 86 Measurements 38 F. Macroscopic measurements and A. Methods 38 systematics 86 1. Total delayed neutron yields and G. New recommendations of group constants 87 time-dependent parameters 38 2. Delayed neutron integral spectra 41 Acknowledgments 87 B. New measurements and compilation 41 1. Energy dependence of the relative References 88 abundances and half-lives of delayed neutrons 42 2. Energy dependence of total delayed I. INTRODUCTION neutron yield 42 3. Integral energy spectra of delayed Neutron-rich nuclei can emit neutrons after β- neutrons from fission of 235U by thermal decay when their decay Q-value is larger than the neutrons 44 (one/two/three. . . ) neutron separation energy: Qβ > 4. ALDEN: new measurements of delayed Sxn. This decay mode is called “β-delayed (βn) neu-
2 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. tron emission” and was discovered in 1939 by Roberts et ities and spectra in the late 1970s. A large number of al. [1,2], shortly after the discovery of fission by Meit- data was published in the following decades and subse- ner, Hahn, and Strassmann in 1938 [3]. Immediately af- quently compiled and evaluated by Rudstam and collab- ter this, Bohr and Wheeler [4] improved their liquid drop orators [15, 16, 20, 21] and Pfeiffer et al. [22]. model to describe these new decay mechanisms. Values from these compilations and evaluations have The β-delayed two- and three-neutron emission (β2n sometimes been incorporated in ENSDF [23], together and β3n) process was discovered 40 years later in 11Li with independent recommendations by the respec- [5,6], and in 1988 the (so far) only β4n emitter, 17B, tive evaluator. Consequently, these values have also has been investigated and a tentative branching ratio of been transferred to evaluated data libraries such as P4n = 0.4(3) has been measured [7]. ENDF/B [24, 25], JENDL [26, 27] and JEFF [28, 29] ‘Delayed’ in this context means, that the neutron is decay-data libraries, as well as in nuclear astrophysics emitted with the β-decay half-life of the precursor AZ, applications. Other independent evaluations have also ranging from a few milliseconds for the most neutron- been performed to address specific needs in applica- rich isotopes up to 80(18) s for 211Tl, the longest- tions, such as the highly-specific evaluations of decay lived β-delayed neutron (βn) precursor observed to date. data for short- to very short-lived radioisotopes, includ- These delayed neutrons have to be distinguished from the ing β-delayed neutron emitters, which were performed prompt neutrons evaporated immediately (in the order of from 1996 to 1998 by Nichols [30] for the UK-based data 10−14 s) after a fission event from a neutron-rich nucleus file. These evaluations were subsequently adopted in the and usually form only a small fraction of the total neu- JEFF decay-data libraries [29]. trons emitted from fission (≈ 1%). Nevertheless, they The first comprehensive compilation and evaluation of play an important role in controlling the dynamic rate of β-delayed neutron emission spectra was published by fission in a chain-reacting assembly and hence in reactor Brady [31] in an effort to combine measured spectra with kinetics and the safe operation of a reactor. statistical model calculations to complete the experimen- This important role of delayed neutrons in determining tal data that was sometimes scattered or limited in en- the reactor kinetics and reactivity was recognized early ergy range. The resulting spectra were adopted in the on and led to several measurements of delayed neutron ENDF/B libraries and were further supplemented with groups [8–12]. In 1954, an extensive delayed neutron model calculations by Kawano et al. [32]. program was initiated at Los Alamos National Labora- Until 2011 only sporadic efforts have been made to tory. One major outcome of this program was the 6- collect data of ‘microscopic’ quantities related to the β- group fitting parameters by Keepin and co-workers (see delayed neutron emission, namely half-lives and neutron- Sect.VA1,[13]). In a least-squares fit of the delayed neu- branching ratios of very neutron-rich nuclei to reveal nu- tron emission curves of 235U, they found that by grouping clear structure properties for nuclei far from stability. the βn emitters in six groups – according to their half- This data has also been used directly as nuclear physics lives – they could reproduce the experimental data. input for astrophysical reaction network calculations for In 1990, Subgroup 6 was established by the NEA- the ”rapid neutron capture” (r) process in explosive as- OECD “Working Party on International Nuclear Data trophysical scenarios like core collapse supernovae and Evaluation Co-operation” (WPEC-SG6) [14] to assess the binary neutron star mergers, or indirectly to benchmark differences between the calculated and measured values of and improve theoretical predictions which are needed for the reactivity scale obtained from reactor kinetics. The extrapolating models towards the yet inaccessible regions observed discrepancies were largely attributable to the of the nuclide chart (“Terra Incognita”). uncertainties in the delayed neutron data used in these Recent experimental efforts relied largely on the exist- calculations. Their effort to address this included interna- ing generation of radioactive beam (RIB) facilities but tional benchmark measurements of the effective delayed with the next generation of these facilities becoming op- neutron fraction, made on fast critical assemblies, which erational within this decade, a plethora of new nuclides aimed to provide high-quality experimental information will become accessible. The vast majority of these new for 235U, 238U, and 239Pu. nuclides on the neutron-rich side will be βn-emitters. One of the outcomes of this working group was the re- According to the Atomic Mass Evaluation in 2016 finement of the Keepin parameters by the laboratories in (AME2016) [33], 651 out of the 3435 known nuclei are Los Alamos and Obninsk from a 6-group fit to an 8-group βn emitters, meaning that their Qβ1n is (within mass fit. This was achieved by separating the three longest- uncertainties) positive. Theoretical models estimate that lived delayed neutron precursors (87Br, 137I, and 88Br) about 4000 additional nuclei should exist between the into single groups. proton- and neutron dripline, most of them are on the Independent of the effort to determine integral delayed neutron-rich side and will be βn emitters. neutron yields for applications, there have also been sev- In 2011, 72 years after the discovery of β-delayed neu- eral attempts to measure, compile and evaluate β-delayed trons , only 216 βn emitters had measured neutron- neutron data of individual precursors. Groups in Swe- branching ratios: 81 nuclei in the lighter mass region den [15, 16], USA [17, 18], and Germany [19] became (Z < 28), 134 nuclei in the fission regions between active in measurements of β-delayed neutron probabil- Z = 29 − 57, and a single (unconfirmed) measurement
3 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
210 of the P1n value of Tl [34] in the whole mass region delayed neutron data are described and compared with above Z > 57. recommended values in Sect.VI. The impact of the new Currently (as of August 2020) 653 nuclei between 8He β-delayed neutron data on integral reactor calculations and 237Ac have been identified as potential β-delayed one- is discussed in Sect.VII. Finally, the systematics of the neutron emitters (Qβ1n > 0 keV). However, as of August measured macroscopic data are discussed in Sect.VIII 2020 only 306 of those nuclei have measured β-delayed and new recommendations for 6- and 8-group kinetic one-neutron branching ratios (P1n): 114 in the light mass parameters for major and minor actinides are given in region Z < 28 [35], 183 in the fission fragment region Sect.IX. (Z = 29 − 57), and 9 in the region above Z > 57 [36]. The structure of the reference database is detailed in A total of 28 of these measurements provided only an Sect.X and a summary of the main conclusions of this upper limit for the P1n value, and another six provided work is provided in Sect.XI. only lower limits. All the microscopic and macroscopic data presented As far as delayed multi-neutron emission is concerned, herein, are available online at http://www-nds.iaea. 300 potential β2n-emitters have been identified, along- org/beta-delayed-neutron/database.html [43]. side 138 β3n-emitters, and 58 β4n-emitters. The number of measured multi-neutron branching ratios is however much smaller: only 31 P2n values have been measured II. MICROSCOPIC DATA: METHODS AND MEASUREMENTS (23 for Z < 28 and eight for Z > 28) and four P3n values (11Li,17B, 19B, and 23N). One single tentative measure- 17 ment for the P4n value of B has been reported. A. Delayed neutron emission probabilities One of the largest and most focused recent efforts to fill this large gap between identified and measured β- The neutron-branching ratio of a nucleus into the β1n delayed neutron branching ratios is the ”Beta-delayed decay channel is defined as neutrons at RIKEN” (BRIKEN) project [37]. Since 2016 this international collaboration has been measuring half- N1n lives and neutron-branching ratios of almost all accessible P1n = . (1) βn emitters. The collaboration will conclude its experi- Ndecays mental campaigns in 2021 and is expected to add in the near future > 250 new β-delayed neutron emitters to where N1n stands for the number of decays (Ndecays) the existing list. The data from the experimental runs that emit one neutron. Similar expressions hold for any in 2016–19 is presently being analyzed, and a small frac- other β-delayed multi-neutron decay channel. This quan- tion has already been published and included in the latest tity can be measured in various ways with different tech- Z > 28 evaluation [36]. An overview of the BRIKEN and niques, either directly by counting the number of emit- future campaigns can be found in Sect.IID. ted neutrons, or indirectly by counting γ-rays or other decay products. Special care has to be taken that the This rapidly developing landscape of measured number of counted events are free of any contaminations, βn emitters and respective application-driven needs for i.e. that the contributions from background activities are reliable βn data necessitates both, a thorough update of subtracted and the random noise is properly corrected the existing βn-data compilations and a review of the for. recommended group parameters. In response to these ur- Eight methods for the determination of P values were gent data needs, the International Atomic Energy Agency 1n extracted from the previous Rudstam evaluation [21] and (IAEA) coordinated an international effort to produce a were reviewed in the IAEA CRP summary reports [39– ‘Reference Database of Beta-Delayed Neutron Emission 41] for their strengths and weaknesses and, in particular, Data’ [39–42]. potential sources of systematic errors and uncertainties. The results of this coordinated effort are reported in this paper. The methods for measurement of microscopic These methods were complemented by newer methods. β-delayed neutron data are described in Sect.II, fol- Note that only a few of these methods are also suited for lowed by the compilation and evaluation efforts (includ- the measurement of multi-neutron emission probabilities. Throughout this article the following notation will be ing new recommendations of standards for P1n measure- A ments) and systematic descriptions in Sect.III. Compar- used: precursor ( Z) or mother nucleus stands for the β- isons of the β-delayed neutron data (half-lives, β-delayed decaying nucleus; β-decay daughter or intermediate nu- cleus is the product from the β decay of the precursor neutron -branching ratios and spectra) with global theo- A retical approaches are presented in Sect.IV. ( Z+1), and β-delayed neutron daughter stands for the final nucleus with A−1Z+1 for the case of one-neutron Macroscopic (integral) delayed neutron yields and A−2 spectra have been used to verify and validate the mi- emission, Z+1 for the case of two-neutron emission, croscopic data. Measurement methods, new data, and etc. compilation efforts are presented in Sect.V. Summation The majority of Pxn values in the literature have been and time-dependent calculations of total delayed neutron determined using methods that involve neutron count- yields and spectra based on the evaluated microscopic β- ing using moderated neutron detectors. These type of counters combine a medium acting as a neutron energy
4 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 1. Chart of nuclides with color-coded Qβn window (in keV, retrieved from NuDat 2.7 [38]). moderator with a detector sensitive to thermal neutrons counters with some additional 3He tubes from RIKEN to maximize detection efficiency. Typical moderators are [37, 51]. made of dense hydrogenated material, like high-density The detection efficiency of these type of detectors ex- polyethylene that has replaced the paraffin wax of earlier hibits a dependency on neutron energy which can be a designs. The first counters used BF3-filled proportional source of systematic uncertainty in the determination of tubes as a thermal neutron detector, but nowadays the Pxn values (see discussion in IIA1). The effect is re- more expensive 3He is preferred because of its advantages lated to the moderation process and depends on the ge- [44] (see also Sect.IIB). This detector can be considered ometry of the counter (number and spatial distribution as an evolution of the neutron long counter [45], where of neutron detectors). Thus the detector design can be very high detection efficiencies are achieved by arranging optimized to minimize the dependency of the efficiency several tubes in a 4π geometry around the source. Apart on neutron energy. An example of this is the LOENIE from the large detection efficiency, another distinct ad- counter [47]: With only 18 3He tubes a very small (17%) vantage of this type of detector is the virtual absence of but flat efficiency up to 2 MeV is achieved. a low-energy threshold. Nowadays the general trend is to minimize the energy Recent examples of this kind of setup applied to the dependency of the efficiency curve while still maximiz- measurement of Pxn are NERO (now HABANERO) ing the overall detection efficiency. This can be achieved at the National Superconduction Laboratory (NSCL) by increasing the number of tubes considerably, as has at Michigan State University (USA) [46], LOENIE at been done for example in the case of the BRIKEN setup Lohengrin-Institute Laue Langevin (France) [47], 3Hen which could use up to 166 3He tubes. The optimization of at HRIBF-Oak Ridge National Laboratory (USA) [48], the geometry incorporating many tubes is a complicated BELEN at the GSI Helmholtz Center for Heavy Ion Re- problem that can be solved by combining Monte Carlo search in Darmstadt (Germany) and at the heavy ion ac- simulations with a parametric approach [37]. In this way celerator laboratory of the University of Jyv¨askyl¨a(Fin- a configuration with an efficiency in excess of 75% below land) [49], and TETRA at ALTO - Institute de Physique 1 MeV was obtained that drops (only) to an efficiency of Nucleaire in Orsay (IPNO, France) [50]. The latest setup 54% at 5 MeV. is BRIKEN at RIKEN Nishina Center in Wako (Japan) The efficiency calibration of these instruments requires which is a temporary merger of the 3Hen and BELEN the use of mono-energetic isotropic neutron sources span-
5 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. ning the range of neutron energies of interest to well above (almost) constant neutron detection efficiency for ener- 5 MeV for exotic nuclei [52]. Photo-neutron (γ, n) and α- gies up to 1 MeV, and then the curve drops slowly. These induced reaction sources (like 241Am/9Be or 24xCm/12C), effects might introduce systematic errors if the calibrating as well as fission sources have been employed [53] but isotopes have very different β- or neutron-energy spectra these are difficult to prepare or have a very broad energy [49]. Examples of efficiency curves for a typical β- and spectrum. Beam-induced (p, n) and (α, n) reactions have moderated neutron detector used in experiments at the also been used [46] although it seems difficult to reach IGISOL facility in Jyv¨askyl¨aare shown in Fig.2[54]. It the required accuracy. can be seen that the highest neutron energies are con- In experiments at radioactive beam facilities β-delayed nected to low-energy β’s with the lowest detection effi- neutron emitters with well know P1n values and neutron ciency due to the detector thresholds, and vice versa. spectra are used for the efficiency calibration of the setup. An additional source of uncertainty is that the energy For sources with a broad neutron energy distribution the dependence (slope) of the neutron efficiency curve is often measured neutron detection efficiency is averaged over ignored. However, if the neutron energy distribution of the respective neutron energy range. However it should the calibrant isotope is very different to the one of the be kept in mind that this is an approximation as different isotope to be measured, this will induce systematic effects neutron spectra – detected with different average efficien- which need to be corrected for. For moderated neutron cies – can have the same average energy. Thus, for these detectors the “average detection efficiency” depends on sources the measurements should be combined with de- the neutron energy spectra of each isotope. Fig.3 shows tailed Monte Carlo simulations, although the accuracy of the neutron spectra of the βn emitters 87Br (green line), the latter in complex geometries is still under discussion. 95Rb (blue line), and 137I (red line) as an example. For The moderation process introduces a delay between comparison, two simulated neutron detection efficiency neutron emission and detection that reaches few hun- curves for two state-of-the-art neutron detectors, BELEN dreds of µs. This poses a problem for conventional trig- [54] and BRIKEN [37], are plotted as dashed lines. gered event-based data acquisition systems. The use of The correct description for the neutron-branching ratio an event gate of this length can lead to a sizable loss of should therefore include detection efficiencies hi which data due to the acquisition dead time. This problem is are averaged over the respective energy ranges: overcome by modern event-less self-triggered data acqui- N sition systems, as in the case of the most recent neutron β β1n P1n = 0 (3) detection setups 3Hen, BELEN, TETRA, and BRIKEN. hβ · ni Nβ In the following, the eight different methods for obtain- 0 In this equation β is the mean β efficiency. hβ · ni is ing P1n values as well as their advantages and disadvan- tages are discussed. the product of the averaged β efficiency above the neutron separation energy weighted by the respective β intensity spectrum and the averaged neutron efficiency weighted 0 1. Beta-neutron coincidences (“β–n”) by the neutron spectrum [49, 54]. If β = hβi, then the simplified Eq.2 is recovered. Even with this simplification, the method requires the In this method β’s and neutrons are counted in coinci- knowledge of the neutron spectra of the respective isotope dence and compared with the number of β counts. Thus to be measured (see Sect.IIB). Since these spectra are the efficiency of the β detector is greatly reduced, while not available for all isotopes, one can use as a guidance the neutron efficiency in absolute terms is needed. This the spectra of calibration isotopes like the ones shown in method was labelled “n/β” in Ref. [21] and has been re- Fig.3 and listed in SectionsIIIC andIIIE, selected to named here to “β–n” to account for the proper sequence cover the expected neutron energy range of the isotopes of emissions of the particles. of interest. For measuring the one-neutron emission probability When analyzing decay curves the signals from all con- (P ), the respective equation is 1n tributing decays need to be considered. This requires a proper knowledge of the different half-lives of all involved 1 Nβ1n nuclei and possible isomeric states. In turn, if a half-life P1n = , (2) in a decay chain is not know, this method can also be used 1n Nβ to extract the unknown half-lives of βn emitters from β- where 1n is the one-neutron detection efficiency, Nβ1n and βn-decay curves. is the number of time-correlated β-1n events, and Nβ is Extensions of the “β–n” method to include multi- the amount of detected β’s. The main assumption in Eq. neutron emission increase the complexity as the neutron- 2 is that the detector efficiency for β decays is constant for emission multiplicity (x) increases. Since the neutron- all energies and thus cancels out. However, β detectors detection efficiency is < 100%, a βxn decay channel will have a lower threshold of typically 50–150 keV, and the contribute to the counts Nβyn with y ≤ x [6] which respective efficiency curve shows a steep drop below ≈1–2 must be taken into account. It is also important to in- MeV (see Fig.2). clude corrections for accidental coincidences with back- In contrast, a moderated 3He-neutron detector has an ground neutrons, in particular when the background rate
6 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 2. Examples of efficiency curves for a β- and a neutron detector [54] in comparison with an example neutron spectra. Note that the highest neutron energies are connected to β’s with the lowest detection efficiency due to the detector thresholds.
neutrons, since the coincidence relations will increase the sensitivity to the neutrons of interest. It is a reliable method if all pitfalls as described above are carefully con- sidered and accounted for. The main sources of system- atic uncertainty are the β and neutron detection efficien- cies that have to be simulated carefully.
2. Neutron and β counting (“n,β”)
In this method the β’s and neutrons are counted si- multaneously but separately (not in coincidence). This method was labelled “n–β” in Ref. [21]. The P1n value can be extracted from the ratio of counted events cor- rected by the ratio of detection efficiencies:
FIG. 3. Examples of neutron spectra for 87Br (green), 95Rb (blue), and 137I (red) in comparison with two (simulated) ef- β N1n P1n = . (4) ficiency curves for state-of-the art neutron detectors (dashed n Nβ and dotted line) used in experiments at the IGISOL facility in Jyv¨askyl¨a[54] and at RIKEN [37]. The “n,β” method requires that either both efficien- cies are determined separately or that their ratio is ob- tained from calibration nuclei. All systematic errors on β and neutron detection efficiencies related to the differ- is large and the Pxn value is small. The importance of a ent intensity distributions and the shape of the efficiency proper background correction has been discussed recently curves discussed in the previous method apply also here in Refs. [51, 54, 55]. As in the P1n case (Eq.3), the cor- and have to be included properly. The β detection needs rect expressions must include variable, isotope- and decay to be checked for potential sources of background, e.g. channel-dependent β- and neutron-detection efficiencies. conversion electrons or changing noise in the detectors. In summary, the “β–n” method is a generic and widely This method is also quite generic and applicable to the employed method that is well-suited for the measurement measurement of multiple neutron emission. The main of single- and multi-neutron-branching ratios . It is the practical limitation is the neutron signal-to-background method of choice if the rate of neutrons from the respec- ratio that can be overcome by the use of β-n coincidences tive decay is small compared to the rate of background (the “β–n” method).
7 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
A 3. Relative neutron counting (“P1n Z”) 5. Neutron and γ counting (“n,γ”)
This is a neutron-only counting method and requires This method was called “γ AZ” in [21] and is renamed that in the decay chain either a descendant or the parent to “n,γ” method to account for the separate counting of A Z is a βn-emitter with well known P1n value which be- neutrons from the precursor and γ’s from the daughter comes the normalization. In addition, all half-lives have AZ decay (AZ+1 in the nomenclature in this paper). It to be known including that of the measured precursor should not be confused with the “n–γ” method in [21] which could be determined independently from β activ- which actually designates a “γ,γ” method (see below). ity curves. It is based on the analysis of the measured Because of the usually much longer half-live of the daugh- neutron activity curve to disentangle the different contri- ter it is convenient to measure first the neutron activity butions. The method has been applied only in a few cases of the sample, and then the γ-activity in a separate ded- after chemical or mass separation, thus activity curves icated setup [58]. abs must be corrected for contaminants. In its simplest form The method requires that the absolute γ-intensity Iγ it assumes that the neutron energy spectrum (and thus of one or more transitions in the decay of the daughter the neutron detection efficiency) for parent and descen- is known. Corrections for growth and decay of activities, dants are identical but eventually corrections as discussed which require a knowledge of half-lives, must be applied before should be applied. before calculating the P1n:
abs P γ · I N 4. Neutron counting per fission (“fiss,n”) 1n = γ n (6) 1 − P1n n Nγ This method was labelled “fiss” in Ref. [21]. Histori- cally this is the oldest method used to measure neutron The γ-ray detection efficiency γ is usually well known, emission probabilities of individual precursors [56]. It is thus the main issue is the availability and accuracy of based on the determination of the number of β-delayed absolute γ intensities per decay. If the daughter decay n activity is also co-produced and present in the sample, neutrons per fission of a given precursor YA,Z with a suit- able neutron detector. This number is compared then to this needs to be accounted for. The method can be used also with γ-rays from subsequent descendants in the β- the fission yield YA,Z of that precursor to obtain the P1n value: decay chain.
n YA,Z 6. Double γ counting (“γ,γ”) P1n = (5) YA,Z This method involves counting of γ-rays emitted by de- Typically fission of 235U at thermal neutron energies scendants in both, the β-decay and the βxn-decay chains. has been used in these measurements. The main source The γ activities can be labeled by their masses, which in of systematic error is the accuracy of the fission yields the case of β1n emitters are A and A-1, respectively. The available at the time of the measurement. For example, in method requires knowledge of the absolute γ-intensities the evaluation of Rudstam [21] the results of earlier stud- per decay for both transitions used [59]. After correct- ies were updated using the fission yields YA from Wahl ing the γ counts by the growth and decay of activities, [57] and the respective charge fraction PZ . However, new the possible branching-out/branching-in from/to the re- mass- and charge-dependent fission yields YA,Z are now spective decay chains and the presence of contaminant available, e.g. from evaluated libraries like ENDF/B- activities in the sample, the P1n can be deduced from: VIII.0 and JEFF-3.1.1. Thus, a renormalization of the older data to the new fission yields standards is advis- P A · Iabs,A N A−1 able. 1n = γ γ γ (7) A−1 abs,A−1 A This method was applied mostly to chemically sepa- 1 − P1n γ · Iγ Nγ rated samples [21], a procedure that requires a certain processing time (in the order of 1 s) and is not suited The main advantage of this technique lies in the rel- for very short-lived β-delayed neutron emitters. It is af- ative simplicity of detecting γ rays as compared to de- fected by systematic uncertainties of the chemical separa- tecting neutrons, but the problem is shifted to the de- tion efficiency as well as the neutron efficiency. For older termination or knowledge of the absolute γ intensities experimental data, the lack of separation of the different per decay. The main drawback can be the low signal- isotopes and the presence of contaminants could have af- to-background ratio requiring that the limited peak effi- fected the determination of the P1n value. Experience ciency of γ-ray detectors is compensated by high count- shows that the reliability of this method can be poor and ing statistics. Tagging with β particles can be used to that should be taken into account when comparing with improve the sensitivity but at the cost of introducing a other methods. possible dependence on the β detection efficiency.
8 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
If β-decaying isomers are present with half-lives similar advances in instrumentation. It is based on counting di- to the ground state, γ rays can provide a uniquely selec- rectly the number of parent nuclei and the number of tive tool to discriminate them, and measure correctly the daughter nuclei produced after neutron emission, circum- Pxn values of isomeric and ground states. venting the direct neutron detection. Fragmentation facilities offer new opportunities for this The method includes several different techniques like technique, as descendants of the nucleus of interest AZ in detecting recoils in a Beta-decay Paul trap (BPT) [64], the β and βxn decay chains can be independently pro- observing recoil tracks with an optical time-projection duced and transmitted with the beam. In this case, ab- chamber (OTPC) [65], or separating and counting ions solute γ intensities for suitable transitions can be deter- with a heavy-ion storage ring [66] or a Multi-Reflection mined in the same experiment using the number of im- Time-Of-Flight (MR-TOF) spectrometer [67]. With the planted ions for normalization [60]. However, the possible first two methods it is also possible to extract neutron existence of isomers should be investigated, and the ac- energy spectra for the one-neutron emitters via the en- tivity should be corrected according to the correlation ergy of the recoil, but not in the case of multi-neutron time. Since ions with different Z are implanted at dif- emitters. More details can be found in Sect.IIC. ferent depths in the respective implantation detector, β efficiency corrections might also be needed. The “γ,γ” method is suitable for the measurement of B. Delayed neutron spectra multi-neutron emission branches. Most of the available data on the spectra of β-delayed neutrons have been obtained using three different tech- 7. Ion and neutron counting (“ion,n”) niques: 3He spectrometers, gaseous proton recoil spec- trometers, and time-of-flight spectrometers. In the fol- This corresponds to the “ion” method as defined in [21]. lowing the principles, instruments, and analysis methods It is based on the direct measurement of the number of of the three techniques are described. parent nuclei rather than the β’s emitted. It has been applied at mass separators by counting the number of se- lected ions. For example, ions can be counted by implan- 1. 3He spectrometers tation on an electron multiplier that is surrounded by the neutron detector [61]. In other arrangements two mea- This method is based on the detection of fast neutrons surements are needed, one for counting the ions, where through the reaction n +3 He → p +3 H, which has a the beam is deflected to the electron multiplier, and one positive Q-value of 764 keV. In a 3He filled proportional for counting neutrons where the beam is implanted in detector the reaction leads to the appearance of a peak a collector in the middle of the neutron detector [62]. in the spectrum for mono-energetic neutrons. The neu- In this case, additional systematic uncertainties must be tron energy E can be determined from the displacement considered. The neutron emission probability is deter- n of the observed peak with respect to the thermal neu- mined as tron peak. For a continuum neutron energy distribution, the measured spectrum must be deconvoluted using the ion N1n energy-dependent neutron detector response. P1n = . (8) Following the introduction of high-pressure gridded n Nion ionization chambers [68] with good energy resolution Here Nion is the number of ions after correction for (10−20 keV) and efficiency, this method has been used to contaminants and background and ion is the ion detec- measure β-delayed neutron spectra of several fission sys- tion efficiency. The determination of the ion detection tems and individual precursors [31]. One disadvantage efficiency of the electron multiplier can be an issue as of this type of detector is the fast drop in efficiency with shown in [61]. neutron energy: whereas at thermal energies the reaction More recently (see [63] as an example), non-destructive cross section is as large as 5316 b it drops to 0.8 b at ion counting with microchannel plates before implanta- 1 MeV. This imposes a practical limit in the maximum tion has been applied to post-accelerated mass-separated measurable energy of about 3 MeV. In addition, the large ion beams. In this example however, γ counting of de- peak due to the thermal-epithermal neutron background scendant decays was used instead of neutron counting to limits the minimum measurable energy to about 100 keV obtain the P1n value (“γ,ion” method). or more. The accurate determination of the detector response is an essential step in this method. This is typically 8. Double ion counting (“ion,ion”) achieved using mono-energetic neutrons from reactions like 7Li(p, n)7Be. The shape of the measured response The “ion,ion” method is possibly the only new method departs from that of a simple peak due to several ef- for the measurement of β-delayed neutron emission prob- fects. Apart from electronic and microphonic noise, the abilities that was developed in recent years, thanks to response includes the effect of γ-ray interactions (includ-
9 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. ing pulse pile-up), a peak due to thermal neutron back- 3. Neutron energy spectroscopy with time-of-flight ground, the effect of neutron-induced recoils in the gas, spectrometers the detector-wall effect, and ballistic effects (incomplete charge integration, see [69] for a detailed discussion). In- The time-of-flight (TOF) methods involve detection formation on signal risetime can be used to eliminate of neutrons using the β-neutron delayed coincidence some of these artifacts [70], at the cost of lowering the method. In radioactive beam experiments, the start sig- efficiency. nal for the neutron time-of-flight measurement is given Nowadays, digital electronics allow a more detailed by a β detector trigger, and the stop signal is provided pulse shape analysis [71]or the use of the full response by a neutron detector, which is typically a fast plastic function [72]. After parametrization of the detector re- or liquid scintillator. The required time resolution for sponse as a function of En, several suitable algorithms the TOF measurements is in the nanosecond range. The can be used to unfold the response from the measured TOF distance and timing resolution determine the total spectrum [20, 69]. The evaluation of systematic uncer- resolution [79]. tainties in this method is not a simple task. Information The key advantages of TOF arrays are the relatively about their magnitude can be obtained by comparison of high intrinsic efficiency for high-energy neutrons, the abil- results obtained by different groups and detector systems ity to measure neutrons without recording their full en- [20, 73]. This issue is discussed in Sect.IIIE. ergy, and the relatively well-understood interaction cross- sections with the scintillator which enables reliable simu- lations. The disadvantage of the TOF technique is the limitation of the achievable resolution for high-energy 2 2. Gaseous proton recoil spectrometers neutrons (> 3 MeV) due to the E ∼ 1/T OF dependence in the typical geometry and susceptibility of the detector This method is based on the detection of protons that response to neutron scattering. Neutron-γ discrimina- are elastically scattered by neutrons. The detectors used tion can be achieved by using liquid scintillators. The are either methane- or hydrogen-filled proportional coun- high cost of the detector material prevented the use of ters at high pressure in which the gas acts at the same solid organic crystals such as stilbene and para-terphenyl time as the target for the reaction and as detecting me- for making large TOF arrays. Recently developed plastic dia. These detectors are insensitive to thermal neutrons scintillators with neutron-γ discrimination capability [80] and compared to 3He spectrometers have a better en- are now available. ergy resolution below ∼ 200 keV. Because of this they Until recently the most successful TOF arrays used for are considered to be superior in this energy range [31]. β-delayed neutron spectroscopy were TONNERRE [81], They have, however, a lower efficiency which limits the which was used in experiments at CERN and at GANIL practical measurable energy range up to about 1 MeV. (France), and the TOF array developed at the NSCL At typical βn energies, struck protons are scattered (USA) [82]. Both were used to measure light-mass nu- clei with distinct neutron transitions [83–85]. The heavi- isotropically and the ideal detector response is a step 51−53 function. Thus the neutron energy distribution is ob- est studied isotopes were K[86]. Both TOF arrays tained by spectrum differentiation [74] rather than by share a similar design with about 1 m flight path, and deconvolution [75]. Fourier filtering can be applied to long and curved plastic scintillator bars. reduce the fluctuations of the extracted distribution [76]. The currently operational TOF arrays for decay stud- However the true detector response departs from the ies are VANDLE [87, 88] (used at various facilities in ideal response and a number of corrections must be ap- the USA, at CERN, and at RIKEN), DESCANT [89, 90] plied [17]. These type of detectors are sensitive to γ- at TRIUMF/ISAC in Vancouver (Canada), and MON- rays that contaminate the low-energy part of the spec- STER which has been built for FAIR and already used trum (< 100 keV), but γ signals are efficiently elimi- at ISOLDE and at the IGISOL facility in Jyv¨asky¨a[91]. nated using pulse shape discrimination. Different dis- VANDLE is a TOF array using a fully digital data ac- crimination methods have been applied [17, 77] and the quisition system [92]. The focus is on a high detection energy-dependent impact on neutron counting must be efficiency with the use of straight bars made of EJ200 evaluated. In addition, the response is altered by the re- plastic scintillators (Eljen Technology, USA). A digital coils of heavier nuclei present in the gas mixture, by wall- distributed trigger enables a low-energy neutron thresh- and end effects and by position-dependent variations of old down to 100 keV. the signal amplitude. The PSNS code [78] allows the MONSTER is a TOF array based on BC501A (St. inclusion of all these corrections and is commonly em- Gobain) and EJ301 (Eljen Technology) cylindrical liq- ployed to analyze the data. In general the corrections are uid scintillators. Both BC501 and EJ301 offer excel- found to be small and sometimes instead of correcting lent neutron/γ-ray pulse shape discrimination capabili- the data they are used to evaluate the size of the related ties, which allow to strongly suppress the competing γ- systematic uncertainties [17]. A comparison of the results ray background and β-neutron detector coincidences due obtained with this technique and with 3He spectrometers to cosmic rays. MONSTER is operated with a fully dig- is presented in Sect.IIIE. ital data acquisition system based on 14 bits resolution
10 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. and 1 Gsample/s digitisers. A total absorption spectrometer is a γ-ray calorimeter DESCANT is a 70-element array of deuterated liquid- made with a large volume of inorganic scintillation ma- scintillator detectors (BC-537, St. Gobain) that can be terial surrounding the source with a ∼ 4π solid-angle used with both the TIGRESS and GRIFFIN γ-ray spec- coverage. Such a detector is capable of absorbing the full trometers at TRIUMF-ISAC for reaction studies and the energy of the γ-cascade following the decay, thus identi- investigation of β-delayed neutron emitters [89, 90]. DES- fying the excitation energy of the populated state. How- CANT forms a close-packed array that replaces the for- ever, due to detector inefficiency and β penetration, this ward “lampshade” of 4 HPGe clover detectors in each of is not true for all events, and a spectrum deconvolution the GRIFFIN and TIGRESS γ-ray spectrometers, pro- procedure is needed to recover the information. viding high detection efficiency for neutrons in the range Due to the large detector size, neutrons emitted after of ≈ 100 keV to 10 MeV. β decay have a large probability to interact with the scin- The decay of 17N, a well-defined standard for neutron- tillation material through elastic collisions, inelastic scat- branching ratio measurements (see TableI), can be used tering, and radiative capture. The latter two processes for energy- and detector-response calibrations in TOF ex- generate γ-rays that are easily detected. In principle, this periments. The energy of the three main neutron lines are is a source of contamination of the spectrum that needs well-known from 3He ion chamber measurements [93–97] to be determined and subtracted [107]. and appear fully resolved in TOF measurements [81, 98]. The neutron reaction cross sections, specific to the scin- Alternatively, if a calibration with neutron transitions tillation material, determine the amount of background of energies higher than 2 MeV is desired, the use of 15B signals. Large differences are observed in (n, γ) cross sec- or 49K is possible. However, there are no high-resolution tions between scintillators, which is the reason for the 3 15 He-detector measurements for the decay of B, and only choice of BaF2 in the spectrometer of Ref. [107] since the one measurement for 49K[99]. Several measurements of neutron-induced background is minimized. the delayed neutron emission of 15B[81, 82, 98] and 49K On the other hand, Ref. [106] proposes to turn the [86, 99] have been published using TOF arrays and show large capture cross section of NaI(Tl) into an advantage, deviations up to 10% in the reported neutron energies. allowing the measurement of neutron-emission probabili- Monte Carlo simulations are a complementary tool ties and spectra. Beta-delayed neutrons captured in the to calibrations since they offer a large flexibility and scintillator produce a bump in the energy spectrum above can reach a high level of accuracy in the simulation of ∼ 6.8 MeV, whose shape depends on the neutron-energy neutron interactions and detector properties [100, 101]. distribution and the height is proportional to the P1n They are used for extending the overall efficiency cali- value. bration of a TOF spectrometer to a broader neutron en- The difficulty of the method is that it relies on the ergy range than available with mono-energetic neutron Monte Carlo simulation of the transport and generation sources. Moreover, Monte Carlo simulations can also be of radiation by neutrons in the detector to extract the applied for determining the time-energy relation of the desired information [108]. Contrary to the simulation of neutron detectors, quantifying the cross-talk and estimat- photon transport, the accuracy of such simulations is not ing the time-correlated neutron and γ-ray background guaranteed a priori and should be investigated by com- due to elastic, inelastic and (n,γ) reactions. Many new parison between different codes and nuclear data libraries measurements have been performed in the past few years [100] and, if possible, by dedicated measurements [109]. and the results are presently being analyzed. Alternatively, a collection of neutron emitters with well-known P1n values and neutron spectra can be mea- sured with TAGS and compared with the results of the standard methods. This was done in Ref. [106] for the C. New methods (for Pn values and energy 137 spectra) case of I, and a good agreement was found with the P1n value recommended in this CRP [36] (see TableI). 1. Total absorption γ-ray spectroscopy However, in a similar work reported in Ref. [108], the 137 P1n value obtained for I was 11% smaller, while the value for 95Rb was found to be 22% larger than the eval- Total absorption γ-ray spectroscopy (TAGS) was first uation [36]. It is clear from this that further work is applied to β-decay studies at ISOLDE-CERN in the needed, both experimentally and from the point of view 1970’s [102] as a tool to measure accurately β-decay in- of simulations, to clarify the situation and establish the tensity distributions. It has been used since then to tackle systematic uncertainties of the method. This is worth the a variety of problems, including the investigation of reac- effort – given the sizable number of βn emitters that have tor decay heat evolution [103], the nuclear level densities been measured recently and will be measured in the near and γ-ray strength functions far from stability [104], and future with the TAGS technique. the spectrum of reactor anti-neutrinos [105]. Recently it was proposed to use this technique to measure P1n val- ues and delayed neutron spectra [106]. This is indeed an exciting possibility, the current status of the method is summarized here.
11 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
2. Ion-recoil methods ciple is the same as for the discovery of the two-proton radioactivity: a gas-filled TPC is used to produce ioniza- The Beta-decay Paul Trap (BPT) [110], an open- tion tracks of the incoming βn-precursor and the recoiling geometry linear radiofrequency quadrupole (RFQ) ion βn-daughter products which are then made visible via a trap designed for precision β-decay studies, has been used CCD camera. Fig.4 shows two different decay events of 8 to determine β-delayed neutron (βn) branching ratios and the βn emitter He. The picture shows only the incoming 8 7 energy spectra from the recoil energy imparted to the He nucleus and the recoiling Li, the β-delayed neutron daughter ions following β decay. This novel way to per- is not visible. In the bottom of Fig.4 the second decay 8 form βn spectroscopy circumvents the challenges associ- possibility of He into an α, a triton, and a neutron is ated with direct neutron detection. shown. The BPT collects and suspends radioactive ions in vac- In principle one can reconstruct the energy of the emit- uum at the center of a radiation-detector array consist- ted neutron from the length of the recoil tracks. This has 8 ing of ∆E-E plastic scintillator telescopes, microchan- been done for the β-delayed triton decay of He via mea- nel plate (MCP) detectors, and high-purity germanium surement of the momenta of the recoiling charged parti- (HPGe) detectors, used to detect the β particles, re- cles (α and triton) [115]. coiling daughter ions, and γ rays, respectively, emitted The main limitation comes from the generally low re- following the decay. Following β decay, the recoil ions coil energy and the finite resolution (diffusion effects) emerge from the trap volume essentially unperturbed by of the OTPC. So the technique is limited to low-mass scattering and the TOF of the recoiling ion to the MCP precursors and reasonably high neutron energies. In the detectors is determined. The βn decays can be identi- present setup with a gas mixture of 95% He and 5% N2 fied because the higher-energy recoil ions characteristic at atmospheric pressure, a lower threshold for neutrons of neutron emission (with kinetic energies up to tens of of ≈ 1.4 MeV corresponding to a minimal track length of keV) have shorter TOFs than the other recoiling daugh- 10 mm [115] is required for the recoil identification. ter ions (which have kinetic energies typically less than The OTPC method is not suitable for measuring P2n 500 eV). values (and higher) since the energy reconstruction is not The βn branching ratios can be deduced by comparing possible with two outgoing (undetectable) neutrons. the number of β-ion coincidences with TOFs character- istic of neutron emission to the total number of β decays from trapped ions obtained from the number of detected: (1) β particles emitted from the trapped precursor, (2) β-ion coincidences with longer TOFs characteristic of β- decay to bound states, and (3) β-delayed γ rays. The agreement obtained for the βn branching ratios when using these three approaches provides additional confidence in the results. This approach was first demon- strated in Ref. [64] and has subsequently been used to study mass-separated ion beams of 137,138,140I, 135,136Sb, and 144,145Cs [111, 112] delivered by the CARIBU facility at Argonne National Laboratory. Additional upgrades to the experimental setup can be implemented to increase the solid angle of the detector array and to improve the ion collection and cooling while minimizing the effect of the electric fields on the recoil ions. Future experiments will also benefit from increases in the intensity and purity of the low-energy beams de- livered by the CARIBU facility [113, 114]. FIG. 4. (Top) β-delayed neutron decay of 8He into 7Li and a The recoil-ion approach is well suited to reconstruct neutron (not visible) in the OTPC. (Bottom) β-delayed triton the kinematics of a single undetected particle through decay of 8He into an α, a triton, and a neutron (not visible) conservation of momentum and energy. The studies of β- in the OTPC. Pictures courtesy of Janas Zenon, University of delayed multi-neutron emission will need to incorporate Warzaw/ Poland. auxiliary neutron detectors.
3. Optical Time Projection Chamber 4. Measurement in a storage ring
The potential of using an optical Time Projection The Pn value could also be measured via ion counting Chamber (OTPC) to measure the neutron branching ra- in a heavy-ion storage ring. This complimentary method tio via ion counting has been shown in [65]. The prin- has been proposed for radioactive beams produced at the
12 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
Fragment Separator (FRS) at the GSI Helmholtz Center nation of masses and half-lives [67, 119]. The setup at the for Heavy Ion Research in Darmstadt/ Germany which Fragment Separator (FRS) at GSI Darmstadt/ Germany, are injected into the Experimental Storage Ring (ESR). the so-called FRS Ion Catcher [120], consists of a Cryo- Similarly, at the future Facility for Antiproton and Ion genic Stopping Cell (CSC) and a MR-TOF spectrometer. Research (FAIR) the βn-emitter will be produced at the The composition of the fast (E > 500 MeV/u) beam is Super-FRS and injected into the Collector Ring (CR). identified with the standard set of detectors at the FRS. While the βn-emitter is circulating in the ring and de- A beam containing only the isotopes of interest is stopped cays, the measurement of the neutron-branching ratio and stored in the CSC. After some decay time the ions could then be carried out via the ion counting method are then extracted from the CSC into the MR-TOF for (see Sect. IIA8). This method has been described in further identification and measurement. The relative in- Ref. [66] but no proof-of-principle measurement has been tensities of βn precursors and respective decay daughters carried out yet. will then allow to extract the neutron-branching ratios. The advantage of this method is that symbiotically the Fig.5 gives a schematic view of the setup and a sim- mass and/or half-life of the stored ion could be measured ulated MR-TOF mass spectrum for the yet unmeasured non-destructively by time-resolved Schottky mass spec- β2n decay of 137Sb. Within the next years this promising trometry [116]. However, Schottky mass spectrometry method will be tested with radioactive beams at the GSI requires “cooling” to reduce the momentum spread of Helmholtz Center in Darmstadt/ Germany. the beam. These beam cooling procedures take several seconds and thus reduce the range of accessible βn pre- cursors. Capacitive Schottky pick-up plates provide the revolu- tion frequencies of the ions in the storage ring and allow their identification. When these ions decay, the change in mass and charge state leads to a change in the respec- tive revolution frequencies. With Schottky spectrometry it is only possible to study the daughter nuclei that re- main within the acceptance of the ring. This restricts the method in the ESR to parent-daughter pairs within an A/q change of ±1.5%, but it can be complemented by particle detectors in pockets behind the dipole magnets to detect also β−/ βn-daughter nuclei. Like all described ion-counting methods this method is independent of the energy of the emitted neutron and the FIG. 5. Schematic view of the MR-TOF technique respective neutron detection efficiency. However, differ- using a Fragment Separator (FRS), a Cryogenic Stop- ent loss mechanisms during the storage in the ring need ping Cell (CSC) and a MR-TOF device to count the to be carefully addressed and investigated. ions. At the bottom a simulated MR-TOF mass spec- trum is shown. Picture reprinted with permission from [67] The decay of the ions can occur anywhere on the ring (https://creativecommons.org/licenses/by/4.0/). orbit (which is ≈ 108 m for the ESR at GSI Darmstadt) but only those decays that occur in the long straight sec- tions can be detected by a particle detector behind the respective dipole section. All other decays outside this D. Outlook for future measurements section will lead to a loss of the ions. For the ESR, a multi-purpose particle detector (CsISiPHOS) has been installed and tested [117], which allows a geometrical ef- There has been a lot of progress in measuring βn emit- ficiency of ≈ 20%. For experiments in the CR at FAIR ters in the last decade, especially with precision bet- two particle detectors will be installed and operated by ter than 5%. The earlier measurement campaigns were the ILIMA collaboration [118]. mainly driven by fission reactor studies and thus fo- cused around the two – low mass and high mass – fis- sion fragment groups. The reactor community still strives for more precise measurements of important fission frag- 5. Ion counting using a MR-TOF ments that are the main contributors to the delayed neu- tron yield νd (see priority list on p. 15 in Ref. [40]) but Multi-Reflection Time-Of-Flight (MR-TOF) spectrom- the general scientific focus has now shifted. eters are becoming more and more popular as devices The nuclear astrophysics community is interested in to measure precisely masses, determine the beam com- measurements of the most neutron-rich, heavy (Z > 26) position, or to manipulate and clean ion beams. This nuclei that can be accessed with the new generation of versatility has also led to the proposal to measure radioactive beam facilities, especially those around the neutron-branching ratios (including multi-neutron emis- N = 50, 82, and 126 shell closures which are the pro- sion branches) with a MR-TOF, parallel to the determi- genitor nuclei that are responsible for the creation of the
13 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. solar r-process abundance peaks at A ≈ 80, 130, and 195. each measurement, from all the available measurements Masses, decay half-lives, neutron-capture cross sections, published in peer-reviewed journals as well as in certain and neutron-branching ratios of these neutron-rich nu- selected secondary references. This was followed by a se- clides are important nuclear physics input parameters for lection of the most reliable set of data for a nuclide which astrophysical network calculations, and in turn also help could be averaged to obtain recommended values. to improve theoretical nuclear models that are required Details of the procedures and policies of evaluations for the prediction of even more exotic nuclei. for this CRP are described in two publications. A de- A detailed review of the r-process nucleosynthesis, the tailed paper on the compilation and evaluation of data astrophysical scenarios, the required nuclear data, how in the Z = 2 – 28 region (cut-off date December 2014) they can be measured at radioactive beam facilities and was published by M. Birch et al. [35] in 2015. The newly with what accuracy, is provided in Ref. [121]. Figs.6 and evaluated data for Z = 29 – 83 (cut-off date August 2020) 7 give an overview of the status quo of βn-emitter as of can be found in J. Liang et al. [36]. August 2020. These pictures include the vast amount of measurements that were performed during the course of the evaluation work that has been published in Refs. [35, B. New evaluated (T1/2,Pn) data 36] (see following Sect.III), as well as the majority of the (yet unpublished) data from the BRIKEN project [37, 51] In total, as of August 2020 data for nuclide identi- which will conclude its campaign in 2021. fication, half-lives and Pn values have been extracted This decade will see the transition to the new genera- and evaluated for 653 neutron-rich known or potential tion of radioactive ion beam (RIB) facilities that are still β-delayed neutron emitters from more than 750 publica- under construction. Two of the main competitors for the tions, mostly in main-stream peer-reviewed journals. study of the most neutron-rich βn-emitter will be the Fa- The evaluations in Refs. [35, 36] included only the very cility for Rare Isotope Beams (FRIB) at Michigan State first batch of data from the aforementioned BRIKEN University/ USA, which is planning to be operational by measuring campaign at RIKEN. The upcoming publica- the end of 2021, and the Facility for Antiproton and Ion tions will be included in evaluations of these mass ranges Research (FAIR) in Darmstadt/ Germany which is aim- in the near future. The cut-off date for the data used ing at producing the first beams in ≈ 2027. Until then, for the calculations of macroscopic data in this CRP pa- the “FAIR Phase-0” program will be operating at the ex- per (see Secs.VI andVII) was March 2020. The data isting GSI facilities, using many of the newly designed used in these calculations can be found in the reference detector setups and upgrades of the accelerator facilities. database [43]. These new facilities will push the limits on what can be Some data that was published in the months between measured on the neutron-rich side of the chart of nuclides March 2020 and the cut-off date of the Z >28 evaluation closer towards the driplines. Two of the key regions for was included in the published version [36]. This data was future r-process experiments at these new facilities are not used for the calculation of macroscopic parameters the so-called “Rare Earth Region” (A = 160 − 170) and in Secs.VI andVII as already mentioned, however, the the nuclei with Z < 75 along the N = 126 shell closure. differences are minute and have a negligible impact on As can be seen in Fig.7 and in the recently published the calculated parameters. evaluation [36], only a small number of β-delayed neu- The most noteworthy addition is the inclusion of data tron emitters has been measured so far for Z > 57. Ex- published by the EURICA collaboration at RIKEN [123], periments in the next decade will allow us to access these measuring half-lives of 55 neutron-rich nuclei between neutron-rich nuclei for the first time and achieve a better 134Sn and 155La. In addition, for 13 neutron-rich nu- understanding of the formation of these two abundance clei the β-decay half-life was measured for the first time: peaks in the solar r-process abundance distribution. 140−142Sb, 139−144Te, 144−146I, and 148Xe. This led, for example, to tiny changes in the recommended half-lives 144,145,146 for the P1n standards Cs: III. COMPILATION AND EVALUATION • 144Cs: 581(6) ms −→ 582(6) ms
A. Compilation and evaluation methodology • 145Cs: 321.7(11) ms −→ 321.7(10) ms
Standard procedures were followed for the evaluation • 146Cs: 229.7(10) ms −→ 229.5(10) ms of nuclear spectroscopic data as defined for the Evalu- ated Nuclear Structure Data File (ENSDF) and discussed at three Research Coordination meetings of the CRP 1. Z = 2 – 28 region [39, 40, 42], with the addition of guidelines for half-life evaluations as formulated by A.L. Nichols and B. Singh Based on Q values deduced from the Atomic Mass Eval- (see p. 31–40 in Ref. [122]). Emphasis has been placed on uation (AME) 2016 for neutron-rich nuclides [33, 124], a making sure that the Pn value and half-life data have been total of 221 experimentally observed nuclides were iden- compiled thoroughly, together with the method used for tified as known or potential β-delayed neutron emitters
14 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 6. Chart of nuclides for Z = 0 – 14 (top), Z = 15 – 28 (middle) and Z = 29 – 43 (bottom) with status for half-life and neutron-branching ratio measurements. Boxes show previously measured neutron-branching ratios, whereas circles indicate isotopes that have recently measured but have not yet been published or included in the latest evaluations [35, 36]. Crosses indicate isotopes for which half-lives have been measured for the first time recently.
15 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 7. Chart of nuclides for Z = 44 – 57 (top) and Z = 58 – 86 (bottom) with status for half-life and neutron-branching ratio measurements. Boxes show previously measured neutron-branching ratios, whereas circles indicate isotopes that have recently measured but have not yet been published or included in the latest evaluations [35, 36]. Crosses indicate isotopes for which half-lives have been measured for the first time recently.
16 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. in Ref. [35]. Of these nuclei, half-lives had been mea- • Measured using reliable methods outlined in sured for 180 nuclides, P1n values for 114, P2n values for Sect.IIA in four or more independent experiments. 21, and P3n values for only three nuclides. There is only 17 Consistent results with overall uncertainty of less one nucleus, B, where a tentative P4n value has been • reported [7]. In addition there were 9 isomeric values for than 5%. P values reported. For about 50% of the cases, P 1n xn Exceptions are 49K, 82Ga, and 146,147Cs, where the over- values are available from only one measurement. all uncertainty is 5–10%. 147Cs has only three indepen- Since the cut-off date for this part of the evaluation was dent measurements, but these measurements are very ro- in 2015, work on an update to include recent publications bust and consistent. has started by the Canadian community. So far, four In summary, the following nuclides are recommended new nuclides have been identified as βn-emitter , and two as P standards, covering different Z regions up to decay half-lives and three new P values were measured. 1n 2n Z=55: 9Li, 16C, 17N, 49K, 82Ga, 87Br, 88Br, 94Rb, 95Rb, 137I, 145Cs, 146Cs, and 147Cs. The recommended half- lives and P1n values for these standards can be found in 2. Z = 29 – 57 region TableI.
The latest evaluation comprised the heavier mass re- TABLE I. Table with P1n standards and their corresponding gion and was recently published [36] with a cut-off date half-lives as recommended from this CRP [35, 36]. of August 2020. A total of 318 nuclei were identified as Nuclide Half-life (s) P1n (%) Ref. β-delayed neutron emitters, based on Q values deduced 9 from the AME2016 [33, 124]. This region is of utmost Li 0.1782(4) 50.5(10) [35] 16C 0.7546(80) 99.28(12) [35] importance due to the production of most of these nu- 17N 4.171(4) 95.1(7) [35] clides in the fission of uranium and other actinides, and 49K 1.263(50) 86(9) [35] therefore of relevance in nuclear power reactors. Out of 82Ga 0.601(2) 22.7(20) [36] these 318 nuclei, half-lives have been measured for 284 87Br 55.64(15) 2.53(10) [36] nuclides, P1n values for 183 nuclides, and P2n values for 88Br 16.29(8) 6.72(27) [36] only eight nuclides. So far no P3n or P4n values have been 94Rb 2.704(15) 10.39(22) [36] determined in this mass region. 95Rb 0.378(2) 8.8(4) [36] The data for isomers are generally poorly determined 137I 24.59(10) 7.63(14) [36] 145 and often the half-lives and Pn values are mixed with Cs 0.582(6) 13.5(6) [36] the data for ground states. For a large number of cases 146Cs 0.3217(10) 14.3(8) [36] 147 P1n values are available from only one measurement, and Cs 0.2295(10) 28.5(20) [36] for several others only upper limits are reported. How- ever the data for the main isotopes contributing to fission yields seem to be in good order.
D. Systematics for P1n values
3. Z = 58 – 87 region In several studies [127–130] correlations between known values of P1n and other β decay gross properties (T1/2, In this mass region 110 experimentally known nuclides Qβ, Sn) have been described to estimate unknown P1n were identified as β-delayed neutron emitters based on values. The basic correlation between P1n and β-decay Q values from the AME2016 [33, 124]. Half-lives are properties is: known for only 54 of these, and P1n values for only 9 nuclides. Measurements for six Tl isotopes (211−216Tl) 210,211 R Qβ and two Hg isotopes ( Hg) were reported only re- S Sβ(E)f(Z,Qβ − E) dE 210 n cently [125, 126], while a tentative value of P for Tl Pn = Q , (9) 1n R β S (E)f(Z,Q − E) dE was reported in a secondary publication in 1961 [34]. 0 β β
where Sβ(E) is the β-decay strength function and f(Z,Qβ − E) is the Fermi integral. The integrations is C. P1n Standards performed over the available β energy window, Qβ − E. One of the earliest semi-empirical expressions origi- The selection of certain radionuclides as standards for nates from Amiel et al. in 1970 [127]. Already then it was P1n follows criteria which are given in Refs. [35, 36]. The known that the neutron emission probability depends on main criteria for “good” standards are: the available Qβn energy window, the level density, and the competition of neutron and γ emission. A plot of the • Easy production at various facilities (large quanti- log P1n vs. log Qβn resulted in a linear dependence with ties, clean beams, no isomers). slope m. This led to the conclusion that the dependency
17 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
of the P1n value on the Qβn window is dominant over any improved systematics was developed by McCutchan et. influence from the level density. al. [129] which considers both, the β-delayed neutron The evaluated half-lives and P1n values in this CRP emission probability as well as the half-life. As the half- were analyzed by the three empirical approaches pro- life is inversely proportional to the decay Q value, the posed by K.L. Kratz and G. Herrmann [128] (Kratz- ratio P1n/T1/2 is considered as a function of the Qβn Herrmann formula, KHF), an improved KHF system- value: atics developed by E.A. McCutchan et al. [129], and a level-density dependent systematics by K. Miernik [130]. These estimates can be used as guidance for future experi- P1n d ∼ c Qβn, (11) ments close to the known region, as well as to cross-check T1/2 if outliers in existing data appear. They are however of limited use for r-process abundance calculations since where c and d are fitted parameters. they do not reach out far enough in A. For Z = 2 – 28 region, global fits with the Kratz- Herrmann and McCutchan et al. systematics were pre- 3. Improved systematics along Z chains sented in Ref. [35]. More detailed systematic plots using all the three approaches are presented for the Z = 28 – While the P1n/T1/2 systematics from McCutchan show 57 region in Ref. [36], where a new approach of fitting the an improvement over the KHF, there is still considerable McCutchan et al. systematics by individual Z has also scatter. ”Isotopic” systematics along a given Z can pro- been presented for most of the elements in this region. It vide a much better picture for unknown P1n values. In appears that this approach has better potential for pre- Ref. [36] this method was used for nuclides in the region dicting unknown P1n values. Unfortunately, there is not Z = 29 − 57 which had three or more isotopes with a enough data in the Z = 58 – 87 region to attempt fits measured value for their P1n and T1/2, respectively. with the systematics. The improved systematics along a given Z shows linear behavior when plotted in a log-log scale. Despite some scatter in this presentation, extreme outliers can indicate 1. Kratz-Herrmann formula isotopes that require further investigation, e.g. the iso- topic chains of As, Pd, and Ba [36]. In the Kratz-Herrmann formula (KHF) [128] the β- strength function is assumed to be constant above a cer- tain cut-off energy C, and to be zero below that energy. 4. Systematics based on the Effective Density Model It also assumes that the Fermi integral might evolve as a function of the Qβn value. This leads to the following The phenomenological effective level density function description: model by Miernik [130] assumes that the β = strength function is proportional to the level density ρ(E) fed by β-decay: b Qβn P1n ∼ a , (10) Qβ − C √ exp(a E) S (E) ∝ ρ(E) = d (12) where a and b are free parameters (fitted to the re- β E2/3 spective data range), and C is the cut-off parameter. The values for the cut-off parameter depend on the pro- The parameter ad is fitted to existing data, so that the ton/neutron number: derived formula can be used across an entire mass surface. In order to calculate Pn using this systematics, one • C = 0 for even-even nuclei. needs to first find the ad for the nuclides of interest. Then √ the β-decay energy Qβ of the parent and the neutron C = 13/ A for even-odd/odd-even nuclei. • separation energy Sn can be used within the integration √ of Eq.9. The Fermi integral can be calculated precisely C = 26/ A for odd-odd nuclei. • with available tables, and the integration of Eq.9 can be For the latest fit parameters and a log-log plot of Pn as carried out numerically. b An attempt to extend this model to multi-neutron a function of Qβn , the reader is referred to Ref. [36]. Qβ −C emission (up to P3n) can be found in [131]. The assump- tion is that the neutrons are emitted sequentially, and that the many-body channels are negligible. The limi- 2. McCutchan systematics tation of this description is also that only Gamow-Teller transitions are considered and that the neutron-γ com- The considerable scatter and large reduced χ2 in the petition is approximated with a ”threshold method”. plots of the KHF systematics questions the reliability of Numerical values for the fit parameters and more de- predictions made employing this systematics. Thus, an tails are given in Refs. [36, 130].
18 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
E. Delayed neutron spectra for individual or inter-comparison purposes. This issue is further dis- precursors cussed below.
Delayed neutron spectra have been measured for less than 50 βn emitters. These include 34 fission products for which good quality data exists which was obtained 1. Digitization of delayed neutron (DN) energy spectra using 3He and gaseous proton recoil spectrometers and measured by several groups at different installations. As Digitization of several β-delayed neutron spectra has it was pointed out in the late 1970’s [19] discrepancies been done using the GSYS2.4.3 code and its later ver- were observed between the different sets of data in par- sion 2.4.7 [133]. This code is routinely and extensively ticular at low neutron energies. This lead to an attempt used for the digitization of cross section plots for the EX- to benchmark the different methods at different labora- FOR database [134]. The quality of the digitization was tories [70] and triggered new measurements at other fa- checked by plotting the extracted data and comparing cilities [17, 73]. the plot with the corresponding figure in the paper. An evaluation of the available neutron spectra for the The extracted datasets of all the figures were sent to 34 individual fission products was carried out during the the IAEA Nuclear Data Section for further checking by thorough and comprehensive work presented in the PhD an expert staff scientist. All the data files were finally Thesis of M. C. Brady in 1989 [31]. The list of isotopes prepared for inclusion into the Reference Database for included in the evaluation is given in Appendix D of delayed neutrons [43]. that work. Most of the experimental data was coming a. Z = 2 – 28 region: Beta-delayed neutron spectra from 3He spectrometer measurements carried out by the from 8He [135], 9Li [136, 137], 11Li [5, 83, 138, 139], 14Be Studsvik group (G. Rudstam et al. ), the Mainz group [138], 15B[81, 98], 16C[140], 17B[141], 17N[81, 98, 140], (K.-L. Kratz et al. ) and the Pacific Northwest Labora- 18N[142], 21N[143, 144], and 27,28Na [145] were digitized. tory (PNL) group (P. Reeder et al. ). Proton recoil data Note that some of these are time-of-flight spectra, e.g. for from the Idaho National Engineering Laboratory (INEL) 9Li from Ref. [136]. group (R. C. Greenwood et al. ) was also included. Addi- b. Z = 29 – 57 region: Spectra are available from tional data was published by this group in 1997 [18] and different groups: thus was not included in the Brady evaluation. The data from the Studsvik group was provided with R.C. Greenwood and A.J. Caffrey [17]: Spectra of uncertainties. The uncertainty for the other data sets • delayed neutrons for the isotope-separated, fission was inferred from generic assessments by the respective product precursors 93−97Rb and 143−145Cs were authors. Brady [31] stated that the Mainz and Studsvik digitized independently by two groups over an en- data are in ”fair agreement” above few hundred keV, and ergy region of ≈ 10 – 1300 keV. The data was trans- are considered to be of better quality than PNL data. formed to an equidistant scale with an energy bin Probably because of this she gave preference to a single of 1 keV for use in the summation method. data set instead of averaging different data sets, when available. In general, due to the smaller uncertainties and broader • P.L. Reeder et al. [73]: Measured energy spectra of energy range, preference was given to 3He data from the delayed neutrons from 93−95Rb and 143Cs were dig- Mainz group. If available, the proton recoil data from itized. INEL replaced the 3He data below 200 keV. The selected data were extended by theoretical model calculations to cover the full Qβn window. This procedure lead to a • M.C. Brady [31]: Several measured delayed neu- renormalization of the experimental data. The recom- tron spectra obtained from different experimental mended data is presented in graphical form (figures 12 to groups (mainly Mainz, Studsvik, and groups at 45 in Ref. [31]). the TRISTAN online separator) were analyzed in The originally evaluated data from Brady was incorpo- this work as raw and adjusted spectra. This work 79−81 85 87−98 129,130 134 rated into the ENDF/B-VI decay data sub-library [132]. includes Ga, As, Rb, In, Sn, 135 136 137−141 141−147 Uncertainties are not included explicitly in the recom- Sb, Te, I, and Cs. Both the raw mended data of [31] neither can they be retrieved from as well as the adjusted spectra from this work were ENDF/B. In addition an evaluation of possible system- digitized. atic errors from the comparison of different data sets was not carried out in Brady’s work. The lack of realistic un- • R.C. Greenwood and K.D. Watts [18]: Delayed neu- certainties for individual precursor spectra diminishes its tron spectra of 87,88,90Br, 137−139I, and 136Te were helpfulness. For example, it does not allow for an eval- presented in this paper. All were digitized, and were uation of uncertainties on aggregate spectra calculated subsequently transformed to equidistant scale with by the summation method. In addition it is difficult to an energy bin of 1 keV for use in the summation define reference β-delayed neutron spectra for calibration method.
19 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
2. Assessment of spectra as a reference for calibration 0.025 purposes 0.02 [a.u.]
n 94 Defining β-delayed neutron reference spectra would be I Rb an important asset in practical work. These could serve 0.015 to benchmark new setups and methods for the measure- 0.01 ment of β-delayed neutron spectra. It will also serve to study the neutron energy dependency of the efficiency in 0.005 devices aiming to determine neutron emission probabili- ties. The lack of an evaluation of uncertainties, in partic- 0 ular systematic uncertainties, introduces ambiguities in 0 0.5 1 1.5 2 2.5 E [MeV] the results. n An attempt was made within the CRP to asses the 0.025 size of systematic uncertainties by comparison of spectral data obtained by different groups and different methods. 0.02 [a.u.]
n 95 However, it was realized that with the available informa- I Rb tion it is not possible to make a quantitative assessment. 0.015 To illustrate the difficulty of the task, the cases of 94Rb 0.01 and 95Rb are shown in Fig.8, which have reasonably large P1n values and are easily produced in fission reactions. 0.005 In the evaluation work of Brady [31], 3He spectrometer data was available for both nuclei, provided directly by 0 the Mainz group and the Studsvik group. Data was avail- 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 E [MeV] able also from the PNL group [73], but was disregarded n because it had lower statistics and a smaller energy range. Preference was given to Mainz data over Studsvik data FIG. 8. Comparison of measured β-delayed neutron spectra for both isotopes, but no information is given about the for 94Rb (top panel) and 95Rb (bottom panel). Red line: Data differences between both data sets. Below 200 keV the retrieved from ENDF/B-VII.1 database. Blue line: Data from Mainz data was replaced by proton-recoil data from INEL Greenwood et al. (INEL) [17]. Green line: Data from Reeder [17]. et al. (PNL) [73] Figure8 displays the resulting evaluated spectra as re- trieved from the ENDF/B-VII.1 library [24] (same spec- tra are available in ENDF/B-VIII.0 [25]). It is compared For practical applications the interested user should be with the full INEL spectrum obtained by digitization of warned about possible large uncertainties in the shape of the graphical representation in [17]. The data from PNL β-delayed neutron spectra retrieved from the ENDF/B- is also included in the figure, and was obtained by dig- VIII.0 [25] and ENDF/B-VII.1 [24] libraries. Whenever itization of figures in [73]. The three data sets are his- possible the effect of spectra obtained from different mea- togramed with a 10 keV bin width and are normalized surements should be studied. It should be noted also that to the counts in the region from 0.2 to 1.2 MeV. As can the experimental spectra in the ENDF/B data base are be seen there are notable differences. When comparing augmented with model calculations to cover the full Qβn values integrated over a range of 100 keV, to minimize energy window. These calculations use FRDM+QRPA the effect of different energy resolution and statistics, the β-strength distributions [147] and Hauser-Feshbach neu- discrepancies reach up to 30–50%. The drop observed in tron branching ratios [32] to obtain the neutron spectra the PNL data at around 0.5 MeV is likely due to a sud- as described in the following sections. The model spectra den rise of their efficiency that the authors quote as of additions are easily identified due to the smooth charac- uncertain origin but was also observed by other groups. ter of their distribution. Later this rise of efficiency has been identified as a spuri- ous effect of neutron resonances excited in iron which is generally present in the experimental setup [146]. IV. THEORY OF MICROSCOPIC DATA It is not easy to arrive to a conclusion regarding the size of systematic errors from this information alone. In In this section we compare the new evaluated CRP particular the loss of the primary information handled (T1/2,Pn) data described in Sect.III with “state-of-the- by Brady is unfortunate. In addition, it seems that the art” microscopic models. The comparisons are limited systematic errors of the different techniques are not fully to global self-consistent β-decay models that take into ac- understood. A proper evaluation work will require new count both the allowed Gamow-Teller (GT) and first for- well-controlled measurements that should profit from the bidden (FF) decays on equal footing and can be used to advances in radioisotope production methods and instru- calculate both half-lives and β-delayed neutron emission mentation made in recent years. probabilities. Self-consistent models have the advantage
20 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. of being founded on first principles. Their parameters are are estimated by considering the competition be- determined from a limited set of sample nuclei and are tween neutron and γ emission in the daughter nu- then kept unchanged for the whole nuclear chart. Thus, cleus using the Hauser-Feshbach statistical model they are considered more reliable for extrapolations to (HF). This approach ‘FRDM12+(Q)RPA+HF’ is extreme N/Z ratios. However, before they are used for different from the sharp cut-off model used in the predicting nuclear properties in unexplored mass regions, previous calculations [147] whereby multi-β-delayed their accuracy has to be tested against existing exper- neutron emission emission was assumed to occur imental data. The aim of this section is to ascertain at energies above the corresponding neutron sepa- whether these global models are robust in the descrip- ration. tion of the evaluated CRP data which include new re- In addition to the above models, which are used in iso- sults from ongoing large-scale experiments at RIB facili- topic chains and global comparisons of both T and P ties worldwide. This will then mean they can be used in 1/2 xn extrapolations to mass regions where experimental data values, we include the results from the finite amplitude are either unavailable or inaccessible. method (FAM) [156] in selected comparisons of half-lives. Beta-decay properties are influenced by various differ- The FAM is a mean-field model used up to now solely for ent factors. First of all, the competition between the calculations of half-lives, however, it provides a break- GT and high-energy FF transitions accelerates the β- through approach for treating deformed nuclei [157–159] decay, while in the region of the ground-state spin inver- with the Skyrme functional SkO’ [160], therefore is worth sion such concurrence may produce isotopic irregularities comparing with the CRP data for deformed nuclei. A few attempts have been made to include complex in the half-lives and Pn values. Second, the effects of many particle-many hole (np-nh) configurations substan- np-nh configurations within the self-consistent approach tially enriches and softens the β-decay strength distribu- in the study of the fine structure of the β-decay strength tion. Third, a local fragmentation of the β-strength due functions and integral β-decay properties. The preferred to deformation may change nuclear β-decay characteris- models are the latest versions of the interacting shell- tics. In transitional nuclei a co-existence of spherical and model (SM) [161, 162] which include the allowed GT deformed shapes may also be of importance. and FF decays. Also “beyond the QRPA approxima- In the comparisons, we use two self-consistent global tion” models based on the EDF approach have been de- models which treat GT and FF transitions consis- veloped and applied to selected nuclei near the closed tently, but do not take into account complex np- shells. Practical schemes of this kind of models are the nh configurations or deformation effects. For com- phonon-phonon coupling model (PPC) within the FRSA pleteness we also include the results of the recently approximation [163–165] and particle-vibration coupling updated microscopic-macroscopic model of Ref. [148] models (PVC) [166–168] of the β-decay strength function. which are taken from the publicly available supplemen- However, these models are still limited to the allowed GT approximation. Among the models with complex con- tary tables (https://t2.lanl.gov/nis/molleretal/ figurations, β-delayed neutron emission branching ra- publications/ADNDT-BETA-2018.html). tios have been tackled so far only in the SM frameworks The three models are listed below: [161, 162] and FRSA [164, 165]. We have included pub- • Fayans Energy Density Functional (EDF) [149] lished results from the SM [161, 162], FRSA [163–165], combined with continuum quasiparticle random- and PVC models [166–168] in the comparison of half- 78 phase approximation (CQRPA). A common abbre- lives of Ni in this report. viation DF+CQRPA is used below for different ver- Finally, we mention the self-consistent models for sions of the density functional, namely DF3 from spherical and deformed nuclei in the GT approximation [150–152] as well as the extended DF3a which was [169–174] which have been used in selected isotopic re- developed within the CRP to compare with the new gions. Ref. [170] in particular is dedicated to the next evaluated data. generation of JENDL evaluated decay data, including β-decay half-lives, β-delayed neutron branching ratios, • Relativistic Hartree–Bogoliubov+QRPA [153] as well as neutron energy spectra. A detailed review based on the D3C∗ functional of the self-consistent models of β decay can be found in (RHB+RQRPA) [154]. This model was also Refs. [175–177]. extended within the CRP to perform global Sect.IVA discusses the extensions to the two self- calculations. consistent models that were performed within the CRP. Comparisons of the models with the evaluated CRP • Microscopic-macroscopic model of Ref. [148] based (T1/2,Pn) data are shown in Sect.IVB for: (i) isotopic on the latest version of the finite-range droplet chains of nuclides in the fission mass region, (ii) isotopic model FRDM12 published in Ref. [155] combined chains of nuclides relevant to nuclear astrophysics and with the BCS+RPA ‘(Q)RPA’ description for the (iii) global comparisons with respect to the (T1/2,Pn) GT decays and a phenomenological approach based data and mass number A. Finally, theoretical DN spec- on the ‘Gross theory’ for first forbidden (FF) de- tra are compared with experimental spectra from Brady’s cays. The β-delayed neutron emission properties thesis [31] in Sect.IVE.
21 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
A. Self-consistent QRPA-based models the same energy-independent quenching factor Q = 2 (gA/GA) . The one-pion component of the residual inter- The DF+CQRPA and RHB+QRPA models have been action is quenched by the same factor Q. The axial-vector extended to perform large-scale calculations of β-decay coupling constant |gA/GA| = 1.0 is used both for GT and half-lives and delayed multi-neutron emission probabil- FF decays, and the amplification of the γ5 component of ities within the CRP. These extensions are briefly de- the spin-dipole operator is accounted for as in Ref. [181]. scribed in the following. The energy-dependent effect of np-nh configurations can be included via the spreading width Γ↓ [182]. In the case of odd-A and odd-odd nuclei (hereafter 1. DF+CQRPA called “odd” nuclei), their degenerate ground states are treated through averaging the corresponding multiplet components via equal filling approximation. The interac- The β-decay strength function is calculated in the ex- tion of the odd nucleon with a “core” and rearrangement tended finite Fermi systems theory (FFST) [178] based on effects are properly taken into account [150]. the energy density functional (EDF) theory which allows In 2016, the DF+CQRPA framework was extended by for fully self-consistent description of the ground state fixing the odd nucleon in a given state. In this way, the properties of nuclei with pairing. Treating the effects of DF3a+CQRPA model can account for possible ground- momentum and energy dependence of the functional and state spin inversion effects which have been shown to effective forces on an equal footing leads one to the Fayans impact the β-decay observables [177, 183]. Though the EDF [149] with the fractional-linear dependence on nor- present version of the model is basically limited by the mal density. The surface and pairing components of the spherical approximation, it should be noted that the de- EDF contain the gradient of the density. Such an ansatz formation caused by the unpaired nucleon is properly is more general than the standard density dependence taken into account in the FFST equations for the odd- of the Skyrme functional that is a result of effective 3N odd and odd-A nuclei. correlations. Within the present CRP the ground-state and β-decay characteristics for more than 200 (quasi-) spherical nuclei were re-calculated with the DF3a func- tional [179] which differs from the previous version DF3 2. Relativistic Density Functional+QRPA by stronger effective tensor-like components. It gives the same quality of description of the ground-state properties Large-scale calculations of β-decay half-lives and β- as the DF3 functional and a better description of β-decay delayed neutron-emission probabilities were performed in of heavy nuclei [180]. Ref. [153] with the main goal of applying the results Based on such an approach a framework that allows for in heavy-element nucleosynthesis and comparing them large-scale continuum QRPA calculations of the GT and with the new CRP (T1/2,Pn) database. The model was FF β-decay properties (DF3a+CQRPA) was developed based on the relativistic nuclear energy density func- in Refs. [151, 152]. For large-scale β-decay calculations tional (RNEDF), which constructs the nuclear ground the effective approximation is used for the spin-isospin ef- state from the nucleon mean-field and a minimal set of fective NN interaction. In the one-particle-one-hole (ph) exchange mesons together with the electromagnetic field. channel it consists of the Landau- Migdal interaction and The ground-state properties were obtained from the rel- one-π and ρ-meson exchange terms (modified for nuclear ativistic Hartree-Bogoliubov model which correctly de- medium effects). For nuclei with pairing correlations, scribes the pairing in open-shell nuclei. In the calculation, the density-dependent, A-dependent, zero-range T = 1 the D3C∗ density functional, which includes momentum- ground-state pairing is used. For the spin-isospin effec- dependent terms [154], was used together with the pairing tive NN interaction in the particle-particle (pp) channel part of the Gogny D1S interaction for the description of (isoscalar T = 0 effective interaction or dynamic pairing) the pairing correlations. For the isoscalar proton-neutron a mass-independent zero-range interaction with constant (pairing-like) NN-interaction, a two-Gaussian form was strength is assumed. used following its introduction in Ref. [184]. Unlike the Allowed and first-forbidden transitions are treated on DF+CQRPA approach, which is free from fitting the pa- equal footing in terms of the full set of multipole oper- rameters to the ‘output’ half-lives, the strength of the ators that depend on the space and spin variables. The interaction was adjusted to reproduce the experimental relativistic vector operator α and axial charge operators half-lives available at the time. γ5 are reduced to their non-relativistic space-dependent The FF transitions were treated at the same level as counter parts via the relations for conservation of vector the allowed GT transitions. The same value of the axial- current (CVC) and partial conservation of axial current vector coupling constant |gA/GA| = 1.0 was used for all (PCAC). An advantage of such a scheme is that it is con- transitions, both allowed and first-forbidden, while the venient to use in the full ph-basis continuum pnQRPA amplification of the γ5 operator was not included. For framework [151, 152]. odd-A and odd-odd nuclei, the ground state energy was The correlations beyond the QRPA are included by computed by constraining the average number of particles re-scaling the spin-dependent multipole responses by to be odd which resulted in an even RHB state. Its energy
22 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
differs from the true RHB state by the energy of the odd increase of the P1n value at N = 54 as seen in the quasi-particle (or both in the case of odd-odd nuclei). FRDM12+(Q)RPA+HF and RHB+RQRPA models. This approach was tested on Qβ-values for even-Z and b. Bromine isotopes The CRP half-lives and P1n odd-Z isotopic chains and was found to be in satisfactory values for Br isotopes (Z = 35, Fig.9) are com- agreement with the data. pared with results obtained from the DF3a+CQRPA and The limitations of the model are: (i) the model con- RHB+RQRPA models. Both frameworks show a be- siders nuclei to be spherical, disregarding the possible havior that is typical of the dominance of the GT de- effects of deformation on the transition spectra, and (ii) cays. The deformation in the Br isotopes is estimated the model determines the ground state and excitations at to be higher than in As isotopes (|β2| ≈ 0.15 – 0.2) zero temperature which has implications for applications [185]. However it is observed that the T1/2 given by the to the r-process nucleosynthesis since the latter takes spherical DF3a+CQRPA and RHB+RQRPA are closer place at finite temperatures. Both approximations are to the experimental data than those from the deformed necessary to reduce the total computational cost of the FRDM12+(Q)RPA model. The latter underestimates large-scale calculations. the T1/2 by up to a factor of 25 at A < 93 and shows a sudden increase of T1/2 at A = 94. At the same time, the N dependence of P1n values in FRDM12+(Q)RPA+HF
B. Comparison with new evaluated (T1/2,Pn) data is in better agreement with the evaluated data, while the for important β-delayed neutron emitters DF3a+CQRPA calculation describes the general trend of experimental half-lives and P1n values. In TableII, the evaluated half-lives and P values for The evaluated CRP (T1/2,Pn) data described in 1n Sect.III are compared with the previously listed self- Br isotopes are compared with those obtained from the consistent global models [151, 153] and microscopic- DF3a+CQRPA (spherical option) and RHB+RQRPA macroscopic results [148] for light and heavy fission prod- models, as well as with the results from the FRDM+RPA ucts. The CRP data can be downloaded from the IAEA model of Ref. [22] that uses the spherical GT approxima- online database [43]. The selected isotopic chains con- tion (“QRPA-2”). The latter model (based on FRDM of tain the most important βn-emitter associated with the Ref. [186]) was used to estimate the total (integral) DN highest fission yields, as recommended by the evalu- yields in the evaluation of Wilson and England [187]. Fi- ated libraries ENDF/B-VIII.0 [25], JEFF-3.1.1 [28] and nally, the deformed FRDM12+(Q)RPA calculations that JENDL/DDF-2015 [27]. In order to assess the quality of include the FF decays within the Gross theory [148] are the description of the β-decay characteristics, we analyze also displayed. As can be seen in the table, DF3a+CQRPA describes both the T1/2 and Pxn values. the general trend of the isotopic distribution of the half- lives and P1n values for N < 58. It can be used for es- 1. Light fission products timating the total (integral) DN yield, as the cumulative fission yields of the DN emitters amplify the contribution of this region (see Sect. VI B 1). a. Arsenic isotopes The β-decay characteristics of c. Molybdenum isotopes The influence of deforma- the As isotopic chain (Z = 33) are compared with the tion on the β-decay properties is illustrated for the Mo calculations from the new DF3a functional within the isotopes (Z = 42) in Fig.9. Here, β is ≈ −0.3 for DF3a+CQRPA model, as well as with the RHB+QRPA 2 110 < A < 118 while for A ≈ 126 it is β2 ≈ 0, as es- [153] and FRDM12+(Q)RPA+HF [148] models in Fig.9. timated from the quadrupole moments [185]. As can be seen in the figure, the DF3a+CQRPA calcu- As can be seen from Fig.9, the T1/2 in Mo iso- lation in the GT approximation overestimates the eval- topes obtained from the spherical DF3a+CQRPA and uated half-lives by up to a factor of 10. Inclusion of RHB+RQRPA model are in agreement with the CRP the FF decays and the experimental ground-state spin π − data. They are also close to the results from the deformed (J = 5/2 ) leads to a slight underestimation (by up to FRDM12+(Q)RPA and FAM models in the N = 72 − 82 19%) of the evaluated half-lives, while the RHB+QRPA region. However, the FRDM12+(Q)RPA shows some and FRDM12+(Q)RPA result in an underestimation by irregularities in the left wing of the isotopic chain cor- up to a factor of 8. The deformation of As isotopes is responding to stronger deformation (β2 ≈ −0.35) for rather moderate and the quadrupole moments are esti- N < 72, as well as in the (quasi-)spherical region N > 82. mated to be β2 ≈ 0.1 − 0.15 [185]. The difference in At the same time, the A-dependence of the half-lives in the results of spherical and deformed calculations of the the spherical models is quite smooth at N > 82. It should half-lives is not significant. be noted that a detailed study in Ref. [188] has shown A sudden decrease of the P1n value when crossing N = some sensitivity of the Mo half-lives due to the oblate, 50 is correlated with the onset of the FF transitions. At prolate or spherical shape. N = 52, 53 a qualitative agreement is observed between d. Rubidium isotopes For the strongly deformed Rb DF3a+CQRPA and FRDM12+(Q)RPA+HF, though the isotopes (Z = 37) with β ≈ 0.35 [185], the spherical experimental P values are up to a factor of 3 higher. 2 1n RHB+RQRPA describes the experimental half-lives ex- However, the experimental data do not support a sharp
23 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
92−94 cept for Rb but does not explain the P1n values curacy comparable to the highly parametrized deformed for the important β-delayed neutron emitters 95,96Rb FRDM12+(Q)RPA framework. (Fig.9). The deformed FRDM12+(Q)RPA model tends In TableII, the DF3a+CQRPA half-lives and P1n val- to underestimate the half-lives in the region of strong de- ues for iodine isotopes are compared with the spherical formation for 92−98Rb, nevertheless it provides a good option of the FRDM+RPA model (“QRPA-2”) [22] and description of the P1n data. the recent deformed option of the FRDM12+(Q)RPA [148]. As can be seen, the FRDM+RPA, which was used in the evaluation of Wilson and England [187] for estimat- 2. Heavy fission products ing the total (integral) DN yields, often overestimates the experimental P1n values. As it does not reproduce the ex- perimental half-lives, the corresponding isotopes are not a. Iodine isotopes In Fig.9, the half-lives and P 1n correctly assigned to the I-VI group of the 6-group model values for the iodine chain (Z = 53) are shown, including that is used to describe the total (integral) DN activity the important β-delayed neutron emitters 137−139I which (see Sect.IX for more on group constants). are the main contributors to the total (integral) DN yields It can be concluded that for fission products a complex for thermal neutron-induced fission of 235U. The spheri- interplay of many factors influences the β-decay prop- cal DF3a+CQRPA and the deformed FRDM12+(Q)RPA erties. The local fragmentation of the β strengths due describe the half-lives reasonably for these moderately de- to deformation is not necessarily the main driving force. formed nuclei with β2 ≈ 0.10 − 0.15 [185]. What is important is an appropriate description of the As the FRDM+BCS+RPA framework adopted in ground-state properties: total energy release and single- Ref. [148] does not retain the SO(8)-symmetry of the particle energies. The competition of the GT and high- QRPA, it infers a notable spurious odd-even staggering energy FF transitions, effects of complex np-nh configu- in half-lives and Pn values. The RHB+RQRPA produces rations and possible ground-state spin inversion also have a very smooth N dependence and systematically under- a strong impact on the results. In transitional nuclei, a estimates the half-lives but interestingly agrees best with co-existence of spherical and deformed shapes may also the experimental P1n values. play an important role. In the DF3a+CQRPA calculation two effects are im- Summarizing, the self-consistent models give a ro- portant: a ground-state spin inversion and contribution bust description of the β-decay properties of fission of the FF decays. The first one leads to a stabilization of 141−143 products. The quality of description is comparable to the half-lives at N = 88 – 90 ( I) and the second that of the highly-parametrized microscopic-macroscopic induces a reduction in the P1n values. Notice that a sta- 144−146 FRDM12+(Q)RPA+HF model. As it will be shown in bilization of the half-lives of I was observed in the the next section, the self-consistent models are preferable recently published RIKEN data [189] (shown in Fig.9 due to their predictive power. and in TableII). b. Cesium isotopes For the isotopic chain of the Cs β-delayed neutron emitters (Z = 55, Fig.9), the de- C. Comparison for isotopes relevant to nuclear formation is β2 ≈ 0.15 − 0.2 [185]. Both, deformed astrophysics FRDM12+(Q)RPA and spherical self-consistent models, underestimate the T1/2 at the beginning of the chain (A In this section, the global models are compared with = 141) where nuclei have a lower deformation. the CRP (T1/2,Pn) data for very neutron-rich isotopes The evaluated half-lives at A ≈ 150 indicate a kind of which are at the edge of the current possibilities of on- plateau like the one observed in the iodine isotopes and going RIB experiments and yet are indispensable for the in self-consistent calculations. The evaluated P1n values modeling of the r-process nucleosynthesis in nuclear as- at A > 146 also show a stabilization effect. This could trophysics. be related to the ground-state spin inversion at N > 90 We focus on two reference semi-magic isotopic chains: [190]. Ni (Z = 28, π1f orbital is filled) near the crossing of As in the iodine case discussed before, the 7/2 the neutron major shell N = 50 and Sn (Z = 50, π1g9/2 RHB+RQRPA seems to reproduce the experimental P1n orbital is filled) around N = 82. The reason for choosing values the best yet underestimates the half-lives. On the these isotopic chains is that β-decay properties of isotopes other hand, the DF3a+CQRPA and FRDM12+(Q)RPA before and after the major neutron shell crossings offer a show a behavior that reflects the competition of GT and sensitive test for nuclear structure models. The theoreti- FF decays in the region up to N = 92. For N > 92 cal results are compared with the available CRP data and (A > 147), the DF3a+CQRPA and RHB+RQRPA de- the half-lives from the latest RIKEN experiment [189]. scribe the experimental P1n values due to increasing con- Throughout the comparisons we show that the mass tribution from the FF decays. Thus, for the I and Cs dependence of total half-lives and Pxn is sensitive to a isotopic chains, deformation is still not a decisive factor structural indicator %FF that represents the ratio of the for the mass dependence of the P1n values. It can be seen FF (or GT) decays over the total rate: that for these isotopes the spherical approaches, namely DF3a+CQRPA and RHB+RQRPA, provide a level of ac- %FF = λFF/λ = (TGT − T )/TGT. (13)
24 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
100 DF3a +CQRPA DF3a +CQRPA (GT only) 80 FRDM12 +(Q)RPA. 10 RHB + RQRPA CRP 60 (%) (s)
1n 40 1/2 1 P T 20 Arsenic (Z=33) 0.1 0 84 85 86 87 88 84 85 86 87 88 A A 100 100
CRP 10 DF3 + CQRPA 80 RHB + RQRPA FRDM12 + (Q)RPA 60 (s) 1 (%) 1n 1/2 40 T P 0.1 20 Bromine (Z=35) 0.01 0 87 88 89 90 91 92 93 94 87 88 89 90 A 91 92 93 94 A 100 100
CRP 10 80 RHB+RQRPA 60
FRDM12+(Q)RPA (s) (%) 1 1n 1/2 40 P T 0.1 20 Rubidium (Z=37) 0.01 0 92 94 96 98 100 102 92 94 96 98 100 102 A A 1 100 CRP
DF3a + CQRPA 80 FRDM12 + (Q)RPA (s) 0.1 RHB + RQRPA
(%) 60 1/2 FAM 1n T 40 P 0.01 20 Molybdenum (Z=42) 0
110 112 114 116 118 120 122 124 126 110 112 114 116 118 120 122 124 126 A A 100 10 Df3a +CQRPA FRDM12 + (Q)RPA 80 RHB + RQRPA 60 1 CRP (%) 1n (s) RIKEN 2020 P
1/2 40 T 0.1 20 Iodine (Z=53) 0 0.01 138 140 142 144 146 148 138 140 142 144 146 148 A A 100 10 FRDM12 + (Q)RPA 80 DF3a + CQRPA RHB + RQRPA
(%) 60 (s)
1 CRP 1n 1/2
P 40 T 0.1 20 Cesium (Z=55) 0.01 0 142 144 146 148 150 152 142 144 146 148 150 152 A A
FIG. 9. (Left panels) Predictions for half-lives from various theoretical models for As, Br, Rb, Mo, I, and Cs isotopes compared with the evaluated CRP data. For 143−146I also the latest RIKEN data from Ref. [189] is shown. (Right panels) Predictions for P1n values for the same isotopes from various theoretical models compared with the evaluated CRP data. For more information, see text. 25 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
TABLE II. Evaluated and theoretical T1/2 and P1n values for bromine (Z= 35) and iodine (Z=55) isotopes. “FRDM+RPA” is the model labelled “QRPA-2” in Ref. [22]. ∗Latest data from Ref. [189]. CRP evaluation DF3a (this work) FRDM + RPA [22] FRDM12 + (Q)RPA + HF [148] RHB + RQRPA [153] T1/2 (s) P1n (%) T1/2 (s) P1n (%) T1/2 (s) P1n (%) T1/2 (s) P1n (%) T1/2 (s) P1n (%) 87Br 55.64(15) 2.53(10) 47.730 0.99 37.25 1.129 6.543 0.0 3.407 5.3 88Br 16.29(8) 6.72(27) 16.550 2.9 104.99 28.599 1.400 1.0 1.292 3.7 89Br 4.338(22) 13.7(6) 6.290 12.6 10.76 41.746 0.177 1.0 0.597 7.8 90Br 1.911(10) 25.6(15) 1.770 17.5 17.33 99.800 0.106 3.0 0.313 5.4 91Br 0.544(10) 29.8(8) 0.980 18.5 0.762 73.365 0.051 9.0 0.175 15.8 92Br 0.334(14) 33.1(25) 0.250 20.1 0.054 8.457 0.034 11.0 0.108 11.7 93Br 0.152(8) 64(7) 0.120 33.4 0.221 100 0.040 27.0 0.072 48.5 94Br 0.070(20) 30(10) 0.100 50.1 0.034 14.221 0.108 53.0 0.052 46.9 137I 24.59(10) 7.63(14) 29.28 11.8 3365.424 98.980 18.099 3.0 0.948 5.5 138I 6.251(31) 5.30(21) 5.1 1.5 9020.949 69.824 11.038 5.0 0.496 4.6 139I 2.280(11) 9.74(33) 0.87 3.95 58.145 99.988 2.133 13.0 0.299 10.3 140I 0.86(4) 7.88(43) 0.79 3.69 17.216 99.728 0.433 13.0 0.194 8.5 141I 0.43(2) 21.2(30) 0.46 19.5 2.347 100 0.428 51.0 0.131 29.3 142I 0.235(11)∗ - 0.401 5.61 1.400 99.952 0.074 19.0 0.096 27.0 143I 0.182(8)∗ - 0.319 11.1 0.150 77.120 0.128 70.0 0.072 37.2 144I 0.094(8)∗ - 0.198 10.4 0.058 29.379 0.049 36.0 0.055 34.9 145I 0.0897(93)∗ - 0.126 61.9 - - 0.076 82.0 0.043 39.1 146I 0.094(26) ∗ - 0.0762 62.7 - - 0.037 54.0 0.034 34.6
26 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
Here, λ and λ are the total rate and rate of the FF 1 FF 10 RHB+RQRPA transitions only, respectively, T and TGT are the total FAM half-life and the half-life of the GT transitions only, re- 0 DF3+CQRPA (s) 10 spectively. Naturally, %FF reflects selection rules and, FRDM+(Q)RPA in particular, a reduction of the unique FF compared to 1/2 -1 CRP the non-unique FF. T 10 10-2 Ni
78 1. β-decay properties in the Ni region 72 74 76 78 80 82 84 86 88
For the reference doubly-magic nucleus 78Ni (N= 50, 100 Z= 28) and the isotopes in the region (both with Z < 28 and Z ≥ 28), the high-energy (low Ex) GT transitions 10 mostly define the total half-lives at N < 50, while the FF decays give less-important contribution to the total % FF 1 β-decay rates. Eventually, for nuclei with occupied pro- Ni ton 1f7/2 orbital beyond the N = 50 neutron shell, the high transition-energy FF decays start to compete with 72 74 76 78 80 82 84 86 88 the GT decays, although the former continue to make A a substantial contribution to the total half-life as their transition energies increase with N − Z. FIG. 10. a) (Upper panel) Half-lives for the Ni iso- The total half-lives for the Ni chain obtained topes (Z= 28) calculated from DF3a+CQRPA (present from DF3a+CQRPA [151, 177], RHB+RQRPA [153], work), RHB+RQRPA [153], FAM (even-even only) [157], and FRDM12+(Q)RPA [148], and FAM [157] are compared FRDM12+(Q)RPA [148] compared with the evaluated half- with evaluated data in Fig. 10. The FRDM12+(Q)RPA life data [35]; b) (Lower panel) Contributions of the first- and DF3a+CQRPA results are close to the experimental forbidden transitions to the total decay rate (%FF value). data, though an odd-even staggering is observed in the FRDM12+(Q)RPA. For the spherical RHB+RQRPA as well as for the de- multi-neutron emission thresholds, and quasi-particle lev- formed FAM, the deviation from the evaluated half-lives els predicted by the different models are very important. at N < 50 is quite pronounced. For the doubly-magic nu- The different phase-(sub)-spaces defined by these quanti- cleus 78Ni, they predict 258 ms and 605 ms, respectively, ties are reflected in the β-strength distributions and are compared to the experimental value of T1/2=122.2(51) ms translated into the half-lives and delayed multi-neutron [35, 191]. emission probabilities. In order to understand the origin of these deviations, The Pxn for the Ni chain obtained from it is necessary to analyze the contribution of the FF (or the DF3a+CQRPA, RHB+RQRPA, and GT) decays to the total rate (%FF value, Fig. 10). For Ni FRDM12+(Q)RPA+HF models are compared in Fig. 11. isotopes with neutron numbers N ≤ 50, as the π1f7/2 or- It can be seen that in the RHB+RQRPA calculations, bital is filled, the GT decays (ν1f5/2, π1f5/2), (ν1p1/2,3/2, the peak of the one-neutron emission probability occurs 78 π1p1/2,3/2) overwhelm the hindered FF-unique decay at Ni, two mass units earlier compared to the other two (ν1g9/2, π1f5/2). Naturally, the %FF values obtained model calculations, and that the two-neutron emission from the DF3a+CQRPA for A < 78 are low and the probability mass behavior also differs. However, the 75−77 drop of the %FF values at A = 78 is mostly related to three experimental P1n data points for Ni are very the competition with the GT transitions. For the FAM nicely reproduced by the FRDM12+(Q)RPA+HF and and RHB+QRPA models, the calculated %FF values for the RHB+RQRPA model, whereas the DF3a+CQRPA A ≤ 78 are higher by 30 and 50%, respectively. However, model seems to underpredict these isotopes. high %FF values for A < 78 are not supported by the Finally, in Fig. 12 we compare the half-lives for the selection rules and the available decay schemes. Ni chain with the recently developed “beyond the pn- Thus, an overestimation of the half-lives for A < 78 by QRPA” models, such as the FRSA+PPC [163] which RHB+RQRPA and FAM may stem from the predicted takes into account the phonon-phonon coupling, the balance between the GT and FF strengths. For N > PVA models with the particle-vibration coupling (PVC) 50, the non-unique decays (2d5/2, 2s1/2→2p1/2,3/2, 1f5/2) [166, 167, 192], and the “quasi-particle time blocking ap- contribute to the decay process, yet the %FF values cal- proximation” (QTBA) [168]. Hybrid shell-model results culated from DF3a+CQRPA reach only 25% for A = 86, from [193] are also included in the figure. which is quite close to the FAM prediction. For A > The calculations in Fig. 12 have been performed within 78, all calculations give similar half-lives although the the GT approximation which is more justified for A < RHB+QRPA gives different %FF values. 78. In the pnQRPA+PVC developed in the recent pa- It turns out that the differences in the total Qβ values, per [167], the experimental half-lives for the Ni isotopes
27 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
80 P FRDM12 + (Q)RPA 1n QTBA 60 P 2n (%) SM
4n P 40 3n 0 DF3+spread P 10 4n - P SM+FF
1n P CRP 20 1n FRSA-PPC P
0 (s) HFB+QPVC CRP 72 74 76 78 80 82 84 86 1/2 T
80 -1 P 10 1n DF3 + CQRPA 60 P 2n (%) P 4n 40 3n P 4n
- P P CRP 20 1n
1n Ni P 0 10-2
72 74 76 78 80 82 84 86 72 74 76 78 80 82 84 80 A P 1n RHB +RQRPA
(%) 60 FIG. 12. Half-lives for Ni isotopes calculated from QTBA
P 2n
4n P [168], FRSA+QRPA+PPC [165] and HFB+QRPA+PVC 40 3n P [167] in comparison with with the evaluated data from 4n - P P CRP Ref. [35]. Also shown are calculations from the SM [161], 1n 20 1n NuShellX [193] and DF3a+CQRPA with spreading [177]. P 0 72 74 76 78 80 82 84 86 A brid” models, i.e. DF3a+CQRPA including the complex configurations through the spreading width [177] and the FIG. 11. Delayed multi-neutron emission probabili- shell-model NuShellX with the FF decays added from the ties (P1n–P4n) for Ni (Z = 28) isotopes calculated DF3a+CQRPA [193], can give an overall reasonable de- with FRDM12+(Q)RPA+HF [148], DF3a+CQRPA (present work), and RHB+RQRPA [153]. These models are compared scription of the experimental half-lives in Ni for A > 78. 75−77 Only a simultaneous account of the FF decays and np- to the evaluated P1n data for Ni as published in Ref. [35]. nh configurations in these models allows one to describe the “sudden shortening” of the half-lives after crossing are reasonably well described assuming a rather strong the N = 50 shell that was observed in the RIKEN T = 0 pairing (except for A = 80). The correspond- data [191]. It is therefore reasonable to expect that ing 1p-1h QRPA calculation in the same framework has the self-consistent ‘beyond pnQRPA’ models have the given a factor of about 10-15 longer half-lives [167]. Thus potential to describe all these effects. In this respect, the renormalization due to the PVC turns out to be very it is important to include the contribution of the first- strong leading to QRPA+PVC half-lives that are only by forbidden decays in the present (spherical) ‘beyond pn- a factor 1.5 shorter for 72Ni and by a factor 1.5 longer for QRPA’ schemes. 78Ni. On the other hand, an effort to include the contribu- A similar strong renormalization due to the PVC ef- tion of the complex configurations into the deformed Fi- fect is obtained in the QTBA approach [168]. Again, nite Amplitude Method (pnFAM) [157] would result in the pure pnQRPA half-lives obtained without assum- a universal approach that could be applied to the whole ing phonon-phonon coupling in the QTBA are much nuclear chart. higher than the corresponding quantities obtained from the FRSA [165] which includes phonon-phonon coupling. 132 Thus, the renormalization factor due to the effective in- 2. β-decay properties in the Sn region teractions in these approaches differs significantly (by about 1–2 orders of magnitude). Most probably this dif- For the isotopes beyond the major neutron shell N = ference is due to the fact that the QTBA model uses ex- 82 in the 132Sn region, the concurrence of GT and FF 78 perimental Qβ values instead of self-consistently derived decays has much in common with the Ni region. We ones. show that the contribution of GT and FF transitions to The above discussion is important if one considers the the β-decay rates differs for the Z < 50 and Z ≥ 50 iso- potential of these models: once the contribution of FF topes. The intensive GT decays in the Z < 50 nuclei decays is included, they could be implemented in global mostly correspond to the (ν1g7/2, π1g9/2) configuration. calculations of β-decay half-lives and β-delayed neutron The high-energy GT decays contribute strongly to the branching ratios. It has already been shown that “hy- total half-lives of the nuclei with Z < 50, N < 82. After
28 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. crossing the N = 82 neutron shell, the high-energy FF 100 decays which are mainly related to the (ν1h11/2, π1g9/2) Tin configuration are active. In contrast to what is observed in nuclei in the Z ≥ 28 region, the high-energy GT transi- 80 tions for Z > 50 nuclei are hindered, as the p1g9/2 orbital is blocked. At the same time, the high energy FF transi- tions (ν2f7/2, π1g7/2), (ν1f7/2, π2d5/2) come into effect 60 and mostly define the total half-life. Fig. 13 shows the comparison of the half-lives obtained from DF3a+CQRPA, RHB+RQRPA, % FF 40 FRDM12+(Q)RPA, and FAM with the evaluated CRP data. For completeness, we include 132Sn in the plot even though it is not a βn-emitter . Within the 20 DF3a + CQRPA GT-only approximation, the DF3a half-lives are up to FAM a factor 5 longer than the data. Taking into account RHB+RQRPA the FF decays sets the balance between the GT and FF 0 strengths as shown in Fig. 14 and allows for an excellent 132 134 136 138 140 description of the experimental half-lives for the Sn A isotopic chain. A rather high percentage of FF decays of 60–75% calcu- FIG. 14. Contributions of the FF transitions to the to- lated from the RHB+RQRPA may be the reason for the tal decay rate (%FF values) for Sn isotopes obtained from strong underestimation of the half-lives of 133−140Sn by a DF3a+CQRPA (present work), RHB+RQRPA [153], and factor of 4–5. On the other hand the FRDM12+(Q)RPA FAM [157]. half-lives overestimate the experimental data by a factor of ≈ 2 while the DF3a+CQRPA (GT+FF) calculations are quite close to the experimental values. observing 99% intensity for the decay to the GT state at 100 1.325 MeV with a log(ft) = 4.2. Accordingly, the DF3a DF3a+CQRPA (GT+FF) calculates a half-life of 35.7 s, compared to the evaluated DF3a +CQRPA (GT only) value of 39.7(5) s from ENSDF. In contrast, the RHB+ 10 RHB+RQRPA FRDM12+(Q)RPA RQRPA predicts 610 s and the FAM a value of 820 s. CRP The corresponding P1n values (Fig. 13) obtained with (s) 1 RHB+RQRPA and DF3a+CQRPA agree in general with 1/2
T the mass dependence of the %FF values predicted by 0.1 Tin these models, however they underestimate the CRP data for A > 133. The FRDM12+(Q)RPA+HF [148] overes- 0.01 timates the P1n values considerably which is consistent 100 with a gradually increasing share of the FF decays at N > 82 due to the increasing contribution of the ∆J = 0 80 g.s.-to-g.s. transitions at N = 83 − 86. The stabilization 60 of the P1n values at N > 86 predicted by RHB+RQRPA (%) and DF3a+CQRPA reflects a decrease of the %FF val- 1n 40 ues for N = 87 − 90 (137−140Sn) due to the competition P from the open high-transition energy GT decays. 20 As was already mentioned, the β-decay characteris- 0 tics of the semi-magic Ni and Sn isotopic chains are 132 134 136 138 140 highly sensitive to the details of the nuclear structure A aspects of the models, in particular the energy-density functional used. As we move further away from the
FIG. 13. Half-lives (upper panel) and P1n values (lower closed shells, these characteristics become less sensitive panel) for Sn isotopes (Z = 50). Calculated values from to the details of the models. For example in the Pd DF3a+CQRPA [177], RHB+RQRPA [153], FAM [157], and and Cd isotopes that have pairing in both neutron and FRDM12+(Q)RPA [148] are compared with the evaluated proton sectors there is good agreement between the self- 132 (T1/2,P1n) data. The half-life for Sn (no βn emitter, consistent DF3a+CQRPA and RHB+RQRPA calcula- 39.7(8) s) was taken from ENSDF. tions [153, 176, 177]. To summarize, from the comparison of theory and Note that for the reference doubly-magic nucleus 132Sn evaluated data for a fairly broad range of nuclei shown the very low %FF value predicted by DF3a+CQRPA in this section, the recommended theoretical approach calculations (Fig. 14) is consistent with the decay scheme is the quasiparticle random phase approximation (pn-
29 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
QRPA) based on the energy density functional (EDF) ap- ongoing large-scale experimental campaigns at RIB facil- proach. Recently developed self-consistent models have ities worldwide (see, e.g. Refs. [125, 189, 191, 194]). been successfully used in this CRP to calculate the half- In Fig. 15, the ratio of half-lives and P1n values from lives and β-delayed neutron rates of fission products FRDM12+(Q)RPA [148] and the evaluated data are plot- and nuclei near the major closed shells that are impor- ted with respect to the evaluated values. It is obvi- tant for r-process modeling. Compared to the old semi- ous that the half-lives and neutron-emission probabilities empirical models (“Gross theory” versions) and the semi- of almost all βn-emitter considered here are reproduced microscopic FRDM-based (Q)RPA framework for the GT within one order of magnitude. decays that is augmented by a statistical model for the In the case of half-lives, the majority is reproduced FF decays, the self-consistent β-decay models have the within a factor of two or better, denoted by the black obvious advantage of being well-founded on first princi- dashed lines. It is also noticeable that the scatter around ples. By definition, the parameters of the self-consistent the ratio of 1 decreases with decreasing experimental half- global models should be kept the same when extrapolat- life, which is a common pattern. Shorter half-lives are ing to different mass regions across the nuclear chart. characterized by larger Q-values, and a larger portion of Global models based on the energy density functional the decay strength is included within the Qβ window. approach are not only used in β-decay studies but also for As a result the half-lives of shorter-lived nuclei are less a consistent description of a range of nuclear properties sensitive to the deviations of the model from the true in a broad mass region, such as binding energies, radii, strength function. magnetic moments, fission barriers as well as nuclear re- In the case of the comparison with experimental P1n sponses to different probes. Consequently, there is a wide values in the bottom part of Fig. 15 several nuclei can be variety of data available to constrain the EDF parame- found aligned on diagonal lines. This surprising feature ters and one can avoid using the “output” half-lives and is due to theoretical P1n values that appear only as inte- Pn values for such a purpose. Self-consistent global mod- ger values in the FRDM12+(Q)RPA+HF model. Lines els are thus a reliable and universal instrument for the for P1n= 1,2,3,4, and 5% are shown in red to guide the description of the middle-heavy, heavy and superheavy eye. For cases where the experimental values are very nuclei which are in the scope of the present paper. small, e.g. lower than 1%, any deviation of the model Finally, it is worth mentioning the rapidly developing predictions from the true energies and strengths of the “beyond the pnQRPA” models such as FRSA+phonon- transitions will, as a consequence, result in relatively large phonon coupling [163, 164], the PVA models with deviations of the total emission probabilities from the ob- particle-vibration coupling [192] and the “quasi-particle served values. time blocking approximation” (QTBA) [168]. These ap- Figure 16 shows the same as Fig. 15, but for theo- proaches can potentially describe the β-decay strength retical results obtained using the RHB+RQRPA model function over a wide range of masses across the nu- [153] which was extended for nucleosynthesis calcula- clear chart with a quality comparable to the multi- tions and for use in comparisons in this CRP. The configurational shell-model. For example, the “sudden general behaviour of the model is quite similar to the acceleration” of the β-decay after crossing the N = 50 FRDM12+(Q)RPA+HF shown above, but with some no- shell, that was recently discovered in the Ni isotopic chain table differences. In the case of half-lives, the relativistic at RIKEN [191], has been described by assuming np-nh approach manages to provide, on average, a somewhat configurations within the DF3a+CQRPA approach using better description for half-lives <200 ms as indicated by quasi-particle spreading. This confirms that such a phe- the smaller scatter of ratios around the value of 1. How- nomenon can only be explained by considering FF decays ever, for longer half-lives >10 s the model seems to sys- and complex configurations simultaneously. tematically underestimate the experimental half-lives. In the case of neutron emission probabilities, the same diag- onal trend as for the FRDM12+(Q)RPA+HF model in D. Global comparisons of theoretical results Fig. 15 can be seen, due to the fact that many nuclei have similar theoretical P1n values around 1–2%. Figures 17 and 18 present the same results as in the pre- Global comparisons between the CRP evaluated vious two figures, but now the ratios are plotted versus data and the self-consistent models DF3a+CQRPA, the atomic mass number A, in order to explore the impact RHB+QRPA, as well as the microscopic-macroscopic of shells on the results. The FRDM12+(Q)RPA+HF for- model FRDM12+(Q)RPA(+HF) are presented in this malism provides a very robust description of the data, section. Large-scale calculations of the β-decay energy without showing any visible strong dependence of the releases (Qβ, Qβn), half-lives (T1/2) and delayed neutron model predictions on the atomic mass number (except (DN) emission probabilities (Pn) for hundreds or even for a notable dip in the deformed A ≈ 90 − 95 region). thousands of spherical, near-spherical and deformed nu- This is true both for the half-lives and for the one-neutron clei have been performed so far within several models. emission probabilities. Nowadays these models provide a reliable data input for However, the results obtained with the relativistic the modeling of abundances in the astrophysical r pro- RHB+RQRPA model (Fig. 18) show significant effects cess. They have also been used as predictions for the
30 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
100 FRDM12 + (Q)RPA + HF
p 10 x e 2 / 1 T
/ 1 o e h t 2 / 1
T 0.1
0.01 -3 -2 -1 0 1 2 10 10 10 exp 10 10 10 T (s) 1/2 100 P theo= 5% 1n P theo= 1% 1n 10 p x e n 1 1 P / o e h t n
1 0.1 P FRDM12 + (Q)RPA + HF 0.01 0.01 0.1 1 10 100 P exp (%) 1n
FIG. 15. Ratios of theoretical and experimental results for the half-lives (upper panel) and P1n values (lower panel) calculated using the FRDM12+(Q)RPA+HF formalism [148], plotted versus the corresponding experimental values. Diagonally aligned nuclei in the lower plot are nuclei that have the same theoretical P1n value. Lines for P1n= 1, 2, 3, 4, and 5% are shown to guide the eye. of the nuclear shell structure on the quality of the re- (near-) spherical nuclei with Z = 25−35 and Z = 44−55 sults in the description of the β-decay half-lives. These regions with T1/2 ranging from 10 ms to 10 s. This results display a clear arch between two mass regions, constraint is a consequence of the strength function for- from A ≈ 90 to A ≈ 140, which is consistent with the malism which assumes that QRPA results are better for fact that the model does not take nuclear deformations short-lived nuclei with relatively large Qβ energy window. into account. Therefore, the model produces deviations Also, the limitation of the comparison to spherical and in the regions of the nuclear chart between closed shells near-spherical nuclei makes it easier to understand the where nuclei have non-negligible ground-state deforma- origin of the deviations from the experimental half-lives. tions. However, the β-decay half-lives are well reproduced Fig. 19 shows that most of the half-lives are in agreement within a factor of 2 (although longer than measured) for within a factor of two with the evaluated data. the mass region A = 100 − 130. In the lower panel of Fig. 19, the ratio of the theoret- In the upper panel of Fig. 19, the half-lives calculated ical and experimental half-lives are plotted as a function within the spherical DF3a+CQRPA approach are plotted of mass number A. The results indicate that the data is with respect to the evaluated values. The figure includes in general reproduced within a factor of 2 over the whole
31 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
100 p x
e 10 2 / 1 T / o
e 1 h t 2 / 1 T 0.1
RHB+RQRPA 0.01 10-3 10-2 10-1 100 101 102 T exp (s) 1/2 100 P theo= 1% 1n
p 10 x
e theo
n P = 0.2% 1n 1 P
/ 1 o e h t n 1
P 0.1 RHB+RQRPA 0.01 0.01 0.1 1 10 100 exp P (%) 1n
FIG. 16. Ratios of theoretical and experimental results for the half-lives (upper panel) and P1n values (lower panel) calculated using the relativistic RHB+RQRPA formalism [153], plotted versus the corresponding experimental values. Diagonally aligned nuclei in the lower plot are nuclei that have the same theoretical P1n value. Lines for P1n= 0.2 and 1% are shown to guide the eye.
mass region, with no obvious discrepancies at shell clo- underestimation of the P1n values. sures or for deformed nuclei around A ≈85 and 120. It is worth stressing that the β-decay strength function In a number of cases, however, such as in K and Sb can be modified by including the quasiparticle-phonon isotopes, the deviation has been shown to be related to coupling. The sensitivity of the neutron emission prob- ground-state spin inversion rather than to deformation abilities to the strength distribution in the vicinity of [177, 183]. the neutron emission thresholds Sxn is extremely high. Similarly, the ratio of neutron emission probabilities of Consequently, in the one-particle-one-hole pnQRPA ap- (near-)spherical nuclei is shown in Fig. 20 as a function of proach, the accuracy of global predictions for the P1n A. Overall, the experimental data are reproduced within values is lower than for half-lives. the same factor of two. Again the exceptions to this trend In Fig. 21, the ratio of DF3a+CQRPA predictions to are weakly deformed nuclei in the Z = 31 − 35, A ≈ 85 experimental half-lives [125, 126] is shown for the heav- and Z = 44−45, A ≈ 120 regions which have rather small ier Os to Bi nuclei (Z = 76–83). As one approaches the P1n values. In these cases, neglecting the spreading of N = 126 shell closure, nuclei are either spherical or have the β strength due to deformation leads to a systematic a small ground-state deformation. For the bulk of the
32 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
100 FRDM12 + (Q)RPA + HF
p 10 x e 2 / 1 T
/ 1 o e h t 2 / 1
T 0.1
0.01 20 40 60 80 100 120 140 160 180 200 220 240 A 100
p 10 x e n 1 P
/ 1 o e h t n 1
P 0.1
FRDM12 + (Q)RPA + HF 0.01 20 40 60 80 100 120 140 160 180 200 220 240 A
FIG. 17. Ratios of theoretical and experimental results for the half-lives (upper panel) and total P1n values (lower panel) calculated using the FRDM12 +(Q)RPA +HF formalism [148], plotted versus the atomic mass number. nuclei included in the DF3 calculations, we observe a decay half-lives have also been performed within the de- larger spread of the half-lives of up to a factor of 5. For formed self-consistent QRPA model of Ref. [174]. The Hg isotopes around A = 210 with relatively low Qβ val- model uses the GT approximation and a global D1M ues, the DF3a+CQRPA underestimates the half-lives by Gogny interaction. Only the half-lives of 145 spherical one and RHB+RQRPA by two orders of magnitude, re- and deformed even-even nuclei for which the experimen- spectively (see Ref. [125]). Both DF3a+CQRPA and tal half-lives are known were included in the comparison. theo exp RHB+RQRPA underestimate the experimental half-life The calculated ratio T1/2 /T1/2 is quite homogeneous of 206Au [195, 196] by a factor of about 30 – 60. with respect to the mass number A and, more partic- A comparison of the existing eight experimental P1n ularly, the quadrupole deformation parameter β2. The values for nuclei beyond the N = 126 shell closure [125, treatment of deformation reduces the overall deviation 196] is shown in Fig. 22. Note, however, that these are the from the data to basically the same factor of 10 which was heaviest nuclei for which P1n values have been measured obtained for spherical nuclei in Refs. [151, 153]. Including so far. As can be seen, the experimental uncertainties are deformation in the model is therefore expected to lead to still too large for any detailed comparison of the different a marked improvement in the RHB+RQRPA half-lives model performances. for deformed nuclei and similar levels of accuracy that It is worth mentioning that global comparisons of β- were obtained for spherical nuclei in Refs. [151, 153].
33 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
100 RHB+RQRPA p x
e 10 2 / 1 T /
o 1 e h t 2 / 1
T 0.1
0.01
20 40 60 80 100 120 140 160 180 200 220 240
A 100 p
x 10 e n 1 P
/ 1 o e h t n 1 0.1 P
RHB+RQRPA 0.01
20 40 60 80 100 120 140 160 180 200 220 240 A
FIG. 18. Ratios of theoretical and experimental results for the half-lives (upper panel) and P1n values (lower panel) calculated using the relativistic RHB+RQRPA formalism [153], plotted versus the atomic mass number.
As far as the comparison of different microscopic mod- From the global comparisons we deduce that els with the CRP data is concerned, we observe that re- DF3a+CQRPA [151], RHB+RQRPA [153], and gardless of the treatment of deformation all models using FRDM+(Q)RPA(+HF) [148] models give overall a a concept of the β-decay strength function have difficulty comparable description of the β-decay properties. Note describing nuclei close to the valley of β stability with that a somewhat smaller spread in the DF3a+CQRPA low Qβ values. The RHB+RQRPA calculations, for in- half-life calculations (Fig. 19) with respect to the stance, tend to underestimate the half-lives for very heavy RHB+RQRPA (Fig. 18) is partially a result of the nuclei, around A ≈ 210 by more than an order of magni- limited comparison of (quasi-)spherical-only nuclei used tude [125, 126]. This generic feature is also true for the in the DF3a calculations. An important source of microscopic-macroscopic model and “Gross theory”. deviations from the experimental data may come from Neutron-rich nuclei near the closed shells are the main the description of the ground-state properties. playground for developments in the microscopic theories. In particular, the underestimation of the total decay The theoretical predictions of their β-decay properties – phase-space Qβ, as well as Qβxn sub-spaces may have a while being a crucial test for any nuclear structure theory destructive impact on the predicted β-decay observables. – are also critical for planning experiments at RIB facil- This is especially true for the Pxn-values, as these ob- ities. Additionally, these nuclei are also key for a better servables are more sensitive to the different ”systematic” understanding of the astrophysical r process. uncertainties of the theoretical strength functions.
34 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
10
Z=18 - 55, T1/2<10 s exp 1/2
T 1 /
theo 1/2 T 0.1 0.01 0.1 1 10 T exp (s) 1/2 10 exp 1/2 T
/ 1
theo
1/2 T 0.1 40 60 80 100 120 140 160 A
FIG. 19. Ratio of theoretical and experimental values of half-lives of (quasi-) spherical nuclei with Z = 18 − 21, 25 − 35, and 44−55, calculated using the DF3a+CQRPA formalism and plotted versus the evaluated experimental half-lives in the top panel and versus mass A in the bottom panel.
Finally, with regards to the predictive power of the tential [199] and the generalized Lorentzian E1 strength models, if we consider the smaller number of parameters function [200] were adopted, respectively, with a parame- involved in the self-consistent models, as well as their uni- ter set obtained from the Reference Input Parameter Li- versality, we can conclude that the predictions of these brary (RIPL-2) [201]. For the nuclear level densities, they models are more reliable than those of the “Gross the- used RIPL-2 for discrete levels and the Gilbert-Cameron ory” and FRDM-based models whose parameters are de- level density formula [202] at high energies. Spin and par- termined near the stability line. ity selection rules for the decay of the precursor nuclei to daughter nuclei were observed. The microscopic-macroscopic approach E. Delayed neutron spectra calculations (FRDM+(Q)RPA+HF) to calculating β-delayed neutron spectra is important for various applications as it can predict the spectra for nuclei located on the neutron-rich The microscopic-macroscopic approach described in side of the nuclear chart for which experimental data the previous sections was also used in the calculations are not available. Furthermore, as has been discussed of delayed neutron (DN) spectra of fission products in in Sect.IV, it provides an improved description of the Ref. [32]. The FRDM95 [186] was combined with the decay process leading to DN emission as it considers (Q)RPA model [147, 198] to estimate β-strength func- nuclear structure effects such as nuclear shell effects, tions which were then used as input in the Hauser- level structures, and the spin-parity selection rules. The Feshbach statistical model (HF) to obtain DN branch- β-delayed neutron spectra calculated by [32] are now ing ratios and spectra. The HF calculations require sev- adopted in the ENDF/B-VIII.-0 decay-data library [25]. eral nuclear ingredients. For neutron and γ-ray transmis- Within the CRP, we performed similar calculations of sion coefficients, the Koning-Delaroche global optical po-
35 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. 10
Z=18 - 55, T1/2<10 s exp
n 1 / P 1
theo 1n P 0.1 0.01 0.1 1 10 100 exp P1n (%) 10 exp
1n / P 1
theo 1n P 0.1 40 60 80 100 120 140 160 A
FIG. 20. Ratio of theoretical and experimental values of β-delayed one-neutron emission probabilities (P1n) of (quasi-)spherical nuclei with Z = 18 − 21, 25 − 35 and 44 − 55, calculated using the DF3a+CQRPA formalism and plotted versus the evaluated experimental P1n values in the top panel and the mass A in the bottom panel.
DN spectra by using a Hartree-Fock-BCS model (SHF- considered. The HF statistical model implemented in the BCS) – instead of the FRDM – combined with QRPA code CCONE [204] was used to calculate the competition to obtain the β-decay properties. The DN emission between DN and γ decays. The global neutron optical po- part was calculated using the HF model [170]. In tential of Ref. [199] and enhanced generalized Lorentzian this SHFBCS+QRPA+HF model, the pairing correla- (EGLO) model [200, 205] were used to calculate the neu- tion is treated by the BCS approximation with a zero- tron and γ transmission coefficients, respectively. For range volume-type force. The ground and excited states level densities, we used RIPL-3 [205] for discrete levels were self-consistently calculated with the SkO’ Skyrme and the formula of Gilbert and Cameron with Mengoni– force [203] in the cylindrical coordinate, where the ax- Nakajima parameter set [206] for higher energies. ial symmetry was assumed. A zero-range volume-type Figure 23 shows the β-delayed neutron spec- pairing was used in the pp-channel (T = 0 pairing), and tra obtained for Rb and Cs isotopes from SHF- the strength parameters were adjusted to obtain the best BCS+QRPA+HF and FRDM+(Q)RPA+HF. The re- possible description of available half-life data for each iso- sults are compared with the evaluated data of ENDF/B- topic chain. Only allowed GT transitions are considered VIII.0 [25] and the experimental data obtained from in the β-decay channel in Ref. [170]. Odd-mass nuclei Ref. [31] that were described in Sect. IIIE1. were treated in the same way as the even-mass nuclei, For FRDM+(Q)RPA+HF, the β-delayed neutron spec- and the blocking and polarization effects caused by cou- tra of 94Rb and 145Cs deviate from experimental data, pling between valence particle(s) and core nuclei were not while those of 95Rb and 143,144,146Cs are reproduced ap-
36 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
Z= 76 - 83, T < 500 s 10 1/2 exp
1/2 1 T
/ th 0.1 1/2 T 0.01 1 10 100 1000 T exp (s) 1/2
exp 10 1/2 T
/ 1 th
1/2
T 0.1
0.01 200 210 220 A
FIG. 21. Same as in Fig. 19 but for heavy (Z = 76 − 83, Os to Bi) spherical isotopes calculated using the DF3a+CQRPA formalism [125, 126, 197]. The experimental data is taken from the recent evaluation [36].
50 while those of the other nuclei are approximately repro- 45 duced. It is worth noting that some of the β-delayed DF3+CQRPA neutron spectra are reproduced reasonably well by these 40 RHB+RQPA FRDM+(Q)RPA theoretical models, even though QRPA approaches are 35 CRP unable to estimate the nuclear excited states within an 30 uncertainty of a few hundreds keV. The deviations from
(%) experimental data are mainly attributed to the calculated 25 β-strength function since both approaches use similar HF 1n
P 20 implementations. For example, the β-delayed neutron spectrum of 145Cs is reasonably reproduced with SHF- 15 BCS+QRPA+HF but not with FRDM+(Q)RPA+HF, 10 which implies that the β-strength function of SHF-
5 BCS+QRPA+HF is more suitable for this nucleus than the FRDM+(Q)RPA+HF strength function. This con- 0 210Hg 211Hg 211Tl 212Tl 213Tl 214Tl 215Tl 216Tl clusion is based on the assumption that the particle evap- oration process is correctly described. FIG. 22. Theoretical and experimental values for the β- Similar systematic detailed comparisons of experimen- delayed one-neutron emission probabilities of heavy (quasi-) tal and calculated DN spectra that duly consider the un- spherical Hg and Tl isotopes (Z = 80, 81), calculated using certainties in the measured spectra and the ingredients the DF3a+CQRPA [125, 126, 197], RHB+RQRPA [153], and of the models are expected to lead to further improve- FRDM12+(Q)RPA+HF models [148]. The experimental data ments of the theoretical β-decay models in the future. is taken from the evaluation in Ref. [36]. However, as already mentioned in Sect. IIIE1, it is still not clear how to quantify the experimental uncertainties in the measured spectra of Ref. [31], therefore, until that proximately. The SHFBCS+QRPA+HF model cannot is possible, caution should be exercised when comparing describe the β-delayed neutron spectra of 94Rb and 144Cs, with these experimental DN spectra.
37 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
The published and digitized experimental DN spectra V. MACROSCOPIC DATA: METHODS AND are available on the IAEA online database [43]. MEASUREMENTS
) In the following sections we deal with macroscopic 4 Exp. -1 ENDF/B-VIII.-0 delayed neutron data produced as a result of neutron- 3 SHFBCS+QRPA+HF induced fission of major and minor actinides. These data FRDM+(Q)RPA+HF are of utmost importance for studies on kinetic response 2 and safe operation of fission reactors.
1
DN Spectrum (MeV (a) Rb-94 A. Methods
) 4 -1
3 1. Total delayed neutron yields and time-dependent parameters 2 The direct methods used for measuring the macro- 1 scopic delayed-neutron characteristics from neutron-
DN Spectrum (MeV (b) Rb-95 induced fission of heavy elements are based on irradiat-
) 4 ing the investigated sample by neutrons and subsequently -1 measuring the appropriate delayed-neutron observables: 3 the time dependence of delayed-neutron activities or the delayed-neutron energy spectrum. Neutron sources used 2 for this purpose are mainly the neutron flux in nuclear 1 reactors [207–214], critical assemblies [13, 215, 216], or electrostatic accelerator-based nuclear reactions: Li(p,n),
DN Spectrum (MeV (c) Cs-143 T(p,n), T(d,n), D(d,n) [217–223]. All experimental meth- ) 4 -1 ods are based on the assumption that the delayed-neutron 3 activity arising from the irradiation of the fissionable sample can be represented by the linear superposition 2 of the exponential decay modes, each with its own decay constant 1 m
DN Spectrum (MeV (d) Cs-144 dN(t) X −λi·tirr −λi·t = A (1 − e ) · ai · e , (14) ) 4 dt -1 i=1
3 A = n · σf · φ · Nf · νd ,
2 where A is the saturation delayed-neutron activity, ai and λi the relative abundance and decay constant of the i-th 1 delayed neutron group, respectively, tirr the irradiation
DN Spectrum (MeV (e) Cs-145 time, m the number of delayed neutron groups each of
) 4 which is associated with a group of delayed-neutron pre- -1 cursors with similar decay constants, n the efficiency of 3 the neutron detector, νd the total delayed-neutron yield 2 2 per one fission, φ the neutron flux (1/cm s), σf the fis- sion cross section (barn), Nf the number of atoms in a 1 fissionable sample. (f) Cs-146 A comprehensive study of the delayed-neutron emis- DN Spectrum (MeV 0 sion properties conducted by Keepin et al. [13] on 232Th, 0.0 0.5 1.0 1.5 2.0 2.5 233U, 235U, 238U, 239Pu, and 240Pu lead to the devel- DN Energy [MeV] opment of a 6-group model of the temporary delayed- neutron characteristics (ai,Ti), where Ti = ln2/λi. This 94,95 FIG. 23. β-delayed neutron spectra for Rb and model has been successfully used by the reactor commu- 143−146 Cs calculated with SHFBCS+QRPA+HF [170] and nity for many years. Only recently, an alternative 8-group FRDM+(Q)RPA+HF [32]. The results are compared with model based on a set of decay constants that are univer- experimental data available in [31] and the evaluated spectra sal for all fissioning systems was proposed and justified in ENDF/B-VIII.0 [24]. by Spriggs et al. [224, 225]. In the case of instantaneous irradiation (prompt burst
38 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
irradiation, λitirr << 1), Eq. (14) becomes by the integration of Eq. (15)
dN(t) m X −λi·t Z tk+1 = A · tirr λi · ai · e . (15) dN(t) dN(tk) dt dt = · ∆tk = i=1 tk dt dt m From Eqs. (14) and (15) it becomes clear that the ac- X −λi·∆tk −λi·tk A · tirr ai · (1 − e ) · e + B · ∆tk . curacy of the delayed-neutron parameters νd and (ai,Ti) i=1 extracted from measured data depends to a large extent (17) on parameters such as the accumulated statistics, the fis- sion rate in the sample Rs = σf · φ · Nf , the transporta- tion time of the sample from the irradiation position to a It is this type of irradiation that was used by Keepin neutron detector and the neutron background. The qual- et al. [13] for the measurement of the decay curves which ity of the obtained delayed-neutron data is also related were used for the estimation of short-lived group param- to the incident neutron energy because the temporary eters (from 3 through 6) and the determination of the delayed-neutron parameters (ai,Ti) are known to have a total delayed neutron yield. The Godiva reactor (a bare noticeable energy dependence when expressed in terms spherical 235U assembly) was used as neutron source. In of the average half-life of delayed neutron precursors hT i addition, the total delayed-neutron yield was determined ([226], [227]). with the help of the following approximation [13]. a. Nuclear reactor-based methods. The most impor- tant advantage of the nuclear reactor-based method is t2 X dN(tk) the high intensity of the neutron flux. With such a high- · ∆t − B · (t − t ) = intensity flux it is possible to use very small samples of dt k 2 1 t1 fissionable nuclides (up to micrograms) as was demon- m strated by Waldo et al. [207]. However, the long time X −λi·t1 −λi·t2 A · tirr ai · e − e = it takes to transport the sample in the reactor measure- i=1 ments does not allow one to resolve all 6 groups of delayed m X neutrons. Furthermore, given that the energy of the pri- A · tirr ai = νd · · Rs · tirr , (18) mary neutrons can take a wide range of values or maybe i=1 is unspecified in many of the reactor experiments, it is not straightforward to use these data in the evaluation of where t is the time of the neutron pulse, R · t the delayed-neutron parameters. irr s irr total amount of fissions in a sample, and the sum of the As a rule, a single irradiation procedure is used in the relative abundances is equal to 1. nuclear reactor-based experiments. The time dependence b. Accelerator-based delayed neutron experiments. of the delayed neutron activity registered by a multichan- The neutron sources used in this method are the monoen- nel analyzer can be obtained by integration of Eq. (14) ergetic neutron fluxes from the accelerator-based nuclear ([222, 226]. reactions Li(p,n), T(p,n), T(d,n)4He, D(d,n)3He. The Z tk+1 dN(t) dN(t ) main advantage of the accelerator-based method is the dt = k · ∆t = possibility to investigate the energy dependence of de- dt dt k tk layed neutron properties in a wide range of incident neu- m a tron energies from thermal to 20 MeV [217–223, 231]. X −λi·tirr i −λi·∆tk −λi·tk A · (1 − e ) · · (1 − e ) · e Eqs. (14) and (15) reflect one single cycle of measure- λi i=1 ment which includes the irradiation of the sample, the de- + B · ∆tk , (16) cay, and the counting of delayed-neutron activity. How- ever, in accelerator-based experiments a sample under- where (dN(tk)/dt) · ∆tk is the number of counts reg- goes a number of cycles with the purpose of increasing istered in the time channel tk with width ∆tk, B the the statistics. The general equation for the determina- intensity of neutron background. This decay curve is tion of the total delayed-neutron yields and temporal used to estimate the temporal delayed-neutron parame- delayed-neutron characteristics (ai,Ti) on the basis of the ters (ai,Ti), the saturation activity A, the intensity of the delayed-neutron decay curve accumulated during n cycles background B, and the covariance matrix of the estimated is given by ([222], [232], [233], [226], [234], [230]) parameters with the help of the least squares method (LSM). The obtained delayed-neutron decay curve, the dN(t ) values of the estimated parameters (ai,Ti) and the back- k · ∆t = ground B are used to calculate the total delayed neutron dt k yield [212, 226, 228–230]. m X ai A · T · · (1 − e−λi·∆tk ) · e−λi·tk + B · ∆t , One possible type of irradiation is the instantaneous hi λ k (prompt burst) irradiation (see Eq. (15)). In this case, a i=1 i measured delayed-neutron decay curve can be obtained (19)
39 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
where Thi determines the history of irradiation the equation for the non-cumulative DN decay curve measured after the n-th cycle of irradiation-counting se- −λi·tirr Thi = (1 − e )× quences ! n 1 − e−n·λi·T −λi·T − e · , dN(tk) 1 − e−λi·T (1 − e−λi·T )2 · ∆t = dt k A = · σ · φ · N · ν , m f f d X ai A · T · · (1 − e−λi·∆tk ) · e−λi·tk + B · ∆t , hi λ k with n the number of cycles, T the period of one cycle i=1 i (irradiation, decay, counting, delay). (20) This equation is used to estimate the values of the sat- uration activity A, the temporal delayed-neutron param- −n·λ ·T eters (ai,Ti(i = 1, ..., m)) and the background B on the 1 − e i m −λi·tirr Thi = (1 − e ) · basis of the decay curves measured with different irra- 1 − e−λi·Tm diation times tirr. The A and B values and group pa- rameters of delayed neutrons are estimated using the it- where the saturation activity is A = n · σf · φ · Nf · νd , erative least-squares method [226]. The decay curves ob- and Tm the modulation period. After reaching the equi- tained with short irradiation times were used to estimate librium (n → ∞) in the short time periodic irradiation the parameters of the short-lived groups of delayed neu- (λiTm << 1) with the period of modulation Tm, the sum- trons, and the decay curves corresponding to long expo- mation of delayed-neutron counts in the counting time sure times were used to estimate the parameters of long- (Tm − tc) can be expressed as follows: lived groups. The result of the analyses of data from one run of measurements comprising n-cycles of irradiation t2 X dN(tk) and counting is the set of the relative abundances of de- · ∆t − B · (t − t ) = dt k 2 1 layed neutrons ai, the half-lives of their precursor Ti, and t1 m the covariance matrix of the group parameters [235]. −λi·t1 −λi·t2 X ai (e − e ) c. Modulated accelerator beam technique. In delayed A · (1 − e−λi·tirr ) · · = λ (1 − e−λi·Tm ) neutron measurements, a short irradiation time (see i=1 i Eq. (15)) simplifies the interpretation of the experimen- Tm − tc A · tirr (21) tal data obtained in the total delayed yield measure- Tm ments. This idea helped to develop an efficient method for the determination of the energy dependence of the where t1 is the time when the delayed-neutron count- total delayed-neutron yield from fission of several heavy ing tc starts, and t2 is the period of modulation Tm (see nuclei [217]. This technique is based on a modulated Fig. 24). Another advantage of this method as compared neutron-source generated by a Cockcroft-Walton acceler- to the nuclear reactor neutron-source techniques is a min- ator which utilizes the D(d,n) reaction for producing 3.1 imum of scattering and thermalizing material present in MeV and the T(d,n) reaction for 14.9 MeV neutrons. Fig- the vicinity of the sample. The total delayed-neutron ure 24 shows the time sequences used for the modulated yield was measured from the fission of 232Th, 233U, 235U, neutron-source technique. 238U, and 239Pu [217] nuclides. Later, Krick and Evans [218], using the same methodology, measured the total delayed-neutron yield as a function of the neutron energy from the fission of 233U, 235U, 238U, 239Pu, and 242Pu in the energy range 0.1 - 6.5 MeV. It should be noted that the on-beam arrangement of a neutron detector can lead to a degradation of its counting characteristics. Special measures should be taken to find and eliminate this ef- fect (see paragraph devoted to measurements with the T(d,n)4He neutron source). d. The neutron detector in delayed-neutron measure- ments. The neutron detector most commonly used in delayed-neutron measurements is an assembly of 3He or FIG. 24. Time sequences in the modulated neutron-source boron counters distributed in a polyethylene or paraffin techniques [217]. Delayed neutron counting is made in time moderator. A moderator allows one to shift the neu- intervals (Tm − tc). tron energy to the region of high registration efficiency. The efficiency of the neutron detector is determined using A long-counter neutron detector operating in anti- calibrated neutron sources (Pu-Be, Am-Li, 252Cf) or the synchronism with the accelerator beam registers delayed monoenergetic neutrons produced from 7Li(p,n), D(d,n), 51 neutrons in the time interval (Tm − tc). Interpretation V(p,n) nuclear reactions measured at a charged-particle of the experimental data for this method is based on accelerator [236].
40 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
e. Fission rate determination in accelerator and measured decay curves of the delayed neutron activity in reactor-based methods. The main method used for this the measurements of the relative abundances and peri- purpose in nuclear reactor experiments is spectroscopy of ods from fission of 232Th, 233U, 235U, 236U, 238U, 239Pu, gamma rays emitted from fission products. Keepin et al. 237Np, and 241Am by neutrons in the energy range from [13] determined the total number of fission in the samples 14 to 18 MeV. Details of the procedure are presented in by standard counting of the 67 h β-activity from 99Mo. [239]. Benedetti et al. [211] determined the fission rate in the samples by measuring the gamma activity of fission prod- ucts (103Ru, 131I, 140Ba/140La) that were induced in the 2. Delayed neutron integral spectra samples after irradiation. In the accelerator-based tech- nique, the fission rate in the sample is determined by mea- There are two main methods used for the investiga- suring the neutron flux with the help of a fission chamber tion of integral delayed-neutron spectra. These are neu- that is placed in the immediate vicinity of the sample. In tron spectrometry using 3He [240] and proton-recoil de- the modulated neutron-source method, the fission rate tectors [17, 18, 223] and the time-of-flight method [241]. in the fissionable sample is determined by measuring the The time-of-flight method gives an excellent resolution count rates in two fission chambers sandwiching the sam- at low energy, however this resolution worsens rapidly at ple ([217], [218]). In the case of accelerator-based meth- energies above a few hundred keV. As an example, the ods with off-beam arrangement of the neutron detector, full width at half maximum (FWHM) of a TOF system the fission chambers are placed in front of and behind the with flight path 50 cm and time resolution 3 ns is 1.6, sample along the beam line of the accelerator ([222],[228], 17.3, and 82.6 keV at energies 50, 300 and 900 keV, re- [237], [229]). spectively. The proton-recoil spectrometer resembles the f. The effect of the concomitant neutron source time-of-flight method: the resolution is good at low ener- 3 D(d,n) He in delayed- neutron experiments with an accel- gies, while deteriorating at high energies. The full width 4 erator neutron source T(d,n) He. In experiments using at half maximum of a proton-recoil spectrometer with 2.5 the T(d,n)4He reaction as the neutron source, one has atm of H2 for the same energies mentioned above, is 4.2, to take into account the effect of a concomitant neutron 12.6 and 25.6, respectively. The 3He-spectrometer seems source which originates from the implantation of deuteron to be the optimal tool over the whole energy range. The ions in the backing of the tritium target. As the number resolution of this type of spectrometer is about 12-30 keV of implanted ions grows, the intensity of the concomi- for energies ranging from 10 keV to 3 MeV. A detailed 3 tant neutron source D(d,n) He increases. This feature description of the techniques used for measurements of 4 of the T(d,n) He reaction makes it difficult to interpret delayed-neutron spectra including unfolding procedures the experimental data obtained in the respective experi- are given in a comprehensive review by Das [242]. Since ments. A special method has been developed to measure the publication of this review, experimental work on com- 3 the intensity of the D(d,n) He neutron source in correla- posite delayed-neutron spectra has been carried out only tion with the ion charge accumulated on a tritium target at IPPE using an accelerator-based neutron beam [243]. [238]. This method allows one to account for the effect By applying the least-squares method to the composite of the concomitant neutron source and has been used spectra measured in several delayed time intervals follow- in measurements of the relative abundances and periods ing fission [244], one can unfold the group energy spectra 232 of delayed neutrons for neutron-induced fission of Th, χ (E ). The main purpose of obtaining the experimen- 233 235 236 238 239 237 241 i n U, U, U, U, Pu, Np, and Am in the tal group spectra is to use them to calculate the effective energy range from 14 to 18 MeV. Details of the procedure fraction of delayed neutrons βeff [245]. are presented in ([238], [239]). g. The effect of degradation of the neutron detec- tor counting rate characteristic in an intense field of B. New measurements and compilation high energy neutrons. The delayed-neutron decay curve recorded in measurements of relative abundances and pe- riods of delayed neutrons from the fission of heavy nu- Since the last evaluation of macroscopic delayed- clides induced by neutrons from the T(d,n)4He reaction neutron parameters was performed in the framework of was found to undergo noticeable degradation in the first WPEC-SG6 [14, 246], experimental activities related to several seconds after switching off the ion beam. A special macroscopic delayed neutron data have mainly been car- experiment has been conducted to measure the counting ried out at IPPE (Obninsk). These activities have been property of the neutron detector after its irradiation by devoted to the measurement of the energy dependence neutrons from the T(d,n)4He reaction. In this experi- of the total delayed neutron yield and the relative abun- ment, the time dependence of the count rates from the dances of delayed neutrons and half-lives of their pre- cursors for neutron-induced fission of 232Th, 233U, 235U, Am-Li neutron source was measured immediately after 236 238 239 237 241 the irradiation of the neutron detector by a T(d,n)4He U, U, Pu, Np, and Am in the energy range neutron flux of different intensities in the energy range from thermal to 18 MeV [226, 228–230]. In addition, the 14-18 MeV. The obtained data were used to correct the integral energy spectra of delayed neutrons from ther- mal neutron-induced fission of 235U have been measured
41 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
[243]. The relative abundances and half-lives of delayed neutrons from the relativistic proton-induced fission of TABLE III. Results of the regression analysis of the energy dependence hT (E )i = A + B · E . 238U have also been measured [247]. All the data re- n n lated to neutron-induced fission have been compiled in Isotope Intercept value A Slope value B the macroscopic section of the database developed by this (s) (s/MeV) CRP [40] under the guidance of the Nuclear Data Section. 233U 12.63±0.06 -0.107±0.028 The database is an extension of the summary of delayed 235U 9.04±0.01 -0.086±0.005 236 neutron parameters by Spriggs and Campbell [248] and U 8.08±0.08 -0.114±0.020 238 a the compilation by Tuttle [249] and comprises the to- U 5.56±0.05 -0.070±0.012 239 tal delayed neutron yields, and the relative abundances Pu 10.54±0.04 -0.200±0.014 and half-lives of delayed neutrons in the 6- and 8-group a The energy range for analyses of 238U data was 3.2-5 MeV. models with covariance and correlation data. The main features of the new IPPE macroscopic delayed-neutron data and the earlier measured data are presented in the half-life of delayed-neutron precursors for all fissioning following sections. systems under consideration decreases with increasing en- ergy of primary neutrons in the energy range from ther- mal to 5 MeV. The real scale of the variation in the rel- 1. Energy dependence of the relative abundances and ative abundances and half-lives of delayed-neutrons can half-lives of delayed neutrons be estimated in terms of the average half-life of delayed- neutron precursors with the help of the regression anal- ysis of the IPPE data shown in Fig. 27 for the ura- The energy dependence of the relative abundances a i nium isotopes 233U, 235U, 236U, 238U and 239Pu. Re- and half-lives T of delayed neutrons have been mea- i sults of the regression analysis of the energy dependence sured at the IPPE electrostatic accelerators for neutron- hT (E )i = A + B · E are presented in TableIII. induced fission of 232Th, 233U, 235U,236U, 238U, 239Pu, n n The temporal parameters of delayed neutrons pre- 237Np, and 241Am in the energy range from thermal to sented in this paragraph have been used in the evalua- 18 MeV. Monoenergetic neutron beams were generated tion of the relative abundances of delayed neutrons and by means of the nuclear reactions T(p,n), D(d,n), T(d,n). half-lives of their precursors (see SectionIX). The method used in these experiments is described in SectionVA1. Efforts were made to improve the exper- imental and data processing procedures by a) shorten- ing the transportation time, b) extending the delayed- 2. Energy dependence of total delayed neutron yield neutron counting time, and c) averaging the temporal delayed-neutron parameters obtained in different experi- The energy dependence of the total delayed-neutron mental runs to improve the accuracy of the delayed neu- yield from neutron-induced fission of 232Th, 233U, 235U, 236 238 239 237 241 tron temporal parameters (ai,Ti). The IPPE measure- U, U, Pu, Np and Am in the energy range ments in the high energy range (14 – 18 MeV) revealed 0.3 – 5 MeV measured at the IPPE accelerator facility is two plausible reasons for the large discrepancies observed presented in Figs. 28 and 29. in the delayed-neutron temporal parameters: (i) the dete- Most of the data have been obtained after the last rec- rioration of the counting efficiency of the neutron detector ommendations were made for the major nuclides 235U, in a high-intensity neutron flux and ii) the contamination 238U and 239Pu [246]. The IPPE data are compared of the accelerator target by neutrons produced by the with both the earlier measured data and the evalu- concomitant D(d,n)3He neutron source in measurements ated data from the ENDF/B-VII.1 [24], JEFF-3.1.1 [28], using the T(d,n)4He reaction. and JENDL-4.0 [258] data libraries. The total delayed- The comparison of the delayed-neutron parameters neutron yields obtained on the basis of the systematics (ai,Ti) obtained in the IPPE measurements with other and correlation properties of delayed neutrons νd(En) = b published data is shown in Figs. 25 and 26. The quan- a · hT (En)i [227] are also presented in Figs. 28 and 29. tity that is compared is the average half-life hT i of the The IPPE data [227–230, 256, 257] are in agreement delayed-neutron precursors introduced by Piksaikin et with the corresponding data measured by Krick et al. al. [227]. Since hT i is the combination of the delayed [218] for target nuclides 233U, 235U, 238U and 239Pu. The neutron parameters (ai,Ti) obtained in the least-squares energy dependence of both sets of data for the nuclides fitting procedure and represents an average property of under consideration is well reproduced by the data anal- the system, it is not affected by the correlations which in- ysis on the basis of the correlation properties of delayed b evitably exist between the delayed-neutron parameters. neutrons νd(En) = a · hT (En)i at least in the energy It should be noted that the hT i value depends neither range up to 3-4 MeV [227]. According to the WPEC- on whether one considers N as individual precursors or SG6 recommendations [246] based mainly on the mea- delayed-neutron groups of precursors, nor on the time surements of the effective delayed neutron fraction βeff boundaries of the delayed neutron groups. at thermal and fast assemblies, the energy dependence of It can be seen from Figs. 25 and 26 that the average the total delayed-neutron yield for 235U and 239Pu does
42 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 25. The energy dependence of the average half-life of delayed-neutron precursors from neutron-induced fission of 233U, 235U, 236U, 238U. The references of the data can be taken from the compilation by Spriggs and Campbell [248] except for the IPPE data [226, 230, 235, 239, 250, 251].
not exceed 1%. The energy dependence of νd in a 1 MeV interval calculated on the basis of Krick’s data [218] shows TABLE IV. Comparison of the energy dependence of the total delayed neutron yield from fission of 233U, 235U, 239Pu. the following relative increase: 233U - 1.98% (±0.6%), 235U - 0.52 % (±1%), 239Pu - 2.06 % (±0.5%). In Ta- dνd/νd/1 MeV (%) bleIV these data are compared with the corresponding Isotope Krick et al. Keepin IPPE [227] b data by Keepin [13] and values calculated on the basis [218] [13] from νd(En) = a · hT (En)i b 233 of the systematics νd(En) = a · hT (En)i [227]. As can U 1.98±0.6 4.80 1.84±0.62 be seen in the table both the experimental data and the 235U 0.52±1.0 2.47 2.04±1.31 239 data obtained on the basis of the systematics show note- Pu 2.06±0.5 3.13 6.26±1.69 worthy energy dependence of the total delayed neutron yield. The difference between the thermal and fast energy 235 239 νd values obtained for U and Pu by the summation calculation on the basis of the JEFF-3.1.1 fission yield data [28] is 24.5% and 21.2%, respectively [260]. Thus, this strong indication of an energy dependence 233 235 239 of νd for U, U, and Pu observed in direct mea- surements is in contrast to the results obtained with help way to do this is to conduct relative measurements with of the βeff methods [246]. More work is needed to deter- accuracy ±1% [261]. Apart from this, the reason why the mine the energy dependence of the total delayed neutron summation method gives such a large difference between yields reliably. The most efficient and straightforward thermal and fast νd values also needs to be explored [260].
43 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 26. The energy dependence of the average half-life of delayed-neutron precursors from neutron-induced fission of 232Th, 237Np, 239Pu, 241Am. The references of the data can be taken from the compilation by Spriggs and Campbell [248] except for the IPPE data [226, 230, 232–234, 252–255].
3. Integral energy spectra of delayed neutrons from fission s after the end of irradiation. The second set of mea- of 235U by thermal neutrons surements was made with an irradiation time interval of 20 s and measurements of delayed neutron spectra in the A comprehensive analysis of the composite delayed time intervals: 0.12-1 s, 1-2 s, 2-3 s, 3-4 s, and 4-44 s after neutron (DN) spectra measurements can be found in the the end of irradiation. About 3000 irradiation-counting review by Das [242]. This review includes also the en- cycles have been made within these measurements. The ergy spectra from individual precursors measured by [19], time sequence of irradiation and counting intervals used [15, 16] and [17, 18], and the spectra of 235 nuclides cal- made it possible to consider some of the measured spec- culated with the help of the evaporation model. Recently tra as quasi-equilibrium. The comparison of these quasi- the integral energy spectra were measured from the ep- equilibrium spectra with the corresponding spectra from ithermal neutron induced fission of 235U at the IPPE ac- other experiments are presented in Fig. 30. celerator based neutron beams. The measurements were It can be seen from Fig. 30 that the IPPE [243] near- made to emphasize particular groups of delayed neutron equilibrium spectrum has more pronounced peak struc- precursors using the high-resolution 3He-spectrometer of ture compared to the corresponding experimental data FNS-1 type and applying different irradiation and count- obtained by other authors. The peak structure observed ing time intervals. Two sets of measurements were per- in the data of Evans et al. [240] is similar to the peak formed. The first one was made with an irradiation time structure of the IPPE data but the low-energy part of interval of 120 s, which was followed by measurements of the spectra (below 300 keV) is about 30% lower than delayed neutron spectra in a sequence of time intervals both the IPPE data and the data obtained by the sum- as follows: 0.12-2 s, 2-12 s, 12-22 s, 22-32 s, and 32-152 mation method using the microscopic spectra from indi-
44 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
that is independent of the incident neutron energy [262]. LOENIE is composed of a cylindrical polyethylene block with a central hole for the fissile material and 16 smaller holes for the 3He counters. An airtight fission chamber containing a deposit of fissile material was used as target in the experiment. With the fission chamber one esti- mates the fission rate during the irradiation phase, which is then used to normalize the delayed-neutron activity curve. During the decay phase, the 3He counters register delayed neutrons as a function of time. This time- dependent activity curve was used to derive the groups abundances in the 8-group model. Several irradiation cycles were performed, each with a different irradiation duration in order to saturate a specific group of precur- sors. FIG. 27. The energy dependence of the average half-life The result obtained for thermal neutron-induced fission of delayed-neutron precursors from neutron-induced fission of 235U from this first experimental campaign in 2018 is 233 235 236 238 239 of U, U, U, U and Pu based on the data compared with other experimental data and evaluated from IPPE [226, 230] and Keepin et al. [13]. The data of libraries in Fig. 31. More ALDEN measurements for the Maksyutenko et al. [248] presented in Figs. 25-26 were not other major actinides are planned to take place in the considered. near future. vidual precursors. This could be explained by possible errors in the determination of the efficiency of the 3He- C. Estimation of energy dependence of total spectrometer or the correction that accounts for the neu- delayed neutron yields tron flux attenuation in the shielding used for rejection of gamma rays [243]. Another important feature of the As discussed in SectionVB2, the incident-energy de- present near-equilibrium spectrum is the lower intensity pendence of the total delayed neutron yield has been mea- of DN below 200 keV as compared to the earlier exper- sured extensively by the IPPE group in the past twenty iments with the exception of the data of Evans et al. years [222, 228–230, 236]. On the other hand, the theoret- [240]. Comparison of the IPPE integral energy spectra ical study of the energy-dependence of the total delayed measured at different time intervals after the end of the neutron yield has not been pursued that actively. The irradiation with the spectra calculated on the basis of evaluated libraries, such as ENDF/B-VIII [25], JEFF- microscopic DN data is presented in SectionVIC. More 3.1.1 [28] and JENDL-4.0 [258], either perform an inter- details of the experimental procedure employed in the polation of the existing data or use the systematics pro- IPPE experiments can be found in [243]. posed by Tuttle [263], Benedetti [211], Waldo [207], and Manero [264]. As a result, the evaluated values lie on a structureless curve. 4. ALDEN: new measurements of delayed neutron data The dependence of total delayed neutron yields on the incident neutron energy results from two factors. One is Apart from the experimental activities at IPPE (Ob- the dependence of fission product yields on the incident ninsk), recently new experiments have been conducted at neutron energy. Fission yield distributions are usually the Institut Laue-Langevin (ILL) to measure aggregate determined by the potential energy surface of fissioning delayed-neutron data. The ALDEN (Average Lifetime of nuclei, which varies with incident neutron energy and ex- DElayed Neutrons) project consists of a series of experi- citation energy of the compound nuclei. The other factor ments designed to measure delayed neutron group abun- is the competition with open channels other than (n, f). dances with the aim of improving the data and provid- In particular, (n, n0f) and (n, 2nf) channels have a large ing realistic uncertainties and correlations for the most impact on the energy dependence of the total delayed important fissioning systems involved in reactor applica- neutron yields. Similarly, other channels accompanying tions (235U, 238U, 239Pu, 241Pu) [262]. multiple particle emissions like (n, 3nf) also affect the The first experimental campaign of ALDEN was held energy dependence of delayed neutron emission. at the end of 2018 at ILL (Institut Laue-Langevin) [262]. In Ref. [265], a new method for evaluating the inci- The ILL reactor produces a very intense and well- dent neutron energy dependence of total delayed neutron characterized thermal neutron flux. The experiment con- yields using least square fitting was presented. The model sisted of irradiating the fissile target and letting it decay adopted the incident neutron energy dependence of fission naturally. Delayed neutrons were measured by LOENIE, product yields, the multi-chance fission, and the evalu- a 3He detector, optimally designed to have an efficiency ated decay data of Ref. [35](Z = 2 – 28). As mentioned
45 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 28. The energy dependence of the total delayed-neutron yield from neutron induced fission of 233U, 235U, 236U, 238U. The references can be taken from the evaluation by Tuttle [249] except for the IPPE data by Piksaikin et al. [227, 228, 230, 256, 257] and Roshchenko et al. [229].
in SectionIII, new evaluated ( T1/2,Pn data for Z > 28 odd nuclei, respectively. The renormalization is carried nuclei were produced within the CRP [36] in addition out after applying the multiplication factors. The energy to Ref. [35]. New calculations of the incident neutron en- dependence of the odd-even effect is approximated as ergy dependence have thus been performed using the new CRP data of sectionIII, and the results are presented in 0 2 Ck(En) = Ck (23) the following. 1 + exp (ξEn) In the model of [265], the most probable charge Z is p where C0 (k = 1, 2) are taken from Madland and Eng- expressed by the following equations; k land’s evaluation [269]. The values of α, β, γ in Eq. (22) 0 2 and ξ in Eq. (23) are determined by a least squares fit- Zp(A, En) = αZp (A) + βc(A)En + γEn. (22) ting to experimental data on total delayed neutron yields The function c(A) accounts for the fact that the energy for 233,235,236,238U and 239−242Pu. More details about the dependence of Zp changes smoothly for light fragments method can be found in Ref. [265]. and relatively rapidly for heavy fragments [266, 267]. The For 240,241,242Pu, the parameter γ was set equal to zero 0 parameter Zp (A) is taken from the evaluation of England because there are not enough available experimental data and Rider [268]. points for all the parameters to be determined uniquely. To take into account the odd-even effect observed in The least square fitting was thus performed only for α, β, fission yield distributions, the method proposed by Mad- and ξ for the above-mentioned actinides. The parame- land and England [269] was adopted, where (1 + C1), ters obtained from fitting the new evaluated CRP data (1 − C1), (1 − C2), and (1 + C2) are multiplied by the are listed in TableV, The differences observed between fission yields of even-even, odd-odd, odd-even, and even- the values in this table and those in Table 3 of [265] re-
46 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 29. The energy dependence of the total delayed-neutron yield from neutron-induced fission of 232Th, 237Np, 239Pu, 241Am. The references can be taken from the evaluation by Tuttle [249] except for the IPPE data of Piksaikin et al. [227, 230, 256, 257, 259].
flect the impact of the new evaluated CRP data (T1/2,Pn) from SectionIII. It should be noted that in this work, the TABLE V. Parameters determined by the least square fitting to the experimental data. total delayed neutron yield of 234U (i.e. from 234U(n, f)) taken from JENDL-4.0 is multiplied by a factor of 1.2. Nuclide α β γ ξ −4 This increase improves the results obtained for the de- U-233 1.00380 0.855 −5.083 × 10 1.272 −4 235 236 U-235 0.99994 1.520 −8.631 × 10 1.124 layed neutron yields of U and U at energies where −3 235 0 U-236 1.00028 −0.195 1.596 × 10 1.924 the second-chance fission U(n, n f) or the third-chance −3 fission 236U(n, 2nf) occurs. This prescription is validated U-238 1.00411 −1.922 4.101 × 10 0.948 Pu-239 1.00013 −0.010 6.816 × 10−4 0.530 by the large ambiguity affecting the original total delayed 234 Pu-240 1.00166 0.326 0 0.379 neutron yield value of U of JENDL-4.0 as it is obtained Pu-241 0.99951 0.794 0 0.741 from Tuttle’s systematics [263]. Pu-242 1.00373 0.414 0 1.554 The total delayed neutron yields for uranium isotopes of A = 233, 235, 236, 238 are shown in Fig. 32 together with the experimental data and the evaluated values from JENDL-4.0 [258], ENDF/B-VII.1 [24] and JEFF- new results give a smooth curve up to about 4 MeV, 3.1.1 [28] libraries. The energy dependence of the exper- while the evaluated data give a straight line. Second, imental data for 233,235,236,238U is reproduced reasonably they are systematically higher than the evaluated data in well with the new parameterization. the energy region from 8 to 10 MeV. At energies above The present results agree roughly with the evaluated 15 MeV, the new results are lower than the evaluated data, however, several differences are noted: First, the data for 236,238U. At these high energies, measurements
47 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
closest to JENDL-4.0 up to around 12 MeV, but tend to be lower above 14 MeV except for 239Pu. However, there are fewer available experimental data for the plutonium isotopes 240,241,242Pu compared to the uranium isotopes. As in the case of the uranium isotopes, measurements of total delayed neutron yields at energies close to the second-chance fission threshold are required, in order to make a reliable extrapolation up to 20 MeV.
1.0 This Work JENDL-4.0 0.9 U-233 ENDF/B-VII.1,JEFF-3.2 Keepin,1957 0.8 Brunso,1955 Masters,1969 0.7 Conant,1971 Rose,1957 0.6 Notea,1969 Piksaikin,2006,2013 0.5 Krick,1970 0.4 0.3 FIG. 30. Quasi-equilibrium spectrum of delayed neutrons (neutron/100 fissions)
235 d 0.2 from the fission of U by epi-thermal neutrons [243]. IN ν 0.1 IPPE spectrum (1): counting time is 1.88 s after irradiation 0.0 This Work of 120 s. In IPPE spectrum (2): counting time is 0.88 s after 2.6 JENDL-4.0 U-235 ENDF/B-VII.1 irradiation of 20 s. Quasi-equilibrium spectra by other au- 2.4 JEFF-3.2 2.2 Piksaikin,1999 (absolute value) thors have been calculated with the help of the group spectra Krick,1970 2.0 Clifford,1972 obtained by decomposition of appropriate aggregate spectra Cox,1974 1.8 Synetos,1979 in different time windows. The data references can be found 1.6 Piksaikin,1999 in the review of Das [242]. 1.4 1.2 1.0 0.8
(neutron/100 fissions) 0.6 Masters,1969 d Besan,1977 ν 0.4 Keepin,1957 Conant,1971 0.2 Loaiza,1998 0.0 1.8 This Work 2.8 JENDL-4.0 U-236 ENDF/B-VII.1 2.6 JEFF-3.2 2.4 Bobkov,1989
) Brady,1989 2.2 Tuttle,1979 1.6 2.0 Gudkov,1989 1.8 Roshchenko,2006; Piksaikin,2013 1.6
/100fissions 1.4 1.4 1.2 ENDF/B-VII.1 Keepin, 1957 1.0 JEFF-3.1.2 Conant, 1971 0.8 (neutron/100 fissions) (neutrons d JENDL-4.0 Besant, 1977 0.6 1.2 ν Krick, Evans, 1970 Loaiza, 1998 0.4 0.2 Clifford, 1972 Piksaikin, 1999 0.0 Cox, 1974 Piksaikin, 1997 (absolute value) This Work 6.0 JENDL-4.0 Synetos, 1979 Foligno, 2019 U-238 ENDF/B-VII.1 1.0 5.5 JEFF-3.2 5.0 Keepin,1957,1969 0.0 0.5 1.0 1.5 2.0 Rose,1957 E (MeV) 4.5 Cox,1970 n Cox,1974 4.0 Besant,1977 3.5 Brunson,1955 3.0 FIG. 31. The energy dependence of the total delayed-neutron Masters,1969 235 2.5 Krick,1970 yield from neutron induced fission of U. The red cross Maksyutenko,1959 2.0 Clifford,1972 Herrmann,1966 comes from the 2018 ALDEN campaign [262]. The other ex- 1.5 Notea,1969 (neutron/100 fissions)
d Benedict,1970 perimental data are taken from the evaluation of Tuttle [249] ν 1.0 Bucko,1966 0.5 Brown,1971 except for the IPPE data which are from by Piksaikin et al. Piksaikin,2002 0.0 [256, 257]. 0 2 4 6 8 10 12 14 16 18 20 En (MeV) of total delayed neutron yields are required to improve FIG. 32. Delayed neutron yield of uranium isotopes with A = the evaluated data. 233 to 238. The data are taken from Bobkov et al. [270], Figure 33 shows the total delayed neutron yields for Loaiza et al. [216], Brady et al. [132], Gudkov et al. [214], and the evaluation by Tuttle [249, 263], while the IPPE data plutonium isotopes of A = 239, 240, 241, 242. As seen are taken from Piksaikin et al. [227, 228, 230, 256, 257] and in the uranium isotopes, the energy dependence of the Roshchenko et al. [229]. total delayed neutron yields is well reproduced. When compared with the evaluated data, the present results are
48 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
2011) [26]. For 235U, the present calculations and 1.0 This Work JENDL/FPY,FPD-2011 reproduce Keepin’s data rea- 0.9 Pu-239 JENDL-4.0 ENDF/B-VII.1 sonably. However, JENDL/FPY,FPD-2011 overesti- 0.8 JEFF-3.2 Piksaikin,2006,2013 mates the delayed neutron activity at times t < 5 sec. 0.7 Masters,1969 Keepin,1957 For 239Pu and 233U, the delayed neutron activities of 0.6 JENDL/FPY,FPD-2011 overestimate those of Keepin’s 0.5 six group evaluation [13] significantly, while the present 0.4 Conant,1970 calculations are rather close to Keepin’s results within 0.3 Besant,1977 Notea,1969
(neutron/100 fissions) errors.
d 0.2 Rose,1957
ν Krick,1970 0.1 Brunson,1955 Figure 35 shows the delayed neutron activities for the Cox,1974 235 238 233 239 240 0.0 fast neutron fission of U, U, U, Pu and Pu. 1.1 Pu-240 Similar to the results observed for thermal neutron fis- 1.0 sion, the delayed neutron activity obtained with the new 0.9 0.8 model calculations is close to Keepin’s six group evalu- 0.7 ations within errors. JENDL/FPY,FPD-2011 overesti- 0.6 mates Keepin’s data significantly except for 233U. 0.5 0.4 This Work 0.3 JENDL-4.0 0.0045 (neutron/100 fissions)
d ENDF/B-VII.1,JEFF-3.2 Keepin
ν 0.2 0.0040 Benedetti,1982 U-235 JENDL/FPY,FPD-2011 0.1 Cesana,1980 0.0035 This Work Keeipn,1957,1969 0.0 0.0030 This Work 0.0025 1.8 JENDL-4.0 Pu-241 ENDF/B-VII.1 0.0020 1.6 JEFF-3.2 Cox,1961 0.0015 Meadows,1976 0.0010 1.4 Cesana,1980 0.0005 1.2 Benedetti,1982 DN act. x time (s / fission) Keepin,1969 0.0000 1.0 0.0018 Keepin Pu-239 JENDL/FPY,FPD-2011 0.8 0.0016 This Work 0.0014 0.6 0.0012
(neutron/100 fissions) 0.0010
d 0.4
ν 0.0008 0.2 0.0006 0.0 0.0004 This Work 2.2 Pu-242 JENDL-4.0,ENDF/B-VII.1 DN act. x time (s / fission) 0.0002 2.0 JEFF-3.2 0.0000 Krick,1972 0.0018 Keepin 1.8 Bobkov,1989 JENDL/FPY,FPD-2011 0.0016 U-233 This Work 1.6 0.0014 1.4 0.0012 1.2 0.0010 1.0 0.0008 0.8 0.0006 0.6 0.0004 (neutron/100 fissions) DN act. x time (s / fission) 0.0002 d 0.4 ν 0.0000 0.2 10-3 10-2 10-1 100 101 102 103 104 0.0 0 2 4 6 8 10 12 14 16 18 20 Time after fission burst (s) En (MeV) FIG. 34. Delayed neutron activities multiplied by time t after neutron radiation for thermal neutron fission of 235U (top), FIG. 33. Delayed neutron yield of plutonium isotopes with 239Pu (middle), and 233U (bottom). Instant neutron radi- A = 239 to 242. The data are taken from Bobkov et al. [270], ation is assumed. The present calculation indicated by the Benedetti et al. [211], Cesana et al. [271, 272], and the eval- dashed lines are compared with the summation calculation of uation by Tuttle [249], while the IPPE data are taken from JENDL/FPY-2011 and FPD-2011 (JENDL/FPY,FPD-2011) Piksaikin et al. [230]. [26] (solid line) and the results of Keepin’s six group [13] (er- rors are shown by the shade).
Figure 34 shows the delayed neutron activities for To conclude, the incident neutron energy dependence thermal neutron fission of 235U, 239Pu and 233U. of the delayed neutron yields of uranium and plutonium The activities are plotted as a function of time af- isotopes was re-estimated with the model of Ref. [265] us- ter fission burst, t, and the figures compare the re- ing the new evaluated CRP data introduced in Sect.III. sults obtained with the current model with Keepin’s The present model reproduces the experimental total de- six group constants [13] and the summation cal- layed neutron yields reasonably well. Moreover, the de- culations using the fission product yields from the layed neutron activities at thermal and fast neutron fis- JENDL/FPY-2011 library and the radioactive decay data sions evaluated by Keepin’s six group model are also well from the JENDL/FPD-2011 library (JENDL/FPY,FPD- reproduced.
49 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
0.0020 Keepin A. A numerical (T1/2,Pn) data file for practical 0.0018 JENDL/FPY,FPD-2011 0.0016 U-233 This Work applications 0.0014 0.0012 0.0010 A numerical file of definitive T1/2 and Pn values has 0.0008 been created for use in summation and other calcula- 0.0006 tions for applications spanning fission reactor technolo- 0.0004 DN act. x time (s / fission) 0.0002 gies, anti-neutrino spectra and nuclear astrophysics. The 0.0000 file is based exclusively on the evaluated data produced 0.0045 Keepin U-235 JENDL/FPY,FPD-2011 by the CRP as described in Sect.III[43]. These data are 0.0040 This Work 0.0035 purely experimental without any input from systematics 0.0030 or theory. It is left to the user to decide how to treat 0.0025 0.0020 discrepant values (when compared with theory or sys- 0.0015 tematics) or missing values due to lack of experimental 0.0010
DN act. x time (s / fission) 0.0005 data. 0.0000 0.0110 Keepin Furthermore, the “non-numerical” features of the JENDL/FPY,FPD-2011 0.0100 U-238 (T ,Pn) tables, such as limits or approximate values, 0.0090 This Work 1/2 0.0080 have been replaced with purely numerical values, and 0.0070 0.0060 asymmetric limits have been symmetrized. A list of all 0.0050 0.0040 the modifications implemented in the numerical file with 0.0030 respect to the tables available online is given below: 0.0020
DN act. x time (s / fission) 0.0010 0.0000 Keepin • upper limits were symmetrized, i.e. if P1n < 0.2 0.0016 Pu-239 JENDL/FPY,FPD-2011 0.0014 This Work then the adopted value in the numerical file would 0.0012 be Pn = 0.1 ± 0.1. 0.0010 0.0008 lower limits are ’symmetrized’ as follows: P > x 0.0006 • n 0.0004 becomes Pn = y ± z so that x = y − z and 100 =
DN act. x time (s / fission) 0.0002 y +z. When the latter cannot be enforced, then the 0.0000 10-3 10-2 10-1 100 101 102 103 104 lower limit is doubled with an uncertainty of 50%, Time after fission burst (s) i.e. y = 2x and z = x, with y + z = 3x ≤ 100.
FIG. 35. Same as Fig. 34, but for fast neutron fission of 235U, • asymmetric uncertainties were symmetrized follow- 238U, 233U, and 239Pu. ing the prescription:
a − b a + b x(a) → y±z; y = x+ , z = b 2 2 However, new evaluations of fission yield data which • approximate values were assigned 50 % uncertainty, are able to reproduce both decay heat and delayed- i.e. P1n ≈ 0.2 was replaced by P1n = 0.2 ± 0.1 neutron emission data simultaneously, using the summa- tion method, are urgently needed. New measurements of • when multiple activities are measured but cannot total delayed neutron yields at different incident neutron be distinguished, i.e. when a Pn is given for both energies are also important for the evaluation of both fis- the ground state and isomeric states, then the latter sion yields as well as delayed neutron yields. is split equally among the states while the uncer- tainty for each state remains the same as the orig- inal recommended one. As an example, one takes for 80Ga a recommended value for both, the ground state and the isomeric state , of P1n = 0.90(7)% . In the numerical file, the value P1n = 0.45(7)% was adopted for the ground state and isomeric state, re- VI. SUMMATION CALCULATIONS spectively. If the original uncertainty is larger than the split P1n value, one can adopt a 100% uncer- Summation and time-dependent calculations of total tainty for the split value. For example, in the case 116 delayed neutron yields, delayed neutron spectra and de- of Rh, where P1n < 2.1% for both ground state layed neutron activities using the CRP evaluated micro- and isomeric state, the limit is symmetrized to give scopic (T1/2,Pn) data are compared with recommended P1n = 1.05 ± 1.05 % and the resulting Pn is split data and evaluated libraries in this section. The aim is among the two states with 100% uncertainty, i.e. to assess the new CRP data against a broad range of as P1n = 0.525 ± 0.525 % for the ground state and available macroscopic data. isomeric state, respectively.
50 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
• there are a few cases (38,39Al, 43P) where the re- for equilibrium irradiation where the production and de- ported measured Pn value can take unphysical val- cay rates of the precursor nuclides are proportional to ues within the reported uncertainty range. In such their cumulative fission yields. cases the uncertainty was adjusted so that Pn lies Assuming the Nn and CFY quantities in the above within the physical limits (0 ≤ Pn ≤ 100 %). expression are independent, the uncertainty of the total DN yield is given by • zero values for Qβxn and Pxn have the following meaning: (i) there are no measurements available, X X (ii) the decay mode is forbidden. σ2 = N 2 · σ2 + σ2 · CFY 2 (25) νd ni CFYi Nni i i i • nuclides that have only been identified in measure-
ments and have no measured T1/2 or Pn value are where σNni and σCFYi are the uncertainties of the de- not included in the numerical file. layed neutron multiplicities Nni and cumulative fission yields CFYi of the individual precursors. The format of the file is as follows: In the following calculations, the cumulative fission Nuclide name, Z, A, liso (index indicating order of states, yields CFYi are taken from the JEFF-3.1.1 fission yield li- i.e. for ground state liso = 0, for 1st isomeric state liso brary [28] which has been widely tested in a broad range = 1, etc.), Level Energy, ∆ Level Energy, %β−, ∆%β−, of applications, including in summation calculations of Qβn1, ∆Qβ1n, Qβn1, ∆Qβ2n, Qβn3, ∆Qβ3n, T1/2, ∆T1/2, DN yields, and has proved to perform well. The β−decay P1n, ∆P1n, P2n, ∆P2n, P3n, ∆P3n. and β-delayed neutron data, (T1/2,Pn), are taken from The numerical file has been used in the summation cal- the new CRP tables. Nuclides for which there are no culations described in the following sections with the aim Pn values in the new CRP tables are not considered as of verifying and validating the new evaluated Pn tables β-delayed neutron emitters, i.e. Pn = 0. The remain- produced by the CRP. It is available for downloading at ing decay data (γ decays-IT), when used in the basic the IAEA Beta-delayed neutron database [43]. A more summation calculations or time-dependent calculations, detailed description of the contents of the IAEA Beta- are taken from ENDF/B-VIII.0 [25]. This combination delayed neutron database is given in SectionX. of input data is hereafter referred to as CRP+ENDF/B- VIII.0. The calculations are compared with total DN yields obtained by combining CFYi from JEFF-3.1.1 and B. Total delayed neutron yields decay data from the ENDF/B-VIII.0 decay data sub- library. This comparison enables us to observe the di- 1. Basic summation calculations rect impact of the new CRP data (T1/2,Pn) on the to- tal DN yields as all the other data inputs are the same. In this section, summation calculations of the total It is also interesting to explore the differences between delayed neutron yield ν performed with the CRP data the new improved CRP+ENDF/B-VIII.0 decay data sub- d library and the other available evaluated decay-data sub- (T1/2,Pn) are compared with other decay data libraries and with recommended values. libraries. Therefore, additional calculations were per- a. Comparison with evaluated libraries. The first formed by combining the CFYi from JEFF-3.1.1 [28] with step in verifying and testing the new evaluated CRP the decay data from JEFF-3.1.1 [28] and JENDL/DDF- 2015 [27] sub-libraries. Note that the release of JEFF-3.3 (T1/2,Pn) values using the numerical file was to perform a basic inter-comparison of summation calculations of the library occurred after most of the work carried out by the delayed neutron yields among the CRP participants. To- CRP was completed. Since JEFF-3.3 and JEFF-3.1.1 do tal delayed neutron yields (DN) were calculated for ther- not differ with respect to the DN data, it was decided mal and fast neutron-induced fission of major and mi- that for the purposes of the CRP work it was sufficient to nor actinides by four independent groups, namely CEA- present the results that had already been obtained with Cadarache (Foligno), CIEMAT (Cano), JAEA (Minato) JEFF-3.1.1 as they would in practice be identical with and SUBATECH-Nantes (Fallot). those obtained with JEFF-3.3. The results of the summation calculations for the major The calculations are based on the summation method 235 238 239,241 whereby the total delayed neutron yield ν is defined as actinides U, U, Pu are given in TableVI. Per- d fect agreement is observed among the four different sum- mation calculations which confirms that the calculations X performed by the different groups within this project are νd = Nni · CFYi (24) i equivalent and hence reliable. We do not need to report the results of all four groups hereafter, but just the re- where Nni is the average delayed neutron multiplic- sults of the CRP. The values in TableVI are considered as ity emitted by the precursor i and is given by Nn = reference values for the given combination of decay data 1 · P1n + 2 · P2n + 3 · P3n..., Pxn are the delayed neutron and fission product yields library. emission probabilities Pxn and CFYi are the cumulative The results displayed in TableVI show that overall fission yields for precursor i. The above equation holds the CRP Pn data lead to larger values of the total DN
51 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. CEA0.01799 (5%) 0.01799 (5%)0.04492 (3%) 0.01799 (5%) 0.04492 (3%) 0.017990.00722 CIEMAT (5%) (7%) 0.04492 (3%) 0.00722 (7%) 0.044920.01407 (3%) (5%) 0.00722 (7%) 0.01407 JAEA (5%) 0.00722 (7%) 0.01407 (5%) 0.01407 (5%) NANTES CEA0.01883 (5%) 0.01883 (5%)0.04604 (3%) 0.01883 (5%) 0.04604 (3%) 0.018830.00756 CIEMAT (5%) (7%) 0.04604 (3%) 0.00756 (7%) 0.046040.01462(5%) (3%) 0.00756 (7%) 0.01462(5%) JAEA 0.00756 (7%) 0.01462(5%) 0.01462(5%) NANTES CEA0.01698 (5%) 0.01698 (5%)0.04037 (3%) 0.01698 (5%) 0.04037 (3%) 0.016980.00675 CIEMAT (5%) (7%) 0.04037 (3%) 0.00675 (7%) 0.040370.0129 (3%) (5%) 0.00675 (7%) JAEA 0.00675 0.0129 (7%) (5%) 0.0129 (5%) 0.0129(5%) NANTES CEA0.01849 (6%) 0.01849 (6%)0.04511 (4%) 0.01849 (6%) 0.04511 (4%)0.00739 CIEMAT (7%) 0.04511 0.01849 (4%) 0.00739 (6%) (7%)0.01436 (5%) 0.00739 0.04511 (7%)3 0.01436 (4%) JAEA (5%) 0.00739 (7%) 0.01436 (5%) 0.01436 (5%) NANTES JEFF-3.1.1 decay data ENDF/B-VIII.0 decay data JENDL/DDF-2015 decay data CRP+ENDF/B-VIII.0 decay data thermalthermalthermal fast fast fast thermal fast ) for thermal and fast neutron-induced fission of major actinides obtained from the summation method using three different d ν CEA CIEMAT JAEA NANTES CEA CIEMAT JAEA NANTES CEA CIEMAT JAEA NANTES CEA CIEMAT JAEA NANTES U Pu 0.00680 (7%)Pu 0.00680 (7%) 0.01397 (4%) 0.00680 (7%) 0.01397 (4%) 0.00680 (7%) 0.01397 (4%) 0.01397 (4%) U Pu 0.00605 (7%)Pu 0.00605 (7%) 0.01232 (5%) 0.00605 (7%) 0.01232 (5%) 0.00605 (7%) 0.01232 (5%) 0.01232 (5%) U 0.01673 (5%) 0.01673 (5%) 0.01673 (5%) 0.01673 (5%) U 0.01477 (5%) 0.01477 (5%) 0.01477 (5%) 0.01477 (5%) U Pu 0.0066 (7%)Pu 0.01374 0.0066 (5%) (6%) 0.01374 (5%) 0.0066(7%) 0.01374 (5%) 0.01374 0.0066 (5%) (7%) U 0.01609 (5%) 0.01609 (5.%) 0.01609 (5%) 0.01609 (5%) U Pu 0.00666 (7%)Pu 0.00666 (7%) 0.01378 (5%) 0.00666 (7%) 0.01378 (5%) 0.00666 (7%) 0.01378 (5%) 0.01378 (5%) U 0.01637 (6%) 0.01637 (6%) 0.01637 (6%) 0.01637 (6%) d d d d 238 239 241 238 239 241 ν 235 ν 235 238 239 241 ν 235 238 239 241 ν 235 TABLE VI. Total delayed neutron yields ( decay data libraries (CRP+ENDF/B-VIII.0, ENDF/B-VIII.0,Four groups JEFF-3.1.1, JENDL/DDF-2015) participated in in combinationrelative the with uncertainties the inter-comparison, neglecting cumulative namely fission correlations. CEA yields CFY (France), from CIEMAT JEFF-3.1.1. (Spain), JAEA (Japan), and NANTES (France). The numbers in brackets are the
52 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
yields (νd) compared to the other evaluated decay data adopted the value of Pn = 1.311(5) from ENSDF. On libraries ENDF/B-VIII.0, JEFF-3.1.1 and JENDL/DDF- the other hand, this precursor does not have a delayed- 2015. The observed differences are larger with respect neutron branch in the JEFF-3.1.1 decay data library. to JEFF-3.1.1, which generally gives the lowest values The contribution of 89Br to the CRP list is lower by of total DN yields among all three evaluated libraries. more than 13% in comparison to JEFF-3.1.1. This differ- In the case of fast fission of 238U and thermal fission of ence is observed across several of the major actinides, and 241 Pu, the CRP+ENDF/B-VIII.0 total DN yields are in is related to the small difference in Pn values, with Pn = agreement with the ENDF/B-VIII.0 and JENDL/DDF- 13.7(6) being the CRP value and Pn = 14.1(4) [30] the 2015 results. value in JEFF-3.1.1. ENDF/B-VIII.0 and JENDL/DDF- The 20 most important contributors to the total DN 2015 contributions are identical which is due to the fact yields (νd) displayed in TableVI are listed in TableVII. that they have adopted the same value Pn = 13.8(4) from The first column includes the percentage contribution ENSDF. of the precursors resulting from the CRP (T1/2,Pn) The above-mentioned differences in the individual pre- data, while the neighboring columns include the con- cursor Pn values may have a non-negligible impact on the tributions of the same precursors obtained when us- DN group spectra which are estimated from summations ing the ENDF/B-VIII.0, JEFF-3.1.1 and JENDL/DDF- over the energy spectra of individual precursors weighted 2015 decay data, respectively. Both ENDF/B-VIII.0 by their relative contributions to the total DN yield, in and JENDL/DDF-2015 have adopted Pn values from addition to the differences observed in the total delayed ENSDF [23], however, depending on the date they were neutron yields. retrieved from the ENSDF database, the values may dif- The uncertainties associated with the new total DN fer among each other and with respect to the current yields in TableVI are fairly small compared to the ENSDF data value. JEFF-3.1.1 uses Pn values from an large uncertainties affecting the cumulative fission yields evaluation by Nichols [30] as well as from ENSDF. CFYs. However, this is expected considering that the From TableVII we see that overall, the differences be- estimation of the uncertainties ignores the effect of cor- tween the contributions of the top 20 precursors among relations which are particularly strong between the CFY the various libraries are less than 10%. Some notable uncertainties. A more rigorous treatment of the propa- cases are discussed below: gation of uncertainties including correlations would yield 137 I is the main contributor to νd for nearly all the different results. major actinides and incident neutron energies. The 137I b. Comparison with recommended values. The total contribution is enhanced (on average by a few %) when DN yields presented in TableVI are compared with rec- using the CRP data and is mainly responsible for the ommended DN yields from [246] and the recommended small increase in the total DN yields observed in TableVI. DN yield data from JEFF-3.1.1, ENDF/B-VIII.0 and The increased relative contribution of 137I is related to a JENDL-4.0 [258] libraries in Figs. 36-38. It should be small increase in the CRP Pn value for this precursor, noted that the total DN yields from JEFF-3.1.1 are iden- which is now Pn = 7.63(14)%, compared to the value tical to the JEFF-3.3 data therefore the results in Fig. 36 of 7.14(23)% used in ENDF/B-VIII.0 and JENDL/DDF- are valid for both libraries. As shown in the figures, the 2015. The increase is more pronounced with respect to total DN yields obtained with the CRP (T1/2,Pn) data for 235 239 JEFF-3.1.1 which has an adopted value of Pn = 6.5(4)% thermal fission of U and Pu are in good agreement from [30]. with the recommended data [246] and the data from the The precursor 85As is among the top twelve contribu- other libraries except for the thermal value of 239Pu from 238 tors in the CRP list with Pn = 62.5(10)%. It also fig- JENDL-4.0. The results for the fast fission of U are ures high in the ENDF/B-VIII.0 and JENDL/DDF-2015 substantially closer to the recommended data compared contributors lists due to the large value Pn = 59.4(24)% with the previous summation calculations of [187]. The adopted in both files, from ENSDF. The difference be- differences between the calculated and recommended val- tween the CRP contribution and JEFF-3.1.1 is consider- ues of the total DN yields for the fast neutron-induced ably larger (over 100%) due to the relatively lower value fission of 239Pu and 235U and for both fast and ther- 241 of Pn = 22(3)% adopted in JEFF-3.1.1 [30]. mal fission of Pu are more than 10%. These differ- 91Br is among the top 20 contributors in the CRP list ences highlight the need for better treatment of uncer- with Pn = 29.8(8)%, which is similar to the JEFF-3.1.1 tainties and further improvement of the microscopic data value Pn = 20(2)% taken from [30] and JENDL/DDF- for these nuclides, mainly the fission yield data. Another 2015 Pn = 20(3) taken from ENSDF. In ENDF/B-VIII.0, important feature of the CRP results presented in Ta- however, this precursor does not have a delayed-neutron bleVI is the clear evidence of the energy dependence of decay branch, therefore it does not contribute to the de- the total DN yields for fission of 235U and 239Pu. This is layed neutron activity. in agreement with the latest experimental data presented 136Te appears in the top 20 list of the CRP calculations in this paper (see SectionVC) and supports the need for although it has a relatively small β-delayed neutron frac- further investigation of this dependency. tion of Pn = 1.37(5)%. Similar results are obtained with ENDF/B-VIII.0 and JENDL/DDF-2015 which have both
53 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
FIG. 36. Comparison of the total DN yields obtained from FIG. 38. Same as in Fig. 36 but for the JENDL-4.0 decay the summation method using the CRP (T1/2,Pn) data and sub-library (bar chart). the JEFF-3.1.1 decay-data sub-library, respectively, with the recommended DN yields [246]. The total DN yields for the fission of various nuclides (listed on the abscissa axis) are plot- the new CRP data. ted with symbols. The right ordinate shows the percentage a. Time-dependent calculations during the irradiation difference between the total DN yield data calculated using and decay phase. The integral DN activity experiment the CRP data and the corresponding data from JEFF-3.1.1 is simulated by summing up the microscopic contribution (bar chart). of each of the precursors produced during the fission pro- cess. The Bateman solver [262] is a computer program that uses the decay data to reconstruct the family tree of the precursors and then solves the Bateman equations to obtain the precursors’ concentration during both the irradiation and the decay phase. The concentrations-in- time are then multiplied by the effective delayed-neutron emission probability of the precursor (Eq. 26).
N X nd(t) = Ci(t) P n eff i, (26) i=1
4 P where P n eff i = x Pxn is the effective delayed- x=1 neutron emission probability of the precursor i, which takes into account multiple decay branches leading to the neutron emission. The following decay-processes are con- − − − − − sidered in the code [262]: IT , β , βn , β2n, β3n and β4n. The Bateman solver has been validated through compar- FIG. 37. Same as in Fig. 36 but for the ENDF/B-VIII.0 decay sub-library. isons with the DARWIN code [273]. Figure 39 shows the DN activity for thermal neutron- induced fission of 235U computed by combining the 2. Time-dependent calculations CRP+ENDF/B-VIII.0 and ENDF/B-VIII.0 decay data with JEFF-3.1.1 fission yields. The figure shows the two Apart from summation calculations that use cumula- DN activities as w ell as their ratio, as a function of time tive fission yields to account for the abundances of the after the irradiation. fission products in equilibrium, time-dependent calcula- The ratio does not start from 1 because the νd is tions that start from independent fission yields and follow different (1.60 neutrons/100 fissions for ENDF/B-VIII.0 the evolution of the system into the decay phase, can also and 1.66 neutrons/100 fissions for CRP+ENDF/B-VIII.0 235 shed light on the impact of the new CRP (T1/2,Pn) data for U). TableVIII displays the 20 most important on integral quantities. In this section, we describe two precursors according to the CRP calculations (see also different time-dependent calculations of the integral DN TableVII), together with their partial yield νd,i, the activity and total DN yields that were performed with difference with respect to the partial yields calculated
54 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al. Pu using 241 Pu and 239 Pu fast U, 241 238 U, 235 Y 5.9 6.0 6.6 5.9 I 21.9I 21.3 6.3 21.2 6.9 20.9 I 7.1 3.9 6.7 Nb 4.2 2.4Sb 4.4 2.5 2.3 4.2 CsTe 2.7 2.6 1.6 1.4Nb 2.5 1.8 2.1 1.2 1.4 2.6 1.9 1.3 0.0 1.7 1.4 1.4 1.2 Rb 10.4Br 10.9 6.2Rb 11.5Br 4.6 6.4 10.7 4.5As 4.8 7.2Y 4.6 3.4Br 6.3 5.1 3.3 3.4 4.9 3.1Rb 4.7 2.9 4.5 2.4 3.2 1.4Br 3.2 3.3 2.2 3.5 2.0 3.2 Rb 3.1 2.4 0.0 1.4Br 2.1 1.5 1.4 1.2 1.4 1.5 1.3 1.3 1.4 1.3 m Nuclide CRP137 ENDF/B JEFF JENDL 94 138 89 98 95 90 139 85 99 88 105 93 135 91 143 136 96 106 87 U fast 238 Y 2.0 2.0 2.2 2.0 I 9.3 8.9ISb 9.0 5.4 5.0 8.9 I 5.7 5.7 4.6Cs 6.2 4.5 5.0 2.4 5.8 5.6 5.3 2.7 4.9 Sb 2.9I 1.7I 2.7 1.7 1.1 1.6 2.1 1.3 1.6 2.2 1.1 1.8 2.0 1.6 RbBr 7.6Br 6.7As 7.9 6.3 6.7 5.5 8.4 6.5Br 5.4 7.3 7.8 7.4 4.8Rb 6.7 2.2Br 6.4 4.2 0.0 5.3 2.8Rb 4.3 3.7Y 2.2 2.8 3.3 4.7Rb 2.0 2.1 3.2 4.2 1.9 1.8 2.8 2.4 2.0Rb 2.0 2.1 2.1 1.6 2.0 1.9 1.4 1.6 1.4 m Nuclide CRP137 ENDF/B JEFF JENDL 94 90 89 85 139 135 91 138 95 88 145 96 99 98 97 136 140 141 93 Pu thermal 241 Y 5.2 5.3 6.0 5.3 II 24.4I 8.4 23.2 9.0 23.6 5.8Sb 23.1 6.0 9.6 3.9 9.0 6.6 4.4 6.2 3.5Te 4.4 2.3CsNb 2.3 2.0 1.6Te 2.1 0.0 1.6 1.2Cs 2.3 2.4 1.8 1.3 1.0 2.1 1.6 1.1 1.4 1.3 1.2 1.1 Rb 8.0Br 8.2 4.9Y 8.8Br 5.0Rb 3.4 8.2 3.3Br 3.2 5.7As 3.0 3.3 3.1 3.3 5.0 2.5Rb 3.3 3.1 3.6 3.6 2.2 2.4 3.3 3.3 3.5Br 3.2 1.9 1.0 3.1 1.2Br 2.4 2.2 0.0 1.2 1.9 1.3 0.9 0.8 1.4 1.3 m Nuclide CRP137 ENDF/B JEFF JENDL 138 94 139 98 89 135 99 90 95 88 85 136 93 143 105 91 137 87 145 U fast 235 Y 5.9 6.1 6.5 5.9 I 12.9 12.7 12.2 12.3 I 3.9 4.3ICs 4.3 2.2 1.4 4.2 Sb 2.4 1.5Te 0.6 2.5 0.6 1.6 0.7 2.4 1.5 0.6 0.6 0.0 0.7 0.6 Rb 13.5BrBr 14.2 11.5Br 14.5 8.3 12.1 8.2As 13.9 8.5 13.1Rb 8.4 5.0 11.8 5.0 8.8Rb 5.0 9.0Br 5.1 8.3 3.6Y 8.1 2.0 3.1Br 3.2 5.4 3.0 4.8 3.3 2.5 5.0 3.5Rb 2.7 0.0 3.4As 3.1 1.3 2.9 3.2 0.7 1.8 1.3 3.0 Rb 1.7 0.8 1.4 0.5 0.8 1.3 0.6 0.7 0.6 0.5 m Nuclide CRP94 ENDF/B JEFF JENDL 137 89 90 88 98 85 95 138 93 87 99 91 139 143 96 86 135 136 97 for thermal and fast neutron-induced fission of the major actinides Pu fast d ν 239 Y 10.8 11.2 12.1 10.9 I 17.7 17.4 16.9 17.0 INb 2.7I 1.3 3.0Cs 1.1 1.3 3.1 0.7Te 1.1 1.4 2.9 Sb 0.8 0.6 1.3 1.2 0.5 0.6 0.8 1.1 0.6 0.7 0.0 0.4 0.6 0.5 Rb 13.6BrBr 14.4 9.2As 14.8 7.0Br 9.7 5.2Y 14.1 7.1 5.1Rb 10.6 5.2Rb 4.6 4.6 5.2 7.8Br 9.5 4.1 2.1 4.2 4.2 7.0 5.5 2.8Br 5.1 4.3 4.5 5.1 4.5 3.0 1.6 4.5 4.6 Rb 4.1 0.0 3.1 4.2 0.9As 2.9 1.2 0.9 0.6 1.1 1.0 0.7 0.9 0.7 0.6 m Nuclide CRP137 ENDF/B JEFF JENDL 94 98 89 88 85 90 99 93 95 87 138 91 105 139 96 143 86 136 135 U thermal 235 Y 4.1 4.2 4.6 4.1 II 4.7 5.1 3.5 3.7 5.3 5.0 Sb 4.0Te 2.1 3.8 Cs 1.6 2.4 1.6 1.6Te 1.7 1.9 0.0 0.8 2.4 1.8 1.6 0.9 1.6 1.0 0.9 I 16.3 15.8 15.7 15.6 BrRb 11.1Br 9.3Br 11.6 7.4As 9.8 13.0 7.3 7.6 5.3 11.4 10.2 7.4 5.3 8.1Rb 9.6 8.2Rb 7.5 3.5 2.1Br 7.3 3.4Br 5.2 3.6 3.2Y 3.1 2.7 3.8 3.5 2.2 3.4 0.0 3.5 3.6As 2.0 3.0 2.1Rb 3.4 0.9 2.2 1.9 0.9 1.0 2.2 0.8 1.0 0.9 0.9 0.8 m 89 94 90 88 85 138 98 139 95 93 87 91 99 135 136 143 86 96 137 137 Nuclide CRP ENDF/B JEFF JENDL Pu thermal 239 Y 9.4 9.7 10.6 9.5 Cs 1.4Nb 1.2 1.5Te 1.2 1.6 0.5 1.5 1.3 0.6 1.2 0.6 0.5 I 5.3 5.7I 5.9 5.6 2.8Sb 3.0 1.5Te 3.2 1.7 1.4 3.1 1.4 1.4 1.7 0.0 1.3 I 25.8 25.0 24.7 24.6 Rb 10.6Br 11.1 6.2 11.6 6.5 11.0 Br 7.2 3.6 6.4 3.7Br 3.7 2.5 3.6 2.6Br 2.8 1.1 2.6 0.0 0.9 0.8 RbY 4.0Rb 3.9 3.6 3.4As 3.4 3.9 3.5 2.8 3.7 3.5 3.7 2.8 3.8 3.4 1.1Rb 2.7 0.8 0.7 0.8 0.7 Br 5.0 5.0 5.6 5.0 m 94 89 90 87 143 105 91 137 98 138 93 99 95 139 85 135 136 96 137 Nuclide CRP ENDF/B JEFF JENDL 88 TABLE VII. The twenty most important contributors to the CRP+ENDF/B-VIII.0 decayENDF/B-VIII.0, data JEFF-3.1.1 and (CRP). JENDL/DDF-2015 In decay-data libraries, columns respectively. ’ENDF/B’, All the ’JEFF’ contributions and are given ’JENDL’, in we %. give the contributions of these precursors obtained when using
55 Development of a Reference Database ... NUCLEAR DATA SHEETS P. Dimitriou et al.
TABLE VIII. Detailed comparison of the microscopic data of the 20 most important precursors for thermal neutron-induced 235 1 fission of U: atomic number (Z), mass (A), isomeric num- ber (index) (I), the partial delayed neutron yield computed 1 by using CRP+ENDF/B-VIII.0 decay data (νd,i), the differ- ence in the partial yield when using CRP+ENDF/B-VIII.0 1 rather than ENDF/B-VIII.0 (∆νd,i) and the change in mi- croscopic Pn from ENDF/B-VIII.0 to CRP+ENDF/B-VIII.0