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.

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

We study the emergence of prethermal Floquet Time Crystal (pFTC) in disordered chiral multiferroic chains. The model is an extension of the usual periodically driven nearest-neighbor disordered Heisenberg chain, with additional next-nearest-neighbor Heisenberg couplings and DMI interactions due to external magnetic and electric couplings. We derive the phase diagram of the model, characterizing the magnetization, entanglement, and coherence dynamics of the system along the extended interactions. In addition, we explore the application of the pFTC as quantum sensors of AC fields. The sensor performance to estimate small AC fields is quantified through the quantum Fisher information (QFI) measure. The sensor offers several advantages as compared to those composed of non-interacting spins due to its intrinsic robustness, long coherent interrogation time, and many-body correlations. Specifically, the sensor can overcome the standard quantum limit ($\rm{SQL} \sim N t^2$) during the prethermal regime, reaching an optimum performance at the pFTC lifetime $t^*$, where the $\rm{QFI}/Nt^{*^2} \sim N^\alpha$ with $\alpha > 0$, scaling superlinarly with the number of spins. Different from \text{full} FTCs, the prethermal lifetime does not diverge in the thermodynamic limit, nevertheless it can be exponentially long with tuning system parameters., Comment: 9 pages, 5 figures

Intrinsic bulk viscosity of the one-component plasma (OCP) is computed and analyzed using equilibrium molecular dynamics simulations and the Green-Kubo formalism. It is found that bulk viscosity exhibits a maximum at $\Gamma \approx 1$, corresponding to the condition that the average kinetic energy of particles equals the potential energy at the average inter-particle spacing. The weakly coupled and strongly coupled limits are analyzed and used to construct a model that captures the full range of coupling strengths simulated: $\Gamma \approx 10^{-2} - 10^2$. Simulations are also run of the Yukawa one-component plasma (YOCP) in order to understand the impact of electron screening. It is found that electron screening leads to a smaller bulk viscosity due to a reduction in the excess heat capacity of the system. Bulk viscosity is shown to be at least an order of magnitude smaller than shear viscosity in both the OCP and YOCP. The generalized frequency-dependent bulk viscosity coefficient is also analyzed. This is found to exhibit a peak near twice the plasma frequency in strongly coupled conditions, which is associated with the oscillatory decay observed in the bulk viscosity autocorrelation function. The generalized shear and bulk viscosity coefficients are found to have a similar magnitude for $\omega \gtrsim 2\omega_p$ at strongly coupled conditions.

CDF-S XT1 and XT2 are considered two "canonical" extragalactic fast X-ray transients (FXTs). In this work, we report new constraints on both FXTs, based on recent JWST NIRCam and MIRI photometry, as well as NIRspec spectroscopy for CDF-S XT2 that allow us to improve our understanding of their distances, energetics, and host galaxy properties compared to the pre-JWST era. We use the available HST and JWST archival data to determine the host properties and constrain the energetics of each FXT based on spectral energy distribution (SED) photometric fitting. The host of CDF-S XT1 is now constrained to lie at $z_{phot}{=}2.76_{-0.13}^{+0.21}$, implying a host absolute magnitude $M_{R}=-19.14$ mag, stellar mass $M_{*}$=2.8e8 $M_\odot$, and star formation rate SFR=0.62 $M_\odot$~yr$^{-1}$. These properties lie at the upper end of previous estimates, leaving CDF-S XT1 with a peak X-ray luminosity of $L_{X,peak}$=2.8e47 erg s$^{-1}$. We argue that the best progenitor scenario for XT1 is a low-luminosity gamma-ray burst (GRB), although we do not fully rule out a proto-magnetar association or a jetted tidal disruption event involving a white dwarf and an intermediate-massive black hole. In the case of CDF-S XT2, JWST imaging reveals a new highly obscured component of the host galaxy, previously missed by HST, while NIRspec spectroscopy securely places the host at $z_{spec}{=}3.4598{\pm}0.0022$. The new redshift implies a host with $M_{R}=-21.76$ mag, $M_{*}$=5.5e10 $M_\odot$, SFR=160 $M_\odot$ yr$^{-1}$, and FXT $L_{X,peak}$=1.4e47 erg s$^{-1}$. The revised energetics, similarity to X-ray flash event light curves, small host offset, and high host SFR favor a low-luminosity collapsar progenitor for CDF-S XT2. While these HST and JWST observations shed light on the host galaxies of XT1 and XT2, and by extension, on the nature of FXTs, a unique explanation for both sources remains elusive., Comment: The manuscript was submitted to Astronomy & Astrophysics in August 2024

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

GRBs produced by the collapse of massive stars are usually found near the most prominent star-forming regions of star-forming galaxies. GRB 171205A happened in the outskirts of a spiral galaxy, a peculiar location in an atypical GRB host. In this paper we present a highly-resolved study of the molecular gas of this host, with CO(1-0) observations from ALMA. We compare with GMRT atomic HI observations, and with data at other wavelengths to provide a broad-band view of the galaxy. The ALMA observations have a spatial resolution of 0.2" and a spectral resolution of 10 km/s, observed when the afterglow had a flux density of ~53 mJy. This allowed a molecular study both in emission and absorption. The HI observations allowed to study the host galaxy and its extended environment. The CO emission shows an undisturbed spiral structure with a central bar, and no significant emission at the location of the GRB. Our CO spectrum does not reveal any CO absorption, with a column density limit of < 10^15 cm^-2. This argues against the progenitor forming in a massive molecular cloud. The molecular gas traces the galaxy arms with higher concentration in the regions dominated by dust. The HI gas does not follow the stellar light or the molecular gas and is concentrated in two blobs, with no emission towards the centre of the galaxy, and is slightly displaced towards the southwest of the galaxy, where the GRB exploded. Within the extended neighbourhood of the host galaxy, we identify another prominent HI source at the same redshift, at a projected distance of 188 kpc. Our observations show that the progenitor of this GRB is not associated to a massive molecular cloud, but more likely related to low-metallicity atomic gas. The distortion in the HI gas field is indicator of an odd environment that could have triggered star formation and could be linked to a past interaction with the companion galaxy., Comment: 13 pages, 10 figures, 8 tables, A&A submitted after 1st referee review

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

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

Bulk viscosity of a plasma consisting of strongly coupled diatomic ions is computed using molecular dynamics simulations. The simulations are based on the rigid rotor one-component plasma, which is introduced as a model system that adds two degrees of molecular rotation to the traditional one-component plasma. It is characterized by two parameters: the Coulomb coupling parameter, $\Gamma$, and the bond length parameter, $\Omega$. Results show that the long-range nature of the Coulomb potential can lead to long rotational relaxation times, which in turn yield large values for bulk viscosity. The bulk-to-shear viscosity ratio is found to span from small to large values depending on the values of $\Gamma$ and $\Omega$. Although bulk viscosity is often neglected in plasma modeling, these results motivate that it can be large in molecular plasmas with rotational degrees of freedom.

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

The Epoch of Reionization (EoR) is a key period of cosmological history in which the intergalactic medium (IGM) underwent a major phase change from being neutral to almost completely ionized. Gamma-ray bursts (GRBs) are luminous and unique probes of their environments that can be used to study the timeline for the progression of the EoR. Here we present a detailed analysis of the ESO Very Large Telescope X-shooter spectrum of GRB 210905A, which resides at a redshift of z ~ 6.3. We focus on estimating the fraction of neutral hydrogen, xHI, on the line of sight to the host galaxy of GRB 210905A by fitting the shape of the Lyman-alpha damping wing of the afterglow spectrum. The X-shooter spectrum has a high signal to noise ratio, but the complex velocity structure of the host galaxy limits the precision of our conclusions. The statistically preferred model suggests a low neutral fraction with an 3-sigma upper limit of xHI < 0.15, indicating that the IGM around the GRB host galaxy is mostly ionized. We discuss complications in current analyses and potential avenues for future studies of the progression of the EoR and its evolution with redshift., Comment: 12 pages, 11 figures, submitted to Monthly Notices of the Royal Astronomical Society

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

Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secure example of this phenomenon to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising (t ~ 2 hr), luminous (M_UV,peak ~ -24.4 mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing strong limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of Gamma_0 ~ 30-50, significantly lower than in any well-observed GRB and consistent with the theoretically-predicted "dirty fireball" scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-Gamma, off-axis from a classical high-Gamma jet core. This event represents a milestone in orphan afterglow searches, demonstrating that luminous afterglows with weak or no detectable gamma-ray radiation exist in nature and can be discovered by high-cadence optical surveys., Comment: Submitted to MNRAS

Recent work has shown that ions are strongly coupled in atmospheric pressure plasmas when the ionization fraction is sufficiently large, leading to a temperature increase from disorder-induced heating that is not accounted for in standard modelling techniques. Here, we extend this study to molecular plasmas. A main finding is that the energy gained by ions in disorder-induced heating gets spread over both translational and rotational degrees of freedom on a nanosecond timescale, causing the final ion and neutral gas temperatures to be lower in the molecular case than in the atomic case. A model is developed for the equilibrium temperature that agrees well with molecular dynamics simulations. The model and simulations are also applied to pressures up to ten atmospheres. We conclude that disorder-induced heating is a significant and predictable phenomenon in molecular atmospheric pressure plasmas.

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!

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

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

Future detection of high-redshift gamma-ray bursts (GRBs) will be an important tool for studying the early Universe. Fast and accurate redshift estimation for detected GRBs is key for encouraging rapid follow-up observations by ground- and space-based telescopes. Low-redshift dusty interlopers pose the biggest challenge for GRB redshift estimation using broad photometric bands, as their high extinction can mimic a high-redshift GRB. To assess false alarms of high-redshift GRB photometric measurements, we simulate and fit a variety of GRBs using phozzy, a simulation code developed to estimate GRB photometric redshifts, and test the ability to distinguish between high- and low-redshift GRBs when using simultaneously observed photometric bands. We run the code with the wavelength bands and instrument parameters for the Photo-z Infrared Telescope (PIRT), an instrument designed for the Gamow mission concept. We explore various distributions of host galaxy extinction as a function of redshift, and their effect on the completeness and purity of a high-redshift GRB search with the PIRT. We find that for assumptions based on current observations, the completeness and purity range from $\sim 82$ to $88\%$ and from $\sim 84$ to $>99\%$, respectively. For the priors optimized to reduce false positives, only $\sim 0.6\%$ of low-redshift GRBs will be mistaken as a high-redshift one, corresponding to $\sim 1$ false alarm per 500 detected GRBs., Comment: 14 pages, 15 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 chemical enrichment of dust and metals in the interstellar medium (ISM) of galaxies throughout cosmic time is one of the key driving processes of galaxy evolution. Here we study the evolution of the gas-phase metallicities, dust-to-gas (DTG), and dust-to-metal (DTM) ratios of 36 star-forming galaxies at $1.7 < z < 6.3$ probed by gamma-ray bursts (GRBs). We compile all GRB-selected galaxies with intermediate (R=7000) to high (R>40,000) resolution spectroscopic data for which at least one refractory (e.g. Fe) and one volatile (e.g. S or Zn) element have been detected at S/N>3. This is to ensure that accurate abundances and dust depletion patterns can be obtained. We first derive the redshift evolution of the dust-corrected, absorption-line based gas-phase metallicity [M/H]$_{\rm tot}$ in these galaxies, for which we determine a linear relation with redshift ${\rm [M/H]_{tot}}(z) = (-0.21\pm 0.04)z -(0.47\pm 0.14)$. We then examine the DTG and DTM ratios as a function of redshift and through three orders of magnitude in metallicity, quantifying the relative dust abundance both through the direct line-of-sight visual extinction $A_V$ and the derived depletion level. We use a novel method to derive the DTG and DTM mass ratios for each GRB sightline, summing up the mass of all the depleted elements in the dust-phase. We find that the DTG and DTM mass ratios are both strongly correlated with the gas-phase metallicity and show a mild evolution with redshift as well. While these results are subject to a variety of caveats related to the physical environments and the narrow pencil-beam sightlines through the ISM probed by the GRBs, they provide strong implications for studies of dust masses to infer the gas and metal content of high-redshift galaxies, and particularly demonstrate the large offset from the average Galactic value in the low-metallicity, high-redshift regime., Comment: Accepted in A&A

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