• ### 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.
• ### Hydrogen-Poor Superluminous Supernovae from the Pan-STARRS1 Medium Deep Survey(1708.01619)

Jan. 12, 2018 astro-ph.HE
We present light curves and classification spectra of 17 hydrogen-poor superluminous supernovae (SLSNe) from the Pan-STARRS1 Medium Deep Survey (PS1 MDS). Our sample contains all objects from the PS1 MDS sample with spectroscopic classification that are similar to either of the prototypes SN2005ap or SN2007bi, without an explicit limit on luminosity. With a redshift range $0.3 < z < 1.6$, PS1MDS is the first SLSN sample primarily probing the high-redshift population; our multi-filter PS1 light curves probe the rest-frame UV emission, and hence the peak of the spectral energy distribution. We measure the temperature evolution and construct bolometric light curves, and find peak luminosities of $(0.5-5) \times 10^{44}$ erg s$^{-1}$ and lower limits on the total radiated energies of $(0.3-2) \times 10^{51}$ erg. The light curve shapes are diverse, with both rise- and decline times spanning a factor of $\sim 5$, and several examples of double-peaked light curves. When correcting for the flux-limited nature of our survey, we find a median peak luminosity at 4000 {\AA} of $M_{\rm 4000} = -21.1$ mag, and a spread of $\sigma = 0.7$ mag.
• ### 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.
• ### 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.
• ### An Empirical Study of Contamination in Deep, Rapid, and Wide-Field Optical Follow-Up of Gravitational Wave Events(1710.02144)

Oct. 9, 2017 astro-ph.HE
We present an empirical study of contamination in deep, rapid, and wide-field optical follow-up searches of GW sources from aLIGO. We utilize dedicated observations during four nights of imaging with DECam. Our search covered $\sim56$ deg$^2$, with two visits per night separated by $\approx 3$~hours, in $i$- and $z$-band, followed by an additional set of $griz$ images three weeks later to serve as reference images for subtraction, and for the purpose of identifying galaxy and stellar counterparts for any transient sources. We achieve $5\sigma$ point-source limiting magnitudes of $i \approx 23.5$ and $z \approx 22.4$ mag in the coadded single-epoch images. We conduct a search for transient objects that can mimic the $i-z$ color behavior of both red ($i-z > 0.5$~mag) and blue ($i-z < 0$~mag) kilonova emission, finding 11 and 10 contaminants, respectively. Independent of color, we identify 48 transients of interest. Additionally, we leverage the rapid cadence of our observations to search for sources with characteristic timescales of $\approx1$ day and $\approx3$ hours, finding no potential contaminants. We assess the efficiency of our pipeline and search methodology with injected point sources, finding that we are 90\% (60\%) efficient when searching for red (blue) kilonova-like sources to a limiting magnitude of $i \lesssim 22.5$ mag. Applying these efficiencies, we derive sky rates for kilonova contaminants in the red and blue regimes of $\mathcal{R}_{\rm red} \approx 0.16$ deg$^{-2}$ and $\mathcal{R}_{\rm blue} \approx 0.80$ deg$^{-2}$. The total contamination rate, independent of color, is $\mathcal{R}_{\rm all} \approx 1.79$ deg$^{-2}$. We compare our derived results to optical follow-up searches of the GW events GW150914 and GW151226 and comment on the outlook for GW follow-up searches as additional GW detectors (e.g., KAGRA, LIGO India) come online in the next decade.
• ### A Reverse Shock and Unusual Radio Properties in GRB 160625B(1705.08455)

Sept. 6, 2017 astro-ph.HE
We present multi-wavelength observations and modeling of the exceptionally bright long $\gamma$-ray burst GRB 160625B. The optical and X-ray data are well-fit by synchrotron emission from a collimated blastwave with an opening angle of $\theta_j\approx 3.6^\circ$ and kinetic energy of $E_K\approx 2\times10^{51}$ erg, propagating into a low density ($n\approx 5\times10^{-5}$ cm$^{-3}$) medium with a uniform profile. The forward shock is sub-dominant in the radio band; instead, the radio emission is dominated by two additional components. The first component is consistent with emission from a reverse shock, indicating an initial Lorentz factor of $\Gamma_0\gtrsim 100$ and an ejecta magnetization of $R_B\approx 1-100$. The second component exhibits peculiar spectral and temporal evolution and is most likely the result of scattering of the radio emission by the turbulent Milky Way interstellar medium (ISM). Such scattering is expected in any sufficiently compact extragalactic source and has been seen in GRBs before, but the large amplitude and long duration of the variability seen here are qualitatively more similar to extreme scattering events previously observed in quasars, rather than normal interstellar scintillation effects. High-cadence, broadband radio observations of future GRBs are needed to fully characterize such effects, which can sensitively probe the properties of the ISM and must be taken into account before variability intrinsic to the GRB can be interpreted correctly.
• ### The properties of GRB 120923A at a spectroscopic redshift of z=7.8(1703.09052)

March 27, 2017 astro-ph.CO, astro-ph.HE
Gamma-ray bursts (GRBs) are powerful probes of early stars and galaxies, during and potentially even before the era of reionization. Although the number of GRBs identified at z>6 remains small, they provide a unique window on typical star-forming galaxies at that time, and thus are complementary to deep field observations. We report the identification of the optical drop-out afterglow of Swift GRB 120923A in near-infrared Gemini-North imaging, and derive a redshift of z=7.84_{-0.12}^{+0.06} from VLT/X-shooter spectroscopy. At this redshift the peak 15-150 keV luminosity of the burst was 3.2x10^52 erg/s, and in fact the burst was close to the Swift/BAT detection threshold. The X-ray and near-infrared afterglow were also faint, and in this sense it was a rather typical long-duration GRB in terms of rest-frame luminosity. We present ground- and space-based follow-up observations spanning from X-ray to radio, and find that a standard external shock model with a constant-density circumburst environment with density, n~4x10^-2 cm^-3 gives a good fit to the data. The near-infrared light curve exhibits a sharp break at t~3.4 days in the observer frame, which if interpreted as being due to a jet corresponds to an opening angle of ~5 degrees. The beaming corrected gamma-ray energy is then E_gamma~2x10^50 erg, while the beaming-corrected kinetic energy is lower, E_K~10^49 erg, suggesting that GRB 120923A was a comparatively low kinetic energy event. We discuss the implications of this event for our understanding of the high-redshift population of GRBs and their identification.
• ### PS1-14bj: A Hydrogen-Poor Superluminous Supernova With a Long Rise and Slow Decay(1605.05235)

Oct. 10, 2016 astro-ph.HE
We present photometry and spectroscopy of PS1-14bj, a hydrogen-poor superluminous supernova (SLSN) at redshift $z=0.5215$ discovered in the last months of the Pan-STARRS1 Medium Deep Survey. PS1-14bj stands out by its extremely slow evolution, with an observed rise of $\gtrsim 125$ rest-frame days, and exponential decline out to $\sim 250$ days past peak at a measured rate of $0.01~{\rm mag~day}^{-1}$, consistent with fully-trapped $^{56}$Co decay. This is the longest rise time measured in a SLSN to date, and the first SLSN to show a rise time consistent with pair-instability supernova (PISN) models. Compared to other slowly-evolving SLSNe, it is spectroscopically similar to the prototype SN2007bi at maximum light, though lower in luminosity ($L_{\rm peak} \simeq 4.6 \times 10^{43} {\rm erg s}^{-1}$) and with a flatter peak than previous events. PS1-14bj shows a number of peculiar properties, including a near-constant color temperature for $>200$ days past peak, and strong emission lines from [O III] $\lambda$5007 and [O III] $\lambda$4363 with a velocity width of $\sim$3400 km/s, in its late-time spectra. These both suggest there is a sustained source of heating over very long timescales, and are incompatible with a simple $^{56}$Ni-powered/PISN interpretation. A modified magnetar model including emission leakage at late times can reproduce the light curve, in which case the blue continuum and [O III] features are interpreted as material heated and ionized by the inner pulsar wind nebula becoming visible at late times. Alternatively, the late-time heating could be due to interaction with a shell of H-poor circumstellar material.
• ### A DECam Search for an Optical Counterpart to the LIGO Gravitational Wave Event GW151226(1606.04538)

Sept. 8, 2016 astro-ph.CO, astro-ph.HE
We report the results of a Dark Energy Camera (DECam) optical follow-up of the gravitational wave (GW) event GW151226, discovered by the Advanced LIGO detectors. Our observations cover 28.8 deg$^2$ of the localization region in the $i$ and $z$ bands (containing 3% of the BAYESTAR localization probability), starting 10 hours after the event was announced and spanning four epochs at $2-24$ days after the GW detection. We achieve $5\sigma$ point-source limiting magnitudes of $i\approx21.7$ and $z\approx21.5$, with a scatter of $0.4$ mag, in our difference images. Given the two day delay, we search this area for a rapidly declining optical counterpart with $\gtrsim 3\sigma$ significance steady decline between the first and final observations. We recover four sources that pass our selection criteria, of which three are cataloged AGN. The fourth source is offset by $5.8$ arcsec from the center of a galaxy at a distance of 187 Mpc, exhibits a rapid decline by $0.5$ mag over $4$ days, and has a red color of $i-z\approx 0.3$ mag. These properties roughly match the expectations for a kilonova. However, this source was detected several times, starting $94$ days prior to GW151226, in the Pan-STARRS Survey for Transients (dubbed as PS15cdi) and is therefore unrelated to the GW event. Given its long-term behavior, PS15cdi is likely a Type IIP supernova that transitioned out of its plateau phase during our observations, mimicking a kilonova-like behavior. We comment on the implications of this detection for contamination in future optical follow-up observations.
• ### The Double-Peaked SN2013ge: a Type Ib/c SN with an Asymmetric Mass Ejection or an Extended Progenitor Envelope(1507.02694)

We present extensive observations of the Type Ib/c SN2013ge from -13 to +457 days, including spectra and Swift UV-optical photometry beginning 2-4 days post-explosion. This data set makes SN2013ge one of the best observed normal Type Ib/c SN at early times---when the light curve is particularly sensitive to the progenitor configuration and mixing of radioactive elements---and reveals two distinct light curve components in the UV bands. The first component rises over 4-5 days and is visible for the first week post-explosion. Spectra of the first component have blue continua and show a plethora of high velocity (~15,000 km/s) but narrow (~3500 km/s) features, indicating that the line-forming region is restricted. The explosion parameters estimated for the bulk explosion are standard for Type Ib/c SN, and there is evidence for weak He features at early times. In addition, SN2013ge exploded in a low metallicity environment and we have obtained some of the deepest radio and X-ray limits for a Type Ib/c SN to date, which constrain the progenitor mass-loss rate. We are left with two distinct progenitor scenarios for SN2013ge, depending on our interpretation of the early emission. If the first component is cooling envelope emission, then the progenitor of SN2013ge either possessed a low-mass extended envelope or ejected a portion of its envelope in the final <1 year before core-collapse. Alternatively, if the first component is due to outwardly mixed Ni-56, then our observations are consistent with the asymmetric ejection of a distinct clump of nickel-rich material at high velocities. Current models for the collision of a SN shock with a binary companion cannot reproduce both the timescale and luminosity of the early emission in SN2013ge. Finally, the spectra of the first component of SN2013ge are similar to those of the rapidly-declining SN2002bj.
• ### A Dark Energy Camera Search for an Optical Counterpart to the First Advanced LIGO Gravitational Wave Event GW150914(1602.04198)

Feb. 18, 2016 astro-ph.CO, astro-ph.IM
We report initial results of a deep search for an optical counterpart to the gravitational wave event GW150914, the first trigger from the Advanced LIGO gravitational wave detectors. We used the Dark Energy Camera (DECam) to image a 102 deg$^2$ area, corresponding to 38% of the initial trigger high-probability sky region and to 11% of the revised high-probability region. We observed in i and z bands at 4-5, 7, and 24 days after the trigger. The median $5\sigma$ point-source limiting magnitudes of our search images are i=22.5 and z=21.8 mag. We processed the images through a difference-imaging pipeline using templates from pre-existing Dark Energy Survey data and publicly available DECam data. Due to missing template observations and other losses, our effective search area subtends 40 deg$^{2}$, corresponding to 12% total probability in the initial map and 3% of the final map. In this area, we search for objects that decline significantly between days 4-5 and day 7, and are undetectable by day 24, finding none to typical magnitude limits of i= 21.5,21.1,20.1 for object colors (i-z)=1,0,-1, respectively. Our search demonstrates the feasibility of a dedicated search program with DECam and bodes well for future research in this emerging field.
• ### A Dark Energy Camera Search for Missing Supergiants in the LMC After the Advanced LIGO Gravitational Wave Event GW150914(1602.04199)

Feb. 18, 2016 astro-ph.CO, astro-ph.HE
The collapse of the core of a star is expected to produce gravitational radiation. While this process will usually produce a luminous supernova, the optical signatue could be subluminous and a direct collapse to a black hole, with the star just disappearing, is possible. The gravitational wave event GW150914 reported by the LIGO Virgo Collaboration (LVC) on 2015 September 16, was detected by a burst analysis and whose high probability spatial localization included the Large Magellanic Cloud. Shortly after the announcement of the event, we used the Dark Energy Camera to observe 102 deg$^2$ of the localization area, including a 38 deg$^2$ area centered on the LMC. Using a catalog of 152 LMC luminous red supergiants, candidates to undergo a core collapse without a visible supernova, we find that the positions of 144 of these are inside our images, and that all are detected - none have disappeared. There are other classes of candidates: we searched existing catalogs of red supergiants, yellow supergiants, Wolf-Rayet stars, and luminous blue variable stars, recovering all that were inside the imaging area. Based on our observations, we conclude that it is unlikely that GW150914 was caused by the core collapse of a supergiant in the LMC, consistent with the LIGO Collaboration analyses of the gravitational wave form as best described by a binary black hole merger. We discuss how to generalize this search for future very nearby core collapse candidates.
• ### Ejection of the massive Hydrogen-rich envelope timed with the collapse of the stripped SN2014C(1601.06806)

Jan. 25, 2016 astro-ph.HE
We present multi-wavelength observations of SN2014C during the first 500 days. These observations represent the first solid detection of a young extragalactic stripped-envelope SN out to high-energy X-rays. SN2014C was the explosion of an H-stripped progenitor star with ordinary explosion parameters. However, over the time scale of ~1yr, SN2014C experienced a complete metamorphosis and evolved from an ordinary H-poor supernova of type Ib into a strongly interacting, H-rich supernova of type IIn. Signatures of the SN shock interacting with a dense medium are observed across the spectrum. Coordinated observations with Swift, Chandra and NuSTAR have captured the evolution in detail and revealed the presence of a massive shell of ~1 Msun of hydrogen-rich material at ~6d16 cm from the explosion site. We estimate that the shell was ejected by the progenitor star in the decades to centuries before core collapse. This result poses significant challenges to current theories of massive star evolution, as it requires a physical mechanism responsible for the ejection of the deepest hydrogen layer of H-poor SN progenitors synchronized with the onset of stellar collapse. Theoretical investigations point at binary interactions and/or instabilities during the last stages of nuclear burning in massive stars as potential triggers of the time-dependent mass loss. We constrain these scenarios utilizing the sample of 183 SNe Ib/c with public radio observations. Our analysis identifies SN2014C-like signatures in ~10% of SNe with constraining radio data. This fraction is somewhat larger but reasonably consistent with the expectation from the theory of recent envelope ejection due to binary evolution IF the ejected material can survive in the close environment for 1000-10000 yrs. Alternatively, nuclear burning instabilities extending all the way to the core C-burning phase might also play a critical role.
• ### The Swift Gamma-Ray Burst Host Galaxy Legacy Survey - I. Sample Selection and Redshift Distribution(1504.02482)

Jan. 20, 2016 astro-ph.GA, astro-ph.HE
We introduce the Swift Gamma-Ray Burst Host Galaxy Legacy Survey ("SHOALS"), a multi-observatory high-redshift galaxy survey targeting the largest unbiased sample of long-duration gamma-ray burst hosts yet assembled (119 in total). We describe the motivations of the survey and the development of our selection criteria, including an assessment of the impact of various observability metrics on the success rate of afterglow-based redshift measurement. We briefly outline our host-galaxy observational program, consisting of deep Spitzer/IRAC imaging of every field supplemented by similarly-deep, multi-color optical/NIR photometry, plus spectroscopy of events without pre-existing redshifts. Our optimized selection cuts combined with host-galaxy follow-up have so far enabled redshift measurements for 110 targets (92%) and placed upper limits on all but one of the remainder. About 20% of GRBs in the sample are heavily dust-obscured, and at most 2% originate from z>5.5. Using this sample we estimate the redshift-dependent GRB rate density, showing it to peak at z~2.5 and fall by about an order of magnitude towards low (z=0) redshift, while declining more gradually towards high (z~7) redshift. This behavior is consistent with a progenitor whose formation efficiency varies modestly over cosmic history. Our survey will permit the most detailed examination to date of the connection between the GRB host population and general star-forming galaxies, directly measure evolution in the host population over cosmic time and discern its causes, and provide new constraints on the fraction of cosmic star-formation occurring in undetectable galaxies at all redshifts.