Nuclear Physics

Large-Scale Calculation of Fission Barrier Parameters for 5254 Nuclei with 171 ≤ A ≤ 330

Peter Möller, Arnold J. Sierk, T-16

n previous Theoretical Division Nuclear the calculated potential- surfaces. At Weapons Highlights issues, we have this stage we have extracted the heights of discussed calculations of fission the first and second barrier peaks and the potential-energy surfaces as functions fission-isomeric state for all 5254 nuclei. Iof up to five different nuclear shape coordinates as the system evolves from a Figure 1 shows a comparison between ground-state shape to two separated fission calculated and experimental barrier fragments. The five shape coordinates we parameters for a sequence of consider beyond the second minimum in the nuclei. The experimental data are from a fission barrier (fission isomer) are elongation, review by Madland [1]. Sometimes the neck radius, spheroidal deformations of the results of several experiments are plotted for emerging left and right fragments, and left- the same . We have made such right asymmetry. For less deformed comparisons for isotope chains of all shapes between the spherical shape and the elements from Th to Es (Z = 90 to Z = 99). second minimum in the barrier we consider Such detailed knowledge about the barrier three shape coordinates in the Nilsson structure and other ε ε perturbed-spheroid parameterization: 2 properties are necessary to model (n,f), ε γ (quadrupole), 4 (hexadecapole), and (axial (n,2n), and many other cross sections. The asymmetry). Important parameters of the potential-energy calculations took 30000 CPU calculated nuclear potential-energy surface, hours and the subsequent analysis 20000 such as the local minima, saddle points CPU hours (so far) on the T-16 cluster between all pairs of minima, valleys as “nuclei.” functions of elongation, and the height of the ridges that separate the valleys are extracted To model many astrophysical scenarios; for by similar immersion techniques to those example to model the end of the rapid- used in geography. Some illustrative results -capture process (r-process) in which for a few nuclei have been presented in many heavy elements are formed in stars, it Nuclear Weapons Highlights reports in is necessary to know fission-barrier heights previous years. of a large number of nuclei. When a neutron is captured in the r-process it is energetically This year we have used our models to possible for the nucleus to fission if the undertake production calculations so that neutron of the compound we obtain detailed descriptions of the fission system is larger than the fission-barrier potential-energy for a large number of nuclei. height. We display in Fig. 2 the difference Specifically we calculated 5-D potential between the calculated fission-barrier height energy surfaces for 5254 nuclei between the and the one-neutron separation energy. and neutron drip lines from A = 171 to When this quantity is negative it is A = 330. For each of these nuclei the energy energetically possible for the nucleus to is calculated for 5,009,325 different . shapes. In addition we studied the effect of triaxiality at the first peak in the barrier by The duration of the large neutron flux that calculating 3-D potential-energy surfaces for powers the r-process is thought to be of the these same nuclei in the ε parameterization. order of 1 s. Thermonuclear explosions also We have developed highly automated scripts generate large fluxes of , but these that use our immersion codes to extract are of shorter duration, and relevant barrier structure parameters from

Assocate Drectorate for Theory, Simulaton, and Computaton (ADTSC) in this environment is therefore referred to as and the threshold saddle-point for prompt neutron capture. However, this these modes. process can be thought of as an approximation to the r-process on earth [2], For more information contact Peter Möller at least over a limited range of nuclei. Several at [email protected]. nuclear weapons tests between 1952 and 1969 involved studies of the prompt capture [1] D.G. Madland, Los Alamos National process. In Fig. 2, we have indicated where Laboratory, personal communication, 2000. capture chains starting with 232Th, 238U, and [2] S.A. Becker, Carnegie Observatories Astrophysics 242Pu terminate [2]. We note that the U and Series, Vol. 4, A. McWilliam and M. Rauch, Eds. Pu chains end very near the dark green area, (Pasadena: Carnegie Observatories, http://www. ociw.edu/ociw/symposia/series/symposium4/ where nuclei would fission immediately proceedings.html). upon neutron capture. Thus, our results are very consistent with the experimental Funding Acknowledgements observations, and a very encouraging sign NNSA’s Advanced Simulation and that our models produce reliable fission- Computing (ASC), Materials and Physics Fig. 1. barrier parameters far from stability. Program. Calculated height Consequently, another exciting application is of the first E and to use our calculated database in modeling A second EB peak in the the end of the r-process, which runs far from fission-barrier and stability in regions of experimentally calculated energy EII of inaccessible nuclei where the neutron the fission isomer for separation energy is 1–2 MeV or so. The Th U , compared does not reach the dark green to experimental data area. However, it was observed [2] that this where available. The capture chain is severely blocked at calculated energies are all relative to the A = 242 and A = 244 where the neutron calculated energy of the capture cross section becomes very small. ground-state minimum. A more detailed interpretation of the observations on these capture chains requires more elaborate calculations of the reaction paths using realistic neutron spectra, which are not all thermal.

In the near future we plan to configure our immersion codes to establish when the calculated potential-energy surfaces exhibit more than one fission mode, the character (mass symmetric division or mass Fig. 2. asymmetric division) of the different modes, Difference between calculated fission- barrier height and neu- tron separation energy. When this quantity becomes negative (dark green area) fission is energetically possible. We indicate in the fig- ure the beginning and ending of these capture chains observed in weapons tests [2]. Two of the chains (U and Pu) end next to the dark green area where fission becomes ener- getically possible.

Nuclear Weapons Highlights 2007 LALP-07-041