The merger of two neutron stars is predicted to give rise to three major detectable phenomena: a short burst of γ-rays, a gravitational-wave signal, and a transient optical–near-infrared source powered by the synthesis of large amounts of very heavy elements via rapid neutron capture (the r-process)^1,2,3. Such transients, named ‘macronovae’ or ‘kilonovae’^4,5,6,7, are believed to be centres of production of rare elements such as gold and platinum⁸. The most compelling evidence so far for a kilonova was a very faint near-infrared rebrightening in the afterglow of a short γ-ray burst^9,10 at redshift z = 0.356, although findings indicating bluer events have been reported¹¹. Here we report the spectral identification and describe the physical properties of a bright kilonova associated with the gravitational-wave source¹² GW170817 and γ-ray burst^13,14 GRB 170817A associated with a galaxy at a distance of 40 megaparsecs from Earth. Using a series of spectra from ground-based observatories covering the wavelength range from the ultraviolet to the near-infrared, we find that the kilonova is characterized by rapidly expanding ejecta with spectral features similar to those predicted by current models^15,16. The ejecta is optically thick early on, with a velocity of about 0.2 times light speed, and reaches a radius of about 50 astronomical units in only 1.5 days. As the ejecta expands, broad absorption-like lines appear on the spectral continuum, indicating atomic species produced by nucleosynthesis that occurs in the post-merger fast-moving dynamical ejecta and in two slower (0.05 times light speed) wind regions. Comparison with spectral models suggests that the merger ejected 0.03 to 0.05 solar masses of material, including high-opacity lanthanides.

Long-duration (>2 s) γ-ray bursts that are believed to originate from the death of massive stars are expected to be accompanied by supernovae. GRB 060614, that lasted 102 s, lacks a supernova-like emission down to very stringent limits and its physical origin is still debated. Here we report the discovery of near-infrared bump that is significantly above the regular decaying afterglow. This red bump is inconsistent with even the weakest known supernova. However, it can arise from a Li-Paczyński macronova—the radioactive decay of debris following a compact binary merger. If this interpretation is correct, GRB 060614 arose from a compact binary merger rather than from the death of a massive star and it was a site of a significant production of heavy r-process elements. The significant ejected mass favours a black hole–neutron star merger but a double neutron star merger cannot be ruled out.

GRB 050709 was the first short Gamma-ray Burst (sGRB) with an identified optical counterpart. Here we report a reanalysis of the publicly available data of this event and the discovery of a Li-Paczynski macronova/kilonova that dominates the optical/infrared signal at t>2.5 days. Such a signal would arise from 0.05 "Equation missing" r-process material launched by a compact binary merger. The implied mass ejection supports the suggestion that compact binary mergers are significant and possibly main sites of heavy r-process nucleosynthesis. Furthermore, we have reanalysed all afterglow data from nearby short and hybrid GRBs (shGRBs). A statistical study of shGRB/macronova connection reveals that macronova may have taken place in all these GRBs, although the fraction as low as 0.18 cannot be ruled out. The identification of two of the three macronova candidates in the I-band implies a more promising detection prospect for ground-based surveys.

The jet breaks in the afterglow light curves of short gamma-ray bursts (SGRBs), rarely detected so far, are crucial for estimating the half-opening angles of the ejecta (θ _j) and hence the neutron star merger rate. In this work, we report the detection of jet decline behaviors in GRB 150424A and GRB 160821B, and find θ _j ∼ 0.1 rad. Together with five events reported before 2015 and three others “identified” recently (GRB 050709, GRB 060614, and GRB 140903A), we have a sample consisting of nine SGRBs and one long-short GRB with reasonably estimated θ _j. In particular, three Swift bursts in the sample have redshifts z ≤ 0.2, with which we estimate the local neutron star merger rate density to be or if the narrowly beamed GRB 061201 is excluded. Inspired by the typical θ _j ∼ 0.1 rad found currently, we further investigate whether the off-beam GRBs (in the uniform jet model) or the off-axis events (in the structured jet model) can significantly enhance the GRB/GW association. For the former, the enhancement is at most moderate, while for the latter the enhancement can be much greater and a high GRB/GW association probability of ∼10% is possible. We also show that the data of GRB 160821B may contain a macronova/kilonova emission component with a temperature of ∼3100 K at ∼3.6 days after the burst and more data are needed to ultimately clarify.

The historic first joint detection of both gravitational-wave and electromagnetic emission from a binary neutron star merger cemented the association between short gamma-ray bursts (SGRBs) and compact object mergers, as well as providing a well-sampled multi-wavelength light curve of a radioactive kilonova (KN) for the first time. Here, we compare the optical and near-infrared light curves of this KN, AT 2017gfo, to the counterparts of a sample of nearby (z < 0.5) SGRBs to characterize their diversity in terms of their brightness distribution. Although at similar epochs AT 2017gfo appears fainter than every SGRB-associated KN claimed so far, we find three bursts (GRBs 050509B, 061201, and 080905A) where, if the reported redshifts are correct, deep upper limits rule out the presence of a KN similar to AT 2017gfo by several magnitudes. Combined with the properties of previously claimed KNe in SGRBs this suggests considerable diversity in the properties of KN drawn from compact object mergers, despite the similar physical conditions that are expected in many NS–NS mergers. We find that observer angle alone is not able to explain this diversity, which is likely a product of the merger type (NS–NS versus NS–BH) and the detailed properties of the binary (mass ratio, spins etc.). Ultimately disentangling these properties should be possible through observations of SGRBs and gravitational-wave sources, providing direct measurements of heavy element enrichment throughout the universe.