Search

Mergers of black hole-neutron star (BHNS) binaries have now been observed by GW detectors with the recent announcement of GW200105 and GW200115. Such observations not only provide confirmation that these systems exist, but will also give unique insights into the death of massive stars, the evolution of binary systems and their possible association with gamma-ray bursts, $r$-process enrichment and kilonovae. Here we perform binary population synthesis of isolated BHNS systems in order to present their merger rate and characteristics for ground-based GW observatories. We present the results for 420 different model permutations that explore key uncertainties in our assumptions about massive binary star evolution (e.g. mass transfer, common-envelope evolution, supernovae), and the metallicity-specific star formation rate density, and characterize their relative impacts on our predictions. We find intrinsic local BHNS merger rates spanning $\mathcal{R}_{\rm{m}}^0 \approx 4$-$830\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ for our full range of assumptions. This encompasses the rate inferred from recent BHNS GW detections, and would yield detection rates of $\mathcal{R}_{\rm{det}} \approx 1$-$180\, \rm{yr}^{-1}$ for a GW network consisting of LIGO, Virgo and KAGRA at design sensitivity. We find that the binary evolution and metallicity-specific star formation rate density each impact the predicted merger rates by order $\mathcal{O}(10)$. We also present predictions for the GW detected BHNS merger properties and find that all 420 model variations predict that $\lesssim 5\%$ of the BHNS mergers have BH masses $\gtrsim 18\,M_{\odot}$, total masses $ \gtrsim 20\,M_{\odot}$, chirp masses $\gtrsim 5.5\,M_{\odot}$, mass ratios $ \gtrsim 12$ or $\lesssim 2$. Moreover, we find that massive NSs $\gtrsim 2\,M_{\odot}$ are expected to be commonly detected in BHNS mergers in almost all our model variations., 38 pages, 18 figures, accepted to MNRAS. The authors welcome suggestions and feedback. All data and code to reproduce the results in this paper are publicly available

We present X-ray and multiband optical observations of the afterglow and host galaxy of GRB 180418A, discovered by Swift/BAT and Fermi/GBM. We present a reanalysis of the GBM and BAT data deriving durations of the prompt emission of T 90 ≈ 2.56 and 1.90 s, respectively. Modeling the Fermi/GBM catalog of 1405 bursts (2008–2014) in the hardness–T 90 plane, we obtain a probability of ≈60% that GRB 180418A is a short-hard burst. From a combination of Swift/XRT and Chandra observations, the X-ray afterglow is detected to ≈38.5 days after the burst and exhibits a single power-law decline with F X ∝ t −0.98. Late-time Gemini observations reveal a faint r ≈ 25.69 mag host galaxy at an angular offset of ≈0.″16. At the likely redshift range of z ≈ 1–2.25, we find that the X-ray afterglow luminosity of GRB 180418A is intermediate between short and long gamma-ray bursts (GRBs) at all epochs during which there are contemporaneous data and that GRB 180418A lies closer to the E γ,peak–E γ,iso correlation for short GRBs. Modeling the multiwavelength afterglow with the standard synchrotron model, we derive the burst explosion properties and find a jet opening angle of θ j ≳ 9°–14°. If GRB 180418A is a short GRB that originated from a neutron star merger, it has one of the brightest and longest-lived afterglows along with an extremely faint host galaxy. If, instead, the event is a long GRB that originated from a massive star collapse, it has among the lowest-luminosity afterglows and lies in a peculiar space in terms of the hardness–T 90 and E γ,peak–E γ,iso planes.

We report on VLBI observations of the fast and blue optical transient (FBOT), AT 2018cow. At ~62 Mpc, AT 2018cow is the first relatively nearby FBOT. The nature of AT 2018cow is not clear, although various hypotheses from a tidal disruption event to different kinds of supernovae have been suggested. It had a very fast rise time (3.5 d) and an almost featureless blue spectrum although high photospheric velocities (40,000 km s$^{-1}$) were suggested early on. The X-ray luminosity was very high, ~$1.4 \times 10^{43}$ erg s$^{-1}$, larger than those of ordinary SNe, and more consistent with those of SNe associated with gamma-ray bursts. Variable hard X-ray emission hints at a long-lived "central engine." It was also fairly radio luminous, with a peak 8.4-GHz spectral luminosity of ~$4 \times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$, allowing us to make VLBI observations at ages between 22 and 287 d. We do not resolve AT 2018cow. Assuming a circularly symmetric source, our observations constrain the average apparent expansion velocity to be, Accepted for publication in MNRAS 3 figurs]es

The number of long gamma-ray bursts (GRBs) known to have occurred in the distant Universe (z greater than 5) is small (approx 15), however these events provide a powerful way of probing star formation at the onset of galaxy evolution. In this paper, we present the case for GRB100205A being a largely overlooked high-redshift event. While initially noted as a high-z candidate, this event and its host galaxy have not been explored in detail. By combining optical and near-infrared Gemini afterglow imaging (at t less than 1.3 days since burst) with deep late-time limits on host emission from the Hubble Space Telescope, we show that the most likely scenario is that GRB100205A arose in the redshift range 4-8. GRB100205A is an example of a burst whose afterglow, even at 1 hour post-burst, could only be identified by 8m class IR observations, and suggests that such observations of all optically dark bursts may be necessary to significantly enhance the number of high-redshift GRBs known., Accepted for publication in MNRAS

There is a shortage of multi-wavelength and spectroscopic followup capabilities given the number of transient and variable astrophysical events discovered through wide-field, optical surveys such as the upcoming Vera C. Rubin Observatory. From the haystack of potential science targets, astronomers must allocate scarce resources to study a selection of needles in real time. Here we present a variational recurrent autoencoder neural network to encode simulated Rubin Observatory extragalactic transient events using 1% of the PLAsTiCC dataset to train the autoencoder. Our unsupervised method uniquely works with unlabeled, real time, multivariate and aperiodic data. We rank 1,129,184 events based on an anomaly score estimated using an isolation forest. We find that our pipeline successfully ranks rarer classes of transients as more anomalous. Using simple cuts in anomaly score and uncertainty, we identify a pure (~95% pure) sample of rare transients (i.e., transients other than Type Ia, Type II and Type Ibc supernovae) including superluminous and pair-instability supernovae. Finally, our algorithm is able to identify these transients as anomalous well before peak, enabling real-time follow up studies in the era of the Rubin Observatory., 13 pages,9 figures, submitted to AAS Journals. Comments welcome!

We present detailed radio observations of the tidal disruption event (TDE) AT2019dsg, obtained with the Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA), and spanning $55-560$ days post-disruption. We find that the peak brightness of the radio emission increases until ~200 days and subsequently begins to decrease steadily. Using the standard equipartition analysis, including the effects of synchrotron cooling as determined by the joint VLA-ALMA spectral energy distributions, we find that the outflow powering the radio emission is in roughly free expansion with a velocity of $\approx 0.07c$, while its kinetic energy increases by a factor of about 5 from 55 to 200 days and plateaus at $\approx 5\times 10^{48}$ erg thereafter. The ambient density traced by the outflow declines as $\approx R^{-1.6}$ on a scale of $\approx (1-4)\times 10^{16}$ cm ($\approx 6300-25000$ $R_s$), followed by a steeper decline to $\approx 6\times 10^{16}$ cm ($\approx 37500$ $R_s$). Allowing for a collimated geometry, we find that to reach even mildly relativistic velocities ($\Gamma=2$) the outflow requires an opening angle of $\theta_j\approx 2^\circ$, which is narrow even by the standards of GRB jets; a truly relativistic outflow requires an unphysically narrow jet. The outflow velocity and kinetic energy in AT2019dsg are typical of previous non-relativistic TDEs, and comparable to those from Type Ib/c supernovae, raising doubts about the claimed association with a high-energy neutrino event., Comment: Accepted to the Astrophysical Journal

We present optical photometry and spectroscopy of SN\,2019stc (=ZTF19acbonaa), an unusual Type Ic supernova (SN Ic) at a redshift of $z=0.117$. SN\,2019stc exhibits a broad double-peaked light curve, with the first peak having an absolute magnitude of $M_r=-20.0$ mag, and the second peak, about 80 rest-frame days later, $M_r=-19.2$ mag. The total radiated energy is large, $E_{\rm rad}\approx 2.5\times 10^{50}$ erg. Despite its large luminosity, approaching those of Type I superluminous supernovae (SLSNe), SN\,2019stc exhibits a typical SN Ic spectrum, bridging the gap between SLSNe and SNe Ic. The spectra indicate the presence of Fe-peak elements, but modeling of the first light curve peak with radioactive heating alone leads to an unusually high nickel mass fraction of $f_{\rm Ni}\approx 31\%$ ($M_{\rm Ni}\approx 3.2$ M$_\odot$). Instead, if we model the first peak with a combined magnetar spin-down and radioactive heating model we find a better match with $M_{\rm ej}\approx 4$ M$_\odot$, a magnetar spin period of $P_{\rm spin}\approx 7.2$ ms and magnetic field of $B\approx 10^{14}$ G, and $f_{\rm Ni}\lesssim 0.2$ (consistent with SNe Ic). The prominent second peak cannot be naturally accommodated with radioactive heating or magnetar spin-down, but instead can be explained as circumstellar interaction with $\approx 0.7$ $M_\odot$ of hydrogen-free material located $\approx 400$ AU from the progenitor. Including the remnant mass leads to a CO core mass prior to explosion of $\approx 6.5$ M$_\odot$. The host galaxy has a metallicity of $\approx 0.26$ Z$_\odot$, low for SNe Ic but consistent with SLSNe. Overall, we find that SN\,2019stc is a transition object between normal SNe Ic and SLSNe., 14 pages, 13 figures, Accepted to ApJ

Superluminous supernovae (SLSNe) are luminous transients that can be detected to high redshifts with upcoming optical time-domain surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). An interesting open question is whether the properties of SLSNe evolve through cosmic time. To address this question, in this paper we model the multi-color light curves of all 21 Type I SLSNe from the Dark Energy Survey (DES) with a magnetar spin-down engine, implemented in the Modular Open Source Fitter for Transients (MOSFiT). With redshifts up to $z\approx 2$, this sample includes some of the highest-redshift SLSNe. We find that the DES SLSNe span a similar range of ejecta and magnetar engine parameters to previous samples of mostly lower-redshift SLSNe (spin period $P\approx 0.79-13.61$ ms, magnetic field $B\approx (0.03-7.33)\times10^{14}$ G, ejecta mass $M_{\rm ej}\approx 1.54-30.32$ M$_{\odot}$, and ejecta velocity $v_{\rm ej}\approx (0.55-1.45)\times 10^4$ km s$^{-1}$). The DES SLSN sample by itself exhibits the previously found negative correlation between $M_{\rm ej}$ and $P$, with a pronounced absence of SLSNe with low ejecta mass and rapid spin. Combining our results for the DES SLSNe with 60 previous SLSNe modeled in the same way, we find no evidence for redshift evolution in any of the key physical parameters., Comment: 12 pages, 6 figures; matches published version. Minor changes following referee comments; conclusions unchanged

The discovery of GW170817 and GRB 170817A in tandem with AT 2017gfo cemented the connection between neutron star mergers, short gamma-ray bursts (GRBs), and kilonovae. To investigate short GRB observations in the context of diverse kilonova behavior, we present a comprehensive optical and near-infrared (NIR) catalog of 85 bursts discovered over 2005-2020 on timescales of $\lesssim12$ days. The sample includes previously unpublished observations of 23 bursts, and encompasses both detections and deep upper limits. We identify 11.8% and 15.3% of short GRBs in our catalog with upper limits that probe luminosities lower than those of AT 2017gfo and a fiducial NSBH kilonovae model (for pole-on orientations), respectively. We quantify the ejecta masses allowed by the deepest limits in our catalog, constraining blue and `extremely blue' kilonova components of 14.1% of bursts to $M_{\rm ej}\lesssim0.01-0.1 M_{\odot}$. The sample of short GRBs is not particularly constraining for red kilonova components. Motivated by the large catalog as well as model predictions of diverse kilonova behavior, we investigate modified search strategies for future follow-up to short GRBs. We find that ground-based optical and NIR observations on timescales of $\gtrsim 2$ days can play a significant role in constraining more diverse outcomes. We expect future short GRB follow up efforts, such as from the {\it James Webb Space Telescope}, to expand the reach of kilonova detectability to redshifts of $z\approx 1$., Comment: 32 pages, 7 figures, 2 tables. Submitted to ApJ

In this work we study the formation of the first two black hole-neutron star (BHNS) mergers detected in gravitational waves (GW200115 and GW200105) from massive stars in wide isolated binary systems - the isolated binary evolution channel. We use 560 BHNS binary population synthesis model realizations from Broekgaarden et al. (2021a) and show that the system properties (chirp mass, component masses and mass ratios) of both GW200115 and GW200105 match predictions from the isolated binary evolution channel. We also show that most model realizations can account for the local BHNS merger rate densities inferred by LIGO-Virgo. However, to simultaneously also match the inferred local merger rate densities for BHBH and NSNS systems we find we need models with moderate kick velocities ($\sigma\lesssim 10^2\,\rm{km}\,\rm{s}^{-1}$) or high common-envelope efficiencies ($\alpha_{\rm{CE}}\gtrsim 2$) within our model explorations. We conclude that the first two observed BHNS mergers can be explained from the isolated binary evolution channel for reasonable model realizations., Comment: Accepted to ApJ Letters, all code and data is publicly available at: https://github.com/FloorBroekgaarden/NSBH_GW200105_and_GW200115

We present continued radio and X-ray observations of the previously relativistic tidal disruption event (TDE) Swift J164449.3+573451 (\sw) extending to about 9.4 years post disruption, as part of ongoing campaigns with the Jansky Very Large Array (VLA) and the \textit{Chandra} X-ray observatory. We find that the X-ray emission has faded below detectable levels, with an upper limit of $\lesssim 3.5\times 10^{-15}$ erg cm$^{-2}$ s$^{-1}$ in a 100 ks observation, while the radio emission continues to be detected and steadily fade. Both are consistent with forward shock emission from a non-relativistic outflow, although we find that the radio spectral energy distribution is better fit at these late times with an electron power law index of $p\approx 3$ (as opposed to $p\approx 2.5$ at earlier times). With the revised spectral index we find $\epsilon_B\approx 0.01$ using the radio and X-ray data, and a density of $\approx 0.04$ cm$^{3}$ at a radius of $R\approx 0.65$ pc ($R_{\rm sch}\approx 2\times 10^6$ R$_\odot$) from the black hole. The energy scale of the blastwave is $\approx 10^{52}$ erg. We also report detections of \sw\ at 3 GHz from the first two epochs of the VLA Sky Survey (VLASS), and find that $\sim 10^2$ off-axis \sw-like events to $z\sim 0.5$ may be present in the VLASS data. Finally, we find that \sw\ itself will remain detectable for decades at radio frequencies, although observations at sub-GHz frequencies will become increasingly important to characterize its dynamical evolution., Comment: Submitted to the Astrophysical Journal 10 pages, 7 figures

We thank the anonymous referee for their many suggestions that improved this paper. We thank Miguel Perez-Torres for helpful discussions. MN is supported by a Royal Astronomical Society Research Fellowship. TW is funded in part by European Research Council grant 320360 and by European Commission grant 730980. PGJ and GC acknowledge support from European Research Council Consolidator Grant 647208. GL and PC are supported by a research grant (19054) from Villum Fonden. MG is supported by the Polish NCN MAESTRO Grant 2014/14/A/ST9/00121. NI is partially supported by Polish NCN DAINA Grant 2017/27/L/ST9/03221. IA is a CIFAR Azrieli Global Scholar in the Gravity and the Extreme Universe Program and acknowledges support from that program, from the Israel Science Foundation (grant numbers 2108/18 and 2752/19), from the United States-Israel Binational Science Foundation (BSF), and from the Israeli Council forHigher Education Alon Fellowship. JB, DH, and CP were supported by NASA Grant 80NSSC18K0577. TWCacknowledges the EU Funding underMarie Sklodowska-Curie grant agreement No. 842471. LG was funded by the European Union's Horizon 2020 Framework Programme under theMarie Sklodowska-Curie grant agreement No. 839090. Thiswork has been partially supported by the Spanish grant PGC2018-095317-B-C21 within the European Funds for Regional Development (FEDER). SGGacknowledges support by FCT under Projects CRISP PTDC/FIS-AST-31546 and UIDB/00099/2020. IM is a recipient of the Australian Research Council Future Fellowship FT190100574. TMBwas funded by the CONICYT PFCHA/DOCTORADOBECAS CHILE/2017-72180113. KDA acknowledges support provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51403.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. This work is based on data collected at the European Organisation for Astronomical Research in the Southern Hemisphere, Chile, under ESO programmes 1103.D-0328 and 0104.B-0709 and as part of ePESSTO+(the advanced Public ESO Spectroscopic Survey for Transient Objects Survey), observations from the Las Cumbres Observatory network, the 6.5 meterMagellan Telescopes located at Las Campanas Observatory, Chile, the MMT Observatory, a joint facility of the University of Arizona and the Smithsonian Institution and the WilliamHerschel Telescope (programme W19B/P7). Swift data were supplied by the UK Swift Science Data Centre at the University of Leicester. The Liverpool Telescope andWilliam Herschel Telescope are operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias with financial support from the UK Science and Technology Facilities Council., At 66 Mpc, AT2019qiz is the closest optical tidal disruption event (TDE) to date, with a luminosity intermediate between the bulk of the population and the faint-and-fast event iPTF16fnl. Its proximity allowed a very early detection and triggering of multiwavelength and spectroscopic follow-up well before maximum light. The velocity dispersion of the host galaxy and fits to the TDE light curve indicate a black hole mass ≈106M , disrupting a star of ≈1M . By analysing our comprehensive UV, optical, and X-ray data, we show that the early optical emission is dominated by an outflow, with a luminosity evolution L ∝ t2, consistent with a photosphere expanding at constant velocity ( 2000 km s−1), and a line-forming region producing initially blueshifted H and He II profiles with v = 3000–10 000 km s−1. The fastest optical ejecta approach the velocity inferred from radio detections (modelled in a forthcoming companion paper from K. D. Alexander et al.), thus the same outflow may be responsible for both the fast optical rise and the radio emission – the first time this connection has been observed in a TDE. The light-curve rise begins 29 ± 2 d before maximum light, peaking when the photosphere reaches the radius where optical photons can escape. The photosphere then undergoes a sudden transition, first cooling at constant radius then contracting at constant temperature. At the same time, the blueshifts disappear from the spectrum and Bowen fluorescence lines (N III) become prominent, implying a source of far-UV photons, while the X-ray light curve peaks at ≈1041 erg s−1. Assuming that these X-rays are from prompt accretion, the size and mass of the outflow are consistent with the reprocessing layer needed to explain the large optical to X-ray ratio in this and other optical TDEs, possibly favouring accretion-powered over collision-powered outflow models., Royal Astronomical Society Research Fellowship, European Research Council (ERC) 320360 647208, European Commission European Commission Joint Research Centre 730980, Villum Fonden 19054, Polish NCN MAESTRO Grant 2014/14/A/ST9/00121, Polish NCN DAINA Grant 2017/27/L/ST9/03221, Israel Science Foundation 2108/18 2752/19, US-Israel Binational Science Foundation, Israeli Council forHigher Education Alon Fellowship, National Aeronautics & Space Administration (NASA) 80NSSC18K0577 NAS5-26555, European Union (EU) 842471 839090, European Union (EU) PGC2018-095317-B-C21, Portuguese Foundation for Science and Technology CRISP PTDC/FIS-AST-31546 UIDB/00099/2020, Australian Research Council FT190100574, CONICYT PFCHA/DOCTORADOBECAS CHILE 2017-72180113, NASA through the NASA Hubble Fellowship grant - Space Telescope Science Institute HST-HF2-51403.001-A, European Organisation for Astronomical Research in the Southern Hemisphere, Chile, under ESO programmes 1103.D-0328 0104.B-0709, WilliamHerschel Telescope W19B/P7, Science & Technology Facilities Council (STFC)

The classification of supernovae (SNe) and its impact on our understanding of the explosion physics and progenitors have traditionally been based on the presence or absence of certain spectral features. However, current and upcoming wide-field time-domain surveys have increased the transient discovery rate far beyond our capacity to obtain even a single spectrum of each new event. We must therefore rely heavily on photometric classification, connecting SN light curves back to their spectroscopically defined classes. Here we present Superphot, an open-source Python implementation of the machine-learning classification algorithm of Villar et al., and apply it to 2315 previously unclassified transients from the Pan-STARRS1 Medium Deep Survey for which we obtained spectroscopic host-galaxy redshifts. Our classifier achieves an overall accuracy of 82%, with completenesses and purities of >80% for the best classes (SNe Ia and superluminous SNe). For the worst performing SN class (SNe Ibc), the completeness and purity fall to 37% and 21%, respectively. Our classifier provides 1257 newly classified SNe Ia, 521 SNe II, 298 SNe Ibc, 181 SNe IIn, and 58 SLSNe. These are among the largest uniformly observed samples of SNe available in the literature and will enable a wide range of statistical studies of each class., Updated to match published version in ApJ. Companion paper to Villar et al. (2020). Full machine-readable tables available as ancillary files at right

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$��$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r, 51 pages, 27 figures. Accepted for publication in ApJ

Gravitational wave (GW) events detectable by LIGO and Virgo have several possible progenitors, including black hole mergers, neutron star mergers, black hole–neutron star mergers, supernovae, and cosmic string cusps. A subset of GW events is expected to produce electromagnetic (EM) emission that, once detected, will provide complementary information about their astrophysical context. To that end, the LIGO–Virgo Collaboration (LVC) sends GW candidate alerts to the astronomical community so that searches for their EM counterparts can be pursued. The DESGW group, consisting of members of the Dark Energy Survey (DES), the LVC, and other members of the astronomical community, uses the Dark Energy Camera (DECam) to perform a search and discovery program for optical signatures of LVC GW events. DESGW aims to use a sample of GW events as standard sirens for cosmology. Due to the short decay timescale of the expected EM counterparts and the need to quickly eliminate survey areas with no counterpart candidates, it is critical to complete the initial analysis of each night's images as quickly as possible. We discuss our search area determination, imaging pipeline, and candidate selection processes. We review results from the DESGW program during the first two LIGO–Virgo observing campaigns and introduce other science applications that our pipeline enables., Funding for the DES Projects has been provided by the DOE and NSF, MEC/MICINN/MINECO (Spain), STFC, HEFCE, NCSA (UIUC), KICP (U. Chicago), CCAPP (Ohio State), MIFPA (Texas A&M), CNPQ, FAPERJ, FINEP, DFG and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne Lab, UC Santa Cruz, University of Cambridge, CIEMAT-Madrid, University of Chicago, University College London, DES-Brazil Consortium, University of Edinburgh, ETH Zürich, Fermilab, University of Illinois, ICE (IEEC-CSIC), IFAE Barcelona, Lawrence Berkeley Lab, LMU München and the associated Excellence Cluster Universe, University of Michigan, NOAO, University of Nottingham, Ohio State University, University of Pennsylvania, University of Portsmouth, SLAC National Lab, Stanford University, University of Sussex, Texas A&M University, and the OzDES Membership Consortium. Based in part on observations at Cerro Tololo Inter-American Observatory at NSF’s NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The DES Data Management System is supported by the NSF under Grant Numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2015-71825, ESP2015-88861, FPA2015-68048, and Centro de Excelencia SEV-2016-0588, SEV-2016-0597 and MDM-2015-0509. Research leading to these results has received funding from the ERC under the EU’s 7th Framework Programme including grants ERC 240672, 291329 and 306478. We acknowledge support from the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO) , through project number CE110001020. This research uses services or data provided by the NOAO Science Archive. NOAO is operated by the Association of Universities for Research in Astronomy (AURA), Inc. under a cooperative agreement with the National Science Foundation. This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics. The U.S. Government retains and the publisher, by accepting the article for publication, acknowledges that the U.S. Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for U.S. Government purposes.

We present a re-analysis of the Swift Burst Alert Telescope (BAT) data around the radio detection of FRB 131104. Possible evidence of a gamma-ray counterpart was presented by DeLaunay et al. However, based on our analysis using all the available BAT data, no significant emission is found in either the temporal or the image domain. We place a 5~sigma fluence upper limit of 3.3 $\times$ 10$^{-6}$ erg cm$^{-2}$ and 2.7 $\times$ 10$^{-6}$ erg cm$^{-2}$ (15-150 keV) with an integration time of 300 s assuming a simple power-law spectrum with photon index of -1.2 and -2.0, respectively. Our result does not support the association of this FRB with a high-energy counterpart, in agreement with growing observational evidence that most FRBs are not associated to catastrophic events such as gamma-ray bursts., Comment: 12 pages, 5 figures, accepted for publication in ApJ

AT 2018hyz (=ASASSN-18zj) is a tidal disruption event (TDE) located in the nucleus of a quiescent E+A galaxy at a redshift of $z = 0.04573$, first detected by the All-Sky Automated Survey for Supernovae (ASAS-SN). We present optical+UV photometry of the transient, as well as an X-ray spectrum and radio upper limits. The bolometric light curve of AT 2018hyz is comparable to other known TDEs and declines at a rate consistent with a $t^{-5/3}$ at early times, emitting a total radiated energy of $E = 9\times10^{50}$ erg. An excess bump appears in the UV light curve about 50 days after bolometric peak, followed by a flattening beyond 250 days. The light curve shows an excess bump in the UV about 50 days after bolometric peak lasting for at least 100 days, which may be related to an outflow. We detect a constant X-ray source present for at least 86 days. The X-ray spectrum shows a total unabsorbed flux of $\sim 4\times10^{-14}$ erg cm$^{-2}$ s$^{-1}$ and is best fit by a blackbody plus power-law model with a photon index of $\Gamma = 0.8$. A thermal X-ray model is unable to account for photons $> 1$ keV, while the radio non-detection favors inverse-Compton scattering rather than a jet for the non-thermal component. We model the optical and UV light curves using the Modular Open-Source Fitter for Transients (MOSFiT) and find a best fit for a black hole of $5.2\times10^6$ M$_\odot$ partially disrupting a $0.1$ M$_\odot$ star (stripping a mass of $\sim 0.01$ M$_\odot$ for the inferred impact parameter, $\beta=0.6$). The low optical depth implied by the small debris mass may explain how we are able to see hydrogen emission with disk-like line profiles in the spectra of AT 2018hyz (see our companion paper, Short et al.~2020)., Comment: 10 pages, 10 figures, published in MNRAS

Over the past decade wide-field optical time-domain surveys have increased the discovery rate of transients to the point that $\lesssim 10\%$ are being spectroscopically classified. Despite this, these surveys have enabled the discovery of new and rare types of transients, most notably the class of hydrogen-poor superluminous supernovae (SLSN-I), with about 150 events confirmed to date. Here we present a machine-learning classification algorithm targeted at rapid identification of a pure sample of SLSN-I to enable spectroscopic and multi-wavelength follow-up. This algorithm is part of the FLEET (Finding Luminous and Exotic Extragalactic Transients) observational strategy. It utilizes both light curve and contextual information, but without the need for a redshift, to assign each newly-discovered transient a probability of being a SLSN-I. This classifier can achieve a maximum purity of about 85\% (with 20\% completeness) when observing a selection of SLSN-I candidates. Additionally, we present two alternative classifiers that use either redshifts or complete light curves and can achieve an even higher purity and completeness. At the current discovery rate, the FLEET algorithm can provide about $20$ SLSN-I candidates per year for spectroscopic follow-up with 85\% purity; with the Legacy Survey of Space and Time we anticipate this will rise to more than $\sim 10^3$ events per year., Comment: 17 pages, 12 figures, submitted to ApJ

We present optical observations of the Swift short-duration gamma-ray burst (GRB) GRB 161104A and its host galaxy at $z=0.793 \pm 0.003$. We model the multiband photometry and spectroscopy with the stellar population inference code Prospector, and explore the posterior using nested sampling. We find that the mass-weighted age $t_m = 2.12^{+0.23}_{-0.21}$~Gyr, stellar mass $\log{(M/M_\odot)} = 10.21 \pm 0.04$, metallicity $\log{(Z/Z_\odot)} = 0.08^{+0.05}_{-0.06}$, dust extinction $A_V = 0.08^{+0.08}_{-0.05}$ mag, and the star formation rate $\text{SFR} = 9.9 \times 10^{-2} M_\odot$~yr$^{-1}$. These properties, along with a prominent 4000 Angstrom break and optical absorption lines classify this host as an early-type, quiescent galaxy. Using Dark Energy Survey galaxy catalogues, we demonstrate that the host of GRB 161104A resides on the outskirts of a galaxy cluster at $z\approx 0.8$, situated $\approx 1$ Mpc from the likely brightest cluster galaxy. We also present new modeling for 20 additional short GRB hosts ($\approx33\%$ of which are early-type galaxies), finding population medians of $\log(M/M_\odot) = 9.94^{+0.88}_{-0.98}$ and $t_m = 1.07^{+1.98}_{-0.67}$~Gyr ($68\%$ confidence). We further find that the host of GRB 161104A is more distant, less massive, and younger than the four other short GRB hosts known to be associated with galaxy clusters. Cluster short GRBs have faint afterglows, in the lower $\approx 11\%$ ($\approx 30\%$) of observed X-ray (optical) luminosities. We place a lower limit on the fraction of short GRBs in galaxy clusters versus those in the field of $\approx 5-13\%$, consistent with the fraction of stellar mass $\approx 10-20\%$ in galaxy clusters at redshifts $0.1 \leq z \leq 0.8$., Comment: 20 pages, 9 figures, ApJ: Vol. 904, No. 1

We present the discovery of the radio afterglow and near-infrared (NIR) counterpart of the Swift short GRB 200522A, located at a small projected offset of $\approx 1$ kpc from the center of a young, star-forming host galaxy at $z=0.5536$. The radio and X-ray luminosities of the afterglow are consistent with those of on-axis cosmological short GRBs. The NIR counterpart, revealed by our HST observations at a rest-frame time of $\approx2.3$ days, has a luminosity of $\approx (1.3-1.7) \times 10^{42}$ erg s$^{-1}$. This is substantially lower than on-axis short GRB afterglow detections, but is a factor of $\approx 8$-$17$ more luminous than the kilonova of GW170817, and significantly more luminous than any kilonova candidate for which comparable observations exist. The combination of the counterpart's color ($i-y = -0.08\pm 0.21$; rest-frame) and luminosity cannot be explained by standard radioactive heating alone. We present two scenarios to interpret the broad-band behavior of GRB 200522A: a synchrotron forward shock with a luminous kilonova (potentially boosted by magnetar energy deposition), or forward and reverse shocks from a $\approx14^{\circ}$, relativistic ($\Gamma_0 \gtrsim 80$) jet. Models which include a combination of enhanced radioactive heating rates, low-lanthanide mass fractions, or additional sources of heating from late-time central engine activity may provide viable alternate explanations. If a stable magnetar was indeed produced in GRB 200522A, we predict that late-time radio emission will be detectable starting $\approx 0.3$-$6$ years after the burst for a deposited energy of $\approx 10^{53}$ erg. Counterparts of similar luminosity to GRB 200522A associated with gravitational wave events will be detectable with current optical searches to $\approx\!250$ Mpc., Comment: 33 pages, 13 figures, 5 tables. Submitted to AAS Journals

Automated classification of supernovae (SNe) based on optical photometric light curve information is essential in the upcoming era of wide-field time domain surveys, such as the Legacy Survey of Space and Time (LSST) conducted by the Rubin Observatory. Photometric classification can enable real-time identification of interesting events for extended multi-wavelength follow-up, as well as archival population studies. Here we present the complete sample of 5,243 "SN-like" light curves (in griz) from the Pan-STARRS1 Medium-Deep Survey (PS1-MDS). The PS1-MDS is similar to the planned LSST Wide-Fast-Deep survey in terms of cadence, filters and depth, making this a useful training set for the community. Using this dataset, we train a novel semi-supervised machine learning algorithm to photometrically classify 2,315 new SN-like light curves with host galaxy spectroscopic redshifts. Our algorithm consists of a random forest supervised classification step and a novel unsupervised step in which we introduce a recurrent autoencoder neural network (RAENN). Our final pipeline, dubbed SuperRAENN, has an accuracy of 87% across five SN classes (Type Ia, Ibc, II, IIn, SLSN-I). We find the highest accuracy rates for Type Ia SNe and SLSNe and the lowest for Type Ibc SNe. Our complete spectroscopically- and photometrically-classified samples break down into: 62.0% Type Ia (1839 objects), 19.8% Type II (553 objects), 4.8% Type IIn (136 objects), 11.7% Type Ibc (291 objects), and 1.6% Type I SLSNe (54 objects). Finally, we discuss how this algorithm can be modified for online LSST data streams., Comment: Submitted to ApJ; Companion paper to Hosseinzadeh et al.; Tables 1 and 2 available on ashleyvillar.com prior to publication

We present results from spectroscopic observations of AT 2018hyz, a transient discovered by the ASAS-SN survey at an absolute magnitude of $M_V\sim -20.2$ mag, in the nucleus of a quiescent galaxy with strong Balmer absorption lines. AT 2018hyz shows a blue spectral continuum and broad emission lines, consistent with previous TDE candidates. High cadence follow-up spectra show broad Balmer lines and He I in early spectra, with He II making an appearance after $\sim70-100$ days. The Balmer lines evolve from a smooth broad profile, through a boxy, asymmetric double-peaked phase consistent with accretion disc emission, and back to smooth at late times. The Balmer lines are unlike typical AGN in that they show a flat Balmer decrement (H$\alpha$/H$\beta\sim1.5$), suggesting the lines are collisionally excited rather than being produced via photo-ionisation. The flat Balmer decrement together with the complex profiles suggest that the emission lines originate in a disc chromosphere, analogous to those seen in cataclysmic variables. The low optical depth of material due to a possible partial disruption may be what allows us to observe these double-peaked, collisionally excited lines. The late appearance of He II may be due to an expanding photosphere or outflow, or late-time shocks in debris collisions., Comment: Published in MNRAS. Accompanied by companion paper Gomez et al. (2020)

We present and analyse a new tidal disruption event (TDE), AT2017eqx at redshift z=0.1089, discovered by Pan-STARRS and ATLAS. The position of the transient is consistent with the nucleus of its host galaxy; it peaks at a luminosity of $L \approx 10^{44}$ erg s$^{-1}$; and the spectrum shows a persistent blackbody temperature $T \gtrsim 20,000$ K with broad H I and He II emission. The lines are initially centered at zero velocity, but by 100 days the H I lines disappear while the He II develops a blueshift of $\gtrsim 5,000$ km s$^{-1}$. Both the early- and late-time morphologies have been seen in other TDEs, but the complete transition between them is unprecedented. The evolution can be explained by combining an extended atmosphere, undergoing slow contraction, with a wind in the polar direction becoming visible at late times. Our observations confirm that a lack of hydrogen a TDE spectrum does not indicate a stripped star, while the proposed model implies that much of the diversity in TDEs may be due to the observer viewing angle. Modelling the light curve suggests AT2017eqx resulted from the complete disruption of a solar-mass star by a black hole of $\sim 10^{6.3} M_\odot$. The host is another quiescent, Balmer-strong galaxy, though fainter and less centrally concentrated than most TDE hosts. Radio limits rule out a relativistic jet, while X-ray limits at 500 days are among the deepest for a TDE at this phase., Comment: Updated to match published version

We investigate the ability of current and third-generation gravitational wave (GW) detectors to determine the delay time distribution (DTD) of binary neutron stars (BNS) through a direct measurement of the BNS merger rate as a function of redshift. We assume that the DTD follows a power law distribution with a slope $\Gamma$ and a minimum merger time $t_{\rm min}$, and also allow the overall BNS formation efficiency per unit stellar mass to vary. By convolving the DTD and mass efficiency with the cosmic star formation history, and then with the GW detector capabilities, we explore two relevant regimes. First, for the current generation of GW detectors, which are only sensitive to the local universe, but can lead to precise redshift determinations via the identification of electromagnetic counterparts and host galaxies, we show that the DTD parameters are strongly degenerate with the unknown mass efficiency and therefore cannot be determined uniquely. Second, for third-generation detectors such as Einstein Telescope (ET) and Cosmic Explorer (CE), which will detect BNS mergers at cosmological distances, but with a redshift uncertainty inherent to GW-only detections ($\delta(z)/z\approx 0.1z$), we show that the DTD and mass efficiency can be well-constrained to better than 10\% with a year of observations. This long-term approach to determining the DTD through a direct mapping of the BNS merger redshift distribution will be supplemented by more near term studies of the DTD through the properties of BNS merger host galaxies at $z\approx 0$ (Safarzadeh & Berger 2019)., Comment: 10 pages, Accepted to ApJ Letters

We present a nuclear transient event, PS1-13cbe, that was first discovered in the Pan-STARRS1 survey in 2013. The outburst occurred in the nucleus of the galaxy SDSS J222153.87+003054.2 at z = 0.123 55, which was classified as a Seyfert 2 in a pre-outburst archival Sloan Digital Sky Survey (SDSS) spectrum. PS1-13cbe showed the appearance of strong broad H α and H β emission lines and a non-stellar continuum in a Magellan spectrum taken 57 d after the peak of the outburst that resembled the characteristics of a Seyfert 1. These broad lines were not present in the SDSS spectrum taken a decade earlier and faded away within 2 yr, as observed in several late-time MDM spectra. We argue that the dramatic appearance and disappearance of the broad lines and a factor of ∼8 increase in the optical continuum are most likely caused by variability in the pre-existing accretion disc than a tidal disruption event, supernova, or variable obscuration. The time-scale for the turn-on of the optical emission of ∼70 d observed in this transient is among the shortest observed in a ‘changing-look’ active galactic nucleus.

Because core-collapse supernovae are the explosions of massive stars, which have relatively short lifetimes, they occur almost exclusively in galaxies with active star formation. On the other hand, the Type Ibn supernova PS1-12sk exploded in an environment much more typical of thermonuclear (Type Ia) supernovae: on the outskirts of the brightest elliptical galaxy in a galaxy cluster. The lack of any obvious star formation at that location presented a challenge to models of Type Ibn supernovae as the explosions of very massive Wolf-Rayet stars. Here we present a supplementary search for star formation at the site of PS1-12sk, now that the supernova has faded, via deep ultraviolet imaging of the host cluster with the Hubble Space Telescope. We do not detect any ultraviolet emission within 1 kpc of the supernova location, which allows us deepen the limit on star formation rate by an order of magnitude compared to the original study on this event. In light of this new limit, we discuss whether the progenitors of Type Ibn supernovae can be massive stars and what reasonable alternatives have been proposed.

Recently born magnetars are promising candidates for the engines powering fast radio bursts (FRBs). The focus thus far has been placed on millisecond magnetars born in rare core-collapse explosions, motivated by the star forming dwarf host galaxy of the repeating FRB 121102, which is remarkably similar to the hosts of superluminous supernovae (SLSNe) and long gamma-ray bursts (LGRB). However, long-lived magnetars may also be created in binary neutron star (BNS) mergers, in the small subset of cases with a sufficiently low total mass for the remnant to avoid collapse to a black hole, or in the accretion-induced collapse (AIC) of a white dwarf. A BNS FRB channel will be characterized by distinct host galaxy and spatial offset distributions than the SLSNe/LGRB channel; we anticipate a similar host population, although possibly different offset distribution for AIC events. We show that both the BNS and AIC channels are consistent with the recently reported FRB 180924, localized by ASKAP to a massive quiescent host galaxy with an offset of about 1.4 effective radii. Using models calibrated to FRB 121102, we make predictions for the dispersion measure, rotation measure, and persistent radio emission from magnetar FRB sources born in BNS mergers or AIC, and show these are consistent with upper limits from FRB 180924 for reasonable parameters. Depending on the rate of AIC, and the fraction of BNS mergers leaving long-lived stable magnetars, the birth rate of repeating FRB sources associated with older stellar populations could be comparable to that of the core-collapse channel. We also discuss potential differences in the repetition properties of these channels, as a result of differences in the characteristic masses and magnetic fields of the magnetars., Comment: submitted to ApJ; comments welcome!

On 2019 April 25.346 and 26.640 UT the LIGO and Virgo gravitational wave (GW) observatories announced the detection of the first candidate events in Observing Run 3 that contain at least one neutron star. S190425z is a likely binary neutron star (BNS) merger at $d_L = 156 \pm 41$ Mpc, while S190426c is possibly the first NS-BH merger ever detected, at $d_L = 377 \pm 100$ Mpc, although with marginal statistical significance. Here we report our optical follow-up observations for both events using the MMT 6.5-m telescope, as well as our spectroscopic follow-up of candidate counterparts (which turned out to be unrelated) with the 4.1-m SOAR telescope. We compare to publicly reported searches, explore the overall areal coverage and depth, and evaluate those in relation to the optical/NIR kilonova emission from the BNS merger GW170817, to theoretical kilonova models, and to short GRB afterglows. We find that for a GW170817-like kilonova, the partial volume covered spans up to about 40% for S190425z and 60% for S190426c. For an on-axis jet typical of short GRBs, the search effective volume is larger, but such a configuration is expected in at most a few percent of mergers. We further find that wide-field $\gamma$-ray and X-ray limits rule out luminous on-axis SGRBs, for a large fraction of the localization regions, although these searches are not sufficiently deep in the context of the $\gamma$-ray emission from GW170817 or off-axis SGRB afterglows. The results indicate that some optical follow-up searches are sufficiently deep for counterpart identification to about 300 Mpc, but that localizations better than 1000 deg$^2$ are likely essential., Comment: The first six authors contributed equally to this work

On 2019 August 14 the Laser Interferometer Gravitational Wave Observatory (LIGO) and the Virgo gravitational wave interferometer announced the detection of a binary merger, S190814bv, with a low false alarm rate (FAR) of about 1 in $1.6\times 10^{25}$ years, a distance of $267\pm 52$ Mpc, a 90\% (50\%) localization region of about 23 (5) deg$^2$, and a probability of being a neutron star--black hole (NS-BH) merger of $>99\%$. The LIGO/Virgo Collaboration (LVC) defines NS-BH such that the lighter binary member has a mass of $5$ M$_\odot$, and this classification is in principle consistent with a BH-BH merger depending on the actual upper mass cut-off for neutron stars. Additionally, the LVC designated a probability that the merger led to matter outside the final BH remnant of $, 10 pages, 1 table, 4 figures, accepted to ApJL

Whether stars could have driven the reionization of the intergalactic medium depends critically on the proportion of ionizing radiation that escapes the galaxies in which it is produced. Spectroscopy of gamma-ray burst (GRB) afterglows can be used to estimate the opacity to extreme ultraviolet (EUV) radiation along the lines-of-sight to the bursts. Assuming that long-duration GRBs trace the locations of the massive stars dominating EUV production, the average escape fraction of ionizing radiation can be calculated independently of galaxy size or luminosity. Here we present a compilation of H i column density (N HI ) measures for 140 GRBs in the range 1.6 < z < 6.7. Although the sample is heterogeneous, in terms of spectral resolution and signal-to-noise ratio, fits to the Ly α absorption line provide robust constraints on N HI , even for spectra of insufficient quality for other purposes. Thus we establish an escape fraction at the Lyman limit of (f esc ) ≈ 0.005, with a 98 per cent confidence upper limit of (f esc ) ≈ 0.015. This analysis suggests that stars provide a small contribution to the ionizing radiation budget at z < 5. At higher redshifts firm conclusions are limited by the small size of the GRB sample (7/140), but any decline in average H i column density seems to be modest. We also find no significant correlation of N HI with galaxy UV luminosity or host stellar mass. We discuss in some detail potential biases and argue that, while not negligible, systematic errors in f esc are unlikely to be more than a factor ~2 in either direction, and so would not affect the primary conclusions. Given that many GRB hosts are low-metallicity dwarf galaxies with high specific star-formation rates, these results present a particular problem for the hypothesis that such galaxies dominated the reionization of the Universe. © 2018 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society., JJ acknowledges support from Nederlandse Onderzoekschool Voor Astronomie(NOVA) and The Netherlands Organisation for Scientific Research (NWO) - The Sao Paulo Research Foundation (FAPESP) grant for advanced instrumentation in astronomy. KEH acknowledges support by a Project Grant (162948-051) from The Icelandic Research Fund. DAK acknowledges support from the Spanish research project AYA2014-58381-P and Juan de la Cierva Incorporacion IJCI-2015-26153. AJL and ERS acknowledge Science and Technology Facilities Council (STFC) consolidated grant ST/L000733/1. NRT and KW acknowledge Science and Technology Facilities Council (STFC) consolidated grant ST/N000757/1. AC acknowledges National Aeronautics and Space Administration (NASA) grant NNX15AP95A. AdUP acknowledges support from a Ramon y Cajal fellowship (RyC-2012-09975), a 2016 Banco Bilbao Vizcaya Argentaria (BBVA) Foundation Grant for Researchers and Cultural Creators, and from the Spanish research project AYA 2014-58381-P. RC acknowledges partial support from National Aeronautics and Space Administration (NASA) Swift grant NNX16AB04G.This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 725246).

We present comprehensive multiwavelength radio to X-ray observations of GRB 181201A spanning from $\approx150$ s to $\approx163$ days after the burst, comprising the first joint ALMA-VLA-GMRT observations of a gamma-ray burst (GRB) afterglow. The radio and mm-band data reveal a distinct signature at $\approx3.9$ days, which we interpret as reverse shock (RS) emission. Our observations present the first time that a single radio-frequency spectral energy distribution can be decomposed directly into RS and forward shock (FS) components. We perform detailed modeling of the full multiwavelength data set, using Markov Chain Monte Carlo sampling to construct the joint posterior density function of the underlying physical parameters describing the RS and FS synchrotron emission. We uncover and account for all degeneracies in the model parameters. The joint RS-FS modeling reveals a weakly magnetized ($\sigma\approx3\times10^{-3}$), mildly relativistic RS, from which we derive an initial bulk Lorentz factor of $\Gamma_0\approx103$ for the GRB jet. Our results support the hypothesis that low-density environments are conducive to the observability of RS emission. We compare our observations to other events with strong RS detections, and find a likely observational bias selecting for longer lasting, non-relativistic reverse shocks. We present and begin to address new challenges in modeling posed by the present generation of comprehensive, multi-frequency data sets., Comment: Accepted version

We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_{\nu}\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_{\nu}\propto \nu^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $\theta_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $\Gamma_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $\theta_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $\epsilon_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $\epsilon_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5\sigma$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (\Gamma\beta)^{-\alpha}$ into the environment, finding that shallow stratification indexes $\alpha\le6$ are disfavored., Comment: version accepted for publication in ApJL, 13 pages, 6 figures

We investigate secular evolution in hierarchical quadruple systems as a formation channel of mass-gap black holes (with masses of about $3-5$ $M_{\odot}$) in systems that will eventually lead to binary black hole mergers detectable by ground-based gravitational wave detectors (LIGO/Virgo). We show that in a 3+1 hierarchical system, two episodes of induced mergers would first cause two neutron stars to merge and form a mass-gap black hole, which will subsequently merge with another (more massive) black hole through a second induced merger. We demonstrate that such systems are stable to flybys, and their formation would predict a high mass ratio and eccentric merger of a mass-gap black hole with a more massive black hole companion. Such a formation channel may explain observed gravitational wave events such as the recently-discovered LIGO/Virgo events S190814bv and S190924h., Comment: Accepted for publication at ApJ Letters

We present optical photometry and spectroscopy of SN 2016iet, an unprecedented Type I supernova (SN) at $z=0.0676$ with no obvious analog in the existing literature. The peculiar light curve has two roughly equal brightness peaks ($\approx -19$ mag) separated by 100 days, and a subsequent slow decline by 5 mag in 650 rest-frame days. The spectra are dominated by emission lines of calcium and oxygen, with a width of only $3400$ km s$^{-1}$, superposed on a strong blue continuum in the first year, and with a large ratio of $L_{\rm [Ca\,II]}/L_{\rm [O\,I]}\approx 4$ at late times. There is no clear evidence for hydrogen or helium associated with the SN at any phase. We model the light curves with several potential energy sources: radioactive decay, central engine, and circumstellar medium (CSM) interaction. Regardless of the model, the inferred progenitor mass near the end of its life (i.e., CO core mass) is $\gtrsim 55$ M$_\odot$ and up to $120$ M$_\odot$, placing the event in the regime of pulsational pair instability supernovae (PPISNe) or pair instability supernovae (PISNe). The models of CSM interaction provide the most consistent explanation for the light curves and spectra, and require a CSM mass of $\approx 35$ M$_\odot$ ejected in the final decade before explosion. We further find that SN 2016iet is located at an unusually large offset ($16.5$ kpc) from its low metallicity dwarf host galaxy ($Z\approx 0.1$ Z$_\odot$, $M\approx 10^{8.5}$ M$_\odot$), supporting the PPISN/PISN interpretation. In the final spectrum, we detect narrow H$\alpha$ emission at the SN location, likely due to a dim underlying galaxy host or an H II region. Despite the overall consistency of the SN and its unusual environment with PPISNe and PISNe, we find that the inferred properties of SN\,2016iet challenge existing models of such events., Comment: 26 Pages, 17 Figures, Submitted to ApJ

This corrects the article DOI: 10.1103/PhysRevLett.121.091102.

We present the largest and deepest late-time radio and millimeter survey to date of superluminous supernovae (SLSNe) and long-duration gamma-ray bursts (LGRBs) to search for associated nonthermal synchrotron emission. Using the Karl G. Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA), we observed 43 sources at 6 and 100 GHz on a timescale of ∼ 1–19 yr post-explosion. We do not detect radio/millimeter emission from any of the sources, with the exception of a 6 GHz detection of PTF10hgi, as well as the detection of 6 GHz emission near the location of the SLSN PTF12dam, which we associate with its host galaxy. We use our data to place constraints on central engine emission due to magnetar wind nebulae and off-axis relativistic jets. We also explore nonrelativistic emission from the SN ejecta, and place constraints on obscured star formation in the host galaxies. In addition, we conduct a search for fast radio bursts (FRBs) from some of the sources using VLA phased-array observations; no FRBs are detected to a limit of 16 mJy (7σ; 10 ms duration) in about 40 minutes on source per event. A comparison to theoretical models suggests that continued radio monitoring may lead to detections of persistent radio emission on timescales of ≳ a decade.

One of the key challenges of real-time detection and parameter estimation of gravitational waves from compact binary mergers is the computational cost of conventional matched-filtering and Bayesian inference approaches. In particular, the application of these methods to the full signal parameter space available to the gravitational-wave detectors, and/or real-time parameter estimation is computationally prohibitive. On the other hand, rapid detection and inference are critical for prompt follow-up of the electromagnetic and astro-particle counterparts accompanying important transients, such as binary neutron-star and black-hole neutron-star mergers. Training deep neural networks to identify specific signals and learn a computationally efficient representation of the mapping between gravitational-wave signals and their parameters allows both detection and inference to be done quickly and reliably, with high sensitivity and accuracy. In this work we apply a deep-learning approach to rapidly identify and characterize transient gravitational-wave signals from binary neutron-star mergers in real LIGO data. We show for the first time that artificial neural networks can promptly detect and characterize binary neutron star gravitational-wave signals in real LIGO data, and distinguish them from noise and signals from coalescing black-hole binaries. We illustrate this key result by demonstrating that our deep-learning framework classifies correctly all gravitational-wave events from the Gravitational-Wave Transient Catalog, GWTC-1 [Phys. Rev. X 9 (2019), 031040]. These results emphasize the importance of using realistic gravitational-wave detector data in machine learning approaches, and represent a step towards achieving real-time detection and inference of gravitational waves., 10 pages, 8 figures, 1 table. Updated grant information. Published in Physics Letters B

The most promising variation of the standard siren technique combines gravitational-wave (GW) data for binary neutron star (BNS) mergers with redshift measurements enabled by their electromagnetic (EM) counterparts, to constrain cosmological parameters such as $H_0$, $\Omega_m$, and $w_0$. Here we evaluate the near- and long-term prospects of multi-messenger cosmology in the era of future GW observatories: Advanced LIGO Plus (A+, 2025), Voyager-like detectors (2030s), and Cosmic Explorer-like detectors (CE, 2035 and beyond). We show that the BNS horizon distance of $\approx 700$ Mpc for A+ is well-matched to the sensitivity of the Vera C. Rubin Observatory (VRO) for kilonova detections. We find that one year of joint A+ and VRO observations will constrain the value of $H_0$ to percent-level precision, given a small investment of VRO time dedicated to target-of-opportunity GW follow-up. In the Voyager era, the BNS-kilonova observations begin to constrain $\Omega_m$ with an investment of a few percent of VRO time. With the larger BNS horizon distance in the Cosmic Explorer era, on-axis short gamma-ray bursts (SGRBs) and their afterglows (though accompanying only some of the GW-detected mergers) supplant kilonovae as the most promising counterparts for redshift identification. We show that five years of joint observations with Cosmic Explorer-like facilities and a next-generation gamma-ray satellite with localization capabilities similar to that presently possible with Swift could constrain both $\Omega_m$ and $w_0$ to $15-20\%$. We therefore advocate for a robust target-of-opportunity (ToO) program with VRO, and a wide-field gamma-ray satellite with improved sensitivity in the 2030s, to enable standard siren cosmology with next-generation gravitational wave facilities., Comment: The Astrophysical Journal Letters accepted version

Massive, rapidly spinning magnetar remnants produced as a result of binary neutron-star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a synchrotron radio signal from the interaction of the ejecta with the circumburst medium. Here, we present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (GRBs; z σ limits on radio continuum emission of F ν ≲ 6–20 μJy at the burst positions, or L ν ≲ (0.6–8.3) × 1028 erg s−1 Hz−1. Comparing these limits with new light-curve modeling that properly incorporates relativistic effects, we obtain limits on the energy deposited into the ejecta of E ej ≲ (0.6–6.7) × 1052 erg ( erg) for an ejecta mass of 0.03 M ⊙ (0.1 M ⊙). We present a uniform reanalysis of 27 short GRBs with 5.5–6.0 GHz observations, and find that ≳50% of short GRBs did not form stable magnetar remnants in their mergers. Assuming short GRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction f GW190425 of cases, we place constraints on the maximum mass of a nonrotating neutron star (NS; Tolman–Oppenheimer–Volkoff mass; M TOV), finding for f GW190425 = 0.4; this limit increases for larger values of f GW190425. The detection (or lack thereof) of radio remnants in untargeted surveys such as the VLA Sky Survey (VLASS) could provide more stringent constraints on the fraction of mergers that produce stable remnants. If radio remnants are discovered in VLASS, this suggests that short GRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce.

The LIGO/Virgo Scientific Collaboration (LVC) recently announced the detection of a compact object binary merger, GW190425, with a total mass of $3.4^{+0.3}_{-0.1}$ M$_{\odot}$, and individual component masses in the range of about 1.1 to 2.5 $M_{\odot}$. If the constituent compact objects are neutron stars, then the total mass is five standard deviations higher than the mean of $2.66\pm 0.12$ M$_{\odot}$ for Galactic binary neutron stars (BNSs). LVC suggests such massive BNS systems are born from a fast-merging channel, and therefore, their non-detection in the Galaxy to be due to a selection effect. However, we are unable to reconcile the inferred formation efficiency from the reported merger rate, $\mathcal{R}_{\rm GW190425}=460^{+1050}_{-390}$ yr$^{-1}$ Gpc$^{-3}$, with predictions from our own study for fast-merging BNS systems. Moreover, the comparable merger rates of GW190425, and GW170817 are possibly in tension with our results for two reasons: (i) more massive systems are expected to have a lower formation rate, and (ii) fast merging channels should constitute $\lesssim 10\%$ of the total BNS systems if case BB unstable mass transfer is permitted to take place as a formation pathway. We argue that to account for the high merger rate of GW190425 as a BNS system requires: (i) revisiting our understanding of NS formation in supernova explosions, or (ii): that more massive NSs need to be preferentially born with either very weak or very high magnetic fields so that they would be undetectable in the radio surveys. Perhaps massive NSs detected in NS-white dwarf binaries are our clues to the formation path of GW190425 systems., Comment: Accepted for publication in ApJ

Despite indications that superluminous supernovae (SLSNe) originate from massive progenitors, the lack of a uniformly analyzed statistical sample has so far prevented a detailed view of the progenitor mass distribution. Here we present and analyze the pre-explosion mass distribution of hydrogen-poor SLSN progenitors as determined from uniformly modelled light curves of 62 events. We construct the distribution by summing the ejecta mass posteriors of each event, using magnetar light curve models presented in our previous works (and using a nominal neutron star remnant mass). The resulting distribution spans $3.6-40$ M$_{\odot}$, with a sharp decline at lower masses, and is best fit by a broken power law described by ${\rm d}N/{\rm dlog}M \propto M^{-0.41 \pm 0.06}$ at $3.6-8.6$ M$_{\odot}$ and $\propto M^{-1.26 \pm 0.06}$ at $8.6-40$ M$_{\odot}$. We find that observational selection effects cannot account for the shape of the distribution. Relative to Type Ib/c SNe, the SLSN mass distribution extends to much larger masses and has a different power-law shape, likely indicating that the formation of a magnetar allows more massive stars to explode as some of the rotational energy accelerates the ejecta. Comparing the SLSN distribution with predictions from single and binary star evolution models, we find that binary models for a metallicity of $Z\lesssim 1/3$ Z$_{\odot}$ are best able to reproduce its broad shape, in agreement with the preference of SLSNe for low metallicity environments. Finally, we uncover a correlation between the pre-explosion mass and the magnetar initial spin period, where SLSNe with low masses have slower spins, a trend broadly consistent with the effects of angular momentum transport evident in models of rapidly-rotating carbon-oxygen stars., 18 pages, 11 figures, Submitted to ApJ

The discovery of a persistent radio source coincident with the first repeating fast radio burst, FRB 121102, and offset from the center of its dwarf host galaxy has been used as evidence for a link with young millisecond magnetars born in superluminous supernovae (SLSNe) or long-duration gamma-ray bursts (LGRBs). A prediction of this scenario is that compact radio sources offset from the centers of dwarf galaxies may serve as signposts for at least some FRBs. Recently, Reines et al. 2019 presented the discovery of 20 such radio sources in nearby ($z\lesssim 0.055$) dwarf galaxies, and argued that these cannot be explained by emission from HII regions, normal supernova remnants, or normal radio supernovae. Instead, they attribute the emission to accreting wandering massive black holes. Here, we explore the alternative possibility that these sources are analogs of FRB 121102. We compare their properties -- radio luminosities, spectral energy distributions, light curves, ratios of radio-to-optical flux, and spatial offsets -- to FRB 121102, a few other well-localized FRBs, and potentially related systems, and find that these are all consistent as arising from the same population. We further compare their properties to the magnetar nebula model used to explain FRB 121102, as well as to theoretical off-axis LGRB light curves, and find overall consistency. Finally, we find a consistent occurrence rate relative to repeating FRBs and LGRBs. We outline key follow-up observations to further test these possible connections., Comment: 11 pages, 6 figures; replaced with referee revisions including minor changes to figures and text

et al., arXiv:1907.01981v3, The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution (≈2farcs5), sensitivity (a 1σ goal of 70 μJy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2–4 GHz). The first observations began in 2017 September, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hr of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (decl. > −40°), a total of 33 885 deg${}^{2}$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations., The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. CIRADA is funded by a grant from the Canada Foundation for Innovation 2017 Innovation Fund (Project 35999), as well as by the Provinces of Ontario, British Columbia, Alberta, Manitoba and Quebec. This research made use of NASA’s Astrophysics Data System (ADS) Abstract Service. Part of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. S.C. acknowledges support from the NSF (AAG 1815242). S.C. and T.J.W.L. are members of the NANOGrav Physics Frontier Center, which is supported by the NSF (award number 1430284. C.J.L. is supported by NSF award 1611606. Partial support for the work of L. Rudnick comes from NSF grants AST-1211595 and 1714205 to the University of Minnesota. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant No. ACI1548562. H.A. benefited from grant CIIC 218/2019 of the University of Guanajuato, Mexico. B.R.K. acknowledges support from the NRAO NINE program and Office of Diversity and Inclusion. G.R.S. acknowledges support from NSERC Discovery Grant RGPIN-06569-2016. S.B.S. is supported by NSF award #1458952. S.P.O. acknowledges financial support from the Deutsche Forschungsgemeinschaft (DFG) under grant BR2026/23. L.L. acknowledges the financial support of DGAPA, UNAM (IN112417), and CONACyT, Mexico.

We present simulated observations to assess the ability of LSST and the WFD survey to detect and characterize kilonovae - the optical emission associated with binary neutron star (and possibly black hole - neutron star) mergers. We expand on previous studies in several critical ways by exploring a range of kilonova models and several choices of cadence, as well as by evaluating the information content of the resulting light curves. We find that, depending on the precise choice of cadence, the WFD survey will achieve an average kilonova detection efficiency of $\approx 1.6-2.5\%$ and detect only $\approx 3-6$ kilonovae per year. The detected kilonovae will be within the detection volume of Advanced LIGO/Virgo (ALV). By refitting the best resulting LSST light curves with the same model used to generate them we find the model parameters are generally weakly constrained, and are accurate to at best a factor of $2-3$. Motivated by the finding that the WFD will yield a small number of kilonova detections, with poor light curves and marginal information content, and that the detections are in any case inside the ALV volume, we argue that target-of-opportunity follow-up of gravitational wave triggers is a much more effective approach for kilonova studies. We outline the qualitative foundation for such a program with the goal of minimizing the impact on LSST operations. We argue that observations in the $gz$-bands with a total time investment per event of $\approx 1.5$ hour per 10 deg$^2$ of search area is sufficient to rapidly detect and identify kilonovae with $\gtrsim 90\%$ efficiency. For an estimated event rate of $\sim20$ per year visible to LSST, this accounts for $\sim1.5\%$ of the total survey time. In this regime, LSST has the potential to be a powerful tool for kilonovae discovery, with detected events handed off to other narrow-field facilities for further monitoring., 17 Pages, 11 Figures, 1 Table; Submitted to ApJ

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, 12 pages of main text

We study the spectroscopic evolution of superluminous supernovae (SLSNe) later than 100 days after maximum light. We present new data for Gaia16apd and SN2017egm, and analyse these with a larger sample comprising 41 spectra of 12 events. The spectra become nebular within 2-4 $e$-folding times after light curve peak, with the rate of spectroscopic evolution correlated to the light curve timescale. Emission lines are identified with well-known transitions of oxygen, calcium, magnesium, sodium and iron. SLSNe are differentiated from other Type Ic SNe by a prominent O I $\lambda$7774 line and higher-ionisation states of oxygen. The iron-dominated region around 5000 \AA\ is more similar to broad-lined SNe Ic than to normal SNe Ic. Principal Component Analysis shows that 5 `eigenspectra' capture 75% of the variance, while a clustering analysis shows no clear evidence for multiple SLSN sub-classes. Line velocities are 5000--8000 km/s, and show stratification of the ejecta. O I $\lambda$7774 likely arises in a dense inner region that also produces calcium emission, while [O I] $\lambda$6300 comes from further out until 300--400 days. The luminosities of O I $\lambda$7774 and Ca II suggest significant clumping, in agreement with previous studies. Ratios of [Ca II]$\lambda$7300/[O I]$\lambda$6300 favour progenitors with relatively massive helium cores, likely $\gtrsim 6$ M$_\odot$, though more modelling is required here. SLSNe with broad light curves show the strongest [O I] $\lambda$6300, suggesting larger ejecta masses. We show how the inferred velocity, density and ionisation structure point to a central power source., Comment: Accepted in ApJ, updated to match accepted version

We present detailed multi-wavelength observations of GRB 161219B at $z=0.1475$, spanning the radio to X-ray regimes, and the first ALMA light curve of a GRB afterglow. The cm- and mm-band observations before $8.5$ d require emission in excess of that produced by the afterglow forward shock (FS). These data are consistent with radiation from a refreshed reverse shock (RS) produced by the injection of energy into the FS, signatures of which are also present in the X-ray and optical light curves. We infer a constant-density circumburst environment with an extremely low density, $n_0\approx 3\times10^{-4}$ cm$^{-3}$ and show that this is a characteristic of all strong RS detections to date. The VLA observations exhibit unexpected rapid variability on $\sim$ minute timescales, indicative of strong interstellar scintillation. The X-ray, ALMA, and VLA observations together constrain the jet break time, $t_{\rm jet}\approx32$ day, yielding a wide jet opening angle of $\theta_{\rm jet}\approx13^{\circ}$, implying beaming corrected $\gamma$-ray and kinetic energies of $E_{\gamma}\approx4.9\times10^{48}$ erg and $E_{\rm K}\approx1.3\times10^{50}$ erg, respectively. Comparing the RS and FS emission, we show that the ejecta are only weakly magnetized, with relative magnetization, $R_{\rm B}\approx1$, compared to the FS. These direct, multi-frequency measurements of a refreshed RS spanning the optical to radio bands highlight the impact of radio and millimeter data in probing the production and nature of GRB jets., Comment: 23 pages, 17 figures, 9 tables. Published in ApJ