No Arabic abstract
Material ejected during (or immediately following) the merger of two neutron stars may assemble into heavy elements by the r-process. The subsequent radioactive decay of the nuclei can power electromagnetic emission similar to, but significantly dimmer than, an ordinary supernova. Identifying such events is an important goal of future transient surveys, offering new perspectives on the origin of r-process nuclei and the astrophysical sources of gravitational waves. Predictions of the transient light curves and spectra, however, have suffered from the uncertain optical properties of heavy ions. Here we consider the opacity of expanding r-process material and argue that it is dominated by line transitions from those ions with the most complex valence electron structure, namely the lanthanides. For a few representative ions, we run atomic structure models to calculate radiative data for tens of millions of lines. We find that the resulting r-process opacities are orders of magnitude larger than that of ordinary (e.g., iron-rich) supernova ejecta. Radiative transport calculations using these new opacities indicate that the transient emission should be dimmer and redder than previously thought. The spectra appear pseudo-blackbody, with broad absorption features, and peak in the infrared (~1 micron). We discuss uncertainties in the opacities and attempt to quantify their impact on the spectral predictions. The results have important implications for observational strategies to find and study the radioactively powered electromagnetic counterparts to compact object mergers.
We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron richness at the high temperatures predicted by the hydrodynamic simulation. Circumventing the complexities of modelling three-dimensional neutrino transport, (anti)neutrino captures are parameterized in terms of prescribed neutrino luminosities and mean energies, guided by published results and assumed as constant in time. Depending sensitively on the adopted neutrino-antineutrino luminosity ratio, neutrino processes increase Y_e to values between 0.25 and 0.40, still allowing for a successful r-process compatible with the observed solar abundance distribution and a significant fraction of the ejecta consisting of r-process nuclei. If the electron neutrino luminosities and mean energies are relatively large compared to the antineutrino properties, the mean Y_e might reach values >0.40 so that neutrino captures seriously compromise the success of the r-process. In this case, the r-abundances remain compatible with the solar distribution, but the total amount of ejected r-material is reduced to a few percent, because the production of iron-peak elements is favored. Proper neutrino physics, in particular also neutrino absorption, have to be included in BNSM simulations before final conclusions can be drawn concerning r-processing in this environment and concerning observational consequences like kilonovae, whose peak brightness and color temperature are sensitive to the composition-dependent opacity of the ejecta.
Star-to-star dispersion of r-process elements has been observed in a significant number of old, metal-poor globular clusters. We investigate early-time neutron-star mergers as the mechanism for this enrichment. Through both numerical modeling and analytical arguments, we show that neutron-star mergers cannot be induced through dynamical interactions early in the history of the cluster, even when the most liberal assumptions about neutron-star segregation are assumed. Therefore, if neutron-star mergers are the primary mechanism for r-process dispersion in globular clusters, they likely result from the evolution of isolated, primordial binaries in the clusters. Through population modeling, we find that moderate fractions of GCs with enrichment are only possible when a significant number of double neutron-star progenitors proceed through Case BB mass transfer --- under various assumptions for the initial properties of globular clusters, a neutron-star merger with the potential for enrichment will occur in ~15-60% (~30-90%) of globular clusters if this mass transfer proceeds stably (unstably). The strong anti-correlation between the pre-supernova orbital separation and post-supernova systemic velocity due to mass loss in the supernova leads to efficient ejection of most enrichment candidates from their host clusters. Thus, most enrichment events occur shortly after the double neutron stars are born. This requires star-forming gas that can absorb the r-process ejecta to be present in the globular cluster 30-50 Myr after the initial burst of star formation. If scenarios for redistributing gas in globular clusters cannot act on these timescales, the number of neutron-star merger enrichment candidates drops severely, and it is likely that another mechanism, such as r-process enrichment from collapsars, is at play.
The observation of a radioactively powered kilonova AT~2017gfo associated with the gravitational wave-event GW170817 from binary neutron star merger proves that these events are ideal sites for the production of heavy $r$-process elements. The gamma-ray photons produced by the radioactive decay of heavy elements are unique probes for the detailed nuclide compositions. Basing on the detailed $r$-process nucleosynthesis calculations and considering radiative transport calculations for the gamma-rays in different shells, we study the gamma-ray emission in a merger ejecta on a timescale of a few days. It is found that the total gamma-ray energy generation rate evolution is roughly depicted as $dot{E}propto t^{-1.3}$. For the dynamical ejecta with a low electron fraction ($Y_{rm e}lesssim0.20$), the dominant contributors of gamma-ray energy are the nuclides around the second $r$-process peak ($Asim130$), and the decay chain of $^{132}$Te ($t_{1/2}=3.21$~days) $rightarrow$ $^{132}$I ($t_{1/2}=0.10$~days) $rightarrow$ $^{132}$Xe produces gamma-ray lines at $228$ keV, $668$ keV, and $773$ keV. For the case of a wind ejecta with $Y_{rm e}gtrsim0.30$, the dominant contributors of gamma-ray energy are the nuclides around the first $r$-process peak ($Asim80$), and the decay chain of $^{72}$Zn ($t_{1/2}=1.93$~days) $rightarrow$ $^{72}$Ga ($t_{1/2}=0.59$~days) $rightarrow$ $^{72}$Ge produces gamma-ray lines at $145$ keV, $834$ keV, $2202$ keV, and $2508$ keV. The peak fluxes of these lines are $10^{-9}sim 10^{-7}$~ph~cm$^{-2}$ s$^{-1}$, which are marginally detectable with the next-generation MeV gamma-ray detector emph{ETCC} if the source is at a distance of $40$~Mpc.
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)
Heavy elements like gold, platinum or uranium are produced in the r-process, which needs neutron-rich and explosive environments. Neutron star mergers are a promising candidate for an r-process site. They exhibit three different channels for matter ejection fulfilling these conditions: dynamic ejecta due to tidal torques, neutrino-driven winds and evaporating matter from the accretion disk. We present a first study of the integrated nucleosynthesis for a neutrino-driven wind from a neutron star merger with a hyper-massive neutron star. Trajectories from a recent hydrodynamical simulation are divided into four different angle regions and post-processed with a reaction network. We find that the electron fraction varies around $Y_e approx 0.1 - 0.4$, but its distribution differs for every angle of ejection. Hence, the wind ejecta do not undergo a robust r-process, but rather possess distinct nucleosynthesis yields depending on the angle range. Compared to the dynamic ejecta, a smaller amount of neutron-rich matter gets unbound, but the production of lighter heavy elements with $A lesssim 130$ in the neutrino-driven wind can complement the strong r-process of the dynamic ejecta.