Nucleosynthesis Of Heavy Elements Research Articles

Kilonova Spectral Inverse Modelling with Simulation-based Inference: An Amortized Neural Posterior Estimation Analysis

Kilonovae represent a category of astrophysical transients, identifiable as the electromagnetic (EM) counterparts associated with the coalescence events of binary systems comprising neutron stars and neutron star–black hole pairs. They act as probes for heavy-element nucleosynthesis in astrophysical environments. These studies rely on an inference of the physical parameters (e.g., ejecta mass, velocity, composition) that describe kilonovae-based on EM observations. This is a complex inverse problem typically addressed with sampling-based methods such as Markov Chain Monte Carlo or nested sampling algorithms. However, repeated inferences can be computationally expensive, due to the sequential nature of these methods. This poses a significant challenge to ensuring the reliability and statistical validity of the posterior approximations and, thus, the inferred kilonova parameters themselves. We present a novel approach: simulation-based inference using simulations produced by KilonovaNet. Our method employs an ensemble of amortized neural posterior estimation (ANPE) with an embedding network to directly predict posterior distributions from simulated spectral energy distributions. We take advantage of the quasi-instantaneous inference time of ANPE to demonstrate the reliability of our posterior approximations using diagnostics tools, including coverage diagnostic and posterior predictive checks. We further test our model with real observations from AT 2017gfo, the only kilonova with multimessenger data, demonstrating agreement with previous likelihood-based methods while reducing inference time down to a few seconds. The inference results produced by ANPE appear to be conservative and reliable, paving the way for testable and more efficient kilonova parameter inference.

Nuclear Uncertainties Associated with the Nucleosynthesis in Ejecta of a Black Hole Accretion Disk

The simulation of heavy element nucleosynthesis requires input from yet-to-be-measured nuclear properties. The uncertainty in the values of these off-stability nuclear properties propagates to uncertainties in the predictions of elemental and isotopic abundances. However, for any given astrophysical explosion, there are many different trajectories, i.e., temperature and density histories, experienced by outflowing material, and thus different nuclear properties can come into play. We consider combined nucleosynthesis results from 460,000 trajectories from a black hole accretion disk and find the spread in elemental predictions due solely to unknown nuclear properties to be a factor of a few. We analyze this relative spread in model predictions due to nuclear variations and conclude that the uncertainties can be attributed to a combination of properties in a given region of the abundance pattern. We calculate a cross-correlation between mass changes and abundance changes to show how variations among the properties of participating nuclei may be explored. Our results provide further impetus for measurements of multiple quantities on individual short-lived neutron-rich isotopes at modern experimental facilities.

A Hubble Space Telescope Search for r-Process Nucleosynthesis in Gamma-Ray Burst Supernovae

The existence of a secondary (in addition to compact object mergers) source of heavy element (r-process) nucleosynthesis, the core-collapse of rapidly rotating and highly magnetized massive stars, has been suggested by both simulations and indirect observational evidence. Here, we probe a predicted signature of r-process enrichment, a late-time (≳40 days post-burst) distinct red color, in observations of gamma-ray burst supernovae (GRB-SNe), which are linked to these massive star progenitors. We present optical to near-IR color measurements of four GRB-SNe at z ≲ 0.4, extending out to >500 days post-burst, obtained with the Hubble Space Telescope and large-aperture ground-based telescopes. Comparison of our observations to models indicates that GRBs 030329, 100316D, and 130427A are consistent with both no enrichment and producing 0.01–0.15 M ⊙ of r-process material if there is a low amount of mixing between the inner r-process ejecta and outer supernova (SN) layers. GRB 190829A is not consistent with any models with r-process enrichment ≥0.01 M ⊙. Taken together the sample of GRB-SNe indicates color diversity at late times. Our derived yields from GRB-SNe may be underestimated due to r-process material hidden in the SN ejecta (potentially due to low mixing fractions) or the limits of current models in measuring r-process mass. We conclude with recommendations for future search strategies to observe and probe the full distribution of r-process produced by GRB-SNe.

Collective Neutrino Oscillations and Heavy-element Nucleosynthesis in Supernovae: Exploring Potential Effects of Many-body Neutrino Correlations

In high-energy astrophysical processes involving compact objects, such as core-collapse supernovae or binary neutron star mergers, neutrinos play an important role in the synthesis of nuclides. Neutrinos in these environments can experience collective flavor oscillations driven by neutrino–neutrino interactions, including coherent forward scattering and incoherent (collisional) effects. Recently, there has been interest in exploring potential novel behaviors in collective oscillations of neutrinos by going beyond the one-particle effective or “mean-field” treatments. Here, we seek to explore implications of collective neutrino oscillations, in the mean-field treatment and beyond, for the nucleosynthesis yields in supernova environments with different astrophysical conditions and neutrino inputs. We find that collective oscillations can impact the operation of the ν p-process and r-process nucleosynthesis in supernovae. The potential impact is particularly strong in high-entropy, proton-rich conditions, where we find that neutrino interactions can nudge an initial ν p-process neutron-rich, resulting in a unique combination of proton-rich low-mass nuclei as well as neutron-rich high-mass nuclei. We describe this neutrino-induced neutron-capture process as the “ν i-process.” In addition, nontrivial quantum correlations among neutrinos, if present significantly, could lead to different nuclide yields compared to the corresponding mean-field oscillation treatments, by virtue of modifying the evolution of the relevant one-body neutrino observables.

Modeling the Light Curve of HD 222925's Progenitor Kilonova

The leading theory for heavy element nucleosynthesis is that a majority of isotopes heavier than iron are formed in neutron star mergers (NSMs), visible through the kilonovae that proceed them. In this work, we derive the potential kilonova light curve associated with the elemental abundances of the “gold standard” metal-poor star HD 222925, which is thought to be enriched by a single NSM. Previous work explores which initial kilonova conditions best match the elemental abundances in HD 222925 and finds two main solutions: either the ejecta is highly constrained, or the underlying nuclear microphysics model is incomplete. We construct the light curves for these two cases in hopes of using kilonova to break this degeneracy. We find that our merger models produce light curves that are observationally too alike to robustly distinguish between the two valid models, offering little hope for constraining kilonova ejecta or the behavior of neutron-rich nuclei.

The Progenitor and Central Engine of a Peculiar GRB 230307A

Recently, a lack of supernova-associated with long-duration gamma-ray burst (GRB 230307A) at such a low redshift z = 0.065, but associated with a possible kilonova emission, has attracted great attention. Its heavy element nucleosynthesis and the characteristic of soft X-ray emission suggest that the central engine of GRB 230307A is a magnetar that is originated from a binary compact star merger. The calculated lower value of ε ∼ 0.05 suggests that GRB 230307A seems to have an ambiguous progenitor. The lower value of f eff = 1.23 implies that GRB 230307A is not likely to be from the effect of “tip of iceberg.” We adopt the magnetar central engine model to fit the observed soft X-ray emission with varying efficiency and find that the parameter constraints of the magnetar falls into a reasonable range, i.e., B < 9.4 × 1015 G and P < 2.5 ms for Γsat = 103, and B < 3.6 × 1015 G and P < 1.05 ms for Γsat = 104. Whether the progenitor of GBR 230307A is from the mergers of neutron star–white dwarf (NS–WD) or neutron star–neutron star (NS–NS) remains unknown. The difference of GW radiation between NS–NS merger and NS–WD merger may be a probe to distinguish the progenitor of GRB 230307A-like events in the future.

Heavy-element production in a compact object merger observed by JWST.

The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs)1, sources of high-frequency gravitational waves (GWs)2 and likely production sites for heavy-element nucleosynthesis by means of rapid neutron capture (the r-process)3. Here we present observations of the exceptionally bright GRB 230307A. We show that GRB 230307A belongs to the class of long-duration GRBs associated with compact object mergers4-6 and contains a kilonova similar to AT2017gfo, associated with the GW merger GW170817 (refs. 7-12). We obtained James Webb Space Telescope (JWST) mid-infrared imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns, which we interpret as tellurium (atomic mass A = 130) and a very red source, emitting most of its light in the mid-infrared owing to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r-process elements across a broad atomic mass range and play a central role in heavy-element nucleosynthesis across the Universe.

Three-dimensional General-relativistic Simulations of Neutrino-driven Winds from Magnetized Proto–Neutron Stars

Formed in the aftermath of a core-collapse supernova or neutron star merger, a hot proto–neutron star (PNS) launches an outflow driven by neutrino heating lasting for up to tens of seconds. Though such winds are considered potential sites for the nucleosynthesis of heavy elements via the rapid neutron capture process (r-process), previous work has shown that unmagnetized PNS winds fail to achieve the necessary combination of high entropy and/or short dynamical timescale in the seed nucleus formation region. We present three-dimensional general-relativistic magnetohydrodynamical simulations of PNS winds which include the effects of a dynamically strong (B ≳ 1015 G) dipole magnetic field. After initializing the magnetic field, the wind quickly develops a helmet-streamer configuration, characterized by outflows along open polar magnetic field lines and a “closed” zone of trapped plasma at lower latitudes. Neutrino heating within the closed zone causes the thermal pressure of the trapped material to rise in time compared to the polar outflow regions, ultimately leading to the expulsion of this matter from the closed zone on a timescale of ∼60 ms, consistent with the predictions of Thompson. The high entropies of these transient ejecta are still growing at the end of our simulations and are sufficient to enable a successful second-peak r-process in at least a modest ≳1% of the equatorial wind ejecta.

Cross section measurement of the Sm144(α,n)Gd147 reaction for studying the α -nucleus optical potential at astrophysical energies

Nuclear reactions involving alpha particles play an important role in various astrophysical processes such as the gamma-process of heavy element nucleosynthesis. The poorly known low-energy alpha-nucleus optical (AOMP) potential is a key parameter to estimate the rates of these reactions. The AOMP can be tested by measuring the cross section of alpha-scattering as well as alpha-induced reactions. Low energy elastic alpha-scattering on 144Sm has recently been measured with high precision. The aim of the present work was to complement that work by measuring the (a,n) cross sections on 144Sm at low energies. The experimental data shall be used to constrain the AOMP. From this potential the 144Sm(a,g)148Gd reaction rate can be derived with reduced uncertainties. The 144Sm(a,n)147Gd reaction was studied by bombarding Sm targets with alpha-beams provided by the cyclotron accelerator of Atomki. The cross section was determined using the activation method. The gamma-radiation following the decay of the 147Gd reaction product was measured with a HPGe detector. The experimental data are analyzed within the statistical model. The cross section was measured in the alpha-energy range between 13 and 20 MeV in 1 MeV steps, i.e., from close above the (a,n) threshold. The results were compared with statistical model calculations using various approaches and parametrizations for the AOMP, and excellent agreement was obtained for two recent potentials. However, these potentials cannot reproduce literature data for the 144Sm(a,g)148Gd reaction with the same accuracy. Constraints for the AOMP were derived from an analysis of the new 144Sm(a,n)147Gd data and literature data for 144Sm(a,g)148Gd. These constraints enable a determination of the reaction rate of the 144Sm(a,g)148Gd reaction with significantly reduced uncertainties of less than a factor of two.

On the impact of compact binary merger ejecta opacity on Kilonova transient signals

Returns of gravitational wave astronomy will largely benefit from the detection and identification of electromagnetic (EM) signatures to gravitationalwave sources. Kilonovae (KNe) are promising EM counterparts to compact binary mergers, offering to astronomers and nuclear astrophysicists a unique window to advance knowledge on the heavy-element nucleosynthesis and mergerdriven mass ejection. However, extremely heterogeneous post-merging ejecta composition of both light- and heavy-r process nuclei, implies strong effects on the KNe light-curve identification due to the varying opacity of the system. Hence, large uncertainties on the r-process final abundance via spectroscopic analysis of KNe signals are still present, hardly fixed by theoretical models. Here we will present some peculiar features of KNe studies, focusing on the opacity issue, from the atomic and plasma physics perspectives. In this view, efforts have been made recently at INFN-LNS, trying to put constraints on plasma opacity of interest for early-stage KNe emission. We will present the experimental progress on the problem, including instruments and methods which could open an interdisciplinary approach to tackle astrophysical problems in laboratory plasmas.

Measurement of the 160Gd(n, γ) cross section at n_TOF and its medical implications

Neutron-capture reactions on gadolinium isotopes play an important role in several fields of physics,in particular in nuclear Astrophysics for the understanding of the nucleosynthesis of heavy elements (beyond iron) in stars via the s- and r-processes [1] and in nuclear technology. Another important application of gadolinium is linked to the production of terbium, that offers a set of clinically interesting isotopes for theranostics, characterized by complementary physical decay characteristics. In particular, the low-energy β −emitter terbium-161 is very similar to lutetium-177 in terms of half-life (6.89 d), β − energy and chemical properties. Being a significant emitter of conversion/Auger electrons, greater therapeutic effect can therefore be expected in comparison to Lu-177 [2, 3]. For this reason, in the last decade, the study of the neutron capture reaction 160Gd(n,γ)161 Gd and the subsequent β − decay in terbium-161 is getting particular attention. As the nuclear data on the Gd-160 neutron capture reaction are quite scarce and inconsistent, a new measurement of the capture cross section of Gd-160 at the CERN neutron Time-Of-Flight facilty was performed in order to provide high resolution, high-accuracy data on this important reaction, in the energy range from thermal to hundreds of keV. In this contribution, the preliminary results of the n_TOF measurement are presented.

The ASTRAL database for neutron-capture nucleosynthesis studies

Present nuclear reaction network computations for astrophysical simulations involve many different types of rates, including neutron-capture reactions of interest for the modeling of heavy-element nucleosynthesis. While for many of them we still have to rely on theoretical calculations, an increasing number of experimentally-determined cross sections have now become available. In this contribution, we present “ASTrophysical Rate and rAw data Library” (ASTRAL), a new online database for neutron-capture cross sections based on experimental results, mainly obtained through activation and timeof-flight measurements. For the evaluation process, cross sections were re-calculated starting from raw data and by considering recent changes in physical properties of the involved isotopes (e.g., half-life and γ-ray intensities). We show the current status of the database, the techniques adopted to derive the recommended Maxwellian-averaged cross sections, and future developments.

Impact of Direct URCA Process on Heavy-Element Nucleosynthesis in X-ray Burst

Type-I X-ray bursts are the most powerful sites for the origin of heavy proton-rich nuclei, which are synthesized by the rapid-proton capture (rp) process. Qualitatively, heavy nuclei up to 107Te are synthe-sized due to the SnSbTe cycle, but how heavy nuclei the rp process can synthesize actually depends on many model parameters of X-ray bursts. In this paper, we focus on neutrino cooling inside NSs, which decreases the temperature of whole NSs. In particular, we examine the effects of nucleon direct Urca process and neu-tron superfluidity on final products in X-ray bursts. We find that most of the final products are insensitive to neutrino cooling properties, except for heavier nuclei more than 104Sn.

Model-independent determination of the dipole response of Zn66 using quasimonoenergetic and linearly polarized photon beams

[Background] Photon strength functions are an important ingredient in calculations relevant for the nucleosynthesis of heavy elements. The relation to the photoabsorption cross section allows to experimentally constrain photon strength functions by investigating the photo-response of atomic nuclei. [Purpose] We determine the photoresponse of $^{66}$Zn in the energy region of 5.6 MeV to 9.9 MeV and analyze the contribution of the "elastic" decay channel back to the ground state. In addition, for the elastic channel electric and magnetic dipole transitions were separated. [Methods] Nuclear resonance fluorescence experiments were performed using a linearly-polarized quasi-monoenergetic photon beam at the High Intensity $\gamma$-ray Source. Photon beam energies from 5.6 to 9.9 MeV with an energy spread of about 3% were selected in steps of 200-300 keV. Two High Purity Germanium detectors were used for the subsequent $\gamma$-ray spectroscopy. [Results] Full photoabsorption cross sections are extracted from the data making use of the monoenergetic character of the photon beam. For the ground-state decay channel, the average contribution of electric and magnetic dipole strengths is disentangled. The average branching ratio back to the ground state is determined as well. [Conclusions] The new results indicate lower cross sections when compared to the values extracted from a former experiment using bremsstrahlung on $^{66}$Zn. In the latter, the average branching ratio to the ground state is estimated from statistical-model calculations in order to analyze the data. Corresponding estimates from statistical-model calculations underestimate this branching ratio compared to the values extracted from the present analysis, which would partly explain the high cross sections determined from the bremsstrahlung data.

Weak process for Iron isotopes chain used as seed for heavy elements nucleosynthesis

The weak processes are a key to understand the nucleosynthesis of heavy elements in the universe. In this work we have studied weak processes such as β-decay and neutrino capture that are very important for the r-process. This process takes place in a region between neutron drip and β-stability line. For nuclei far from β-stability line the experimental data are very scarce. To study the nucleosynthesis of heavy elements in astrophysics sites like Binary Neutron Star merger (BNS) and Supernova explosions, we have employed the Gross Theory of Beta Decay (GTBD) to evaluate the neutrino-nucleus cross sections and β-decay half-lives. Our 56Fe(νe , e −)56Co cross section show a excelent agreement with microscopic models and the half-lives well reproduce the experimental data with a satisfactory accuracy. We have fitted the neutrino-nucleus cross sections obtained by GTBD using a four degree polynomial dependence of the incident neutrino energy. This dependence allows us to systematically calculate the neutrino-nucleus cross sections for a large amount of nuclei, as occurs in the r-process nucleosynthesis calculation (more than 4000 nuclei). Our fitting and interpolation procedures are shown to be in a good agreement and the GTBD cross sections are successfully reproduced.

On the synthesis of heavy nuclei in protomagnetar outflows and implications for ultra-high energy cosmic rays

ABSTRACT It has been suggested that strongly magnetized and rapidly rotating protoneutron stars (PNSs) may produce long duration gamma-ray bursts (GRBs) originating from stellar core collapse. We explore the steady-state properties and heavy element nucleosynthesis in neutrino-driven winds from such PNSs whose magnetic axis is generally misaligned with the axis of rotation. We consider a wide variety of central engine properties such as surface dipole field strength, initial rotation period, and magnetic obliquity to show that heavy element nuclei can be synthesized in the radially expanding wind. This process is facilitated provided the outflow is Poynting-flux dominated such that its low entropy and fast expansion time-scale enables heavy nuclei to form in a more efficient manner as compared to the equivalent thermal GRB outflows. We also examine the acceleration and survival of these heavy nuclei and show that they can reach sufficiently high energies ≳ 1020 eV within the same physical regions that are also responsible for powering gamma-ray emission, primarily through magnetic dissipation processes. Although these magnetized outflows generally fail to achieve the production of elements heavier than lanthanides for our explored electron fraction range 0.4–0.6, we show that they are more than capable of synthesizing nuclei near and beyond iron peak elements.

The R-process Alliance: A Nearly Complete R-process Abundance Template Derived from Ultraviolet Spectroscopy of the R-process-enhanced Metal-poor Star HD 222925* * Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555. These observations are primarily

We present a nearly complete rapid neutron-capture process (r-process) chemical inventory of the metal-poor ([Fe/H] = −1.46 ± 0.10) r-process-enhanced ([Eu/Fe] = +1.32 ± 0.08) halo star HD 222925. This abundance set is the most complete for any object beyond the solar system, with a total of 63 metals detected and seven with upper limits. It comprises 42 elements from 31 ≤ Z ≤ 90, including elements rarely detected in r-process-enhanced stars, such as Ga, Ge, As, Se, Cd, In, Sn, Sb, Te, W, Re, Os, Ir, Pt, and Au. We derive these abundances from an analysis of 404 absorption lines in ultraviolet spectra collected using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope and previously analyzed optical spectra. A series of appendices discusses the atomic data and quality of fits for these lines. The r-process elements from Ba to Pb, including all elements at the third r-process peak, exhibit remarkable agreement with the solar r-process residuals, with a standard deviation of the differences of only 0.08 dex (17%). In contrast, deviations among the lighter elements from Ga to Te span nearly 1.4 dex, and they show distinct trends from Ga to Se, Nb through Cd, and In through Te. The r-process contribution to Ga, Ge, and As is small, and Se is the lightest element whose production is dominated by the r-process. The lanthanide fraction, log X La = −1.39 ± 0.09, is typical for r-process-enhanced stars and higher than that of the kilonova from the GW170817 neutron-star merger event. We advocate adopting this pattern as an alternative to the solar r-process-element residuals when confronting future theoretical models of heavy-element nucleosynthesis with observations.

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 &amp;lt; 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.

The Type Icn SN 2021csp: Implications for the Origins of the Fastest Supernovae and the Fates of Wolf–Rayet Stars

Abstract We present observations of SN 2021csp, the second example of a newly identified type of supernova (SN) hallmarked by strong, narrow, P Cygni carbon features at early times (Type Icn). The SN appears as a fast and luminous blue transient at early times, reaching a peak absolute magnitude of −20 within 3 days due to strong interaction between fast SN ejecta (v ≈ 30,000 km s−1) and a massive, dense, fast-moving C/O wind shed by the WC-like progenitor months before explosion. The narrow-line features disappear from the spectrum 10–20 days after explosion and are replaced by a blue continuum dominated by broad Fe features, reminiscent of Type Ibn and IIn supernovae and indicative of weaker interaction with more extended H/He-poor material. The transient then abruptly fades ∼60 days post-explosion when interaction ceases. Deep limits at later phases suggest minimal heavy-element nucleosynthesis, a low ejecta mass, or both, and imply an origin distinct from that of classical Type Ic SNe. We place SN 2021csp in context with other fast-evolving interacting transients, and discuss various progenitor scenarios: an ultrastripped progenitor star, a pulsational pair-instability eruption, or a jet-driven fallback SN from a Wolf–Rayet (W-R) star. The fallback scenario would naturally explain the similarity between these events and radio-loud fast transients, and suggests a picture in which most stars massive enough to undergo a W-R phase collapse directly to black holes at the end of their lives.

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.