Synthesis Of Heavy Elements Research Articles

Nucleosynthesis in the Innermost Ejecta of Magnetorotational Supernova Explosions in Three Dimensions

Core-collapse supernova (CCSN) explosions powered by rotation and magnetic fields present an interesting astrophysical site for nucleosynthesis that potentially contributes to the production of r-process elements. Here we present yields of the innermost ejecta in 3D magnetorotational CCSN models simulated using the CoCoNuT-FMT code. Strong magnetic fields tap the rotational energy of the proto−neutron star and lead to earlier and more energetic (∼3 × 1051 erg) explosions than typical neutrino-driven CCSNe. Compared to a reference nonmagnetic model, the ejecta in the magnetorotational models have much more neutron-rich components with Y e down to ∼0.25. Our post-processing calculations with the reaction network SkyNet show significant production of weak r-process elements up to mass number ∼130. We find negligible differences in the synthesis of heavy elements between two magnetorotational models with different initial field strengths of 1010 and 1012 G, in accord with their similar explosion dynamics. The magnetorotational models produce about ∼0.19 and 0.14 M ☉ of radioactive 56Ni, on the low end of inferred hypernova nickel masses. The yields are publicly available at Zenodo (doi: 10.5281/zenodo.10578981) for comparison with stellar abundance patterns, inclusion in modeling galactic chemical evolution, and comparison with other yield calculations. Our results add to the yet-restricted corpus of nucleosynthesis yields from 3D magnetorotational supernova simulations and will help quantify yield uncertainties.

Inferences Regarding a Stellar Metallicity Distribution Following a Single Primordial Chemical Enrichment Event

An idealized scenario is presented according to which a power-law distribution in metallicity might have characterized the first metal-poor stars. The scenario assumes an initial synthesis of heavy elements within a first source of metals, followed by a subsequent ejection of these metals and their dilution within a metal-free gas environment.

Neutron star mergers as the dominant contributor to the production of heavy r-process elements

ABSTRACT The discovery of the radioactively powered kilonova AT2017gfo, associated with the short-duration gamma-ray burst GRB 170817A and the gravitational wave source GW170817, has provided the first direct evidence supporting binary neutron star mergers as crucial astrophysical sites for the synthesis of heavy elements beyond iron through r-process nucleosysthesis in the universe. However, recent identification of kilonovae following long-duration gamma-ray bursts, such as GRB 211211A and GRB 230307A, has sparked discussions about the potential of neutron star–white dwarf mergers to also produce neutron-rich ejecta and contribute to the production of heavy r-process elements. In this work, we estimate the contribution of binary neutron star mergers to the total mass of r-process elements in the Milky Way and investigate the possibility of neutron star–white dwarf mergers as alternative astrophysical sites for r-process nucleosynthesis through an analysis of the total mass of the r-process elements in the Milky Way. Our results reveal that binary neutron star mergers can sufficiently account for the Galactic heavy r-process elements, suggesting that these events are the dominant contributor to the production of heavy r-process elements in the Milky Way. Considering the total mass of r-process elements in the Milky Way and the higher occurrence rate of neutron star–white dwarf mergers, it is unlikely that such mergers can produce a significant amount of neutron-rich ejecta, with the generated mass of r-process elements being lower than $0.005\, {\rm M}_{\odot }$.

Multimessenger Diagnostics of the Engine behind Core-collapse Supernovae

Core-collapse supernova explosions play a wide role in astrophysics by producing compact remnants (neutron stars or black holes) and the synthesis and injection of many heavy elements into their host galaxy. Because they are produced in some of the most extreme conditions in the universe, they can also probe physics in extreme conditions (matter at nuclear densities and extreme temperatures and magnetic fields). To quantify the impact of supernovae on both fundamental physics and our understanding of the universe, we must leverage a broad set of observables of this engine. In this paper, we study a subset of these probes using a suite of one-dimensional, parameterized mixing models: ejecta remnants from supernovae, ultraviolet, optical and infrared light curves, and transient gamma-ray emission. We review the other diagnostics and show how the different probes tie together to provide a more clear picture of the supernova engine. Join us in improving and evolving this document through active community engagement. Instructions are provided at this link: https://github.com/clfryer/MM-SNe.

Collapse of Rotating Massive Stars Leading to Black Hole Formation and Energetic Supernovae

We explore a possible explosion scenario resulting from core collapses of rotating massive stars that leave a black hole by performing radiation-viscous-hydrodynamics simulations in numerical relativity. We take moderately and rapidly rotating compact pre-collapse stellar models with zero-age main-sequence masses of 9M ⊙ and 20M ⊙ based on stellar evolution calculations as the initial conditions. We find that viscous heating in the disk formed around the central black hole is the power source for an outflow. The moderately rotating models predict a small ejecta mass of the order of 0.1M ⊙ and an explosion energy of ≲1051 erg. Due to the small ejecta mass, these models may predict a short-timescale transient with a rise time of 3–5 days. This can lead to a bright (∼1044 erg s−1) transient, like superluminous supernovae in the presence of a dense massive circumstellar medium. For hypothetically rapidly rotating models that have a high mass-infall rate onto the disk, the explosion energy is ≳3 × 1051 erg, which is comparable to or larger than that of typical stripped-envelope supernovae, indicating that a fraction of such supernovae may be explosions powered by black hole accretion disks. The explosion energy is still increasing at the end of the simulations with a rate of >1050 erg s−1, and thus, it may reach ∼1052 erg. A nucleosynthesis calculation shows that the mass of 56Ni amounts to ≳0.1M ⊙, which, together with the high explosion energy, may satisfy the required amount for broad-lined type Ic supernovae. Irrespective of the models, the lowest value of the electron fraction of the ejecta is ≳0.4; thus, synthesis of heavy r-process elements is not found in our models.

Hydrodynamic Mixing of Accretion Disk Outflows in Collapsars: Implications for r-process Signatures

The astrophysical environments capable of triggering heavy-element synthesis via rapid neutron capture (the r-process) remain uncertain. While binary neutron star mergers (NSMs) are known to forge r-process elements, certain rare supernovae (SNe) have been theorized to supplement—or even dominate—r-production by NSMs. However, the most direct evidence for such SNe, unusual reddening of the emission caused by the high opacities of r-process elements, has not been observed. Recent work identified the distribution of r-process material within the SN ejecta as a key predictor of the ease with which signals associated with r-process enrichment could be discerned. Though this distribution results from hydrodynamic processes at play during the SN explosion, thus far it has been treated only in a parameterized way. We use hydrodynamic simulations to model how disk winds—the alleged locus of r-production in rare SNe—mix with initially r-process-free ejecta. We study mixing as a function of the wind mass, wind duration, and the initial SN explosion energy, and find that it increases with the first two of these and decreases with the third. This suggests that SNe accompanying the longest long-duration gamma-ray bursts are promising places to search for signs of r-process enrichment. We use semianalytic radiation transport to connect hydrodynamics to electromagnetic observables, allowing us to assess the mixing level at which the presence of r-process material can be diagnosed from SN light curves. Analytic arguments constructed atop this foundation imply that a wind-driven r-process-enriched SN model is unlikely to explain standard energetic SNe.

Superheavy Elements in Kilonovae

As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first “kilonova” associated with a binary neutron star merger (NSM) suggests that these sites are hosts of the rapid neutron capture (“r”) process. However, it is unknown just how robust elemental production can be in mergers. Identifying signposts of the production of particular nuclei is critical for fully understanding merger-driven heavy-element synthesis. In this study, we investigate the properties of very neutron-rich nuclei for which superheavy elements (Z ≥ 104) can be produced in NSMs and whether they can similarly imprint a unique signature on kilonova light-curve evolution. A superheavy-element signature in kilonovae represents a route to establishing a lower limit on heavy-element production in NSMs as well as possibly being the first evidence of superheavy-element synthesis in nature. Favorable NSM conditions yield a mass fraction of superheavy elements X Z≥104 ≈ 3 × 10−2 at 7.5 hr post-merger. With this mass fraction of superheavy elements, we find that the component of kilonova light curves possibly containing superheavy elements may appear similar to those arising from lanthanide-poor ejecta. Therefore, photometric characterizations of superheavy-element rich kilonova may possibly misidentify them as lanthanide-poor events.

Competition between fusion and quasifission in the angular momentum dependent dynamics of heavy element synthesis reactions

Background: Mass and angle distribution measurements have illuminated many aspects of the physical variables controlling quasifission. However, mapping from detection angle to reaction time is clouded by the wide range of contributing angular momenta $L\ensuremath{\hbar}$, ranging from $0\ensuremath{\hbar}$ to the maximum of the reaction ($⪆100\ensuremath{\hbar}$), which complicates the mapping, and thus limits our understanding of the reaction dynamics.Purpose: To investigate the angular momentum dependence of the reaction dynamics in quasifission and determine the fission fragment mass evolutions in the reaction.Method: The mass and angular distributions of products of the reactions $^{52}\mathrm{Cr}+^{198}\mathrm{Pt}$ and $^{54}\mathrm{Cr}+^{196}\mathrm{Pt}$ were measured. The distributions were compared with quasifission simulation results.Results: Mass angle distributions are reproduced by utilizing a new quasifission mass evolution model, and including a mass-symmetric component associated with low $L$. The latter increases in yield at higher beam energy.Conclusions: The symmetric component represents the total of slow quasifission and fusion-fission. The increase in contribution with energy above the Coulomb barrier suggests that an extra-push energy is required to achieve fusion.

S-process Nuclear Reaction Rates

In stars the 13C(α, n)16O and 22Ne(α, n)25Mg reactions are the two main sources of neutrons for the so-called slow neutron capture process (s-process), which is one of the main mechanisms for the stellar synthesis of heavy elements. About 13C(α, n)16O, in despite of many efforts in measuring its cross section at the lower energies, only high uncertainty data above the s-process Gamow window (150 keV < Ecm < 230 keV) were available, due mostly to the difficulties on suppress the natural background. Indeed, only recently the LUNA collaboration performed high precision underground measurements of the reaction cross section inside the Gamow window, improving the accuracy of its extrapolation at the lower energies. Again due to natural background, only upper limits for the 22Ne(α, n)25Mg reaction cross section are currently known in the s-process Gamow window (450 keV < Ecm < 750 keV). For this, the ERC founded project SHADES (Unina/INFN) aims to perform high precision and high sensitivity measurements of the 22Ne(α, n)25Mg reaction cross section down to the neutron threshold. A sensitivity improvement of at least two orders of magnitude over the state of the art is expected thanks to the low natural background environment of the INFN-LNGS laboratory in Italy, the high beam current of the new LUNAMV accelerator and the Beam Induced Background events suppression performed by SHADES hybrid detectors array. In this paper I will present the LUNA efforts to estimate nuclear reaction rates for 13C(α, n)16O, with a focus on the R-Matrix analysis performed with the code AZURE2 to extrapolate nuclear reaction rates at stellar energies and the estimate of their uncertainty through Monte Carlo methods. I will also present an overview of the SHADES project to measure 22Ne(α, n)25Mg in the Gamow window and the first results on the setup commissioning.

Nucleon-nucleon interaction and the neutrino surface of neutron star mergers

Neutrinos act as an important cooling mechanism in an array of explosive astrophysical events. Modeling of their fluxes and spectra is key for the interpretation of their detection and the understanding of the synthesis of heavy elements. The neutrino surface is key in the modeling of neutrino spectra. In previous work, the points of last neutrino scattering were found assuming no interaction between nucleons. In this work we incorporate this interaction using a few different models and observe differences in where the neutrinos decouple from the system.

Indirect measurements of neutron-induced reaction cross sections at storage rings

Neutron-induced reaction cross sections of unstable nuclei are essential for understanding the synthesis of heavy elements in stars. However, their measurement is very difficult due to the radioactivity of the targets involved. We propose to circumvent this problem by using for the first time the surrogate reaction method in inverse kinematics at heavy-ion storage rings. In this contribution, we describe the developments we have done to perform surrogate-reaction studies at the storage rings of GSI/FAIR. In particular, we present the first results of the proof of principle experiment, which we conducted recently at the Experimental Storage Ring (ESR).

Magnetorotational supernovae: a nucleosynthetic analysis of sophisticated 3D models

ABSTRACT Magnetorotational supernovae are a rare type of core-collapse supernovae where the magnetic field and rotation play a central role in the dynamics of the explosion. We present the post-processed nucleosynthesis of state-of-the-art neutrino-MHD supernova models that follow the post explosion evolution for few seconds. We find three different dynamical mechanisms to produce heavy r-process elements: (i) a prompt ejection of matter right after core bounce, (ii) neutron-rich matter that is ejected at late times due to a reconfiguration of the protoneutronstar shape, (iii) small amount of mass ejected with high entropies in the centre of the jet. We investigate total ejecta yields, including the ones of unstable nuclei such as 26Al, 44Ti, 56Ni, and 60Fe. The obtained 56Ni masses vary between $0.01\!-\!1\, \mathrm{M_\odot }$. The latter maximum is compatible with hypernova observations. Furthermore, all of our models synthesize Zn masses in agreement with observations of old metal-poor stars. We calculate simplified light curves to investigate whether our models can be candidates for superluminous supernovae. The peak luminosities obtained from taking into account only nuclear heating reach up to a few $\sim 10^{43} \, \mathrm{erg\, s^{-1}}$. Under certain conditions, we find a significant impact of the 66Ni decay chain that can raise the peak luminosity up to $\sim 38{{\ \rm per \, cent}}$ compared to models including only the 56Ni decay chain. This work reinforces the theoretical evidence on the critical role of magnetorotational supernovae to understand the occurrence of hypernovae, superluminous supernovae, and the synthesis of heavy elements.

Development of Micro-Pattern Gaseous Detectors for Nuclear Reaction Studies

One of the frontiers of today’s nuclear physics research is the synthesis of Super Heavy Elements (SHE). Fusion-fission dynamics, namely the competition between quasi fission and fusion is one of the key challenges to optimize the SHE. To have an insight into the dynamics, one requires the study of fission fragment mass and angular distribution near barrier energies for heavy-ion induced fission reactions. Recent successful installation of linear accelerators in India offers a unique opportunity to study the dynamics of nuclear reactions and formation process of SHE. For the effective utilization of these current, as well as upcoming facilities, development of novel detectors to study reaction dynamics, formation process of SHE with heavier projectiles and higher beam energies is needed. Gaseous detectors have undergone a rapid improvement in terms of spatial, temporal and energy resolution, rate capability, radiation hardness, ion feedback etc., ushering in a new genre of micro-structured devices based on semi-conductor technology, commonly known as Micro-Pattern Gaseous Detectors (MPGDs). Although many of the MPGD structures were primarily developed for high-rate tracking of charged particles in high energy physics experiments, stability of operation, simplicity of construction and relatively low cost make these detectors suitable for other applications, such as low-energy nuclear physics experiments. The present activities encompass a detailed evaluation of the operational conditions of Micromesh-Multi Wire and THGEM-Multi Wire hybrid detector operated in low-pressure isobutane gas with a view to optimizing their use in the detection of charged particles and fission fragments.

Planetary Nebulae as Sources of Chemical Enrichment of the Galaxy

During the late asymptotic giant branch (AGB) phase of stellar evolution, the element carbon is created through the triple-α nuclear reaction. This is followed by the synthesis of other heavy elements via neutron capture in the core and the dredge-up of these elements to the surface by convection. Simple molecules and solid-state minerals begin to form in the upper photosphere. These molecules and solids are ejected into the interstellar medium by a strong stellar wind. During the subsequent post-AGB evolution, complex organics with aromatic and aliphatic structures are synthesized in the circumstellar envelope. Planetary nebulae, formed by swept-up circumstellar material, are found to show strong spectral signatures of ions, atoms, molecules, and inorganic and organic solids. These ejected materials of planetary nebulae serve as primordial ingredients of formation of new stars and planetary systems. Stellar synthesized organic solids may survive their journeys through the interstellar medium and evidence for remnants of such organics can be found in our own Solar System. In this paper, we summarize the recent observations of circumstellar synthesis of molecules and solids and discuss the implications of these ejected circumstellar materials on the chemical enrichment of the Galaxy and planetary systems.

Systematic exploration of heavy element nucleosynthesis in protomagnetar outflows

ABSTRACTWe study the nucleosynthesis products in neutrino-driven winds from rapidly rotating, highly magnetized and misaligned protomagnetars using the nuclear reaction network SkyNet. We adopt a semi-analytic parametrized model for the protomagnetar and systematically study the capabilities of its neutrino-driven wind for synthesizing nuclei and eventually producing ultra-high energy cosmic rays (UHECRs). We find that for neutron-rich outflows (Ye < 0.5), synthesis of heavy elements ($\overline{A}\sim 20-65$) is possible during the first $\sim 10\, {\rm s}$ of the outflow, but these nuclei are subjected to composition-altering photodisintegration during the epoch of particle acceleration at the dissipation radii. However, after the first $\sim 10\, {\rm s}$ of the outflow, nucleosynthesis reaches lighter elements ($\overline{A}\sim 10-50$) that are not subjected to subsequent photodisintegration. For proton-rich (Ye ≥ 0.5) outflows, synthesis is more limited ($\overline{A}\sim 4-15$). These suggest that while protomagnetars typically do not synthesize nuclei heavier than second r-process peak elements, they are intriguing sources of intermediate/heavy mass UHECRs. For all configurations, the most rapidly rotating protomagnetars are more conducive for nucleosynthesis with a weaker dependence on the magnetic field strength.

Mass Equilibration and Fluctuations in the Angular Momentum Dependent Dynamics of Heavy Element Synthesis Reactions.

Mass and angle distributions for the ^{52}Cr+^{198}Pt and ^{54}Cr+^{196}Pt reactions (both forming ^{250}No) were measured and subtracted, giving new information on fast quasifission mass evolution, and the first direct determination of the dependence of sticking times on angular momentum. TDHF calculations showed good agreement with average experimental values, but experimental mass distributions unexpectedly extended to symmetric splits while the peak yield remained close to the initial masses. This implies a strong role of fluctuations in mass division early in the collision, giving insights into the transition from fast energy dissipative deep-inelastic collisions to quasifission.

Measurement of the 72Ge(n,γ) cross section over a wide neutron energy range at the CERN n_TOF facility

The Ge72(n,γ) cross section was measured for neutron energies up to 300keV at the neutron time-of-flight facility n_TOF (CERN), Geneva, for the first time covering energies relevant to heavy-element synthesis in stars. The measurement was performed at the high-resolution beamline EAR-1, using an isotopically enriched GeO272 sample. The prompt capture γ rays were detected with four liquid scintillation detectors, optimized for low neutron sensitivity. We determined resonance capture kernels up to a neutron energy of 43keV, and averaged cross sections from 43 to 300keV. Maxwellian-averaged cross section values were calculated from kT=5 to 100keV, with uncertainties between 3.2% and 7.1%. The new results significantly reduce uncertainties of abundances produced in the slow neutron capture process in massive stars.

The Study of Key Reactions Shaping the Post-Main Sequence Evolution of Massive Stars in Underground Facilities

The chemical evolution of the Universe and several phases of stellar life are regulated by minute nuclear reactions. The key point for each of these reactions is the value of cross-sections at the energies at which they take place in stellar environments. Direct cross-section measurements are mainly hampered by the very low counting rate and by cosmic background; nevertheless, they have become possible by combining the best experimental techniques with the cosmic silence of an underground laboratory. In the nineties, the LUNA (Laboratory for Underground Nuclear Astrophysics) collaboration opened the era of underground nuclear astrophysics, installing first a homemade 50 kV and, later on, a second 400 kV accelerator under the Gran Sasso mountain in Italy: in 25 years of experimental activity, important reactions responsible for hydrogen burning could have been studied down to the relevant energies thanks to the high current proton and helium beams provided by the machines. The interest in the next and warmer stages of star evolution (i.e., post-main sequence and helium and carbon burning) drove a new project based on an ion accelerator in the MV range called LUNA-MV, able to deliver proton, helium, and carbon beams. The present contribution is aimed to discuss the state of the art for some selected key processes of post-main sequence stellar phases:12C(α,γ)16O and12C+12C are fundamental for helium and carbon burning phases, and13C(α,n)16O and22Ne(α,n)25Mg are relevant to the synthesis of heavy elements in AGB stars. The perspectives opened by an underground MV facility will be highlighted.

The origin of r-process elements

Nuclear Astrophysics Theoretical models predict that the synthesis of heavy elements by the rapid neutron capture process (r-process) occurs in extreme astrophysical environments such as neutron star mergers or some types of supernovae. Testing those predictions by comparing them with the isotopic record has been difficult. Cote et al. examined two r-process isotopes, iodine-129 and curium-247, both of which have half-lives of 15.6 million years. Therefore, their ratio remains constant even long after the nucleosynthesis event. The ratio of those isotopes at the time of Solar System formation is recorded in meteorites. Comparing this value with nuclear astrophysics calculations shows that the most likely source was moderately neutron-rich material ejected from a binary neutron star merger. Science , this issue p. [945][1] [1]: /lookup/doi/10.1126/science.aba1111

Nucleosynthesis in magneto-rotational supernovae

Abstract We present the nucleosynthesis of magneto-rotational supernovae (MR-SNe) including neutrino-driven and magneto-rotational-driven ejecta based, for the first time, on two-dimensional simulations with accurate neutrino transport. The models analysed here have different rotation and magnetic fields, allowing us to explore the impact of these two key ingredients. The accurate neutrino transport of the simulations is critical to analyse the slightly neutron rich and proton rich ejecta that are similar to the, also neutrino-driven, ejecta in standard supernovae. In the model with strong magnetic field, the r-process produces heavy elements up to the third r-process peak (A ∼ 195), in agreement with previous works. This model presents a jet-like explosion with proton-rich jets surrounded by neutron rich material where the r-process occurs. We have estimated a lower limit for 56Ni of 2.5 × 10−2M⊙, which is still well below the expected hypernova value. Longer simulations including the accretion disk evolution are required to get a final prediction. In addition, we have found that the late evolution is critical in a model with weak magnetic field in which lately ejected neutron rich matter produces elements up to the second r-process peak. Even if we cannot yet provide conclusions for hypernova nucleosynthesis, our results agree with observations of old stars and radioactive isotopes in supernova remnants. This makes MR-SNe a good additional scenario to neutron star mergers for the synthesis of heavy elements and brings us closer to understand their origin and the role of MR-SNe in the early galaxy nucleosynthesis.