S250818k: Sub-threshold GW Candidate

• S250818k is a sub-threshold gravitational-wave event potentially arising from a binary neutron star merger involving a sub-solar-mass neutron star.
• Extensive electromagnetic follow-ups across optical, radio, and X-ray bands identified SN 2025ulz, a Type IIb supernova initially mimicking kilonova emission.
• The event highlights challenges in multi-messenger candidate identification and underscores the need for advanced ranking methods in crowded transient searches.

S250818k is the designation for a sub-threshold gravitational-wave (GW) candidate reported by the LIGO-Virgo-KAGRA (LVK) collaboration on 2025 August 18. The event was consistent with a binary neutron star (BNS) merger involving a potentially sub-solar-mass neutron star. Prompted by its astrophysical and observational significance—particularly the pursuit of electromagnetic (EM) counterparts to rare low-mass compact object mergers—S250818k triggered extensive multi-messenger follow-up across optical, radio, and X-ray bands. These searches led to the identification and intensive investigation of the transient AT2025ulz (subsequently classified as SN 2025ulz), whose early-time behavior mimicked kilonova emission but was ultimately recognized as a Type IIb supernova. The association between S250818k and SN 2025ulz remains unproven, yet the episode exemplifies the challenges of multi-wavelength counterpart identification, the astrophysical implications of sub-solar-mass GW events, and the landscape of future searches for "superkilonovae"—energetic transients postulated to be powered by hierarchical neutron-star mergers within core-collapse supernovae (Hall et al., 27 Oct 2025, Kasliwal et al., 27 Oct 2025, Franz et al., 20 Oct 2025, Hall et al., 28 Oct 2025, O'Dwyer et al., 6 Apr 2026, Ackley et al., 4 May 2026).

1. Gravitational-Wave Detection of S250818k

S250818k was issued as a low-significance public alert, with a false-alarm probability of 71% (i.e., 29% probability of astrophysical origin) and an estimated binary neutron star merger probability $P(\mathrm{BNS}) \approx 0.29$. The GW parameters included a source-frame chirp mass constrained to $\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$, making the event potentially a sub-solar-mass neutron star merger (Hall et al., 27 Oct 2025, Kasliwal et al., 27 Oct 2025, Franz et al., 20 Oct 2025). The sky localization encompassed $786$–$1500$ deg$^2$ (updated $90\%$ credible area: $949$ deg$^2$), with a luminosity distance posterior centered at $237\pm62$ Mpc (BAYESTAR/LALInference), and a $2\sigma$ host-compatible redshift of $\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$0 ($\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$1 Mpc) (Hall et al., 27 Oct 2025, Yang et al., 21 Oct 2025, Ackley et al., 4 May 2026). The low chirp mass is unprecedented among known BNS systems and motivated theoretical discussion of origin channels such as core fission and accretion-disk fragmentation (collapsar disk) (Kasliwal et al., 27 Oct 2025, Wu et al., 29 Apr 2026).

2. Electromagnetic Follow-Up and Identification of SN 2025ulz

In response to S250818k, wide-field surveys (ZTF, Pan-STARRS, ATLAS) launched a multi-band optical search, rapidly uncovering new transients (Gillanders et al., 1 Oct 2025). ZTF detected AT2025ulz roughly 3 hours post-trigger; the transient exhibited a fast, blue decline, resembling kilonova expectations (by analogy to GW170817/AT2017gfo)—declining by up to $\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$2 mag d$\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$3 in $\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$4 and $\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$5 mag d$\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$6 in $\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$7 over the first two days (Yang et al., 21 Oct 2025, Franz et al., 20 Oct 2025). Early color evolution ($\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$8 to $\mathcal{M}_{\mathrm{chirp}} < 0.87\,M_\odot$9 mag over 48 hr) and brightness (factor $786$0 over AT2017gfo in $786$1) were only marginally compatible with kilonova models, requiring implausibly large ejecta masses ($786$2 at $786$3–$786$4) (Yang et al., 21 Oct 2025, Hall et al., 28 Oct 2025). Spectroscopic monitoring (GTC, Gemini-N, VLT/X-shooter, MUSE, FTW/3KK, HST) tracked the transient’s evolution, revealing by $786$5 d a rebrightening and emergence of broad P-Cygni H$786$6 features at velocities of $786$7–$786$8 km s$786$9, establishing its nature as a young, stripped-envelope Type IIb supernova (hereafter SN 2025ulz) (Yang et al., 21 Oct 2025, Hall et al., 27 Oct 2025, Hall et al., 28 Oct 2025, Ackley et al., 4 May 2026).

Ancillary photometry (HST, VLT, CTIO/DECam, NOT/ALFOSC, MMT/Binospec, SOAR/Goodman, T80N-Cam, CFHT) mapped the shock-cooling tail, plateau, and radioactive peak (Franz et al., 20 Oct 2025, Ackley et al., 4 May 2026). Light-curve modeling with Arnett-type $1500$0Ni + shock cooling fits yielded $1500$1, $1500$2 erg, and $1500$3, incompatible with standard kilonovae (Ackley et al., 4 May 2026). The explosion epoch was inferred to precede the GW trigger by approximately $1500$4 days, ruling out causality for kilonova association (Ackley et al., 4 May 2026).

3. Statistical Association and Host-Galaxy Properties

Comprehensive host-galaxy studies exploited DESI’s pre-existing spectroscopic catalog to rapidly measure the redshift directly beneath SN 2025ulz as $1500$5, corresponding to $1500$6 Mpc (Hall et al., 27 Oct 2025). Volumetric comparison with the GW posterior localization yielded an integral overlap $1500$7–$1500$8, while the reference event GW170817/GRB170817A provides $1500$9 (Hall et al., 27 Oct 2025). Morphology and spectral energy distribution (SED) analysis indicated a star-forming, dusty host of mass $^2$0, $^2$1 mag, metallicity $^2$2, and SFR $^2$3, consistent with the broader class of core-collapse SN and short GRB hosts but lacking specificity for a GW counterpart (Hall et al., 27 Oct 2025, Yang et al., 21 Oct 2025). Statistical assessment, based on the volumetric rate of Type IIb supernovae and shock-cooling phase duration, yields a $^2$4–$^2$5\% chance of observing at least one interloping SN IIb in a comparable GW search volume and time window (Kasliwal et al., 27 Oct 2025). This rate is much higher than for BNS mergers ($^2$6–$^2$7 Mpc$^2$8 yr$^2$9), making supernovae a dominant contamination channel in kilonova searches (Ackley et al., 4 May 2026).

4. Multi-Wavelength Constraints: Radio and X-ray Follow-Up

Deep radio follow-up (VLA/JAGWAR, uGMRT, MeerKAT, eMERLIN) and X-ray campaigns (Swift, XMM-Newton, Chandra) set stringent constraints on non-thermal emission possibly expected from a BNS merger, off-axis jet, or a "superkilonova" (O'Dwyer et al., 6 Apr 2026, O'Connor et al., 27 Oct 2025, Ackley et al., 4 May 2026). Key findings:

• VLA S/C/X-band coverage at 3–210 d yielded a faint, but statistically significant, SN counterpart at 6–10 GHz ($90\%$0Jy at day 89 in C-band, $90\%$1Jy combining X-band B-array epochs), with marginal evidence of a rise near 74 d (O'Dwyer et al., 6 Apr 2026).
• uGMRT at 8.5 d set F$90\%$263 $90\%$3Jy (3$90\%$4) upper limits.
• MeerKAT detections were dominated by host-galaxy flux.
• X-ray limits at 19.41 d: Chandra/ACIS-S $90\%$5 erg cm$90\%$6 s$90\%$7 (0.3–10 keV), constraining afterglow models.
• Combined, these limits rule out a GW170817-like afterglow for viewing angles $90\%$8 at 400 Mpc (O'Connor et al., 27 Oct 2025).

Radio emission is consistent with two scenarios: (1) a standard CSM-interacting SN IIb with shock velocity $90\%$9 km s$949$0 and mass-loss rate $949$1, or (2) an off-axis relativistic jet (requiring $949$2 erg, $949$3, $949$4), with model degeneracies unresolved due to sparse light-curve sampling (O'Dwyer et al., 6 Apr 2026). However, the absence of expected afterglow emission and plausible modeling favor the CSM interaction scenario for SN 2025ulz and argue against a physical connection between S250818k and a relativistic jet-powered transient (O'Dwyer et al., 6 Apr 2026, O'Connor et al., 27 Oct 2025).

5. Theoretical Frameworks: Superkilonova Scenarios and Hierarchical Mergers

The unusual GW parameter space (sub-solar $949$5) motivated consideration of non-standard progenitor scenarios. Chief among these are collapsar-disk fragmentation and core fission, both of which can result in sub-solar-mass neutron stars and associated supernovae with atypical EM signatures (Kasliwal et al., 27 Oct 2025, Wu et al., 29 Apr 2026).

Hierarchical mergers in a collapsar disk, as explored in Wu et al., predict multi-body dynamics where fragments coalesce to form a $949$6 neutron star, before a final eccentric merger with a central black hole. Numerical relativity simulations demonstrate that the final merger retains large eccentricity, $949$7 at merger, a strong discriminant compared to $949$8 expected from isolated binaries (Wu et al., 29 Apr 2026). GW waveforms are characterized by pericenter bursts and phase modulation, observable in the 50–300 Hz band with SNR$949$910.

Predicted EM counterparts in this framework include a blue, fast shock-cooling transient, rebrightening from $^2$0Ni decay, extended red NIR emission from r-process elements, and possible late-time radio/X-ray flares tied to pericenter passages of a non-circular merger (Wu et al., 29 Apr 2026, Kasliwal et al., 27 Oct 2025). Direct multi-messenger signatures, such as pericenter GW arrival times correlated with radio/optical flares, would provide compelling evidence for a hierarchical origin.

6. Implications for Kilonova Searches and Candidate Ranking Methodologies

The S250818k episode highlighted the diagnostic ambiguity of early-time photometry: for $^2$1 d post-merger, light curves and colors could plausibly match both kilonova and shock-cooling SN models (Hall et al., 28 Oct 2025, Franz et al., 20 Oct 2025). Only subsequent evolution (plateau, rebrightening, and spectroscopic emergence of H$^2$2 P Cygni) securely established the Type IIb nature.

To improve candidate triage, quantitative ranking frameworks have been developed. The scoring algorithm presented by (Franz et al., 20 Oct 2025), now implemented in the "TROVE" system, assesses transients via a composite score $^2$3, combining sky and distance overlap, artifact rejection, and photometric properties. SN 2025ulz rapidly dropped in rank as its photometric and spectroscopic behavior diverged from kilonova expectations. At all times after discovery, there existed $^2$4 better candidates by this metric, underscoring the need for real-time, multi-criteria evaluation to minimize resource expenditure on contaminating SNe (Franz et al., 20 Oct 2025).

7. Lessons Learned and Prospects for Future Multi-Messenger Campaigns

The S250818k/AT2025ulz campaign provides several critical lessons for GW-EM follow-up:

• Fast-evolving Type IIb SNe with extended envelopes ("shock cooling tails") are a major source of kilonova impostors, with much greater volumetric rates than BNS mergers (Ackley et al., 4 May 2026).
• Robust identification of kilonovae demands high-cadence, multi-band photometry—including NIR, early spectroscopic follow-up to $^2$5–23 mag, and quantitative modeling frameworks (Hall et al., 28 Oct 2025).
• Deep, multi-epoch radio and X-ray campaigns are necessary to rule out engine-driven outflow and afterglow scenarios, particularly as viewing angles and circum-merger environments differ from GW170817 (O'Dwyer et al., 6 Apr 2026, O'Connor et al., 27 Oct 2025).
• Spectroscopic host-galaxy catalogs (e.g., DESI) are essential for rapid volumetric association/rejection and facilitate host subtraction crucial for isolating transient features (Hall et al., 27 Oct 2025).
• For GW events with marginal significance and large localization volumes, follow-up strategies that combine rapid scoring with targeted observations toward high $^2$6 hosts and iterative ranking are effective (Franz et al., 20 Oct 2025).
• Identification of hierarchical merger channels and "superkilonovae" requires coordinated GW analysis using eccentric templates, EM counterpart searches for pericenter–flare correlation, and enhanced multi-messenger models (Wu et al., 29 Apr 2026, Kasliwal et al., 27 Oct 2025).

In summary, S250818k illustrates both the scientific potential and operational complexity of low-significance GW triggers: it motivates improved observational protocols, theoretical modeling of unusual merger channels, and deploys a new generation of candidate-vetting methodologies for future multi-messenger astrophysics.