THE ROLE OF FISSION IN NEUTRON STAR MERGERS AND ITS IMPACT ON THEr-PROCESS PEAKS

β-Delayed One and Two Neutron Emission Probabilities Southeast of ^{132}Sn and the Odd-Even Systematics in r-Process Nuclide Abundances

The β-delayed one- and two-neutron emission probabilities (P_{1n} and P_{2n}) of 20 neutron-rich nuclei with N≥82 have been measured at the RIBF facility of the RIKEN Nishina Center. P_{1n} of ^{130,131}Ag, ^{133,134}Cd, ^{135,136}In, and ^{138,139}Sn were determined for the first time, and stringent upper limits were placed on P_{2n} for nearly all cases. β-delayed two-neutron emission (β2n) was unambiguously identified in ^{133}Cd and ^{135,136}In, and their P_{2n} were measured. Weak β2n was also detected from ^{137,138}Sn. Our results highlight the effect of the N=82 and Z=50 shell closures on β-delayed neutron emission probability and provide stringent benchmarks for newly developed macroscopic-microscopic and self-consistent global models with the inclusion of a statistical treatment of neutron and γ emission. The impact of our measurements on r-process nucleosynthesis was studied in a neutron star merger scenario. Our P_{1n} and P_{2n} have a direct impact on the odd-even staggering of the final abundance, improving the agreement between calculated and observed Solar System abundances. The odd isotope fraction of Ba in r-process-enhanced (r-II) stars is also better reproduced using our new data.

Beyond-mean-field calculations of allowed and first-forbidden β− decays of r-process waiting-point nuclei

β-decay rates of neutron-rich nuclei, in particular those located at neutron shell closures, play a central role in simulations of the heavy-element nucleosynthesis and resulting abundance distributions. We present β-decay half-lives of even-even N = 82 and N = 126 r-process waiting-point nuclei calculated in the approach based on relativistic quasiparticle random phase approximation with quasiparticle-vibration coupling. The calculations include both allowed and first-forbidden transitions. In the N = 82 chain, the quasiparticlevibration coupling has an important impact close to stability, as it increases the contribution of Gamow-Teller modes and improves the agreement with the available data. In the N = 126 chain, we find the decay to proceed dominantly via first-forbidden transitions, even when the coupling to vibrations is included.

129I and 247Cm in meteorites constrain the last astrophysical source of solar r-process elements

The composition of the early Solar System can be inferred from meteorites. Many elements heavier than iron were formed by the rapid neutron capture process (r-process), but the astrophysical sources where this occurred remain poorly understood. We demonstrate that the near-identical half-lives [Formula: see text] of the radioactive r-process nuclei iodine-129 and curium-247 preserve their ratio, irrespective of the time between production and incorporation into the Solar System. We constrain the last r-process source by comparing the measured meteoritic ratio ¹²⁹I/²⁴⁷Cm = 438 ± 184 with nucleosynthesis calculations based on neutron star merger and magneto-rotational supernova simulations. Moderately neutron-rich conditions, often found in merger disk ejecta simulations, are most consistent with the meteoritic value. Uncertain nuclear physics data limit our confidence in this conclusion.

Following Fission Products in Explosive Astrophysical Environments

The astrophysical process by which the heaviest elements are formed in the universe is known as the rapid neutron capture process, or r process, of nucleosynthesis. The r process is characterized by the neutron capture and β− decay of short-lived, neutron-rich atomic nuclei; in suitably extreme environments, nuclear fission can also play a major role in determining the ensuing nucleosynthesis. In this work, we present the application of our recently developed nucleosynthesis tracing framework to precisely quantify the impact that neutron-induced and β− -delayed fission processes have in r-process environments that produce fissioning nuclei.

Probing the fission properties of neutron-rich actinides with the astrophysical r process

We review recent work examining the influence of fission in rapid neutron capture (r-process) nucleosynthesis which can take place in astrophysical environments. We briefly discuss the impact of uncertain fission barriers and fission rates on the population of heavy actinide species. We demonstrate the influence of the fission fragment distributions for neutron-rich nuclei and discuss currently available treatments, including recent macroscopic-microscopic calculations. We conclude by comparing our nucleosynthesis results directly with stellar data for metal-poor stars rich in r-process elements to consider whether fission plays a role in the so-called ‘universality’ of r-process abundances observed from star to star.

Kilonovae

The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergoes rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a "kilonova", which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ∼ hour-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.

Bayesian Evaluation of Incomplete Fission Yields

Fission product yields are key infrastructure data for nuclear applications in many aspects. It is a challenge both experimentally and theoretically to obtain accurate and complete energy-dependent fission yields. We apply the Bayesian neural network (BNN) approach to learn existing fission yields and predict unknowns with uncertainty quantification. We demonstrated that the BNN is particularly useful for evaluations of fission yields when incomplete experimental data are available. The BNN evaluation results are quite satisfactory on distribution positions and energy dependencies of fission yields.

Review on the progress in nuclear fission-experimental methods and theoretical descriptions

An overview is given on some of the main advances in the experimental methods, experimental results, theoretical models and ideas of the last few years in the field of nuclear fission. New approaches have considerably extended the availability of fissioning systems for the experimental study of nuclear fission, and have provided a full identification of all fission products in A and Z for the first time. In particular, the transition from symmetric to asymmetric fission around ²²⁶Th, some unexpected structures in the mass distributions in the fission of systems around Z  =  80-84, and an extended systematics of the odd-even effect in the fission fragment Z distributions have all been measured (Andreyev et al 2018 Rep. Prog. Phys. 81 016301). Three classes of model descriptions of fission presently appear to be the most promising or the most successful. Self-consistent quantum-mechanical models fully consider the quantum-mechanical features of the fission process. Intense efforts are presently being made to develop suitable theoretical tools (Schunck and Robledo 2016 Rep. Prog. Phys. 79 116301) for modeling the non-equilibrium, large-amplitude collective motion leading to fission. Stochastic models provide a fully developed technical framework. The main features of the fission-fragment mass distribution have been well reproduced from mercury to fermium and beyond (Möller and Randrup 2015 Phys. Rev. C 91 044316). However, limited computer resources still impose restrictions, for example, on the number of collective coordinates and on an elaborate description of the fission dynamics. In an alternative semi-empirical approach (Schmidt et al 2016 Nucl. Data Sheets 131 107), considerable progress in describing the fission observables has been achieved by combining several theoretical ideas, which are essentially well known. This approach exploits (i) the topological properties of a continuous function in multidimensional space, (ii) the separability of the influence of fragment shells and the macroscopic properties of the compound nucleus, (iii) the properties of a quantum oscillator coupled to a heat bath of other nuclear degrees of freedom, (iv) an early freeze-out of collective motion, and (v) the application of statistical mechanics for describing the thermalization of intrinsic excitations in the nascent fragments. This new approach reveals a high degree of regularity and allows the calculation of high-quality data that is relevant to nuclear technology without specifically adjusting the empirical data of individual systems.

Primordial Black Holes and r-Process Nucleosynthesis

We show that some or all of the inventory of r-process nucleosynthesis can be produced in interactions of primordial black holes (PBHs) with neutron stars (NSs) if PBHs with masses 10^{-14}  M_{⊙}<M_{PBH}<10^{-8}  M_{⊙} make up a few percent or more of dark matter. A PBH captured by a NS sinks to the center of the NS and consumes it from the inside. When this occurs in a rotating millisecond-period NS, the resulting spin-up ejects ∼0.1  M_{⊙}-0.5  M_{⊙} of relatively cold neutron-rich material. This ejection process and the accompanying decompression and decay of nuclear matter can produce electromagnetic transients, such as a kilonova-type afterglow and fast radio bursts. These transients are not accompanied by significant gravitational radiation or neutrinos, allowing such events to be differentiated from compact object mergers occurring within the distance sensitivity limits of gravitational-wave observatories. The PBH-NS destruction scenario is consistent with pulsar and NS statistics, the dark-matter content, and spatial distributions in the Galaxy and ultrafaint dwarfs, as well as with the r-process content and evolution histories in these sites. Ejected matter is heated by beta decay, which leads to emission of positrons in an amount consistent with the observed 511-keV line from the Galactic center.

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.

Kilonovae

The mergers of double neutron star (NS-NS) and black hole (BH)-NS binaries are promising gravitational wave (GW) sources for Advanced LIGO and future GW detectors. The neutron-rich ejecta from such merger events undergoes rapid neutron capture (r-process) nucleosynthesis, enriching our Galaxy with rare heavy elements like gold and platinum. The radioactive decay of these unstable nuclei also powers a rapidly evolving, supernova-like transient known as a "kilonova" (also known as "macronova"). Kilonovae are an approximately isotropic electromagnetic counterpart to the GW signal, which also provides a unique and direct probe of an important, if not dominant, r-process site. I review the history and physics of kilonovae, leading to the current paradigm of week-long emission with a spectral peak at near-infrared wavelengths. Using a simple light curve model to illustrate the basic physics, I introduce potentially important variations on this canonical picture, including: [Formula: see text]day-long optical ("blue") emission from lanthanide-free components of the ejecta; [Formula: see text]hour-long precursor UV/blue emission, powered by the decay of free neutrons in the outermost ejecta layers; and enhanced emission due to energy input from a long-lived central engine, such as an accreting BH or millisecond magnetar. I assess the prospects of kilonova detection following future GW detections of NS-NS/BH-NS mergers in light of the recent follow-up campaign of the LIGO binary BH-BH mergers.

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.

Impact of Nuclear Mass Uncertainties on the r Process

Nuclear masses play a fundamental role in understanding how the heaviest elements in the Universe are created in the r process. We predict r-process nucleosynthesis yields using neutron capture and photodissociation rates that are based on the nuclear density functional theory. Using six Skyrme energy density functionals based on different optimization protocols, we determine for the first time systematic uncertainty bands-related to mass modeling-for r-process abundances in realistic astrophysical scenarios. We find that features of the underlying microphysics make an imprint on abundances especially in the vicinity of neutron shell closures: Abundance peaks and troughs are reflected in trends of neutron separation energy. Further advances in the nuclear theory and experiments, when linked to observations, will help in the understanding of astrophysical conditions in extreme r-process sites.