Search

Quasi-periodic eruptions (QPEs) are an extreme X-ray variability phenomenon associated with low-mass supermassive black holes. First discovered in the nucleus of the galaxy GSN 069, they have been so far securely detected in five other galaxies, including RX J1301.9+2747. When detected, the out-of-QPE emission (quiescence) is consistent with the high-energy tail of thermal emission from an accretion disk. We present the X-ray and radio properties of RX J1301.9+2747, both in quiescence and during QPEs. We analyse X-ray data taken during five XMM-Newton observations between 2000 and 2022. The last three observations were taken in coordination with radio observations with the Karl G. Jansky Very Large Array. We also make use of EXOSAT, ROSAT, and Chandra archival observations taken between 1983 and 2009. XMM-Newton detected 34 QPEs of which 8 have significantly lower amplitudes than the others. No correlated radio/X-ray variability was observed during QPEs. In terms of timing properties, the QPEs in RX J1301.9+2747 do not exhibit the striking regularity observed in the discovery source GSN 069. In fact there is no clear repetition pattern between QPEs: the average time separation between their peaks is about four hours, but it can be as short as one, and as long as six hours. The QPE spectral properties of RX J1301.9+2747 as a function of energy are however very similar to those of GSN 069 and of other QPE sources. The quiescent emission of RX J1301.9+2747 is more complex than that of GSN 069, as it requires a soft X-ray excess-like component in addition to the thermal emission from the accretion disk. Its long-term X-ray quiescent flux variations are of low-amplitude and not strictly monotonic, with a general decay over $\sim 22$ years. We discuss our observational results in terms of some of the ideas and models that have been proposed so far for the physical origin of QPEs., Comment: Accepted for publication in Astronomy & Astrophysics. 22 pages and 17 figures in the main text; 11 pages, 5 tables, and 4 figures in the Appendix. Abstract shortened to comply with arXiv requirements

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 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

Galactic nuclei showing recurrent phases of activity and quiescence have recently been discovered, with recurrence times as short as a few hours to a day -- known as quasi-periodic X-ray eruption (QPE) sources -- to as long as hundreds to a thousand days for repeating nuclear transients (RNTs). Here we present a multi-wavelength overview of Swift J023017.0+283603 (hereafter Swift J0230+28), a source that exhibits repeating and quasi-periodic X-ray flares from the nucleus of a previously unremarkable galaxy at $\sim$ 165 Mpc, with a recurrence time of approximately 22 days, an intermediary timescale between known RNTs and QPE sources. The source also shows transient radio emission, likely associated with the X-ray emission. Such recurrent soft X-ray eruptions, with no accompanying UV/optical emission, are strikingly similar to QPE sources. However, in addition to having a recurrence time that is $\sim 25$ times longer than the longest-known QPE source, Swift J0230+28's eruptions exhibit somewhat distinct shapes and temperature evolution than the known QPE sources. Scenarios involving extreme mass ratio inspirals are favored over disk instability models. The source reveals an unexplored timescale for repeating extragalactic transients and highlights the need for a wide-field, time-domain X-ray mission to explore the parameter space of recurring X-ray transients., Comment: Final version, appeared on Nature Astronomy on 12 January 2024

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

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

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 report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-bolometric luminosity of log (L$_{\text{max}} / [\text{erg}/\text{s}]$) = 45.7 from slowly fading ztf-$\it{g}$ and ztf-$r$ light curves spanning over 1000 observer-frame days. The host galaxy is not detected in archival Pan-STARRS observations ($g > 23.3$ mag), implying a lower limit to the outburst amplitude of more than 5 mag relative to the quiescent host galaxy. Optical spectra from Lick and Keck Observatories exhibit strong emission lines with narrow cores from the H Balmer series and ultraviolet semi-forbidden lines of Si III] $\lambda$1892, C III] $\lambda$1909, and C II] $\lambda$2325. Typical nebular lines in AGN spectra from ions such as [O II] and [O III] are not detected. These spectral features, along with the smooth light curve that is unlike most AGN flaring activity, and the luminosity that exceeds any observed or theorized supernova, lead us to conclude that AT 2021lwx is most likely an extreme tidal disruption event (TDE). Modeling of ZTF photometry with MOSFiT suggests that the TDE was between a $\approx 14 M_{\odot}$ star and a supermassive black hole of mass $M_{\text{BH}} \sim$ $10^{8} M_{\odot}$. Continued monitoring of the still-evolving light curve along with deep imaging of the field once AT 2021lwx has faded can test this hypothesis and potentially detect the host galaxy., Comment: 15 pages, 4 figures, 1 Table; Version as published in The Astrophysical Journal Letters. Observations of AT 2021lwx published in the paper can be found at https://bsubraya.github.io/research/

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 photometric and spectroscopic observations of the extraordinary gamma-ray burst (GRB) 221009A in search of an associated supernova. Some past GRBs have shown bumps in the optical light curve that coincide with the emergence of supernova spectral features, but we do not detect any significant light curve features in GRB~221009A, nor do we detect any clear sign of supernova spectral features. Using two well-studied GRB-associated supernovae (SN~2013dx, $M_{r,max} = -19.54$; SN~2016jca, $M_{r,max} = -19.04$) at a similar redshift as GRB~221009A ($z=0.151$), we modeled how the emergence of a supernova would affect the light curve. If we assume the GRB afterglow to decay at the same rate as the X-ray data, the combination of afterglow and a supernova component is fainter than the observed GRB brightness. For the case where we assume the best-fit power law to the optical data as the GRB afterglow component, a supernova contribution should have created a clear bump in the light curve, assuming only extinction from the Milky Way. If we assume a higher extinction of $E(B-V)$=$1.74$ mag (as has been suggested elsewhere), the supernova contribution would have been hard to detect, with a limit on the associated supernova of $M_{r,max} \approx-$19.54. We do not observe any clear supernova features in our spectra, which were taken around the time of expected maximum light. The lack of a bright supernova associated with GRB~221009A may indicate that the energy from the explosion is mostly concentrated in the jet, leaving a lower energy budget available for the supernova., Comment: 19 pages, 8 figures, 4 tables, accepted for publication in ApJL

A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to understand the radiative mechanisms operating in super-Eddington jets. Here, we present multi-wavelength (X-ray, UV, optical, and radio) observations of the optically discovered transient \target at $z=1.193$. Its unusual X-ray properties, including a peak observed luminosity of $\gtrsim$10$^{48}$ erg s$^{-1}$, systematic variability on timescales as short as 1000 seconds, and overall duration lasting more than 30 days in the rest-frame are traits associated with relativistic TDEs. The X-ray to radio spectral energy distributions spanning 5-50 days after discovery can be explained as synchrotron emission from a relativistic jet (radio), synchrotron self-Compton (X-rays), and thermal emission similar to that seen in low-redshift TDEs (UV/optical). Our modeling implies a beamed, highly relativistic jet akin to blazars but requires extreme matter-domination, i.e, high ratio of electron-to-magnetic field energy densities in the jet, and challenges our theoretical understanding of jets., Comment: To appear in Nature Astronomy on 30th November 2022. Also see here for an animation explaining the result: https://youtu.be/MQHdSbxuznY

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 X-ray, UV, optical, and radio observations of the nearby ($\approx78$ Mpc) tidal disruption event (TDE) AT2021ehb/ZTF21aanxhjv during its first 430 days of evolution. AT2021ehb occurs in the nucleus of a galaxy hosting a $\approx 10^{7}\,M_\odot$ black hole ($M_{\rm BH}$ inferred from host galaxy scaling relations). High-cadence Swift and NICER monitoring reveals a delayed X-ray brightening. The spectrum first undergoes a gradual ${\rm soft }\rightarrow{\rm hard}$ transition and then suddenly turns soft again within 3 days at $\delta t\approx 272$ days during which the X-ray flux drops by a factor of ten. In the joint NICER+NuSTAR observation ($\delta t =264$ days, harder state), we observe a prominent non-thermal component up to 30 keV and an extremely broad emission line in the iron K band. The bolometric luminosity of AT2021ehb reaches a maximum of $6.0^{+10.4}_{-3.8}\% L_{\rm Edd}$ when the X-ray spectrum is the hardest. During the dramatic X-ray evolution, no radio emission is detected, the UV/optical luminosity stays relatively constant, and the optical spectra are featureless. We propose the following interpretations: (i) the ${\rm soft }\rightarrow{\rm hard}$ transition may be caused by the gradual formation of a magnetically dominated corona; (ii) hard X-ray photons escape from the system along solid angles with low scattering optical depth ($\sim\,$a few) whereas the UV/optical emission is likely generated by reprocessing materials with much larger column density -- the system is highly aspherical; (iii) the abrupt X-ray flux drop may be triggered by the thermal-viscous instability in the inner accretion flow leading to a much thinner disk., Comment: 35 pages, 21 figures, accepted by 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 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

The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of neutron star mergers and other gravitational wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving neutron stars (about tens per year) out to distances of several hundred Mpc. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of neutron star and other compact object mergers, and yet unknown classes of gravitational wave events., Comment: Accepted for publication in ApJS. arXiv admin note: text overlap with arXiv:1812.04051

We present deep X-ray and radio observations of the Fast Blue Optical Transient (FBOT) AT2020xnd/ZTF20acigmel at $z=0.2433$ from $13$d to $269$d after explosion. AT2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT2018cow. AT2020xnd shows luminous radio emission reaching $L_{\nu}\approx8\times10^{29}$ergs$^{-1}$Hz$^{-1}$ at 20GHz and $75$d post explosion, accompanied by luminous and rapidly fading soft X-ray emission peaking at $L_{X}\approx6\times10^{42}$ergs$^{-1}$. Interpreting the radio emission in the context of synchrotron radiation from the explosion's shock interaction with the environment we find that AT2020xnd launched a high-velocity outflow ($v\sim$0.1-0.2$c$) propagating into a dense circumstellar medium (effective $\dot M\approx10^{-3}M_{\rm{sol}}$yr$^{-1}$ for an assumed wind velocity of $v_w=1000$kms$^{-1}$). Similar to AT2018cow, the detected X-ray emission is in excess compared to the extrapolated synchrotron spectrum and constitutes a different emission component, possibly powered by accretion onto a newly formed black hole or neutron star. These properties make AT2020xnd a high-redshift analog to AT2018cow, and establish AT2020xnd as the fourth member of the class of optically-luminous FBOTs with luminous multi-wavelength counterparts., Comment: 21 pages, 6 figures, 6 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

Massive, rapidly-spinning magnetar remnants produced as a result of binary neutron star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a synchrotron radio signal from the interaction of the ejecta with the circumburst medium. We present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (SGRBs; $z<0.5$) on rest-frame timescales of $\approx2.4-13.9$ yr following the bursts. We place $3\sigma$ limits on radio continuum emission of $F_{\nu}\lesssim6-20\,\mu$Jy at the burst positions, or $L_{\nu}\lesssim(0.6-8.3)\times10^{28}$erg s$^{-1}$Hz$^{-1}$. Comparing these limits with new light curve modeling which properly incorporates relativistic effects, we obtain limits on the energy deposited into the ejecta of $E_{ej}\lesssim(0.6-6.7)\times 10^{52}$erg ($E_{ej}\lesssim(1.8-17.6)\times10^{52}$erg) for an ejecta mass of $0.03\,M_{\odot}$ ($0.1\,M_{\odot}$). We present a uniform re-analysis of 27 SGRBs with $5.5-6.0$ GHz observations, and find that $\gtrsim50\%$ of SGRBs did not form stable magnetar remnants in their mergers. Assuming SGRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction $f_{GW190425}$ of cases, we place constraints on the maximum mass of a non-rotating neutron star (NS) ($M_{TOV}$), finding $M_{TOV}\lesssim2.23\,M_{\odot}$ for $f_{GW190425}=0.4$; this limit increases for larger values of $f_{GW190425}$. The detection (or lack thereof) of radio remnants in untargeted surveys such as the VLA Sky Survey (VLASS) could provide more stringent constraints on the fraction of mergers that produce stable remnants. If $\gtrsim30-300$ radio remnants are discovered in VLASS, this suggests that SGRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce., Comment: 21 pages, 7 figures, submitted to ApJ

Radio observations of tidal disruption events (TDEs) probe material ejected by the disruption of stars by supermassive black holes (SMBHs), uniquely tracing the formation and evolution of jets and outflows, revealing details of the disruption hydrodynamics, and illuminating the environments around previously-dormant SMBHs. To date, observations reveal a surprisingly diverse population. A small fraction of TDEs (at most a few percent) have been observed to produce radio-luminous mildly relativistic jets. The remainder of the population are radio quiet, producing less luminous jets, non-relativistic outflows or, possibly, no radio emission at all. Here, we review the radio observations that have been made of TDEs to date and discuss possible explanations for their properties, focusing on detected sources and, in particular, on the two best-studied events: Sw J1644+57 and ASASSN-14li. We also discuss what we have learned about the host galaxies of TDEs from radio observations and review constraints on the rates of bright and faint radio outflows in TDEs. Upcoming X-ray, optical, near-IR, and radio surveys will greatly expand the sample of TDEs, and technological advances open the exciting possibility of discovering a sample of TDEs in the radio band unbiased by host galaxy extinction., Comment: Resubmitted for publication in Springer Space Science Reviews following referee comments. Chapter in ISSI review "The Tidal Disruption of Stars by Massive Black Holes" vol. 79. Table 2 is available in machine-readable format upon request

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

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 <0.49$c$ by t = 98 d (3$\sigma$ limit). We also constrain the proper motion of AT 2018cow to be <0.51$c$. Since the radio emission generally traces the fastest ejecta, our observations make the presence of a long-lived relativistic jet with a lifetime of more than one month very unlikely., Comment: Accepted for publication in MNRAS 3 figurs]es

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 ALMA 97.5 GHz total intensity and linear polarization observations of the mm-band afterglow of GRB 190114C spanning 2.2 to 5.2 hours after the burst. We detect linear polarization at the $\approx 5\,\sigma$ level, decreasing from $\Pi=(0.87\pm0.13)\%$ to $(0.60\pm0.19)\%$, and evolving in polarization position angle from $(10\pm5)^\circ$ to $(-44\pm12)^\circ$ during the course of the observations. This represents the first detection of polarized millimeter emission in a $\gamma$-ray burst. We show that the optical and X-ray observations between $0.03$ days and $\sim0.3$ days are consistent with a fast cooling forward shock expanding into a wind environment. However, the optical observations at $\lesssim0.03$ days, as well as the radio and millimeter observations arise from a separate component, which we interpret as emission from the reverse-shocked ejecta. Using the measured linear polarization, we constrain the coherence scale of tangled magnetic fields in the ejecta to an angular size of $\theta_{\rm B} \approx10^{-3}$ radian, while the rotation of the polarization angle rules out the presence of large scale, ordered axisymmetric magnetic fields, and in particular a large scale toroidal field, in the jet., Comment: Published in ApJ Letters

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network., Comment: White paper submitted to the Astro2020 Decadal Survey

We present the first observations of a Type I superluminous supernova (SLSN) at $\gtrsim 1000$ days after maximum light. We observed SN 2015bn using the Hubble Space Telescope Advanced Camera for Surveys in the F475W, F625W and F775W filters at 721 days and 1068 days. SN 2015bn is clearly detected and resolved from its compact host, allowing reliable photometry. A galaxy template constructed from these data further enables us to isolate the SLSN flux in deep ground-based imaging. We measure a light curve decline rate at $>700$ days of $0.19 \pm 0.03$ mag (100 d)$^{-1}$, much shallower than the earlier evolution, and slower than previous SLSNe (at any phase) or the decay rate of $^{56}$Co. Neither additional radioactive isotopes nor a light echo can consistently account for the slow decline. A spectrum at 1083 days shows the same [O I] and [Ca II] lines as seen at $\sim300-400$ days, with no new features to indicate strong circumstellar interaction. Radio limits with the Very Large Array rule out an extended wind for mass-loss rates $10^{-2.7} \lesssim \dot{M}/v_{10} \lesssim 10^{-1.1}$ M$_\odot$ yr$^{-1}$ (where $v_{10}$ is the wind velocity in units of 10 km s$^{-1}$). The optical light curve is consistent with $L \propto t^{-4}$, which we show is expected for magnetar spin-down with inefficient trapping; furthermore, the evolution matches predictions from earlier magnetar model fits. The opacity to magnetar radiation is constrained at $\sim 0.01$ cm$^2$ g$^{-1}$, consistent with photon-matter pair-production over a broad $\sim$GeV-TeV range. This suggests the magnetar spectral energy distribution, and hence the 'missing energy' leaking from the ejecta, may peak in this range., Comment: Accepted for publication in ApJL, updated to match accepted version

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

We present Karl G. Jansky Very Large Array radio observations of the long gamma-ray burst GRB 161219B ($z=0.147$) spanning $1-37$ GHz. The data exhibit unusual behavior, including sharp spectral peaks and minutes-timescale large-amplitude variability centered at $20$ GHz and spanning the full frequency range. We attribute this behavior to scattering of the radio emission by the turbulent ionized Galactic interstellar medium (ISM), including both diffractive and refractive scintillation. However, the scintillation is much stronger than predicted by a model of the Galactic electron density distribution (NE2001); from the measured variability timescale and decorrelation bandwidth we infer a scattering measure of $SM\approx {(8-70)\times 10^{-4}}$ kpc m$^{-20/3}$ (up to ${25}$ times larger than predicted in NE2001) and a scattering screen distance of $d_{\rm scr}\approx {0.2-3}$ kpc. We infer an emission region size of $\theta_s \approx {0.9-4}$ $\mu$as ($\approx {(1-4)}\times 10^{16}$ cm) at $\approx4$ days, and find that prior to 8 days the source size is an order of magnitude smaller than model predictions for a uniformly illuminated disk or limb-brightened ring, indicating a slightly off-axis viewing angle or significant substructure in the emission region. Simultaneous multi-hour broadband radio observations of future GRB afterglows will allow us to characterize the scintillation more completely, and hence to probe the observer viewing angle, the evolution of the jet Lorentz factor, the structure of the afterglow emission regions, and ISM turbulence at high Galactic latitudes., Comment: 14 pages, 7 figures, 1 table. Accepted to ApJ. The data used to create Figures 2-6 is available with the published article, or by request to kate.alexander@northwestern.edu

We present the first effort to aggregate, homogenize, and uniformly model the combined ultraviolet, optical, and near-infrared dataset for the electromagnetic counterpart of the binary neutron star merger GW170817. By assembling all of the available data from 18 different papers and 46 different instruments, we are able to identify and mitigate systematic offsets between individual datasets, and to identify clear outlying measurements, with the resulting pruned and adjusted dataset offering an opportunity to expand the study of the kilonova. The unified dataset includes 647 individual flux measurements, spanning 0.45 to 29.4 days post-merger, and thus has greater constraining power for physical models than any single dataset. We test a number of semi-analytical models and find that the data are well modeled with a three-component kilonova model: a "blue" lanthanide-poor component with Mej~0.020 Msol and vej~0.27c; an intermediate opacity "purple" component with Mej~0.047 Msol and vej~0.15c; and a "red" lanthanide-rich component with Mej~0.011 Msol and vej~0.14c. We further explore the possibility of ejecta asymmetry and its impact on the estimated parameters. From the inferred parameters we draw conclusions about the physical mechanisms responsible for the various ejecta components, the properties of the neutron stars, and, combined with an up-to-date merger rate, the implications for r-process enrichment via this channel. To facilitate future studies of this keystone event we make the unified dataset and our modeling code public., Comment: 35 pages; 5 figures; 3 tables. Significant changes made to Table 3. Accepted to ApJL

We present detailed multi-frequency, multi-epoch radio observations of GRB 140304A at $z=5.283$ from 1 to 86 GHz and 0.45 d to 89 d. The radio and mm data exhibit unusual multiple spectral components, which cannot be simply explained by standard forward and reverse shock scenarios. Through detailed multi-wavelength analysis spanning radio to X-rays, we constrain the forward shock parameters to $E_{\rm K, iso}\approx4.9\times10^{54}\,$erg, $A_* \approx 2.6\times10^{-2}$, $\epsilon_{\rm e}\approx2.5\times10^{-2}$, $\epsilon_{\rm B}\approx5.9\times10^{-2}$, $p\approx2.6$, and $\theta_{\rm jet}\approx1.1^{\circ}$, yielding a beaming corrected $\gamma$-ray and kinetic energy, $E_{\gamma}\approx2.3\times10^{49}\,$erg and $E_{\rm K}\approx9.5\times10^{50}\,$erg, respectively. We model the excess radio emission as due to a combination of a late-time reverse shock (RS) launched by a shell collision, which also produces a re-brightening in the X-rays at $\approx0.26\,$d, and either a standard RS or diffractive interstellar scintillation. Under the standard RS interpretation, we invoke consistency arguments between the forward and reverse shocks to derive a deceleration time, $t_{\rm dec}\approx100\,$s, the ejecta Lorentz factor, $\Gamma(t_{\rm dec})\approx300$, and a low RS magnetization, $R_{\rm B}\approx0.6$. Our observations highlight both the power of radio observations in capturing RS emission and thus constraining the properties of GRB ejecta and central engines, and the challenge presented by interstellar scintillation in conclusively identifying RS emission in GRB radio afterglows., Comment: 24 pages, 17 figures, 7 tables; submitted to ApJ

While the Atacama Large Millimeter/submillimeter Array (ALMA) is a uniquely powerful telescope, its impact in certain fields of astrophysics has been limited by observatory policies rather than the telescope's innate technical capabilities. In particular, several observatory policies present challenges for observations of variable, mobile, and/or transient sources --- collectively referred to here as "time-domain" observations. In this whitepaper we identify some of these policies, describe the scientific applications they impair, and suggest changes that would increase ALMA's science impact in Cycle 6 and beyond. Parties interested in time-domain science with ALMA are encouraged to join the ALMA Time-domain Special Interest Group (ATSIG) by signing up for the ATSIG mailing list at https://groups.google.com/group/alma-td-sig ., Comment: 9 pages; whitepaper submitted to the ALMA Science Advisory Council; corresponding author P. K. G. Williams (pwilliams@cfa.harvard.edu)

We present the second multi-frequency radio detection of a reverse shock in a $\gamma$-ray burst. By combining our extensive radio observations of the Fermi-LAT GRB 160509A at $z = 1.17$ up to $20$ days after the burst with Swift X-ray observations and ground-based optical and near-infrared data, we show that the afterglow emission comprises distinct reverse shock and forward shock contributions: the reverse shock emission dominates in the radio band at $\lesssim10~$days, while the forward shock emission dominates in the X-ray, optical, and near-infrared bands. Through multi-wavelength modeling, we determine a circumburst density of $n_0\approx10^{-3}~$cm$^{-3}$, supporting our previous suggestion that a low-density circumburst environment is conducive to the production of long-lasting reverse shock radiation in the radio band. We infer the presence of a large excess X-ray absorption column, $N_{\rm H} \approx 1.5\times10^{22}~$cm$^{-2}$, and a high rest-frame optical extinction, $A_{\rm V}\approx3.4~$mag. We identify a jet break in the X-ray light curve at $t_{\rm jet}\approx6~$d, and thus derive a jet opening angle of $\theta_{\rm jet}\approx4~$deg, yielding a beaming-corrected kinetic energy and radiated $\gamma$-ray energy of $E_{\rm K}\approx4\times10^{50}~$erg and $E_{\gamma}\approx1.3\times10^{51}~$erg ($1$-$10^4~$keV, rest frame), respectively. Consistency arguments connecting the forward and reverse shocks suggest a deceleration time of $t_{\rm dec} \approx 460~$s$~\approx T_{90}$, a Lorentz factor of $\Gamma(t_{\rm dec})\approx330$, and a reverse shock to forward shock fractional magnetic energy density ratio of $R_{\rm B}\equiv\epsilon_{\rm B,RS}/\epsilon_{\rm B,FS}\approx8$., Comment: 8 pages, 2 tables, 4 figures. Submitted to ApJ

We report the discovery of transient radio emission from the nearby optically-discovered TDE ASASSN-14li (distance of 90 Mpc), making it the first typical TDE detected in the radio, and unambiguously pointing to the formation of a non-relativistic outflow with a kinetic energy of $\approx (4-10)\times10^{47}$ erg, a velocity of $\approx 12,000-36,000$ km s$^{-1}$, and a mass of $\approx 3\times10^{-5}-7\times10^{-4}$ M$_{\odot}$. We show that the outflow was ejected on 2014 August 11--25, in agreement with an independent estimate of the timing of super-Eddington accretion based on the optical, UV, and X-ray observations, and that the ejected mass corresponds to about $1-10\%$ of the mass accreted in the super-Eddington phase. The temporal evolution of the radio emission also uncovers the circumnuclear density profile, $\rho(R)\propto R^{-2.5}$ on a scale of about 0.01 pc, a scale that cannot be probed via direct measurements even in the nearest SMBHs. Our discovery of radio emission from the nearest well-studied TDE to date, with a radio luminosity lower than all previous limits, indicates that non-relativistic outflows are ubiquitous in TDEs, and that future, more sensitive, radio surveys will uncover similar events., Comment: 10 pages, 6 figures, 2 tables; published in ApJ Letters. Fixes an error in our density calculation and expands discussion, major conclusions unchanged from previous version