R-process Nuclei Research Articles

ALMA and GMRT Constraints on the Off-axis Gamma-Ray Burst 170817A from the Binary Neutron Star Merger GW170817

Abstract Binary neutron-star mergers (BNSMs) are among the most readily detectable gravitational-wave (GW) sources with the Laser Interferometer Gravitational-wave Observatory (LIGO). They are also thought to produce short γ-ray bursts (SGRBs) and kilonovae that are powered by r-process nuclei. Detecting these phenomena simultaneously would provide an unprecedented view of the physics during and after the merger of two compact objects. Such a Rosetta Stone event was detected by LIGO/Virgo on 2017 August 17 at a distance of ∼44 Mpc. We monitored the position of the BNSM with Atacama Large Millimeter/submillimeter Array (ALMA) at 338.5 GHz and the Giant Metrewave Radio Telescope (GMRT) at 1.4 GHz, from 1.4 to 44 days after the merger. Our observations rule out any afterglow more luminous than in these bands, probing >2–4 dex fainter than previous SGRB limits. We match these limits, in conjunction with public data announcing the appearance of X-ray and radio emission in the weeks after the GW event, to templates of off-axis afterglows. Our broadband modeling suggests that GW170817 was accompanied by an SGRB and that the γ-ray burst (GRB) jet, powered by erg, had a half-opening angle of , and was misaligned by from our line of sight. The data are also consistent with a more collimated jet: erg, . This is the most conclusive detection of an off-axis GRB afterglow and the first associated with a BNSM-GW event to date. We use the viewing angle estimates to infer the initial bulk Lorentz factor and true energy release of the burst.

Kilonova from post-merger ejecta as an optical and near-Infrared counterpart of GW170817

Abstract Recent detection of gravitational waves from a neutron star (NS) merger event GW170817 and identification of an electromagnetic counterpart provide a unique opportunity to study the physical processes in NS mergers. To derive properties of ejected material from the NS merger, we perform radiative transfer simulations of kilonova, optical and near-infrared emissions powered by radioactive decays of r-process nuclei synthesized in the merger. We find that the observed near-infrared emission lasting for >10 d is explained by 0.03 M⊙ of ejecta containing lanthanide elements. However, the blue optical component observed at the initial phases requires an ejecta component with a relatively high electron fraction (Ye). We show that both optical and near-infrared emissions are simultaneously reproduced by the ejecta with a medium Ye of ∼0.25. We suggest that a dominant component powering the emission is post-merger ejecta, which exhibits that the mass ejection after the first dynamical ejection is quite efficient. Our results indicate that NS mergers synthesize a wide range of r-process elements and strengthen the hypothesis that NS mergers are the origin of r-process elements in the Universe.

The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. III. Optical and UV Spectra of a Blue Kilonova from Fast Polar Ejecta

Abstract We present optical and ultraviolet spectra of the first electromagnetic counterpart to a gravitational-wave (GW) source, the binary neutron star merger GW170817. Spectra were obtained nightly between 1.5 and 9.5 days post-merger, using the Southern Astrophysical Research and Magellan telescopes; the UV spectrum was obtained with the Hubble Space Telescope at 5.5 days. Our data reveal a rapidly fading blue component ( T ≈ 5500 K at 1.5 days) that quickly reddens; spectra later than ≳ 4.5 days peak beyond the optical regime. The spectra are mostly featureless, although we identify a possible weak emission line at ∼7900 Å at t ≲ 4.5 days. The colors, rapid evolution, and featureless spectrum are consistent with a “blue” kilonova from polar ejecta comprised mainly of light r-process nuclei with atomic mass number A ≲ 140 . This indicates a sightline within θ obs ≲ 45 ° of the orbital axis. Comparison to models suggests ∼0.03 M ⊙ of blue ejecta, with a velocity of ∼ 0.3 c . The required lanthanide fraction is ∼ 10 − 4 , but this drops to < 10 − 5 in the outermost ejecta. The large velocities point to a dynamical origin, rather than a disk wind, for this blue component, suggesting that both binary constituents are neutron stars (as opposed to a binary consisting of a neutron star and a black hole). For dynamical ejecta, the high mass favors a small neutron star radius of ≲ 12 km. This mass also supports the idea that neutron star mergers are a major contributor to r-process nucleosynthesis.

J-GEM observations of an electromagnetic counterpart to the neutron star merger GW170817

Abstract GW170817 is the first detected gravitational wave source from a neutron star merger. We present the Japanese collaboration for gravitational-wave electro-magnetic (J-GEM) follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5 mag decline in the z band during the period between 1.7 and 7.7 d after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder during later epochs: the z − H color has changed by approximately 2.5 mag during the period between 0.7 and 7.7 d. The rapid evolutions of both the color and the optical brightness are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized r-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that r-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with an ejecta mass of 0.01 M⊙, which suggests a more intense mass ejection (∼0.03 M⊙) or possibly an additional energy source.

Multi-messenger Observations of a Binary Neutron Star Merger*

Abstract On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ∼ 1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40 − 8 + 8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 M ⊙ . An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ∼ 40 Mpc ) less than 11 hours after the merger by the One-Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ∼10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ∼ 9 and ∼ 16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC 4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta.

Neutron Star Mergers and Nucleosynthesis of Heavy Elements

The existence of neutron star mergers has been supported since the discovery of the binary pulsar and the observation of its orbital energy loss, consistent with General Relativity. They are considered nucleosynthesis sites of the rapid neutron-capture process (r-process), which is responsible for creating approximately half of all heavy elements beyond Fe and is the only source of elements beyond Pb and Bi. Detailed nucleosynthesis calculations based on the decompression of neutron star matter are consistent with solar r-process abundances of heavy nuclei. Neutron star mergers have also been identified with short-duration [Formula: see text]-ray bursts via their IR afterglow. The high neutron densities in ejected matter permit a violent r-process, leading to fission cycling of the heaviest nuclei in regions far from (nuclear) stability. Uncertainties in several nuclear properties affect the abundance distributions. The modeling of astrophysical events also depends on the hydrodynamic treatment, the occurrence of a neutrino wind after the merger and before the possible emergence of a black hole, and the properties of black hole accretion disks. We discuss the effect of nuclear and modeling uncertainties and conclude that binary compact mergers are probably a (or the) dominant site of the production of r-process nuclei in our Galaxy.

Refining mass formulas for astrophysical applications: A Bayesian neural network approach

Exotic nuclei, particularly those near the driplines, are at the core of one of the fundamental questions driving nuclear structure and astrophysics today: what are the limits of nuclear binding? Exotic nuclei play a critical role in both informing theoretical models as well as in our understanding of the origin of the heavy elements. Our purpose is to refine existing mass models through the training of an artificial neural network that will mitigate the large model discrepancies far away from stability. The basic paradigm of our two-pronged approach is an existing mass model that captures as much as possible of the underlying physics followed by the implementation of a Bayesian Neural Network (BNN) refinement to account for the missing physics. Bayesian inference is employed to determine the parameters of the neural network so that model predictions may be accompanied by theoretical uncertainties. Despite the undeniable quality of the mass models adopted in this work, we observe a significant improvement (of about 40%) after the BNN refinement is implemented. Indeed, in the specific case of the Duflo-Zuker mass formula, we find that the rms deviation relative to experiment is reduced from rms =0.503MeV to rms=0.286 MeV. These newly refined mass tables are used to map the neutron drip lines (or rather "drip bands") and to study a few critical r-process nuclei. The BNN approach is highly successful in refining the predictions of existing mass models. In particular, the large discrepancy displayed by the original "bare" models in regions where experimental data is unavailable is considerably quenched after the BNN refinement. This lends credence to our approach and has motivated us to publish refined mass tables that we trust will be helpful for future astrophysical applications.

Empirical formula for β −-decay half-lives of r-process nuclei

Experimental data of $\beta^{-}$-decay half-lives of nuclei with atomic number between 20 and 190 are investigated. A systematic formula has been proposed to calculate $\beta^{-}$-decay half-lives of neutron-rich nuclei, with a particular consideration on shell and pair effects, the decay energy $Q$ as well as the nucleon numbers ($Z$, $N$). Although the formula has relatively few parameters, it reproduces the experimental $\beta^{-}$-decay half-lives of neutron-rich nuclei very well. The predicted half-lives for the r-process relevant nuclei obtained with the current formula serve as reliable input in the r-process model calculations.

Neutrino-heated winds from millisecond protomagnetars as sources of the weak r-process

We explore heavy element nucleosynthesis in neutrino-driven winds from rapidly-rotating, strongly magnetized proto-neutron stars for which the magnetic dipole is aligned with the rotation axis, and the field is assumed to be a static force-free configuration. We process the proto-magnetar wind trajectories calculated by Vlasov et al 2014 through the r-process nuclear reaction network SkyNet using contemporary models for the evolution of the wind electron fraction during the proto-neutron star cooling phase. Although we do not find a successful second or third peak r-process for any rotation period P, we show that proto-magnetars with P around 1-5 ms produce heavy element abundance distributions that extend to higher nuclear mass number than from otherwise equivalent spherical winds (with the mass fractions of some elements enhanced by factors of 100-1000). The heaviest elements are synthesized by outflows emerging along flux tubes which graze the closed zone and pass near the equatorial plane outside the light cylinder. Due to dependence of the nucleosynthesis pattern on the magnetic field strength and rotation rate of the proto-neutron star, natural variations in these quantities between core collapse events could contribute to the observed diversity of the abundances of weak r-process nuclei in metal-poor stars. Further diversity, including possibly even a successful third-peak r-process, could be achieved for misaligned rotators with non-zero magnetic inclination with respect to the rotation axis. If proto-magnetars are central engines for GRBs, their relativistic jets should contain a high mass fraction of heavy nuclei of characteristic average mass number A of order 100, providing a possible source for ultra-high energy cosmic rays comprised of heavy nuclei with an energy spectrum that extends beyond the nominal GZK cut-off for protons or iron nuclei.

The Intermediate r-process in Core-collapse Supernovae Driven by the Magneto-rotational Instability

Abstract We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability (MRI). A series of axisymmetric magneto-hydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the MRI. Neutrino-heating dominated explosions, enhanced by magnetic fields, showed mildly neutron-rich ejecta producing nuclei up to (i.e., the weak r-process), while explosion models with stronger magnetic fields reproduce a solar-like r-process pattern. More commonly seen abundance patterns in our models are in between the weak and regular r-process, producing lighter and intermediate-mass nuclei. These intermediate r-processes exhibit a variety of abundance distributions, compatible with several abundance patterns in r-process-enhanced metal-poor stars. The amount of Eu ejecta in magnetically driven jets agrees with predicted values in the chemical evolution of early galaxies. In contrast, neutrino-heating dominated explosions have a significant amount of Fe ( ) and Zn, comparable to regular supernovae and hypernovae, respectively. These results indicate magneto-rotational supernovae can produce a wide range of heavy nuclei from iron-group to r-process elements, depending on the explosion dynamics.

Dynamical ejecta from precessing neutron star-black hole mergers with a hot, nuclear-theory based equation of state

Neutron star-black hole binaries are among the strongest sources of gravitational waves detectable by current observatories. They can also power bright electromagnetic signals (gamma-ray bursts, kilonovae), and may be a significant source of production of r-process nuclei. A misalignment of the black hole spin with respect to the orbital angular momentum leads to precession of that spin and of the orbital plane, and has a significant effect on the properties of the post-merger remnant and of the material ejected by the merger. We present a first set of simulations of precessing neutron star-black hole mergers using a hot, composition dependent, nuclear-theory based equation of state (DD2). We show that the mass of the remnant and of the dynamical ejecta are broadly consistent with the result of simulations using simpler equations of state, while differences arise when considering the dynamics of the merger and the velocity of the ejecta. We show that the latter can easily be understood from assumptions about the composition of low-density, cold material in the different equations of state, and propose an updated estimate for the ejecta velocity which takes those effects into account. We also present an updated mesh-refinement algorithm which allows us to improve the numerical resolution used to evolve neutron star-black hole mergers.

Chemical Evolution of 244Pu in the Solar Vicinity and Its Implications for the Properties of r-process Production

Abstract Meteoritic abundances of r-process elements are analyzed to deduce the history of chemical enrichment by the r-process, from the beginning of disk formation to the present time in the solar vicinity. Our analysis combines the abundance information from short-lived radioactive nuclei such as 244Pu with the abundance information from stable r-process nuclei such as Eu. These two types of nuclei can be associated with one r-process event and an accumulation of events until the formation of the solar system, respectively. With the help of the observed local star formation (SF) history, we deduce the chemical evolution of 244Pu and obtain three main results: (i) the last r-process event occurred 130–140 Myr before the formation of the solar system; (ii) the present-day low 244Pu abundance as measured in deep-sea reservoirs results from the low recent SF rate compared to ∼4.5−5 Gyr ago; and (iii) there were ∼15 r-process events in the solar vicinity from the formation of the Galaxy to the time of solar system’s formation and ∼30 r-process events to the present time. Then, adopting the hypothesis that a neutron star (NS) merger is the r-process production site, we find that the ejected r-process elements are extensively spread out and mixed with interstellar matter, with a mass of M ⊙, which is about 100 times larger than that for supernova ejecta. In addition, the event frequency of r-process production is estimated to be 1 per ~1400 core-collapse supernovae, which is identical to the frequency of NS mergers estimated from the analysis of stellar abundances.

Comparison of yields of neutron-rich nuclei in proton- and photon-induced U238 fission

A comparative study of fission of actinides specially $^{238}$U, by proton and bremsstrahlung photon is performed. Relative mass distribution of $^{238}$U fission fragments have been explored theoretically for both proton and photon induced fission. The integrated yield along with charge distribution of the products are calculated to find out the neutron richness in comparison to the nuclei produced by r-process in nucleosynthesis. Some r-process nuclei in intermediate mass range for symmetric fission mode are found to be produced almost two order of magnitude more for proton induced fission than photofission, although rest of the neutron rich nuclei in the asymmetric mode are produced in comparable proportion for both the processes.

The influence of neutrinos on r-process nucleosynthesis in the ejecta of black hole–neutron star mergers

During the merger of a black hole and a neutron star, baryonic mass can become unbound from the system. Because the ejected material is extremely neutron-rich, the r-process rapidly synthesizes heavy nuclides as the material expands and cools. In this work, we map general relativistic models of black hole-neutron star (BHNS) mergers into a Newtonian smoothed particle hydrodynamics (SPH) code and follow the evolution of the thermodynamics and morphology of the ejecta until the outflows become homologous. We investigate how the subsequent evolution depends on our mapping procedure and find that the results are robust. Using thermodynamic histories from the SPH particles, we then calculate the expected nucleosynthesis in these outflows while varying the level of neutrino irradiation coming from the postmerger accretion disk. We find that the ejected material robustly produces r-process nucleosynthesis even for unrealistically high neutrino luminosities, due to the rapid velocities of the outflow. Nonetheless, we find that neutrinos can have an impact on the detailed pat- tern of the r-process nucleosynthesis. Electron neutrinos are captured by neutrons to produce protons while neutron capture is occurring. The produced protons rapidly form low mass seed nuclei for the r-process. These low mass seeds are eventually incorporated into the first r-process peak at A~78, producing mainly Ge and Se. We consider the mechanism of this process in detail and discuss if it can impact galactic chemical evolution of the first peak r-process nuclei.

TRANSPORT AND MIXING OF r-PROCESS ELEMENTS IN NEUTRON STAR BINARY MERGER BLAST WAVES

ABSTRACT The r-process nuclei are robustly synthesized in the material ejected during neutron star binary mergers (NSBMs). If NSBMs are indeed solely responsible for the solar system r-process abundances, a galaxy like our own would be required to host a few NSBMs per million years, with each event ejecting, on average, about 5 × 10−2 M ⊙ of r-process material. Because the ejecta velocities in the tidal tail are significantly larger than those in ordinary supernovae, NSBMs deposit a comparable amount of energy into the ISM. In contrast to extensive efforts studying spherical models for supernova remnant evolution, calculations quantifying the impact of NSBM ejecta in the ISM have been lacking. To better understand their evolution, we perform a suite of three-dimensional hydrodynamic simulations of isolated NSBM ejecta expanding in environments with conditions adopted from Milky-Way-like galaxy simulations. Although the remnant morphology is highly complex at early times, the subsequent radiative evolution is remarkably similar to that of a standard supernova. This implies that sub-resolution supernova feedback models can be used in galaxy-scale simulations that are unable to resolve the key evolutionary phases of NSBMs. Among other quantities, we examine the radius, mass, and kinetic energy content of the remnant at shell formation. We find that the shell formation epoch is attained when the swept-up mass is about 103(n H/1 cm−3)−2/7 M ⊙; at this point, the mass fraction of r-process material is enhanced up to two orders of magnitude in relation to a solar metallicity ISM.

On the robustness of the r-process in neutron-star mergers against variations of nuclear masses

r-process calculations have been performed for matter ejected dynamically in neutron star mergers (NSM), such calculations are based on a complete set of trajectories from a three-dimensional relativistic smoothed particle hydrodynamic (SPH) simulation. Our calculations consider an extended nuclear reaction network, including spontaneous, β- and neutron-induced fission and adopting fission yield distributions from the ABLA code. In this contribution we have studied the sensitivity of the r-process abundances to nuclear masses by using diferent mass models for the calculation of neutron capture cross sections via the statistical model. Most of the trajectories, corresponding to 90% of the ejected mass, follow a relatively slow expansion allowing for all neutrons to be captured. The resulting abundances are very similar to each other and reproduce the general features of the observed r-process abundance (the second and third peaks, the rare-earth peak and the lead peak) for all mass models as they are mainly determined by the fission yields. We find distinct differences in the predictions of the mass models at and just above the third peak, which can be traced back to different predictions of neutron separation energies for r-process nuclei around neutron number N = 130.

MODELS OF KILONOVA/MACRONOVA EMISSION FROM BLACK HOLE–NEUTRON STAR MERGERS

ABSTRACT Black hole–neutron star (BH–NS) mergers are among the most promising gravitational-wave sources for ground-based detectors, and gravitational waves from BH–NS mergers are expected to be detected in the next few years. The simultaneous detection of electromagnetic counterparts with gravitational waves would provide rich information about merger events. Among the possible electromagnetic counterparts from BH–NS mergers is the so-called kilonova/macronova, emission powered by the decay of radioactive r-process nuclei, which is one of the best targets for follow-up observations. We derive fitting formulas for the mass and the velocity of ejecta from a generic BH–NS merger based on recently performed numerical-relativity simulations. We combine these fitting formulas with a new semi-analytic model for a BH–NS kilonova/macronova lightcurve, which reproduces the results of radiation-transfer simulations. Specifically, the semi-analytic model reproduces the results of each band magnitude obtained by the previous radiation-transfer simulations within ∼1 mag. By using this semi-analytic model we found that, at 400 Mpc, the kilonova/macronova is as bright as 22–24 mag for cases with a small chirp mass and a high black hole spin, and >28 mag for a large chirp mass and a low black hole spin. We also apply our model to GRB 130603B as an illustration, and show that a BH–NS merger with a rapidly spinning black hole and a large neutron star radius is favored.

Impact of low-energy multipole excitations and pygmy resonances on radiative nucleon captures

Nuclear structure theory is considered in the framework of the development of a microscopic model for nucleon-capture astrophysical implementations. In particular, microscopically obtained strength functions from a theoretical method incorporating density functional theory and quasiparticle-phonon model are used as an input in a statistical reaction model. The approach is applied in systematic investigations of the impact of low-energy multipole excitations and pygmy resonances on dipole photoabsorption and radiative neutronand proton-capture cross sections of key s- and r-process nuclei which is discussed in comparison with the experiment. For the cases of the short-lived isotopes 89 Zr and 91 Mo theoretical predictions are made.

Nuclear robustness of therprocess in neutron-star mergers

We have performed r-process calculations for matter ejected dynamically in neutron star mergers based on a complete set of trajectories from a three-dimensional relativistic smoothed particle hydrodynamic simulation. Our calculations consider an extended nuclear network, including spontaneous, $\beta$- and neutron-induced fission and adopting fission yield distributions from the ABLA code. We have studied the sensitivity of the r-process abundances to nuclear masses by using different models. Most of the trajectories, corresponding to 90% of the ejected mass, follow a relatively slow expansion allowing for all neutrons to be captured. The resulting abundances are very similar to each other and reproduce the general features of the observed r-process abundance (the second and third peaks, the rare-earth peak and the lead peak) for all mass models as they are mainly determined by the fission yields. We find distinct differences in the abundance yields at and just above the third peak, which can be traced back to different predictions of neutron separation energies for r-process nuclei around neutron number $N=130$. The remaining trajectories, which contribute 10% by mass to the total integrated abundances, follow such a fast expansion that the r process does not use all the neutrons. This also leads to a larger variation of abundances among trajectories as fission does not dominate the r-process dynamics. The total integrated abundances are dominated by contributions from the slow abundances and hence reproduce the general features of the observed r-process abundances. We find that at timescales of weeks relevant for kilonova light curve calculations, the abundance of actinides is larger than the one of lanthanides. Hence actinides can be even more important than lanthanides to determine the photon opacities under kilonova conditions. (Abridged)

NEUTRINO-DRIVEN WINDS IN THE AFTERMATH OF A NEUTRON STAR MERGER: NUCLEOSYNTHESIS AND ELECTROMAGNETIC TRANSIENTS

We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of a binary neutron star merger. Our focus is the initial remnant phase when a massive central neutron star is present. Using tracers from a recent hydrodynamical simulation, we determine total masses and integrated abundances to characterize the composition of unbound matter. We find that the nucleosynthetic yields depend sensitively on both the life time of the massive neutron star and the polar angle. Matter in excess of up to $9 \cdot 10^{-3} M_\odot$ becomes unbound until $\sim 200~{\rm ms}$. Due to electron fractions of $Y_{\rm e} \approx 0.2 - 0.4$ mainly nuclei with mass numbers $A < 130$ are synthesized, complementing the yields from the earlier dynamic ejecta. Mixing scenarios with these two types of ejecta can explain the abundance pattern in r-process enriched metal-poor stars. Additionally, we calculate heating rates for the decay of the freshly produced radioactive isotopes. The resulting light curve peaks in the blue band after about $4~{\rm h}$. Furthermore, high opacities due to heavy r-process nuclei in the dynamic ejecta lead to a second peak in the infrared after $3-4~{\rm d}$.