Observations of gamma-ray burst (GRB) afterglows have begun to readily reveal contamination from a kilonova or a supernova. This contamination presents significant challenges towards traditional methods of inferring the properties of these phenomena from observations. Given current knowledge of kilonova and afterglow modelling, observations – as expected – with near-infrared bands and at early observing times provide the greatest diagnostic power for both observing the presence of a kilonova and inferences on its properties in GRB afterglows. However, contemporaneous observations in radio and X-ray are critical for reducing the afterglow parameter space and for more efficient parameter estimation. We compare different methods for fitting joint kilonova and afterglow observations under different scenarios. We find that ignoring the contribution of one source (even in scenarios where the source is sub-dominant) can lead to significantly biased estimated parameters but could still produce great light-curve fits that do not raise suspicion. This bias is also present for analyses that fit data where one source is ‘subtracted’. In most scenarios, the bias is smaller than the systematic uncertainty inherent to kilonova models but significant for afterglow parameters, particularly in the absence of high-quality radio and X-ray observations. Instead, we show that the most reliable method for inference in any scenario where contamination cannot be confidently dismissed is to jointly fit for both an afterglow and kilonova/supernova, and showcase a Bayesian framework to make this joint analysis computationally feasible.

Context. Gravitationally lensed type Ia supernovae (glSNe Ia) are unique astronomical tools that can be used to study cosmological parameters, distributions of dark matter, the astrophysics of the supernovae, and the intervening lensing galaxies themselves. A small number of highly magnified glSNe Ia have been discovered by ground-based telescopes such as the Zwicky Transient Facility (ZTF), but simulations predict that a fainter, undetected population may also exist. Aims. We present a systematic search for glSNe Ia in the ZTF archive of alerts distributed from June 1 2019 to September 1 2022. Methods. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry (with constraints on the peak absolute magnitude and the distance to a catalogue-matched galaxy, as examples) before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. The candidates were also cross-matched with the DESI spectroscopic catalogue. Results. Seven transients were identified that passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. Although superluminous supernovae (SLSNe) cannot be fully rejected as contaminants, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 ± 3 and 34 ± 1 days for the two events, respectively, which suggests that we have uncovered a population of glSNe Ia with longer time delays. Conclusions. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, such as the Rubin Observatory Legacy Survey of Space and Time.

Upgrades to existing gravitational-wave observatories have the potential to simultaneously constrain the nuclear equation of state and Hubble’s constant H₀ to per cent level with merging neutron star binaries. In practice, performing simultaneous inference of H₀ and the equation of state is limited computationally by the requirement to solve the equations of general-relativistic hydrostatic equilibrium millions of times. We develop a machine-learning model to solve the Tolman–Oppenheimer–Volkoff equations in less than a millisecond, and demonstrate its utility by performing direct inference of both equation of state and Hubble’s constant for synthetic neutron star merger signals with LIGO–Virgo–KAGRA operating at A+ sensitivities. We show that a population of 15 mergers observed with A+ allows for the radius of a 1.4 M_⊙ neutron star and H₀ to be constrained to km and at 90 per cent credible interval and 68 per cent credible interval, respectively. These constraints utilize only the gravitational-wave information to infer cosmological parameters; such numbers will be further improved with the addition of electromagnetic counterparts and/or galaxy catalogues.

With a small sample of fast X-ray transients (FXTs) with multiwavelength counterparts discovered to date, their progenitors and connections to γ-ray bursts (GRBs) and supernovae (SNe) remain ambiguous. Here, we present photometric and spectroscopic observations of SN 2025kg, the SN counterpart to the FXT EP 250108a. At z = 0.17641, this is the closest known SN discovered following an Einstein Probe (EP) FXT. We show that SN 2025kg’s optical spectra reveal the hallmark features of a broad-lined Type Ic SN. Its light-curve evolution and expansion velocities are comparable to those of GRB-SNe, including SN 1998bw, and two past FXT-SNe. We present JWST/NIRSpec spectroscopy taken around SN 2025kg’s maximum light, and find weak absorption due to He I 1.0830 μm and 2.0581 μm and a broad, unidentified emission feature at ∼4-4.5 μm. Further, we observe broadened Hα in optical data at 42.5 days that is not detected at other epochs, indicating interaction with H-rich material. From its light curve, we derive a 56Ni mass of 0.2-0.6 M⊙. Together with our companion Letter, our broadband data are consistent with a trapped or low-energy (≲1051 erg) jet-driven explosion from a collapsar with a zero-age main-sequence mass of 15-30 M⊙. Finally, we show that the sample of EP FXT-SNe supports past estimates that low-luminosity jets seen through FXTs are more common than successful (GRB) jets, and that similar FXT-like signatures are likely present in at least a few percent of the brightest Type Ic-BL SNe.

We present multiband observations and analysis of EP240801a, a low-energy, extremely soft gamma-ray burst (GRB) discovered on 2024 August 1 by the Einstein Probe (EP) satellite with a weak contemporaneous signal also detected by the Fermi Gamma-ray Burst Monitor (GBM). Optical spectroscopy of the afterglow, obtained by Gran Telescopio Canarias and Keck, identified the redshift of z = 1.6734. EP240801a exhibits a burst duration of 148 s in X-rays and 22.3 s in gamma rays, with X-rays leading by 80.61 s. Spectral lag analysis indicates that the gamma-ray signal arrived 8.3 s earlier than the X-rays. Joint spectral fitting of EP Wide-field X-ray Telescope and Fermi/GBM data yields an isotropic energy E γ , iso = ( 5.5 7 − 0.50 + 0.54 ) × 1 0 51 erg , a peak energy E peak = 14.9 0 − 4.71 + 7.08 keV , and a fluence ratio S(25-50 keV)/S ( 50 - 100 keV ) = 1.6 7 − 0.46 + 0.74 , classifying EP240801a as an X-ray flash (XRF). The host-galaxy continuum spectrum, inferred using Prospector, was used to correct its contribution for the observed outburst optical data. Unusual early R-band behavior and EP Follow-up X-ray Telescope observations suggest multiple components in the afterglow. Three models are considered: a two-component jet model, a forward-reverse shock model, and a forward shock model with energy injection. All three provide reasonable explanations. The two-component jet model and the energy injection model imply a relatively small initial energy and velocity of the jet in the line of sight, while the forward-reverse shock model remains typical. Under the two-component jet model, EP240801a may resemble GRB 221009A (BOAT) if the bright narrow beam is viewed on-axis. Therefore, EP240801a can be interpreted as an off-beam (narrow) jet or an intrinsically weak GRB jet. Our findings provide crucial clues for uncovering the origin of XRFs.

We present optical, radio, and X-ray observations of EP250108a/SN 2025kg, a broad-line Type Ic supernova (SN Ic-BL) accompanying an Einstein Probe (EP) fast X-ray transient at z = 0.176. EP250108a/SN 2025kg possesses a double-peaked optical light curve, and its spectrum transitions from a blue underlying continuum to a typical SN Ic-BL spectrum over time. We fit a radioactive decay model to the second peak of the optical light curve and find SN parameters that are consistent with the SN Ic-BL population, while its X-ray and radio properties are consistent with those of low-luminosity GRB (LLGRB) 060218/SN 2006aj. We explore three scenarios to understand the system’s multiwavelength emission: (a) SN ejecta interacting with an extended circumstellar medium (CSM), (b) the shocked cocoon of a collapsar-driven jet choked in its stellar envelope, and (c) the shocked cocoon of a collapsar-driven jet choked in an extended CSM. Models (b) and (c) can explain the optical light curve and are also consistent with the radio and X-ray observations. We favor model (c) because it can self-consistently explain both the X-ray prompt emission and first optical peak, but we do not rule out model (b). From the properties of the first peak in model (c), we find evidence that EP250108a/SN 2025kg interacts with an extended CSM and infer an envelope mass M_e ∼ 0.1 M_⊙ and radius R_e ∼ 4 × 10¹³ cm. EP250108a/SN 2025kg’s multiwavelength properties make it a close analog to LLGRB 060218/SN 2006aj and highlight the power of early follow-up observations in mapping the environments of massive stars prior to core collapse.

We present photometric and spectroscopic observations of SN 2020xga and SN 2022xgc, two hydrogen-poor superluminous supernovae (SLSNe-I) at z=-0.4296 and z = 0.3103, respectively, which show an additional set of broad Mg II absorption lines, blueshifted by a few thousands kilometer second-1 with respect to the host galaxy absorption system. Previous work interpreted this as due to resonance line scattering of the SLSN continuum by rapidly expanding circumstellar material (CSM) expelled shortly before the explosion. The peak rest-frame g-band magnitude of SN 2020xga is -22.30 ± 0.04 mag and of SN 2022xgc is -21.97 ± 0.05 mag, placing them among the brightest SLSNe-I. We used high-quality spectra from ultraviolet to near-infrared wavelengths to model the Mg II line profiles and infer the properties of the CSM shells. We find that the CSM shell of SN 2020xga resides at ∼1.3×1016 cm, moving with a maximum velocity of 4275 km s-1, and the shell of SN 2022xgc is located at ∼0.8×1016 cm, reaching up to 4400 km s-1. These shells were expelled ∼11 and ∼5 months before the explosions of SN 2020xga and SN 2022xgc, respectively, possibly as a result of luminous-blue-variable-like eruptions or pulsational pair instability (PPI) mass loss. We also analyzed optical photometric data and modeled the light curves, considering powering from the magnetar spin-down mechanism. The results support very energetic magnetars, approaching the mass-shedding limit, powering these SNe with ejecta masses of ∼7-9M⊙. The ejecta masses inferred from the magnetar modeling are not consistent with the PPI scenario pointing toward stars > 50M⊙ He-core; hence, alternative scenarios such as fallback accretion and CSM interaction are discussed. Modeling the spectral energy distribution of the host galaxy of SN 2020xga reveals a host mass of 107.8 M⊙, a star formation rate of 0.96-0.26+0.47 M⊙ yr-1, and a metallicity of ∼0.2 Z⊙

Stars are initially powered by the fusion of hydrogen to helium. These ashes serve as fuel in a series of stages1, 2–3, transforming massive stars into a structure of shells. These are composed of natal hydrogen on the outside and consecutively heavier compositions inside, predicted to be dominated by He, C/O, O/Ne/Mg and O/Si/S (refs. 4,5). Silicon and sulfur are fused into iron, leading to the collapse of the core and either a supernova explosion or the formation of a black hole6, 7, 8–9. Stripped stars, in which the outer hydrogen layer has been removed and the internal He-rich or even the C/O layer below it is exposed10, provide evidence for this shell structure and the cosmic element production mechanism it reflects. The supernova types that arise from stripped stars embedded in shells of circumstellar material (CSM) confirm this scenario11, 12, 13, 14–15. However, direct evidence for the most interior shells, which are responsible for producing elements heavier than oxygen, is lacking. Here we report the discovery of the supernova (SN) 2021yfj resulting from a star stripped to its O/Si/S-rich layer. We directly observe a thick, massive Si/S-rich shell, expelled by the progenitor shortly before the supernova explosion. Exposing such an inner stellar layer is theoretically challenging and probably requires a rarely observed mass-loss mechanism. This rare supernova event reveals advanced stages of stellar evolution, forming heavier elements, including silicon, sulfur and argon, than those detected on the surface of any known class of massive stars.

The ejection of neutron-rich matter is one of the most important consequences of a neutron star merger. While the bulk of the matter is ejected at fast, but non-relativistic velocities (⁠∼0.2c⁠), a small amount of mildly relativistic dynamic ejecta have been seen in a number of numerical simulations. Such ejecta can have far-reaching observational consequences ranging from the shock breakout burst of gamma-rays promptly after the merger, to an early (⁠∼1 h post-merger) blue kilonova precursor signal, to synchrotron emission years after the merger (‘kilonova afterglow’). These all potentially carry the imprint of the binary system parameters and the equation of state. By analysing Lagrangian simulations in full general relativity, performed with the code sphincs_bssn, we identify two ejection mechanisms for fast ejecta: (i) about 30 per cent of the ejecta with v>0.4c are ‘sprayed out’ from the shear interface between the merging stars and escape along the orbital plane and (ii) the remaining ∼70 per cent of the fast ejecta result from the central object ‘bouncing back’ after strong, general-relativistic compression. This ‘bounce component’ is ejected in a rather isotropic way and reaches larger velocities (by ∼0.1c⁠) so that its faster parts can catch up with and shock slower parts of the spray ejecta. Even for a case that promptly collapses to a black hole, we find fast ejecta with similar properties to the non-collapsing case, while slower matter parts are swallowed by the forming black hole. We discuss observational implications of these fast ejecta, including shock breakout and kilonova afterglow.

The nature of the minute-to-hour-long fast X-ray transients localized 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 1,600-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 and directly detect leaking ionizing Lyman continuum. The observed properties are consistent with EP240315a being a long-duration gamma-ray burst, and these observations suggest a possible interpretation in which a substantial fraction of the X-ray-transient population are lower-luminosity examples of similar events. If correct, then sensitive narrow-field searches could be a powerful complementary probe to traditional wide-field transient detection in the identification of samples of gamma-ray-burst-like events into the epoch of reionization.