Ejecta Of Neutron Star Mergers Research Articles

Signatures of r-process elements in kilonova spectra

AbstractBinary neutron star (NS) mergers have been expected to synthesize r-process elements and cause electromagnetic radiation called kilonovae. Although r-process nucleosynthesis was confirmed by the observations of GW170817/AT2017gfo, individual elements have not been identified except for strontium. Toward identification of elements in kilonova spectra, we perform radiative transfer simulations in NS merger ejecta. We find that Sr II triplet lines appear in the spectrum, which is consistent with the absorption feature observed in GW170817/AT2017gfo. The synthetic spectrum also shows the strong Ca II triplet lines. Absence of the Ca II line features in GW170817/AT2017gfo implies that the Ca/Sr ratio is < 0.002 in mass fraction, which is consistent with nucleosynthesis for electron fraction ≥ 0.40 and entropy per nucleon (in units of Boltzmann constant) ≥ 25. Identification of absorption lines in near-infrared wavelengths which have not yet been decoded may lead to clarify the abundances synthesized in NS merger ejecta.

Modeling kilonova emission from neutron star mergers

AbstractCoalescence of binary neutron stars gives rise to kilonova, thermal emission powered by radioactive decays of newly synthesized r-process nuclei. Observational properties of kilonova are largely affected by bound-bound opacities of r-process elements. It is, thus, important to understand atomic properties of heavy elements to link the observed signals with nucleosynthesis of neutron star mergers. In this paper, we introduce the latest status of kilonova modeling by focusing on the aspects of atomic physics. We perform systematic atomic structure calculations of r-process elements to understand element-to-element variation in the opacities. We demonstrate that the properties of the atomic structure of heavy elements are imprinted in the opacities of the neutron star merger ejecta and consequently in the kilonova light curves and spectra. Using this latest opacity dataset, we briefly discuss implications for GW170817, expected diversity of kilonova emission, and prospects for element identification in kilonova spectra.

Radioactive Heating Rate of r-process Elements and Macronova Light Curve

Abstract We study the heating rate of r-process nuclei and thermalization of decay products in neutron star merger ejecta and macronova (kilonova) light curves. Thermalization of charged decay products, i.e., electrons, α-particles, and fission fragments, is calculated according to their injection energy. The γ-ray thermalization processes are also properly calculated by taking the γ-ray spectrum of each decay into account. We show that the β-decay heating rate at later times approaches a power-law decline as ∝t −2.8, which agrees with the result of Waxman et al. We present a new analytic model to calculate macronova light curves, in which the density structure of the ejecta is accounted for. We demonstrate that the observed bolometric light curve and temperature evolution of the macronova associated with GW170817 are reproduced well by the β-decay heating rate with the solar r-process abundance pattern. We interpret the break in the observed bolometric light curve around a week as a result of the diffusion wave crossing a significant part of the ejecta rather than a thermalization break. We also show that the time-weighted integral of the bolometric light curve (Katz integral) is useful to provide an estimate of the total r-process mass from the observed data, which is independent of the highly uncertain radiative transfer. For the macronova in GW170817, the ejecta mass is robustly estimated as ≈0.05 M ⊙ for A min ≤ 72 and 85 ≤ A min ≤ 130 with the solar r-process abundance pattern. The code for computation of the heating rate and light curve for given initial nuclear abundances is publicly available.

A line-binned treatment of opacities for the spectra and light curves from neutron star mergers

ABSTRACT The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy r-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound–bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion–opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and thus eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion–opacity formalism. The agreement between the line-binned and expansion–opacity results is explained as arising from the similarity in their opacities in the limit of low optical depth, where radiation transport is important in the ejecta. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semirelativistic approaches.

Physical Conditions for the r-process. I. Radioactive Energy Sources of Kilonovae

Abstract Radioactive energies from unstable nuclei made in the ejecta of neutron star mergers play principal roles in powering kilonovae. In previous studies, power-law-type heating rates (e.g., ) have frequently been used, which may be inadequate if the ejecta are dominated by nuclei other than the A ∼ 130 region. We consider, therefore, two reference abundance distributions that match the r-process residuals to the solar abundances for A ≥ 69 (light trans-iron plus r-process elements) and A ≥ 90 (r-process elements). Nucleosynthetic abundances are obtained by using free-expansion models with three parameters: expansion velocity, entropy, and electron fraction. Radioactive energies are calculated as an ensemble of weighted free-expansion models that reproduce the reference abundance patterns. The results are compared with the bolometric luminosity (> a few days since merger) of the kilonova associated with GW170817. We find that the former case (fitted for A ≥ 69) with an ejecta mass 0.06 M ⊙ reproduces the light curve remarkably well, including its steepening at ≳7 days, in which the mass of r-process elements is ≈0.01 M ⊙. Two β-decay chains are identified: 66Ni 66Cu 66Zn and 72Zn 72Ga 72Ge with similar halflives of parent isotopes (≈2 days), which leads to an exponential-like evolution of heating rates during 1–15 days. The light curve at late times (>40 days) is consistent with additional contributions from the spontaneous fission of 254Cf and a few Fm isotopes. If this is the case, the GW170817 event is best explained by the production of both light trans-iron and r-process elements that originate from dynamical ejecta and subsequent disk outflows from the neutron star merger.

Synchrotron Radiation from the Fast Tail of Dynamical Ejecta of Neutron Star Mergers

Abstract We find, using high-resolution numerical relativistic simulations, that the tail of the dynamical ejecta of neutron star mergers extends to mildly relativistic velocities faster than 0.7c. The kinetic energy of this fast tail is ∼1047–1049 erg, depending on the neutron star equation of state and on the binary masses. The synchrotron flare arising from the interaction of this fast tail with the surrounding interstellar medium (ISM) can power the observed nonthermal emission that followed GW170817, provided that the ISM density is , the two neutron stars had roughly equal masses and the neutron star equation of state is soft (small neutron star radii). One of the generic predictions of this scenario is that the cooling frequency crosses the X-ray band on a timescale of a few months to a year, leading to a cooling break in the X-ray light curve. While the recent observation of the superluminal motion resolved by very long baseline interferometry (VLBI) rules out the dynamical ejecta scenario, the model described in this paper is generic and can be applied for future neutron star merger events.

Shock Acceleration of Electrons and Synchrotron Emission from the Dynamical Ejecta of Neutron Star Mergers

Abstract Neutron star mergers (NSMs) eject energetic subrelativistic dynamical ejecta into circumbinary media. Analogous to supernovae and supernova remnants, the NSM dynamical ejecta are expected to produce nonthermal emission by electrons accelerated at a shock wave. In this paper, we present the expected radio and X-ray signals by this mechanism, taking into account nonlinear diffusive shock acceleration (DSA) and magnetic field amplification. We suggest that the NSM is unique as a DSA site, where the seed relativistic electrons are abundantly provided by the decays of r-process elements. The signal is predicted to peak at a few 100–1000 days after the merger, determined by the balance between the decrease of the number of seed electrons and the increase of the dissipated kinetic energy, due to the shock expansion. While the resulting flux can ideally reach the maximum flux expected from near-equipartition, the available kinetic energy dissipation rate of the NSM ejecta limits the detectability of such a signal. It is likely that the radio and X-ray emission are overwhelmed by other mechanisms (e.g., an off-axis jet) for an observer placed in a jet direction (i.e., for GW170817). However, for an off-axis observer, to be discovered once a number of NSMs are identified, the dynamical ejecta component is predicted to dominate the nonthermal emission. While the detection of this signal is challenging even with near-future facilities, this potentially provides a robust probe of the creation of r-process elements in NSMs.

Properties of Kilonovae from Dynamical and Post-merger Ejecta of Neutron Star Mergers

Abstract Ejected material from neutron star mergers gives rise to electromagnetic emission powered by radioactive decays of r-process nuclei, the so-called kilonova or macronova. While properties of the emission are largely affected by opacities in the ejected material, available atomic data for r-process elements are still limited. We perform atomic structure calculations for r-process elements: Se (Z = 34), Ru (Z = 44), Te (Z = 52), Ba (Z = 56), Nd (Z = 60), and Er (Z = 68). We confirm that the opacities from bound–bound transitions of open f-shell, lanthanide elements (Nd and Er) are higher than those of the other elements over a wide wavelength range. The opacities of open s-shell (Ba), p-shell (Se and Te), and d-shell (Ru) elements are lower than those of open f-shell elements, and their transitions are concentrated in the ultraviolet and optical wavelengths. We show that the optical brightness can be different by mag depending on the element abundances in the ejecta such that post-merger, lanthanide-free ejecta produce brighter and bluer optical emission. Such blue emission from post-merger ejecta can be observed from the polar directions if the mass of the preceding dynamical ejecta in these regions is small. For the ejecta mass of 0.01 , observed magnitudes of the blue emission will reach 21.0 mag (100 Mpc) and 22.5 mag (200 Mpc) in the g and r bands within a few days after the merger, which are detectable with 1 m or 2 m class telescopes.

Synthesis of Heavy Elements in the Ejecta of Neutron Star Mergers

We present a nucleosynthesis study on the neutrino-driven wind after neutron star merger. Post-processing tracers from a hydrodynamical simulation, we determine nucleosynthesis yields that depend on both life time of the massive neutron star and polar angle. In comparison with the dynamic ejecta we find a complementary nucleosynthesis, which can explain the r-process pattern beyond the first r-process peak. Additionally, we predict possible forms of light curves for the detection of a kilonova.

Dynamics, nucleosynthesis, and kilonova signature of black hole—neutron star merger ejecta

We investigate the ejecta from black hole—neutron star mergers by modeling the formation and interaction of mass ejected in a tidal tail and a disk wind. The outflows are neutron-rich, giving rise to optical/infrared emission powered by the radioactive decay of r-process elements (a kilonova). Here we perform an end-to-end study of this phenomenon, where we start from the output of a fully-relativistic merger simulation, calculate the post-merger hydrodynamical evolution of the ejecta and disk winds including neutrino physics, determine the final nucleosynthetic yields using post-processing nuclear reaction network calculations, and compute the kilonova emission with a radiative transfer code. We study the effects of the tail-to-disk mass ratio by scaling the tail density. A larger initial tail mass results in fallback matter becoming mixed into the disk and ejected in the subsequent disk wind. Relative to the case of a disk without dynamical ejecta, the combined outflow has lower mean electron fraction, faster speed, larger total mass, and larger absolute mass free of high-opacity Lanthanides or Actinides. In most cases, the nucleosynthetic yield is dominated by the heavy r-process contribution from the unbound part of the dynamical ejecta. A Solar-like abundance distribution can however be obtained when the total mass of the dynamical ejecta is comparable to the mass of the disk outflows. The kilonova has a characteristic duration of 1 week and a luminosity of ∼ erg , with orientation effects leading to variations of a factor ∼2 in brightness. At early times (1 d) the emission includes an optical component from the (hot) Lanthanide-rich material, but the spectrum evolves quickly to the infrared thereafter.

Analytic heating rate of neutron star merger ejecta derived from Fermi's theory of beta decay

Macronovae (kilonovae) that arise in binary neutron star mergers are powered by radioactive beta decay of hundreds of $r$-process nuclides. We derive, using Fermi's theory of beta decay, an analytic estimate of the nuclear heating rate. We show that the heating rate evolves as a power law ranging between $t^{-6/5}$ to $t^{-4/3}$. The overall magnitude of the heating rate is determined by the mean values of nuclear quantities, e.g., the nuclear matrix elements of beta decay. These values are specified by using nuclear experimental data. We discuss the role of higher order beta transitions and the robustness of the power law. The robust and simple form of the heating rate suggests that observations of the late-time bolometric light curve $\propto t^{-\frac{4}{3}}$ would be a direct evidence of a $r$-process driven macronova. Such observations could also enable us to estimate the total amount of $r$-process nuclei produced in the merger.

Radioactive decay products in neutron star merger ejecta: heating efficiency and γ-ray emission

The radioactive decay of the freshly synthesized $r$-process nuclei ejected in compact binary mergers power optical/infrared macronovae (kilonovae) that follow these events. The light curves depend critically on the energy partition among the different products of the radioactive decay and this plays an important role in estimates of the amount of ejected $r$-process elements from a given observed signal. We study the energy partition and $\gamma$-ray emission of the radioactive decay. We show that $20$-$50\%$ of the total radioactive energy is released in $\gamma$-rays on timescales from hours to a month. The number of emitted $\gamma$-rays per unit energy interval has roughly a flat spectrum between a few dozen keV and $1$ MeV so that most of this energy is carried by $\sim 1$ MeV $\gamma$-rays. However at the peak of macronova emission the optical depth of the $\gamma$-rays is $\sim 0.02$ and most of the $\gamma$-rays escape. The loss of these $\gamma$-rays reduces the heat deposition into the ejecta and hence reduces the expected macronova signals if those are lanthanides dominated. This implies that the ejected mass is larger by a factor of $2$-$3$ than what was previously estimated. Spontaneous fission heats up the ejecta and the heating rate can increase if a sufficient amount of transuranic nuclei are synthesized. Direct measurements of these escaping $\gamma$-rays may provide the ultimate proof for the macronova mechanisms and an identification of the $r$-process nucleosynthesis sites. However, the chances to detect these signals are slim with current X-ray and $\gamma$-ray missions. New detectors, more sensitive by at least a factor of ten, are needed for a realistic detection rate.

The long-term evolution of neutron star merger remnants - I. The impact of r-process nucleosynthesis

We follow the longterm evolution of the dynamic ejecta of neutron star mergers for up to 100 years and over a density range of roughly 40 orders of magnitude. We include the nuclear energy input from the freshly synthesized, radioactively decaying nuclei in our simulations and study its effects on the remnant dynamics. Although the nuclear heating substantially alters the longterm evolution, we find that running nuclear networks over purely hydrodynamic simulations (i.e. without heating) yields actually acceptable nucleosynthesis results. The main dynamic effect of the radioactive heating is to quickly smooth out inhomogeneities in the initial mass distribution, subsequently the evolution proceeds self-similarly and after 100 years the remnant still carries the memory of the initial binary mass ratio. We also explore the nucleosynthetic yields for two mass ejection channels. The dynamic ejecta very robustly produce "strong" r-process elements with $A > 130$ with a pattern that is essentially independent of the details of the merging system. From a simple model we find that neutrino-driven winds yield "weak" r-process contributions with $50 < A < 130$ whose abundance patterns vary substantially between different merger cases. This is because their electron fraction, set by the ratio of neutrino luminosities, varies considerably from case to case. Such winds do not produce any $^{56}{\rm Ni}$, but a range of radioactive isotopes that are long-lived enough to produce a second, radioactively powered electromagnetic transient in addition to the "macronova" from the dynamic ejecta. While our wind model is very simple, it nevertheless demonstrates the potential of such neutrino-driven winds for electromagnetic transients and it motivates further, more detailed neutrino-hydrodynamic studies. The properties of the mentioned transients are discussed in more detail in a companion paper.

The long-term evolution of neutron star merger remnants – II. Radioactively powered transients

We use 3D hydrodynamic simulations of the long-term evolution of neutron star merger ejecta to predict the light curves of electromagnetic transients that are powered by the decay of freshly produced r-process nuclei. For the dynamic ejecta that are launched by tidal and hydrodynamic interaction, we adopt grey opacities of 10 cm$^2$/g, as suggested by recent studies. For our reference case of a 1.3-1.4 $M_\odot$ merger, we find a broad IR peak 2-4 d after the merger. The peak luminosity is $\approx 2\times 10^{40}$ erg/s for an average orientation, but increased by up to a factor of 4 for more favourable binary parameters and viewing angles. These signals are rather weak and hardly detectable within the large error box (~100 deg$^2$) of a gravitational wave trigger. A second electromagnetic transient results from neutrino-driven winds. These winds produce `weak' r-process material with $50 < A < 130$ and abundance patterns that vary substantially between different merger cases. For an adopted opacity of 1 cm$^2$/g, the resulting transients peak in the UV/optical about 6 h after the merger with a luminosity of $\approx 10^{41}$ erg/s (for a wind of 0.01 $M_\odot$) These signals are marginally detectable in deep follow-up searches (e.g. using Hypersuprime camera on Subaru). A subsequent detection of the weaker but longer lasting IR signal would allow an identification of the merger event. We briefly discuss the implications of our results to the recent detection of an nIR transient accompanying GRB 130603B.

Electromagnetic counterparts of compact object mergers powered by the radioactive decay of r-process nuclei

The most promising astrophysical sources of kHz gravitational waves (GWs) are the inspiral and merger of binary neutron star(NS)/black hole systems. Maximizing the scientific return of a GW detection will require identifying a coincident electro-magnetic (EM) counterpart. One of the most likely sources of isotropic EM emission from compact object mergers is a supernova-like transient powered by the radioactive decay of heavy elements synthesized in ejecta from the merger. We present the first calculations of the optical transients from compact object mergers that self-consistently determine the radioactive heating by means of a nuclear reaction network; using this heating rate, we model the light curve with a one dimensional Monte Carlo radiation transfer calculation. For an ejecta mass ~1e-2 M_sun[1e-3 M_sun] the resulting light curve peaks on a timescale ~ 1 day at a V-band luminosity nu L_nu ~ 3e41[1e41] ergs/s (M_V = -15[-14]); this corresponds to an effective "f" parameter ~3e-6 in the Li-Paczynski toy model. We argue that these results are relatively insensitive to uncertainties in the relevant nuclear physics and to the precise early-time dynamics and ejecta composition. Due to the rapid evolution and low luminosity of NS merger transients, EM counterpart searches triggered by GW detections will require close collaboration between the GW and astronomical communities. NS merger transients may also be detectable following a short-duration Gamma-Ray Burst or "blindly" with present or upcoming optical transient surveys. Because the emission produced by NS merger ejecta is powered by the formation of rare r-process elements, current optical transient surveys can directly constrain the unknown origin of the heaviest elements in the Universe.

Nucleosynthesis calculations for the ejecta of neutron star coalescences

We present the results of fully dynamical r-process network calculations for the ejecta of neutron star mergers (NSMs). The late stages of the inspiral and the final violent coalescence of a neutron star binary have been calculated in detail using a 3D hydrodynamics code (Newtonian gravity plus backreaction forces emerging from the emission of gravitational waves) and a realistic nuclear equation of state. The found trajectories for the ejecta serve as input for dynamical r-process calculations where all relevant nuclear reactions (including beta-decays depositing nuclear energy in the expanding material) are followed. We find that all the ejected material undergoes r-process. For an initial Ye close to 0.1 the abundance distributions reproduce very accurately the solar r-process pattern for nuclei with A above 130. For lighter nuclei strongly underabundant (as compared to solar) distributions are encountered. We show that this behaviour is consistent with the latest observations of very old, metal-poor stars, despite simplistic arguments that have recently been raised against the possibility of NSM as possible sources of Galactic r-process material.