Search

We report spectroscopic and photometric follow-up of a dormant black hole (BH) candidate from Gaia DR3. The system, which we call Gaia BH2, contains a $\sim 1M_{\odot}$ red giant and a dark companion with mass $M_2 = 8.9\pm 0.3\,M_{\odot}$ that is very likely a BH. The orbital period, $P_{\rm orb} = 1277$ days, is much longer than that of any previously studied BH binary. Our radial velocity (RV) follow-up over a 7-month period spans more than 90% of the orbit's dynamic range in RV and is in excellent agreement with predictions of the Gaia solution. UV imaging and high-resolution optical spectra rule out all plausible luminous companions that could explain the orbit. The star is a bright ($G=12.3$), slightly metal-poor ($\rm [Fe/H]=-0.22$) low-luminosity giant ($T_{\rm eff}=4600\,\rm K$; $R = 7.8\,R_{\odot}$; $\log\left[g/\left({\rm cm\,s^{-2}}\right)\right] = 2.6$). The binary's orbit is moderately eccentric ($e=0.52$). The giant is strongly enhanced in $\alpha-$elements, with $\rm [\alpha/Fe] = +0.26$, but the system's Galactocentric orbit is typical of the thin disk. We obtained X-ray and radio nondetections of the source near periastron, which support BH accretion models in which the net accretion rate at the horizon is much lower than the Bondi-Hoyle-Lyttleton rate. At a distance of 1.16 kpc, Gaia BH2 is the second-nearest known BH, after Gaia BH1. Its orbit -- like that of Gaia BH1 -- seems too wide to have formed through common envelope evolution. Gaia BH1 and BH2 have orbital periods at opposite edges of the Gaia DR3 sensitivity curve, perhaps hinting at a bimodal intrinsic period distribution for wide BH binaries. Dormant BH binaries like Gaia BH1 and Gaia BH2 likely significantly outnumber their close, X-ray bright cousins, but their formation pathways remain uncertain., Comment: Accepted to MNRAS

We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly-modeled sample to-date. Using the Prospector stellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of $z=0.64^{+0.83}_{-0.32}$ ($68\%$ confidence), including 10 new or revised photometric redshifts at $z\gtrsim1$. We further find a median mass-weighted age of $t_m=0.8^{+2.71}_{-0.53}$Gyr, stellar mass of $\log(M_*/M_\odot)=9.69^{+0.75}_{-0.65}$, star formation rate of SFR=$1.44^{+9.37}_{-1.35}M_\odot$yr$^{-1}$, stellar metallicity of $\log(Z_*/Z_\odot)=-0.38^{+0.44}_{-0.42}$, and dust attenuation of $A_V=0.43^{+0.85}_{-0.36}$~mag (68\% confidence). Overall, the majority of short GRB hosts are star-forming ($\approx84\%$), with small fractions that are either transitioning ($\approx6\%$) or quiescent ($\approx10\%$); however, we observe a much larger fraction ($\approx40\%$) of quiescent and transitioning hosts at $z\lesssim0.25$, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution, with a fast-merging channel at $z>1$ and a decreased BNS formation efficiency at lower redshifts. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the BRIGHT website., 32 pages, 15 figures, 3 tables, published with ApJ

Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, $Z$, and redshift $z, \mathcal{S}(Z,z)$. Considering these uncertainties we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal{S}(Z,z)$ variations can impact the predicted intrinsic and detectable merger rates by factors $10^2$-$10^4$. We find that BHBH rates are dominantly impacted by $\mathcal{S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal{S}(Z,z)$ for BHBH, BHNS and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95% of BHBH detections to contain a BH $\gtrsim 8\,\rm{M}_{\odot}$ and have mass ratios $\lesssim 4$. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties., Comment: 25 pages of which 12 are (rainbow) figures. Main results: Fig 2, 4, 5 and 6. Code/results publicly available at https://github.com/FloorBroekgaarden/Double-Compact-Object-Mergers. Comments welcome

GRB 221009A ($z=0.151$) is one of the closest known long $\gamma$-ray bursts (GRBs). Its extreme brightness across all electromagnetic wavelengths provides an unprecedented opportunity to study a member of this still-mysterious class of transients in exquisite detail. We present multi-wavelength observations of this extraordinary event, spanning 15 orders of magnitude in photon energy from radio to $\gamma$-rays. We find that the data can be partially explained by a forward shock (FS) from a highly-collimated relativistic jet interacting with a low-density wind-like medium. Under this model, the jet's beaming-corrected kinetic energy ($E_K \sim 4\times10^{50}$ erg) is typical for the GRB population. The radio and mm data provide strong limiting constraints on the FS model, but require the presence of an additional emission component. From equipartition arguments, we find that the radio emission is likely produced by a small amount of mass ($\lesssim6\times10^{-7} M_\odot$) moving relativistically ($\Gamma\gtrsim9$) with a large kinetic energy ($\gtrsim10^{49}$ erg). However, the temporal evolution of this component does not follow prescriptions for synchrotron radiation from a single power-law distribution of electrons (e.g. in a reverse shock or two-component jet), or a thermal electron population, perhaps suggesting that one of the standard assumptions of afterglow theory is violated. GRB 221009A will likely remain detectable with radio telescopes for years to come, providing a valuable opportunity to track the full lifecycle of a powerful relativistic jet., Comment: Accepted for publication in the Astrophysical Journal Letters

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

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited., We present centimeter-band and millimeter-band afterglow observations of five long-duration γ-ray bursts (GRBs; GRB 130131A, 130420B, 130609A, 131229A, 140713A) with dust-obscured optical afterglow emission, known as “dark” GRBs. We detect the radio afterglow of two of the dark GRBs (GRB 130131A and 140713A), along with a tentative detection of a third (GRB 131229A) with the Karl G. Jansky Very Large Array (VLA). Supplemented by three additional VLA-detected dark GRBs from the literature, we present uniform modeling of their broadband afterglows. We derive high line-of-sight dust extinctions of AV,GRB ≈ 2.2– ≳ 10.6 mag. Additionally, we model the host galaxies of the six bursts in our sample, and derive host galaxy dust extinctions of AV,Host ≈ 0.3–4.7 mag. Across all tested γ-ray (fluence and duration) and afterglow properties (energy scales, geometries, and circumburst densities), we find dark GRBs to be representative of more typical unobscured long GRBs, except in fluence, for which observational biases and inconsistent classification may influence the dark GRB distribution. Additionally, we find that AV,GRB is not related to a uniform distribution of dust throughout the host, nor to the extremely local environment of the burst, indicating that a larger-scale patchy dust distribution is the cause of the high line-of-sight extinction. Since radio observations are invaluable to revealing heavily dust-obscured GRBs, we make predictions for the detection of radio emission from host star formation with the next-generation VLA. © 2022. The Author(s). Published by the American Astronomical Society., G.S. acknowledges support for this work was provided by the NSF through Student Observing Support award SOSP20B-001 from the NRAO. The Fong Group at Northwestern acknowledges support by the National Science Foundation under grant No. AST-1814782 and CAREER grant No. AST-2047919. W.F. gratefully acknowledges support by the David and Lucile Packard Foundation. T.L. acknowledges support from the Radboud Excellence Initiative. A.J.C.T. acknowledges support from the Spanish Ministry Project PID2020-118491GB-I00, Junta de Andalucía Project P20_01068 and the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). K.D.A. acknowledges support provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51403.001 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.

The Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo Interferometer Collaborations have now detected all three classes of compact binary mergers: binary black hole (BBH), binary neutron star (BNS), and neutron star-black hole (NSBH). For coalescences involving neutron stars, the simultaneous observation of gravitational and electromagnetic radiation produced by an event, has broader potential to enhance our understanding of these events, and also to probe the equation of state (EOS) of dense matter. However, electromagnetic follow-up to gravitational wave (GW) events requires rapid real-time detection and classification of GW signals, and conventional detection approaches are computationally prohibitive for the anticipated rate of detection of next-generation GW detectors. In this work, we present the first deep learning based results of classification of GW signals from NSBH mergers in \textit{real} LIGO data. We show for the first time that a deep neural network can successfully distinguish all three classes of compact binary mergers and separate them from detector noise. Specifically, we train a convolutional neural network (CNN) on $\sim 500,000$ data samples of real LIGO noise with injected BBH, BNS, and NSBH GW signals, and we show that our network has high sensitivity and accuracy. Most importantly, we successfully recover the two confirmed NSBH events to-date (GW200105 and GW200115) and the two confirmed BNS mergers to-date (GW170817 and GW190425), together with $\approx 90\%$ of all BBH candidate events from the third Gravitational Wave Transient Catalog, GWTC-3. These results are an important step towards low-latency real-time GW detection, enabling multi-messenger astronomy., 9 pages, 5 figures. Accepted for publication in Physics Letters B

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 present the discovery of the first millimeter afterglow of a short-duration $\gamma$-ray burst (SGRB) and the first confirmed afterglow of an SGRB localized by the GUANO system on Swift. Our Atacama Large Millimeter/Sub-millimeter Array (ALMA) detection of SGRB 211106A establishes an origin in a faint host galaxy detected in Hubble Space Telescope (HST) imaging at $0.7\lesssim z\lesssim1.4$. From the lack of a detectable optical afterglow, coupled with the bright millimeter counterpart, we infer a high extinction, $A_{\rm V}\gtrsim2.6$ mag along the line of sight, making this the one of the most highly dust-extincted SGRBs known to date. The millimeter-band light curve captures the passage of the synchrotron peak from the afterglow forward shock and reveals a jet break at $t_{\rm jet}=29.2^{+4.5}_{-4.0}$~days. For a presumed redshift of $z=1$, we infer an opening angle, $\theta_{\rm jet}=(15.5\pm1.4)$~degrees, and beaming-corrected kinetic energy of $\log(E_{\rm K}/{\rm erg})=51.8\pm0.3$, making this one of the widest and most energetic SGRB jets known to date. Combining all published millimeter-band upper limits in conjunction with the energetics for a large sample of SGRBs, we find that energetic outflows in high density environments are more likely to have detectable millimeter counterparts. Concerted afterglow searches with ALMA should yield detection fractions of 24-40% on timescales of $\gtrsim2$~days at rates $\approx0.8$-1.6 per year, outpacing the historical discovery rate of SGRB centimeter-band afterglows., Comment: 20 pages, 9 figures, 4 tables. Version accepted for publication in ApJ Letters

Full list of authors: Rastinejad, Jillian C.; Gompertz, Benjamin P.; Levan, Andrew J.; Fong, Wen-Fai; Nicholl, Matt; Lamb, Gavin P.; Malesani, Daniele B.; Nugent, Anya E.; Oates, Samantha R.; Tanvir, Nial R.; Postigo, Antonio de Ugarte; Kilpatrick, Charles D.; Moore, Christopher J.; Metzger, Brian D.; Ravasio, Maria Edvige; Rossi, Andrea; Schroeder, Genevieve; Jencson, Jacob; Sand, David J.; Smith, Nathan; Fernandez, Jose Feliciano Agui; Berger, Edo; Blanchard, Peter K.; Chornock, Ryan; Cobb, Bethany E.; De Pasquale, Massimiliano; Fynbo, Johan P. U.; Izzo, Luca; Kann, D. Alexander; Laskar, Tanmoy; Marini, Ester; Paterson, Kerry; Escorial, Alicia Rouco; Sears, Huei M.; Thone, Christina C., Gamma-ray bursts (GRBs) are divided into two populations1,2; long GRBs that derive from the core collapse of massive stars (for example, ref. 3) and short GRBs that form in the merger of two compact objects4,5. Although it is common to divide the two populations at a gamma-ray duration of 2 s, classification based on duration does not always map to the progenitor. Notably, GRBs with short (≲2 s) spikes of prompt gamma-ray emission followed by prolonged, spectrally softer extended emission (EE-SGRBs) have been suggested to arise from compact object mergers6,7,8. Compact object mergers are of great astrophysical importance as the only confirmed site of rapid neutron capture (r-process) nucleosynthesis, observed in the form of so-called kilonovae9,10,11,12,13,14. Here we report the discovery of a possible kilonova associated with the nearby (350 Mpc), minute-duration GRB 211211A. The kilonova implies that the progenitor is a compact object merger, suggesting that GRBs with long, complex light curves can be spawned from merger events. The kilonova of GRB 211211A has a similar luminosity, duration and colour to that which accompanied the gravitational wave (GW)-detected binary neutron star (BNS) merger GW170817 (ref. 4). Further searches for GW signals coincident with long GRBs are a promising route for future multi-messenger astronomy. © 2022, The Author(s), under exclusive licence to Springer Nature Limited., The Fong group at Northwestern acknowledges support by the National Science Foundation under grant nos. AST-1814782 and AST-1909358 and CAREER grant no. AST-2047919. W.F. gratefully acknowledges support by the David and Lucile Packard Foundation. A.J.L. and D.B.M. are supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 725246). M.N. and B.P.G. are supported by the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 948381). M.N. acknowledges a Turing Fellowship. G.P.L. is supported by the UK Science and Technology Facilities Council grant ST/S000453/1. A.R. and E.M. acknowledge support from the INAF research project ‘LBT - Supporto Arizona Italia’. J.F.A.F. acknowledges support from the Spanish Ministerio de Ciencia, Innovación y Universidades through the grant PRE2018-086507. D.A.K. and J.F.A.F. acknowledge support from Spanish National Research Project RTI2018-098104-J-I00 (GRBPhot). W. M. Keck Observatory and MMT Observatory access was supported by Northwestern University and the Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA). Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration (NASA). The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. We wish to recognize and acknowledge the very important cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. Observations reported here were obtained at the MMT Observatory, a joint facility of the University of Arizona and the Smithsonian Institution. On the basis of observations obtained at the international Gemini Observatory (programme ID GN2021B-Q-109), a programme of NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation on behalf of the Gemini Observatory partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Ministerio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovações e Comunicações (Brazil) and Korea Astronomy and Space Science Institute (Republic of Korea). Processed using the Gemini IRAF package and DRAGONS (Data Reduction for Astronomy from Gemini Observatory North and South). This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This research is based on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute, which is operated by the AURA, Inc., under NASA contract NAS 5-26555. These observations are associated with programme no. 16923. This work is partly based on observations made with the Gran Telescopio Canarias, installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, on the island of La Palma. Partly based on observations collected at the Calar Alto Astronomical Observatory, operated jointly by Instituto de Astrofísica de Andalucía (CSIC) and Junta de Andalucía. Partly based on observations made with the Nordic Optical Telescope, under programme 64-502, owned in collaboration by the University of Turku and Aarhus University, and operated jointly by Aarhus University, the University of Turku and the University of Oslo, representing Denmark, Finland and Norway, respectively, the University of Iceland and Stockholm University at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofísica de Canarias. The LBT is an international collaboration among institutions in the United States, Italy and Germany. LBT Corporation partners are: The University of Arizona on behalf of the Arizona Board of Regents; Istituto Nazionale di Astrofisica, Italy; LBT Beteiligungsgesellschaft, Germany, representing the Max Planck Society, The Leibniz Institute for Astrophysics Potsdam and Heidelberg University; The Ohio State University, representing OSU, University of Notre Dame, University of Minnesota and University of Virginia.

We present a comprehensive optical and near-infrared census of the fields of 90 short gamma-ray bursts (GRBs) discovered in 2005-2021, constituting all short GRBs for which host galaxy associations are feasible ($\approx$ 60% of the total Swift short GRB population). We contribute 245 new multi-band imaging observations across 49 distinct GRBs and 25 spectra of their host galaxies. Supplemented by literature and archival survey data, the catalog contains 335 photometric and 40 spectroscopic data sets. The photometric catalog reaches $3\sigma$ depths of $\gtrsim 24-27$ mag and $\gtrsim 23-26$ mag for the optical and near-infrared bands, respectively. We identify host galaxies for 84 bursts, in which the most robust associations make up 54% (49/90) of events, while only a small fraction, 6.7%, have inconclusive host associations. Based on new spectroscopy, we determine 17 host spectroscopic redshifts with a range of $z\approx 0.15-1.6$ and find that $\approx$ 25-44% of Swift short GRBs originate from $z>1$. We also present the galactocentric offset catalog for 83 short GRBs. Taking into account the large range of individual measurement uncertainties, we find a median of projected offset of $\approx 7.9$ kpc, for which the bursts with the most robust associations have a smaller median of $\approx 4.9$ kpc. Our catalog captures more high-redshift and low-luminosity hosts, and more highly-offset bursts than previously found, thereby diversifying the population of known short GRB hosts and properties. In terms of locations and host luminosities, the populations of short GRBs with and without detectable extended emission are statistically indistinguishable. This suggests that they arise from the same progenitors, or from multiple progenitors which form and evolve in similar environments. All of the data products are available on the BRIGHT website., Comment: 53 pages, 9 figures, 6 tables, submitted

With the upcoming Vera C.~Rubin Observatory Legacy Survey of Space and Time (LSST), it is expected that only $\sim 0.1\%$ of all transients will be classified spectroscopically. To conduct studies of rare transients, such as Type I superluminous supernovae (SLSNe), we must instead rely on photometric classification. In this vein, here we carry out a pilot study of SLSNe from the Pan-STARRS1 Medium-Deep Survey (PS1-MDS) classified photometrically with our SuperRAENN and Superphot algorithms. We first construct a sub-sample of the photometric sample using a list of simple selection metrics designed to minimize contamination and ensure sufficient data quality for modeling. We then fit the multi-band light curves with a magnetar spin-down model using the Modular Open-Source Fitter for Transients (MOSFiT). Comparing the magnetar engine and ejecta parameter distributions of the photometric sample to those of the PS1-MDS spectroscopic sample and a larger literature spectroscopic sample, we find that these samples are overall consistent, but that the photometric sample extends to slower spins and lower ejecta masses, which correspond to lower luminosity events, as expected for photometric selection. While our PS1-MDS photometric sample is still smaller than the overall SLSN spectroscopic sample, our methodology paves the way to an orders-of-magnitude increase in the SLSN sample in the LSST era through photometric selection and study., Comment: 13 pages, 6 figures, submitted to ApJ

Recent work has revealed that the light curves of hydrogen-poor (Type I) superluminous supernovae (SLSNe), thought to be powered by magnetar central engines, do not always follow the smooth decline predicted by a simple magnetar spin-down model. Here we present the first systematic study of the prevalence and properties of "bumps" in the post-peak light curves of 34 SLSNe. We find that the majority (44-76%) of events cannot be explained by a smooth magnetar model alone. We do not find any difference in supernova properties between events with and without bumps. By fitting a simple Gaussian model to the light-curve residuals, we characterize each bump with an amplitude, temperature, phase, and duration. We find that most bumps correspond with an increase in the photospheric temperature of the ejecta, although we do not see drastic changes in spectroscopic features during the bump. We also find a moderate correlation ($\rho\approx0.5$; $p\approx0.01$) between the phase of the bumps and the rise time, implying that such bumps tend to happen at a certain "evolutionary phase," $(3.7\pm1.4)t_\mathrm{rise}$. Most bumps are consistent with having diffused from a central source of variable luminosity, although sources further out in the ejecta are not excluded. With this evidence, we explore whether the cause of these bumps is intrinsic to the supernova (e.g., a variable central engine) or extrinsic (e.g., circumstellar interaction). Both cases are plausible, requiring low-level variability in the magnetar input luminosity, small decreases in the ejecta opacity, or a thin circumstellar shell or disk., Comment: updated to match accepted version

We present optical imaging and spectroscopy of SN 2018lfe, which we classify as a Type I superluminous supernova (SLSN-I) at a redshift of z = 0.3501 ± 0.0004 with a peak absolute magnitude of M r = −22.1 ± 0.1 mag, one of the brightest SLSNe discovered. SN 2018lfe was identified for follow-up using our FLEET machine-learning pipeline. Both the light curve and the spectra of SN 2018lfe are consistent with the broad population of SLSNe. We fit the light curve with a magnetar central engine model and find an ejecta mass of M ej ≈ 3.8 M ⊙, a magnetar spin period of P ≈ 2.9 ms, and a magnetic field strength of B ⊥ ≈ 2.8 × 1014 G. The magnetic field strength is near the top of the distribution for SLSNe, while the spin period and ejecta mass are near the median values of the distribution for SLSNe. From late-time imaging and spectroscopy we find that the host galaxy of SN 2018lfe has an absolute magnitude of M r = −17.85 ± 0.24, (L B = 0.029 ± 0.007L*), and an inferred metallicity of Z ≈ 0.3 Z ⊙ and star formation rate of ≈0.8 M ⊙ yr−1.

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

We present the complete set of {\it Hubble Space Telescope} imaging of the binary neutron star merger GW170817 and its optical counterpart AT 2017gfo. Including deep template imaging in F814W, F110W, F140W, and F160W at 3.4 years post-merger, we re-analyze the full light curve of AT 2017gfo across 12 bands from 5--1273 rest-frame days after merger. We obtain four new detections of the short $\gamma$-ray burst (GRB) 170817A afterglow from 109--170 rest-frame days post-merger. These detections are consistent with the previously observed $\beta=-0.6$ spectral index in the afterglow light curve with no evidence for spectral evolution. We also analyze our limits in the case of novel late-time optical and IR emission signatures, such as a kilonova afterglow or infrared dust echo, but find our limits are not constraining in these contexts. We use the new data to construct deep optical and infrared stacks, reaching limits of $M=-6.3$ to $-4.6$ mag, to analyze the local environment around AT 2017gfo and low surface brightness features in its host galaxy NGC 4993. We rule out the presence of any globular cluster at the position of AT 2017gfo to $2.3 \times 10^{4} L_{\odot}$, including those with the reddest $V-H$ colors. Finally, we analyze the substructure of NGC 4993 in deep residual imaging, and find shell features which extend up to 71.8\arcsec\ (14.2 kpc) from the center of the galaxy. We find that the shells have a cumulative stellar mass of $6.3\times10^{8} M_{\odot}$, roughly 2% the total stellar mass of NGC 4993, and mass-weighted ages of $>$3 Gyr. We conclude that it was unlikely the GW170817 progenitor system formed in the galaxy merger, which based on dynamical signatures and the stellar population in the shells mostly likely occurred 220--685 Myr ago., Comment: 21 pages, 9 figures, 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

It was proposed that a remnant stable magnetar could be formed in a binary neutron-star merger, leading to a fast X-ray transient (FXT) that can last for thousands of seconds. Recently, Xue et al. suggested that CDF-S XT2 was exactly such a kind of source. If confirmed, such emission can be used to search for electromagnetic counterparts to gravitational wave events from binary neutron-star mergers that have short gamma-ray bursts and the corresponding afterglows seen off-axis and thus too weak to be detected. Here we report the discovery of three new FXTs, XRT 170901, XRT 030511, and XRT 110919, from preliminary search over Chandra archival data. Similar to CDF-S XT2, these new FXTs had a very fast rise (less than a few ten seconds) and a plateau of X-ray flux $\sim$$1.0\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$ lasting for 1-2 ks, followed by a steep decay. Their optical/IR counterparts, if present, are very weak, arguing against a stellar flare origin for these FXTs. For XRT 170901, we identified a faint host galaxy with the source at the outskirts, very similar to CDF-S XT2. Therefore, our newly discovered FXTs are also strong candidates for magnetar-powered X-ray transients resulting from binary neutron star mergers., 13 pages, 9 figures, submitted version after addressing the referee's comments

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 provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

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.