Search

The Legacy Survey of Space and Time (LSST) at Vera C. Rubin Observatory is planned to begin in the Fall of 2025. The LSST survey cadence has been designed via a community-driven process regulated by the Survey Cadence Optimization Committee (SCOC), which recommended up to 3% of the observing time to carry out Target of Opportunity (ToO) observations. Experts from the scientific community, Rubin Observatory personnel, and members of the SCOC were brought together to deliver a recommendation for the implementation of the ToO program during a workshop held in March 2024. Four main science cases were identified: gravitational wave multi-messenger astronomy, high energy neutrinos, Galactic supernovae, and small potentially hazardous asteroids possible impactors. Additional science cases were identified and briefly addressed in the documents, including lensed or poorly localized gamma-ray bursts and twilight discoveries. Trigger prioritization, automated response, and detailed strategies were discussed for each science case. This document represents the outcome of the Rubin ToO 2024 workshop, with additional contributions from members of the Rubin Science Collaborations. The implementation of the selection criteria and strategies presented in this document has been endorsed in the SCOC Phase 3 Recommendations document (PSTN-056). Although the ToO program is still to be finalized, this document serves as a baseline plan for ToO observations with the Rubin Observatory.

We present the radio afterglow of short gamma-ray burst (GRB) 230217A, which was detected less than 1 day after the gamma-ray prompt emission with the Australia Telescope Compact Array (ATCA) and the Karl G. Jansky Very Large Array (VLA). The ATCA rapid-response system automatically triggered an observation of GRB 230217A following its detection by the Neil Gehrels Swift Observatory and began observing the event just 32 minutes post-burst at 5.5 and 9 GHz for 7 hours. Dividing the 7-hour observation into three time-binned images allowed us to obtain radio detections with logarithmic central times of 1, 2.8 and 5.2 hours post-burst, the first of which represents the earliest radio detection of any GRB to date. The decline of the light curve is consistent with reverse shock emission if the observing bands are below the spectral peak and not affected by synchrotron self-absorption. This makes GRB 230217A the fifth short GRB with radio detections attributed to a reverse shock at early times ($<1$ day post-burst). Following brightness temperature arguments, we have used our early radio detections to place the highest minimum Lorentz factor (${\Gamma}_{min} > 50$ at $\sim1$ hour) constraints on a GRB in the radio band. Our results demonstrate the importance of rapid radio follow-up observations with long integrations and good sensitivity for detecting the fast-evolving radio emission from short GRBs and probing their reverse shocks., Comment: 11 pages, 3 figures, submitted to ApJL

Quasi-periodic Eruptions (QPEs) are luminous bursts of soft X-rays from the nuclei of galaxies, repeating on timescales of hours to weeks. The mechanism behind these rare systems is uncertain, but most theories involve accretion disks around supermassive black holes (SMBHs), undergoing instabilities or interacting with a stellar object in a close orbit. It has been suggested that this disk could be created when the SMBH disrupts a passing star, implying that many QPEs should be preceded by observable tidal disruption events (TDEs). Two known QPE sources show long-term decays in quiescent luminosity consistent with TDEs, and two observed TDEs have exhibited X-ray flares consistent with individual eruptions. TDEs and QPEs also occur preferentially in similar galaxies. However, no confirmed repeating QPEs have been associated with a spectroscopically confirmed TDE or an optical TDE observed at peak brightness. Here we report the detection of nine X-ray QPEs with a mean recurrence time of approximately 48 hours from AT2019qiz, a nearby and extensively studied optically-selected TDE. We detect and model the X-ray, ultraviolet and optical emission from the accretion disk, and show that an orbiting body colliding with this disk provides a plausible explanation for the QPEs.

Time-domain astrophysics continues to grow rapidly, with the inception of new surveys drastically increasing data volumes. Democratised, distributed approaches to training sets for machine learning classifiers are crucial to make the most of this torrent of discovery -- with citizen science approaches proving effective at meeting these requirements. In this paper, we describe the creation of and the initial results from the $\textit{Kilonova Seekers}$ citizen science project, built to find transient phenomena from the GOTO telescopes in near real-time. $\textit{Kilonova Seekers}$ launched in July 2023 and received over 600,000 classifications from approximately 2,000 volunteers over the course of the LIGO-Virgo-KAGRA O4a observing run. During this time, the project has yielded 20 discoveries, generated a `gold-standard' training set of 17,682 detections for augmenting deep-learned classifiers, and measured the performance and biases of Zooniverse volunteers on real-bogus classification. This project will continue throughout the lifetime of GOTO, pushing candidates at ever-greater cadence, and directly facilitate the next-generation classification algorithms currently in development., Comment: 20 pages, 15 figures. Accepted in MNRAS

The nature of the minute-to-hour long Fast X-ray Transients (FXTs) localised by telescopes such as Chandra, Swift, and XMM-Newton remains mysterious, with numerous models suggested for the events. Here, we report multi-wavelength observations of EP240315a, a 1600 s long transient detected by the Einstein Probe, showing it to have a redshift of z=4.859. We measure a low column density of neutral hydrogen, indicating that the event is embedded in a low-density environment, further supported by direct detection of leaking ionising Lyman-continuum. The observed properties are consistent with EP240315a being a long-duration gamma-ray burst, and these observations support an interpretation in which a significant fraction of the FXT population are lower-luminosity examples of similar events. Such transients are detectable at high redshifts by the Einstein Probe and, in the (near) future, out to even larger distances by SVOM, THESEUS, and Athena, providing samples of events into the epoch of reionisation., Comment: 41 pages, 7 figures, submitted

We present the semi-analytical light curve modelling of 13 supernovae associated with gamma-ray bursts (GRB-SNe) along with two relativistic broad-lined (Ic-BL) SNe without GRBs association (SNe 2009bb and 2012ap), considering millisecond magnetars as central-engine-based power sources for these events. The bolometric light curves of all 15 SNe in our sample are well-regenerated utilising a $\chi^2-$minimisation code, $\texttt{MINIM}$, and numerous parameters are constrained. The median values of ejecta mass ($M_{\textrm{ej}}$), magnetar's initial spin period ($P_\textrm{i}$) and magnetic field ($B$) for GRB-SNe are determined to be $\approx$ 5.2 M$_\odot$, 20.5 ms and 20.1 $\times$ 10$^{14}$ G, respectively. We leverage machine learning (ML) algorithms to comprehensively compare the 3-dimensional parameter space encompassing $M_{\textrm{ej}}$, $P_\textrm{i}$, and $B$ for GRB-SNe determined herein to those of H-deficient superluminous SNe (SLSNe-I), fast blue optical transients (FBOTs), long GRBs (LGRBs), and short GRBs (SGRBs) obtained from the literature. The application of unsupervised ML clustering algorithms on the parameters $M_{\textrm{ej}}$, $P_\textrm{i}$, and $B$ for GRB-SNe, SLSNe-I, and FBOTs yields a classification accuracy of $\sim$95%. Extending these methods to classify GRB-SNe, SLSNe-I, LGRBs, and SGRBs based on $P_\textrm{i}$ and $B$ values results in an accuracy of $\sim$84%. Our investigations show that GRB-SNe and relativistic Ic-BL SNe presented in this study occupy different parameter spaces for $M_{\textrm{ej}}$, $P_\textrm{i}$, and $B$ than those of SLSNe-I, FBOTs, LGRBs and SGRBs. This indicates that magnetars with different $P_\textrm{i}$ and $B$ can give birth to distinct types of transients., Comment: 13 pages, 7 figures, and 3 tables (including appendix). Accepted for publication in MNRAS

Much of what is known of the chemical composition of the universe is based on emission line spectra from star forming galaxies. Emission-based inferences are, nevertheless, model-dependent and they are dominated by light from luminous star forming regions. An alternative and sensitive probe of the metallicity of galaxies is through absorption lines imprinted on the luminous afterglow spectra of long gamma ray bursts (GRBs) from neutral material within their host galaxy. We present results from a JWST/NIRSpec programme to investigate for the first time the relation between the metallicity of neutral gas probed in absorption by GRB afterglows and the metallicity of the star forming regions for the same host galaxy sample. Using an initial sample of eight GRB host galaxies at z=2.1-4.7, we find a tight relation between absorption and emission line metallicities when using the recently proposed $\hat{R}$ metallicity diagnostic (+/-0.2dex). This agreement implies a relatively chemically-homogeneous multi-phase interstellar medium, and indicates that absorption and emission line probes can be directly compared. However, the relation is less clear when using other diagnostics, such as R23 and R3. We also find possible evidence of an elevated N/O ratio in the host galaxy of GRB090323 at z=3.58, consistent with what has been seen in other $z>4$ galaxies. Ultimate confirmation of an enhanced N/O ratio and of the relation between absorption and emission line metallicities will require a more direct determination of the emission line metallicity via the detection of temperature-sensitive auroral lines in our GRB host galaxy sample., Comment: Accepted for publication in MNRAS; 24 pages, 15 figures

Correlations between intrinsic properties of gamma-ray burst (GRB) light curves provide clues to the nature of the central engine, the jet, and a possible means to standardise GRBs for cosmological use. Here we report on the discovery of a correlation between the intrinsic early time luminosity, $L_{G,\rm 10s}$, measured at rest frame 10s, and the average decay rate measured from rest frame 10s onward, $\alpha_{G,\rm avg>10s}$, in a sample of 13 Fermi Large Array Telescope (LAT) long GRB light curves. We note that our selection criteria, in particular the requirement for a redshift to construct luminosity light curves, naturally limits our sample to energetic GRBs. A Spearman's rank correlation gives a coefficient of -0.74, corresponding to a confidence level of 99.6%, indicating that brighter afterglows decay faster than less luminous ones. Assuming a linear relation with $\log(L_{G,\rm 10s})$, we find $\alpha_{G,\rm avg>10s} = -0.31_{-0.09}^{+0.12}\log(L_{G,\rm 10s}) + 14.43_{-5.97}^{+4.55}$. The slope of -0.31 is consistent at $1\sigma$ with previously identified correlations in the optical/UV and X-ray light curves. We speculate that differences in the rate at which energy is released by the central engine or differences in observer viewing angle may be responsible for the correlation., Comment: 8 pages, 4 figures, accepted for publication in MNRAS

Context: GRB 211106A and GRB 211227A are recent gamma-ray bursts (GRBs) with initial X-ray positions suggesting associations with nearby galaxies (z < 0.7). Their prompt emission characteristics indicate GRB 211106A is a short-duration GRB and GRB 211227A is a short GRB with extended emission, likely originating from compact binary mergers. However, classifying solely based on prompt emission can be misleading. Aims: These short GRBs in the local Universe offer opportunities to search for associated kilonova (KN) emission and study host galaxy properties in detail. Methods: We conducted deep optical and NIR follow-up using ESO-VLT FORS2, HAWK-I, and MUSE for GRB 211106A, and ESO-VLT FORS2 and X-Shooter for GRB 211227A, starting shortly after the X-ray afterglow detection. We performed photometric analysis to look for afterglow and KN emissions associated with the bursts, along with host galaxy imaging and spectroscopy. Optical/NIR results were compared with Swift X-Ray Telescope (XRT) and other high-energy data. Results: For both GRBs we placed deep limits to the optical/NIR afterglow and KN emission. Host galaxies were identified: GRB 211106A at photometric z = 0.64 and GRB 211227A at spectroscopic z = 0.228. Host galaxy properties aligned with typical short GRB hosts. We also compared the properties of the bursts with the S-BAT4 sample to further examined the nature of these events. Conclusions: Study of prompt and afterglow phases, along with host galaxy analysis, confirms GRB 211106A as a short GRB and GRB 211227A as a short GRB with extended emission. The absence of optical/NIR counterparts is likely due to local extinction for GRB 211106A and a faint kilonova for GRB 211227A., Comment: Accepted to A&A on 08 August 2023, 21 pages, 24 figures

We present the discovery and extensive follow-up of a remarkable fast-evolving optical transient, AT2022aedm, detected by the Asteroid Terrestrial impact Last Alert Survey (ATLAS). AT2022aedm exhibited a rise time of $9\pm1$ days in the ATLAS $o$-band, reaching a luminous peak with $M_g\approx-22$ mag. It faded by 2 magnitudes in $g$-band during the next 15 days. These timescales are consistent with other rapidly evolving transients, though the luminosity is extreme. Most surprisingly, the host galaxy is a massive elliptical with negligible current star formation. X-ray and radio observations rule out a relativistic AT2018cow-like explosion. A spectrum in the first few days after explosion showed short-lived He II emission resembling young core-collapse supernovae, but obvious broad supernova features never developed; later spectra showed only a fast-cooling continuum and narrow, blue-shifted absorption lines, possibly arising in a wind with $v\approx2700$ km s$^{-1}$. We identify two further transients in the literature (Dougie in particular, as well as AT2020bot) that share similarities in their luminosities, timescales, colour evolution and largely featureless spectra, and propose that these may constitute a new class of transients: luminous fast-coolers (LFCs). All three events occurred in passive galaxies at offsets of $\sim4-10$ kpc from the nucleus, posing a challenge for progenitor models involving massive stars or massive black holes. The light curves and spectra appear to be consistent with shock breakout emission, though usually this mechanism is associated with core-collapse supernovae. The encounter of a star with a stellar mass black hole may provide a promising alternative explanation., Comment: Accepted in ApJL

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

The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs), sources of high-frequency gravitational waves and likely production sites for heavy element nucleosynthesis via rapid neutron capture (the r-process). These heavy elements include some of great geophysical, biological and cultural importance, such as thorium, iodine and gold. Here we present observations of the exceptionally bright gamma-ray burst GRB 230307A. We show that GRB 230307A belongs to the class of long-duration gamma-ray bursts associated with compact object mergers, and contains a kilonova similar to AT2017gfo, associated with the gravitational-wave merger GW170817. We obtained James Webb Space Telescope mid-infrared (mid-IR) imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns which we interpret as tellurium (atomic mass A=130), and a very red source, emitting most of its light in the mid-IR due to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r-process elements across a broad atomic mass range and play a central role in heavy element nucleosynthesis across the Universe., Comment: Submitted. Comments welcome! Nature (2023)

Luminous Fast Blue Optical Transients (LFBOTs) - the prototypical example being AT2018cow - are a rare class of events whose origins are poorly understood. They are characterised by rapid evolution, featureless blue spectra at early times, and luminous X-ray and radio emission. LFBOTs thus far have been found exclusively at small projected offsets from star-forming host galaxies. We present Hubble Space Telescope, Gemini, Chandra and Very Large Array observations of a new LFBOT, AT2023fhn. The Hubble Space Telescope data reveal a large offset (greater than 3.5 half-light radii) from the two closest galaxies, both at a redshift of 0.24. The location of AT2023fhn is in stark contrast with previous events, and demonstrates that LFBOTs can occur in a range of galactic environments., Comment: Accepted for publication in MNRASL. 7 pages, 4 figures, 2 tables

GRB 201015A is a peculiarly low luminosity, spectrally soft gamma-ray burst (GRB), with $T_{\rm 90} = 9.8 \pm 3.5$ s (time interval of detection of 90\% of photons from the GRB), and an associated supernova (likely to be type Ic or Ic-BL). GRB 201015A has an isotropic energy $E_{\gamma,\rm iso} = 1.75 ^{+0.60} _{-0.53} \times 10^{50}$ erg, and photon index $\Gamma = 3.00 ^{+0.50} _{-0.42}$ (15-150 keV). It follows the Amati relation, a correlation between $E_{\gamma,\rm iso}$ and spectral peak energy $E_{\rm p}$ followed by long GRBs. It appears exceptionally soft based on $\Gamma$, the hardness ratio of HR = $0.47 \pm 0.24$, and low-$E_{\rm p}$, so we have compared it to other GRBs sharing these properties. These events can be explained by shock breakout, poorly collimated jets, and off-axis viewing. Follow-up observations of the afterglow taken in the X-ray, optical, and radio, reveal a surprisingly late flattening in the X-ray from $t = (2.61 \pm 1.27)\times 10^4$ s to $t = 1.67 ^{+1.14} _{-0.65} \times 10^6$ s. We fit the data to closure relations describing the synchrotron emission, finding the electron spectral index to be $p = 2.42 ^{+0.44} _{-0.30}$, and evidence of late-time energy injection with coefficient $q = 0.24 ^{+0.24} _{-0.18}$. The jet half opening angle lower limit ($\theta_{j} \ge 16^{\circ}$) is inferred from the non-detection of a jet break. The launch of SVOM and Einstein Probe in 2023, should enable detection of more low luminosity events like this, providing a fuller picture of the variety of GRBs., Comment: 15 pages, 4 figures

We present a semi-analytic model for predicting kilonova light curves from the mergers of neutron stars with black holes (NSBH). The model is integrated into the MOSFiT platform, and can generate light curves from input binary properties and nuclear equation-of-state considerations, or incorporate measurements from gravitational wave (GW) detectors to perform multi-messenger parameter estimation. The rapid framework enables the generation of NSBH kilonova distributions from binary populations, light curve predictions from GW data, and statistically meaningful comparisons with an equivalent BNS model in MOSFiT. We investigate a sample of kilonova candidates associated with cosmological short gamma-ray bursts, and demonstrate that they are broadly consistent with being driven by NSBH systems, though most have limited data. We also perform fits to the very well sampled GW170817, and show that the inability of an NSBH merger to produce lanthanide-poor ejecta results in a significant underestimate of the early (< 2 days) optical emission. Our model indicates that NSBH-driven kilonovae may peak up to a week after merger at optical wavelengths for some observer angles. This demonstrates the need for early coverage of emergent kilonovae in cases where the GW signal is either ambiguous or absent; they likely cannot be distinguished from BNS mergers by the light curves alone from ~2 days after the merger. We also discuss the detectability of our model kilonovae with the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST)., Comment: 14 pages, 6 figures, 2 tables. Accepted for publication in MNRAS. This is the author's final submitted version. The model code is available through MOSFiT at https://github.com/guillochon/MOSFiT

Long and short gamma-ray bursts (GRBs), canonically separated at around 2 seconds duration, are associated with different progenitors: the collapse of a massive star and the merger of two compact objects, respectively. GRB 191019A was a long GRB ($T_{90}\sim64$ s). Despite the relatively small redshift z=0.248 and HST followup observations, an accompanying supernova was not detected. In addition, the host galaxy did not have significant star formation activity. Here we propose that GRB 191019A was produced by a binary compact merger, whose prompt emission was stretched in time by the interaction with a dense external medium. This would be expected if the burst progenitor was located in the disk of an active galactic nucleus, as supported by the burst localization close to the center of its host galaxy. We show that the light curve of GRB 191019A can be well modeled by a burst of intrinsic duration t=1.1 s and of energy $E_{\rm{iso}}=10^{51}$ erg seen moderately off-axis, exploding in a medium of density $10^7-10^8$ cm$^{-3}$. The double-peaked light curve carries the telltale features predicted for GRBs in high-density media, where the first peak is produced by the photosphere, and the second by the overlap of reverse shocks that take place before the internal shocks could happen. This would make GRB 191019A the first confirmed stellar explosion from within an accretion disk, with important implications for the formation and evolution of stars in accretion flows and for gravitational waves source populations., Comment: 6 pages, 2 figures. Accepted for publication in ApJL

The majority of long duration ($>2$ s) gamma-ray bursts (GRBs) are believed to arise from the collapse of massive stars \cite{Hjorth+03}, with a small proportion created from the merger of compact objects. Most of these systems are likely formed via standard stellar evolution pathways. However, it has long been thought that a fraction of GRBs may instead be an outcome of dynamical interactions in dense environments, channels which could also contribute significantly to the samples of compact object mergers detected as gravitational wave sources. Here we report the case of GRB 191019A, a long GRB (T_90 = 64.4 +/- 4.5 s) which we pinpoint close (<100 pc projected) to the nucleus of an ancient (>1~Gyr old) host galaxy at z=0.248. The lack of evidence for star formation and deep limits on any supernova emission make a massive star origin difficult to reconcile with observations, while the timescales of the emission rule out a direct interaction with the supermassive black hole in the nucleus of the galaxy, We suggest that the most likely route for progenitor formation is via dynamical interactions in the dense nucleus of the host, consistent with the centres of such galaxies exhibiting interaction rates up to two orders of magnitude larger than typical field galaxies. The burst properties could naturally be explained via compact object mergers involving white dwarfs (WD), neutron stars (NS) or black holes (BH). These may form dynamically in dense stellar clusters, or originate in a gaseous disc around the supermassive black hole. Future electromagnetic and gravitational-wave observations in tandem thus offer a route to probe the dynamical fraction and the details of dynamical interactions in galactic nuclei and other high density stellar systems., Comment: Accepted to Nature Astronomy. This is the submitted version and will differ from the published version due to modifications in the refereeing process

Context: The extreme luminosity of gamma-ray bursts (GRBs) makes them powerful beacons for studies of the distant Universe. The most luminous bursts are typically detected at moderate/high redshift, where the volume for seeing such rare events is maximized and the star-formation activity is greater than at z = 0. For distant events, not all observations are feasible, such as at TeV energies. Aims: Here we present a spectroscopic redshift measurement for the exceptional GRB 221009A, the brightest GRB observed to date with emission extending well into the TeV regime. Methods: We used the X-shooter spectrograph at the ESO Very Large Telescope (VLT) to obtain simultaneous optical to near-IR spectroscopy of the burst afterglow 0.5 days after the explosion. Results: The spectra exhibit both absorption and emission lines from material in a host galaxy at z = 0.151. Thus GRB 221009A was a relatively nearby burst with a luminosity distance of 745 Mpc. Its host galaxy properties (star-formation rate and metallicity) are consistent with those of LGRB hosts at low redshift. This redshift measurement yields information on the energy of the burst. The inferred isotropic energy release, $E_{\rm iso} > 5 \times 10^{54}$ erg, lies at the high end of the distribution, making GRB 221009A one of the nearest and also most energetic GRBs observed to date. We estimate that such a combination (nearby as well as intrinsically bright) occurs between once every few decades to once per millennium., Comment: 9 pages, 4 figures, Astronomy & Astrophysics, accepted

We present JWST and Hubble Space Telescope (HST) observations of the afterglow of GRB 221009A, the brightest gamma-ray burst (GRB) ever observed. This includes the first mid-IR spectra of any GRB, obtained with JWST/NIRSPEC (0.6-5.5 micron) and MIRI (5-12 micron), 12 days after the burst. Assuming that the intrinsic spectral slope is a single power-law, with $F_{\nu} \propto \nu^{-\beta}$, we obtain $\beta \approx 0.35$, modified by substantial dust extinction with $A_V = 4.9$. This suggests extinction above the notional Galactic value, possibly due to patchy extinction within the Milky Way or dust in the GRB host galaxy. It further implies that the X-ray and optical/IR regimes are not on the same segment of the synchrotron spectrum of the afterglow. If the cooling break lies between the X-ray and optical/IR, then the temporal decay rates would only match a post jet-break model, with electron index $p<2$, and with the jet expanding into a uniform ISM medium. The shape of the JWST spectrum is near-identical in the optical/nIR to X-shooter spectroscopy obtained at 0.5 days and to later time observations with HST. The lack of spectral evolution suggests that any accompanying supernova (SN) is either substantially fainter or bluer than SN 1998bw, the proto-type GRB-SN. Our HST observations also reveal a disc-like host galaxy, viewed close to edge-on, that further complicates the isolation of any supernova component. The host galaxy appears rather typical amongst long-GRB hosts and suggests that the extreme properties of GRB 221009A are not directly tied to its galaxy-scale environment., Comment: Accepted for publication to the Astrophysical Journal Letters for the GRB 221009A Special Issue. The results of this paper are under press embargo until March 28, 18 UT. 19 pages, 8 figures, 2 tables

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

Developing an effective automatic classifier to separate genuine sources from artifacts is essential for transient follow-ups in wide-field optical surveys. The identification of transient detections from the subtraction artifacts after the image differencing process is a key step in such classifiers, known as real-bogus classification problem. We apply a self-supervised machine learning model, the deep-embedded self-organizing map (DESOM) to this "real-bogus" classification problem. DESOM combines an autoencoder and a self-organizing map to perform clustering in order to distinguish between real and bogus detections, based on their dimensionality-reduced representations. We use 32x32 normalized detection thumbnails as the input of DESOM. We demonstrate different model training approaches, and find that our best DESOM classifier shows a missed detection rate of 6.6% with a false positive rate of 1.5%. DESOM offers a more nuanced way to fine-tune the decision boundary identifying likely real detections when used in combination with other types of classifiers, for example built on neural networks or decision trees. We also discuss other potential usages of DESOM and its limitations.

We present the results from multi-wavelength observations of a transient discovered during the follow-up of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN2019wxt, a young transient in a galaxy whose sky position (in the 80\% GW contour) and distance ($\sim$150\,Mpc) were plausibly compatible with the localisation uncertainty of the GW event. Initially, the transient's tightly constrained age, its relatively faint peak magnitude ($M_i \sim -16.7$\,mag) and the $r-$band decline rate of $\sim 1$\,mag per 5\,days appeared suggestive of a compact binary merger. However, SN2019wxt spectroscopically resembled a type Ib supernova, and analysis of the optical-near-infrared evolution rapidly led to the conclusion that while it could not be associated with S191213g, it nevertheless represented an extreme outcome of stellar evolution. By modelling the light curve, we estimated an ejecta mass of $\sim 0.1\,M_\odot$, with $^{56}$Ni comprising $\sim 20\%$ of this. We were broadly able to reproduce its spectral evolution with a composition dominated by helium and oxygen, with trace amounts of calcium. We considered various progenitors that could give rise to the observed properties of SN2019wxt, and concluded that an ultra-stripped origin in a binary system is the most likely explanation. Disentangling electromagnetic counterparts to GW events from transients such as SN2019wxt is challenging: in a bid to characterise the level of contamination, we estimated the rate of events with properties comparable to those of SN2019wxt and found that $\sim 1$ such event per week can occur within the typical GW localisation area of O4 alerts out to a luminosity distance of 500\,Mpc, beyond which it would become fainter than the typical depth of current electromagnetic follow-up campaigns., Comment: By the ENGRAVE collaboration (engrave-eso.org). 35 pages, 20 figures, final version accepted by A&A

A black-hole neutron-star binary merger can lead to an electromagnetic counterpart called a kilonova if the neutron star is disrupted prior to merger. The observability of a kilonova depends on the amount of neutron star ejecta, which is sensitive to the aligned component of the black hole spin. We explore the dependence of the ejected mass on two main mechanisms that provide high black hole spin in isolated stellar binaries. When the black hole inherits a high spin from a Wolf-Rayet star that was born with least $\sim$ 10% of its breakup spin under weak stellar core-envelope coupling, relevant for all formation pathways, the median of the ejected mass is $\gtrsim$ $10^{-2}$ M$_{\odot}$. Though only possible for certain formation pathways, similar ejected mass results when the black hole accretes $\gtrsim$ 20% of its companion's envelope to gain a high spin. Together, these signatures suggest that a population analysis of black-hole neutron-star binary mergers with observed kilonovae may help distinguish between mechanisms for spin and possible formation pathways. We show that these kilonovae will be difficult to detect with current capabilities, but that future facilities, such as the Vera Rubin Observatory, can do so even if the aligned dimensionless spin of the black hole is as low as $\sim$ 0.2. Our model predicts kilonovae as bright as $M_i$ $\sim$ -14.5 for an aligned black hole spin of $\sim$ 0.9 and mass ratio Q = 3.6.

The temporal and spectral evolution of gamma-ray burst (GRB) afterglows can be used to infer the density and density profile of the medium through which the shock is propagating. In long-duration (core-collapse) GRBs, the circumstellar medium (CSM) is expected to resemble a wind-blown bubble, with a termination shock separating the stellar wind and the interstellar medium (ISM). A long standing problem is that flat density profiles, indicative of the ISM, are often found at lower radii than expected for a massive star progenitor. Furthermore, the presence of both wind-like environments at high radii and ISM-like environments at low radii remains a mystery. In this paper, we perform a 'CSM population synthesis' with long GRB progenitor stellar evolution models. Analytic results for the evolution of wind blown bubbles are adjusted through comparison with a grid of 2D hydrodynamical simulations. Predictions for the emission radii, ratio of ISM to wind-like environments, wind and ISM densities are compared with the largest sample of afterglow-derived parameters yet compiled, which we make available for the community. We find that high ISM densities around 1000/cm3 best reproduce observations. If long GRBs instead occur in typical ISM densities of approximately 1/cm3, then the discrepancy between theory and observations is shown to persist at a population level. We discuss possible explanations for the origin of variety in long GRB afterglows, and for the overall trend of CSM modelling to over-predict the termination shock radius., Comment: 22 pages, 15 figures, 4 tables. Accepted for publication in MNRAS

GRB 201015A is a long-duration Gamma-Ray Burst (GRB) which was detected at very high energies (> 100 GeV) using the MAGIC telescopes. If confirmed, this would be the fifth and least luminous GRB ever detected at this energies. We performed a radio follow-up of GRB 201015A over twelve different epochs, from 1.4 to 117 days post-burst, with the Karl G. Jansky Very Large Array, e-MERLIN and the European VLBI Network. We included optical and X-rays observations, performed with the Multiple Mirror Telescope and the Chandra X-ray Observatory respectively, together with publicly available data. We detected a point-like transient, consistent with the position of GRB 201015A until 23 and 47 days post-burst at 1.5 and 5 GHz, respectively. The source was detected also in both optical (1.4 and 2.2 days post-burst) and X-ray (8.4 and 13.6 days post-burst) observations. The multi-wavelength afterglow light curves can be explained with the standard model for a GRB seen on-axis, which expands and decelerates into a medium with a homogeneous density, while a circumburst medium with a wind-like profile is disfavoured. Notwithstanding the high resolution provided by the VLBI, we could not pinpoint any expansion or centroid displacement of the outflow. If the GRB is seen at the viewing angle which maximises the apparent velocity, we estimate that the Lorentz factor for the possible proper motion is $\Gamma_{\alpha}$ < 40 in right ascension and $\Gamma_{\delta}$ < 61 in declination. On the other hand, if the GRB is seen on-axis, the size of the afterglow is <5 pc and <16 pc at 25 and 47 days. Finally, the early peak in the optical light curve suggests the presence of a reverse shock component before 0.01 days from the burst., Comment: Accepted for publication in Astronomy & Astrophysics, 11 pages, 3 figures

For decades, gamma-ray bursts (GRBs) have been broadly divided into `long'- and `short'-duration bursts, lasting more or less than 2s, respectively. However, this dichotomy does not map perfectly to the two progenitor channels that are known to produce GRBs -- the merger of compact objects (merger-GRBs) or the collapse of massive stars (collapsar-GRBs). In particular, the merger-GRBs population may also include bursts with a short, hard $\lesssim$2s spike and subsequent longer, softer extended emission (EE). The recent discovery of a kilonova -- the radioactive glow of heavy elements made in neutron star mergers -- in the 50s-duration GRB 211211A further demonstrates that mergers can drive long, complex GRBs that mimic the collapsar population. Here we present a detailed temporal and spectral analysis of the high-energy emission of GRB 211211A. We demonstrate that the emission has a purely synchrotron origin, with both the peak and cooling frequencies moving through the $\gamma$-ray band down to the X-rays, and that the rapidly-evolving spectrum drives the EE signature at late times. The identification of such spectral evolution in a merger-GRB opens avenues for diagnostics of the progenitor type., Comment: Author's final submitted version. 6 figures, 5 tables. The Supplementary Information .tex file is included

Here, we report the discovery of a kilonova associated with the nearby (350 Mpc) minute-duration GRB 211211A. In tandem with deep optical limits that rule out the presence of an accompanying supernova to $M_I > -13$ mag at 17.7 days post-burst, the identification of a kilonova confirms that this burst's progenitor was a compact object merger. While the spectrally softer tail in GRB 211211A's gamma-ray light curve is reminiscent of previous extended emission short GRBs (EE-SGRBs), its prompt, bright spikes last $\gtrsim 12$ s, separating it from past EE-SGRBs. GRB 211211A's kilonova has a similar luminosity, duration and color to AT2017gfo, the kilonova found in association with the gravitational wave (GW)-detected binary neutron star (BNS) merger GW170817. We find that the merger ejected $\approx 0.04 M_{\odot}$ of r-process-rich material, and is consistent with the merger of two neutron stars (NSs) with masses close to the canonical $1.4 M_{\odot}$. This discovery implies that GRBs with long, complex light curves can be spawned from compact object merger events and that a population of kilonovae following GRBs with durations $\gg 2$ s should be accounted for in calculations of the NS merger r-process contribution and rate. At 350 Mpc, the current network of GW interferometers at design sensitivity would have detected the merger precipitating GRB 211211A, had it been operating at the time of the event. Further searches for GW signals coincident with long GRBs are therefore a promising route for future multi-messenger astronomy., Comment: Submitted. 69 pages, 11 figures, 3 tables

The Gravitational-wave Optical Transient Observer (GOTO) is an array of wide-field optical telescopes, designed to exploit new discoveries from the next generation of gravitational wave detectors (LIGO, Virgo, KAGRA), study rapidly evolving transients, and exploit multi-messenger opportunities arising from neutrino and very high energy gamma-ray triggers. In addition to a rapid response mode, the array will also perform a sensitive, all-sky transient survey with few day cadence. The facility features a novel, modular design with multiple 40-cm wide-field reflectors on a single mount. In June 2017 the GOTO collaboration deployed the initial project prototype, with 4 telescope units, at the Roque de los Muchachos Observatory (ORM), La Palma, Canary Islands. Here we describe the deployment, commissioning, and performance of the prototype hardware, and discuss the impact of these findings on the final GOTO design. We also offer an initial assessment of the science prospects for the full GOTO facility that employs 32 telescope units across two sites., Comment: 19 pages, 16 Figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Compact binaries such as double neutron stars or a neutron star paired with a black-hole, are strong sources of gravitational waves during coalescence and also the likely progenitors of various electromagnetic phenomena, notably short-duration gamma-ray bursts (SGRBs), and kilonovae. In this work, we generate populations of synthetic binaries and place them in galaxies from the large-scale hydrodynamical galaxy evolution simulation EAGLE. With our zELDA code, binaries are seeded in proportion to star formation rate, and we follow their evolution to merger using both the BPASS and COSMIC binary stellar evolution codes. We track their dynamical evolution within their host galaxy potential, to estimate the galactocentric distance at the time of the merger. Finally, we apply observational selection criteria to allow comparison of this model population with the legacy sample of SGRBs. We find a reasonable agreement with the redshift distribution (peaking at $0.526$)., Comment: Accepted for publication in MNRAS, 21 Pages (6 Tables, 14 Figures), 14 Pages Appendix (4 Tables, 16 Figures)

The typical detection rate of $\sim1$ gamma-ray burst (GRB) per day by the \emph{Fermi} Gamma-ray Burst Monitor (GBM) provides a valuable opportunity to further our understanding of GRB physics. However, the large uncertainty of the \emph{Fermi} localization typically prevents rapid identification of multi-wavelength counterparts. We report the follow-up of 93 \emph{Fermi} GRBs with the Gravitational-wave Optical Transient Observer (GOTO) prototype on La Palma. We selected 53 events (based on favourable observing conditions) for detailed analysis, and to demonstrate our strategy of searching for optical counterparts. We apply a filtering process consisting of both automated and manual steps to 60\,085 candidates initially, rejecting all but 29, arising from 15 events. With $\approx3$ GRB afterglows expected to be detectable with GOTO from our sample, most of the candidates are unlikely to be related to the GRBs. Since we did not have multiple observations for those candidates, we cannot confidently confirm the association between the transients and the GRBs. Our results show that GOTO can effectively search for GRB optical counterparts thanks to its large field of view of $\approx40$ square degrees and its depth of $\approx20$ mag. We also detail several methods to improve our overall performance for future follow-up programs of \emph{Fermi} GRBs.

The identification of the first confirmed neutron star - black hole (NS-BH) binary mergers by the LIGO, Virgo and KAGRA collaboration provides the opportunity to investigate the properties of the early sample of confirmed and candidate events. Here, we focus primarily on the tilt angle of the black hole's spin relative to the orbital angular momentum vector of the binary, and the implications for the physical processes that determine this tilt. The posterior tilt distributions of GW200115 and the candidate events GW190426_152155 and GW190917_114630 peak at significantly anti-aligned orientations (though display wide distributions). Producing these tilts through isolated binary evolution would require stronger natal kicks than are typically considered (and preferentially-polar kicks would be ruled out), and/or an additional source of tilt such as stable mass transfer. The early sample of NS-BH events are less massive than expected for classical formation channels, and may provide evidence for efficient mass transfer that results in the merger of more massive NS-BH binaries before their evolution to the compact phase is complete. We predict that future gravitational-wave detections of NS-BH events will continue to display total binary masses of $\approx 7$ M$_{\odot}$ and mass ratios of $q \sim 3$ if this interpretation is correct. Conversely, the high mass of the candidate GW191219_163120 suggests a dynamical capture origin. Large tilts in a significant fraction of merging NS-BH systems would weaken the prospects for electromagnetic detection. However, EM observations, including non-detections, can significantly tighten the constraints on spin and mass ratio., Comment: 8 pages, 3 figures. Accepted for publication in MNRAS

The advent of wide-field sky surveys has led to the growth of transient and variable source discoveries. The data deluge produced by these surveys has necessitated the use of machine learning (ML) and deep learning (DL) algorithms to sift through the vast incoming data stream. A problem that arises in real-world applications of learning algorithms for classification is imbalanced data, where a class of objects within the data is underrepresented, leading to a bias for over-represented classes in the ML and DL classifiers. We present a recurrent neural network (RNN) classifier that takes in photometric time-series data and additional contextual information (such as distance to nearby galaxies and on-sky position) to produce real-time classification of objects observed by the Gravitational-wave Optical Transient Observer (GOTO), and use an algorithm-level approach for handling imbalance with a focal loss function. The classifier is able to achieve an Area Under the Curve (AUC) score of 0.972 when using all available photometric observations to classify variable stars, supernovae, and active galactic nuclei. The RNN architecture allows us to classify incomplete light curves, and measure how performance improves as more observations are included. We also investigate the role that contextual information plays in producing reliable object classification., Comment: 16 pages, 12 figures, to be published in Monthly Notices of the Royal Astronomical Society

We have adapted the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Science Pipelines to process data from the Gravitational-Wave Optical Transient Observer (GOTO) prototype. In this paper, we describe how we used the Rubin Observatory LSST Science Pipelines to conduct forced photometry measurements on nightly GOTO data. By comparing the photometry measurements of sources taken on multiple nights, we find that the precision of our photometry is typically better than 20~mmag for sources brighter than 16 mag. We also compare our photometry measurements against colour-corrected PanSTARRS photometry, and find that the two agree to within 10~mmag (1$\sigma$) for bright (i.e., $\sim14^{\rm th}$~mag) sources to 200~mmag for faint (i.e., $\sim18^{\rm th}$~mag) sources. Additionally, we compare our results to those obtained by GOTO's own in-house pipeline, {\sc gotophoto}, and obtain similar results. Based on repeatability measurements, we measure a $5\sigma$ L-band survey depth of between 19 and 20 magnitudes, depending on observing conditions. We assess, using repeated observations of non-varying standard SDSS stars, the accuracy of our uncertainties, which we find are typically overestimated by roughly a factor of two for bright sources (i.e., $<15^{\rm th}$~mag), but slightly underestimated (by roughly a factor of 1.25) for fainter sources ($>17^{\rm th}$~mag). Finally, we present lightcurves for a selection of variable sources, and compare them to those obtained with the Zwicky Transient Factory and GAIA. Despite the Rubin Observatory LSST Science Pipelines still undergoing active development, our results show that they are already delivering robust forced photometry measurements from GOTO data., Comment: Accepted for publication in PASA

Large-scale sky surveys have played a transformative role in our understanding of astrophysical transients, only made possible by increasingly powerful machine learning-based filtering to accurately sift through the vast quantities of incoming data generated. In this paper, we present a new real-bogus classifier based on a Bayesian convolutional neural network that provides nuanced, uncertainty-aware classification of transient candidates in difference imaging, and demonstrate its application to the datastream from the GOTO wide-field optical survey. Not only are candidates assigned a well-calibrated probability of being real, but also an associated confidence that can be used to prioritise human vetting efforts and inform future model optimisation via active learning. To fully realise the potential of this architecture, we present a fully-automated training set generation method which requires no human labelling, incorporating a novel data-driven augmentation method to significantly improve the recovery of faint and nuclear transient sources. We achieve competitive classification accuracy (FPR and FNR both below 1%) compared against classifiers trained with fully human-labelled datasets, whilst being significantly quicker and less labour-intensive to build. This data-driven approach is uniquely scalable to the upcoming challenges and data needs of next-generation transient surveys. We make our data generation and model training codes available to the community., Comment: 17 pages, 12 figures, resubmitted to MNRAS following reviewer comments

The past few decades have seen the burgeoning of wide field, high cadence surveys, the most formidable of which will be the Legacy Survey of Space and Time (LSST) to be conducted by the Vera C. Rubin Observatory. So new is the field of systematic time-domain survey astronomy, however, that major scientific insights will continue to be obtained using smaller, more flexible systems than the LSST. One such example is the Gravitational-wave Optical Transient Observer (GOTO), whose primary science objective is the optical follow-up of Gravitational Wave events. The amount and rate of data production by GOTO and other wide-area, high-cadence surveys presents a significant challenge to data processing pipelines which need to operate in near real-time to fully exploit the time-domain. In this study, we adapt the Rubin Observatory LSST Science Pipelines to process GOTO data, thereby exploring the feasibility of using this "off-the-shelf" pipeline to process data from other wide-area, high-cadence surveys. In this paper, we describe how we use the LSST Science Pipelines to process raw GOTO frames to ultimately produce calibrated coadded images and photometric source catalogues. After comparing the measured astrometry and photometry to those of matched sources from PanSTARRS DR1, we find that measured source positions are typically accurate to sub-pixel levels, and that measured L-band photometries are accurate to $\sim50$ mmag at $m_L\sim16$ and $\sim200$ mmag at $m_L\sim18$. These values compare favourably to those obtained using GOTO's primary, in-house pipeline, GOTOPHOTO, in spite of both pipelines having undergone further development and improvement beyond the implementations used in this study. Finally, we release a generic "obs package" that others can build-upon should they wish to use the LSST Science Pipelines to process data from other facilities., Comment: 17 pages, 13 figures. Accepted for publication in PASA

The amount of observational data produced by time-domain astronomy is exponentially in-creasing. Human inspection alone is not an effective way to identify genuine transients fromthe data. An automatic real-bogus classifier is needed and machine learning techniques are commonly used to achieve this goal. Building a training set with a sufficiently large number of verified transients is challenging, due to the requirement of human verification. We presentan approach for creating a training set by using all detections in the science images to be thesample of real detections and all detections in the difference images, which are generated by the process of difference imaging to detect transients, to be the samples of bogus detections. This strategy effectively minimizes the labour involved in the data labelling for supervised machine learning methods. We demonstrate the utility of the training set by using it to train several classifiers utilizing as the feature representation the normalized pixel values in 21-by-21pixel stamps centered at the detection position, observed with the Gravitational-wave Optical Transient Observer (GOTO) prototype. The real-bogus classifier trained with this strategy can provide up to 95% prediction accuracy on the real detections at a false alarm rate of 1%., Comment: 9 pages, 8 figures

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200~Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$~Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma \& Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers., Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in MNRAS. Author's final submitted version

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event., Comment: 52 pages, revised version now accepted for publication in A&A. Abstract abridged to meet arXiv requirements

Short gamma-ray bursts (SGRBs) are now known to be the product of the merger of two compact objects. However, two possible formation channels exist: neutron star -- neutron star (NS -- NS) or NS -- black hole (BH). The landmark SGRB 170817A provided evidence for the NS -- NS channel, thanks to analysis of its gravitational wave signal. We investigate the complete population of SGRBs with an associated redshift (39 events), and search for any divisions that may indicate that a NS -- BH formation channel also contributes. Though no conclusive dichotomy is found, we find several lines of evidence that tentatively support the hypothesis that SGRBs with extended emission (EE; 7 events) constitute the missing merger population: they are unique in the large energy band-sensitivity of their durations, and have statistically distinct energies and host galaxy offsets when compared to regular (non-EE) SGRBs. If this is borne out via future gravitational wave detections it will conclusively disprove the magnetar model for SGRBs. Furthermore, we identify the first statistically significant anti-correlation between the offsets of SGRBs from their host galaxies and their prompt emission energies., Comment: 23 pages, 9 figures, 2 tables. 1 appendix (2 pages, 1 table). Author's final accepted version, to be published in ApJ

A small fraction of candidate tidal disruption events (TDEs) show evidence of powerful relativistic jets, which are particularly pronounced at radio wavelengths, and likely contribute non-thermal emission at a wide range of wavelengths. A non-thermal emission component can be diagnosed using linear polarimetry, even when the total received light is dominated by emission from an accretion disk or disk outflow. In this paper we present Very Large Telescope (VLT) measurements of the linear polarisation of the optical light of jetted TDE Swift J2058+0516. This is the second jetted TDE studied in this manner, after Swift J1644+57. We find evidence of non-zero optical linear polarisation, P_V ~ 8%, a level very similar to the near-infrared polarimetry of Swift J1644+57. These detections provide an independent test of the emission mechanisms of the multiwavelength emission of jetted tidal disruption events., Comment: 7 pages, 4 figures, MNRAS

We present X-ray and optical observations of the short duration gamma-ray burst GRB 071227 and its host at $z=0.381$, obtained using \textit{Swift}, Gemini South and the Very Large Telescope. We identify a short-lived and moderately bright optical transient, with flux significantly in excess of that expected from a simple extrapolation of the X-ray spectrum at 0.2-0.3 days after burst. We fit the SED with afterglow models allowing for high extinction and thermal emission models that approximate a kilonova to assess the excess' origins. While some kilonova contribution is plausible, it is not favoured due to the low temperature and high luminosity required, implying superluminal expansion and a large ejecta mass of $\sim 0.1$ M$_{\odot}$. We find, instead, that the transient is broadly consistent with power-law spectra with additional dust extinction of $E(B-V)\sim0.4$ mag, although a possibly thermal excess remains in the \textit{z}-band. We investigate the host, a spiral galaxy with an edge-on orientation, resolving its spectrum along its major axis to construct the galaxy rotation curve and analyse the star formation and chemical properties. The integrated host emission shows evidence for high extinction, consistent with the afterglow findings. The metallicity and extinction are consistent with previous studies of this host and indicate the galaxy is a typical, but dusty, late-type SGRB host., Comment: 16 pages, 13 figures, accepted to MNRAS

The number of long gamma-ray bursts (GRBs) known to have occurred in the distant Universe (z greater than 5) is small (approx 15), however these events provide a powerful way of probing star formation at the onset of galaxy evolution. In this paper, we present the case for GRB100205A being a largely overlooked high-redshift event. While initially noted as a high-z candidate, this event and its host galaxy have not been explored in detail. By combining optical and near-infrared Gemini afterglow imaging (at t less than 1.3 days since burst) with deep late-time limits on host emission from the Hubble Space Telescope, we show that the most likely scenario is that GRB100205A arose in the redshift range 4-8. GRB100205A is an example of a burst whose afterglow, even at 1 hour post-burst, could only be identified by 8m class IR observations, and suggests that such observations of all optically dark bursts may be necessary to significantly enhance the number of high-redshift GRBs known., Comment: Accepted for publication in MNRAS

We present post-jet-break \textit{HST}, VLA and \textit{Chandra} observations of the afterglow of the long $\gamma$-ray bursts GRB 160625B (between 69 and 209 days) and GRB 160509A (between 35 and 80 days). We calculate the post-jet-break decline rates of the light curves, and find the afterglow of GRB 160625B inconsistent with a simple $t^{-3/4}$ steepening over the break, expected from the geometric effect of the jet edge entering our line of sight. However, the favored optical post-break decline ($f_{\nu} \propto t^{-1.96 \pm 0.07}$) is also inconsistent with the $f_{\nu} \propto t^{-p}$ decline (where $p \approx 2.3$ from the pre-break light curve), which is expected from exponential lateral expansion of the jet; perhaps suggesting lateral expansion that only affects a fraction of the jet. The post-break decline of GRB 160509A is consistent with both the $t^{-3/4}$ steepening and with $f_{\nu} \propto t^{-p}$. We also use {\sc boxfit} to fit afterglow models to both light curves and find both to be energetically consistent with a millisecond magnetar central engine, although the magnetar parameters need to be extreme (i.e. $E \sim 3 \times 10^{52}$ erg). Finally, the late-time radio light curves of both afterglows are not reproduced well by {\sc boxfit} and are inconsistent with predictions from the standard jet model; instead both are well represented by a single power law decline (roughly $f_{\nu} \propto t^{-1}$) with no breaks. This requires a highly chromatic jet break ($t_{j,\mathrm{radio}} > 10 \times t_{j,\mathrm{optical}}$) and possibly a two-component jet for both bursts., Comment: 19 pages, 9 figures. Revised version; accepted for publication in ApJ

We report our identification of the optical afterglow and host galaxy of the short-duration gamma-ray burst GRB 160821B. The spectroscopic redshift of the host is $z=0.162$, making it one of the lowest redshift sGRBs identified by Swift. Our intensive follow-up campaign using a range of ground-based facilities as well as HST, XMM and Swift, shows evidence for a late-time excess of optical and near-infrared emission in addition to a complex afterglow. The afterglow light-curve at X-ray frequencies reveals a narrow jet, $\theta_j\sim1.9^{+0.10}_{-0.03}$ deg, that is refreshed at $>1$ day post-burst by a slower outflow with significantly more energy than the initial outflow that produced the main GRB. Observations of the 5 GHz radio afterglow shows a reverse shock into a mildly magnetised shell. The optical and near-infrared excess is fainter than AT2017gfo associated with GW170817, and is well explained by a kilonova with dynamic ejecta mass $M_{\rm dyn}=(1.0\pm0.6)\times10^{-3}$ M$_{\odot}$ and a secular (postmerger) ejecta mass with $M_{\rm pm}=(1.0\pm0.6)\times10^{-2}$ M$_\odot$, consistent with a binary neutron star merger resulting in a short-lived massive neutron star. This optical and near-infrared dataset provides the best-sampled kilonova light-curve without a gravitational wave trigger to date., Comment: 17 pages, 6 figures, Version accepted by ApJ