We present the discovery of the first millimeter afterglow of a short-duration γ-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 imaging at 0.7 ≲ z ≲ 1.4. From the lack of a detectable optical afterglow, coupled with the bright millimeter counterpart, we infer a high extinction, A _V ≳ 2.6 mag along the line of sight, making this 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}_{\mathrm{jet}}={29.2}_{-4.0}^{+4.5}$ days. For a presumed redshift of z = 1, we infer an opening angle, θ _jet = (15.°5 ± 1.°4), and beaming-corrected kinetic energy of \mathrm{log}({E}_{{\rm{K}}}/\mathrm{erg})=51.8\pm 0.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 ≳2 days at rates of ≍0.8-1.6 per year, outpacing the historical discovery rate of SGRB centimeter-band afterglows.

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

We present the discovery of the radio afterglow and near-infrared (NIR) counterpart of the Swift short gamma-ray burst (GRB) GRB 200522A, located at a small projected offset of ≍1 kpc from the center of a young, star-forming host galaxy at z = 0.5536. The radio and X-ray luminosities of the afterglow are consistent with those of on-axis cosmological short GRBs. The NIR counterpart, revealed by our Hubble Space Telescope observations at a rest-frame time of ≍2.3 days, has a luminosity of ≍(1.3-1.7) × 10⁴² erg s^-1. This is substantially lower than on-axis short GRB afterglow detections but is a factor of ≍8-17 more luminous than the kilonova of GW170817 and significantly more luminous than any kilonova candidate for which comparable observations exist. The combination of the counterpart's color (i - y = -0.08 ± 0.21; rest frame) and luminosity cannot be explained by standard radioactive heating alone. We present two scenarios to interpret the broadband behavior of GRB 200522A: a synchrotron forward shock with a luminous kilonova (potentially boosted by magnetar energy deposition), or forward and reverse shocks from a ≍14°, relativistic (Γ₀ ≳ 10) jet. Models that include a combination of enhanced radioactive heating rates, low-lanthanide mass fractions, or additional sources of heating from late-time central engine activity may provide viable alternate explanations. If a stable magnetar was indeed produced in GRB 200522A, we predict that late-time radio emission will be detectable starting ≍0.3-6 yr after the burst for a deposited energy of ≍10⁵³ erg. Counterparts of similar luminosity to GRB 200522A associated with gravitational wave events will be detectable with current optical searches to ≍250 Mpc.

The Fermi-LAT collaboration presented the second gamma-ray burst (GRB) catalog covering its first 10 years of operations. A significant fraction of afterglow-phase light curves in this catalog cannot be explained by the closure relations of the standard synchrotron forward-shock model, suggesting that there could be an important contribution from another process. In view of the above, we derive the synchrotron self-Compton (SSC) light curves from the reverse shock in the thick- and thin-shell regime for a uniform-density medium. We show that this emission could explain the GeV flares exhibited in some LAT light curves. Additionally, we demonstrate that the passage of the forward shock synchrotron cooling break through the LAT band from jets expanding in a uniform-density environment may be responsible for the late time (≍10² s) steepening of LAT GRB afterglow light curves. As a particular case, we model the LAT light curve of GRB 160509A that exhibited a GeV flare together with a break in the long-lasting emission, and also two very high energy photons with energies of 51.9 and 41.5 GeV observed 76.5 and 242 s after the onset of the burst, respectively. Constraining the microphysical parameters and the circumburst density from the afterglow observations, we show that the GeV flare is consistent with an SSC reverse-shock model, the break in the long-lasting emission with the passage of the synchrotron cooling break through the Fermi-LAT band, and the very energetic photons with SSC emission from the forward shock, when the outflow carries a significant magnetic field (R_B ≃ 30) and it decelerates in a uniform-density medium with a very low density ( 

Massive, rapidly spinning magnetar remnants produced as a result of binary neutron-star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a synchrotron radio signal from the interaction of the ejecta with the circumburst medium. Here, we present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (GRBs; z < 0.5) on rest-frame timescales of ≍2.4-13.9 yr following the bursts. We place 3σ limits on radio continuum emission of Fν ≲ 6-20 μJy at the burst positions, or Lν ≲ (0.6-8.3) × 1028 erg s-1 Hz-1. Comparing these limits with new light-curve modeling that properly incorporates relativistic effects, we obtain limits on the energy deposited into the ejecta of Eej ≲ (0.6-6.7) × 1052 erg (  Eej≲(1.8--17.6)×1052Eej≲(1.8--17.6)×1052  erg) for an ejecta mass of 0.03 M☉ (0.1 M☉). We present a uniform reanalysis of 27 short GRBs with 5.5-6.0 GHz observations, and find that ≳50% of short GRBs did not form stable magnetar remnants in their mergers. Assuming short GRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction fGW190425 of cases, we place constraints on the maximum mass of a nonrotating neutron star (NS; Tolman-Oppenheimer-Volkoff mass; MTOV), finding  MTOV≲2.23M⊙MTOV≲2.23M⊙  for fGW190425 = 0.4; this limit increases for larger values of fGW190425. 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  ≳30--300≳30--300  radio remnants are discovered in VLASS, this suggests that short GRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce.

Context. The identification of the sources that reionized the Universe and their specific contribution to this process are key missing pieces of our knowledge of the early Universe. Faint star-forming galaxies may be the main contributors to the ionizing photon budget during the epoch of reionization, but their escaping photons cannot be detected directly due to inter-galactic medium opacity. Hence, it is essential to characterize the properties of faint galaxies with significant Lyman continuum (LyC) photon leakage up to z ∼ 4 to define indirect indicators allowing analogs to be found at the highest redshift.
Aims: Long gamma-ray bursts (LGRBs) typically explode in star-forming regions of faint, star-forming galaxies. Through LGRB afterglow spectroscopy it is possible to detect directly LyC photons. Our aim is to use LGRBs as tools to study LyC leakage from faint, star-forming galaxies at high redshift.
Methods: Here we present the observations of LyC emission in the afterglow spectra of GRB 191004B at z = 3.5055, together with those of the other two previously known LyC-leaking LGRB host galaxies (GRB 050908 at z = 3.3467, and GRB 060607A at z = 3.0749), to determine their LyC escape fraction and compare their properties.
Results: From the afterglow spectrum of GRB 191004B we determine a neutral hydrogen column density at the LGRB redshift of log(NH I/cm-2) = 17.2 ± 0.15, and negligible extinction (AV = 0.03 ± 0.02 mag). The only metal absorption lines detected are C IV and Si IV. In contrast to GRB 050908 and GRB 060607A, the host galaxy of GRB 191004B displays significant Lyman-alpha (Lyα) emission. From its Lyα emission and the non-detection of Balmer emission lines we constrain its star-formation rate (SFR) to 1 ≤ SFR ≤ 4.7 M☉ yr-1. We fit the Lyα emission with a shell model and find parameters values consistent with the observed ones. The absolute (relative) LyC escape fractions we find for GRB 191004B, GRB 050908 and GRB 060607A are of 0.35-0.11+0.10 (0.43-0.13+0.12), 0.08-0.04+0.05 (0.08-0.04+0.05) and 0.20-0.05+0.05 (0.45-0.15+0.15), respectively. We compare the LyC escape fraction of LGRBs to the values of other LyC emitters found from the literature, showing that LGRB afterglows can be powerful tools to study LyC escape for faint high-redshift star-forming galaxies. Indeed we could push LyC leakage studies to much higher absolute magnitudes. The host galaxies of the three LGRBs presented here have all M1600 > -19.5 mag, with the GRB 060607A host at M1600 > -16 mag. LGRB hosts may therefore be particularly suitable for exploring the ionizing escape fraction in galaxies that are too faint or distant for conventional techniques. Furthermore, the time involved is minimal compared to galaxy studies.

Results based on observations carried out at ESO Observatory, Paranal, Chile, by the Stargate consortium under Program ID: 0104.D-0600, P.I.: N. Tanvir.

We present the discovery of the optical afterglow and host galaxy of the Swift short-duration gamma-ray burst (SGRB) GRB 181123B. Observations with Gemini-North starting ≍9.1 hr after the burst reveal a faint optical afterglow with i ≍ 25.1 mag at an angular offset of 0"59 ± 0"16 from its host galaxy. Using grizYJHK observations, we measure a photometric redshift of the host galaxy of z=1.77+0.30−0.17z=1.77−0.17+0.30 . From a combination of Gemini and Keck spectroscopy of the host galaxy spanning 4500-18000 Å, we detect a single emission line at 13390 Å, inferred as Hβ at z = 1.754 ± 0.001 and corroborating the photometric redshift. The host galaxy properties of GRB 181123B are typical of those of other SGRB hosts, with an inferred stellar mass of ≍9.1 × 10⁹ M_☉, a mass-weighted age of ≍0.9 Gyr, and an optical luminosity of ≍0.9L*. At z = 1.754, GRB 181123B is the most distant secure SGRB with an optical afterglow detection and one of only three at z > 1.5. Motivated by a growing number of high-z SGRBs, we explore the effects of a missing z > 1.5 SGRB population among the current Swift sample on delay time distribution (DTD) models. We find that lognormal models with mean delay times of ≍4-6 Gyr are consistent with the observed distribution but can be ruled out to 95% confidence, with an additional ≍one to five Swift SGRBs recovered at z > 1.5. In contrast, power-law models with ∝t^-1 are consistent with the redshift distribution and can accommodate up to ≍30 SGRBs at these redshifts. Under this model, we predict that ≍1/3 of the current Swift population of SGRBs is at z > 1. The future discovery or recovery of existing high-z SGRBs will provide significant discriminating power on their DTDs and thus their formation channels.

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 × 1050 erg. An excess bump appears in the UV light curve about 50 d after bolometric peak, followed by a flattening beyond 250 d. We detect a constant X-ray source present for at least 86 d. The X-ray spectrum shows a total unabsorbed flux of ∼4 × 10-14 erg cm-2 s-1 and is best fit by a blackbody plus power-law model with a photon index of Γ = 0.8. A thermal X-ray model is unable to account for photons >1 keV, while a radio non-detection favours 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 × 106 M☉ disrupting a 0.1 M☉ star; the model suggests the star was likely only partially disrupted, based on the derived impact parameter of β = 0.6. The low optical depth implied by the small debris mass may explain how we are able to see hydrogen emission with disc-like line profiles in the spectra of AT 2018hyz (see our companion paper).

We present X-ray and radio observations of the Fast Blue Optical Transient CRTS-CSS161010 J045834−081803 (CSS161010 hereafter) at t = 69–531 days. CSS161010 shows luminous X-ray (L _x ~ 5 × 10³⁹ erg s⁻¹) and radio (L _ν  ~ 10²⁹ erg s⁻¹ Hz⁻¹) emission. The radio emission peaked at ~100 days post-transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blast wave. CSS161010 launched a mildly relativistic outflow with velocity Γβc ≥ 0.55c at ~100 days. This is faster than the non-relativistic AT 2018cow (Γβc ~ 0.1c) and closer to ZTF18abvkwla (Γβc ≥ 0.3c at 63 days). The inferred initial kinetic energy of CSS161010 (E _k   10⁵¹ erg) is comparable to that of long gamma-ray bursts, but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger (). This is consistent with the lack of observed γ-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass M _* ~ 10⁷M _⊙ and specific star formation rate sSFR ~ 0.3 Gyr⁻¹. This mass is among the lowest inferred for host galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally located intermediate-mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate < 0.4% of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

We present Atacama Large Millimeter/submillimeter Array 97.5 GHz total intensity and linear polarization observations of the mm-band afterglow of GRB 190114C spanning 2.2-5.2 hr after the burst. We detect linear polarization at the ≈5σ level, decreasing from Π = (0.87 ± 0.13)% to (0.60 ± 0.19)%, and evolving in polarization position angle from (10 ± 5)° to (-44 ± 12)° during the course of the observations. This represents the first detection and measurement of the temporal evolution of polarized radio/millimeter emission in a γ-ray burst. We show that the optical and X-ray observations between 0.03 days and ̃0.3 days are consistent with a fast-cooling forward shock expanding into a wind environment. However, the optical observations at ≲0.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 θ _B ≈ 10^-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.

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

We present a revised and complete optical afterglow light curve of the binary neutron star merger GW170817, enabled by deep Hubble Space Telescope (HST) F606W observations at ≈584 days post-merger, which provide a robust optical template. The light curve spans ≈110-362 days, and is fully consistent with emission from a relativistic structured jet viewed off-axis, as previously indicated by radio and X-ray data. Combined with contemporaneous radio and X-ray observations, we find no spectral evolution, with a weighted average spectral index of < β > =-0.583+/- 0.013, demonstrating that no synchrotron break frequencies evolve between the radio and X-ray bands over these timescales. We find that an extrapolation of the post-peak temporal slope of GW170817 to the luminosities of cosmological short gamma-ray bursts matches their observed jet break times, suggesting that their explosion properties are similar, and that the primary difference in GW170817 is viewing angle. Additionally, we place a deep limit on the luminosity and mass of an underlying globular cluster (GC) of L ≲ 6.7 × 10³ L _☉, or M ≲ 1.3 × 10⁴ M _☉, at least 4 standard deviations below the peak of the GC mass function of the host galaxy, NGC 4993. This limit provides a direct and strong constraint that GW170817 did not form and merge in a GC. As highlighted here, HST (and soon the James Webb Space Telescope) enables critical observations of the optical emission from neutron star merger jets and outflows.

We present comprehensive multiwavelength radio to X-ray observations of GRB 181201A spanning from ≈150 s to ≈163 days after the burst, comprising the first joint ALMA-VLA-GMRT observations of a gamma-ray burst (GRB) afterglow. The radio and millimeter-band data reveal a distinct signature at ≈3.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 discovered degeneracies in the model parameters. The joint RS-FS modeling reveals a weakly magnetized (σ ≈ 3 × 10^-3), mildly relativistic RS, from which we derive an initial bulk Lorentz factor of Γ₀ ≈ 103 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, nonrelativistic RSs. We present and begin to address new challenges in modeling posed by the present generation of comprehensive, multifrequency data sets.

In the era of multi-messenger astronomy the exploration of the early emission from transients is key for understanding the encoded physics. At the same time, current generation networks of fully-robotic telescopes provide new opportunities in terms of fast followup and sky coverage. This work describes our pipeline designed for robotic optical followup of gamma-ray bursts with the Las Cumbrés Observatory network. We designed a Python code to promptly submit observation requests to the Las Cumbrés Observatory network within 3 minutes of the receipt of the socket notice. Via Telegram the pipeline keeps the users informed, allowing them to take control upon request. Our group was able to track the early phases of the evolution of the optical output from gamma-ray bursts with a fully-robotic procedure and here we report the case of GRB180720B as an example. The developed pipeline represents a key ingredient for any reliable and rapid (minutes timescale) robotic telescope system. While successfully utilized and adapted for LCO, it can also be adapted to any other robotic facilities.

We present Chandra and Very Large Array observations of GW170817 at ∼521─743 days post-merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time nonthermal emission follows the expected evolution of an off-axis relativistic jet, with a steep temporal decay {F}ν \propto {t}-1.95+/- 0.15 and power-law spectrum {F}ν \propto {ν }-0.575+/- 0.007. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and find n≤slant 9.6× {10}-3 {cm}}-3. This measurement is independent from inferences based on jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy {E}0={1.5}-1.1+3.6× {10}49 {erg} ({E}iso}={2.1}-1.5+6.4× {10}52 {erg}) and jet opening angle {θ }0={5.9}-0.7+1.0 \deg with characteristic Lorentz factor {{{Γ }}}j={163}-43+23, expanding in a low-density medium with {n}0={2.5}-1.9+4.1× {10}-3 {cm}}-3 and viewed {θ }obs}={30.4}-3.4+4.0 \deg off-axis. The synchrotron emission originates from a power-law distribution of electrons with index p={2.15}-0.02+0.01. The shock microphysics parameters are constrained to {ε }{{e}}={0.18}-0.13+0.30 and {ε }{{B}}={2.3}-2.2+16.0× {10}-3. Furthermore, we investigate the presence of X-ray flares and find no statistically significant evidence of ≥2.5σ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy {E}kKN}\propto {({{Γ }}β )}-α into the environment, finding that shallow stratification indexes α ≤ 6 are disfavored. Future radio and X-ray observations will refine our inferences on the fastest kilonova ejecta properties.

We present the first extensive radio to γ-ray observations of a fast-rising blue optical transient, AT 2018cow, over its first ∼100 days. AT 2018cow rose over a few days to a peak luminosity L _pk ∼ 4 × 10⁴⁴ erg s^‑1, exceeding that of superluminous supernovae (SNe), before declining as L ∝ t ^‑2. Initial spectra at δt ≲ 15 days were mostly featureless and indicated large expansion velocities v ∼ 0.1c and temperatures reaching T ∼ 3 × 10⁴ K. Later spectra revealed a persistent optically thick photosphere and the emergence of H and He emission features with v ∼ 4000 km s^‑1 with no evidence for ejecta cooling. Our broadband monitoring revealed a hard X-ray spectral component at E ≥ 10 keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT 2018cow showed bright radio emission consistent with the interaction of a blast wave with v _sh ∼ 0.1c with a dense environment (\dot{M}∼ {10}^-3-{10}^-4 {M}_ȯ {yr}}^-1 for v _w = 1000 km s^‑1). While these properties exclude ⁵⁶Ni-powered transients, our multiwavelength analysis instead indicates that AT 2018cow harbored a “central engine,” either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released ∼10⁵⁰–10^51.5 erg over ∼10³–10⁵ s and resides within low-mass fast-moving material with equatorial–polar density asymmetry (M _ej,fast ≲ 0.3 M _⊙). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black holes are disfavored by the large environmental density probed by the radio observations.

Whether stars could have driven the reionization of the intergalactic medium depends critically on the proportion of ionizing radiation that escapes the galaxies in which it is produced. Spectroscopy of gamma-ray burst (GRB) afterglows can be used to estimate the opacity to extreme ultraviolet (EUV) radiation along the lines-of-sight to the bursts. Assuming that long-duration GRBs trace the locations of the massive stars dominating EUV production, the average escape fraction of ionizing radiation can be calculated independently of galaxy size or luminosity. Here we present a compilation of H I column density (N_{H I}) measures for 140 GRBs in the range 1.6 < z < 6.7. Although the sample is heterogeneous, in terms of spectral resolution and signal-to-noise ratio, fits to the Ly α absorption line provide robust constraints on N_{H I}, even for spectra of insufficient quality for other purposes. Thus we establish an escape fraction at the Lyman limit of <f_esc> ≈ 0.005, with a 98 per cent confidence upper limit of <f_esc> ≈ 0.015. This analysis suggests that stars provide a small contribution to the ionizing radiation budget at z < 5. At higher redshifts firm conclusions are limited by the small size of the GRB sample (7/140), but any decline in average H I column density seems to be modest. We also find no significant correlation of N_{H I} with galaxy UV luminosity or host stellar mass. We discuss in some detail potential biases and argue that, while not negligible, systematic errors in f_esc are unlikely to be more than a factor ˜2 in either direction, and so would not affect the primary conclusions. Given that many GRB hosts are low-metallicity dwarf galaxies with high specific star-formation rates, these results present a particular problem for the hypothesis that such galaxies dominated the reionization of the Universe.

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 ≈ (8–70) × 10^‑4 kpc m^‑20/3 (up to 25 times larger than predicted in NE2001) and a scattering screen distance of d _scr ≈ 0.2–3 kpc. We infer an emission region size of {θ }_s≈ 0.9{--}4 μas (≈ (1{--}4)× {10}¹⁶ cm) at ≈4 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.