First Gogny-Hartree-Fock-Bogoliubov nuclear mass model

Scission Deformation of the ^{120}Cd/^{132}Sn Neutronless Fragmentation in ^{252}Cf(sf)

We report on a study of the radiative decay of fission fragments populated via neutronless fission of ^{252}Cf(sf). Applying the double-energy method a perfect mass identification is achieved for these rare events. In the specific case of the ^{120}Cd/^{132}Sn cold fragmentation, we find that ^{132}Sn is produced in its ground state. We can therefore directly measure the excitation energy of the complementary fragment, ^{120}Cd. The reproduction of the γ-ray spectrum, measured in coincidence with the neutronless fission events, is sensitive to the angular momentum distribution of the studied primary fragment. The latter estimated using a time-dependent collective Hamiltonian model, allows us to constrain for the first time the deformation (β_{2}≃0.4) of the studied fission fragment at scission. The present work demonstrates the high potential of the understudied neutronless fission channel for extracting detailed information on both fission fragments and process.

Deformation versus Sphericity in the Ground States of the Lightest Gold Isotopes

The changes in mean-squared charge radii of neutron-deficient gold nuclei have been determined using the in-source, resonance-ionization laser spectroscopy technique, at the ISOLDE facility (CERN). From these new data, nuclear deformations are inferred, revealing a competition between deformed and spherical configurations. The isotopes ^{180,181,182}Au are observed to possess well-deformed ground states and, when moving to lighter masses, a sudden transition to near-spherical shapes is seen in the extremely neutron-deficient nuclides, ^{176,177,179}Au. A case of shape coexistence and shape staggering is identified in ^{178}Au which has a ground and isomeric state with different deformations. These new data reveal a pattern in ground-state deformation unique to the gold isotopes, whereby, when moving from the heavy to light masses, a plateau of well-deformed isotopes exists around the neutron midshell, flanked by near-spherical shapes in the heavier and lighter isotopes-a trend hitherto unseen elsewhere in the nuclear chart. The experimental charge radii are compared to those from Hartree-Fock-Bogoliubov calculations using the D1M Gogny interaction and configuration mixing between states of different deformation. The calculations are constrained by the known spins, parities, and magnetic moments of the ground states in gold nuclei and show a good agreement with the experimental results.

Expanding the limits of nuclear stability at finite temperature

Properties of nuclei in hot stellar environments such as supernovae or neutron star mergers are largely unexplored. Since it is poorly understood how many protons and neutrons can be bound together in hot nuclei, we investigate the limits of nuclear existence (drip lines) at finite temperature. Here, we present mapping of nuclear drip lines at temperatures up to around 20 billion kelvins using the relativistic energy density functional theory (REDF), including treatment of thermal scattering of nucleons in the continuum. With extensive computational effort, the drip lines are determined using several REDFs with different underlying interactions, demonstrating considerable alterations of the neutron drip line with temperature increase, especially near the magic numbers. At temperatures T ≲ 12 billion kelvins, the interplay between the properties of nuclear effective interaction, pairing, and temperature effects determines the nuclear binding. At higher temperatures, we find a surprizing result that the total number of bound nuclei increases with temperature due to thermal shell quenching. Our findings provide insight into nuclear landscape for hot nuclei, revealing that the nuclear drip lines should be viewed as limits that change dynamically with temperature.

Evidence of Hexadecapole Deformation in Uranium-238 at the Relativistic Heavy Ion Collider

State-of-the-art hydrodynamic simulations of the quark-gluon plasma are unable to reproduce the elliptic flow of particles observed at the BNL Relativistic Heavy Ion Collider (RHIC) in relativistic ^{238}U+^{238}U collisions when they rely on information obtained from low-energy experiments for the implementation of deformation in the colliding ^{238}U ions. We show that this is due to an inappropriate treatment of well-deformed nuclei in the modeling of the initial conditions of the quark-gluon plasma. Past studies have identified the deformation of the nuclear surface with that of the nuclear volume, though these are different concepts. In particular, a volume quadrupole moment can be generated by both a surface hexadecapole and a surface quadrupole moment. This feature was so far neglected in the modeling of heavy-ion collisions, and is particularly relevant for nuclei like ^{238}U, which is both quadrupole deformed and hexadecapole deformed. With rigorous input from Skyrme density functional calculations, we show that correcting for such effects in the implementation of nuclear deformations in hydrodynamic simulations restores agreement with BNL RHIC data. This brings consistency to the results of nuclear experiments across energy scales, and demonstrates the impact of the hexadecapole deformation of ^{238}U on high-energy collisions.

Unified Equation of State for Neutron Stars Based on the Gogny Interaction

The effective Gogny interactions of the D1 family were established by D. Gogny more than forty years ago with the aim to describe simultaneously the mean field and the pairing field corresponding to the nuclear interaction. The most popular Gogny parametrizations, namely D1S, D1N and D1M, describe accurately the ground-state properties of spherical and deformed finite nuclei all across the mass table obtained with Hartree–Fock–Bogoliubov (HFB) calculations. However, these forces produce a rather soft equation of state (EoS) in neutron matter, which leads to predict maximum masses of neutron stars well below the observed value of two solar masses. To remove this limitation, we built new Gogny parametrizations by modifying the density dependence of the symmetry energy predicted by the force in such a way that they can be applied to the neutron star domain and can also reproduce the properties of finite nuclei as good as their predecessors. These new parametrizations allow us to obtain stiffer EoS’s based on the Gogny interactions, which predict maximum masses of neutron stars around two solar masses. Moreover, other global properties of the star, such as the moment of inertia and the tidal deformability, are in harmony with those obtained with other well tested EoSs based on the SLy4 Skyrme force or the Barcelona–Catania–Paris–Madrid (BCPM) energy density functional. Properties of the core-crust transition predicted by these Gogny EoSs are also analyzed. Using these new Gogny forces, the EoS in the inner crust is obtained with the Wigner–Seitz approximation in the Variational Wigner–Kirkwood approach along with the Strutinsky integral method, which allows one to estimate in a perturbative way the proton shell and pairing corrections. For the outer crust, the EoS is determined basically by the nuclear masses, which are taken from the experiments, wherever they are available, or by HFB calculations performed with these new forces if the experimental masses are not known.

Taming Nuclear Complexity with a Committee of Multilayer Neural Networks

We demonstrate that a committee of deep neural networks is capable of predicting the ground-state and excited energies of more than 1800 atomic nuclei with an accuracy akin to the one achieved by state-of-the-art nuclear energy density functionals (EDFs) and with significantly less computational cost. An active learning strategy is proposed to train this algorithm with a minimal set of 210 nuclei. This approach enables future fast studies of the influence of EDF parametrizations on structure properties over the whole nuclear chart and suggests that for the first time a machine learning framework successfully encoded several correlated aspects of nuclear deformation.

78Ni revealed as a doubly magic stronghold against nuclear deformation

Nuclear magic numbers correspond to fully occupied energy shells of protons or neutrons inside atomic nuclei. Doubly magic nuclei, with magic numbers for both protons and neutrons, are spherical and extremely rare across the nuclear landscape. Although the sequence of magic numbers is well established for stable nuclei, experimental evidence has revealed modifications for nuclei with a large asymmetry between proton and neutron numbers. Here we provide a spectroscopic study of the doubly magic nucleus ⁷⁸Ni, which contains fourteen neutrons more than the heaviest stable nickel isotope. We provide direct evidence of its doubly magic nature, which is also predicted by ab initio calculations based on chiral effective-field theory interactions and the quasi-particle random-phase approximation. Our results also indicate the breakdown of the neutron magic number 50 and proton magic number 28 beyond this stronghold, caused by a competing deformed structure. State-of-the-art phenomenological shell-model calculations reproduce this shape coexistence, predicting a rapid transition from spherical to deformed ground states, with ⁷⁸Ni as the turning point.

Discovery of ^{60}Ca and Implications For the Stability of ^{70}Ca

The discovery of the important neutron-rich nucleus _{20}^{60}Ca_{40} and seven others near the limits of nuclear stability is reported from the fragmentation of a 345  MeV/u ^{70}Zn projectile beam on ^{9}Be targets at the radioactive ion-beam factory of the RIKEN Nishina Center. The produced fragments were analyzed and unambiguously identified using the BigRIPS two-stage in-flight separator. The eight new neutron-rich nuclei discovered, ^{47}P, ^{49}S, ^{52}Cl, ^{54}Ar, ^{57}K, ^{59,60}Ca, and ^{62}Sc, are the most neutron-rich isotopes of the respective elements. In addition, one event consistent with ^{59}K was registered. The results are compared with the drip lines predicted by a variety of mass models and it is found that the models in best agreement with the observed limits of existence in the explored region tend to predict the even-mass Ca isotopes to be bound out to at least ^{70}Ca.

The Gogny-HFB+QRPA dipole strength function and its application to radiative neutron capture cross section

Valuable theoretical predictions of nuclear dipole excitations in the whole chart are of great interest for different nuclear applications, including in particular nuclear astrophysics. Here we extend our large-scale calculations of the E1 and M1 absorption γ-ray strength function obtained in the framework of the axially-symmetric deformed quasiparticle random phase approximation (QRPA) based on the finite-range D1M Gogny force to the determination of the de-excitation strength function. To do so, shell-model calculations of the de-excitation dipole strength function as well as experimental data are considered to provide insight in the low-energy limit and to complement the QRPA estimate phenomenologically. We compare our final prediction of the E1 and M1 strengths with available experimental data at low energies and show that a relatively good agreement can be obtained. Its impact on the average radiative width as well as radiative neutron capture cross section is discussed.

Impact of new data for neutron-rich heavy nuclei on theoretical models for r-process nucleosynthesis

Current models for the r process are summarized with an emphasis on the key constraints from both nuclear physics measurements and astronomical observations. In particular, we analyze the importance of nuclear physics input such as beta-decay rates; nuclear masses; neutron-capture cross sections; beta-delayed neutron emission; probability of spontaneous fission, beta- and neutron-induced fission, fission fragment mass distributions; neutrino-induced reaction cross sections, etc. We highlight the effects on models for r-process nucleosynthesis of newly measured β-decay half-lives, masses, and spectroscopy of neutron-rich nuclei near the r-process path. We overview r-process nucleosynthesis in the neutrino driven wind above the proto-neutron star in core collapse supernovae along with the possibility of magneto-hydrodynamic jets from rotating supernova explosion models. We also consider the possibility of neutron star mergers as an r-process environment. A key outcome of newly measured nuclear properties far from stability is the degree of shell quenching for neutron rich isotopes near the closed neutron shells. This leads to important constraints on the sites for r-process nucleosynthesis in which freezeout occurs on a rapid timescale.

Microscopic theory of nuclear fission: a review

This article reviews how nuclear fission is described within nuclear density functional theory. A distinction should be made between spontaneous fission, where half-lives are the main observables and quantum tunnelling the essential concept, and induced fission, where the focus is on fragment properties and explicitly time-dependent approaches are often invoked. Overall, the cornerstone of the density functional theory approach to fission is the energy density functional formalism. The basic tenets of this method, including some well-known tools such as the Hartree-Fock-Bogoliubov (HFB) theory, effective two-body nuclear potentials such as the Skyrme and Gogny force, finite-temperature extensions and beyond mean-field corrections, are presented succinctly. The energy density functional approach is often combined with the hypothesis that the time-scale of the large amplitude collective motion driving the system to fission is slow compared to typical time-scales of nucleons inside the nucleus. In practice, this hypothesis of adiabaticity is implemented by introducing (a few) collective variables and mapping out the many-body Schrödinger equation into a collective Schrödinger-like equation for the nuclear wave-packet. The region of the collective space where the system transitions from one nucleus to two (or more) fragments defines what are called the scission configurations. The inertia tensor that enters the kinetic energy term of the collective Schrödinger-like equation is one of the most essential ingredients of the theory, since it includes the response of the system to small changes in the collective variables. For this reason, the two main approximations used to compute this inertia tensor, the adiabatic time-dependent HFB and the generator coordinate method, are presented in detail, both in their general formulation and in their most common approximations. The collective inertia tensor enters also the Wentzel-Kramers-Brillouin (WKB) formula used to extract spontaneous fission half-lives from multi-dimensional quantum tunnelling probabilities (For the sake of completeness, other approaches to tunnelling based on functional integrals are also briefly discussed, although there are very few applications.) It is also an important component of some of the time-dependent methods that have been used in fission studies. Concerning the latter, both the semi-classical approaches to time-dependent nuclear dynamics and more microscopic theories involving explicit quantum-many-body methods are presented. One of the hallmarks of the microscopic theory of fission is the tremendous amount of computing needed for practical applications. In particular, the successful implementation of the theories presented in this article requires a very precise numerical resolution of the HFB equations for large values of the collective variables. This aspect is often overlooked, and several sections are devoted to discussing the resolution of the HFB equations, especially in the context of very deformed nuclear shapes. In particular, the numerical precision and iterative methods employed to obtain the HFB solution are documented in detail. Finally, a selection of the most recent and representative results obtained for both spontaneous and induced fission is presented, with the goal of emphasizing the coherence of the microscopic approaches employed. Although impressive progress has been achieved over the last two decades to understand fission microscopically, much work remains to be done. Several possible lines of research are outlined in the conclusion.

First Measurement of Several β-Delayed Neutron Emitting Isotopes Beyond N=126

The β-delayed neutron emission probabilities of neutron rich Hg and Tl nuclei have been measured together with β-decay half-lives for 20 isotopes of Au, Hg, Tl, Pb, and Bi in the mass region N≳126. These are the heaviest species where neutron emission has been observed so far. These measurements provide key information to evaluate the performance of nuclear microscopic and phenomenological models in reproducing the high-energy part of the β-decay strength distribution. This provides important constraints on global theoretical models currently used in r-process nucleosynthesis.

Spectroscopic Quadrupole Moments in {96,98}Sr: Evidence for Shape Coexistence in Neutron-Rich Strontium Isotopes at N=60

Neutron-rich {96,98}Sr isotopes have been investigated by safe Coulomb excitation of radioactive beams at the REX-ISOLDE facility. Reduced transition probabilities and spectroscopic quadrupole moments have been extracted from the differential Coulomb excitation cross sections. These results allow, for the first time, the drawing of definite conclusions about the shape coexistence of highly deformed prolate and spherical configurations. In particular, a very small mixing between the coexisting states is observed, contrary to other mass regions where strong mixing is present. Experimental results have been compared to beyond-mean-field calculations using the Gogny D1S interaction in a five-dimensional collective Hamiltonian formalism, which reproduce the shape change at N=60.

Uncertainty quantification for nuclear density functional theory and information content of new measurements

Statistical tools of uncertainty quantification can be used to assess the information content of measured observables with respect to present-day theoretical models, to estimate model errors and thereby improve predictive capability, to extrapolate beyond the regions reached by experiment, and to provide meaningful input to applications and planned measurements. To showcase new opportunities offered by such tools, we make a rigorous analysis of theoretical statistical uncertainties in nuclear density functional theory using Bayesian inference methods. By considering the recent mass measurements from the Canadian Penning Trap at Argonne National Laboratory, we demonstrate how the Bayesian analysis and a direct least-squares optimization, combined with high-performance computing, can be used to assess the information content of the new data with respect to a model based on the Skyrme energy density functional approach. Employing the posterior probability distribution computed with a Gaussian process emulator, we apply the Bayesian framework to propagate theoretical statistical uncertainties in predictions of nuclear masses, two-neutron dripline, and fission barriers. Overall, we find that the new mass measurements do not impose a constraint that is strong enough to lead to significant changes in the model parameters. The example discussed in this study sets the stage for quantifying and maximizing the impact of new measurements with respect to current modeling and guiding future experimental efforts, thus enhancing the experiment-theory cycle in the scientific method.

New fission fragment distributions and r-process origin of the rare-earth elements

Neutron star (NS) merger ejecta offer a viable site for the production of heavy r-process elements with nuclear mass numbers A≳140. The crucial role of fission recycling is responsible for the robustness of this site against many astrophysical uncertainties, but calculations sensitively depend on nuclear physics. In particular, the fission fragment yields determine the creation of 110≲A≲170 nuclei. Here, we apply a new scission-point model, called SPY, to derive the fission fragment distribution (FFD) of all relevant neutron-rich, fissioning nuclei. The model predicts a doubly asymmetric FFD in the abundant A≃278 mass region that is responsible for the final recycling of the fissioning material. Using ejecta conditions based on relativistic NS merger calculations, we show that this specific FFD leads to a production of the A≃165 rare-earth peak that is nicely compatible with the abundance patterns in the Sun and metal-poor stars. This new finding further strengthens the case of NS mergers as possible dominant origin of r nuclei with A≳140.

Nuclear correlations and the r process

We show that long-range correlations for nuclear masses have a significant effect on the synthesis of heavy elements by the r process. As calculated by Delaroche et al. [Phys. Rev. C 81, 014303 (2010)], these correlations suppress magic number effects associated with minor shells. This impacts the calculated abundances before the third r-process peak (at mass number A≈195), where the abundances are low and form a trough. This trough and the position of the third abundance peak are strongly affected by the masses of nuclei in the transition region between deformed and spherical. Based on different astrophysical environments, our results demonstrate that a microscopic theory of nuclear masses including correlations naturally smoothens the separation energies, thus reducing the trough and improving the agreement with observed solar system abundances.

Evidence for a smooth onset of deformation in the neutron-rich Kr isotopes

The neutron-rich nuclei 94,96Kr were studied via projectile Coulomb excitation at the REX-ISOLDE facility at CERN. Level energies of the first excited 2(+) states and their absolute E2 transition strengths to the ground state are determined and discussed in the context of the E(2(1)(+)) and B(E2;2(1)(+)→0(1)(+)) systematics of the krypton chain. Contrary to previously published results no sudden onset of deformation is observed. This experimental result is supported by a new proton-neutron interacting boson model calculation based on the constrained Hartree-Fock-Bogoliubov approach using the microscopic Gogny-D1M energy density functional.