• ### A Decline in the X-ray through Radio Emission from GW170817 Continues to Support an Off-Axis Structured Jet(1805.02870)

May 8, 2018 astro-ph.HE
We present new observations of the binary neutron star merger GW170817 at $\Delta t\approx 220-260$ days post-merger, at radio (Karl G.~Jansky Very Large Array; VLA), X-ray ({\it Chandra X-ray Observatory}) and optical ({\it Hubble Space Telescope}; HST) wavelengths. These observations provide the first evidence for a turnover in the X-ray light curve, mirroring a decline in the radio emission at $\gtrsim5\sigma$ significance. The radio to X-ray spectral energy distribution exhibits no evolution into the declining phase. Our full multi-wavelength dataset is consistent with the predicted behavior of our previously-published models of a structured jet expanding into a low density circumbinary medium. If future observations continue to track our predictions, we expect the radio and X-ray emission to remain detectable until $\sim 1000$ days post-merger.
• ### A Precise Distance to the Host Galaxy of the Binary Neutron Star Merger GW170817 Using Surface Brightness Fluctuations(1801.06080)

Feb. 27, 2018 astro-ph.CO, astro-ph.GA
The joint detection of gravitational waves and electromagnetic radiation from the binary neutron star (BNS) merger GW170817 has provided unprecedented insight into a wide range of physical processes: heavy element synthesis via the $r$-process; the production of relativistic ejecta; the equation of state of neutron stars and the nature of the merger remnant; the binary coalescence timescale; and a measurement of the Hubble constant via the "standard siren" technique. In detail, all of these results depend on the distance to the host galaxy of the merger event, NGC4993. In this paper we measure the surface brightness fluctuation (SBF) distance to NGC4993 in the F110W and F160W passbands of the Wide Field Camera 3 Infrared Channel on the Hubble Space Telescope (HST). For the preferred F110W passband we derive a distance modulus of $m{-}M=33.05\pm0.08\pm0.10$ mag, or a linear distance $d=40.7\pm1.4\pm1.9$ Mpc (random and systematic errors, respectively); a virtually identical result is obtained from the F160W data. This is the most precise distance to NGC4993 available to date. Combining our distance measurement with the corrected recession velocity of NGC4993 implies a Hubble constant $H_0=71.9\pm 7.1~km~s^{-1}~Mpc^{-1}$. A comparison of our result to the GW-inferred value of $H_0$ indicates a binary orbital inclination of $i\,{\gtrsim}\,137~\deg$. The SBF technique can be applied to early-type host galaxies of BNS mergers to ${\sim\,}100$ Mpc with HST and possibly as far as ${\sim\,}300$ Mpc with the James Webb Space Telescope, thereby helping to break the inherent distance-inclination degeneracy of the GW signals at distances where many future BNS mergers are likely to be detected.
• ### Jets in Hydrogen-poor Super-luminous Supernovae: Constraints from a Comprehensive Analysis of Radio Observations(1711.03428)

Feb. 26, 2018 astro-ph.HE
The energy source powering the extreme optical luminosity of hydrogen-stripped Superluminous Supernovae (SLSNe-I) is not known, but recent studies have highlighted the case for a central engine. Radio and/or X-ray observations are best placed to track the fastest ejecta and probe the presence of outflows from a central engine. We compile all the published radio observations of SLSNe-I to date and present three new observations of two new SLSNe-I. None were detected. Through modeling the radio emission, we constrain the sub-parsec environments and possible outflows in SLSNe-I. In this sample we rule out on-axis collimated relativistic jets of the kind detected in Gamma-Ray Bursts (GRBs). We constrain off-axis jets with opening angles of 5\arcdeg\ (30\arcdeg) to energies of $\rm{E_k<4\times10^{50}\,erg}$ ($\rm{E_k<10^{50}\,erg}$) in environments shaped by progenitors with mass-loss rates of $\dot{M}<10^{-4}\,M_{\odot}\,{\rm yr}^{-1}$ ($\dot{M}<10^{-5}\,M_{\odot}\,{\rm yr}^{-1}$) for all off-axis angles, assuming fiducial values $\epsilon_e=0.1$ and $\epsilon_B=0.01$. The deepest limits rule out emission of the kind seen in faint un-collimated GRBs (with the exception of GRB\,060218), and from relativistic supernovae. Finally, for the closest SLSN-I SN 2017egm we constrained the energy of an uncollimated non-relativistic outflow like those observed in normal SNe to $E_{\rm k}\lesssim10^{48}$ erg.
• ### The Binary Neutron Star event LIGO/VIRGO GW170817 a hundred and sixty days after merger: synchrotron emission across the electromagnetic spectrum(1801.03531)

Feb. 25, 2018 astro-ph.HE
We report deep Chandra, HST and VLA observations of the binary neutron star event GW170817 at $t<160$ d 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 $\nu_c$ is above the X-ray band and the synchrotron frequency $\nu_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 $N(\gamma)\propto \gamma^{-p}$ accelerated by a shock launched by a NS-NS merger to date. We find $p=2.17\pm0.01$, which indicates that radiation from ejecta with $\Gamma\sim3-10$ dominates the observed emission. While constraining the nature of the emission process, these observations do \emph{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, \emph{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 SGRB directed away from our line of sight. Observations at $t\le 200$ days are unlikely to settle the debate as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.
• ### On the nature of Hydrogen-rich Superluminous Supernovae(1604.01226)

We present two hydrogen-rich superluminous supernovae (SLSNe), namely SN2013hx and PS15br. These objects, together with SN2008es are the only SLSNe showing a distinct, broad Halpha feature during the photospheric phase and also do not show any sign of strong interaction between fast-moving ejecta and circumstellar shells in their early spectra. Despite PS15br peak luminosity is fainter than the other two objects, the spectrophotometric evolution is similar to SN2013hx and different than any other supernova in a similar luminosity space. We group all of them as SLSNe II and hence distinct from the known class of SLSN IIn. Both transients show a strong, multi-component Halpha emission after 200 days past maximum which we interpret as an indication of interaction of the ejecta with an asymmetric, clumpy circumstellar material. The spectra and photometric evolution of the two objects are similar to type II supernovae, although they have much higher luminosity and evolve on slower timescales. This is qualitatively similar to how SLSNe I compare with normal type Ic in that the former are brighter and evolve more slowly. We apply a magnetar and an interaction semi-analytical codes to fit the light curves of our two objects and SN2008es. The overall observational dataset would tend to favour the magnetar, or central engine, model as the source of the peak luminosity although the clear signature of late-time interaction indicates that interaction can play a role in the luminosity evolution of SLSNe II at some phases.
• ### Improved constraints on H0 from a combined analysis of gravitational-wave and electromagnetic emission from GW170817(1710.06426)

Oct. 19, 2017 astro-ph.CO
The luminosity distance measurement of GW170817 derived from GW analysis in Abbott et al. 2017 (here, A17:H0) is highly correlated with the measured inclination of the NS-NS system. To improve the precision of the distance measurement, we attempt to constrain the inclination by modeling the broad-band X-ray-to-radio emission from GW170817, which is dominated by the interaction of the jet with the environment. We update our previous analysis and we consider the radio and X-ray data obtained at $t<40$ days since merger. We find that the afterglow emission from GW170817 is consistent with an off-axis relativistic jet with energy $10^{48}\,\rm{erg}<E_{k}\le 3\times 10^{50} \,\rm{erg}$ propagating into an environment with density $n\sim10^{-2}-10^{-4} \,\rm{cm^{-3}}$, with preference for wider jets (opening angle $\theta_j=15$ deg). For these jets, our modeling indicates an off-axis angle $\theta_{\rm obs}\sim25-50$ deg. We combine our constraints on $\theta_{\rm obs}$ with the joint distance-inclination constraint from LIGO. Using the same $\sim 170$ km/sec peculiar velocity uncertainty assumed in A17:H0 but with an inclination constraint from the afterglow data, we get a value of $H_0=$$74.0 \pm \frac{11.5}{7.5} \mbox{km/s/Mpc}, which is higher than the value of H_0=$$70.0 \pm \frac{12.0}{8.0}$ $\mbox{km/s/Mpc}$ found in A17:H0. Further, using a more realistic peculiar velocity uncertainty of 250 km/sec derived from previous work, we find $H_0=$$75.5 \pm \frac{11.6}{9.6}$ km/s/Mpc for H0 from this system. We note that this is in modestly better agreement with the local distance ladder than the Planck CMB, though a significant such discrimination will require $\sim 50$ such events. Future measurements at $t>100$ days of the X-ray and radio emission will lead to tighter constraints.
• ### How Many Kilonovae Can Be Found in Past, Present, and Future Survey Datasets?(1710.05845)

Oct. 16, 2017 astro-ph.IM
The discovery of a kilonova (KN) associated with the Advanced LIGO (aLIGO)/Virgo event GW170817 opens up new avenues of multi-messenger astrophysics. Here, using realistic simulations, we provide estimates of the number of KNe that could be found in data from past, present and future surveys without a gravitational-wave trigger. For the simulation, we construct a spectral time-series model based on the DES-GW multi-band light-curve from the single known KN event, and we use an average of BNS rates from past studies of $10^3 \rm{Gpc}^{-3}/\rm{year}$, consistent with the $1$ event found so far. Examining past and current datasets from transient surveys, the number of KNe we expect to find for ASAS-SN, SDSS, PS1, SNLS, DES, and SMT is between 0 and $0.3$. We predict the number of detections per future survey to be: 8.3 from ATLAS, 10.6 from ZTF, 5.5/69 from LSST (the Deep Drilling / Wide Fast Deep), and 16.0 from WFIRST. The maximum redshift of KNe discovered for each survey is z = 0.8 for WFIRST, z = 0.25 for LSST and z = 0.04 for ZTF and ATLAS. For the LSST survey, we also provide contamination estimates from Type Ia and Core-collapse supernovae: after light-curve and template-matching requirements, we estimate a background of just 2 events. More broadly, we stress that future transient surveys should consider how to optimize their search strategies to improve their detection efficiency, and to consider similar analyses for GW follow-up programs.
• ### The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. V. Rising X-ray Emission from an Off-Axis Jet(1710.05431)

Oct. 16, 2017 astro-ph.HE
We report the discovery of rising X-ray emission from the binary neutron star (BNS) merger event GW170817. This is the first detection of X-ray emission from a gravitational-wave 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.2x10^38 erg/s (isotropic-equivalent). Continued monitoring revealed the presence of an X-ray source that brightened with time, reaching L_x\sim 9x10^39 erg/s 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 broad-band X-ray to radio observations are consistent with emission from a relativistic jet with kinetic energy E_k~10^49-10^50 erg, viewed off-axis with theta_obs~ 20-40 deg. Our models favor a circumbinary density n~ 0.0001-0.01 cm-3, depending on the value of the microphysical parameter epsilon_B=10^{-4}-10^{-2}. A central-engine origin of the X-ray emission is unlikely. Future X-ray observations at $t\gtrsim 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 theta_obs are testable with gravitational wave information on GW170817 from Advanced LIGO/Virgo on the binary inclination.
• ### The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VIII. A Comparison to Cosmological Short-duration Gamma-ray Bursts(1710.05438)

Oct. 16, 2017 astro-ph.HE
We present a comprehensive comparison of the properties of the radio through X-ray counterpart of GW170817 and the properties of short-duration gamma-ray bursts (GRBs). For this effort, we utilize a sample of 36 short GRBs spanning a redshift range of $z \approx 0.12-2.6$ discovered over 2004-2017. We find that the counterpart to GW170817 has an isotropic-equivalent luminosity that is $\approx 3000$ times less than the median value of on-axis short GRB X-ray afterglows, and $\gtrsim10^{4}$ times less than that for detected short GRB radio afterglows. Moreover, the allowed jet energies and particle densities inferred from the radio and X-ray counterparts to GW170817 and on-axis short GRB afterglows are remarkably similar, suggesting that viewing angle effects are the dominant, and perhaps only, difference in their observed radio and X-ray behavior. From comparison to previous claimed kilonovae following short GRBs, we find that the optical and near-IR counterpart to GW170817 is comparatively under-luminous by a factor of $\approx 3-5$, indicating a range of kilonova luminosities and timescales. A comparison of the optical limits following short GRBs on $\lesssim 1$ day timescales also rules out a "blue" kilonova of comparable optical isotropic-equivalent luminosity in one previous short GRB. Finally, we investigate the host galaxy of GW170817, NGC4993, in the context of short GRB host galaxy stellar population properties. We find that NGC4993 is superlative in terms of its large luminosity, old stellar population age, and low star formation rate compared to previous short GRB hosts. Additional events within the Advanced LIGO/VIRGO volume will be crucial in delineating the properties of the host galaxies of NS-NS mergers, and connecting them to their cosmological counterparts.
• ### The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. IV. Detection of Near-infrared Signatures of r-process Nucleosynthesis with Gemini-South(1710.05454)

We present a near-infrared spectral sequence of the electromagnetic counterpart to the binary neutron star merger GW170817 detected by Advanced LIGO/Virgo. Our dataset comprises seven epochs of J+H spectra taken with FLAMINGOS-2 on Gemini-South between 1.5 and 10.5 days after the merger. In the initial epoch, the spectrum is dominated by a smooth blue continuum due to a high-velocity, lanthanide-poor blue kilonova component. Starting the following night, all subsequent spectra instead show features that are similar to those predicted in model spectra of material with a high concentration of lanthanides, including spectral peaks near 1.07 and 1.55 microns. Our fiducial model with 0.04 M_sun of ejecta, an ejection velocity of v=0.1c, and a lanthanide concentration of X_lan=1e-2 provides a good match to the spectra taken in the first five days, although it over-predicts the late-time fluxes. We also explore models with multiple fitting components, in each case finding that a significant abundance of lanthanide elements is necessary to match the broad spectral peaks that we observe starting at 2.5 d after the merger. These data provide direct evidence that binary neutron star mergers are significant production sites of even the heaviest r-process elements.
• ### The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. III. Optical and UV Spectra of a Blue Kilonova From Fast Polar Ejecta(1710.05456)

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 SOAR and Magellan telescopes; the UV spectrum was obtained with the \textit{Hubble Space Telescope} at 5.5 days. Our data reveal a rapidly-fading blue component ($T\approx5500$ K at 1.5 days) that quickly reddens; spectra later than $\gtrsim 4.5$ days peak beyond the optical regime. The spectra are mostly featureless, although we identify a possible weak emission line at $\sim 7900$ \AA\ at $t\lesssim 4.5$ days. The colours, 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\lesssim 140$. This indicates a sight-line within $\theta_{\rm obs}\lesssim 45^{\circ}$ of the orbital axis. Comparison to models suggests $\sim0.03$ M$_\odot$ of blue ejecta, with a velocity of $\sim 0.3c$. The required lanthanide fraction is $\sim 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 $\lesssim 12$ km. This mass also supports the idea that neutron star mergers are a major contributor to $r$-process nucleosynthesis.
• ### The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VI. Radio Constraints on a Relativistic Jet and Predictions for Late-Time Emission from the Kilonova Ejecta(1710.05457)

Oct. 16, 2017 astro-ph.HE
We present Very Large Array (VLA) and Atacama Large Millimeter/sub-millimeter Array ALMA radio observations of GW\,170817, 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$ hours 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 d. 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 $\gtrsim 10^{48}$ erg. For fiducial SGRB parameters, our limits require an observer viewer angle of $\gtrsim 20^{\circ}$. The radio and X-ray data can be jointly explained as the afterglow emission from an SGRB with a jet energy of $\sim 10^{49}-10^{50}$ erg that exploded in a uniform density environment with $n\sim 10^{-4}-10^{-2}$ cm$^{-3}$, viewed at an angle of $\sim 20^{\circ}-40^{\circ}$ 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 $\sim 5-10$ years that will remain detectable for decades with next-generation radio facilities, making GW\,170817 a compelling target for long-term radio monitoring.
• ### The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VII. Properties of the Host Galaxy and Constraints on the Merger Timescale(1710.05458)

Oct. 16, 2017 astro-ph.HE
We present the properties of NGC 4993, the host galaxy of GW170817, the first gravitational wave (GW) event from the merger of a binary neutron star (BNS) system and the first with an electromagnetic (EM) counterpart. We use both archival photometry and new optical/near-IR imaging and spectroscopy, together with stellar population synthesis models to infer the global properties of the host galaxy. We infer a star formation history peaked at $\gtrsim 10$ Gyr ago, with subsequent exponential decline leading to a low current star formation rate of 0.01 M$_{\odot}$ yr$^{-1}$, which we convert into a binary merger timescale probability distribution. We find a median merger timescale of $11.2^{+0.7}_{-1.4}$ Gyr, with a 90% confidence range of $6.8-13.6$ Gyr. This in turn indicates an initial binary separation of $\approx 4.5$ R$_{\odot}$, comparable to the inferred values for Galactic BNS systems. We also use new and archival $Hubble$ $Space$ $Telescope$ images to measure a projected offset of the optical counterpart of $2.1$ kpc (0.64$r_{e}$) from the center of NGC 4993 and to place a limit of $M_{r} \gtrsim -7.2$ mag on any pre-existing emission, which rules out the brighter half of the globular cluster luminosity function. Finally, the age and offset of the system indicates it experienced a modest natal kick with an upper limit of $\sim 200$ km s$^{-1}$. Future GW$-$EM observations of BNS mergers will enable measurement of their population delay time distribution, which will directly inform their viability as the dominant source of $r$-process enrichment in the Universe.
• ### The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. I. Dark Energy Camera Discovery of the Optical Counterpart(1710.05459)

Oct. 16, 2017 astro-ph.HE
We present the Dark Energy Camera (DECam) discovery of the optical counterpart of the first binary neutron star merger detected through gravitational wave emission, GW170817. Our observations commenced 10.5 hours post-merger, as soon as the localization region became accessible from Chile. We imaged 70 deg$^2$ in the $i$ and $z$ bands, covering 93\% of the initial integrated localization probability, to a depth necessary to identify likely optical counterparts (e.g., a kilonova). At 11.4 hours post-merger we detected a bright optical transient located $10.6''$ from the nucleus of NGC\,4993 at redshift $z=0.0098$, consistent (for $H_0 = 70$\, km s$^{-1}$ Mpc$^{-1}$) with the distance of $40 \pm 8$\, Mpc reported by the LIGO Scientific Collaboration and the Virgo Collaboration (LVC). At detection the transient had magnitudes $i\approx 17.30$ and $z\approx 17.45$, and thus an absolute magnitude of $M_i = -15.7$, in the luminosity range expected for a kilonova. We identified 1,500 potential transient candidates. Applying simple selection criteria aimed at rejecting background events such as supernovae, we find the transient associated with NGC\,4993 as the only remaining plausible counterpart, and reject chance coincidence at the 99.5\% confidence level. We therefore conclude that the optical counterpart we have identified near NGC\,4993 is associated with GW170817. This discovery ushers in the era of multi-messenger astronomy with gravitational waves, and demonstrates the power of DECam to identify the optical counterparts of gravitational-wave sources.
• ### The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. II. UV, Optical, and Near-IR Light Curves and Comparison to Kilonova Models(1710.05840)

Oct. 16, 2017 astro-ph.HE
We present UV, optical, and NIR photometry of the first electromagnetic counterpart to a gravitational wave source from Advanced LIGO/Virgo, the binary neutron star merger GW170817. Our data set extends from the discovery of the optical counterpart at $0.47$ days to $18.5$ days post-merger, and includes observations with the Dark Energy Camera (DECam), Gemini-South/FLAMINGOS-2 (GS/F2), and the {\it Hubble Space Telescope} ({\it HST}). The spectral energy distribution (SED) inferred from this photometry at $0.6$ days is well described by a blackbody model with $T\approx 8300$ K, a radius of $R\approx 4.5\times 10^{14}$ cm (corresponding to an expansion velocity of $v\approx 0.3c$), and a bolometric luminosity of $L_{\rm bol}\approx 5\times10^{41}$ erg s$^{-1}$. At $1.5$ days we find a multi-component SED across the optical and NIR, and subsequently we observe rapid fading in the UV and blue optical bands and significant reddening of the optical/NIR colors. Modeling the entire data set we find that models with heating from radioactive decay of $^{56}$Ni, or those with only a single component of opacity from $r$-process elements, fail to capture the rapid optical decline and red optical/NIR colors. Instead, models with two components consistent with lanthanide-poor and lanthanide-rich ejecta provide a good fit to the data, the resulting "blue" component has $M_\mathrm{ej}^\mathrm{blue}\approx 0.01$ M$_\odot$ and $v_\mathrm{ej}^\mathrm{blue}\approx 0.3$c, and the "red" component has $M_\mathrm{ej}^\mathrm{red}\approx 0.04$ M$_\odot$ and $v_\mathrm{ej}^\mathrm{red}\approx 0.1$c. These ejecta masses are broadly consistent with the estimated $r$-process production rate required to explain the Milky Way $r$-process abundances, providing the first evidence that BNS mergers can be a dominant site of $r$-process enrichment.
• The detection of GW170817 in both gravitational waves and electromagnetic waves heralds the age of gravitational-wave multi-messenger astronomy. On 17 August 2017 the Advanced LIGO and Virgo detectors observed GW170817, a strong signal from the merger of a binary neutron-star system. Less than 2 seconds after the merger, a gamma-ray burst (GRB 170817A) was detected within a region of the sky consistent with the LIGO-Virgo-derived location of the gravitational-wave source. This sky region was subsequently observed by optical astronomy facilities, resulting in the identification of an optical transient signal within $\sim 10$ arcsec of the galaxy NGC 4993. These multi-messenger observations allow us to use GW170817 as a standard siren, the gravitational-wave analog of an astronomical standard candle, to measure the Hubble constant. This quantity, which represents the local expansion rate of the Universe, sets the overall scale of the Universe and is of fundamental importance to cosmology. Our measurement combines the distance to the source inferred purely from the gravitational-wave signal with the recession velocity inferred from measurements of the redshift using electromagnetic data. This approach does not require any form of cosmic "distance ladder;" the gravitational wave analysis can be used to estimate the luminosity distance out to cosmological scales directly, without the use of intermediate astronomical distance measurements. We determine the Hubble constant to be $70.0^{+12.0}_{-8.0} \, \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ (maximum a posteriori and 68% credible interval). This is consistent with existing measurements, while being completely independent of them. Additional standard-siren measurements from future gravitational-wave sources will provide precision constraints of this important cosmological parameter.
• ### Hydrogen-rich supernovae beyond the neutrino-driven core-collapse paradigm(1709.10475)

Sept. 29, 2017 astro-ph.SR
We present our study of OGLE-2014-SN-073, one of the brightest Type II SN ever discovered, with an unusually broad lightcurve combined with high ejecta velocities. From our hydrodynamical modelling we infer a remarkable ejecta mass of $60^{+42}_{-16}$~M$_\odot$, and a relatively high explosion energy of $12.4^{+13.0}_{-5.9} \times10^{51}$~erg. We show that this object belongs, with a very small number of other hydrogen-rich SNe, to an energy regime that is not explained by standard core-collapse (CC) neutrino-driven explosions. We compare the quantities inferred by the hydrodynamical modelling with the expectations of various exploding scenarios, trying to explain the high energy and luminosity released. We find some qualitative similarities with pair-instabilities SNe, although a prompt injection of energy by a magnetar seems also a viable alternative to explain such extreme event.
• ### Empirical constraints on the origin of fast radio bursts: volumetric rates and host galaxy demographics as a test of millisecond magnetar connection(1704.00022)

July 11, 2017 astro-ph.CO, astro-ph.HE
The localization of the repeating FRB 121102 to a low-metallicity dwarf galaxy at $z=0.193$, and its association with a quiescent radio source, suggests the possibility that FRBs originate from magnetars, formed by the unusual supernovae in such galaxies. We investigate this via a comparison of magnetar birth rates, the FRB volumetric rate, and host galaxy demographics. We calculate average volumetric rates of possible millisecond magnetar production channels such as superluminous supernovae (SLSNe), long and short gamma-ray bursts (GRBs), and general magnetar production via core-collapse supernovae. For each channel we also explore the expected host galaxy demographics using their known properties. We determine for the first time the number density of FRB emitters (the product of their volumetric birthrate and lifetime), $R_{\rm FRB}\tau\approx 10^4$Gpc$^{-3}$, assuming that FRBs are predominantly emitted from repetitive sources similar to FRB 121102 and adopting a beaming factor of 0.1. By comparing rates we find that production via rare channels (SLSNe, GRBs) implies a typical FRB lifetime of $\approx$30-300 yr, in good agreement with other lines of argument. The total energy emitted over this time is consistent with the available energy stored in the magnetic field. On the other hand, any relation to magnetars produced via normal core-collapse supernovae leads to a very short lifetime of $\approx$0.5yr, in conflict with both theory and observation. We demonstrate that due to the diverse host galaxy distributions of the different progenitor channels, many possible sources of FRB birth can be ruled out with $\lesssim 10$ host galaxy identifications. Conversely, targeted searches of galaxies that have previously hosted decades-old SLSNe and GRBs may be a fruitful strategy for discovering new FRBs and related quiescent radio sources, and determining the nature of their progenitors.
• ### Superluminous supernova progenitors have a half-solar metallicity threshold(1605.04925)

May 23, 2017 astro-ph.GA
Host galaxy properties provide strong constraints on the stellar progenitors of superluminous supernovae. By comparing a sample of 19 low-redshift (z < 0.3) superluminous supernova hosts to galaxy populations in the local Universe, we show that sub-solar metallicities seem to be a requirement. All superluminous supernovae in hosts with high measured gas-phase metallicities are found to explode at large galactocentric radii, indicating that the metallicity at the explosion site is likely lower than the integrated host value. We found that superluminous supernovae hosts do not always have star-formation rates higher than typical star-forming galaxies of the same mass. However, we confirm that high absolute specific star-formation rates are a feature of superluminous supernova host galaxies, but interpret this as simply a consequence of the anti-correlation between gas-phase metallicity and specific star-formation rate and the requirement of on-going star formation to produce young, massive stars greater than ~ 10-20 M_sol. Based on our sample, we propose an upper limit of ~ 0.5 Z_sol for forming superluminous supernova progenitors (assuming an N2 metallicity diagnostic and a solar oxygen abundance of 8.69). Finally, we show that if magnetar powering is the source of the extreme luminosity then the required initial spins appear to be correlated with metallicity of the host galaxy. This correlation needs further work, but if it holds it is a powerful link between the supernova parameters and nature of the progenitor population.
• ### Complexity in the light curves and spectra of slow-evolving superluminous supernovae(1701.00941)

April 25, 2017 astro-ph.CO, astro-ph.HE
A small group of the newly discovered superluminous supernovae show broad and slowly evolving light curves. Here we present extensive observational data for the slow-evolving superluminous supernova LSQ14an, which brings this group of transients to four in total in the low redshift Universe (z$<$0.2; SN 2007bi, PTF12dam, SN 2015bn). We particularly focus on the optical and near-infrared evolution during the period from 50 days up to 400 days from peak, showing that they are all fairly similar in their light curve and spectral evolution. LSQ14an shows broad, blue-shifted [O III] $\lambda\lambda$4959, 5007 lines, as well as a blue-shifted [O II] $\lambda\lambda$7320, 7330 and [Ca II] $\lambda\lambda$7291, 7323. Furthermore, the sample of these four objects shows common features. Semi-forbidden and forbidden emission lines appear surprisingly early at 50-70 days and remain visible with almost no variation up to 400 days. The spectra remain blue out to 400 days. There are small, but discernible light curve fluctuations in all of them. The light curve of each shows a faster decline than $^{56}$Co after 150 days and it further steepens after 300 days. We also expand our analysis presenting X-ray limits for LSQ14an and SN2015bn and discuss their diagnostic power in interpreting their common features. These features are quite distinct from the faster evolving superluminous supernovae and are not easily explained in terms of only a variation in ejecta mass. While a central engine is still the most likely luminosity source, it appears that the ejecta structure is complex, with multiple emitting zones and at least some interaction between the expanding ejecta and surrounding material.
• ### Results from a systematic survey of X-ray emission from Hydrogen-poor Superluminous Supernovae(1704.05865)

April 19, 2017 astro-ph.HE
We present the results from a sensitive X-ray survey of 26 nearby hydrogen-poor superluminous supernovae (SLSNe-I) with Swift, Chandra and XMM. This dataset constrains the SLSN evolution from a few days until ~2000 days after explosion, reaching a luminosity L_x~10^40 erg/s and revealing the presence of significant X-ray emission at the location of PTF12dam. No SLSN-I is detected above L_x~10^41 erg/s, suggesting that the luminous X-ray emission L_x~10^45 erg/s associated with SCP60F6 is not common among SLSNe-I. We constrain the presence of off-axis GRB jets, ionization breakouts from magnetar central engines and the density in the sub-pc environments of SLSNe-I through Inverse Compton emission. The deepest limits rule out the weakest uncollimated GRB outflows, suggesting that IF the similarity of SLSNe-I with GRB/SNe extends to their fastest ejecta, then SLSNe-I are either powered by energetic jets pointed far away from our line of sight theta>30 deg, or harbor failed jets that do not successfully break through the stellar envelope. Furthermore, IF a magnetar central engine is responsible for the exceptional luminosity of SLSNe-I, our X-ray analysis favors large magnetic fields B>2x10^(14) G and ejecta masses M_ej>3 Msun in agreement with optical/UV studies. Finally, we constrain the pre-explosion mass-loss rate of stellar progenitors of SLSNe-I. For PTF12dam we infer Mdot<2x10^(-5) Msun/yr, suggesting that the SN shock interaction with the CSM is unlikely to supply the main source of energy powering the optical transient and that some SLSN-I progenitors end their life as compact stars surrounded by a low-density medium similar to long GRBs and Type Ib/c SNe.
• ### The evolution of superluminous supernova LSQ14mo and its interacting host galaxy system(1611.09910)

We present and analyse an extensive dataset of the superluminous supernova (SLSN) LSQ14mo (z = 0.256), consisting of a multi-colour lightcurve from -30 d to +70 d in the rest-frame and a series of 6 spectra from PESSTO covering -7 d to +50 d. This is among the densest spectroscopic coverage, and best-constrained rising lightcurve, for a fast-declining hydrogen-poor SLSN. The bolometric lightcurve can be reproduced with a millisecond magnetar model with ~ 4 M_sol ejecta mass, and the temperature and velocity evolution is also suggestive of a magnetar as the power source. Spectral modelling indicates that the SN ejected ~ 6 M_sol of CO-rich material with a kinetic energy of ~ 7 x 10^51 erg, and suggests a partially thermalised additional source of luminosity between -2 d and +22 d. This may be due to interaction with a shell of material originating from pre-explosion mass loss. We further present a detailed analysis of the host galaxy system of LSQ14mo. PESSTO and GROND imaging show three spatially resolved bright regions, and we used the VLT and FORS2 to obtain a deep (five-hour exposure) spectra of the SN position and the three star-forming regions, which are at a similar redshift. The FORS spectrum at +300 days shows no trace of SN emission lines and we place limits on the strength of [O I] from comparisons with other Ic SNe. The deep spectra provides a unique chance to investigate spatial variations in the host star-formation activity and metallicity. The specific star-formation rate is similar in all three components, as is the presence of a young stellar population. However, the position of LSQ14mo exhibits a lower metallicity, with 12 + log(O/H) = 8.2 in both the R23 and N2 scales (corresponding to ~ 0.3 Z_sol). We propose that the three bright regions in the host system are interacting, which thus triggers star-formation and forms young stellar populations.
• ### Long-duration superluminous supernovae at late times(1608.02994)

Jan. 25, 2017 astro-ph.SR, astro-ph.HE
Nebular-phase observations and spectral models of Type Ic superluminous supernovae are presented. LSQ14an and SN 2015bn both display late-time spectra similar to galaxy-subtracted spectra of SN 2007bi, and the class shows strong similarity with broad-lined Type Ic SNe such as SN 1998bw. Near-infrared observations of SN 2015bn show a strong Ca II triplet, O I 9263, O I 1.13 um and Mg I 1.50 um, but no distinct He, Si, or S emission. The high Ca II NIR/[Ca II] 7291,7323 ratio of ~2 indicates a high electron density of n_e >~ 10^8 cm^{-3}. Spectral models of oxygen-zone emission are investigated to put constraints on the emitting region. Models require M(O-zone) >~ 10 Msun to produce enough [O I] 6300,6364 luminosity, irrespective of the powering situation and the density. The high oxygen-zone mass, supported by high estimated magnesium masses, points to explosions of massive CO cores, requiring M_ZAMS >~ 40 Msun. Collisions of pair-instability pulsations do not provide enough mass to account for the emission. [O II] and [O III] lines emerge naturally in many models, which strengthens the identification of broad [O II] 7320,7330, [O III] 4363, and [O III] 4959,5007 in some spectra. A small filling factor f <~ 0.01 for the O/Mg zone is needed to produce enough luminosity in Mg I] 4571, Mg I 1.504 um, and O I recombination lines, which shows that the ejecta is clumped. We review the constraints from the nebular spectral modelling in the context of the various scenarios proposed for superluminous supernovae.
• ### An ultraviolet excess in the superluminous supernova Gaia16apd reveals a powerful central engine(1611.06993)

Since the discovery of superluminous supernovae (SLSNe) in the last decade, it has been known that these events exhibit bluer spectral energy distributions than other supernova subtypes, with significant output in the ultraviolet. However, the event Gaia16apd seems to outshine even the other SLSNe at rest-frame wavelengths below $\sim 3000$ \AA. Yan et al (2016) have recently presented HST UV spectra and attributed the UV flux to low metallicity and hence reduced line blanketing. Here we present UV and optical light curves over a longer baseline in time, revealing a rapid decline at UV wavelengths despite a typical optical evolution. Combining the published UV spectra with our own optical data, we demonstrate that Gaia16apd has a much hotter continuum than virtually any SLSN at maximum light, but it cools rapidly thereafter and is indistinguishable from the others by $\sim 10$-15 days after peak. Comparing the equivalent widths of UV absorption lines with those of other events, we show that the excess UV continuum is a result of a more powerful central power source, rather than a lack of UV absorption relative to other SLSNe or an additional component from interaction with the surrounding medium. These findings strongly support the central-engine hypothesis for hydrogen-poor SLSNe. An explosion ejecting $M_{\rm ej} = 4 (0.2/\kappa)$ M$_\odot$, where $\kappa$ is the opacity in cm$^2$g$^{-1}$, and forming a magnetar with spin period $P=2$ ms, and $B=2\times10^{14}$ G (lower than other SLSNe with comparable rise-times) can consistently explain the light curve evolution and high temperature at peak. The host metallicity, $Z=0.18$ Z$_\odot$, is comparable to other SLSNe.
• ### The Superluminous Transient ASASSN-15lh as a Tidal Disruption Event from a Kerr Black Hole(1609.02927)

When a star passes within the tidal radius of a supermassive black hole, it will be torn apart. For a star with the mass of the Sun ($M_\odot$) and a non-spinning black hole with a mass $<10^8 M_\odot$, the tidal radius lies outside the black hole event horizon and the disruption results in a luminous flare. Here we report observations over a period of 10 months of a transient, hitherto interpreted as a superluminous supernova. Our data show that the transient rebrightened substantially in the ultraviolet and that the spectrum went through three different spectroscopic phases without ever becoming nebular. Our observations are more consistent with a tidal disruption event than a superluminous supernova because of the temperature evolution, the presence of highly ionised CNO gas in the line of sight and our improved localisation of the transient in the nucleus of a passive galaxy, where the presence of massive stars is highly unlikely. While the supermassive black hole has a mass $> 10^8 M_\odot$, a star with the same mass as the Sun could be disrupted outside the event horizon if the black hole were spinning rapidly. The rapid spin and high black hole mass can explain the high luminosity of this event.