Kilonovae

Physical mechanism of core-collapse supernovae that neutrinos drive

The current understanding of the mechanism of core-collapse supernovae (CCSNe), one of the most energetic events in the universe associated with the death of massive stars and the main formation channel of compact objects such as neutron stars and black holes, is reviewed for broad readers from different disciplines of science who may not be familiar with the object. Therefore, we emphasize the physical aspects than the results of individual model simulations, although large-scale high-fidelity simulations have played the most important roles in the progress we have witnessed in the past few decades. It is now believed that neutrinos are the most important agent in producing the commonest type of CCSNe. The so-called neutrino-heating mechanism will be the focus of this review and its crucial ingredients in micro- and macrophysics and in numerics will be explained one by one. We will also try to elucidate the remaining issues.

Magnetohydrodynamic Simulations of Self-Consistent Rotating Neutron Stars with Mixed Poloidal and Toroidal Magnetic Fields

We perform the first magnetohydrodynamic simulations in full general relativity of self-consistent rotating neutron stars (NSs) with ultrastrong mixed poloidal and toroidal magnetic fields. The initial uniformly rotating NS models are computed assuming perfect conductivity, stationarity, and axisymmetry. Although the specific geometry of the mixed field configuration can delay or accelerate the development of various instabilities known from analytic perturbative studies, all our models finally succumb to them. Differential rotation is developed spontaneously in the cores of our magnetars which, after sufficient time, is converted back to uniform rotation. The rapidly rotating magnetars show a significant amount of ejecta, which can be responsible for transient kilonova signatures. However, no highly collimated, helical magnetic fields or incipient jets, which are necessary for γ-ray bursts, arise at the poles of these magnetars by the time our simulations are terminated.

PIC methods in astrophysics: simulations of relativistic jets and kinetic physics in astrophysical systems

The Particle-In-Cell (PIC) method has been developed by Oscar Buneman, Charles Birdsall, Roger W. Hockney, and John Dawson in the 1950s and, with the advances of computing power, has been further developed for several fields such as astrophysical, magnetospheric as well as solar plasmas and recently also for atmospheric and laser-plasma physics. Currently more than 15 semi-public PIC codes are available which we discuss in this review. Its applications have grown extensively with increasing computing power available on high performance computing facilities around the world. These systems allow the study of various topics of astrophysical plasmas, such as magnetic reconnection, pulsars and black hole magnetosphere, non-relativistic and relativistic shocks, relativistic jets, and laser-plasma physics. We review a plethora of astrophysical phenomena such as relativistic jets, instabilities, magnetic reconnection, pulsars, as well as PIC simulations of laser-plasma physics (until 2021) emphasizing the physics involved in the simulations. Finally, we give an outlook of the future simulations of jets associated to neutron stars, black holes and their merging and discuss the future of PIC simulations in the light of petascale and exascale computing.

The first 5 years of gravitational-wave astrophysics

Gravitational waves are ripples in spacetime generated by the acceleration of astrophysical objects; a direct consequence of general relativity, they were first directly observed in 2015. Here, I review the first 5 years of gravitational-wave detections. More than 50 gravitational-wave events have been found, emitted by pairs of merging compact objects such as neutron stars and black holes. These signals yield insights into the formation of compact objects and their progenitor stars, enable stringent tests of general relativity, and constrain the behavior of matter at densities higher than that of an atomic nucleus. Mergers that emit both gravitational and electromagnetic waves probe the formation of short gamma-ray bursts and the nucleosynthesis of heavy elements, and they measure the local expansion rate of the Universe.

Multimessenger Binary Mergers Containing Neutron Stars: Gravitational Waves, Jets, and γ-Ray Bursts

Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate magnetic fields ten orders of magnitude larger than those currently constructed on earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger of a binary neutron star system, an event that generated a gravitational wave (GW) signal, a short γ -ray burst (sGRB), and a kilonova. This event serves as the highlight so far of the era of multimessenger astronomy. In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general relativistic magnetohydrodynamic simulations. Such mergers can lead to events like the one on August 2017, GW170817, and its EM counterparts, GRB 170817 and AT 2017gfo. In addition to exploring the GW emission from binary black hole-neutron star and neutron star-neutron star mergers, we also focus on their counterpart EM signals. In particular, we are interested in identifying the conditions under which a relativistic jet can be launched following these mergers. Such a jet is an essential feature of most sGRB models and provides the main conduit of energy from the central object to the outer radiation regions. Jet properties, including their lifetimes and Poynting luminosities, the effects of the initial magnetic field geometries and spins of the coalescing NSs, as well as their governing equation of state, are discussed. Lastly, we present our current understanding of how the Blandford-Znajek mechanism arises from merger remnants as the trigger for launching jets, if, when and how a horizon is necessary for this mechanism, and the possibility that it can turn on in magnetized neutron ergostars, which contain ergoregions, but no horizons.

Black hole-neutron star coalescence: Effects of the neutron star spin on jet launching and dynamical ejecta mass

Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (SGRBs) and kilonovas, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin (a BH/M BH = 0 and 0.75), the NS spin (a NS/M NS = -0.17, 0, 0.23, and 0.33), and the binary mass ratio (q = M BH:M NS = 3:1 and 5:1). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time Δt ~ 3500-5500M ~ 88-138(M NS/1.4 M ⊙) ms after the peak gravitational-wave amplitude, a magnetically driven jet is launched only for q = 3:1 regardless of the initial NS spin. The lifetime of the jets [Δt ~ 0.5-0.8(M NS/1.4 M ⊙) s] and their outgoing Poynting luminosity [L Poyn ~ 1051.5±0.5 erg/s] are consistent with typical SGRBs' luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic-field collimation for the other systems we consider. The mass range of dynamically ejected matter is 10-4.5-10-2(M NS/1.4 M ⊙) M ⊙, which can power kilonovas with peak bolometric luminosities L knova ~ 1040-1041.4 erg/s with rise times ≲6.5 h and potentially detectable by the LSST.

Systematic Uncertainty of Standard Sirens from the Viewing Angle of Binary Neutron Star Inspirals

The independent measurement of the Hubble constant with gravitational-wave standard sirens will potentially shed light on the tension between the local distance ladders and Planck experiments. Therefore, thorough understanding of the sources of systematic uncertainty for the standard siren method is crucial. In this Letter, we focus on two scenarios that will potentially dominate the systematic uncertainty of standard sirens. First, simulations of electromagnetic counterparts of binary neutron star mergers suggest aspherical emissions, so the binaries available for the standard siren method can be selected by their viewing angles. This selection effect can lead to ≳2% bias in Hubble constant measurement even with mild selection. Second, if the binary viewing angles are constrained by the electromagnetic counterpart observations but the bias of the constraints is not controlled under ∼10°, the resulting systematic uncertainty in the Hubble constant will be >3%. In addition, we find that both of the systematics cannot be properly removed by the viewing angle measurement from gravitational-wave observations. Comparing to the known dominant systematic uncertainty for standard sirens, the ≤2% gravitational-wave calibration uncertainty, the effects from the viewing angle appear to be more significant. Therefore, the systematic uncertainty from the viewing angle might be a major challenge before the standard sirens can resolve the tension in the Hubble constant, which is currently ∼9%.

Measuring the Hubble constant with a sample of kilonovae

Kilonovae produced by the coalescence of compact binaries with at least one neutron star are promising standard sirens for an independent measurement of the Hubble constant (H₀). Through their detection via follow-up of gravitational-wave (GW), short gamma-ray bursts (sGRBs) or optical surveys, a large sample of kilonovae (even without GW data) can be used for H₀ contraints. Here, we show measurement of H₀ using light curves associated with four sGRBs, assuming these are attributable to kilonovae, combined with GW170817. Including a systematic uncertainty on the models that is as large as the statistical ones, we find \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${H}_{0}=73.{8}_{-5.8}^{+6.3}\ {\rm{km}}\ {{\rm{s}}}^{-1}\ {{\rm{Mpc}}}^{-1}$$\end{document}H0=73.8−5.8+6.3kms−1Mpc−1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${H}_{0}=71.{2}_{-3.1}^{+3.2}\ {\rm{km}}\ {{\rm{s}}}^{-1}\ {{\rm{Mpc}}}^{-1}$$\end{document}H0=71.2−3.1+3.2kms−1Mpc−1 for two different kilonova models that are consistent with the local and inverse-distance ladder measurements. For a given model, this measurement is about a factor of 2-3 more precise than the standard-siren measurement for GW170817 using only GWs.

Kilonovae observations can be used to out constraints on the Hubble constant (H0). Here, the authors show H0 measurements by combining light curves of four short gamma-ray burts with GW170817 are about a factor of 2-3 more precise than the standard-siren measurements using only gravitational-waves.

Current Status and Developments of the Atomic Database on Rare-Earths at Mons University (DREAM)

The main purpose of the Database on Rare Earths At Mons University (DREAM) is to provide the scientific community with updated spectroscopic parameters related to lanthanide atoms (Z = 57–71) in their lowest ionization stages. The radiative parameters (oscillator strengths and transitions probabilities) listed in the database have been obtained over the past 20 years by the Atomic Physics and Astrophysics group of Mons University, Belgium, thanks to a systematic and extensive use of the pseudo-relativistic Hartree-Fock (HFR) method modified for taking core-polarization and core-penetration effects into account. Most of these theoretical results have been validated by the good agreement obtained when comparing computed radiative lifetimes and accurate experimental values measured by the time-resolved laser-induced fluorescence technique. In the present paper, we report on the current status and developments of the database that gathers radiative parameters for more than 72,000 spectral lines in neutral, singly-, doubly-, and triply-ionized lanthanides.

Magnetohydrodynamic simulations of binary neutron star mergers in general relativity: Effects of magnetic field orientation on jet launching

Binary neutron star mergers can be sources of gravitational waves coincident with electromagnetic counterpart emission across the spectrum. To solidify their role as multimessenger sources, we present fully 3D, general relativistic, magnetohydrodynamic simulations of highly spinning binary neutrons stars initially on quasicircular orbits that merge and undergo delayed collapse to a black hole. The binaries consist of two identical stars modeled as Γ = 2 polytropes with spin χ NS = 0.36 aligned along the direction of the total orbital angular momentum L. Each star is initially threaded by a dynamical unimportant interior dipole magnetic field. The field is extended into the exterior where a nearly force-free magnetosphere resembles that of a pulsar. The magnetic dipole moment μ is either aligned or perpendicular to L and has the same initial magnitude for each orientation. For comparison, we also impose symmetry across the orbital plane in one case where μ in both stars is aligned along L. We find that the lifetime of the transient hypermassive neutron star remnant, the jet launching time, and the ejecta (which can give rise to a detectable kilonova) are very sensitive to the magnetic field orientation. By contrast, the physical properties of the black hole + disk remnant, such as the mass and spin of the black hole, the accretion rate, and the electromagnetic (Poynting) luminosity, are roughly independent of the initial magnetic field orientation. In addition, we find imposing symmetry across the orbital plane does not play a significant role in the final outcome of the mergers. Our results suggest that, as in the black hole-neutron star merger scenario, an incipient jet emerges only when the seed magnetic field has a sufficiently large-scale poloidal component aligned to the initial orbital angular momentum. The lifetime [Δt ≳ 140(M NS/1.625 M ⊙) ms] and Poynting luminosities [L EM ≃ 1052 erg/s] of the jet, when it forms, are consistent with typical short gamma-ray bursts, as well as with the Blandford-Znajek mechanism for launching jets.

Prospects for observing and localizing gravitational-wave transients with Advanced LIGO, Advanced Virgo and KAGRA

We present our current best estimate of the plausible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next several years, with the intention of providing information to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals for the third (O3), fourth (O4) and fifth observing (O5) runs, including the planned upgrades of the Advanced LIGO and Advanced Virgo detectors. We study the capability of the network to determine the sky location of the source for gravitational-wave signals from the inspiral of binary systems of compact objects, that is binary neutron star, neutron star-black hole, and binary black hole systems. The ability to localize the sources is given as a sky-area probability, luminosity distance, and comoving volume. The median sky localization area (90% credible region) is expected to be a few hundreds of square degrees for all types of binary systems during O3 with the Advanced LIGO and Virgo (HLV) network. The median sky localization area will improve to a few tens of square degrees during O4 with the Advanced LIGO, Virgo, and KAGRA (HLVK) network. During O3, the median localization volume (90% credible region) is expected to be on the order of 10 5 , 10 6 , 10 7 Mpc 3 for binary neutron star, neutron star-black hole, and binary black hole systems, respectively. The localization volume in O4 is expected to be about a factor two smaller than in O3. We predict a detection count of 1 - 1 + 12 ( 10 - 10 + 52 ) for binary neutron star mergers, of 0 - 0 + 19 ( 1 - 1 + 91 ) for neutron star-black hole mergers, and 17 - 11 + 22 ( 79 - 44 + 89 ) for binary black hole mergers in a one-calendar-year observing run of the HLV network during O3 (HLVK network during O4). We evaluate sensitivity and localization expectations for unmodeled signal searches, including the search for intermediate mass black hole binary mergers.

Kilonovae

The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergoes rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a "kilonova", which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ∼ hour-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.

A nearby neutron-star merger explains the actinide abundances in the early Solar System

A growing body of evidence indicates that binary neutron-star mergers are the primary origin of heavy elements produced exclusively through rapid neutron capture^1-4 (the 'r-process'). As neutron-star mergers occur infrequently, their deposition of radioactive isotopes into the pre-solar nebula could have been dominated by a few nearby events. Although short-lived r-process isotopes-with half-lives shorter than 100 million years-are no longer present in the Solar System, their abundances in the early Solar System are known because their daughter products were preserved in high-temperature condensates found in meteorites⁵. Here we report that abundances of short-lived r-process isotopes in the early Solar System point to their origin in neutron-star mergers, and indicate substantial deposition by a single nearby merger event. By comparing numerical simulations with the early Solar System abundance ratios of actinides produced exclusively through the r-process, we constrain the rate of occurrence of their Galactic production sites to within about 1-100 per million years. This is consistent with observational estimates of neutron-star merger rates^6-8, but rules out supernovae and stellar sources. We further find that there was probably a single nearby merger that produced much of the curium and a substantial fraction of the plutonium present in the early Solar System. Such an event may have occurred about 300 parsecs away from the pre-solar nebula, approximately 80 million years before the formation of the Solar System.

Positrons and 511 keV Radiation as Tracers of Recent Binary Neutron Star Mergers

Neutron-rich material ejected from neutron star-neutron star (NS-NS) and neutron star-black-hole (NS-BH) binary mergers is heated by nuclear processes to temperatures of a few hundred keV, resulting in a population of electron-positron pairs. Some of the positrons escape from the outer layers of the ejecta. We show that the population of low-energy positrons produced by NS-NS and NS-BH mergers in the Milky Way can account for the observed 511-keV line from the Galactic center (GC). Moreover, we suggest how positrons and the associated 511-keV emission can be used as tracers of recent mergers. Recent discovery of 511-keV emission from the ultrafaint dwarf galaxy Reticulum II, consistent with a rare NS-NS merger event, provides a smoking-gun signature of our proposal.

High-energy emissions from neutron star mergers

In 2017, LIGO-Virgo collaborations reported detection of the first neutron star merger event, GW170817, which is accompanied by electromagnetic counterparts from radio to gamma rays. Although high-energy neutrinos were not detected from this event, mergers of neutron stars are expected to produce such high-energy particles. Relativistic jets are launched when neutron stars merge. If the jets contain protons, they can emit high-energy neutrinos through photomeson production. In addition, neutron star mergers produce massive and fast ejecta, which can be a source of Galactic high-energy cosmic rays above the knee. We briefly review what we learned from the multi-messenger event, GW170817, and discuss prospects for multi-messenger detections and hadronic cosmic-ray production related to the neutron star mergers.

New Constraints on Radii and Tidal Deformabilities of Neutron Stars from GW170817

We explore in a parameterized manner a very large range of physically plausible equations of state (EOSs) for compact stars for matter that is either purely hadronic or that exhibits a phase transition. In particular, we produce two classes of EOSs with and without phase transitions, each containing one million EOSs. We then impose constraints on the maximum mass (M<2.16  M_{⊙}) and on the dimensionless tidal deformability (Λ[over ˜]<800) deduced from GW170817, together with recent suggestions of lower limits on Λ[over ˜]. Exploiting more than 10^{9} equilibrium models for each class of EOSs, we produce distribution functions of all the stellar properties and determine, among other quantities, the radius that is statistically most probable for any value of the stellar mass. In this way, we deduce that the radius of a purely hadronic neutron star with a representative mass of 1.4  M_{⊙} is constrained to be 12.00<R_{1.4}/km<13.45 at a 2σ confidence level, with a most likely value of R[over ¯]_{1.4}=12.39  km; similarly, the smallest dimensionless tidal deformability is Λ[over ˜]_{1.4}>375, again at a 2σ level. On the other hand, because EOSs with a phase transition allow for very compact stars on the so-called "twin-star" branch, small radii are possible with such EOSs although not probable, i.e., 8.53<R_{1.4}/km<13.74 and R[over ¯]_{1.4}=13.06  km at a 2σ level, with Λ[over ˜]_{1.4}>35.5 at a 3σ level. Finally, since these EOSs exhibit upper limits on Λ[over ˜], the detection of a binary with a total mass of 3.4  M_{⊙} and Λ[over ˜]_{1.7}>461 can rule out twin-star solutions.

High-energy astrophysics and the search for sources of gravitational waves

The dawn of the gravitational-wave (GW) era has sparked a greatly renewed interest into possible links between sources of high-energy radiation and GWs. The most luminous high-energy sources-gamma-ray bursts (GRBs)-have long been considered as very likely sources of GWs, particularly from short-duration GRBs, which are thought to originate from the merger of two compact objects such as binary neutron stars and a neutron star-black hole binary. In this paper, we discuss: (i) the high-energy emission from short-duration GRBs; (ii) what other sources of high-energy radiation may be observed from binary mergers; and (iii) how searches for high-energy electromagnetic counterparts to GW events are performed with current space facilities. While current high-energy facilities, such as Swift and Fermi, play a crucial role in the search for electromagnetic counterparts, new space missions will greatly enhance our capabilities for joint observations. We discuss why such facilities, which incorporate new technology that enables very wide-field X-ray imaging, are required if we are to truly exploit the multi-messenger era.This article is part of a discussion meeting issue 'The promises of gravitational-wave astronomy'.

Prospects for observing and localizing gravitational-wave transients with Advanced LIGO, Advanced Virgo and KAGRA

We present possible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next decade, with the intention of providing information to the astronomy community to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals, and study the capability of the network to determine the sky location of the source. We report our findings for gravitational-wave transients, with particular focus on gravitational-wave signals from the inspiral of binary neutron star systems, which are the most promising targets for multi-messenger astronomy. The ability to localize the sources of the detected signals depends on the geographical distribution of the detectors and their relative sensitivity, and [Formula: see text] credible regions can be as large as thousands of square degrees when only two sensitive detectors are operational. Determining the sky position of a significant fraction of detected signals to areas of 5-[Formula: see text] requires at least three detectors of sensitivity within a factor of [Formula: see text] of each other and with a broad frequency bandwidth. When all detectors, including KAGRA and the third LIGO detector in India, reach design sensitivity, a significant fraction of gravitational-wave signals will be localized to a few square degrees by gravitational-wave observations alone.

Optical observations of LIGO source GW 170817 by the Antarctic Survey Telescopes at Dome A, Antarctica

The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW 170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the i-band magnitude declining from 17.23±0.13 magnitude to 17.72±0.09 magnitude in ~1.8 h. The brightness and time evolution of the optical transient associated with GW 170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about ∼10^-2 solar mass of radioactive material at a speed of up to 30% the speed of light.

Swift and NuSTAR observations of GW170817: Detection of a blue kilonova

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and x-ray observations by Swift and the Nuclear Spectroscopic Telescope Array of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading UV emission indicates a high mass (≈0.03 solar masses) wind-driven outflow with moderate electron fraction (Y_e ≈ 0.27). Combined with the x-ray limits, we favor an observer viewing angle of ≈30° away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultrarelativistic, highly collimated ejecta (a γ-ray burst afterglow).

A kilonova as the electromagnetic counterpart to a gravitational-wave source

Gravitational waves were discovered with the detection of binary black-hole mergers and they should also be detectable from lower-mass neutron-star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal. This signal is luminous at optical and infrared wavelengths and is called a kilonova. The gravitational-wave source GW170817 arose from a binary neutron-star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC 4993, which is spatially coincident with GW170817 and with a weak, short γ-ray burst. The transient has physical parameters that broadly match the theoretical predictions of blue kilonovae from neutron-star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 ± 0.01 solar masses, with an opacity of less than 0.5 square centimetres per gram, at a velocity of 0.2 ± 0.1 times light speed. The power source is constrained to have a power-law slope of -1.2 ± 0.3, consistent with radioactive powering from r-process nuclides. (The r-process is a series of neutron capture reactions that synthesise many of the elements heavier than iron.) We identify line features in the spectra that are consistent with light r-process elements (atomic masses of 90-140). As it fades, the transient rapidly becomes red, and a higher-opacity, lanthanide-rich ejecta component may contribute to the emission. This indicates that neutron-star mergers produce gravitational waves and radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

Early spectra of the gravitational wave source GW170817: Evolution of a neutron star merger

On 17 August 2017, Swope Supernova Survey 2017a (SSS17a) was discovered as the optical counterpart of the binary neutron star gravitational wave event GW170817. We report time-series spectroscopy of SSS17a from 11.75 hours until 8.5 days after the merger. Over the first hour of observations, the ejecta rapidly expanded and cooled. Applying blackbody fits to the spectra, we measured the photosphere cooling from [Formula: see text] to [Formula: see text] kelvin, and determined a photospheric velocity of roughly 30% of the speed of light. The spectra of SSS17a began displaying broad features after 1.46 days and evolved qualitatively over each subsequent day, with distinct blue (early-time) and red (late-time) components. The late-time component is consistent with theoretical models of r-process-enriched neutron star ejecta, whereas the blue component requires high-velocity, lanthanide-free material.