The final inspiral of double neutron star and neutron-star-black-hole binaries are likely to be detected by advanced networks of ground-based gravitational wave (GW) interferometers. Maximizing the ...science returns from such a discovery will require the identification of an electromagnetic counterpart. Here we critically evaluate and compare several possible counterparts, including short-duration gamma-ray bursts (SGRBs), "orphan" optical and radio afterglows, and day-long optical transients powered by the radioactive decay of heavy nuclei synthesized in the merger ejecta ("kilonovae"). We assess the promise of each counterpart in terms of four "Cardinal Virtues": detectability, high fraction, identifiability, and positional accuracy. Taking into account the search strategy for typical error regions of tens of square degrees, we conclude that SGRBs are the most useful to confirm the cosmic origin of a few GW events, and to test the association with neutron star mergers. However, for the more ambitious goal of localizing and obtaining redshifts for a large sample of GW events, kilonovae are instead preferred. Off-axis optical afterglows are detectable for at most tens of percent of events, while radio afterglows are promising only for energetic relativistic ejecta in a high-density medium. Our main recommendations are: (1) an all-sky gamma-ray satellite is essential for temporal coincidence detections, and for GW searches of gamma-ray-triggered events; (2) the Large Synoptic Survey Telescope should adopt a one-day cadence follow-up strategy, ideally with 0.5 hr per pointing to cover GW error regions; and (3) radio searches should focus on the relativistic case, which requires observations for a few months.
The kilonova emission observed following the binary neutron star merger event GW170817 provided the first direct evidence for the synthesis of heavy nuclei through the rapid neutron capture process ...(r process). The late-time transition in the spectral energy distribution to near-infrared wavelengths was interpreted as indicating the production of lanthanide nuclei, with atomic mass number A≳140. However, compelling evidence for the presence of even heavier third-peak (A≈195) r-process elements (e.g., gold, platinum) or translead nuclei remains elusive. At early times (∼days) most of the r-process heating arises from a large statistical ensemble of β decays, which thermalize efficiently while the ejecta is still dense, generating a heating rate that is reasonably approximated by a single power law. However, at later times of weeks to months, the decay energy input can also possibly be dominated by a discrete number of α decays, ^{223}Ra (half-life t_{1/2}=11.43 d), ^{225}Ac (t_{1/2}=10.0 d, following the β decay of ^{225}Ra with t_{1/2}=14.9 d), and the fissioning isotope ^{254}Cf (t_{1/2}=60.5 d), which liberate more energy per decay and thermalize with greater efficiency than β-decay products. Late-time nebular observations of kilonovae which constrain the radioactive power provide the potential to identify signatures of these individual isotopes, thus confirming the production of heavy nuclei. In order to constrain the bolometric light to the required accuracy, multiepoch and wideband observations are required with sensitive instruments like the James Webb Space Telescope. In addition, by comparing the nuclear heating rate obtained with an abundance distribution that follows the solar r abundance pattern, to the bolometric lightcurve of AT2017gfo, we find that the yet-uncertain r abundance of ^{72}Ge plays a decisive role in powering the lightcurve, if one assumes that GW170817 has produced a full range of the solar r abundances down to mass number A∼70.
We construct one-dimensional steady-state models of accretion discs produced by the tidal disruption of a white dwarf (WD) by a neutron star (NS) or stellar mass black hole (BH). At radii r ≲ ...108.5-109 cm the mid-plane density and temperature are sufficiently high to burn the initial WD 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 disc a nuclear-dominated accretion flow (NuDAF). At small radii ≲107 cm iron photodisintegrates into helium and then free nuclei, and in the very innermost disc cooling by neutrinos may be efficient. At the high accretion rates of relevance ∼10−4 to 0.1 M⊙ s−1, most of the disc is radiatively inefficient and prone to outflows powered by viscous dissipation and nuclear burning. Outflow properties are calculated by requiring that material in the mid-plane be marginally bound (Bernoulli constant ≲ 0), due (in part) to cooling by matter escaping the disc. For reasonable assumptions regarding the properties of disc winds, we show that a significant fraction (≳ 50-80 per cent) of the total WD mass is unbound. The composition of the ejecta is predominantly O, C, Si, Mg, Ne, Fe and S He, C, Si, S, Ar and Fe, in the case of C-O pure He WDs, respectively, along with a small quantity ∼10−3 to 10−2 M⊙ of radioactive 56Ni and, potentially, a trace amount of hydrogen. Depending on the pressure dependence of wind cooling, we find that the disc may be thermally unstable to nuclear burning, the likelihood of which increases for higher mass WDs. We use our results to evaluate possible electromagnetic counterparts of WD-NS/BH mergers, including optical transients powered by the radioactive decay of 56Ni 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. Ultimately assessing the fate of these events requires global simulations of the disc evolution, which capture the complex interplay between nuclear burning, convection and outflows.
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.
ABSTRACT
Mergers of compact binaries containing two neutron stars (NS–NS), or a neutron star and a stellar mass black hole (NS–BH), are likely progenitors of short-duration gamma-ray bursts (SGRBs). ...A fraction ${\gtrsim } 20{{\ \rm per\ cent}}$ of SGRBs is followed by temporally extended (≳minute-long), variable X-ray emission, attributed to ongoing activity of the central engine. One source of late-time engine activity is fall-back accretion of bound tidal ejecta; however, observed extended emission light curves do not track the naively anticipated, uninterrupted t−5/3 power-law decay, instead showing a lull or gap in emission typically lasting tens of seconds after the burst. Here, we re-examine the impact of heating due to rapid neutron capture (r-process) nucleosynthesis on the rate of the fall-back accretion, using ejecta properties extracted from numerical relativity simulations of NS–BH mergers. Heating by the r-process has its greatest impact on marginally bound matter, hence its relevance to late-time fall-back. Depending on the electron fraction of the ejecta and the mass of the remnant black hole, r-process heating can imprint a range of fall-back behaviour, ranging from temporal gaps of up to tens of seconds to complete late-time cut-off in the accretion rate. This behaviour is robust to realistic variations in the nuclear heating experienced by different parts of the ejecta. Central black holes with masses ${\lesssim } 3\, \mathrm{M}_{\odot }$ typically experience absolute cut-offs in the fall-back rate, while more massive ${\gtrsim } 6\!-\!8\, \mathrm{M}_{\odot }$ black holes instead show temporal gaps. We thus propose that SGRBs showing extended X-ray emission arise from NS–BH, rather than NS–NS, mergers. Our model implies an NS–BH merger detection rate by LIGO that, in steady state, is comparable to or greater than that of NS–NS mergers.
Supernovae are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their ...explosion. These are thought to appear as type Ib/c and type II supernovae, and are associated with young stellar populations. In contrast, the thermonuclear detonation of a carbon-oxygen white dwarf, whose mass approaches the Chandrasekhar limit, is thought to produce type Ia supernovae. Such supernovae are observed in both young and old stellar environments. Here we report a faint type Ib supernova, SN 2005E, in the halo of the nearby isolated galaxy, NGC 1032. The ‘old’ environment near the supernova location, and the very low derived ejected mass (∼0.3 solar masses), argue strongly against a core-collapse origin. Spectroscopic observations and analysis reveal high ejecta velocities, dominated by helium-burning products, probably excluding this as a subluminous or a regular type Ia supernova. We conclude that it arises from a low-mass, old progenitor, likely to have been a helium-accreting white dwarf in a binary. The ejecta contain more calcium than observed in other types of supernovae and probably large amounts of radioactive 44Ti.
We report the discovery of rising X-ray emission from the binary neutron star merger event GW170817. This is the first detection of X-ray emission from a gravitational-wave (GW) source. Observations ...acquired with the Chandra X-ray Observatory (CXO) at t 2.3 days post-merger reveal no significant emission, with L x 3.2 × 10 38 erg s − 1 (isotropic-equivalent). Continued monitoring revealed the presence of an X-ray source that brightened with time, reaching L x 9 × 10 38 erg s − 1 at 15.1 days post-merger. We interpret these findings in the context of isotropic and collimated relativistic outflows (both on- and off-axis). We find that the broadband X-ray to radio observations are consistent with emission from a relativistic jet with kinetic energy E k ∼ 10 49 − 50 erg , viewed off-axis with θ obs ∼ 20 ° - 40 ° . Our models favor a circumbinary density n ∼ 10 − 4 - 10 − 2 cm − 3 , depending on the value of the microphysical parameter ϵ B = 10 − 4 - 10 − 2 . A central-engine origin of the X-ray emission is unlikely. Future X-ray observations at t 100 days, when the target will be observable again with the CXO, will provide additional constraints to solve the model degeneracies and test our predictions. Our inferences on θ obs are testable with GW information on GW170817 from advanced LIGO/Virgo on the binary inclination.
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 electromagnetic (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 ∼10−2 M⊙ (10−3 M⊙) the resulting light-curve peaks on a time-scale ∼1 d at a V-band luminosity νLν∼ 3 × 1041 (1041) erg s−1MV=−15(−14); this corresponds to an effective ‘f’ parameter ∼3 × 10−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. Since NS merger transients peak at a luminosity that is a factor of ∼103 higher than a typical nova, we propose naming these events ‘kilo-novae’. Because of 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.
We present Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA) radio observations of GW170817, the first Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo ...gravitational wave (GW) event from a binary neutron star merger and the first GW event with an electromagnetic (EM) counterpart. Our data include the first observations following the discovery of the optical transient at both the centimeter (13.7 hr post-merger) and millimeter (2.41 days post-merger) bands. We detect faint emission at 6 GHz at 19.47 and 39.23 days after the merger, but not in an earlier observation at 2.46 days. We do not detect cm/mm emission at the position of the optical counterpart at frequencies of 10-97.5 GHz at times ranging from 0.6 to 30 days post-merger, ruling out an on-axis short gamma-ray burst (SGRB) for energies 10 48 erg. For fiducial SGRB parameters, our limits require an observer viewer angle of 20°. The radio and X-ray data can be jointly explained as the afterglow emission from an SGRB with a jet energy of ∼ 10 49 - 10 50 erg that exploded in a uniform density environment with n ∼ 10 − 4 - 10 − 2 cm−3, viewed at an angle of ∼20°-40° from the jet axis. Using the results of our light curve and spectral modeling, in conjunction with the inference of the circumbinary density, we predict the emergence of late-time radio emission from the deceleration of the kilonova (KN) ejecta on a timescale of ∼5-10 years that will remain detectable for decades with next-generation radio facilities, making GW170817 a compelling target for long-term radio monitoring.
We report deep Chandra X-ray Observatory (CXO), Hubble Space Telescope (HST), and Karl J. Jansky Very Large Array (VLA) observations of the binary neutron star event GW170817 at t < 160 days after ...merger. These observations show that GW170817 has been steadily brightening with time and might have now reached its peak, and constrain the emission process as non-thermal synchrotron emission where the cooling frequency c is above the X-ray band and the synchrotron frequency m is below the radio band. The very simple power-law spectrum extending for eight orders of magnitude in frequency enables the most precise measurement of the index p of the distribution of non-thermal relativistic electrons accelerated by a shock launched by a neutron star (NS)-NS merger to date. We find p = 2.17 0.01, which indicates that radiation from ejecta with Γ ∼ 3-10 dominates the observed emission. While constraining the nature of the emission process, these observations do not constrain the nature of the relativistic ejecta. We employ simulations of explosive outflows launched in NS ejecta clouds to show that the spectral and temporal evolution of the non-thermal emission from GW170817 is consistent with both emission from radially stratified quasi-spherical ejecta traveling at mildly relativistic speeds, and emission from off-axis collimated ejecta characterized by a narrow cone of ultra-relativistic material with slower wings extending to larger angles. In the latter scenario, GW170817 harbored a normal short gamma-ray burst (SGRB) directed away from our line of sight. Observations at t ≤ 200 days are unlikely to settle the debate, as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.
PDF
