Research Article

Electromagnetic evidence that SSS17a is the result of a binary neutron star merger

See allHide authors and affiliations

Science  22 Dec 2017:
Vol. 358, Issue 6370, pp. 1583-1587
DOI: 10.1126/science.aaq0073

Photons from a gravitational wave event

Two neutron stars merging together generate a gravitational wave signal and have also been predicted to emit electromagnetic radiation. When the gravitational wave event GW170817 was detected, astronomers rushed to search for the source using conventional telescopes (see the Introduction by Smith). Coulter et al. describe how the One-Meter Two-Hemispheres (1M2H) collaboration was the first to locate the electromagnetic source. Drout et al. present the 1M2H measurements of its optical and infrared brightness, and Shappee et al. report their spectroscopy of the event, which is unlike previously detected astronomical transient sources. Kilpatrick et al. show how these observations can be explained by an explosion known as a kilonova, which produces large quantities of heavy elements in nuclear reactions.

Science, this issue p. 1556, p. 1570, p. 1574, p. 1583; see also p. 1554


Eleven hours after the detection of gravitational wave source GW170817 by the Laser Interferometer Gravitational-Wave Observatory and Virgo Interferometers, an associated optical transient, SSS17a, was identified in the galaxy NGC 4993. Although the gravitational wave data indicate that GW170817 is consistent with the merger of two compact objects, the electromagnetic observations provide independent constraints on the nature of that system. We synthesize the optical to near-infrared photometry and spectroscopy of SSS17a collected by the One-Meter Two-Hemisphere collaboration, finding that SSS17a is unlike other known transients. The source is best described by theoretical models of a kilonova consisting of radioactive elements produced by rapid neutron capture (the r-process). We conclude that SSS17a was the result of a binary neutron star merger, reinforcing the gravitational wave result.

The gravitational wave event GW170817 (1) was detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) (2) and Virgo Interferometer (Virgo) (3). An electromagnetic counterpart, SSS17a, was subsequently identified 11 hours after the gravitational wave event (4). Combining the gravitational wave and electromagnetic signatures could inform models of the physical mechanisms involved in compact object mergers (5) and the rate at which compact object binaries produce heavy elements (6).

Two seconds after GW170817 and 11 hours before SSS17a was identified, the Fermi Gamma-ray Space Telescope (Fermi) and International Gamma-Ray Astrophysics Laboratory (INTEGRAL) detected a γ-ray burst (GRB), GRB 170817A, coincident with GW170817 (7, 8). It has long been hypothesized that short-duration GRBs (SGRBs) are the result of compact binary mergers involving neutron stars and/or black holes (9). Rapidly fading GRB afterglows and kilonova emission powered by the decay of heavy elements have been hypothesized as potential counterparts to compact binary mergers (5).

Whereas the LIGO-Virgo collaboration analyze GW170817 and the information encoded in the waveform of the gravitational signal, as astrophysicists we seek to place stringent constraints on the nature and evolutionary behavior of the progenitor system by interpreting the electromagnetic properties of SSS17a and its galactic environment. We take an agnostic approach motivated by the question: If we had no constraints from gravitational waves and SSS17a was discovered as part of a typical optical transient survey, what could we conclude about the nature of the progenitor system? Such considerations may be unavoidable in future events given that the LIGO and Virgo interferometers will undergo temporary off-line periods. In addition, deep, wide-field surveys, such as the Large Synoptic Survey Telescope, are expected to be more efficient at detecting SSS17a-like objects than targeted follow-up of LIGO-Virgo detections (10).

SSS17a was identified by the One-Meter Two-Hemisphere (1M2H) collaboration and the Swope Supernova Survey as part of a search for gravitational wave counterparts on 17 August 2017 (4). The discovery image revealed a transient source with i-band magnitude of Embedded Image mag ~10.6 arcsec (″) from the center of the S0-type galaxy NGC 4993 (4). We synthesize optical and near-infrared photometry of the transient source (11, 12), its environment (13), and our theoretical interpretation of these data.

Galactic environment in NGC 4993

We first consider the host galaxy of GW170817, and the location of SSS17a within it. The host galaxy morphology and mass provide information about its age, the populations of stars that it contains, and their dynamical properties. Comparison between different classes of astrophysical transients and their host galaxies is therefore an indirect constraint on the properties of their progenitor systems.

The host galaxy of SSS17a is NGC 4993, an S0-type galaxy with prominent lanes of cosmic dust concentrated at its core and reaching out nearly to the projected position of SSS17a (13). We analyzed the surface brightness of NGC 4993 using a Hubble Space Telescope (HST) image (14) and found that it decreases with projected radius (R) and a scaling constant (k) as exp(–kR1/4). This profile is a de Vaucouleurs profile and is typical for S0-type galaxies. Assuming this profile, we find that NGC 4993 has an effective radius Embedded Image kpc (17″).

Fitting stellar population models to 13-band photometry of NGC 4993 (14), we find a galactic stellar mass of Embedded Image, where M is the stellar mass of NGC 4993 and Embedded Image is the mass of the Sun. This stellar mass is consistent with the host galaxies’ masses of SGRBs (15).

We detect no point source at the position of SSS17a in the pre-trigger HST image in the Embedded Image filter, with a V-band apparent magnitude limit of Embedded Image mag. This is equivalent to a limit on the absolute V-band magnitude Embedded Image mag at Embedded Image Mpc, which is the distance to NGC 4993 (14). This limit rules out the most luminous and massive progenitor stars, but is still consistent with the majority of massive stars (16) and with low-mass stars or compact binary systems, both of which are optically dim.

SSS17a is offset 10.6″ from the center of NGC 4993, or 2.0 kpc in projection. This corresponds to an offset of Embedded Image, which is small compared to the typical offsets seen in population studies of SGRBs (17). The small offset is indicative of either a progenitor velocity less than the escape velocity of the galaxy [Embedded Image350 km sEmbedded Image at the location of the transient (13)] or a system that was recently kicked from the site where it formed. If the latter is true, then SSS17a must be close to the site where it formed.

It is thought that the progenitor systems of SGRBs are formed more frequently in globular clusters, which have high stellar densities and more dynamical interactions per star (6). There are Embedded Image100 likely globular clusters in HST imaging of NGC 4993, the closest offset by 290 pc in projection from SSS17a (13). Even if the progenitor system was ejected from the closest cluster with a relative velocity of 10 km sEmbedded Image, this corresponds to a travel time of 28 million years (My), meaning that we cannot exclude the possibility that the SSS17a progenitor system originated in a globular cluster.

If the SSS17a progenitor system had received a large kick velocity of >350 km s−1, it would have become unbound and exited its host galaxy after a relatively short time (<20 My). Given the lack of recent star formation seen in NGC 4993, it is unlikely that the progenitor system could have been as young as 20 My. In addition, such large kicks are physically unlikely given that they tend to unbind a compact binary (18). Thus, it is more likely that the SSS17a progenitor was bound to its host and merged over a long time scale. This hypothesis also agrees with expectations of the progenitor systems of SGRBs in S0-type galaxies, which are thought to be systematically older and occur at lower redshifts (19).

The optical properties of SSS17a

For many known classes of transients, there is sufficient observational evidence to connect them to specific progenitor systems. For instance, pre-explosion images of some supernovae reveal a massive star located at the precise position of the transient (16). Indirect evidence from optical light curves and spectroscopy also associate other transients (e.g., Type Ia and Ic supernovae) with either white-dwarf or massive-star progenitors (16, 20).

SSS17a has different optical properties from other known astrophysical transient objects (10). Figure 1 shows the photometric light curves, in which optical emission rises in Embedded Image day, then fades rapidly, with a rapid color evolution to the red (11). The spectra shown in Fig. 2 exhibit similarly rapid evolution; 11 hours after the trigger the spectrum is blue and smooth, but it transitions within a few days to a redder spectrum with at least one clear spectral bump near Embedded Image Å (12). These optical spectra lack the numerous absorption features typically seen in the spectra of ordinary supernovae (10, 12).

Fig. 1 Photometry of SSS17a compared to fitted kilonova models.

(A) UV-to-NIR photometry of SSS17a from Embedded Image hours after the BNS merger to +Embedded Image days (11). Overplotted lines are our best-fitting kilonova model in each band. (B) Residuals (in magnitudes) between each photometry measurement and our best-fitting model. (C) The integrated luminosity of our best-fitting kilonova model compared with the total integrated luminosity of SSS17a (11). We also show the luminosity of the individual blue and red components of our kilonova model. (D) The derived temperature of our kilonova model compared with the temperature derived by fitting a blackbody SED to each epoch (11).

Fig. 2 Spectral series of SSS17a and kilonova models.

Our flux-calibrated Magellan spectra of SSS17a (12). Each spectrum is labeled with the epoch in rest-frame days since the BNS merger. We overplot our best-fitting kilonova model at each corresponding epoch.

Nearly every known class of astrophysical transients is inconsistent with at least one of the following properties of SSS17a: the rise time of Embedded Image hours, subsequent fading at Embedded Image mag dayEmbedded Image, the color difference in V- and H-band magnitudes (VH) that transitions from Embedded Image to Embedded Image mag over 4 days, and the nearly featureless optical spectra at all epochs (1012). Among previously observed events dominated by thermal continuum emission, the most similar class is rapidly evolving blue transients (21); however, these events have peak magnitudes that are too luminous, typically occur in star-forming galaxies, do not have the observed color evolution, and do not fade as quickly as SSS17a (10). Nonthermal relativistic sources such as GRB afterglows can produce rapidly fading transients but are not expected to produce the quasi-blackbody spectra that are observed (Fig. 2) (12).

SSS17a’s observational properties suggest a different progenitor system. Various types of short-duration transients have been theorized to come from white dwarf or massive-star explosions, but none of these model predictions resemble SSS17a. For example, models of surface detonations on a white dwarf (22) could explain the peak luminosity and spectra dominated by thermal continuum emission. But this model cannot explain the rapid rise time and color evolution of SSS17a detailed above. In particular, SSS17a’s combination of luminosity, rise time, and red color implies a small ejecta mass (Embedded Image) consisting predominately of high-opacity, radioactive material (10, 11). Therefore, we turn to models of other astrophysical transients to describe our observations.

Comparison with kilonova models

Kilonova models and ultraviolet to near-infrared photometry of SSS17a

The early evolution, colors, and ultraviolet (UV) to near-infrared (IR) luminosities of SSS17a are well constrained by its light curve (11). Observations of the optical transient began roughly 11 hours after the LIGO-Virgo trigger (4), and within 1 day of discovery the spectral energy distribution of the transient had been determined from the far-UV (Swift W2-band; Embedded Image Å) to the near-IR (Magellan/FourStar K-band; Embedded Image Å) (14).

Models of compact-object mergers involving neutron stars predict that optical and IR emission should be produced by radioactively powered thermal emission (a kilonova) (5). This emission, which is produced by subrelativistic ejecta thrown off from the neutron stars during or just after the merger, is thought to be observable in all directions. This physical mechanism contrasts with that of SGRBs, which are beamed along certain lines of sight.

Numerical simulations predict that neutron-rich material will be ejected from the system and assemble into heavy elements via rapid neutron capture (r-process) nucleosynthesis (23). The optical and IR appearance of kilonovae are distinguished by the unusual composition of the ejecta. Heavy r-process ejecta (atomic mass Embedded Image) include a large fraction of lanthanide elements, which have complex f-shell valence electron structures, producing millions of bound-bound line transitions (24). This leads to a high optical opacity κ ≈ 5 to to 10 cm2 g−1, which causes the kilonova spectrum to peak in the IR and produces a red transient lasting for several days (25). Ejecta composed primarily of lighter r-process products (with mass number Embedded Image) can be relatively lanthanide-free and have lower opacities, leading to a bluer transient (26).

The peak luminosity, characteristic time scale, and spectral energy distribution of a kilonova is largely determined by the total ejecta mass (Embedded Image), characteristic ejecta velocity (Embedded Image), and lanthanide fraction in the ejecta (Embedded Image). To test whether our data are consistent with a kilonova and to constrain the model parameters, we analyze the UV, optical, and IR emission of SSS17a (11) using numerical kilonova models (27).

The UV-to-IR observations of SSS17a qualitatively match the numerical kilonova model (Figs. 1 and 2) and indicate that the ejecta contained two distinct components. The longer-duration IR light curves require the high opacities of a lanthanide-rich component, whereas the short-duration UV–optical light curves require a lower-opacity, relatively lanthanide-free component. We model the panchromatic data by summing a red kilonova component with Embedded Image, Embedded Image, Embedded Image, and a blue kilonova with Embedded Image, Embedded Image, and a gradient in the lanthanide fraction spanning Embedded Image to Embedded Image (14). The blue kilonova reproduces both the UV luminosity and decline rate of SSS17a at early times. We do not attempt to tune the model parameters or compositional stratification.

Specific features in the data provide evidence for two-components of kilonova ejecta. The distinctive shape of the bolometric light curve—which declines at Embedded Image mag dayEmbedded Image at times Embedded Image days, then plateaus from t ≈ 3 to 8 days (11)—is naturally reproduced by the sum of blue and red kilonova components. Such a spectral energy distribution can be formed by the superposition of two quasi-blackbody sources of different temperatures and was predicted to be a signature of two-component kilonovae (25).

Around 10 to 15 days after merger, the color temperature of SSS17a asymptotes to a value near 2500 K (11). This behavior can be naturally explained by the recombination of open f-shell lanthanides, which occurs around this temperature (24). However, at times >10 days after merger, the model H-band luminosity is fainter than the observations. This may reflect limitations in the radiation transport model calculations, which, owing to uncertainties in the atomic data inputs, may poorly fit the observed H-band luminosities. This deviation is reflected in both the derived luminosity and temperature for SSS17a (11). Luminosities derived Embedded Image days after merger assume that the temperature remains around 2000 K as we have no multiband constraints beyond this point, although the model appears to match the K-band data.

The characteristics of the blue kilonova are well constrained by the early-time UV and optical detections of SSS17a. In binary neutron star (BNS) systems, it is thought that the merger scenario might involve two stars that collide violently; these dynamically expel hot polar ejecta that is neutrino irradiated and can produce a lighter r-process (lanthanide-free) blue kilonova (28). The early component of SSS17a could broadly match predictions of BNS kilonovae (25).

It is thought that kilonovae from black hole–neutron star (BH-NS) systems are generated in systems with relatively small mass ratios—that is, with neutron stars disrupted by relatively low-mass (Embedded Image), high-spin (Embedded Image) black holes (29). In this range, it is possible for tidal forces to eject a small fraction of the neutron star material in tidal tails and for additional disrupted material to assemble into a disk. If the mass ratio is too high, however, the neutron star would not be disrupted until it had already passed below the black hole’s event horizon. The dynamical ejecta in these systems should produce a lanthanide-rich, red kilonova (23).

The photometry clearly exhibits a rapidly fading (Embedded Image mag dayEmbedded Image), blue kilonova in the far-UV Swift bands, where the transient was only detected within the first day of observations (Fig. 1). We find that SSS17a exhibits rapid color evolution from blue to red over the course of days (11), which contrasts with that of other known optical transients, which have characteristic time scales of weeks or months (10). This rapid transformation in color is consistent, however, with the kilonovae expected from BNS merger models.

We model the blue emission of SSS17a using an abundance gradient in the ejecta such that the lanthanide mass fraction is very low (Embedded Image) in the outermost layers but increases inward to a higher, although still trace, fraction (Embedded Image) (14). A model with such a compositional gradient provides a better fit to the color evolution of SSS17a and is consistent with the physical expectation that the fastest-moving material may have experienced the highest neutrino irradiation.

Spectroscopy of SSS17a: Isolating kilonova features

As a check on the validity of our kilonova models, we overplot spectra of SSS17a (12) with the synthetic kilonova spectra from our model in Fig. 2. The observed spectra have been dereddened to remove the effects of dust obscuration, and the Doppler shift caused by the recessional velocity of the host galaxy has been removed. There is qualitative agreement between the synthetic spectra and the optical emission from SSS17a, with a transition from blue at early times to red at later times (12). Given the high velocities associated with our synthetic spectra (≥0.1c) and the lanthanide features across the UV, optical, and IR, we can reproduce the smooth continuum shape of SSS17a with approximately the observed temperature at all epochs.

The presence of distinct blue and red thermal emission is especially apparent in the day 3 to 5 spectra. The spectra appear to be split into two smooth continua with a transition around 9000 to 10,000 Å, indicating that the spectrum can be divided into two components with T > 3300 K and T < 2900 K, which around day 3 each contribute roughly the same amount of flux.

Our kilonova model, which is the sum of distinct blue and red ejecta components, reproduces the general behavior of the observed color evolution, and reproduces the characteristic double-peaked spectral continuum observed at day 3.46. Quantitatively, however, the model evolution shows several deficiencies. The blue kilonova component does not fade rapidly enough, such that the model overpredicts the optical flux at days 4.51 and 7.45. This suggests that the lanthanide gradient of the blue component of ejecta may be even steeper than assumed here, and the inner layers of ejecta even more lanthanide rich. This would result in a more rapid suppression of the optical flux over time. Thus, the time series of spectra provide insight into the layered compositional structure of the ejecta.

On the origin of r-process elements

Theoretical studies of nucleosynthesis have discerned the thermodynamic conditions required for the r-process (30). The specific astrophysical site, on the other hand, has not been unambiguously identified.

Two main alternatives have been commonly discussed. The first attributed the production of r-process elements to the neutron-rich regions in the outer layers of a nascent neutron star in a core-collapse supernova explosion (31), while the second suggested the ejecta from a neutron star merger as a likely site (23). These two avenues are thought to eject different quantities of r-process material, and the expectation is that such differences should be imprinted in the r-process enrichment pattern of stars in the Galactic halo (32).

Because both models imply that elements are synthesized in an explosive event, a strict lower limit for the dilution of such elements can be derived by calculating the mass swept up once the blast wave has lost most of its energy to radiation. The inferred minimum r-process mass per event as a function of explosion energy is shown in Fig. 3, which was derived using the r-process concentrations observed in low-metallicity stars in the Milky Way (14). Also shown are the constraints on the ejecta mass and kinetic energy that we derive for SSS17a, which demonstrate that the ejecta in this event are consistent with the most stringent constraints on r-process production derived using low-metallicity stars.

Fig. 3 Mass of r-process elements and kinetic energy in dynamical ejecta.

Strict lower limit on the mass (Embedded Image) and kinetic energy (Ek) of r-process ejecta derived from the r-process content of low-metallicity stars (dashed line) (14). This argument provides a range for any viable model to be tested—for example, from magnetars and Type II supernovae (14)—and illustrates that SSS17a-like mergers satisfy the mass and energy ejecta constraints required for viable scenarios.

The r-process abundance patterns in these stars (in particular elements with A > 140) closely follow that of the solar system, suggesting that the rate of enrichment has been relatively constant over long periods of time in Galactic history (33). The large star-to-star r-process dispersions that are seen in Galactic halo stars probably suggest that they formed in an early and nucleosynthetically unmixed Milky Way (32). These confined heavy-element inhomogeneities are then expected to be washed out as more events take place and the merger products mix effectively throughout the host galaxy. This is the case for SSS17a, which will likely deposit all of its newly synthesized heavy-element–rich ejecta into NGC 4993.

The rate of SSS17a-like events in the Galaxy is quite uncertain. Based on the observed properties of SSS17a, a conservative limit of the volumetric rate of SSS17a-like events, Embedded Image MpcEmbedded Image yearEmbedded Image (Embedded Image GpcEmbedded Image yearEmbedded Image), has been derived (10). Here we assume an average rate Embedded Image MpcEmbedded Image yearEmbedded Image over the history of the local universe. This value is consistent with rate limits derived from the first observing run performed by the Advanced LIGO Interferometer, although it is a factor of Embedded Image smaller than predictions for what was necessary to detect a BNS merger during the second observing run (34). However, assuming one Milky Way–like galaxy per (4.4 Mpc)3 (34), the assumed rate implies a Milky Way rate RMW = 25 My–1, which is consistent with predictions from binary population synthesis (18).

Our estimates for SSS17a show that an average total mass of about Embedded Image was ejected in the event. If this is representative of the average for BNS mergers, the cumulative mass of r-process elements ejected from all past SSS17a-like events over the Galactic history Embedded Image is Embedded Image.

This value is comparable to the present r-process mass inventory in stars in the Milky Way. We determine that the total mass of r-process elements per Milky Way–like galaxy is Embedded Image, where Embedded Image is the total Galactic mass in stars and gas (35) and Embedded Image is the total mass fraction of r-process nuclei in the solar system (36). This total mass of r-process elements consists of 78% light r-process elements (A < 140; likely powering the blue kilonova) and 22% in the main heavier component (33) (A > 140; likely powering the red kilonova). This indicates that BNS mergers such as SSS17a can produce sufficient r-process material to be a major source of r-process elements.

The progenitor system of SSS17a

SSS17a has been observed throughout the entire electromagnetic spectrum, providing clues to the nature of the progenitor system. The radiated γ-ray flux varies on time scales of ~0.1 s (7, 8), which demonstrates that the site producing γ-ray photons must be very compact. Estimates of the total emitted energy suggest that the event must have emitted >1050 erg (37). Given the requirements of energy, outflow velocity, and γ-ray compactness, it is unlikely that mass can be converted into outflow energy with efficiencies better than a few percent. As a result, the central source giving rise to GW170817/GRB170817A/SSS17a must process upwards of Embedded Image of material through a region that is not much larger than the size of a neutron star (38).

Based on these electromagnetic constraints, several questions remain regarding the central engine that triggered SSS17a. In particular, does the progenitor system involve a black hole and neutron star or two neutron stars? If so, can we distinguish between the two?

The UV-to-IR observations of SSS17a do not strongly constrain the primary source of this radiated emission; rather, they are signatures of secondary reprocessing by the radioactive decay of r-process elements. From this, we can determine the properties of the binary merger only by a chain of less certain inferences. However, one such inference involves the properties of the dynamical ejecta—namely, its mass and velocity. These properties are predicted to correlate primarily with the mass ratio of the compact binary (39).

Comparing our results with numerical simulations (14) in fig. S1, we estimate the mass ratio of the binary system to be ~0.75, which is consistent with findings from the gravitational wave data alone (1). This conclusion is based on a specific set of simulations, but it is generally consistent with a BNS system progenitor where the mass ratio of the progenitor binary is more likely to be close to 1. However, these constraints on the mass ratio depend sensitively on general-relativistic effects and the equation of state of nuclear matter (40).


SSS17a was a peculiar optical transient with a luminosity, rise time, and rate of decline unlike those of known classes of optical transients (1012). It occurred well within the effective radius of the massive, S0-type host galaxy NGC 4993. Either the progenitor star system was bound to its host galaxy, or there was a <20-My delay time between the formation of the system and the optical transient. Given the morphology of the host and the lack of any obvious signs of star formation, we infer that the former is more likely.

The light curves of SSS17a are fairly well matched by kilonova models. Such models involve the merger of a compact binary system with at least one neutron star, ejecting and synthesizing r-process elements. From these models, we infer that SSS17a had an early blue kilonova component with Embedded Image, vk = 0.25c, and a longer-lasting red kilonova with Embedded Image and vk = 0.15 ± 0.03c. From this mass estimate, we find that the rate of r-process nucleosynthesis inferred from SSS17a is consistent with the production in our Galaxy as inferred from halo stars and the abundance of r-process elements in the solar system.

From the properties of the dynamical ejecta as inferred from the UV, optical, and IR light curves of SSS17a, we infer a mass ratio of ~0.75 for the compact object binary progenitor. We therefore conclude that SSS17a is most consistent with a BNS progenitor system.

Supplementary Materials

Materials and Methods

Supplementary Text

Fig. S1

References (4162)

References and Notes

  1. Materials and methods are available as supplementary materials.
Acknowledgments: We thank J. Mulchaey (Carnegie Observatories director), L. Infante (Las Campanas Observatory director), and the entire Las Campanas staff for their dedication, professionalism, and excitement, all of which were critical in the observations used in this study. We thank I. Thompson and the Carnegie Observatory Time Allocation Committee for approving the Swope Supernova Survey and scheduling our program. We thank the University of Copenhagen, DARK Cosmology Centre, and the Niels Bohr International Academy for hosting D.A.C., R.J.F., A.M.-B., E.R.-R., and M.R.S. during part of this work. R.J.F., A.M.-B., and E.R.-R. were participating in the Kavli Summer Program in Astrophysics, “Astrophysics with gravitational wave detections.” This program was supported by the the Kavli Foundation, Danish National Research Foundation (DNRF), the Niels Bohr International Academy, and the DARK Cosmology Centre. The University of California, Santa Cruz (UCSC), group is supported in part by NSF grant AST–1518052, the Gordon & Betty Moore Foundation, the Heising-Simons Foundation, generous donations from many individuals through a UCSC Giving Day grant, and from fellowships from the Alfred P. Sloan Foundation (R.J.F.), the David and Lucile Packard Foundation (R.J.F. and E.R.-R.) and the Niels Bohr Professorship from the DNRF (E.R.-R.). A.M.B. acknowledges support from a UCMEXUS-CONACYT Doctoral Fellowship. D.K. is supported in part by a U.S. Department of Energy (DOE) Early Career award DE-SC0008067, a DOE Office of Nuclear Physics award DE-SC0017616, and a DOE SciDAC award DE-SC0018297 and by the Director, Office of Energy Research, Office of High Energy and Nuclear Physics, Divisions of Nuclear Physics, of the U.S. DOE under contract no.DE-AC02-05CH11231. M.R.D. is a Hubble and Carnegie-Dunlap Fellow. M.R.D. acknowledges support from the Dunlap Institute at the University of Toronto. B.F.M. is an unpaid visiting scientist at the University of Chicago and a consultant to the NASA-IPAC Extragalactic Database. J.X.P. is an affiliate member of the Institute for Physics and Mathematics of the Universe. B.J.S. acknowledges support from a Carnegie-Princeton Fellowship. M.R.D. and B.J.S. were partially supported by NASA through Hubble Fellowship grant HST-HF-51373.001 awarded by the Space Telescope Science Institute. This paper includes data gathered with the 6.5-m Magellan Telescopes located at Las Campanas Observatory, Chile. This research has made use of the NASA-IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. Based in part on observations made with the NASA-ESA Hubble Space Telescope, obtained from the Data Archive at the Space Telescope Science Institute (, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5–26555. These observations are associated with program GO–14840. All photometry and spectroscopy of SSS17a, kilonova models, and code used in our analysis will be available at
View Abstract


Navigate This Article