Search

Fast X-ray transients (FXTs) are extragalactic bursts of X-rays first identified in archival X-ray data, and now routinely discovered by the Einstein Probe in real time, which is continuously surveying the night sky in the soft ($0.5 - 4$ keV) X-ray regime. In this Letter, we report the discovery of the second optical counterpart (AT2024gsa) to an FXT (EP240414a). EP240414a is located at a projected radial separation of 27 kpc from its likely host galaxy at $z = 0.4018 \pm 0.0010$. The optical light curve of AT2024gsa displays three distinct components. The initial decay from our first observation is followed by a re-brightening episode, displaying a rapid rise in luminosity to an absolute magnitude of $M_r \sim -21$ after two rest-frame days. While the early optical luminosity and decline rate is similar to luminous fast blue optical transients, the colour temperature of AT2024gsa is distinctly red and we show that the peak flux is inconsistent with a thermal origin. The third component peaks at $M_i \sim -19$ at $\gtrsim 16$ rest-frame days post-FXT, and is compatible with an emerging supernova. We fit the $riz$-band data with a series of power laws and find that the decaying components are in agreement with gamma-ray burst afterglow models, and that the re-brightening may originate from refreshed shocks. By considering EP240414a in context with all previously reported known-redshift FXT events, we propose that Einstein Probe FXT discoveries may all result from high-redshift gamma-ray bursts, and thus are distinct from the previously discovered lower redshift, lower luminosity population of FXTs., Comment: To be submitted to ApJL

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

Kilonovae are key to advancing our understanding of r-process nucleosynthesis. To date, only two kilonovae have been spectroscopically observed, AT 2017gfo and AT 2023vfi. Here, we present an analysis of the James Webb Space Telescope spectra obtained +29 and +61 days post-merger for AT 2023vfi. After re-reducing and photometrically flux-calibrating the data, we empirically model the observed X-ray to mid-infrared continua with a power law and a blackbody, to replicate the non-thermal afterglow and apparent thermal continuum $\gtrsim 2 \, \mu$m. We fit Gaussians to the apparent emission features, obtaining line centroids of $20218_{-38}^{+37}$, $21874 \pm 89$ and $44168_{-152}^{+153}$\,\AA, and velocity widths spanning $0.057 - 0.110$\,c. These line centroid constraints facilitated a detailed forbidden line identification search, from which we shortlist a number of r-process species spanning all three r-process peaks. We rule out Ba II and Ra II as candidates and propose Te I-III, Er I-III and W III as the most promising ions for further investigation, as they plausibly produce multiple emission features from one (W III) or multiple (Te I-III, Er I-III) ion stages. We compare to the spectra of AT 2017gfo, which also exhibit prominent emission at $\sim 2.1 \, \mu$m, and conclude that [Te III] $\lambda$21050 remains the most plausible cause of the observed $\sim 2.1 \, \mu$m emission in both kilonovae. However, the observed line centroids are not consistent between both objects, and they are significantly offset from [Te III] $\lambda$21050. The next strongest [Te III] transition at 29290\,\AA\ is not observed, and we quantify its detectability. Further study is required, with particular emphasis on expanding the available atomic data to enable quantitative non-LTE spectral modelling., Comment: Main text: 21 pages, 14 figures, 3 tables. Appendices: 7 pages, 20 figures. Submitted to MNRAS. Comments welcome

Fast X-ray Transients (FXTs) are extragalactic bursts of soft X-rays first identified >10 years ago. Since then, nearly 40 events have been discovered, although almost all of these have been recovered from archival Chandra and XMM-Newton data. To date, optical sky surveys and follow-up searches have not revealed any multi-wavelength counterparts. The Einstein Probe, launched in January 2024, has started surveying the sky in the soft X-ray regime (0.5-4 keV) and will rapidly increase the sample of FXTs discovered in real time. Here, we report the first discovery of both an optical and radio counterpart to a distant FXT, the fourth source publicly released by the Einstein Probe. We discovered a fast-fading optical transient within the 3 arcmin localisation radius of EP240315a with the all-sky optical survey ATLAS, and our follow-up Gemini spectrum provides a redshift, z=4.859+/-0.002. Furthermore, we uncovered a radio counterpart in the S-band (3.0 GHz) with the MeerKAT radio interferometer. The optical (rest-frame UV) and radio luminosities indicate the FXT most likely originates from either a long gamma-ray burst or a relativistic tidal disruption event. This may be a fortuitous early mission detection by the Einstein Probe or may signpost a mode of discovery for high-redshift, high-energy transients through soft X-ray surveys, combined with locating multi-wavelength counterparts., Comment: Updated to match version accepted for publication in ApJL (17 pages, 4 figures, 2 tables)

GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverage of the skymap beginning within 1 hour and 3 hours respectively of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 190425 has been suggested, given their spatial and temporal coincidence. The smaller sky localisation area of FRB 190425 and its dispersion measure have led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 +/- 10 Mpc. Our optical imaging covered the galaxy 6.0 hrs after GW190425 was detected and 3.5 hrs after the FRB 190425. No optical emission was detected and further imaging at +1.2 and +13.2 days also revealed no emission. If the FRB 190425 and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 190425 from a BNS merger involves a supramassive magnetised neutron star spinning down by dipole emission on the timescale of hours. We show that magnetar enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB-GW link for this source., Comment: Submitted to MNRAS, 20th Sept 2023, 9 pages

Binary neutron star mergers are the first confirmed site of rapid neutron capture (r-process) element nucleosynthesis. The kilonova AT2017gfo is the only electromagnetic counterpart of a neutron star merger spectroscopically observed. We analyse the entire spectral sequence of AT2017gfo (from merger to +10.4 days) and identify seven emission-like features. We confirm that the prominent 1.08um feature can be explained by the Sr II near-infrared triplet evolving from a P-Cygni profile through to pure emission. We calculate the expected strength of the [Sr II] doublet and show that its absence requires highly clumped ejecta. Near-infrared features at 1.58 and 2.07um emerge after three days and become more prominent as the spectra evolve. We model these as optically thick P-Cygni profiles and alternatively as pure emission features (with FWHM = 35600 +/- 6600 km/s), and favour the latter interpretation. The profile of the strong 2.07um emission feature is best reproduced with two lines, centred at 2.059 and 2.135um. We search for candidate ions for all prominent features in the spectra. Strong, permitted transitions of La III, Ce III, Gd III, Ra II and Ac I are plausible candidates for the emission features. If any of these features are produced by intrinsically weak, forbidden transitions, we highlight candidate ions spanning the three r-process peaks. The second r-process peak elements Te and I have plausible matches to multiple features. We highlight the need for more detailed and quantitative atomic line transition data., Comment: Main text: 19 pages, 7 figures, 6 tables. Appendices: 10 pages, 3 tables. Submitted to MNRAS. Comments welcome

Long duration gamma-ray bursts (GRBs) are powerful cosmic explosions, signaling the death of massive stars. Among them, GRB 221009A is by far the brightest burst ever observed. Due to its enormous energy ($E_\textrm{iso}\!\approx$10$^{55}$ erg) and proximity ($z\!\approx$0.15), GRB 221009A is an exceptionally rare event that pushes the limits of our theories. We present multi-wavelength observations covering the first three months of its afterglow evolution. The X-ray brightness decays as a power-law with slope $\approx\!t^{-1.66}$, which is not consistent with standard predictions for jetted emission. We attribute this behavior to a shallow energy profile of the relativistic jet. A similar trend is observed in other energetic GRBs, suggesting that the most extreme explosions may be powered by structured jets launched by a common central engine., Comment: Submitted version. 53 pages, 9 figures, 6 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

The kilonova (KN) associated with the binary neutron star (BNS) merger GW170817 is the only known electromagnetic counterpart to a gravitational wave source. Here we produce a sequence of radiative transfer models (using $\textsc{tardis}$) with updated atomic data, and compare them to accurately calibrated spectra. We use element compositions from nuclear network calculations based on a realistic hydrodynamical simulation of a BNS merger. We show that the blue spectrum at +1.4 days after merger requires a nucleosynthetic trajectory with a high electron fraction. Our best-fitting model is composed entirely of first $r$-process peak elements (Sr & Zr) and the strong absorption feature is reproduced well by Sr$\,\textsc{ii}$ absorption. At this epoch, we set an upper limit on the lanthanide mass fraction of $X_{\textsc{ln}} \lesssim 5 \times 10^{-3}$. In contrast, all subsequent spectra from $+2.4 - 6.4$ days require the presence of a modest amount of lanthanide material ($X_{\textsc{ln}} \simeq 0.05^{+0.05}_{-0.02}$), produced by a trajectory with $Y_{\rm e} = 0.29$. This produces lanthanide-induced line blanketing below 6000$\,$\AA, and sufficient light $r$-process elements to explain the persistent strong feature at $\sim 0.7 - 1.0 \, \micron$ (Sr$\,\textsc{ii}$). The composition gives good matches to the observed data, indicating that the strong blue flux deficit results in the near-infrared (NIR) excess. The disjoint in composition between the first epoch and all others indicates either ejecta stratification, or the presence of two distinct components of material. This further supports the `two-component' kilonova model, and constrains the element composition from nucleosynthetic trajectories. The major uncertainties lie in availability of atomic data and the ionisation state of the expanding material., Comment: 20 pages, 14 figures and 3 tables. Accepted for publication in MNRAS

We present multi-wavelength follow-up observations of the ATLAS discovered faint Iax supernova SN 2020kyg that peaked at an absolute magnitude of $M_g \approx -14.9 \pm 0.2$, making it another member of the faint Iax supernova population. The bolometric light curve requires only $\approx 7 \times 10^{-3}$ M$_{\odot}$ of radioactive $^{56}$Ni, with an ejected mass of $M_{\rm ej} \sim 0.4$ M$_{\odot}$ and a low kinetic energy of $E \approx 0.05 \pm 0.02 \times 10^{51}$ erg. We construct a homogeneous volume-limited sample of 902 transients observed by ATLAS within 100 Mpc during a 3.5 year span. Using this sample, we constrain the rates of faint Iax ($M_r \gtrsim -16$) events within 60 Mpc at $12^{+14}_{-8}\%$ of the SN Ia rate. The overall Iax rate, at $15^{+17}_{-9}\%$ of the Ia rate, is dominated by the low-luminosity events, with luminous SNe Iax ($M_r \lesssim -17.5$) like 2002cx and 2005hk accounting for only $0.9^{+1.1}_{-0.5}\%$ of the Ia rate (a 2$\sigma$ upper limit of approximately 3\%). We favour the hybrid CONe WD + He star progenitor channel involving a failed deflagration of a near Chandrasekhar mass white dwarf, expected to leave a bound remnant and a surviving secondary companion, as a candidate explanation for faint Iax explosions. This scenario requires short delay times, consistent with the observed environments of SNe Iax. Furthermore, binary population synthesis calculations have suggested rates of $1-18\%$ of the SN Ia rate for this channel, consistent with our rate estimates., Comment: Accepted to MNRAS after minor revision

Investigations of some type Iax supernovae have led to the suggestion that their ejecta must be layered to some degree. Such an ejecta structure has been argued as inconsistent with the well-mixed composition predicted by pure deflagrations. Based on explosion models, we create toy models in which the ejecta are artificially stratified and progressively mixed until a uniform composition is obtained. We find that models that are heavily mixed, containing burned and unburned material at all velocities, produce reasonably good agreement with SN 2012Z, for which a layered structure has been suggested. We also discuss how existing ejecta compositions determined for type Iax supernovae do not necessarily contradict pure deflagration models and may be consistent with a steeper density profile. We investigate previous claims that differences in line profile shapes may be due to strong blending, by presenting a series of models with different plasma states. These models indicate that blending could indeed explain differences in the observed profiles. Alternatively, stratification could also explain such differences, however all of our models indicate that this does not necessarily require stratification in abundance. Sufficient stratification in ionisation state can be achieved even for a well-mixed model. Based on our analysis, we demonstrate that there is insufficient evidence to suggest the ejecta of type Iax supernovae must be layered and therefore argue the pure deflagration scenario is not ruled out, even for the brightest type Iax supernovae. Our analysis does not indicate the ejecta cannot be layered to some degree, but observations within days of explosion are necessary to determine the extent to which the outer ejecta could be layered., Comment: 19 pages, 12 figures, 8 tables. Accepted for publication in MNRAS

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

Recent studies have shown how the distribution of $^{56}$Ni within the ejecta of type Ia supernovae can have profound consequences on the observed light curves. Observations at early times can therefore provide important details on the explosion physics in thermonuclear supernovae. We present a series of radiative transfer calculations that explore variations in the $^{56}$Ni distribution. Our models also show the importance of the density profile in shaping the light curve, which is often neglected in the literature. Using our model set, we investigate the observations that are necessary to determine the $^{56}$Ni distribution as robustly as possible within the current model set. We find that this includes observations beginning at least $\sim$14 days before $B$-band maximum, extending to approximately maximum light with a high ($\lesssim$3 day) cadence, and in at least one blue and one red band are required (such as $B$ and $R$, or $g$ and $r$). We compare a number of well-observed type Ia supernovae that meet these criteria to our models and find that the light curves of $\sim$70-80\% of objects in our sample are consistent with being produced solely by variations in the $^{56}$Ni distributions. The remaining supernovae show an excess of flux at early times, indicating missing physics that is not accounted for within our model set, such as an interaction or the presence of short-lived radioactive isotopes. Comparing our model light curves and spectra to observations and delayed detonation models demonstrates that while a somewhat extended $^{56}$Ni distribution is necessary to reproduce the observed light curve shape, this does not negatively affect the spectra at maximum light. Investigating current explosion models shows that observations typically require a shallower decrease in the $^{56}$Ni mass towards the outer ejecta than is produced for models of a given $^{56}$Ni mass., Comment: 31 pages, 32 figures, 2 tables. Accepted for publication in A&A. Model light curves available at https://github.com/MarkMageeAstro/TURTLS-Light-curves