We present a systematic numerical relativity study of the mass ejection and the associated electromagnetic transients and nucleosynthesis from binary neutron star (NS) mergers. We find that a few 10−3 M ⊙ of material is ejected dynamically during the mergers. The amount and the properties of these outflows depend on binary parameters and on the NS equation of state (EOS). A small fraction of these ejecta, typically ∼10−6 M ⊙, is accelerated by shocks formed shortly after merger to velocities larger than 0.6c and produces bright radio flares on timescales of weeks, months, or years after merger. Their observation could constrain the strength with which the NSs bounce after merger and, consequently, the EOS of matter at extreme densities. The dynamical ejecta robustly produce second and third r-process peak nuclei with relative isotopic abundances close to solar. The production of light r-process elements is instead sensitive to the binary mass ratio and the neutrino radiation treatment. Accretion disks of up to ∼0.2 M ⊙ are formed after merger, depending on the lifetime of the remnant. In most cases, neutrino- and viscously driven winds from these disks dominate the overall outflow. Finally, we generate synthetic kilonova light curves and find that kilonovae depend on the merger outcome and could be used to constrain the NS EOS.

We present fully general-relativistic simulations of binary neutron star mergers with a temperature and composition dependent nuclear equation of state. We study the dynamical mass ejection from both quasi-circular and dynamical-capture eccentric mergers. We systematically vary the level of our treatment of the microphysics to isolate the effects of neutrino cooling and heating and we compute the nucleosynthetic yields of the ejecta. We find that eccentric binaries can eject significantly more material than quasi-circular binaries and generate bright infrared and radio emission. In all our simulations the outflow is composed of a combination of tidally-and shock-driven ejecta, mostly distributed over a broad ∼ 60 • angle from the orbital plane, and, to a lesser extent, by thermally driven winds at high latitudes. Ejecta from eccentric mergers are typically more neutron rich than those of quasi-circular mergers. We find neutrino cooling and heating to affect, quantitatively and qualitatively, composition, morphology, and total mass of the outflows. This is also reflected in the infrared and radio signatures of the binary. The final nucleosynthetic yields of the ejecta are robust and insensitive to input physics or merger type in the regions of the second and third r-process peaks. The yields for elements on the first peak vary between our simulations, but none of our models is able to explain the Solar abundances of first-peak elements without invoking additional first-peak contributions from either neutrino and viscously-driven winds operating on longer timescales after the mergers, or from core-collapse supernovae.

The possibility that long tidal tails formed during compact object mergers may power optical transients through the decay of freshly synthesized r-process material is investigated. Precise modeling of the merger dynamics allows for a realistic determination of the thermodynamic conditions in the ejected debris. The results of hydrodynamic and full nuclear network calculations are combined to calculate the resultant r-process abundances and the heating of the material by their decays. The subsequent homologous structure is mapped into a radiative transfer code to synthesize emergent model light curves and determine how their properties (variability and color evolution) depend on the mass ratio and orientation of the merging binary. The radiation emanating from the ejected debris, though less spectacular than a typical supernova, should be observable in transient surveys and we estimate the associated detection rates. The case for (or against) compact object mergers as the progenitors of short gamma-ray bursts can be tested if such electromagnetic transients are detected (or not) in coincidence with some bursts, although they may be obscured by on-axis afterglows.

We present a first simulation of the post-merger evolution of a black hole-neutron star binary in full general relativity using an energy-integrated general relativistic truncated moment formalism for neutrino transport. We describe our implementation of the moment formalism and important tests of our code, before studying the formation phase of an accretion disk after a black hole-neutron star merger. We use as initial data an existing general relativistic simulation of the merger of a neutron star of mass 1.4M with a black hole of mass 7M and dimensionless spin χBH = 0.8. Comparing with a simpler leakage scheme for the treatment of the neutrinos, we find noticeable differences in the neutron to proton ratio in and around the disk, and in the neutrino luminosity. We find that the electron neutrino luminosity is much lower in the transport simulations, and that both the disk and the disk outflows are less neutron-rich. The spatial distribution of the neutrinos is significantly affected by relativistic effects, due to large velocities and curvature in the regions of strongest emission. Over the short timescale evolved, we do not observe purely neutrino-driven outflows. However, a small amount of material (3 × 10 −4 M ) is ejected in the polar region during the circularization of the disk. Most of that material is ejected early in the formation of the disk, and is fairly neutron rich (electron fraction Ye ∼ 0.15 − 0.25). Through r-process nucleosynthesis, that material should produce high-opacity lanthanides in the polar region, and could thus affect the lightcurve of radioactively powered electromagnetic transients. We also show that by the end of the simulation, while the bulk of the disk remains neutron-rich (Ye ∼ 0.15 − 0.2 and decreasing), its outer layers have a higher electron fraction: 10% of the remaining mass has Ye > 0.3. As that material would be the first to be unbound by disk outflows on longer timescales, and as composition evolution is slower at later times, the changes in Ye experienced during the formation phase of the disk could have an impact on nucleosynthesis outputs from neutrino-driven and viscously-driven outflows. Finally, we find that the effective viscosity due to momentum transport by neutrinos is unlikely to have a strong effect on the growth of the magnetorotational instability in the post-merger accretion disk.

Neutron star mergers are among the most promising sources of gravitational waves for advanced groundbased detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients powered by r-process nucleosynthesis in neutron-rich material ejected by the merger, and radio emission from the interaction of that ejecta with the interstellar medium. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the Spectral Einstein Code (SpEC) to simulate the merger of low-mass neutron star binaries (two 1.2M neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from recent simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived (i.e. rapid uniform rotation is sufficient to prevent its collapse). We observe strong excitations of l = 2, m = 2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes.

We present the SuperNova Explosion Code (SNEC), an open-source Lagrangian code for the hydrodynamics and equilibrium-diffusion radiation transport in the expanding envelopes of supernovae. Given a model of a progenitor star, an explosion energy, and an amount and distribution of radioactive nickel, SNEC generates the bolometric light curve, as well as the light curves in different broad bands assuming black body emission. As a first application of SNEC, we consider the explosions of a grid of 15 M (at zero-age main sequence) stars whose hydrogen envelopes are stripped to different extents and at different points in their evolution. The resulting light curves exhibit plateaus with durations of ∼20 − 100 days if 1.5 − 2 M of hydrogen-rich material is left and no plateau if less hydrogen-rich material is left. If these shorter plateau lengths are not seen for Type IIP supernovae in nature, it suggests that, at least for zero-age main sequence masses 20 M , hydrogen mass loss occurs as an all or nothing process. This perhaps points to the important role binary interactions play in generating the observed mass-stripped supernovae (i.e., Type Ib/c events). These light curves are also unlike what is typically seen for Type IIL supernovae, arguing that simply varying the amount of mass loss cannot explain these events. The most stripped models begin to show double-peaked light curves similar to what is often seen for Type IIb supernovae, confirming previous work that these supernovae can come from progenitors that have a small amount of hydrogen and a radius of ∼ 500 R .

Magnetohydrodynamic turbulence is important in many high-energy astrophysical systems, where instabilities can amplify the local magnetic field over very short timescales. Specifically, the magnetorotational instability and dynamo action have been suggested as a mechanism for the growth of magnetar-strength magnetic fields (of 10(15) gauss and above) and for powering the explosion of a rotating massive star. Such stars are candidate progenitors of type Ic-bl hypernovae, which make up all supernovae that are connected to long γ-ray bursts. The magnetorotational instability has been studied with local high-resolution shearing-box simulations in three dimensions, and with global two-dimensional simulations, but it is not known whether turbulence driven by this instability can result in the creation of a large-scale, ordered and dynamically relevant field. Here we report results from global, three-dimensional, general-relativistic magnetohydrodynamic turbulence simulations. We show that hydromagnetic turbulence in rapidly rotating protoneutron stars produces an inverse cascade of energy. We find a large-scale, ordered toroidal field that is consistent with the formation of bipolar magnetorotationally driven outflows. Our results demonstrate that rapidly rotating massive stars are plausible progenitors for both type Ic-bl supernovae and long γ-ray bursts, and provide a viable mechanism for the formation of magnetars. Moreover, our findings suggest that rapidly rotating massive stars might lie behind potentially magnetar-powered superluminous supernovae.

r-process nucleosynthesis in material ejected during neutron star mergers may lead to radioactively powered transients called kilonovae. The timescale and peak luminosity of these transients depend on the composition of the ejecta, which determines the local heating rate from nuclear decays and the opacity. Kasen et al. and Tanaka & Hotokezaka pointed out that lanthanides can drastically increase the opacity in these outflows. We use the new general-purpose nuclear reaction network SkyNet to carry out a parameter study of r-process nucleosynthesis for a range of initial electron fractions Y e , initial specific entropies s, and expansion timescales τ. We find that the ejecta is lanthanide-free for Y e  0.22−0.30, depending on s and τ. The heating rate is insensitive to s and τ, but certain, larger values of Y e lead to reduced heating rates, due to individual nuclides dominating the heating. We calculate approximate light curves with a simplified gray radiative transport scheme. The light curves peak at about a day (week) in the lanthanide-free (-rich) cases. The heating rate does not change much as the ejecta becomes lanthanide-free with increasing Y e , but the light-curve peak becomes about an order of magnitude brighter because it peaks much earlier when the heating rate is larger. We also provide parametric fits for the heating rates between 0.1 and 100 days, and we provide a simple fit in Y e , s, and τ to estimate whether or not the ejecta is lanthanide-rich.