Search

Over time, core-collapse supernova (CCSN) spectra become redder due to dust formation and cooling of the SN ejecta. A UV detection of a CCSN at late times thus indicates an additional physical process such as interaction between the SN ejecta and the circumstellar material, or viewing down to the central engine of the explosion. Both these models have been proposed to explain the peculiar transient AT2018cow, a luminous fast blue optical transient that has been detected in the UV 2-4 years after the event with only marginal fading over this time period. To identify if the late-time UV detection of AT2018cow could indicate that it is a CCSN, we investigate if CCSNe are detected in the UV between 2-5 years after the explosion. We use a sample of 51 nearby (z<0.065) CCSNe observed with the Hubble Space Telescope within 2-5 years of discovery. We measure their brightness, or determine an upper limit on the emission through an artificial star experiment if there is no detection. For two CCSNe we detect a point source within the uncertainty region of the SN position. Both have a low chance alignment probability with bright objects within their host galaxies and are thus likely related to the SNe. Comparing the absolute UV magnitude of AT2018cow to the absolute UV magnitudes of the two potential SN detections, there is no evidence that a late-time UV detection of AT2018cow is atypical for interacting SNe. However, when limiting to CCSNe closer than AT2018cow, we see that it is brighter than the upper limits on most non-detections. Combined with a very small late time photospheric radius of AT2018cow, this leads us to conclude that AT2018cow's late-time UV detection was not driven by interaction. It suggests instead that we are possibly viewing the inner region of the explosion. Such properties are naturally expected in tidal disruption models and are less straightforward in supernova scenarios., Comment: 15 pages, 6 figures, 2 tables, accepted for publication in A&A

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.

Aims. The locations of binary neutron star (BNS) mergers within their host galaxies encode the systemic kicks that these systems received in the supernova aftermath. We investigate how the galactic potential and the systemic kicks shape the offset distribution of BNS mergers with a case study of GW 170817 and its host NGC 4993. Methods. We derived dynamical constraints on the host potential from integral field spectroscopy with Jeans anisotropic modelling. We evolved the trajectories of synthetic BNSs from the BPASS code in the galactic potential, using two different kick prescriptions to investigate how the observed offsets might differentiate between these two possibilities. The simulation was repeated after swapping the host potential with that of a dwarf galaxy, to test the effect of the potential on the offsets. Results. The location of GW 170817 is entirely consistent with our predictions, regardless of large or small kicks, because the strong potential of NGC 4993 is only diagnostic of very large kicks. In galaxies of similar or greater mass, large offsets can constrain large kicks, while small offsets do not provide much information. In an old dwarf galaxy, on the other hand, small offsets can constrain small kicks, while large offsets would prevent host association., Comment: 11 pages, 7 figures, accepted for publication in A&A

The nearby long gamma-ray burst (GRB) 190829A was observed using the HST/WFC3/IR grisms about four weeks to 500 days after the burst. We find the spectral features of its associated supernova, SN 2019oyw, are redshifted by several thousands km/s compared to the redshift of the large spiral galaxy on which it is superposed. This velocity offset is seen in several features but most clearly in Ca II NIR triplet $\lambda\lambda$ 8498, 8542, 8662 (CaIR3). We also analyze VLT/FORS and X-shooter spectra of the SN and find strong evolution with time of its P-Cygni features of CaIR3 from the blue to the red. However, comparison with a large sample of Type Ic-BL and Ic SNe shows no other object with the CaIR3 line as red as that of SN 2019oyw were it at the z = 0.0785 redshift of the disk galaxy. This implies that SN 2019oyw is either a highly unusual SN or is moving rapidly with respect to its apparent host. Indeed, using CaIR3 we find the redshift of SN 2019oyw is 0.0944 <= z <= 0.1156. The GRB-SN is superposed on a particularly dusty region of the massive spiral galaxy; therefore, while we see no sign of a small host galaxy behind the spiral, it could be obscured. Our work provides a surprising result on the origins of GRB 190829A, as well as insights into the time evolution of GRB-SNe spectra and a method for directly determining the redshift of a GRB-SN using the evolution of strong spectral features such as CaIR3., Comment: 21 pages, 16 figures

Detections of fast X-ray transients (FXTs) have been accrued over the last few decades. However, their origin has remained mysterious. There is now rapid progress thanks to timely discoveries and localisations with the Einstein Probe mission. Early results indicate that FXTs may frequently, but not always, be associated with gamma-ray bursts (GRBs). Here, we report on the multi-wavelength counterpart of FXT EP240414a, which has no reported gamma-ray counterpart. The transient is located 25.7~kpc in projection from a massive galaxy at $z=0.40$. We perform comprehensive photometric and spectroscopic follow-up. The optical light curve shows at least three distinct emission episodes with timescales of $\sim 1, 4$ and 15 days and peak absolute magnitudes of $M_R \sim -20$, $-21$, and $-19.5$, respectively. The optical spectrum at early times is extremely blue, inconsistent with afterglow emission. It may arise from the interaction of both jet and supernova shock waves with the stellar envelope and a dense circumstellar medium, as has been suggested for some Fast Blue Optical Transients (LFBOTs). At late times, the spectrum evolves to a broad-lined~Type~Ic supernova, similar to those seen in collapsar long-GRBs. This implies that the progenitor of EP240414a is a massive star creating a jet-forming supernova inside a dense envelope, resulting in an X-ray outburst with a luminosity of $\sim 10^{48}$ erg s$^{-1}$, and the complex observed optical/IR light curves. If correct, this argues for a causal link between the progenitors of long-GRBs, FXTs and LFBOTs., Comment: 36 pages, 13 figures, submitted to ApJ

The systems creating binary neutron stars (BNSs) experience systemic kicks when one of the components goes supernova. The combined magnitude of these kicks is still a topic of debate, and has implications for the eventual location of the transient resulting from the merger of the binary. For example, the offsets of short-duration gamma-ray bursts (SGRBs) resulting from BNS mergers depend on the BNS kicks. We investigated Galactic BNSs, and traced their motion through the Galaxy. This enabled us to estimate their kinematic ages and construct a BNS kick distribution, based on their Galactic trajectories. We used the pulsar periods and their derivatives to estimate the characteristic spin-down ages of the binaries. Moreover, we used a Monte Carlo estimation of their present-day velocity vector in order to trace back their trajectory and estimate their kinematic ages. These trajectories, in turn, were used to determine the eccentricity of their Galactic orbit. Based on simulations of kicked objects in the Galactic potential, we investigated the relationship between this eccentricity and kick velocity, in order to constrain the kicks imparted to the binaries at birth. We find that the Galactic BNSs are likely older than $\sim40$ Myr, which means their current (scalar) galactocentric speeds are not representative of their initial kicks. However, we find a close relationship between the eccentricity of a Galactic trajectory and the experienced kick. Using this relation, we constrained the kicks of the Galactic BNSs, depending on the kind of isotropy assumed in estimating their velocity vectors. These kick velocities are well-described by a log-normal distribution peaking around $\sim40-50$ km/s, and coincide with the peculiar velocities of the binaries at their last disc crossing., Comment: 22 pages, 15 figures, accepted for publication by A&A

Fast X-ray transients (FXTs) are X-ray flares lasting minutes to hours. Multi-wavelength counterparts to these FXTs have been proven hard to find. As a result distance measurements are through indirect methods such as host galaxy identification. The three main models proposed for FXTs; supernova shock breakout emission (SN SBO), binary neutron star (BNS) mergers and tidal dirsuption events (TDEs) of an intermediate mass black hole (IMBH) disrupting a white dwarf (WD), have a different associated L$_{X, peak}$. Therefore obtaining the distance to FXTs will be a powerful probe to investigate the nature of these FXTs. We use VLT/MUSE observations of a sample of FXTs to report the redshift of between 13 and 22 galaxies per FXT and use these redshifts to calculate the distance, L$_{X, peak}$ and the projected offsets. We find L$_{X, peak}>10^{44}$ erg s$^{-1}$ if we assume any of the sources with a redshift measurement is the true host galaxy of the corresponding FXT. For XRT 100831 we find a very faint galaxy within the 1$\sigma$ uncertainty region with a chance alignment probability of 0.04. For XRT 060207 we find a candidate host galaxy at z = 0.939 with a low chance alignment probability. However, we also report the detection of a late-type star within the 3$\sigma$ uncertainty region with a similar chance alignment probability. For the remaining FXTs we find no sources within their 3$\sigma$ uncertainty regions. We rule out a SN SBO nature for all FXTs based on L$_{X, peak}$ and the projected offsets. For XRT 100831 we conclude the detected galaxy within the 1$\sigma$ uncertainty position is likely to be the host galaxy of this FXT. From the available information, we are not able to determine if XRT 060207 originated from the galaxy found within 1$\sigma$ of the FXT position or was due to a flare from the late-type star detected within the 3$\sigma$ uncertainty region., Comment: 32 pages (of which 20 appendix), 9 figures, 3 tables. Accepted for publication in A&A

Extragalactic fast X-ray transients (FXTs) are short-duration ($\sim$ ks) X-ray flashes of unknown origin, potentially arising from binary neutron star mergers. We investigate the possible link between FXTs and the afterglows of off-axis merger-induced gamma-ray bursts (GRBs). By modelling the broadband afterglows of 13 GRBs, we make predictions for their X-ray light curve behaviour had they been observed off-axis, considering both a uniform jet with core angle $\theta_{C}$ and a Gaussian-structured jet with truncation angle $\theta_{W} = 2\theta_{C}$. We compare their peak X-ray luminosity, duration, and temporal indices with those of the currently known extragalactic FXTs. Our analysis reveals that a slightly off-axis observing angle of $\theta_{\text{obs}}\approx (2.2-3)\theta_{C}$ and a structured jet are required to explain the shallow temporal indices of the FXT light curves. In the case of a structured jet, the durations of the FXTs are consistent with those of the off-axis afterglows for the same range of observing angles, $\theta_{\text{obs}}\approx (2.2-3)\theta_{C}$. While the off-axis peak X-ray luminosities are consistent only for $\theta_{\text{obs}} = 2.2\theta_{C}$, focussing on individual events reveals that the match of all three properties of the FXTs at the same viewing angle is possible in the range $\theta_{\text{obs}} \sim (2.2-2.6)\theta_{C}$. Despite the small sample of GRBs analysed, these results show that there is a region of the parameter space - although quite limited - where the observational properties of off-axis GRB afterglows can be consistent with those of the newly discovered FXTs. Future observations of FXTs discovered by the recently launched Einstein Probe mission and GRB population studies combined with more complex afterglow models will shed light on this possible GRB-FXT connection, and eventually unveil the progenitors of some FXTs., Comment: 24 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

Long GRB hosts at z<1 are usually low-mass, low metallicity star-forming galaxies. Here we present the until now most detailed, spatially resolved study of the host of GRB 171205A, a grand-design barred spiral galaxy at z=0.036. Our analysis includes MUSE integral field spectroscopy, complemented by high spatial resolution UV/VIS HST imaging and CO(1-0) and HI 21cm data. The GRB is located in a small star-forming region in a spiral arm of the galaxy at a deprojected distance of ~ 8 kpc from the center. The galaxy shows a smooth negative metallicity gradient and the metallicity at the GRB site is half solar, slightly below the mean metallicity at the corresponding distance from the center. Star formation in this galaxy is concentrated in a few HII regions between 5-7 kpc from the center and at the end of the bar, inwards of the GRB region, however, the HII region hosting the GRB is in the top 10% of regions with highest specific star-formation rate. The stellar population at the GRB site has a very young component (< 5 Myr) contributing a significant part of the light. Ionized and molecular gas show only minor deviations at the end of the bar. A parallel study found an asymmetric HI distribution and some additional gas near the position of the GRB, which might explain the star-forming region of the GRB site. Our study shows that long GRBs can occur in many types of star-forming galaxies, however, the actual GRB sites consistently have low metallicity, high star formation and a young population. Furthermore, gas inflow or interactions triggering the star formation producing the GRB progenitor might not be evident in ionized or even molecular gas but only in HI., Comment: 22 pages, 18 figures, Appendix with additional figures, A&A under revision

Luminous Fast Blue Optical Transients (LFBOTs) are a class of extragalactic transients notable for their rapid rise and fade times, blue colour and accompanying luminous X-ray and radio emission. Only a handful have been studied in detail since the prototypical example AT2018cow. Their origins are currently unknown, but ongoing observations of previous and new events are placing ever stronger constraints on their progenitors. We aim to put further constraints on the LFBOT AT2023fhn, and LFBOTs as a class, using information from the multi-wavelength transient light-curve, its host galaxy and local environment. Our primary results are obtained by fitting galaxy models to the spectral energy distribution of AT2023fhn's host and local environment, and by modelling the radio light-curve of AT2023fhn as due to synchrotron self-absorbed emission from an expanding blast-wave in the circumstellar medium. We find that neither the host galaxy nor circumstellar environment of AT2023fhn are unusual compared with previous LFBOTs, but that AT2023fhn has a much lower X-ray to ultraviolet luminosity ratio than previous events. We argue that the variety in ultraviolet-optical to X-ray luminosity ratios among LFBOTs is likely due to viewing angle differences, and that the diffuse, yet young local environment of AT2023fhn - combined with a similar circumstellar medium to previous events - favours a progenitor system containing a massive star with strong winds. Plausible progenitor models in this interpretation therefore include black hole/Wolf-Rayet mergers or failed supernovae., Comment: Accepted for publication in A&A. 15 pages, 12 figures, 8 tables

Various stellar objects experience a velocity kick at some point in their evolution. These include neutron stars and black holes at their birth or binary systems when one of the two components goes supernova. For most of these objects, the magnitude of the kick and its impact on the object dynamics remains a topic of debate. We investigate how kicks alter the velocity distribution of objects born in the Milky Way disc, both immediately after the kick and at later times, and whether these kicks are encoded in the observed population of Galactic neutron stars. We simulate the Galactic trajectories of point masses on circular orbits in the disc after being perturbed by an isotropic kick, with a Maxwellian distribution of magnitudes with $\sigma=265$ km/s. Then, we simulate the motion of these point masses for $200$ Myr. These trajectories are then evaluated, either for the Milky Way population as a whole or for those passing within two kiloparsecs of the Sun, to get the time evolution of the velocities. During the first $20$ Myr, the bulk velocity of kicked objects becomes temporarily aligned to the cylindrical radius, implying an anisotropy in the velocity orientations. Beyond this age, the velocity distribution shifts toward lower values and settles to a median of $\sim200$ km/s. Around the Sun, the distribution also loses its upper tail, primarily due to unbound objects escaping the Galaxy. We compare this to the velocities of Galactic pulsars and find that pulsars show a similar evolution with characteristic age. The shift of the velocity distribution is due to bound objects spending most of their orbits at larger radii after the kick. They are, therefore, decelerated by the Galactic potential. We find the same deceleration to be predicted for nearby objects and the total population and conclude it is also observed in Galactic pulsars., Comment: 15 pages, 12 figures, accepted for publication by A&A

The tidal disruption event (TDE) AT2018fyk has unusual X-ray, UV, and optical light curves that decay over the first $\sim$600d, rebrighten, and decay again around 1200d. We explain this behavior as a one-off TDE associated with a massive black hole (BH) \emph{binary}. The sharp drop-offs from $t^{-5/3}$ power laws at around 600d naturally arise when one BH interrupts the debris fallback onto the other BH. The BH mass $M_\bullet$ derived from fitting X-ray spectra with a slim disk accretion model and, independently, from fitting the early UV/optical light curves, is smaller by two orders of magnitude than predicted from the $M_\bullet$--$\sigma_*$ host galaxy relation, suggesting that the debris is accreted onto the secondary, with fallback cut off by the primary. Furthermore, if the rebrightening were associated with the primary, it should occur around 5000d, not the observed 1200d. The secondary's mass and dimensionless spin is $M_{\bullet,{\rm s}}=2.7^{+0.5}_{-1.5} \times 10^5 M_\odot$ and $a_{\bullet,{\rm s}}>0.3$ (X-ray spectral fitting), while the primary's mass is $M_{\bullet,{\rm p}}=10^{7.7\pm0.4}M_\odot$ ($M_\bullet$-$\sigma_*$ relation). An intermediate mass BH secondary is consistent with the observed UV/optical light curve decay, i.e., the secondary's outer accretion disk is too faint to produce a detectable emission floor. The time of the first accretion cutoff constrains the binary separation to be $(6.7\pm 1.2) \times 10^{-3}~{\rm pc}$. X-ray spectral fitting and timing analysis indicate that the hard X-rays arise from a corona above the secondary's disk. The early UV/optical emission, suggesting a super-Eddington phase for the secondary, possibly originates from shocks arising from debris circularization., Comment: 5 tables, 7 figures

Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a, whose bright peak was also detected by the Swift Burst Alert Telescope and Konus-Wind through off-line analyses. At a redshift of $z=4.859$, EP240315a showed a much longer and more complicated light curve in the soft X-ray band than in gamma-rays. Benefiting from a large field-of-view ($\sim$3600 deg$^2$) and a high sensitivity, EP-WXT captured the earlier engine activation and extended late engine activity through a continuous detection. With a peak X-ray flux at the faint end of previously known high-$z$ GRBs, the detection of EP240315a demonstrates the great potential for EP to study the early universe via GRBs., Comment: 41 pages, 8 figures, 7 tables

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

Gamma-ray bursts (GRBs) are ideal probes of the Universe at high redshift (z > 5), pinpointing the locations of the earliest star-forming galaxies and providing bright backlights that can be used to spectrally fingerprint the intergalactic medium and host galaxy during the period of reionization. Future missions such as Gamow Explorer are being proposed to unlock this potential by increasing the rate of identification of high-z GRBs to rapidly trigger observations from 6-10 m ground telescopes, JWST, and the Extremely Large Telescopes. Gamow was proposed to the NASA 2021 Medium-Class Explorer (MIDEX) program as a fast-slewing satellite featuring a wide-field lobster-eye X-ray telescope (LEXT) to detect and localize GRBs, and a 30 cm narrow-field multi-channel photo-z infrared telescope (PIRT) to measure their photometric redshifts using the Lyman-alpha dropout technique. To derive the PIRT sensitivity requirement we compiled a complete sample of GRB optical-near-infrared afterglows from 2008 to 2021, adding a total of 66 new afterglows to our earlier sample, including all known high-z GRB afterglows. We performed full light-curve and spectral-energy-distribution analyses of these afterglows to derive their true luminosity at very early times. For all the light curves, where possible, we determined the brightness at the time of the initial finding chart of Gamow, at different high redshifts and in different NIR bands. We then followed the evolution of the luminosity to predict requirements for ground and space-based follow-up. We find that a PIRT sensitivity of 15 micro-Jy (21 mag AB) in a 500 s exposure simultaneously in five NIR bands within 1000s of the GRB trigger will meet the Gamow mission requirement to recover > 80% of all redshifts at z > 5., Comment: 44 pages, 10 figures, 1 table. Accepted for publication in Astronomy and Astrophysics 15 Feb 2024. Abstract abridged for arXiv

Twenty-two extragalactic fast X-ray transients (FXTs) have now been discovered from two decades of Chandra data (analyzing ~259 Ms of data), with 17 associated with distant galaxies (>100 Mpc). Different mechanisms and progenitors have been proposed to explain their properties; nevertheless, after analyzing their timing, spectral parameters, host-galaxy properties, luminosity function, and volumetric rates, their nature remains uncertain. We interpret a sub-sample of nine FXTs that show a plateau or a fast-rise light curve within the framework of a binary neutron star (BNS) merger magnetar model. We fit their light curves and derive magnetar (magnetic field and initial rotational period) and ejecta (ejecta mass and opacity) parameters. This model predicts two zones: an orientation-dependent free zone (where the magnetar spin-down X-ray photons escape freely to the observer) and a trapped zone (where the X-ray photons are initially obscured and only escape freely once the ejecta material becomes optically thin). We argue that six FXTs show properties consistent with the free zone and three FXTs with the trapped zone. This sub-sample of FXTs has a similar distribution of magnetic fields and initial rotation periods to those inferred for short gamma-ray bursts (SGRBs), suggesting a possible association. We compare the predicted ejecta emission fed by the magnetar emission (called merger-nova) to the optical and near-infrared upper limits of two FXTs, XRT 141001 and XRT 210423 where contemporaneous optical observations are available. The non-detections place lower limits on the redshifts of XRT 141001 and XRT 210423 of z>1.5 and >0.1, respectively. If the magnetar remnants lose energy via gravitational waves, it should be possible to detect similar objects with the current advanced LIGO detectors out to a redshift z<0.03, while future GW detectors will be able to detect them out to z=0.5., Comment: The paper was accepted for publication in Astronomy & Astrophysics

Extragalactic fast X-ray transients (FXTs) are a class of soft (0.3-10 keV) X-ray transients lasting a few hundred seconds to several hours. Several progenitor mechanisms have been suggested to produce FXTs, including supernova shock breakouts, binary neutron star mergers, or tidal disruptions involving an intermediate-mass black hole and a white dwarf. We present detailed host studies, including spectroscopic observations of the host galaxies of 7 XMM-Newton-discovered FXTs. The candidate hosts lie at redshifts 0.0928 $< z <$ 0.645 implying peak X-ray luminosities of 10$^{43}$ erg s$^{-1}$ $< L_X <$ 10$^{45}$ erg s$^{-1}$,and physical offsets of 1 kpc < $r_\mathrm{proj}$ < 22 kpc. These observations increase the number of FXTs with a spectroscopic redshift measurement by a factor of 2, although we note that one event is re-identified as a Galactic flare star. We infer host star formation rates and stellar masses by fitting the combined spectroscopic and archival photometric data. We also report on a contemporaneous optical counterpart search to the FXTs in Pan-STARRS and ATLAS by performing forced photometry at the position of the FXTs. We do not find any counterpart in our search. Given our constraints, including peak X-ray luminosities, optical limits, and host properties, we find that XRT 110621 is consistent with a SN SBO event. Spectroscopic redshifts of likely host galaxies for four events imply peak X-ray luminosities that are too high to be consistent with SN SBOs, but we are unable to discard either the BNS or WD-IMBH TDE scenarios for these FXTs., Comment: 17 pages, 10 figures, 6 tables. Accepted for publication in MNRAS

The existence of a secondary (in addition to compact object mergers) source of heavy element ($r$-process) nucleosynthesis, the core-collapse of rapidly-rotating and highly-magnetized massive stars, has been suggested by both simulations and indirect observational evidence. Here, we probe a predicted signature of $r$-process enrichment, a late-time ($\gtrsim 40$ days post-burst) distinct red color, in observations of GRB-supernovae (GRB-SNe) which are linked to these massive star progenitors. We present optical to near-IR color measurements of four GRB-SNe at $z \lesssim 0.4$, extending out to $> 500$ days post-burst, obtained with the Hubble Space Telescope and large-aperture ground-based telescopes. Comparison of our observations to models indicates that GRBs 030329, 100316D and 130427A are consistent with both no enrichment and producing $0.01 - 0.15 M_{\odot}$ of $r$-process material if there is a low amount of mixing between the inner $r$-process ejecta and outer SN layers. GRB 190829A is not consistent with any models with $r$-process enrichment $\geq 0.01 M_{\odot}$. Taken together the sample of GRB-SNe indicates color diversity at late times. Our derived yields from GRB-SNe may be underestimated due to $r$-process material hidden in the SN ejecta (potentially due to low mixing fractions) or the limits of current models in measuring $r$-process mass. We conclude with recommendations for future search strategies to observe and probe the full distribution of $r$-process produced by GRB-SNe., Comment: Resubmission after comments. Accepted to ApJ. 36 pages, 7 figures, 2 tables

The existence of merging black hole-neutron star (BHNS) binaries has been ascertained through the observation of their gravitational wave (GW) signals. However, to date, no definitive electromagnetic (EM) emission has been confidently associated with these mergers. Such an association could help unravel crucial information on these systems, for example, their BH spin distribution, the equation of state (EoS) of NS and the rate of heavy element production. We model the multi-messenger (MM) emission from BHNS mergers detectable during the fourth (O4) and fifth (O5) observing runs of the LIGO-Virgo-KAGRA GW detector network, in order to provide detailed predictions that can help enhance the effectiveness of observational efforts and extract the highest possible scientific information from such remarkable events. Our methodology is based on a population synthesis-approach, which includes the modelling of the signal-to-noise ratio of the GW signal in the detectors, the GW-inferred sky localization of the source, the kilonova (KN) optical and near-infrared light curves, the relativistic jet gamma-ray burst (GRB) prompt emission peak photon flux, and the GRB afterglow light curves in the radio, optical and X-ray bands. The resulting prospects for BHNS MM detections during O4 are not promising, with a GW detection rate of $15.0^{+15.4}_{-8.8}$ yr$^{-1}$, but joint MM rates of $\sim 10^{-1}$ yr$^{-1}$ for the KN and $\sim 10^{-2}$ yr$^{-1}$ for the jet-related emission. In O5 we find an overall increase in expected detection rates by around an order of magnitude, owing to both the enhanced sensitivity of the GW detector network, and the coming online of future EM facilities. Finally, we discuss direct searches for the GRB radio afterglow with large-field-of-view instruments as a new possible follow-up strategy in the context of ever-dimming prospects for KN detection., Comment: Submitted to A&A. 17 pages, 11 figures, 2 tables. Comments are welcome!

The merger locations of binary neutron stars (BNSs) encode their galactic kinematics and provide insights into their connection to short gamma-ray bursts (SGRBs). In this work, we use the sample of Galactic BNSs with measured proper motions to investigate their kinematics and predict their merger locations. Using a synthetic image of the Milky Way and its Galactic potential we analyse the BNS mergers as seen from an extragalactic viewpoint and compare them to the location of SGRBs on and around their host galaxies. We find that the Galactocentric transverse velocities of the BNSs are similar in magnitude and direction to those of their Local Standards of Rest, which implies that the present-day systemic velocities are not isotropically oriented and the peculiar velocities might be as low as those of BNS progenitors. Both systemic and peculiar velocities fit a lognormal distribution, with the peculiar velocities being as low as $\sim 22-157$ km s$^{-1}$. We also find that the observed BNS sample is not representative of the whole Galactic population, but rather of systems born around the Sun's location with small peculiar velocities. When comparing the predicted BNS merger locations to SGRBs, we find that they cover the same range of projected offsets, host-normalized offsets, and fractional light. Therefore, the spread in SGRB locations can be reproduced by mergers of BNSs born in the Galactic disk with small peculiar velocities, although the median offset match is likely a coincidence due to the biased BNS sample., Comment: 13 pages, 7 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

The bright, blue, rapidly evolving AT2018cow is a well-studied peculiar extragalactic transient. Despite an abundance of multi-wavelength data, there still is no consensus on the nature of the event. We present our analysis of three epochs of Hubble Space Telescope (HST) observations spanning the period from 713-1474 days post burst, paying particular attention to uncertainties of the transient photometry introduced by the complex background in which AT2018cow resides. Photometric measurements show evident fading in the UV and more subtle but significant fading in the optical. During the last HST observation, the transient's optical/UV colours were still bluer than those of the substantial population of compact, young, star-forming regions in the host of AT2018cow, suggesting some continued transient contribution to the light. However, a compact source underlying the transient would substantially modify the resulting spectral energy distribution, depending on its contribution in the various bands. In particular, in the optical filters, the complex, diffuse background poses a problem for precise photometry. An underlying cluster is expected for a supernova occurring within a young stellar environment or a tidal-disruption event (TDE) within a dense older one. While many recent works have focused on the supernova interpretation, we note the substantial similarity in UV light-curve morphology between AT2018cow and several tidal disruption events around supermassive black holes. Assuming AT2018cow arises from a TDE-like event, we fit the late-time emission with a disc model and find $M_{BH} = 10^{3.2{\pm}0.8}$ M$_{\odot}$. Further observations are necessary to determine the late-time evolution of the transient and its immediate environment., Comment: 16 pages, 10 figures. Accepted for publication in MNRAS

The Interplanetary Network (IPN) is a detection, localization and alert system that utilizes the arrival time of transient signals in gamma-ray detectors on spacecraft separated by planetary baselines to geometrically locate the origin of these transients. Due to the changing astrophysical landscape and the new emphasis on time domain and multi-messenger astrophysics (TDAMM) from the Pathways to Discovery in Astronomy and Astrophysics for the 2020s, this Gamma-ray Transient Network Science Analysis Group was tasked to understand the role of the IPN and high-energy monitors in this new era. The charge includes describing the science made possible with these facilities, tracing the corresponding requirements and capabilities, and highlighting where improved operations of existing instruments and the IPN would enhance TDAMM science. While this study considers the full multiwavelength and multimessenger context, the findings are specific to space-based high-energy monitors. These facilities are important both for full characterization of these transients as well as facilitating follow-up observations through discovery and localization. The full document reports a brief history of this field, followed by our detailed analyses and findings in some 68 pages, providing a holistic overview of the role of the IPN and high-energy monitors in the coming decades., Comment: Terms of Reference and additional information on the Science Analysis Group are available at https://pcos.gsfc.nasa.gov/sags/gtn-sag.php

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

Extragalactic fast X-ray transients (FXTs) are short flashes of X-ray photons of unknown origin that last a few minutes to hours. We extend the search for extragalactic FXTs from Quirola et al. 2022 (Paper I; based on sources in the Chandra Source Catalog 2.0, CSC2) to further Chandra archival data between 2014-2022. We extract X-ray data using a method similar to that employed by CSC2 and apply identical search criteria as in Paper I. We report the detection of eight FXT candidates, with peak 0.3-10 keV fluxes between 1$\times$10$^{-13}$ to 1$\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ and $T_{90}$ values from 0.3 to 12.1 ks. This sample of FXTs has likely redshifts between 0.7 to 1.8. Three FXT candidates exhibit light curves with a plateau (${\approx}$1-3 ks duration) followed by a power-law decay and X-ray spectral softening, similar to what was observed for a few previously reported FXTs in Paper I. In light of the new, expanded source lists (eight FXTs with known redshifts from Paper I and this work), we update the event sky rates derived in Paper I, finding 36.9$_{-8.3}^{+9.7}$ deg$^{-2}$ yr$^{-1}$ for the extragalactic samples for a limiting flux of ${\gtrsim}$1${\times}$10$^{-13}$ erg cm$^{-2}$ s$^{-1}$, calculate the first FXT X-ray luminosity function, and compare the volumetric density rate between FXTs and other transient classes. Our latest Chandra-detected extragalactic FXT candidates boost the total Chandra sample by $\sim$50 %, and appear to have a similar diversity of possible progenitors., Comment: 37 pages, 19 figures, 11 tables, accepted for publication in A&A. arXiv admin note: text overlap with arXiv:2201.07773

The majority of massive stars are born in binaries, and most unbind upon the first supernova. With precise proper motion surveys such as Gaia, it is possible to trace back the motion of stars in the vicinity of young remnants to search for ejected companions. Establishing the fraction of remnants with an ejected companion, and the photometric and kinematic properties of these stars, offers unique insight into supernova progenitor systems. In this paper, we employ binary population synthesis to produce kinematic and photometric predictions for ejected secondary stars. We demonstrate that the unbound neutron star velocity distribution from supernovae in binaries closely traces the input kicks. Therefore, the observed distribution of neutron star velocities should be representative of their natal kicks. We evaluate the probability for any given filter, magnitude limit, minimum measurable proper motion (as a function of magnitude), temporal baseline, distance and extinction that an unbound companion can be associated with a remnant. We compare our predictions with results from previous companion searches, and demonstrate that the current sample of stars ejected by the supernova of their companion can be increased by a factor of 5-10 with Gaia data release 3. Further progress in this area is achievable by leveraging the absolute astrometric precision of Gaia, and by obtaining multiple epochs of deep, high resolution near-infrared imaging with the Hubble Space Telescope, JWST and next-generation wide-field near-infrared observatories such as Euclid or the Nancy Grace Roman Space Telescope., Comment: Accepted for publication in MNRAS. 19 pages, 17 figures

The highly variable and energetic pulsed emission of a long gamma-ray burst (GRB) is thought to originate from local, rapid dissipation of kinetic or magnetic energy within an ultra-relativistic jet launched by a newborn compact object, formed during the collapse of a massive star. The spectra of GRB pulses are best modelled by power-law segments, indicating the dominance of non-thermal radiation processes. Spectral lines in the X-ray and soft $\gamma$-ray regime for the afterglow have been searched for intensively, but never confirmed. No line features ever been identified in the high energy prompt emission. Here we report the discovery of a highly significant ($> 6 \sigma$) narrow emission feature at around $10$ MeV in the brightest ever GRB 221009A. By modelling its profile with a Gaussian, we find a roughly constant width $\sigma \sim 1$ MeV and temporal evolution both in energy ($\sim 12$ MeV to $\sim 6$ MeV) and luminosity ($\sim 10^{50}$ erg/s to $\sim 2 \times 10^{49}$ erg/s) over 80 seconds. We interpret this feature as a blue-shifted annihilation line of relatively cold ($k_\mathrm{B}T\ll m_\mathrm{e}c^2$) electron-positron pairs, which could have formed within the jet region where the brightest pulses of the GRB were produced. A detailed understanding of the conditions that can give rise to such a feature could shed light on the so far poorly understood GRB jet properties and energy dissipation mechanism., Comment: Submitted

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

Fast X-ray Transients (FXTs) are X-ray flares with a duration ranging from a few hundred seconds to a few hours. Possible origins include the tidal disruption of a white dwarf by an intermediate-mass black hole, a supernova shock breakout, and a binary neutron star merger. We present the X-ray light curve and spectrum, and deep optical imaging of the FXT XRT 210423, which has been suggested to be powered by a magnetar produced in a binary neutron star merger. Our Very Large Telescope and Gran Telescopio Canarias (GTC) observations began on May 6, 2021, thirteen days after the onset of the flare. No transient optical counterpart is found in the 1" (3$\sigma$) X-ray uncertainty region of the source to a depth $g_{s}$=27.0 AB mag. A candidate host lies within the 1" X-ray uncertainty region with a magnitude of 25.9 $\pm$ 0.1 in the GTC/HiPERCAM $g_s$-filter. Due to its faintness, it was not detected in other bands, precluding a photometric redshift determination. We detect two additional candidate host galaxies; one with $z_{\rm spec}=1.5082 \pm 0.0001$ and an offset of 4.2$\pm$1" (37$\pm$9 kpc) from the FXT and another one with $z_{\rm phot}=1.04^{+0.22}_{-0.14}$, at an offset of 3.6$\pm$1" (30$\pm$8 kpc). Based on the properties of all the prospective hosts we favour a binary neutron star merger, as previously suggested in the literature, as explanation for XRT 210423., Comment: 14 pages, 6 figures. Accepted for publication in ApJ

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, submitted to Astronomy & Astrophysics

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

Binary neutron star mergers are one of the ultimate events of massive binary star evolution, and our understanding of their parent system is still in its infancy. Upcoming gravitational wave detections, coupled with multi-wavelength follow-up observations, will allow us to study an increasing number of these events by characterising their neighbouring stellar populations and searching for their progenitors. Stellar evolution simulations are essential to this work but they are also based on numerous assumptions. Additionally, the models used to study the host galaxies differ from those used to characterise the progenitors and are typically based on single star populations. Here we introduce a framework to perform an end-to-end analysis and deploy it to the first binary neutron star merger - GW170817. With the Binary Population And Spectral Synthesis (BPASS) codes we are able to retrieve the physical properties of the host galaxy NGC 4993 as well as infer progenitor candidates. In our simulations there is a >98% chance that GW170817 originated from a stellar population with Z=0.010 born between 5 and 12.5 Gyrs ago. By carefully weighing the stellar genealogies we find that GW170817 most likely came from a binary system born with a 13-24 Msol primary and 10-12 Msol secondary which underwent two or three common envelope events over their lifetime., Comment: Main text (3 Figures, 10 pages) and Supplementary Information (23 pages, 14 Figures). Published in Nature Astronomy

The study of the properties of galaxies in the first billion years after the Big Bang is one of the major topic of current astrophysics. Optical/near-infrared spectroscopy of the afterglows of long Gamma-ray bursts (GRBs) provide a powerful diagnostic tool to probe the interstellar medium (ISM) of their host galaxies and foreground absorbers, even up to the highest redshifts. We analyze the VLT/X-shooter afterglow spectrum of GRB 210905A, triggered by the Swift Neil Gehrels Observatory, and detect neutral-hydrogen, low-ionization, high-ionization, and fine-structure absorption lines from a complex system at z=6.3118, that we associate with the GRB host galaxy. We study the ISM properties of the host system, revealing the metallicity, kinematics and chemical abundance pattern. The total metallicity of the z~6.3 system is [M/H]=-1.72+/-0.13, after correcting for dust-depletion and taking into account alpha-element enhancement. In addition, we determine the overall amount of dust and dust-to-metal mass ratio (DTM) ([Zn/Fe]_fit=0.33+/-0.09, DTM=0.18+/-0.03). We find indications of nucleosynthesis due to massive stars and evidence of peculiar over-abundance of aluminium. From the analysis of fine-structure lines, we determine distances of several kpc for the low-ionization gas clouds closest to the GRB. Those farther distances are possibly due to the high number of ionizing photons. Using the HST/F140W image of the GRB field, we show the GRB host galaxy as well as multiple objects within 2" from the GRB. We discuss the galaxy structure and kinematics that could explain our observations, also taking into account a tentative detection of Lyman-alpha emission. Deep spectroscopic observations with VLT/MUSE and JWST will offer the unique possibility of combining our results with the ionized-gas properties, with the goal of better understanding how galaxies in the reionization era form and evolve., Comment: Accepted Publication (In Press on A&A) - 22 pages, 10 figures, 6 tables - Appendix: 6 figures, 3 tables

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

CSS100217 was a nuclear, rapid and luminous flare in a Narrow-Line Seyfert 1 galaxy, whose initial interpretation as a supernova is now debated between variability of the active galactic nucleus (AGN) or a tidal disruption event (TDE). In this paper, we present and discuss new evidence in favour of a TDE or extreme flaring episode scenario. After the decay of the flare, the galaxy entered a long-term low luminosity state, 0.4 magnitudes lower than the pre-outburst emission in the V band. We attribute this to the creation of a cavity in the accretion disk after the tidal disruption of a star in a retrograde orbit with respect to the accretion disk rotation, making a TDE our favoured interpretation of the flare. We also show how the host galaxy shows a point-like, compact profile, no evidence for an extended component and a relatively low mass, unlike what expected from an AGN host galaxy at z=0.147. A compact host galaxy may result in an increased TDE rate, strengthening our interpretation of the event., Comment: 11 pages, 4 figures, accepted for publication on MNRAS

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

It is well established that magnetars are neutron stars with extreme magnetic fields and young ages, but the evolutionary pathways to their creation are still uncertain. Since most massive stars are in binaries, if magnetars are a frequent result of core-collapse supernovae, some fraction are expected to have a bound companion at the time of observation. In this paper, we utilise literature constraints, including deep Hubble Space Telescope imaging, to search for bound stellar companions to magnetars. The magnitude and colour measurements are interpreted in the context of binary population synthesis predictions. We find two candidates for stellar companions associated with CXOU J171405.7-381031 and SGR 0755-2933, based on their J-H colours and H-band absolute magnitudes. Overall, the proportion of the Galactic magnetar population with a plausibly stellar near-infrared counterpart candidate, based on their magnitudes and colours, is between 5 and 10 per cent. This is consistent with a population synthesis prediction of 5 per cent, for the fraction of core-collapse neutron stars arising from primaries which remain bound to their companion after the supernova. These results are therefore consistent with magnetars being drawn in an unbiased way from the natal core-collapse neutron star population, but some contribution from alternative progenitor channels cannot be ruled out., Comment: 14 pages, 6 figures, 3 tables, accepted for publication in MNRAS

We report the discovery of six new magnetar counterpart candidates from deep near-infrared Hubble Space Telescope imaging. The new candidates are among a sample of nineteen magnetars for which we present HST data obtained between 2018-2020. We confirm the variability of previously established near-infrared counterparts, and newly identify candidates for PSRJ1622-4950, SwiftJ1822.3-1606, CXOUJ171405.7-381031, SwiftJ1833-0832, SwiftJ1834.9-0846 and AXJ1818.8-1559 based on their proximity to X-ray localisations. The new candidates are compared with the existing counterpart population in terms of their colours, magnitudes, and near-infrared to X-ray spectral indices. We find two candidates for AXJ1818.8-1559 which are both consistent with previously established counterparts. The other new candidates are likely to be chance alignments, or otherwise have a different origin for their near-infrared emission not previously seen in magnetar counterparts. Further observations and studies of these candidates are needed to firmly establish their nature., Comment: 11 pages, 5 figures, 3 tables, accepted for publication in MNRAS

We present a detailed follow-up of the very energetic GRB 210905A at a high redshift of z = 6.312 and its luminous X-ray and optical afterglow. We obtained a photometric and spectroscopic follow-up in the optical and near-infrared (NIR), covering both the prompt and afterglow emission from a few minutes up to 20 Ms after burst. With an isotropic gamma-ray energy release of Eiso = 1.27E54 erg, GRB 210905A lies in the top ~7% of gamma-ray bursts (GRBs) in terms of energy released. Its afterglow is among the most luminous ever observed. It starts with a shallow evolution that can be explained by energy injection, and it is followed by a steeper decay, while the spectral energy distribution is in agreement with slow cooling in a constant-density environment within the standard fireball theory. A jet break at ~ 46.2+-16.3 d (~6.3 d rest-frame) has been observed in the X-ray light curve; however, it is hidden in the H band due to the contribution from the likely host galaxy, the fourth GRB host at z > 6 known to date. We derived a half-opening angle of 8.4+-1.0 degrees, which is the highest ever measured for a z>6 burst, but within the range covered by closer events. The resulting collimation-corrected gamma-ray energy release of 1E52 erg is also among the highest ever measured. The moderately large half-opening angle argues against recent claims of an inverse dependence of the half-opening angle on the redshift. The total jet energy is likely too large to be sustained by a standard magnetar, and it suggests that the central engine of this burst was a newly formed black hole. Despite the outstanding energetics and luminosity of both GRB 210905A and its afterglow, we demonstrate that they are consistent with those of less distant bursts, indicating that the powering mechanisms and progenitors do not evolve significantly with redshift., Comment: 19 pages, 12 figures, 5 tables, accepted in Astronomy & Astrophysics

Extragalactic Fast X-ray Transients (FXRTs) are short flashes of X-ray photons spanning a few seconds to hours, with an uncertain origin. Our ignorance about their physical mechanisms and progenitor systems is due in part to the lack of clear multi-wavelength counterparts in most cases because they only have been identified serendipitously. We develop a systematic search of FXRTs using a straightforward X-ray flare search algorithm, in the Chandra Source Catalog (Data Release 2.0; 169.6 Ms over 592.4 deg$^{2}$ using only observations with $|b|>10^{\circ}$ and before 2015), incorporating various multi-wavelength constraints to rule out Galactic contamination and characterize the candidates. We report the detection of 14 FXRT candidates from a parent sample of 214,701 sources. Candidates have peak 0.5-7 keV fluxes between 1$\times$10$^{-13}$ to 2$\times$10$^{-10}$ erg cm$^{-2}$ s$^{-1}$ and $T_{90}$ values from 4 to 48 ks. The sample can be subdivided into two groups: six "nearby" FXRTs that occurred within $d\lesssim$100 Mpc and eight "distant" FXRTs with likely redshifts $\gtrsim$0.1. Three distant FXRT candidates exhibit light curves with a plateau (${\approx}$1-3 ks duration) followed by a power-law decay and X-ray spectral softening, similar to what was observed for the previously reported FXRT CDF-S~XT2, a proposed magnetar-powered binary neutron star merger event. After applying completeness corrections, we calculate event rates for the nearby and distant samples of 53.7$_{-15.1}^{+22.6}$ and 28.2$_{-6.9}^{+9.8}$ deg$^{-2}$ yr$^{-1}$, respectively. This novel sample of Chandra-detected extragalactic FXRT candidates, although modest in size, breaks new ground in terms of characterizing the diverse properties, nature, and possible progenitors of these enigmatic events., Comment: 44 pages, 22 figures, submitted to A&A on Jan 6, 2022, and accepted on Apr 11, 2022

We present deep Hubble Space Telescope near-infrared (NIR) observations of the magnetar SGR 1935+2154 from June 2021, approximately 6 years after the first HST observations, a year after the discovery of fast radio burst like emission from the source, and in a period of exceptional high frequency activity. Although not directly taken during a bursting period the counterpart is a factor of ~1.5 to 2.5 brighter than seen at previous epochs with F140W(AB) = $24.65\pm0.02$ mag. We do not detect significant variations of the NIR counterpart within the course of any one orbit (i.e. on minutes-hour timescales), and contemporaneous X-ray observations show SGR 1935+2154 to be at the quiescent level. With a time baseline of 6 years from the first identification of the counter-part we place stringent limits on the proper motion of the source, with a measured proper motion of ${\mu} = 3.1\pm1.5$ mas/yr. The direction of proper motion indicates an origin of SGR 1935+2154 very close to the geometric centre of SNR G57.2+08, further strengthening their association. At an adopted distance of $6.6\pm0.7$ kpc, the corresponding tangential space velocity is ${\nu_T} = 97\pm48$ km/s (corrected for differential Galactic rotation and peculiar Solar motion), although its formal statistical determination may be compromised owing to few epochs of observation. The current velocity estimate places it at the low end of the kick distribution for pulsars, and makes it among the lowest known magnetar kicks. When collating the few-magnetar kick constraints available, we find full consistency between the magnetar kick distribution and the much larger pulsar kick sample, Comment: Accepted for publication in ApJ

We observed GRB190114C (redshift z = 0.4245), the first GRB ever detected at TeV energies, at optical and near-infrared wavelengths with several ground-based telescopes and the Hubble Space Telescope, with the primary goal of studying its underlying supernova, SN2019jrj. The monitoring spanned the time interval between 1.3 and 370 days after the burst, in the observer frame. We find that the afterglow emission can be modelled with a forward shock propagating in a uniform medium modified by time-variable extinction along the line of sight. A jet break could be present after 7 rest-frame days, and accordingly the maximum luminosity of the underlying SN ranges between that of stripped-envelope corecollapse supernovae (SNe) of intermediate luminosity, and that of the luminous GRB-associated SN2013dx. The observed spectral absorption lines of SN2019jrj are not as broad as in classical GRB-SNe, and are rather more similar to those of less-luminous core-collapse SNe. Taking the broad-lined stripped-envelope core-collapse SN2004aw as an analogue, we tentatively derive the basic physical properties of SN2019jrj. We discuss the possibility that a fraction of the TeV emission of this source might have had a hadronic origin and estimate the expected high-energy neutrino detection level with IceCube., Comment: 9 pages, 4 figures, accepted for publication in A&A

The Gamow Explorer will use Gamma Ray Bursts (GRBs) to: 1) probe the high redshift universe (z > 6) when the first stars were born, galaxies formed and Hydrogen was reionized; and 2) enable multi-messenger astrophysics by rapidly identifying Electro-Magnetic (IR/Optical/X-ray) counterparts to Gravitational Wave (GW) events. GRBs have been detected out to z ~ 9 and their afterglows are a bright beacon lasting a few days that can be used to observe the spectral fingerprints of the host galaxy and intergalactic medium to map the period of reionization and early metal enrichment. Gamow Explorer is optimized to quickly identify high-z events to trigger follow-up observations with JWST and large ground-based telescopes. A wide field of view Lobster Eye X-ray Telescope (LEXT) will search for GRBs and locate them with arc-minute precision. When a GRB is detected, the rapidly slewing spacecraft will point the 5 photometric channel Photo-z Infra-Red Telescope (PIRT) to identify high redshift (z > 6) long GRBs within 100s and send an alert within 1000s of the GRB trigger. An L2 orbit provides > 95% observing efficiency with pointing optimized for follow up by the James Webb Space Telescope (JWST) and ground observatories. The predicted Gamow Explorer high-z rate is >10 times that of the Neil Gehrels Swift Observatory. The instrument and mission capabilities also enable rapid identification of short GRBs and their afterglows associated with GW events. The Gamow Explorer will be proposed to the 2021 NASA MIDEX call and if approved, launched in 2028., Comment: 14 pages, 8 Figures

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