On The Nature of Einstein Probe Transient EP250916a: Insights from X-ray, Optical, and Radio Observations

Abstract: We report multi-wavelength studies of the transient EP250916a, detected by the Einstein Probe on 2025 September 16. Located at low Galactic latitude, the source exhibited a rapid X-ray brightening, reaching an unabsorbed 0.5--10 keV flux of $(6.4 \pm 0.1) \times 10<sup>{-10}$ erg cm$<sup>{-2}$ s$<sup>{-1}$, followed by a plateau and a two-stage decay lasting over 40 days. Swift/XRT monitoring shows a persistently hard spectrum ($Γ\approx 1.6$--2.2) with only modest softening during decay, while a NuSTAR observation confirms a hard-state continuum extending up to 70 keV. Timing analysis of XMM-Newton data reveals a weak quasi-periodic oscillation (QPO) at $\sim$13 Hz. No other coherent pulsations or thermonuclear bursts are detected. Broadband spectral modeling favors a nonthermal power-law continuum with partial-covering absorption, and shows no significant thermal disk component. Optical imaging obtained with NOT/ALFOSC, LCO, and GaiaDR3 identifies two faint sources within the 2 arcsec Swift/XRT positional uncertainty. A MeerKAT observation at 1.28 GHz yielded no radio counterpart, with a 3$σ$ upper limit of 60 $μ$Jy beam$<sup>{-1}$. The combination of a long-lasting outburst, a hard nonthermal X-ray spectrum, a weak QPO detection, the absence of coherent timing features, and faint potential optical counterparts disfavors a stellar-flare or extragalactic origin and supports an accreting compact-object scenario. Comparisons with similar faint, hard-state transients place EP250916a within a growing population of low-luminosity, hard-state black hole X-ray binary candidates.

• The paper demonstrates that EP250916a is a low-luminosity, hard-state black hole candidate through coordinated X-ray, optical, and radio observations.
• It applies broadband spectral modeling and timing analysis, including detection of a weak ~13 Hz Type-C QPO, to infer the accretion dynamics.
• The study highlights the importance of rapid multi-observatory follow-up in revealing faint X-ray binaries that challenge traditional detection paradigms.

Multiwavelength Analysis of the Einstein Probe Transient EP250916a

Overview

This work presents a detailed multiwavelength analysis of the transient X-ray source EP250916a, discovered by the Einstein Probe (EP) on 2025 September 16. The study exploits coordinated observations across X-ray (Swift/XRT, NuSTAR, XMM-Newton), optical (NOT, LCO, Gaia), and radio (MeerKAT) domains to elucidate the physical properties, classification, and evolutionary state of the source. The combination of persistently hard X-ray spectrum, temporal structure, and weak QPO detection is leveraged to argue for an accreting compact-object scenario, with particular emphasis on the identification of EP250916a as a member of the growing population of faint, hard-state black hole X-ray binary (BHB) candidates (2606.08752).

Temporal Evolution and X-ray Timing Properties

The outburst evolution of EP250916a is characterized by an initial rapid X-ray rise, a plateau phase, and a two-stage decay spanning over forty days. Swift/XRT monitoring constrains the spectral hardness and reveals modest softening concurrent with flux decline, but the system remains in the hard X-ray regime throughout. The timing analysis with XMM-Newton identifies a weak QPO at approximately 13 Hz with a fractional rms amplitude of 7% and quality factor $Q \sim 1.6$, most consistent with a Type-C-like QPO commonly observed in intermediate/hard states of BHBs. No coherent pulsations or thermonuclear (Type I) bursts are observed, and NuSTAR places strict upper limits on possible QPO features inconsistent with strong neutron star signatures.

Broadband X-ray Spectral Modeling

The combined Swift/XRT, NuSTAR, and XMM-Newton spectral coverage (0.5–80 keV) is best described by an absorbed power-law continuum of photon index $\Gamma \sim 1.6$–2.2, with significant partial-covering absorption and a narrow Gaussian emission line at 0.71 keV. There is no evidence for a soft thermal disk component; attempts to model the spectrum with thermal or disk-blackbody components yield unphysical parameters (disk normalization implying sub-kilometer inner radius). The high-energy cut-off, derived from broadband fits, is broadly consistent with a Comptonizing corona with $E_{\rm cut} \gtrsim 200$ keV. The hydrogen column density is slightly elevated with respect to the Galactic value, implying additional local absorption.

Multiwavelength Counterpart Constraints

Optical imaging identifies two Gaia DR3 sources within the refined X-ray error circle, both consistent with faint Galactic stars. Astrometric and photometric analysis suggests the most plausible counterparts are mid-K to early-M dwarf stars or (if more distant) somewhat more massive donors. However, the lack of significant optical brightening argues against a luminous early-type companion or active, optically bright accretion disk at the time of observation. Gaia data alone are insufficient to distinguish between low- and intermediate-mass X-ray binary scenarios. MeerKAT radio follow-up yields non-detections with an upper limit of 60 μJy beam$^{-1}$ at 1.28 GHz.

Source Classification and Context

The collective evidence—a long-lived predominantly hard X-ray outburst, the absence of neutron star-specific timing or thermal features, lack of strong optical flaring, and faint/distant potential counterparts—strongly disfavors flare stars or extragalactic background origins and supports a scenario involving accretion onto a compact object within a binary. The persistent hard X-ray spectral state, non-detection of radio jets, low inferred X-ray luminosity ($L_{\rm X,peak} \sim 5 \times 10^{36}$ erg s$^{-1}$ at 8 kpc), and weak Type-C-like QPOs align with emerging characterizations of low-luminosity, hard-state transients, particularly those interpreted as “failed” outbursts of BHBs (e.g., never fully transitioning to disk-dominated states). The system is broadly similar to EP J182730.0-095633 but with a longer decay and less dramatic spectral/temporal variability.

The LR–LX limits derived from MeerKAT and X-ray data place EP250916a in a region overlapping both NS and BH LMXBs, but the overall phenomenology leans toward a BH accretor interpretation. The partial covering absorption and lack of Fe K or reflection features suggest a truncated disk possibly undergoing radiatively inefficient accretion. The observed quadruple broken power-law evolution of the lightcurve is consistent with disk-instability models of transient accretion.

Implications and Future Directions

EP250916a exemplifies the growing population of faint, hard-state X-ray binaries accessible to new-generation wide-field X-ray monitors. Its properties challenge detection and classification via traditional bright-outburst paradigms, emphasizing the importance of coordinated, rapid multiwavelength follow-up. The data reinforce the idea that the Galactic population of transient black hole binaries likely includes many low-luminosity/hard-state members that have thus far avoided detection.

The continuing systematic identification of such sources has significant implications for constraining the BH and NS binary mass functions, population synthesis models, and the low-accretion-rate regime of disk physics. Future developments should prioritize more sensitive dynamical mass and infrared studies to unambiguously confirm the nature of faint XRBs and ascertain their donor types, as well as high-resolution X-ray imaging (e.g., Chandra) to better localize and identify optical/infrared counterparts.

Conclusion

The multiwavelength study of EP250916a demonstrates a robust approach to the classification of faint X-ray binaries in the Galactic plane. The combined evidence points toward a compact-object accretion system, most consistently interpreted as a low-luminosity, hard-state black hole X-ray binary. This work establishes EP250916a as a significant contribution to the sample of faint hard-state transients and underscores the value of time-domain, multi-observatory strategies for the study of accreting compact objects (2606.08752).