Search

We present X-ray to radio frequency observations of the bright long gamma-ray burst GRB 210702A. Our ALMA 97.5 GHz observations show a significant rebrightening by a factor of ~2 beginning at 8.2 days post-burst and rising to peak brightness at 18.1 days before declining again. This is the first such rebrightening seen in a millimeter afterglow light curve. A standard forward shock model in a stellar wind circumburst medium can explain most of our X-ray, optical and millimeter observations prior to the rebrightening, but significantly over-predicts the self-absorbed radio emission, and cannot explain the millimeter rebrightening. We investigate possible explanations for the millimeter rebrightening and find that energy injection or a reverse shock from a late-time shell collision are plausible causes. Similar to other bursts, our radio data may require alternative scenarios such as a thermal electron population or a structured jet to explain the data. Our observations demonstrate that millimeter light curves can exhibit some of the rich features more commonly seen in optical and X-ray afterglow light curves, motivating further millimeter wavelength studies of GRB afterglows., Comment: 22 pages, 10 figures, accepted for publication in ApJ

We present optical/UV observations and the spectroscopic classification of the transient AT2023vto as a tidal disruption event (TDE) at z = 0.4846. The spectrum is dominated by a broad He II $\lambda$4686 emission line, with a width of ~ $3.76 \times 10^4$ km/s and a blueshift of ~ $1.05 \times 10^4$ km/s, classifying it as a member of the TDE-He class. The light curve exhibits a long rise and decline timescale, with a large peak absolute magnitude of M$_g$ ~ -23.6, making it the most luminous of the classical optical TDEs (H, H+He, He) discovered to date by about 2 mag (and ~ 4 mag compared to the mean of the population). The light curve exhibits a persistent blue color of g - r ~ -0.4 mag throughout its evolution, similar to other TDEs, but distinct from supernovae. We identify the host galaxy of AT2023vto in archival Pan-STARRS images and find that the transient is located at the galaxy center, and that its inferred central black hole mass is ~ $10^7~M_{\odot}$. Modeling the light curves of AT2023vto, we find that it resulted from the disruption of a ~ 9 $M_{\odot}$ star by a ~$10^7~M_{\odot}$ supermassive black hole. The star mass is about 5 times larger than the highest star masses previously inferred in TDEs, and the black hole mass is at the high end of the distribution. AT2023vto is comparable in luminosity and timescale to some putative TDEs (with a blue featureless continuum), as well as to the mean of the recently identified population of ambiguous nuclear transients (ANTs), although the latter are spectroscopically distinct and tend to have longer timescales. ANTs have been speculated to arise from tidal disruptions of massive stars, perhaps in active galactic nuclei, and AT2023vto may represent a similar case but in a dormant black hole, thereby bridging the TDE and ANT populations. We anticipate that Rubin Observatory / LSST will uncover similar luminous TDEs to z ~ 3., Comment: 21 pages, 8 figures

We present multiwavelength observations of the Swift short $\gamma$-ray burst GRB 231117A, localized to an underlying galaxy at redshift $z = 0.257$ at a small projected offset ($\sim 2~$kpc). We uncover long-lived X-ray (Chandra) and radio/millimeter (VLA, MeerKAT, and ALMA) afterglow emission, detected to $\sim 37~$days and $\sim 20~$days (rest frame), respectively. We measure a wide jet ($\sim 10.4^\circ$) and relatively high circumburst density ($\sim 0.07~{\rm cm}^{-3}$) compared to the short GRB population. Our data cannot be easily fit with a standard forward shock model, but they are generally well fit with the incorporation of a refreshed forward shock and a reverse shock at $< 1~$day. We incorporate GRB 231117A into a larger sample of 132 X-ray detected events, 71 of which were radio-observed (17 cm-band detections), for a systematic study of the distributions of redshifts, jet and afterglow properties, galactocentric offsets, and local environments of events with and without detected radio afterglows. Compared to the entire short GRB population, the majority of radio-detected GRBs are at relatively low redshifts ($z < 0.6$) and have high circumburst densities ($> 10^{-2}~{\rm cm}^{-3}$), consistent with their smaller ($< 8~$kpc) projected galactocentric offsets. We additionally find that 70% of short GRBs with opening angle measurements were radio-detected, indicating the importance of radio afterglows in jet measurements, especially in the cases of wide ($> 10^\circ$) jets where observational evidence of collimation may only be detectable at radio wavelengths. Owing to improved observing strategies and the emergence of sensitive radio facilities, the number of radio-detected short GRBs has quadrupled in the past decade., Comment: 30 pages, 11 figures, submitted to ApJ

We present the most comprehensive catalog to date of Type I Superluminous Supernovae (SLSNe), a class of stripped envelope supernovae (SNe) characterized by exceptionally high luminosities. We have compiled a sample of 262 SLSNe reported through 2022 December 31. We verified the spectroscopic classification of each SLSN and collated an exhaustive data set of UV, optical and IR photometry from both publicly available data and our own FLEET observational follow-up program, totaling over 30,000 photometric detections. Using these data we derive observational parameters such as the peak absolute magnitudes, rise and decline timescales, as well as bolometric luminosities, temperature and photospheric radius evolution for all SLSNe. Additionally, we model all light curves using a hybrid model that includes contributions from both a magnetar central engine and the radioactive decay of $^{56}$Ni. We explore correlations among various physical and observational parameters, and recover the previously found relation between ejecta mass and magnetar spin, as well as the overall progenitor pre-explosion mass distribution with a peak at $\approx 6.5$ M$_\odot$. We find no significant redshift dependence for any parameter, and no evidence for distinct sub-types of SLSNe. We find that $< 3$\% of SLSNe are best fit with a significant contribution from radioactive decay $\gtrsim 50$\%, representing a set of relatively dim and slowly declining SNe. We provide several analytical tools designed to simulate typical SLSN light curves across a broad range of wavelengths and phases, enabling accurate K-corrections, bolometric scaling calculations, and inclusion of SLSNe in survey simulations or future comparison works. The complete catalog, including all of the photometry, models, and derived parameters, is made available as an open-source resource on GitHub., Comment: 59 pages, 22 Figures, Submitted to MNRAS

Gravitational wave detectors are observing compact object mergers from increasingly far distances, revealing the redshift evolution of the binary black hole (BBH) -- and soon the black hole-neutron star (BHNS) and binary neutron star (BNS) -- merger rate. To help interpret these observations, we investigate the expected redshift evolution of the compact object merger rate from the isolated binary evolution channel. We present a publicly available catalog of compact object mergers and their accompanying cosmological merger rates from population synthesis simulations conducted with the COMPAS software. To explore the impact of uncertainties in stellar and binary evolution, our simulations use two-parameter grids of binary evolution models that vary the common-envelope efficiency with mass transfer accretion efficiency, and supernova remnant mass prescription with supernova natal kick velocity, respectively. We quantify the redshift evolution of our simulated merger rates using the local ($z\sim 0$) rate, the redshift at which the merger rate peaks, and the normalized differential rates (as a proxy for slope). We find that although the local rates span a range of $\sim 10^3$ across our model variations, their redshift-evolutions are remarkably similar for BBHs, BHNSs, and BNSs, with differentials typically within a factor $3$ and peaks of $z\approx 1.2-2.4$ across models. Furthermore, several trends in our simulated rates are correlated with the model parameters we explore. We conclude that future observations of the redshift evolution of the compact object merger rate can help constrain binary models for stellar evolution and gravitational-wave formation channels., Comment: 16 pages, 9 figures, link to code https://github.com/Adam-Boesky/Exploring_Parameter_Space, link to data https://gwlandscape.org.au/compas/publication/Q29tcGFzUHVibGljYXRpb25Ob2RlOjYz/

Gravitational-wave detectors are now making it possible to investigate how the merger rate of binary black holes (BBHs) evolves with redshift. In this study, we examine whether the BBH merger rate of isolated binaries deviates from a scaled star formation rate density (SFRD) -- a frequently used model in state-of-the-art research. To address this question, we conduct population synthesis simulations using COMPAS with a grid of stellar evolution models, calculate their cosmological merger rates, and compare them to a scaled SFRD. We find that our simulated rates deviate by factors up to $3.5\times$ at $z\sim0$ and $5\times$ at $z\sim 9$ due to two main phenomena: (i) The formation efficiency of BBHs is an order of magnitude higher at low metallicities than at solar metallicity; and (ii) BBHs experience a wide range of delays (from a few Myr to many Gyr) between formation and merger. Deviations are similar when comparing to a $\textit{delayed}$ SFRD, and even larger (up to $\sim 10\times$) when comparing to SFRD-based models scaled to the local merger rate. Interestingly, our simulations find that the BBH delay time distribution is redshift-dependent, increasing the complexity of the redshift distribution of mergers. We find similar results for simulated merger rates of BHNSs and BNSs. We conclude that the rate of BBH, BHNS, and BNS mergers from the isolated channel can significantly deviate from a scaled SFRD, and that future measurements of the merger rate will provide insights into the formation pathways of gravitational-wave sources., Comment: 20 pages, 12 figures, link to code https://github.com/Adam-Boesky/Exploring_Parameter_Space, link to data https://gwlandscape.org.au/compas/publication/Q29tcGFzUHVibGljYXRpb25Ob2RlOjYz/

We present results from a search for radio emission in 77 stellar systems hosting 140 exoplanets, predominantly within 17.5 pc using the Very Large Array (VLA) at $4-8$ GHz. This is the largest and most sensitive search to date for radio emission in exoplanetary systems in the GHz frequency range. We obtained new observations of 58 systems, and analyzed archival observations of an additional 19 systems. Our choice of frequency and volume limit are motivated by radio detections of ultracool dwarfs (UCDs), including T dwarfs with masses at the exoplanet threshold of $\sim\!13\,M_J$. Our surveyed exoplanets span a mass range of $\approx\,10^{-3}-10\,M_J$ and semi-major axes of $\approx\,10^{-2}-10\,$AU. We detect a single target - GJ 3323 (M4) hosting two exoplanets with minimum masses of 2 and 2.3$\,M_\oplus$ - with a circular polarization fraction of $\approx\,40\%$; the radio luminosity agrees with its known X-ray luminosity and the G\"udel-Benz relation for stellar activity suggesting a likely stellar origin, but the high circular polarization fraction may also be indicative of star-planet interaction. For the remaining sources our $3\sigma$ upper limits are generally $L_\nu\lesssim\,10^{12.5}\,\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{Hz}^{-1}$, comparable to the lowest radio luminosities in UCDs. Our results are consistent with previous targeted searches of individual systems at GHz frequencies while greatly expanding the sample size. Our sensitivity is comparable to predicted fluxes for some systems considered candidates for detectable star-planet interaction. Observations with future instruments such as the Square Kilometer Array and Next Generation Very Large Array will be necessary to further constrain emission mechanisms from exoplanet systems at GHz frequencies., Comment: Submitted to ApJ, 18 pages, 8 figures

We present detailed radio observations of the tidal disruption event (TDE) ASASSN-19bt/AT2019ahk, obtained with the Australia Telescope Compact Array (ATCA), the Atacama Large Millimeter/submillimeter Array (ALMA), and the MeerKAT radio telescopes, spanning 40 to 1464 days after the onset of the optical flare. We find that ASASSN-19bt displays unusual radio evolution compared to other TDEs, as the peak brightness of its radio emission increases rapidly until 457 days post-optical discovery and then plateaus. Using a generalized approach to standard equipartition techniques, we estimate the energy and corresponding physical parameters for two possible emission geometries: a non-relativistic spherical outflow and a relativistic outflow observed from an arbitrary viewing angle. We find that the non-relativistic solution implies a continuous energy rise in the outflow from $E\sim10^{46}$ erg to $E\sim10^{49}$ erg with $\beta \approx 0.05$, while the off-axis relativistic jet solution instead suggests $E\approx10^{52}$ erg with $\Gamma\sim10$ erg at late times in the maximally off-axis case. We find that neither model provides a holistic explanation for the origin and evolution of the radio emission, emphasizing the need for more complex models. ASASSN-19bt joins the population of TDEs that display unusual radio emission at late times. Conducting long-term radio observations of these TDEs, especially during the later phases, will be crucial for understanding how these types of radio emission in TDEs are produced., Comment: 25 pages. Submitted to ApJ

Rowlinson et al. 2023 recently claimed the detection of a coherent radio flash 76.6 minutes after a short gamma-ray burst. They proposed that the radio emission may be associated with a long-lived neutron star engine. We show through theoretical and observational arguments that the coherent radio emission, if real and indeed associated with GRB 201006A and at the estimated redshift, is unlikely to be due to the collapse of the neutron star, ruling out a blitzar-like mechanism. Instead, we show if a long-lived engine was created, it must have been stable with the radio emission likely linked to the intrinsic magnetar activity. However, we find that the optical upper limits require fine-tuning to be consistent with a magnetar-driven kilonova: we show that neutron-star engines that do satisfy the optical constraints would have produced a bright kilonova afterglow that should already be observable by the VLA or MeerKAT (for ambient densities typical for short GRBs). Given the optical limits and the current lack of a kilonova afterglow, we instead posit that no neutron star survived the merger, and the coherent radio emission was produced far from a black hole central engine via mechanisms such as synchrotron maser or magnetic reconnection in the jet -- a scenario consistent with all observations. We encourage future radio follow-up to probe the engine of this exciting event and continued prompt radio follow-up of short GRBs., Comment: Submitted. 10 pages, 4 figures

Photometric classifications of supernova (SN) light curves have become necessary to utilize the full potential of large samples of observations obtained from wide-field photometric surveys, such as the Zwicky Transient Facility (ZTF) and the Vera C. Rubin Observatory. Here, we present a photometric classifier for SN light curves that does not rely on redshift information and still maintains comparable accuracy to redshift-dependent classifiers. Our new package, Superphot+, uses a parametric model to extract meaningful features from multiband SN light curves. We train a gradient-boosted machine with fit parameters from 6,061 ZTF SNe that pass data quality cuts and are spectroscopically classified as one of five classes: SN Ia, SN II, SN Ib/c, SN IIn, and SLSN-I. Without redshift information, our classifier yields a class-averaged F1-score of 0.61 +/- 0.02 and a total accuracy of 0.83 +/- 0.01. Including redshift information improves these metrics to 0.71 +/- 0.02 and 0.88 +/- 0.01, respectively. We assign new class probabilities to 3,558 ZTF transients that show SN-like characteristics (based on the ALeRCE Broker light curve and stamp classifiers), but lack spectroscopic classifications. Finally, we compare our predicted SN labels with those generated by the ALeRCE light curve classifier, finding that the two classifiers agree on photometric labels for 82 +/- 2% of light curves with spectroscopic labels and 72% of light curves without spectroscopic labels. Superphot+ is currently classifying ZTF SNe in real time via the ANTARES Broker, and is designed for simple adaptation to six-band Rubin light curves in the future., Comment: 37 pages, 25 figures. Submitted to AAS Journals

Astrometry from the Gaia mission was recently used to discover the two nearest known stellar-mass black holes (BHs), Gaia BH1 and Gaia BH2. Both systems contain $\sim 1\,M_{\odot}$ stars in wide orbits ($a\approx$1.4 AU, 4.96 AU) around $\sim9\,M_{\odot}$ BHs. These objects are among the first stellar-mass BHs not discovered via X-rays or gravitational waves. The companion stars -- a solar-type main sequence star in Gaia BH1 and a low-luminosity red giant in Gaia BH2 -- are well within their Roche lobes. However, the BHs are still expected to accrete stellar winds, leading to potentially detectable X-ray or radio emission. Here, we report observations of both systems with the Chandra X-ray Observatory and radio observations with the Very Large Array (for Gaia BH1) and MeerKAT (for Gaia BH2). We did not detect either system, leading to X-ray upper limits of $L_X < 10^{29.4}$ and $L_X < 10^{30.1}\,\rm erg\,s^{-1}$ and radio upper limits of $L_r < 10^{25.2}$ and $L_r < 10^{25.9}\,\rm erg\,s^{-1}$. For Gaia BH2, the non-detection implies that the the accretion rate near the horizon is much lower than the Bondi rate, consistent with recent models for hot accretion flows. We discuss implications of these non-detections for broader BH searches, concluding that it is unlikely that isolated BHs will be detected via ISM accretion in the near future. We also calculate evolutionary models for the binaries' future evolution using Modules for Experiments in Stellar Astrophysics (MESA). We find that Gaia BH1 will be X-ray bright for 5--50 Myr when the star is a red giant, including 5 Myr of stable Roche lobe overflow. Since no symbiotic BH X-ray binaries are known, this implies either that fewer than $\sim 10^4$ Gaia BH1-like binaries exist in the Milky Way, or that they are common but have evaded detection, perhaps due to very long outburst recurrence timescales., Comment: Submitted to PASP

Identifying the astrophysical sites of the $r$-process, one of the primary mechanisms by which heavy elements are formed, is a key goal of modern astrophysics. The discovery of the brightest gamma-ray burst of all time, GRB 221009A, at a relatively nearby redshift, presented the first opportunity to spectroscopically test the idea that $r$-process elements are produced following the collapse of rapidly rotating massive stars. Here we present spectroscopic and photometric $\textit{James Webb Space Telescope}$ (JWST) observations of GRB 221009A obtained $+168$ and $+170$ rest-frame days after the initial gamma-ray trigger, and demonstrate they are well-described by a supernova (SN) and power-law afterglow, with no evidence for an additional component from $r$-process emission, and that the SN component strongly resembles the near-infrared spectra of previous SNe, including SN 1998bw. We further find that the SN associated with GRB 221009A is slightly fainter than the expected brightness of SN 1998bw at this phase, concluding that the SN is therefore not an unusual GRB-SN. We infer a nickel mass of $\approx0.09$ M$_{\odot}$, consistent with the lack of an obvious SN detection in the early-time data. We find that the host galaxy of GRB 221009A has a very low metallicity of $\approx0.12$ Z$_{\odot}$ and our resolved host spectrum shows that GRB 221009A occurred in a unique environment in its host characterized by strong H$_2$ emission lines consistent with recent star formation, which may hint at environmental factors being responsible for its extreme energetics., Comment: 32 pages, 14 Figures, Submitted to Nature Astronomy

We present radio observations of 23 optically discovered tidal disruption events (TDEs) on timescales of 500-3200 days post discovery. We detect nine new TDEs that did not have detectable radio emission at earlier times, indicating a late-time brightening after several hundred (and up to 2300) days; an additional seven TDEs exhibit radio emission whose origin is ambiguous or may be attributed to the host galaxy or an active galactic nucleus. We also report a new rising component in one TDE previously detected in the radio at 10^3 days. While the radio emission in some of the detected TDEs peaked on a timescale 2-4 yr, over half of the sample still show rising emission. The range of luminosities for the sample is 10^37-10^39 erg/s, about 2 orders of magnitude below the radio luminosity of the relativistic TDE Sw J1644+57. Our data set indicates 40% of all optical TDEs are detected in radio hundreds to thousands of days after discovery, and that this is probably more common than early radio emission peaking at 10^2 days. Using an equipartition analysis, we find evidence for a delayed launch of the radio-emitting outflows, with delay timescales of 500-2000 days, inferred velocities of 0.02-0.15c, and kinetic energies of 10^47-10^49 erg. We rule out off axis relativistic jets as a viable explanation for this population, and conclude delayed outflows are a more likely explanation, possibly from delayed disk formation. We conclude late radio emission marks a fairly ubiquitous but heretofore overlooked phase of TDE evolution., Comment: 30 pages. Published in ApJ

We present the discovery of the radio afterglow of the short $\gamma$-ray burst (GRB) 210726A, localized to a galaxy at a photometric redshift of $z\sim 2.4$. While radio observations commenced $\lesssim 1~$day after the burst, no radio emission was detected until $\sim11~$days. The radio afterglow subsequently brightened by a factor of $\sim 3$ in the span of a week, followed by a rapid decay (a "radio flare"). We find that a forward shock afterglow model cannot self-consistently describe the multi-wavelength X-ray and radio data, and underpredicts the flux of the radio flare by a factor of $\approx 5$. We find that the addition of substantial energy injection, which increases the isotropic kinetic energy of the burst by a factor of $\approx 4$, or a reverse shock from a shell collision are viable solutions to match the broad-band behavior. At $z\sim 2.4$, GRB 210726A is among the highest redshift short GRBs discovered to date as well as the most luminous in radio and X-rays. Combining and comparing all previous radio afterglow observations of short GRBs, we find that the majority of published radio searches conclude by $\lesssim 10~$days after the burst, potentially missing these late rising, luminous radio afterglows., Comment: 29 pages, 10 figures, accepted to ApJ

We present an extensive $\textit{Hubble Space Telescope}$ ($\textit{HST}$) rest-frame ultraviolet (UV) imaging study of the locations of Type I superluminous supernovae (SLSNe) within their host galaxies. The sample includes 65 SLSNe with detected host galaxies in the redshift range $z\approx 0.05-2$. Using precise astrometric matching with SN images, we determine the distributions of physical and host-normalized offsets relative to the host centers, as well as the fractional flux distribution relative to the underlying UV light distribution. We find that the host-normalized offsets of SLSNe roughly track an exponential disk profile, but exhibit an overabundance of sources with large offsets of $1.5-4$ times their host half-light radius. The SLSNe normalized offsets are systematically larger than those of long gamma-ray bursts (LGRBs), and even Type Ib/c and II SNe. Furthermore, we find that about 40\% of all SLSNe occur in the dimmest regions of their host galaxies (fractional flux of 0), in stark contrast to LGRBs and Type Ib/c and II SNe. We do not detect any significant trends in the locations of SLSNe as a function of redshift, or as a function of explosion and magnetar engine parameters inferred from modeling of their optical lights curves. The significant difference in SLSN locations compared to LGRBs (and normal core-collapse SNe) suggests that at least some of their progenitors follow a different evolutionary path. We speculate that SLSNe arise from massive runaway stars from disrupted binary systems, with velocities of $\sim 10^2$ km s$^{-1}$., Comment: 31 pages, 14 figures, 5 tables. Submitted to ApJ. Comments welcomed

We present new spectroscopic and photometric follow-up observations of the known sample of extreme coronal line emitting galaxies (ECLEs) identified in the Sloan Digital Sky Survey (SDSS). With these new data, observations of the ECLE sample now span a period of two decades following their initial SDSS detections. We confirm the nonrecurrence of the iron coronal line signatures in five of the seven objects, further supporting their identification as the transient light echoes of tidal disruption events (TDEs). Photometric observations of these objects in optical bands show little overall evolution. In contrast, mid-infrared (MIR) observations show ongoing long-term declines. The remaining two objects had been classified as active galactic nuclei (AGN) with unusually strong coronal lines rather than being TDE related, given the persistence of the coronal lines in earlier follow-up spectra. We confirm this classification, with our spectra continuing to show the presence of strong, unchanged coronal-line features and AGN-like MIR colours and behaviour. We have constructed spectral templates of both subtypes of ECLE to aid in distinguishing the likely origin of newly discovered ECLEs. We highlight the need for higher cadence, and more rapid, follow-up observations of such objects to better constrain their properties and evolution. We also discuss the relationships between ECLEs, TDEs, and other identified transients having significant MIR variability., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. Note the corrected caption of Figure 1 continued, which in this version correctly refers to 'SDSS J124' rather than the erroneous 'SDSS J1341' in the published version. 29 Pages, 14 Figures

We present the discovery of the Type II supernova SN 2023ixf in M101 and follow-up photometric and spectroscopic observations, respectively, in the first month and week of its evolution. Our discovery was made within a day of estimated first light, and the following light curve is characterized by a rapid rise ($\approx5$ days) to a luminous peak ($M_V\approx-18.2$ mag) and plateau ($M_V\approx-17.6$ mag) extending to $30$ days with a fast decline rate of $\approx0.03$ mag day$^{-1}$. During the rising phase, $U-V$ color shows blueward evolution, followed by redward evolution in the plateau phase. Prominent flash features of hydrogen, helium, carbon, and nitrogen dominate the spectra up to $\approx5$ days after first light, with a transition to a higher ionization state in the first $\approx2$ days. Both the $U-V$ color and flash ionization states suggest a rise in the temperature, indicative of a delayed shock breakout inside dense circumstellar material (CSM). From the timescales of CSM interaction, we estimate its compact radial extent of $\sim(3-7)\times10^{14}$ cm. We then construct numerical light-curve models based on both continuous and eruptive mass-loss scenarios shortly before explosion. For the continuous mass-loss scenario, we infer a range of mass-loss history with $0.1-1.0\,M_\odot\,{\rm yr}^{-1}$ in the final $2-1$ yr before explosion, with a potentially decreasing mass loss of $0.01-0.1\,M_\odot\,{\rm yr}^{-1}$ in $\sim0.7-0.4$ yr toward the explosion. For the eruptive mass-loss scenario, we favor eruptions releasing $0.3-1\,M_\odot$ of the envelope at about a year before explosion, which result in CSM with mass and extent similar to the continuous scenario. We discuss the implications of the available multiwavelength constraints obtained thus far on the progenitor candidate and SN 2023ixf to our variable CSM models., Comment: Updated to match the published letter in ApJL, 2023 September 19

We present a detailed study on SN2019szu, a Type I superluminous supernova at $z=0.213$, that displayed unique photometric and spectroscopic properties. Pan-STARRS and ZTF forced photometry shows a pre-explosion plateau lasting $\sim$ 40 days. Unlike other SLSNe that show decreasing photospheric temperatures with time, the optical colours show an apparent temperature increase from $\sim$15000 K to $\sim$20000 K over the first 70 days, likely caused by an additional pseudo-continuum in the spectrum. Remarkably, the spectrum displays a forbidden emission line even during the rising phase of the light curve, inconsistent with an apparently compact photosphere. We show that this early feature is [O II] $\lambda\lambda$7320,7330. We also see evidence for [O III] $\lambda\lambda$4959, 5007, and [O III] $\lambda$4363 further strengthening this line identification. Comparing with models for nebular emission, we find that the oxygen line fluxes and ratios can be reproduced with $\sim$0.25 M$_{\odot}$ of oxygen rich material with a density of $\sim10^{-15} \rm{g cm}^{-3}$. The low density suggests a circumstellar origin, but the early onset of the emission lines requires that this material was ejected within the final months before the terminal explosion, consistent with the timing of the precursor plateau. Interaction with denser material closer to the explosion likely produced the pseudo-continuum bluewards of $\sim$5500 \AA. We suggest that this event is one of the best candidates to date for a pulsational pair-instability ejection, with early pulses providing the low density material needed for the forbidden emission line, and collisions between the final shells of ejected material producing the pre-explosion plateau., Comment: Accepted for publication in MNRAS

Type-I Superluminous Supernovae (SLSNe) are an exotic class of core-collapse SN (CCSN) that can be up to 100 times brighter and more slowly-evolving than normal CCSNe. SLSNe represent the end-stages of the most massive stripped stars, and are thought to be powered by the spin-down energy of a millisecond magnetar. Studying them and measuring their physical parameters can help us to better understand stellar mass-loss, evolution, and explosions. Moreover, thanks to their high luminosities, SLSNe can be seen up to greater distances, allowing us to explore how stellar physics evolves as a function of redshift. The High Latitude Time Domain Survey (HLTDS) will provide us with an exquisite dataset that will discover 100s of SLSNe. Here, we focus on the question of which sets of filters and cadences will allow us to best characterize the physical parameters of these SLSNe. We simulate a set of SLSNe at redshifts ranging from z = 0.1 to z = 5.0, using six different sets of filters, and cadences ranging from 5 to 100 days. We then fit these simulated light curves to attempt to recover the input parameter values for their ejecta mass, ejecta velocity, magnetic field strength, and magnetar spin period. We find that four filters are sufficient to accurately characterize SLSNe at redshifts below $z = 3$, and that cadences faster than 20 days are required to obtain measurements with an uncertainty below 10\%, although a cadence of 70 days is still acceptable under certain conditions. Finally, we find that the nominal survey strategy will not be able to properly characterize the most distant SLSNe at $z = 5$. We find that the addition of 60-day cadence observations for 4 years to the nominal HLTDS survey can greatly improve the prospect of characterizing these most extreme and distant SNe, with only an 8\% increase to the time commitment of the survey., Comment: Submitted as part of the Nancy Grace Roman Space Telescope's Core Community Surveys call for white papers

Gravitational-wave astronomy has revolutionized humanity's view of the universe, a revolution driven by observations that no other field can make. This white paper describes an observatory that builds on decades of investment by the National Science Foundation and that will drive discovery for decades to come: Cosmic Explorer. Major discoveries in astronomy are driven by three related improvements: better sensitivity, higher precision, and opening new observational windows. Cosmic Explorer promises all three and will deliver an order-of-magnitude greater sensitivity than LIGO. Cosmic Explorer will push the gravitational-wave frontier to almost the edge of the observable universe using technologies that have been proven by LIGO during its development. With the unprecedented sensitivity that only a new facility can deliver, Cosmic Explorer will make discoveries that cannot yet be anticipated, especially since gravitational waves are both synergistic with electromagnetic observations and can reach into regions of the universe that electromagnetic observations cannot explore. With Cosmic Explorer, scientists can use the universe as a laboratory to test the laws of physics and study the nature of matter. Cosmic Explorer allows the United States to continue its leading role in gravitational-wave science and the international network of next-generation observatories. With its extraordinary discovery potential, Cosmic Explorer will deliver revolutionary observations across astronomy, physics, and cosmology including: Black Holes and Neutron Stars Throughout Cosmic Time, Multi-Messenger Astrophysics and Dynamics of Dense Matter, New Probes of Extreme Astrophysics, Fundamental Physics and Precision Cosmology, Dark Matter and the Early Universe.

We present 1.3 mm (230 GHz) observations of the recent and nearby Type II supernova, SN2023ixf, obtained with the Submillimeter Array (SMA) at 2.6-18.6 days after explosion. The observations were obtained as part the SMA Large Program POETS (Pursuit of Extragalactic Transients with the SMA). We do not detect any emission at the location of SN2023ixf, with the deepest limits of $L_\nu(230\,{\rm GHz})\lesssim 8.6\times 10^{25}$ erg s$^{-1}$ Hz$^{-1}$ at 2.7 and 7.7 days, and $L_\nu(230\,{\rm GHz})\lesssim 3.4\times 10^{25}$ erg s$^{-1}$ Hz$^{-1}$ at 18.6 days. These limits are about a factor of 2 times dimmer than the mm emission from SN2011dh (IIb), about an order of magnitude dimmer compared to SN1993J (IIb) and SN2018ivc (IIL), and about 30 times dimmer than the most luminous non-relativistic SNe in the mm-band (Type IIb/Ib/Ic). Using these limits in the context of analytical models that include synchrotron self-absorption and free-free absorption we place constraints on the proximate circumstellar medium around the progenitor star, to a scale of $\sim 2\times 10^{15}$ cm, excluding the range $\dot{M}\sim {\rm few}\times 10^{-6}-10^{-2}$ M$_\odot$ yr$^{-1}$ (for a wind velocity, $v_w=115$ km s$^{-1}$, and ejecta velocity, $v_{\rm eje}\sim (1-2)\times 10^4$ km s$^{-1}$). These results are consistent with an inference of the mass loss rate based on optical spectroscopy ($\sim 2\times 10^{-2}$ M$_\odot$ yr$^{-1}$ for $v_w=115$ km s$^{-1}$), but are in tension with the inference from hard X-rays ($\sim 7\times 10^{-4}$ M$_\odot$ yr$^{-1}$ for $v_w=115$ km s$^{-1}$). This tension may be alleviated by a non-homogeneous and confined CSM, consistent with results from high-resolution optical spectroscopy., Comment: Submitted

We present optical photometry and spectroscopy of the Type IIn supernova (SN) 2021qqp. Its unusual light curve is marked by a long precursor for $\approx300$ days, a rapid increase in brightness for $\approx60$ days, and then a sharp increase of $\approx1.6$ mag in only a few days to a first peak of $M_r \approx -19.5$ mag. The light curve then declines rapidly until it re-brightens to a second distinct peak of $M_r \approx -17.3$ mag centered at $\approx335$ days after the first peak. The spectra are dominated by Balmer lines with a complex morphology, including a narrow component with a width of $\approx 1300$ km s$^{-1}$ (first peak) and $\approx 2500$ km s$^{-1}$ (second peak) that we associate with the circumstellar medium (CSM) and a P Cygni component with an absorption velocity of $\approx 8500$ km s$^{-1}$ (first peak) and $\approx 5600$ km s$^{-1}$ (second peak) that we associate with the SN-CSM interaction shell. Using the luminosity and velocity evolution, we construct a flexible analytical model, finding two significant mass-loss episodes with peak mass loss rates of $\approx 10$ and $\approx 5\,M_{\odot}$ yr$^{-1}$ about $0.8$ and $2$ yr before explosion, respectively, with a total CSM mass of $\approx 2-4\,M_{\odot}$. We show that the most recent mass-loss episode could explain the precursor for the year preceding the explosion. The SN ejecta mass is constrained to be $\approx 5-30\,M_{\odot}$ for an explosion energy of $\approx (3-10)\times10^{51}$ erg. We discuss eruptive massive stars (luminous blue variable, pulsational pair instability) and an extreme stellar merger with a compact object as possible progenitor channels., Comment: Updated to match the published article in ApJ, 2024 March 29

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

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

We present the largest compilation to date of optical observations during and following fast radio bursts (FRBs). The data set includes our dedicated simultaneous and follow-up observations, as well as serendipitous archival survey observations, for a sample of 15 well-localized FRBs: eight repeating and seven one-off sources. Our simultaneous (and nearly simultaneous with a $0.4$ s delay) optical observations of 13 (1) bursts from the repeating FRB 20220912A provide the deepest such limits to date for any extragalactic FRB, reaching a luminosity limit of $\nu L_\nu\lesssim 10^{42}$ erg s$^{-1}$ ($\lesssim 2\times10^{41}$ erg s$^{-1}$) with $15-400$ s exposures; an optical-flux-to-radio-fluence ratio of $f_{\rm opt}/F_{\rm radio}\lesssim 10^{-7}$ ms$^{-1}$ ($\lesssim 10^{-8}$ ms$^{-1}$); and flux ratio of $f_{\rm opt}/f_{\rm radio}\lesssim 0.02-\lesssim 2\times 10^{-5}$ ($\lesssim 10^{-6}$) on millisecond to second timescales. These simultaneous limits provide useful constraints in the context of FRB emission models, such as the pulsar magnetosphere and pulsar nebula models. Interpreting all available optical limits in the context of the synchrotron maser model, we find that they constrain the flare energies to $\lesssim 10^{43}-10^{49}$ erg (depending on the distances of the various repeating FRBs, with $\lesssim 10^{39}$ erg for the Galactic SGR 1935+2154). These limits are generally at least an order of magnitude larger than those inferred from the FRBs themselves, although in the case of FRB 20220912A our simultaneous and rapid follow-up observations severely restrict the model parameter space. We conclude by exploring the potential of future simultaneous and rapid-response observations with large optical telescopes., Comment: Updated to match the published letter in ApJL, 2023 April 26

We present a detailed compilation and analysis of the X-ray phase space of low- to intermediate-redshift ($ 0\le z \le 1$) transients that consolidates observed light curves (and theory where necessary) for a large variety of classes of transient/variable phenomena in the 0.3--10 keV energy band. We include gamma-ray burst afterglows, supernovae, supernova shock breakouts and shocks interacting with the environment, tidal disruption events and active galactic nuclei, fast blue optical transients, cataclysmic variables, magnetar flares/outbursts and fast radio bursts, cool stellar flares, X-ray binary outbursts, and ultraluminous X-ray sources. Our overarching goal is to offer a comprehensive resource for the examination of these ephemeral events, extending the X-ray duration-luminosity phase space (DLPS) to show luminosity evolution. We use existing observations (both targeted and serendipitous) to characterize the behavior of various transient/variable populations. Contextualizing transient signals in the larger DLPS serves two primary purposes: to identify areas of interest (i.e., regions in the parameter space where one would expect detections, but in which observations have historically been lacking) and to provide initial qualitative guidance in classifying newly discovered transient signals. We find that while the most luminous (largely extragalactic) and least luminous (largely Galactic) part of the phase space is well-populated at $t > 0.1$ days, intermediate luminosity phenomena (L$_x = 10^{34} - 10^{42}$ erg s$^{-1}$) represent a gap in the phase space. We thus identify L$_x = 10^{34} - 10^{42}$ erg s$^{-1}$ and $t = 10^{-4} - 0.1$ days as a key discovery phase space in transient X-ray astronomy., Comment: 12 figures, 13 tables; version accepted to ApJ

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., Comment: 9 pages, 5 figures. Accepted for publication in Physics Letters B

In November 2019 we began operating FLEET (Finding Luminous and Exotic Extragalactic Transients), a machine learning algorithm designed to photometrically identify Type I superluminous supernovae (SLSNe) in transient alert streams. Using FLEET, we spectroscopically classified 21 of the 50 SLSNe identified worldwide between November 2019 and January 2022. Based on our original algorithm, we anticipated that FLEET would achieve a purity of about 50\% for transients with a probability of being a SLSN, \pslsn$>0.5$; the true on-sky purity we obtained is closer to 80\%. Similarly, we anticipated FLEET could reach a completeness of about 30\%, and we indeed measure an upper limit on the completeness of $\approx 33$\%. Here, we present FLEET 2.0, an updated version of FLEET trained on 4,780 transients (almost 3 times more than in FLEET 1.0). FLEET 2.0 has a similar predicted purity to FLEET 1.0, but outperforms FLEET 1.0 in terms of completeness, which is now closer to $\approx 40$\% for transients with \pslsn$>0.5$. Additionally, we explore possible systematics that might arise from the use of FLEET for target selection. We find that the population of SLSNe recovered by FLEET is mostly indistinguishable from the overall SLSN population, in terms of physical and most observational parameters. We provide FLEET as an open source package on GitHub https://github.com/gmzsebastian/FLEET, Comment: 10 pages, 3 figures, submitted to ApJ

We present an expansion of FLEET, a machine learning algorithm optimized to select transients that are most likely to be tidal disruption events (TDEs). FLEET is based on a random forest algorithm trained on the light curves and host galaxy information of 4,779 spectroscopically classified transients. For transients with a probability of being a TDE, \ptde$>0.5$, we can successfully recover TDEs with a $\approx40$\% completeness and a $\approx30$\% purity when using the first 20 days of photometry, or a similar completeness and $\approx50$\% purity when including 40 days of photometry. We find that the most relevant features for differentiating TDEs from other transients are the normalized host separation, and the light curve $(g-r)$ color during peak. Additionally, we use FLEET to produce a list of the 39 most likely TDE candidates discovered by the Zwicky Transient Facility that remain currently unclassified. We explore the use of FLEET for the Legacy Survey of Space and Time on the Vera C. Rubin Observatory (\textit{Rubin}) and the \textit{Nancy Grace Roman Space Telescope} (\textit{Roman}). We simulate the \textit{Rubin} and \textit{Roman} survey strategies and estimate that $\sim 10^4$ TDEs could be discovered every year by \textit{Rubin}, and $\sim200$ TDEs per year by \textit{Roman}. Finally, we run FLEET on the TDEs in our \textit{Rubin} survey simulation and find that we can recover $\sim 30$\% of those at a redshift $z <0.5$ with \ptde$>0.5$. This translates to $\sim3,000$ TDEs per year that FLEET could uncover from \textit{Rubin}. FLEET is provided as a open source package on GitHub https://github.com/gmzsebastian/FLEET, Comment: 17 pages, 12 figures, submitted to ApJ

We present a comprehensive study of 29 short gamma-ray bursts (SGRBs) observed $\approx 0.8-60$ days post-burst using $Chandra$ and $XMM-Newton$. We provide the inferred distributions of SGRB jet opening angles and true event rates to compare against neutron star merger rates. We perform uniform analysis and modeling of their afterglows, obtaining 10 opening angle measurements and 19 lower limits. We report on two new opening angle measurements (SGRBs 050724A and 200411A) and eight updated values, obtaining a median value of $\langle \theta_{\rm j} \rangle \approx 6.1^{\circ}$ [-3.2$^{\circ}$,+9.3$^{\circ}$] (68\% confidence on the full distribution) from jet measurements alone. For the remaining events, we infer $\theta_{\rm j}\gtrsim 0.5-26^{\circ}$. We uncover a population of SGRBs with wider jets of $\theta_{\rm j} \gtrsim 10^{\circ}$ (including two measurements of $\theta_{\rm j} \gtrsim 15^{\circ}$), representing $\sim 28\%$ of our sample. Coupled with multi-wavelength afterglow information, we derive a total true energy of $\langle E_{\rm true, tot} \rangle \approx 10^{49}-10^{50}$\,erg which is consistent with MHD jet launching mechanisms. Furthermore, we determine a range for the beaming-corrected event rate of $\mathfrak{R}_{\rm true} \approx360-1800$ Gpc$^{-3}$ yr$^{-1}$, set by the inclusion of a population of wide jets on the low end, and the jet measurements alone on the high end. From a comparison with the latest merger rates, our results are consistent with the majority of SGRBs originating from binary neutron star mergers. However, our inferred rates are well above the latest neutron star-black hole merger rates, consistent with at most a small fraction of SGRBs originating from such mergers., Comment: 32 pages, 8 figures, 4 tables (Submitted to ApJ)

We present late-time radio/millimeter (as well as optical/UV and X-ray) detections of the tidal disruption event (TDE) AT2018hyz, spanning $970 - 1300$ d after optical discovery. In conjunction with earlier deeper limits, including at $\approx 700$ d, our observations reveal rapidly rising emission at $0.8-240$ GHz, steeper than $F_\nu\propto t^5$ relative to the time of optical discovery. Such a steep rise cannot be explained in any reasonable scenario of an outflow launched at the time of disruption (e.g., off-axis jet, sudden increase in the ambient density), and instead points to a delayed launch. Our multi-frequency data allow us to directly determine the radius and energy of the radio-emitting outflow, showing that it was launched $\approx 750$ d after optical discovery. The outflow velocity is mildly relativistic, with $\beta\approx 0.25$ and $\approx 0.6$ for a spherical and a $10^\circ$ jet geometry, respectively, and the minimum kinetic energy is $E_K\approx 5.8\times 10^{49}$ and $\approx 6.3\times 10^{49}$ erg, respectively. This is the first definitive evidence for the production of a delayed mildly-relativistic outflow in a TDE; a comparison to the recently-published radio light curve of ASASSN-15oi suggests that the final re-brightening observed in that event (at a single frequency and time) may be due to a similar outflow with a comparable velocity and energy. Finally, we note that the energy and velocity of the delayed outflow in AT2018hyz are intermediate between those of past non-relativistic TDEs (e.g., ASASSN-14li, AT2019dsg) and the relativistic TDE Sw\,J1644+57. We suggest that such delayed outflows may be common in TDEs., Comment: Submitted to ApJ

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

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., Comment: 32 pages, 15 figures, 3 tables, published with ApJ

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

We present cm-band and mm-band afterglow observations of five long-duration $\gamma$-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 $A_{V, \rm GRB} \gtrsim 2.2 - 10.6~{\rm mag}$. Additionally, we model the host galaxies of the six bursts in our sample, and derive host galaxy dust extinctions of $A_{V, \rm Host} \approx 0.3-4.7~{\rm mag}$. Across all tested $\gamma$-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 $A_{V, \rm 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., Comment: 40 pages, 15 figures, submitted to ApJ

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

Stripped-envelope core-collapse supernovae can be divided into two broad classes: the common Type Ib/c supernovae (SNe Ib/c), powered by the radioactive decay of $^{56}$Ni, and the rare superluminous supernovae (SLSNe), most likely powered by the spin-down of a magnetar central engine. Up to now, the intermediate regime between these two populations has remained mostly unexplored. Here, we present a comprehensive study of 40 \textit{luminous supernovae} (LSNe), SNe with peak magnitudes of $M_r = -19$ to $-20$ mag, bound by SLSNe on the bright end and by SNe Ib/c on the dim end. Spectroscopically, LSNe appear to form a continuum between Type Ic SNe and SLSNe. Given their intermediate nature, we model the light curves of all LSNe using a combined magnetar plus radioactive decay model and find that they are indeed intermediate, not only in terms of their peak luminosity and spectra, but also in their rise times, power sources, and physical parameters. We sub-classify LSNe into distinct groups that are either as fast-evolving as SNe Ib/c or as slow-evolving as SLSNe, and appear to be either radioactively or magnetar powered, respectively. Our findings indicate that LSNe are powered by either an over-abundant production of $^{56}$Ni or by weak magnetar engines, and may serve as the missing link between the two populations., Comment: 39 pages, 16 figures, submitted to ApJ

The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community., Comment: 29 pages, 12 figures, Table of Implementation Recommendations, Appendix of Community Science Pitches, AURA-commissioned whitepaper submitted to the Director of STScI (Ken Sembach) and the Director of NOIRLab (Pat McCarthy)

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., Comment: 13 pages, 9 figures, submitted version after addressing the referee's comments

The advent of sensitive gravitational wave (GW) detectors, coupled with wide-field, high cadence optical time-domain surveys, raises the possibility of the first joint GW-electromagnetic (EM) detections of core-collapse supernovae (CCSNe). For targeted searches of GWs from CCSNe optical observations can be used to increase the sensitivity of the search by restricting the relevant time interval, defined here as the GW search window (GSW). The extent of the GSW is a critical factor in determining the achievable false alarm probability (FAP) for a triggered CCSN search. The ability to constrain the GSW from optical observations depends on how early a CCSN is detected, as well as the ability to model the early optical emission. Here we present several approaches to constrain the GSW, ranging in complexity from model-independent analytical fits of the early light curve, model-dependent fits of the rising or entire light curve, and a new data-driven approach using existing well-sampled CCSN light curves from {\it Kepler} and the Transiting Exoplanet Survey Satellite (TESS). We use these approaches to determine the time of core-collapse and its associated uncertainty (i.e., the GSW). We apply our methods to two Type II SNe that occurred during LIGO/Virgo Observing Run 3: SN\,2019fcn and SN\,2019ejj (both in the same galaxy at $d=15.7$ Mpc). Our approach shortens the duration of the GSW and improves the robustness of the GSW compared to techniques used in past GW CCSN searches.

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

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$ with a peak absolute magnitude of $M_r\approx -22.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_{\rm ej}\approx 3.8$ M$_\odot$, a magnetar spin period of $P\approx 2.9$ ms and a magnetic field strength of $B_{\perp}\approx 2.8\times 10^{14}$ 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\approx -17.85$, ($L_B \approx 0.029$ $L^*$), and an inferred metallicity of $Z\approx 0.3$ Z$_\odot$, and star formation rate of $\approx 0.8$ M$_\odot$ yr$^{-1}$.

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

We present the first systematic search for GHz frequency radio emission from directly imaged exoplanets using Very Large Array (VLA) observations of sufficient angular resolution to separate the planets from their host stars. We obtained results for five systems and eight exoplanets located at $\lesssim 50$ pc, through new observations (Ross 458, GU Psc, and 51 Eri) and archival data (GJ 504 and HR 8799). We do not detect radio emission from any of the exoplanets, with $3\sigma$ luminosity upper limits of $(0.9-23)\times10^{21}$ erg s$^{-1}$. These limits are comparable to the level of radio emission detected in several ultracool dwarfs, including T dwarfs, whose masses are only a factor of two times higher than those of the directly-imaged exoplanets. Despite the lack of detections in this pilot study, we highlight the need for continued GHz frequency radio observations of nearby exoplanets at $\mu$Jy-level sensitivity., Comment: Submitted to ApJ 7 pages, 4 figures

The light curve diversity of hydrogen-poor superluminous supernovae (SLSNe) has kept open the possibility that multiple power sources account for the population. Specifically, pair-instability explosions (PISNe), which produce large masses of $^{56}$Ni, have been argued as the origin of some slowly-evolving SLSNe. Here we present detailed observations of SN 2016inl (=PS16fgt), a slowly-evolving SLSN at $z=0.3057$, whose unusually red spectrum matches PS1-14bj, a SLSN with an exceptionally long rise time consistent with a PISN. Ground-based and Hubble Space Telescope data, spanning about 800 rest-frame days, reveal a significant light curve flattening, similar to that seen in SN 2015bn, and much slower than the decline rate expected from radioactive decay of $^{56}$Co. We therefore conclude that despite its slow evolution, SN 2016inl is inconsistent with a PISN. Instead, the light curve evolution matches the expected power-law spin-down of a magnetar central engine, but with a shallower power law ($L\propto t^{-2.8}$) compared to that in SN 2015bn, indicating a possible difference in the $\gamma$-ray opacity between the two events. Analytical modeling indicates typical magnetar engine parameters, but one of the highest ejecta masses ($\approx 20$ M$_{\odot}$) inferred for a SLSN. Our results indicate that monitoring the late-time light curve evolution of SLSNe provides a powerful diagnostic of their energy source., Comment: 16 pages, 10 figures, Submitted 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

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., Comment: 13 pages,9 figures, submitted to AAS Journals. Comments welcome!