Nuclear-dominated accretion and subluminous supernovae from the merger of a white dwarf with a neutron star or black hole

We construct one dimensional steady-state models of accretion disks produced by the tidal disruption of a white dwarf (WD) by a neutron star (NS) or stellar mass black hole (BH). At radii r <~ 1e8.5-1e9 cm the midplane density and temperature are sufficiently high to burn the initial white dwarf material into increasingly heavier elements (e.g. Mg, Si, S, Ca, Fe, and Ni) at sequentially smaller radii. When the energy released by nuclear reactions is comparable to that released gravitationally, we term the disk a nuclear-dominated accretion flow (NuDAF). At small radii <~1e7 cm Fe photo-disintegrates into He and then free nuclei, and cooling by neutrinos may be efficient. At the high accretion rates of relevance ~ 0.1-1e-4 Msun/s, most of the disk is radiatively inefficient and prone to outflows powered by viscous dissipation and nuclear burning. Outflow properties are calculated by requiring that material in the midplane be marginally bound (Bernoulli constant <~ 0), due (in part) to cooling by matter escaping the disk. For reasonable assumptions regarding the properties of disk winds, we show that a significant fraction >50-80 per cent of the total WD mass is unbound. The ejecta composition is predominantly O, C, Si, Mg, Ne, Fe, and S [He, C, Si, S, Ar, and Fe], in the case of C-O [He] WDs, respectively, along with a small quantity ~1e-3-1e-2 Msun of radioactive Ni56 and, potentially, a trace amount of H. We use our results to evaluate possible EM counterparts of WD-NS/BH mergers, including optical transients powered by the radioactive decay of Ni56 and radio transients powered by the interaction of the ejecta with the interstellar medium. We address whether recently discovered subluminous Type I supernovae result from WD-NS/BH mergers. Our results also have implications for accretion following the core collapse of massive stars in collapsar models for gamma-ray bursts.Comment: 13 pages, 7 figures, 2 tables, now accepted to MNRA

High-Energy Neutrinos from Millisecond Magnetars formed from the Merger of Binary Neutron Stars

The merger of a neutron star (NS) binary may result in the formation of a long-lived, or indefinitely stable, millisecond magnetar remnant surrounded by a low-mass ejecta shell. A portion of the magnetar's prodigious rotational energy is deposited behind the ejecta in a pulsar wind nebula, powering luminous optical/X-ray emission for hours to days following the merger. Ions in the pulsar wind may also be accelerated to ultra-high energies, providing a coincident source of high energy cosmic rays and neutrinos. At early times, the cosmic rays experience strong synchrotron losses; however, after a day or so, pion production through photomeson interaction with thermal photons in the nebula comes to dominate, leading to efficient production of high-energy neutrinos. After roughly a week, the density of background photons decreases sufficiently for cosmic rays to escape the source without secondary production. These competing effects result in a neutrino light curve that peaks on a few day timescale near an energy of $\sim10^{18}$ eV. This signal may be detectable for individual mergers out to $\sim$ 10 (100) Mpc by current (next-generation) neutrino telescopes, providing clear evidence for a long-lived NS remnant, the presence of which may otherwise be challenging to identify from the gravitational waves alone. Under the optimistic assumption that a sizable fraction of NS mergers produce long-lived magnetars, the cumulative cosmological neutrino background is estimated to be $\sim 10^{-9}-10^{-8}\,\rm GeV\,cm^{-2}\,s^{-1}\,sr^{-1}$ for a NS merger rate of $10^{-7}\,\rm Mpc^{-3}\,yr^{-1}$, overlapping with IceCube's current sensitivity and within the reach of next-generation neutrino telescopes.Comment: 13 pages, 5 figures, accepted for publication in Ap

Shock-powered light curves of luminous red novae as signatures of pre-dynamical mass loss in stellar mergers

Luminous red novae (LRN) are a class of optical transients believed to originate from the mergers of binary stars, or "common envelope" events. Their light curves often show secondary maxima, which cannot be explained in the previous models of thermal energy diffusion or hydrogen recombination without invoking multiple independent shell ejections. We propose that double-peaked light curves are a natural consequence of a collision between dynamically-ejected fast shell and pre-existing equatorially-focused material, which was shed from the binary over many orbits preceding the dynamical event. The fast shell expands freely in the polar directions, powering the initial optical peak through cooling envelope emission. Radiative shocks from the collision in the equatorial plane power the secondary light curve peak on the radiative diffusion timescale of the deeper layers, similar to luminous Type IIn supernovae and some classical novae. Using a detailed 1D analytic model, informed by complementary 3D hydrodynamical simulations, we show that shock-powered emission can explain the observed range of peak timescales and luminosities of the secondary peaks in LRN for realistic variations in the binary parameters and fraction of the binary mass ejected. The dense shell created by the radiative shocks in the equatorial plane provides an ideal location for dust nucleation consistent with the the inferred aspherical geometry of dust in LRN. For giant stars, the ejecta forms dust when the shock-powered luminosity is still high, which could explain the infrared transients recently discovered by Spitzer. Our results suggest that pre-dynamical mass loss is common if not ubiquitous in stellar mergers, providing insight into the instabilities responsible for driving the binary merger.Comment: 12 pages, 6 figures, submitte

Three-dimensional GRMHD simulations of neutrino-cooled accretion disks from neutron star mergers

Merging binaries consisting of two neutron stars (NSs) or an NS and a stellar-mass black hole typically form a massive accretion torus around the remnant black hole or long-lived NS. Outflows from these neutrino-cooled accretion disks represent an important site for $r$-process nucleosynthesis and the generation of kilonovae. We present the first three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulations including weak interactions and a realistic equation of state of such accretion disks over viscous timescales ($380\,\mathrm{ms}$). We witness the emergence of steady-state MHD turbulence, a magnetic dynamo with an $\sim\!20\,\mathrm{ms}$ cycle, and the generation of a `hot' disk corona that launches powerful thermal outflows aided by the energy released as free nucleons recombine into $\alpha$-particles. We identify a self-regulation mechanism that keeps the midplane electron fraction low ($Y_e\sim0.1$) over viscous timescales. This neutron-rich reservoir, in turn, feeds outflows that retain a sufficiently low value of $Y_e\approx 0.2$ to robustly synthesize third-peak $r$-process elements. The quasi-spherical outflows are projected to unbind $40\%$ of the initial disk mass with typical asymptotic escape velocities of $0.1c$, and may thus represent the dominant mass ejection mechanism in NS-NS mergers. Including neutrino absorption, our findings agree with previous hydrodynamical $\alpha-$disk simulations that the entire range of $r$-process nuclei from the first to the third $r$-process peak can be synthesized in the outflows, in good agreement with observed solar system abundances. The asymptotic escape velocities and the quantity of ejecta, when extrapolated to moderately higher disk masses, are consistent with those needed to explain the red kilonova emission following the NS merger GW170817.Comment: 16 figures, 24 pages; matches published versio

Kilonovae

The mergers of double neutron star (NS-NS) and black hole (BH)-NS binaries are promising gravitational wave (GW) sources for Advanced LIGO and future GW detectors. The neutron-rich ejecta from such merger events undergoes rapid neutron capture (r-process) nucleosynthesis, enriching our Galaxy with rare heavy elements like gold and platinum. The radioactive decay of these unstable nuclei also powers a rapidly evolving, supernova-like transient known as a "kilonova". Kilonovae provide an approximately isotropic electromagnetic counterpart to the GW signal, which also provides a unique and direct probe of an important, if not dominant, r-process site. This handbook reviews the history and physics of kilonovae, leading to the current paradigm of week-long emission with a spectral peak at near-infrared wavelengths. Using a simple light curve model to illustrate the basic physics, I introduce potentially important variations on this canonical picture, including: ~day-long optical ("blue") emission from lanthanide-free components of the ejecta; ~hours-long precursor UV/blue emission, powered by the decay of free neutrons in the outermost ejecta layers; and enhanced emission due to energy input from a long-lived central engine, such as an accreting BH or millisecond magnetar. I assess the prospects of detecting kilonovae following future GW detections of NS-NS/BH-NS mergers in light of the recent follow-up campaign of the LIGO binary BH-BH mergers.Comment: published in Living Reviews in Relativit