We present the discovery of the optical afterglow and host galaxy of the {\it Swift} short-duration gamma-ray burst, GRB\,181123B. Observations with Gemini-North starting at $\approx 9.1$~hr after the burst reveal a faint optical afterglow with $i\approx25.1$~mag, at an angular offset of 0.59 $\pm$ 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.17}$. From a combination of Gemini and Keck spectroscopy of the host galaxy spanning 4500-18000~\AA , we detect a single emission line at 13390~\AA, inferred as H$\beta$ at $z = 1.754 \pm 0.001$ and corroborating the photometric redshift. The host galaxy properties of GRB\,181123B are typical to those of other SGRB hosts, with an inferred stellar mass of $\approx 1.7 \times 10^{10}\,M_{\odot}$, mass-weighted age of $\approx 0.9$~Gyr and optical luminosity of $\approx 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 {\it Swift} sample on delay time distribution models. We find that log-normal models with mean delay times of $\approx 4-6$~Gyr are consistent with the observed distribution, but can be ruled out to $95\%$ confidence with an additional $\approx1-5$~{\it Swift} SGRBs recovered at $z>1.5$. In contrast, power-law models with $\propto$ $t^{-1}$ are consistent with the redshift distribution and can accommodate up to $\approx30$ SGRBs at these redshifts. Under this model, we predict that $\approx 1/3$ of the current {\it 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 delay time distributions, and thus their formation channels.

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 ($\approx10^2$ 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 a 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_{\rm B} \simeq 30$) and it decelerates in a uniform-density medium with a very low density ($n=4.554^{+1.128}_{-1.121}\times 10^{-4}\,{\rm cm^{-3}}$).

We present the observations of Lyman continuum (LyC) emission in the afterglow spectra of GRB 191004B at $z=3.5055$, together with those of the other two previously known LyC-emitting long gamma-ray bursts (LGRB) (GRB 050908 at $z=3.3467$, and GRB 060607A at $z=3.0749$), to determine their LyC escape fraction and compare their properties. From the afterglow spectrum of GRB 191004B we determine a neutral hydrogen column density at the LGRB redshift of $\log(N_{\rm HI}/cm^{-2})= 17.2 \pm 0.15$, and negligible extinction ($A_{\rm V}=0.03 \pm 0.02$ mag). The only metal absorption lines detected are CIV and SiIV. In contrast to GRB 050908 and GRB 060607A, the host galaxy of GRB 191004B displays significant Ly$\alpha$ emission. From its Ly$\alpha$ emission and the non-detection of Balmer emission lines we constrain its star-formation rate (SFR) to $1 \leq$ SFR $\leq 4.7$ M$_{\odot}\ yr^{-1}$. We fit the Ly$\alpha$ emission with a shell model and find parameters values consistent with the observed ones. The absolute LyC escape fractions we find for GRB 191004B, GRB 050908 and GRB 060607A are of $0.35^{+0.10}_{-0.11}$, $0.08^{+0.05}_{-0.04}$ and $0.20^{+0.05}_{-0.05}$, 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 LGRB presented here have all $M_{\rm 1600} > -19.5$ mag, with the GRB 060607A host at $M_{\rm 1600} > -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 investment is very small compared to galaxy studies. [Abridged]

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 $\approx 2.4-13.9$ yr following the bursts. We place $3\sigma$ limits on radio continuum emission of $F_{\nu}\lesssim 6-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_{\rm ej}\lesssim (0.6-6.7)\times 10^{52}$ erg ($E_{\rm ej}\lesssim (1.7-14.1) \times 10^{52}$ erg) for an ejecta mass of $0.03\,M_{\odot}$ ($0.1\,M_{\odot}$). We present a uniform re-analysis of 27 short GRBs with 5.5-6.0 GHz observations, and find that $\gtrsim 50\%$ of short GRBs did not form stable magnetar remnants in their mergers. Assuming short GRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction $f_{\rm GW190425}$ of cases, we place constraints on the maximum mass of a non-rotating neutron star (NS) ($M_{\rm TOV}$), finding $M_{\rm TOV} \lesssim 2.4\,M_{\odot}$ for $f_{\rm GW190425}=0.4$. 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 $\gtrsim 100-1000$ radio remnants are discovered in VLASS, this suggests that short GRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce.

The detection of mJy/sub-mJy point sources is a significant challenge for single-dish radio telescopes. Detection or upper limits on the faint afterglow from GRBs or other sources at cosmological distances are important means of constraining the source modeling. Using the Sardinia Radio Telescope (SRT), we compare the sensitivity and robustness of three methods applied to the detection of faint radio sources from raster maps around a known source position: the smart quick-look method, the source extraction method (typical of high-energy astronomy), and the fit with a 2-D Gaussian. We developed a Python code specific for the analysis of point-like radio sources applied to the SRT C-band (6.9 GHz) observations of both undetected sources (GRB afterglows of 181201A and 190114C) and the detected Galactic X-ray binary GRS 1915+105. Our comparative analysis of the different detection methods made extensive use of simulations as a useful complement to actual radio observations. The best method for the SRT data analysis is the fit with a 2-D Gaussian, as it pushes down the sensitivity limits of single-dish observations -- with respect to more traditional techniques -- to ~ 1.8 mJy, improving by ~ 40 % compared with the initial value. This analysis shows that -- especially for faint sources -- good maps of the scanned region pre- or post-outburst are essential.

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) 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 blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-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_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. 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.

Polarization measurements of gamma-ray burst (GRB) afterglows are a promising means of probing the structure, geometry, and magnetic composition of relativistic GRB jets. However, a precise treatment of instrumental calibration is vital for a robust physical interpretation of polarization data, requiring tests of and validations against potential instrumental systematics. We illustrate this with ALMA Band 3 (97.5 GHz) observations of GRB 171205A taken $\approx5.19$ days after the burst, where a detection of linear polarization was recently claimed. We describe a series of tests for evaluating the stability of polarization measurements with ALMA. Using these tests to re-analyze and evaluate the archival ALMA data, we uncover systematics in the polarization calibration at the $\approx0.09\%$ level. We derive a 3$\sigma$ upper limit on the linearly polarized intensity of $P<97.2~\mu$Jy, corresponding to an upper limit on the linear fractional polarization of $\Pi_{\rm L}<0.30\%$, in contrast to the previously claimed detection. Our upper limit improves upon existing constraints on the intrinsic polarization of GRB radio afterglows by a factor of 3. We discuss this measurement in the context of constraints on the jet magnetic field geometry. We present a compilation of polarization observations of GRB radio afterglows, and demonstrate that a significant improvement in sensitivity is desirable for eventually detecting signals polarized at the $\approx0.1\%$ level from typical radio afterglows.

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

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.

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.

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

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 $\approx\!584$ days post-merger, which provide a robust optical template. The light curve spans $\approx 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 $\langle \beta \rangle = -0.583 \pm 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 GRBs 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 of $L \lesssim 6.7 \times 10^{3}\,L_{\odot}$, or $M \lesssim 1.3 \times 10^{4}\,M_{\odot}$, at least 4 standard deviations below the peak of the globular cluster mass function of the host galaxy, NGC4993. This limit provides a direct and strong constraint that GW170817 did not form and merge in a globular cluster. As highlighted here, HST (and soon JWST) enables critical observations of the optical emission from neutron star merger jets and outflows.

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

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 Cumbres Observatory network. We designed a Python code to promptly submit observation requests to the Las Cumbres 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 represent 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.

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.

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.

The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars ($\sim$10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources.

To date, nearly all multi-wavelength modeling of long-duration gamma-ray bursts has ignored synchrotron radiation from the significant population of electrons expected to pass the shock without acceleration into a power-law distribution. We investigate the effect of including the contribution of thermal, non-accelerated electrons to synchrotron absorption and emission in the standard afterglow model, and show that these thermal electrons provide an additional source of opacity to synchrotron self-absorption, and yield an additional emission component at higher energies. The extra opacity results in an increase in the synchrotron self-absorption frequency by factors of 10--100 for fiducial parameters. The nature of the additional emission depends on the details of the thermal population, but is generally observed to yield a spectral peak in the optical brighter than radiation from the nonthermal population by similar factors a few seconds after the burst, remaining detectable at millimeter and radio frequencies several days later.

We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 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. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow 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 $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). 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.

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.

The part played by stars in the ionization of the intergalactic medium remains an open question. A key issue is the proportion of the stellar ionizing radiation that escapes the galaxies in which it is produced. Spectroscopy of gamma-ray burst afterglows can be used to determine the neutral hydrogen column-density in their host galaxies and hence the opacity to extreme ultra-violet radiation along the lines-of-sight to the bursts. Thus, making the reasonable assumption that long-duration GRB locations are representative of the sites of massive stars that dominate EUV production, one can calculate an average escape fraction of ionizing radiation in a way that is independent of galaxy size, luminosity or underlying spectrum. Here we present a sample of NH measures for 138 GRBs in the range 1.6<z<6.7 and use it to establish an average escape fraction at the Lyman limit of <fesc>~0.005, with a 98% confidence upper limit of ~0.015. This analysis suggests that stars provide a small contribution to the ionizing radiation budget of the IGM at z<5, where the bulk of the bursts lie. At higher redshifts, z>5, firm conclusions are limited by the small size of the GRB sample, but any decline in average HI column-density seems to be modest. We also find no indication of a significant correlation of NH with galaxy UV luminosity or host stellar mass, for the subset of events for which these are available. We discuss in some detail a number of selection effects and potential biases. Drawing on a range of evidence we argue that such effects, while not negligible, are unlikely to produce systematic errors of more than a factor ~2, and so would not affect the primary conclusions. Given that many GRB hosts are low metallicity, high specific star-formation rate, dwarf galaxies, these results present a particular problem for the hypothesis that such galaxies dominated the reionization of the universe.

In the standard synchrotron afterglow model, a power law of electrons is responsible for all aspects of photon production and absorption. Recent numerical work has shown that the vast majority of particles in the downstream medium are actually "thermal" particles, which were shock-heated but did not enter the Fermi acceleration process (the name stands in contrast to the nonthermal high-energy tail, rather than connoting a Maxwellian distribution). There are substantial differences at optical and higher energies when these thermal electrons participate in the afterglow, but early work along these lines ignored the radio end of the electromagnetic spectrum. We report here on an extension of previous Monte Carlo simulations of gamma-ray burst afterglows. The model now includes the synchrotron self-absorption (SSA) process and so can simulate afterglows across the entire EM spectrum, and several orders of magnitude in time. In keeping with earlier work, inclusion of the thermal electrons increases the SSA frequency by a factor of 30, and the radio intensity by a factor of 100. Furthermore, these changes happen with no modification to the late optical or X-ray afterglow. Our results provide very strong evidence that thermal electrons must be considered in any multiwavelength model for afterglows.

We present the second edition of a Best Practices Guide for academic departments and other institutions striving to create more inclusive environments for physicists and astronomers in the LGBT+ community. Our recommendations incorporate new research since the original, 2014 edition, and are designed for anyone who wishes to become aware of -- and help mitigate -- the extra burdens that face members of the LGBT+ community in the physical sciences.

We present light curves and classification spectra of 17 hydrogen-poor superluminous supernovae (SLSNe) from the Pan-STARRS1 Medium Deep Survey (PS1 MDS). Our sample contains all objects from the PS1 MDS sample with spectroscopic classification that are similar to either of the prototypes SN2005ap or SN2007bi, without an explicit limit on luminosity. With a redshift range $0.3 < z < 1.6$, PS1MDS is the first SLSN sample primarily probing the high-redshift population; our multi-filter PS1 light curves probe the rest-frame UV emission, and hence the peak of the spectral energy distribution. We measure the temperature evolution and construct bolometric light curves, and find peak luminosities of $(0.5-5) \times 10^{44}$ erg s$^{-1}$ and lower limits on the total radiated energies of $(0.3-2) \times 10^{51}$ erg. The light curve shapes are diverse, with both rise- and decline times spanning a factor of $\sim 5$, and several examples of double-peaked light curves. When correcting for the flux-limited nature of our survey, we find a median peak luminosity at 4000 {\AA} of $M_{\rm 4000} = -21.1$ mag, and a spread of $\sigma = 0.7$ mag.

Jet breaks in gamma ray burst (GRB) afterglows provide a direct probe of their collimation angle. Modeling a jet break requires an understanding of the "jet spreading" process, whereby the jet transitions from a collimated outflow into the spherical Sedov-Taylor solution at late times. Currently, direct numerical calculations are the most accurate way to capture the deceleration and spreading process, as analytical models have previously given inaccurate descriptions of the dynamics. Here (in paper I) we present a new, semi-analytical model built empirically by performing relativistic numerical jet calculations and measuring the relationship between Lorentz factor and opening angle. We then calculate the Lorentz factor and jet opening angle as a function of shock radius and compare to the numerical solutions. Our analytic model provides an efficient means of computing synthetic GRB afterglow light curves and spectra, which is the focus of paper II.

We present multi-wavelength observations and modeling of the exceptionally bright long $\gamma$-ray burst GRB 160625B. The optical and X-ray data are well-fit by synchrotron emission from a collimated blastwave with an opening angle of $\theta_j\approx 3.6^\circ$ and kinetic energy of $E_K\approx 2\times10^{51}$ erg, propagating into a low density ($n\approx 5\times10^{-5}$ cm$^{-3}$) medium with a uniform profile. The forward shock is sub-dominant in the radio band; instead, the radio emission is dominated by two additional components. The first component is consistent with emission from a reverse shock, indicating an initial Lorentz factor of $\Gamma_0\gtrsim 100$ and an ejecta magnetization of $R_B\approx 1-100$. The second component exhibits peculiar spectral and temporal evolution and is most likely the result of scattering of the radio emission by the turbulent Milky Way interstellar medium (ISM). Such scattering is expected in any sufficiently compact extragalactic source and has been seen in GRBs before, but the large amplitude and long duration of the variability seen here are qualitatively more similar to extreme scattering events previously observed in quasars, rather than normal interstellar scintillation effects. High-cadence, broadband radio observations of future GRBs are needed to fully characterize such effects, which can sensitively probe the properties of the ISM and must be taken into account before variability intrinsic to the GRB can be interpreted correctly.

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.

We present the first results from a recently concluded study of GRBs at $z\gtrsim5$ with the Karl G. Jansky Very Large Array (VLA). Spanning $1$ to $85.5$ GHz and 7 epochs from 1.5 to 82.3 d, our observations of GRB 140311A are the most detailed joint radio and millimeter observations of a GRB afterglow at $z\gtrsim5$ to date. In conjunction with optical/near-IR and X-ray data, the observations can be understood in the framework of radiation from a single blast wave shock with energy $E_{\rm K,iso}\approx8.5\times10^{53}$ erg expanding into a constant density environment with density, $n_0\approx8\,{\rm cm}^{-3}$. The X-ray and radio observations require a jet break at $t_{\rm jet}\approx0.6$ d, yielding an opening angle of $\theta_{\rm jet}\approx4^{\circ}$ and a beaming-corrected blast wave kinetic energy of $E_{\rm K}\approx2.2\times10^{50}$ erg. The results from our radio follow-up and multi-wavelength modeling lend credence to the hypothesis that detected high-redshift GRBs may be more tightly beamed than events at lower redshift. We do not find compelling evidence for reverse shock emission, which may be related to fast cooling driven by the moderately high circumburst density.

Gamma-ray bursts (GRBs) are powerful probes of early stars and galaxies, during and potentially even before the era of reionization. Although the number of GRBs identified at z>6 remains small, they provide a unique window on typical star-forming galaxies at that time, and thus are complementary to deep field observations. We report the identification of the optical drop-out afterglow of Swift GRB 120923A in near-infrared Gemini-North imaging, and derive a redshift of z=7.84_{-0.12}^{+0.06} from VLT/X-shooter spectroscopy. At this redshift the peak 15-150 keV luminosity of the burst was 3.2x10^52 erg/s, and in fact the burst was close to the Swift/BAT detection threshold. The X-ray and near-infrared afterglow were also faint, and in this sense it was a rather typical long-duration GRB in terms of rest-frame luminosity. We present ground- and space-based follow-up observations spanning from X-ray to radio, and find that a standard external shock model with a constant-density circumburst environment with density, n~4x10^-2 cm^-3 gives a good fit to the data. The near-infrared light curve exhibits a sharp break at t~3.4 days in the observer frame, which if interpreted as being due to a jet corresponds to an opening angle of ~5 degrees. The beaming corrected gamma-ray energy is then E_gamma~2x10^50 erg, while the beaming-corrected kinetic energy is lower, E_K~10^49 erg, suggesting that GRB 120923A was a comparatively low kinetic energy event. We discuss the implications of this event for our understanding of the high-redshift population of GRBs and their identification.

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 .

We present optical spectroscopy, ultraviolet to infrared imaging and X-ray observations of the intermediate luminosity optical transient (ILOT) SN 2010da in NGC 300 (d=1.86 Mpc) spanning from -6 to +6 years relative to the time of outburst in 2010. Based on the light curve and multi-epoch SEDs of SN 2010da, we conclude that the progenitor of SN 2010da is a ~10-12 Msol yellow supergiant possibly transitioning into a blue loop phase. During outburst, SN 2010da had a peak absolute magnitude of M<-10.4 mag, dimmer than other ILOTs and supernova impostors. We detect multi-component hydrogen Balmer, Paschen, and Ca II emission lines in our high-resolution spectra, which indicate a dusty and complex circumstellar environment. Since the 2010 eruption, the star has brightened by a factor of ~5 and remains highly variable in the optical. Furthermore, we detect SN 2010da in archival Swift and Chandra observations as an ultraluminous X-ray source (L~6x10^{39} erg/s). We additionally attribute He II 4686 Angstrom and coronal Fe emission lines in addition to a steady X-ray luminosity of ~10^{37} erg/s to the presence of a compact companion.

We present the discovery of the X-ray and optical afterglows of the short-duration GRB 150101B, pinpointing the event to an early-type host galaxy at z=0.1343 +/- 0.0030. This makes GRB 150101B the most nearby short GRB with an early-type host galaxy discovered to date. Fitting the spectral energy distribution of the host galaxy results in an inferred stellar mass of ~7x10^10 M_sol, stellar population age of ~2-2.5 Gyr, and star formation rate of <0.4 M_sol yr^-1. The host of GRB 150101B is one of the largest and most luminous short GRB host galaxies, with a B-band luminosity of ~4.3L* and half-light radius of ~8 kpc. GRB 150101B is located at a projected distance of 7.35 +/- 0.07 kpc from its host center, and lies on a faint region of its host rest-frame optical light. Its location, combined with the lack of associated supernova, is consistent with a NS-NS/NS-BH merger progenitor. From modeling the evolution of the broad-band afterglow, we calculate isotropic-equivalent gamma-ray and kinetic energies of ~1.3x10^49 erg and ~(6-14)x10^51 erg, respectively, a circumburst density of ~(0.8-4)x10^-5 cm^-3, and a jet opening angle of >9 deg. Using observations extending to ~30 days, we place upper limits of <(2-4)x10^41 erg s^-1 on associated kilonova emission. We compare searches following previous short GRBs to existing kilonova models, and demonstrate the difficulty of performing effective kilonova searches from cosmological short GRBs using current ground-based facilities. We show that at the Advanced LIGO/VIRGO horizon distance of 200 Mpc, searches reaching depths of ~23-24 AB mag are necessary to probe a meaningful range of kilonova models.

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

We introduce the Swift Gamma-Ray Burst Host Galaxy Legacy Survey ("SHOALS"), a multi-observatory high-redshift galaxy survey targeting the largest unbiased sample of long-duration gamma-ray burst hosts yet assembled (119 in total). We describe the motivations of the survey and the development of our selection criteria, including an assessment of the impact of various observability metrics on the success rate of afterglow-based redshift measurement. We briefly outline our host-galaxy observational program, consisting of deep Spitzer/IRAC imaging of every field supplemented by similarly-deep, multi-color optical/NIR photometry, plus spectroscopy of events without pre-existing redshifts. Our optimized selection cuts combined with host-galaxy follow-up have so far enabled redshift measurements for 110 targets (92%) and placed upper limits on all but one of the remainder. About 20% of GRBs in the sample are heavily dust-obscured, and at most 2% originate from z>5.5. Using this sample we estimate the redshift-dependent GRB rate density, showing it to peak at z~2.5 and fall by about an order of magnitude towards low (z=0) redshift, while declining more gradually towards high (z~7) redshift. This behavior is consistent with a progenitor whose formation efficiency varies modestly over cosmic history. Our survey will permit the most detailed examination to date of the connection between the GRB host population and general star-forming galaxies, directly measure evolution in the host population over cosmic time and discern its causes, and provide new constraints on the fraction of cosmic star-formation occurring in undetectable galaxies at all redshifts.

We present rest-frame NIR luminosities and stellar masses for a large and uniformly-selected population of GRB host galaxies using deep Spitzer Space Telescope imaging of 119 targets from the Swift GRB Host Galaxy Legacy Survey spanning 0.03 < z < 6.3, and determine the effects of galaxy evolution and chemical enrichment on the mass distribution of the GRB host population across cosmic history. We find strong evolution in the host luminosity distribution between z~0.5 (median absolute NIR AB magnitude ~ -18.5, corresponding to M* ~ 3x10^8 M_sun and z~1.5), but negligible variation between z~1.5 and z~5 (median magnitude ~ -21.2, corresponding to M* ~ 5x10^9 M_sun). Dust-obscured GRBs dominate the massive host population but are only rarely seen associated with low-mass hosts, indicating that massive star-forming galaxies are universally and (to some extent) homogeneously dusty at high-redshift while low-mass star-forming galaxies retain little dust in their ISM. Comparing our luminosity distributions to field surveys and measurements of the high-z mass-metallicity relation, our results have good consistency with a model in which the GRB rate per unit star-formation is constant in galaxies with gas-phase metallicity below approximately the Solar value but heavily suppressed in more metal-rich environments. This model also naturally explains the previously-reported "excess" in the GRB rate beyond z>2; metals stifle GRB production in most galaxies at z<1.5 but have only minor impact at higher redshifts. The metallicity threshold we infer is much higher than predicted by single-star models and favors a binary progenitor. Our observations also constrain the fraction of cosmic star-formation in low-mass galaxies undetectable to Spitzer to be a small minority at most redshifts (~10% at z~2, ~25% at z~3, and ~50% at z=3.5-6.0).

We present multi-wavelength observations and modeling of Gamma-ray Bursts (GRBs) that exhibit a simultaneous re-brightening in their X-ray and optical light curves, and are also detected at radio wavelengths. We show that the re-brightening episodes can be modeled by injection of energy into the blastwave and that in all cases the energy injection rate falls within the theoretical bounds expected for a distribution of energy with ejecta Lorentz factor. Our measured values of the circumburst density, jet opening angle, and beaming corrected kinetic energy are consistent with the distribution of these parameters for long-duration GRBs at both z~1 and z>6, suggesting that the jet launching mechanism and environment of these events are similar to that of GRBs that do not have bumps in their light curves. However, events exhibiting re-brightening episodes have lower radiative efficiencies than average, suggesting that a majority of the kinetic energy of the outflow is carried by slow-moving ejecta, which is further supported by steep measured distributions of the ejecta energy as a function of Lorentz factor. We do not find evidence for reverse shocks over the energy injection period, implying that the onset of energy injection is a gentle process. We further show that GRBs exhibiting simultaneous X-ray and optical re-brightenings are likely the tail of a distribution of events with varying rates of energy injection, forming the most extreme events in their class. Future X-ray observations of GRB afterglows with Swift and its successors will thus likely discover several more such events, while radio follow-up and multi-wavelength modeling of similar events will unveil the role of energy injection in GRB afterglows.

We study the late-time (t>0.5 days) X-ray afterglows of nearby (z<0.5) long Gamma-Ray Bursts (GRB) with Swift and identify a population of explosions with slowly decaying, super-soft (photon index Gamma_x>3) X-ray emission that is inconsistent with forward shock synchrotron radiation associated with the afterglow. These explosions also show larger-than-average intrinsic absorption (NH_x,i >6d21 cm-2) and prompt gamma-ray emission with extremely long duration (T_90>1000 s). Chance association of these three rare properties (i.e. large NH_x,i, super-soft Gamma_x and extreme duration) in the same class of explosions is statistically unlikely. We associate these properties with the turbulent mass-loss history of the progenitor star that enriched and shaped the circum-burst medium. We identify a natural connection between NH_x,i Gamma_x and T_90 in these sources by suggesting that the late-time super-soft X-rays originate from radiation reprocessed by material lost to the environment by the stellar progenitor before exploding, (either in the form of a dust echo or as reprocessed radiation from a long-lived GRB remnant), and that the interaction of the explosion's shock/jet with the complex medium is the source of the extremely long prompt emission. However, current observations do not allow us to exclude the possibility that super-soft X-ray emitters originate from peculiar stellar progenitors with large radii that only form in very dusty environments.

We present rest-frame far-infrared (FIR) and optical observations of the host galaxy of GRB090423 at z=8.23 from the Atacama Large Millimeter Array (ALMA) and the Spitzer Space Telescope, respectively. The host remains undetected to 3-sigma limits of Fnu(222 GHz)<33 microJy and Fnu(3.6 micron)<81 nJy. The FIR limit is about 20 times fainter than the luminosity of the local ULIRG Arp220, and comparable to the local starburst M82. Comparing to model spectral energy distributions we place a limit on the IR luminosity of L_IR(8-1000 micron)<3e10 Lsun, corresponding to a limit on the obscured star formation rate of SFR_IR<5 Msun/yr; for comparison, the limit on the unobscured star formation rate from Hubble Space Telescope rest-frame UV observations is SFR_UV<1 Msun/yr. We also place a limit on the host galaxy stellar mass of <5e7 Msun (for a stellar population age of 100 Myr and constant star formation rate). Finally, we compare our millimeter observations to those of field galaxies at z>4 (Lyman break galaxies, Ly-alpha emitters, and submillimeter galaxies), and find that our limit on the FIR luminosity is the most constraining to date, although the field galaxies have much larger rest-frame UV/optical luminosities than the host of GRB090423 by virtue of their selection techniques. We conclude that GRB host galaxies at z>4, especially those with measured interstellar medium metallicities from afterglow spectroscopy, are an attractive sample for future ALMA studies of high redshift obscured star formation.

We present the discovery and subsequent spectroscopy with Gemini-North of the optical afterglow of the Swift gamma-ray burst (GRB) 140515A. The spectrum exhibits a well-detected continuum at wavelengths longer than 8915 Angs with a steep decrement to zero flux blueward of 8910 Angs due to Ly-alpha absorption at redshift z~6.33. Some transmission through the Lyman-alpha forest is present at 5.2<z<5.733, but none is detected at higher redshift, consistent with previous measurements from quasars and GRB 130606A. We model the red damping wing of Lyman-alpha in three ways that provide equally good fits to the data: (a) a single host galaxy absorber at z=6.327 with log(N_HI)=18.62+/-0.08; (b) pure intergalactic medium (IGM) absorption from z=6.0 to z=6.328 with a constant neutral hydrogen fraction of x_HI=0.056+0.011-0.027; and (c) a hybrid model with a host absorber located within an ionized bubble of radius 10 comoving Mpc in an IGM with x_HI=0.12+/-0.05 (x_HI<0.21 at the 2-sigma level). Regardless of the model, the sharpness of the dropoff in transmission is inconsistent with a substantial neutral fraction in the IGM at this redshift. No narrow absorption lines from the host galaxy are detected, indicating a host metallicity of [Z/H]<~ -0.8. Even if we assume that all of the hydrogen absorption is due to the host galaxy, the column is unusually low for a GRB sightline, similar to two out of the other three highest-redshift bursts with measured log(N_HI). This is possible evidence that the escape fraction of ionizing photons from normal star-forming galaxies increases at z>~6.

We present optical spectroscopy and optical/near-IR photometry of 31 host galaxies of hydrogen-poor superluminous supernovae (SLSNe), including 15 events from the Pan-STARRS1 Medium Deep Survey. Our sample spans the redshift range 0.1 < z < 1.6 and is the first comprehensive host galaxy study of this specific subclass of cosmic explosions. Combining the multi-band photometry and emission-line measurements, we determine the luminosities, stellar masses, star formation rates and metallicities. We find that as a whole, the hosts of SLSNe are a low-luminosity (<M_B> ~ -17.3 mag), low stellar mass (<M_*> ~ 2 x 10^8 M_sun) population, with a high median specific star formation rate (<sSFR> ~ 2 Gyr^-1). The median metallicity of our spectroscopic sample is low, 12 + log(O/H}) ~ 8.35 ~ 0.45 Z_sun, although at least one host galaxy has solar metallicity. The host galaxies of H-poor SLSNe are statistically distinct from the hosts of GOODS core-collapse SNe (which cover a similar redshift range), but resemble the host galaxies of long-duration gamma-ray bursts (LGRBs) in terms of stellar mass, SFR, sSFR and metallicity. This result indicates that the environmental causes leading to massive stars forming either SLSNe or LGRBs are similar, and in particular that SLSNe are more effectively formed in low metallicity environments. We speculate that the key ingredient is large core angular momentum, leading to a rapidly-spinning magnetar in SLSNe and an accreting black hole in LGRBs.

We present radio, optical/NIR, and X-ray observations of the afterglow of the short-duration 130603B, and uncover a break in the radio and optical bands at 0.5 d after the burst, best explained as a jet break with an inferred jet opening angle of 4-8 deg. GRB 130603B is only the third short GRB with a radio afterglow detection to date, and the first time that a jet break is evident in the radio band. We model the temporal evolution of the spectral energy distribution to determine the burst explosion properties and find an isotropic-equivalent kinetic energy of (0.6-1.7) x 10^51 erg and a circumburst density of 5 x 10^-3-30 cm^-3. From the inferred opening angle of GRB 130603B, we calculate beaming-corrected energies of Egamma (0.5-2) x 10^49 erg and EK (0.1-1.6) x 10^49 erg. Along with previous measurements and lower limits we find a median short GRB opening angle of 10 deg. Using the all-sky observed rate of 10 Gpc^-3 yr^-1, this implies a true short GRB rate of 20 yr^-1 within 200 Mpc, the Advanced LIGO/VIRGO sensitivity range for neutron star binary mergers. Finally, we uncover evidence for significant excess emission in the X-ray afterglow of GRB 130603B at >1 d and conclude that the additional energy component could be due to fall-back accretion or spin-down energy from a magnetar formed following the merger.

We present optical, near-infrared, and radio observations of the afterglow of GRB 120521C. By modeling the multi-wavelength dataset, we derive a photometric redshift of z~6.0, which we confirm with a low signal-to-noise ratio spectrum of the afterglow. We find that a model with a constant-density environment provides a good fit to the afterglow data, with an inferred density of 0.05/cm^3. The radio observations reveal the presence of a jet break at 7 d, corresponding to a jet opening angle of ~ 3 deg. The beaming-corrected gamma-ray and kinetic energies are both ~ 3e50 erg. We quantify the uncertainties in our results using a detailed Markov Chain Monte Carlo analysis, which allows us to uncover degeneracies between the physical parameters of the explosion. To compare GRB 120521C to other high-redshift bursts in a uniform manner we re-fit all available afterglow data for the two other bursts at z>6 with radio detections (GRBs 050904 and 090423). We find a jet break at ~ 15 d for GRB 090423, in contrast to previous work. Based on these three events, we find that GRBs at z>6 appear to explode in constant-density environments, and exhibit a wide range of energies and densities that span the range inferred for lower redshift bursts. On the other hand, we find a hint for narrower jets in the z>6 bursts, potentially indicating a larger true event rate at these redshifts. Overall, our results indicate that long GRBs share a common progenitor population at least to z~8.

We present high signal-to-noise ratio Gemini and MMT spectroscopy of the optical afterglow of the gamma-ray burst (GRB) 130606A at redshift z=5.913, discovered by Swift. This is the first high-redshift GRB afterglow to have spectra of comparable quality to those of z~6 quasars. The data exhibit a smooth continuum at near-infrared wavelengths that is sharply cut off blueward of 8410 Angs due to absorption from Ly-alpha at redshift z~5.91, with some flux transmitted through the Ly-alpha forest between 7000-7800 Angs. We use column densities inferred from metal absorption lines to constrain the metallicity of the host galaxy between a lower limit of [Si/H]>-1.7 and an upper limit of [S/H]<-0.5 set by the non-detection of S II absorption. We demonstrate consistency between the dramatic evolution in the transmission fraction of Ly-alpha seen in this spectrum over the redshift range z=4.9 to 5.85 with that previously measured from observations of high-redshift quasars. There is an extended redshift interval of Delta-z=0.12 in the Ly-alpha forest at z=5.77 with no detected transmission, leading to a 3-sigma upper limit on the mean Ly-alpha transmission fraction of <0.2% (or tau_eff(Ly-alpha) > 6.4). This is comparable to the lowest-redshift Gunn-Peterson troughs found in quasar spectra. We set a 2-sigma upper limit of 0.11 on the neutral fraction of the IGM at the redshift of the GRB from the lack of a Ly-alpha red damping wing, assuming a model with a constant neutral density. Some Ly-beta and Ly-gamma transmission is detected in this redshift window, indicating that it is not completely opaque, and hence that the IGM is nonetheless mostly ionized at these redshifts. GRB 130606A thus for the first time realizes the promise of GRBs as probes of the first galaxies and cosmic reionization.

We present extensive radio and millimeter observations of the unusually bright GRB 130427A at z=0.340, spanning 0.67 to 12 days after the burst. Taken in conjunction with detailed multi-band UV, optical, NIR, and X-ray observations we find that the broad-band afterglow emission is composed of distinct reverse shock and forward shock contributions. The reverse shock emission dominates in the radio/millimeter and at <0.1 days in the UV/optical/NIR, while the forward shock emission dominates in the X-rays and at >0.1 days in the UV/optical/NIR. We further find that the optical and X-ray data require a Wind circumburst environment, pointing to a massive star progenitor. Using the combined forward and reverse shock emission we find that the parameters of the burst are an isotropic kinetic energy of E_Kiso~2e53 erg, a mass loss rate of Mdot~3e-8 Msun/yr (for a wind velocity of 1,000 km/s), and a Lorentz factor at the deceleration time of Gamma(200s)~130. Due to the low density and large isotropic energy, the absence of a jet break to ~15 days places only a weak constraint on the opening angle of theta_j>2.5 deg, and therefore a total energy of E_gamma+E_K>1.2e51 erg, similar to other GRBs. The reverse shock emission is detectable in this burst due to the low circumburst density, which leads to a slow cooling shock. We speculate that this is a required property for the detectability of reverse shocks in the radio and millimeter bands. Following on GRB 130427A as a benchmark event, observations of future GRBs with the exquisite sensitivity of VLA and ALMA, coupled with detailed modeling of the reverse and forward shock contributions will test this hypothesis.

We present observations of the afterglows and host galaxies of three short-duration gamma-ray bursts (GRBs): 100625A, 101219A and 110112A. We find that GRB 100625A occurred in a z=0.452 early-type galaxy with a stellar mass of 4.6e9 M_Sun and a stellar population age of 0.7 Gyr, and GRB 101219A originated in a star-forming galaxy at z=0.718 with a stellar mass of 1.4e9 M_Sun, a star formation rate of 16 M_Sun yr^-1, and a stellar population age of 50 Myr. We also report the discovery of the optical afterglow of GRB 110112A, which lacks a coincident host galaxy to i>26 mag and we cannot conclusively identify any field galaxy as a possible host. The bursts have inferred circumburst densities of ~1e-4-1 cm^-3, and isotropic-equivalent gamma-ray and kinetic energies of 1e50-1e51 erg. These events highlight the diversity of galaxies that host short GRBs. To quantify this diversity, we use the sample of 36 Swift short GRBs with robust associations to an environment (~1/2 of 68 short bursts detected by Swift to May 2012) and classify them as follows: late-type (50%), early-type (15%), inconclusive (20%), and host-less (lacking a coincident host galaxy to limits of >26 mag; 15%). To find likely ranges for the true late- and early-type fractions, we assign each of the host-less bursts to the late- or early-type category using probabilistic arguments, and consider the scenario that all hosts in the inconclusive category are early-type galaxies to set an upper bound on the early-type fraction. The most likely ranges for the late- and early-type fractions are ~60-80% and ~20-40%, respectively. We find no clear trend between gamma-ray duration and host type, and no change to the fractions when excluding events recently claimed as possible contaminants from the long GRB/collapsar population. Our reported demographics are consistent with a short GRB rate driven by both stellar mass and star formation.

We present the optical discovery and sub-arcsecond optical and X-ray localization of the afterglow of the short GRB 120804A, as well as optical, near-IR, and radio detections of its host galaxy. X-ray observations with Swift/XRT, Chandra, and XMM-Newton to ~19 d reveal a single power law decline. The optical afterglow is faint, and comparison to the X-ray flux indicates that GRB 120804A is "dark", with a rest-frame extinction of A_V~2.5 mag (at z~1.3). The intrinsic neutral hydrogen column density inferred from the X-ray spectrum, N_H~2x10^22 cm^-2, is commensurate with the large extinction. The host galaxy exhibits red optical/near-IR colors. Equally important, JVLA observations at 0.9-11 d reveal a constant 5.8 GHz flux density and an optically-thin spectrum, unprecedented for GRB afterglows, but suggestive instead of emission from the host galaxy. The optical/near-IR and radio fluxes are well fit with the scaled spectral energy distribution of the local ultra-luminous infrared galaxy (ULIRG) Arp 220 at z~1.3, with a resulting star formation rate of ~300 Msun/yr. The inferred extinction and small projected offset (2.2+/-1.2 kpc) are also consistent with the ULIRG scenario, as is the presence of a companion galaxy at a separation of about 11 kpc. The limits on radio afterglow emission, in conjunction with the observed X-ray and optical emission, require a circumburst density of ~10^-3 cm^-3 an isotropic-equivalent energy scale of E_gamma,iso ~ E_K,iso ~ 7x10^51 erg, and a jet opening angle of >8 deg. The expected fraction of luminous infrared galaxies in the short GRB host sample is ~0.01-0.3 (for pure stellar mass and star formation weighting, respectively). Thus, the observed fraction of 2 events in about 25 hosts (GRBs 120804A and 100206A), provides additional support to our previous conclusion that short GRBs track both stellar mass and star formation activity.

We present X-ray, optical, near-infrared, and radio observations of GRBs 110709B and 111215A, as well as optical and near-IR observations of their host galaxies. The combination of X-ray detections and deep optical/near-infrared limits establish both bursts as "dark". Sub-arcsecond positions enabled by radio detections lead to robust host galaxy associations, with optical detections that indicate z < 4 (110709B) and 1.8 < z < 2.7 (111215A). Using the radio and X-ray data for each burst we find that GRB 110709B requires A_V > 5.3 mag and GRB 111215A requires A_V > 8.5 mag (z=2), among the largest extinction values inferred for dark bursts to-date. The two bursts also exhibit large neutral hydrogen column densities (N_H > 10^22/cm^2; z=2) as inferred from their X-ray spectra, in agreement with the trend for dark GRBs. Finally, we find that for both bursts the afterglow emission is best explained by a collimated outflow with a total beaming-corrected energy of E_gamma+E_K ~ 7-9 x 10^51 erg (z=2) expanding into a wind medium with a high density (n~100-350 cm^-3 at 10^17 cm). While the energy release is typical of long GRBs, the inferred density may be indicative of larger mass loss rates for GRB progenitors in dusty (and hence metal rich) environments. This study establishes the critical role of radio observations in determining the origin and properties of dark GRBs. Observations with the JVLA and ALMA will provide a sample with sub-arcsecond positions and robust host associations that will help shed light on obscured star formation and the role of metallicity in GRB progenitors.

Long-duration gamma-ray bursts (GRBs) provide a premier tool for studying high-redshift star-forming galaxies thanks to their extreme brightness and association with massive stars. Here we use GRBs to study the galaxy stellar mass-metallicity (M*-Z) relation at z~3-5, where conventional direct metallicity measurements are extremely challenging. We use the ISM metallicities of LGRB hosts derived from afterglow absorption spectroscopy (Z~0.01-1 solar), in conjunction with host galaxy stellar masses determined from deep Spitzer 3.6 micron observations of 20 GRB hosts. We detect about 1/4 of the hosts with I-band magnitudes ~ -21.5 to -22.5 AB mag, and place a limit of M > -19 mag on the remaining hosts from a stacking analysis. Using these observations, we present the first rest-frame optical luminosity distribution of long GRB hosts at z>3 and find that it is similar to the distribution of long GRB hosts at z~1. In comparison to Lyman-break galaxies at the same redshift, GRB hosts are generally fainter, but the sample is too small to rule out an overall similar luminosity function. On the other hand, the GRB hosts appear to be more luminous than the population of Lyman-alpha emitters at z~3-4. Using a conservative range of mass-to-light ratios for simple stellar populations (with ages of 70 Myr to ~2 Gyr), we infer the host stellar masses and present mass-metallicity measurements at z~3-5 (<z> ~ 3.5). We find that the detected GRB hosts, with M*~2e10 solar masses, display a wide range of metallicities, but that the mean metallicity at this mass scale, Z~0.1 solar, is lower than measurements at z<3. Combined with stacking of the non-detected hosts with M*< 3e9 solar masses and Z<0.03 solar, we find evidence for the existence of an M*-Z relation at z~3.5 and continued evolution of this relation to systematically lower metallicities from z~2.

Long-duration gamma-ray bursts (GRBs) provide a premier tool for studying high-redshift star-forming galaxies thanks to their extreme brightness and association with massive stars. Here we use GRBs to study the galaxy mass-metallicity (M*-Z) relation at z~3-5, where conventional direct metallicity measurements are extremely challenging. We use the interstellar medium metallicities of long-duration GRB hosts derived from afterglow absorption pectroscopy (Z~0.01-1 solar), in conjunction with host galaxy stellar masses determined from deep Spitzer 3.6 micron observations of 20 GRB hosts. We detect about 1/4 of the hosts with I-band absolute magnitudes, M ~ -21.5 to -22.5 AB mag, and place a limit of M>-19 mag on the remaining hosts from a stacking analysis. Using a conservative range of mass-to-light ratios for simple stellar populations (with ages of 70 Myr to ~ 2 Gyr), we infer the host stellar masses and present the galaxy mass-metallicity measurements at z ~ 3-5 (<z> ~ 3.5). We find that the detected GRB hosts, with M* ~ 2e10 solar masses, display a wide range of metallicities, but that the mean metallicity at this mass scale, Z~ 0.1 solar, is lower than measurements at z < 3. Combined with stacking of the non-detected hosts (with M* < 4e9 solar masses and Z < 0.03 solar), we find evidence for the existence of an M*-Z relation at z ~ 3.5 and continued evolution of this relation to systematically lower metallicities from z ~ 2.

We report XMM-Newton observations of three FR II radio galaxies at redshifts between 0.85 and 1.34, which show extended diffuse X-ray emission within the radio lobes, likely due to inverse-Compton up-scattering of the cosmic microwave background. Under this assumption, through spectrum-fitting together with archival VLA radio observations, we derive an independent estimate of the magnetic field in the radio lobes of 3C 469.1 and compare it with the equipartition value. We find concordance between these two estimates as long as the turnover in the energy distribution of the particles occurs at a Lorentz factor in excess of ~ 250. We determine the total energy in relativistic particles in the radio emitting lobes of all three sources to range between 3e59 and 8e59 erg. The nuclei of these X-ray sources are heavily-absorbed powerful AGN.

We have measured the abundances of Fe, Si, Ni, Ti, and Na in 20 Pleiads with \teff values near solar and with low \vsini using high-resolution, high signal-to-noise echelle spectra. We have validated our procedures by also analyzing 10 field stars of a range of temperatures and metallicities that were observed by \citet{Vale05}. Our result for the Pleiades is [Fe/H] = $+0.03\pm0.02\pm0.05$ (statistical and systematic). The average of published measurements for the Pleiades is $+0.042\pm0.021$.

Pre-main sequence (PMS) binaries and multiples enable critical tests of stellar models if masses, metallicities, and luminosities of the component stars are known. We have analyzed high-resolution, high signal-to-noise echelle spectra of the quadruple-star system HD 98800 and using spectrum synthesis, computed fits to the composite spectrum for a full range of plausible stellar parameters for the components. We consistently find that sub-solar metallicity yields fits with lower $\chi^2$ values, with an overall best-fit of $[M/H] = -0.20\pm0.10$. This metallicity appears to be consistent with PMS evolutionary tracks for the measured masses and luminosities of the components of HD 98800 but additional constraints on the system and modelling are needed.

J2018-556, an X-shaped FR-I radio galaxy provides a unique opportunity to discuss the formation scenarios for the extended structures seen in this class of radio sources. An understanding of the origin of X-shaped source structures has important implications for radio galaxy models as well as black hole-black hole mergers. Guided by our ATCA radio continuum observations of J2018-556 we discuss evidence in support of 'wings' forming as a result of a redirection of backflows by asymmetric gas distribution associated with the host galaxy. Formation models invoking interaction between backflows with asymmetric environments seem the preferred model for X-shaped radio sources despite apparent discoveries of X-shaped FR-I radio sources.

We consider sputtering of dust grains, believed to be formed in cooling supernovae ejecta, under the influence of reverse shocks. In the regime of self-similar evolution of reverse shocks, we can follow the evolution of ejecta density and temperature analytically as a function of time in different parts of the ejecta, and calculate the sputtering rate of graphite and silicate grains embedded in the ejecta as they encounter the reverse shock. Through analytic (1D) calculations, we find that a fraction of dust mass ($ 1\hbox{--}20$% for silicates and %$\le 5$% for graphites) can be sputtered by reverse shocks, the fraction varying with the grain size distribution and the steepness of the density profile of the ejecta mass. It is expected that many more grains will get sputtered in the region between the forward and reverse shocks, so that our analytical results provide a lower limit to the destroyed fraction of dust mass.

The most detailed constraints on the accelerating expansion of the universe and details of nature of dark energy are derived from the high redshift supernova data, assuming that the errors in the measurements are Gaussian in nature. There is a possibility that there are direction dependent systematics in the data, either due to uncorrected, known physical processes or because there are tiny departures from the cosmological principle, making the universe slightly anisotropic. To investigate this possibility we introduce a statistic based on the extreme value theory and apply it to the gold data set from Riess et al. (2004). Our analysis indicates a systematic, direction dependent non-gaussianity at about one sigma level.