2MASX J12571076+2724177 (J1257) X-ray Variability

• The paper demonstrates that J1257, a low-mass SMBH in the Coma cluster, exhibits distinctive 25 ks X-ray variability featuring both quasi-periodic modulations and flare-like events.
• Methodologically, it uses 17 pointed Chandra and XMM-Newton observations over 20 years, applying light curve modeling and Lomb–Scargle periodograms to characterize the variability.
• The analysis reveals spectral evolution with dual flux states and a softer-when-brighter trend, challenging conventional AGN classifications and hinting at possible QPO or QPE mechanisms.

Searching arXiv for the primary paper and closely related variability/QPE context papers. 2MASX J12571076+2724177 (J1257) is an active nucleus in the Coma cluster at redshift $z = 0.02068$ that has repeatedly exhibited peculiar short-term X-ray variability over a baseline of approximately 20 years. In 17 pointed \textit{Chandra} and \textit{XMM-Newton} observations, the source shows repeated variability on typical timescales of $\simeq 20$–25 ks, together with hints of a softer-when-brighter behavior and two well-separated flux states. The available dataset does not permit a definitive classification, and the variability could either represent a quasi-periodic oscillation at particularly low frequency or be associated with quasi-periodic eruptions in an AGN with peculiar spectral properties (Imbrogno et al., 3 Oct 2025).

1. Source identification and nuclear classification

J1257 is coincident with the galactic center in 2MASS and Gaia, confirming it as the SMBH host. The Gaia J2000 coordinates used for all extractions are RA $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$, Dec $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$. It lies in the low-mass tail of the SMBH population, with a broad-line virial mass estimate $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$, corresponding to $M \approx 2 \times 10^6\, M_\odot$ (Imbrogno et al., 3 Oct 2025).

Its optical classification is debated. Several optical studies classify it as a Seyfert 1, whereas infrared/optical SED-based work suggests a Seyfert 2, raising the possibility of changing-look behavior. This ambiguity is important because the X-ray phenomenology is being interpreted in the context of accretion-state changes, obscuration, and transient or quasi-periodic inner-flow structure rather than as a purely stationary Seyfert continuum.

The low black-hole mass is central to the interpretation of the timing signal. The source is discussed as part of a small but growing set of relatively low-mass SMBHs showing extreme short-timescale X-ray variability. In that sense, J1257 occupies a regime where characteristic dynamical and thermal disk timescales are short enough for hour-scale modulation to be astrophysically plausible, but where the observed recurrence frequency is still unusually low when compared with other SMBH timing candidates of similar mass.

2. X-ray observing baseline and reduction methodology

The analysis is based on 17 pointed X-ray observations with \textit{Chandra} and \textit{XMM-Newton} over $\sim 20$ years. The \textit{XMM-Newton} set spans 2000, 2006, 2010, and 2012; the \textit{Chandra} sample spans 2010 and 2020. The authors focus in particular on a three-observation \textit{Chandra} sequence in March 2020, designated C1 and composed of ObsIDs 22648, 22649, and 23182, because it provides the best coverage of the repeating short-term variability (Imbrogno et al., 3 Oct 2025).

Data reduction followed standard procedures. \textit{Chandra} data were reprocessed with CIAO v4.16 and CALDB v4.11.3, with VFAINT mode cleaning enabled (\texttt{checkvfpha=yes}). Source events were extracted in the 0.5–7 keV band from 4\arcsec-radius circular regions, except for off-axis ObsID 12887, for which a $7^{\prime\prime}\times 5^{\prime\prime}$ elliptical region was used. Background regions were nearby, source-free circles of 60\arcsec radius. Spectral products were generated with \texttt{specextract}, ensuring $\ge 95\%$ source contribution, and ObsID 23361 was excluded from spectroscopy due to poor statistics.

\textit{XMM-Newton} data were processed with SAS v21. Event selection used PATTERN $\le 4$ for pn and $\simeq 20$0 for MOS. High-energy ($\simeq 20$1 keV) flares were filtered with typical thresholds of $\simeq 20$2 for pn and $\simeq 20$3 for MOS, with slightly higher thresholds for selected high-background ObsIDs. Source events were extracted in the 0.3–10 keV band from 15\arcsec-radius circular apertures, with background from 40\arcsec circles on the same CCD. Vignetting and background corrections were applied using \texttt{epiclccorr}, and spectra were rebinned to $\simeq 20$4 count per bin. No barycentric corrections were applied because the variability timescales are tens of kiloseconds, and no pile-up issues were reported.

These reduction choices matter because the timing inference is limited by sparse cycles, heterogeneous background conditions, different soft-band responses, and gaps between observations. The paper explicitly notes that this heterogeneity prevents more rigorous regression-based or power-spectral significance assessments in several parts of the analysis.

3. Repeating short-timescale variability

J1257 repeatedly exhibits short-timescale X-ray variability on $\simeq 20$5–25 ks. In the combined 0.5–7 keV \textit{Chandra} light curve of the C1 sequence, a sinusoid plus a constant yields a period $\simeq 20$6 ks, with $\simeq 20$7. The fit quality is poor because of a flare-like deviation near the end of the sequence; adding a Gaussian component to model that flare improves the fit to $\simeq 20$8 (Imbrogno et al., 3 Oct 2025).

Lomb–Scargle periodograms computed with Stingray/hendrics show a broad peak centered at $\simeq 20$9 Hz, consistent with quasi-coherent variability. The paper emphasizes, however, that the low-frequency regime precludes a robust white-noise false-alarm probability and that the limited number of observed cycles, together with the observational gaps, complicates any formal significance assessment. The period inferred from the sinusoidal fit is consistent with the expected relation $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$0: a timescale $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$1–$= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$2 s corresponds to $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$3–$= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$4 Hz, broadly compatible with the periodogram peak.

Phase-folding all 2020 \textit{Chandra} data at $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$5 ks reveals a semi-amplitude $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$6 in 0.5–2 keV and $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$7 in 2–7 keV. This indicates a modestly larger modulation in the soft band, although the difference is only marginally significant at approximately $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$8. The waveform is not strictly sinusoidal, and flare-like events recur near minima.

A distinct flare-like event is seen in \textit{XMM-Newton} ObsID 0691610201. The 0.3–10 keV light curve shows a low state lasting $= 12^\mathrm{h}\ 57^\mathrm{m}\ 10.76^\mathrm{s}$9 ks at $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$0, followed by a rapid $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$1 increase over $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$2 ks and a $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$3 ks high-flux phase before declining. Using $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$4, where $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$5–10 keV and $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$6–1 keV, the hard-to-soft count-rate ratio is $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$7 outside the high-flux phase and $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$8 during the flare, consistent within $= +27^\circ\ 24^\prime\ 17.7^{\prime\prime}$9. Similar flare amplitudes and durations are seen in the late part of \textit{Chandra} ObsID 23182.

The timing phenomenology therefore consists of two linked components: quasi-coherent modulation on $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$0 ks scales and superposed flare-like events. The paper treats both as potentially intrinsic manifestations of a single accretion-flow process rather than as unrelated variability classes.

4. Spectral structure and flux-dependent behavior

In most epochs, J1257 is well described by a power law absorbed by Galactic foreground, modeled as TBabs with $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$1 fixed to $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$2, plus optional components when warranted by residuals (Imbrogno et al., 3 Oct 2025). A soft thermal excess, modeled as a blackbody with $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$3 keV, is detected mainly before 2012; later \textit{Chandra} epochs lack the soft sensitivity needed to constrain it. Intrinsic neutral absorption in excess of the Galactic component is required in several recent observations, with $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$4–$\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$5. The photon index spans $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$6–2.0, indicating genuine spectral evolution rather than unmodeled absorption.

Across the campaign, the absorption-corrected 0.3–10 keV luminosity varies between $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$7 and $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$8, with factor $\log(M_\mathrm{BH}/M_\odot) \simeq 6.3$9 changes on day timescales. Adopting a flat $M \approx 2 \times 10^6\, M_\odot$0CDM cosmology with $M \approx 2 \times 10^6\, M_\odot$1, $M \approx 2 \times 10^6\, M_\odot$2, $M \approx 2 \times 10^6\, M_\odot$3, and a NED-computed luminosity distance $M \approx 2 \times 10^6\, M_\odot$4 Mpc sets the luminosity scale used throughout the study.

The authors’ inspection of $M \approx 2 \times 10^6\, M_\odot$5 versus absorption-corrected 0.3–10 keV flux suggests two regimes. Below $M \approx 2 \times 10^6\, M_\odot$6, the spectral slope spans a wide range, including very hard states with $M \approx 2 \times 10^6\, M_\odot$7. Above that threshold, the source tends to be softer-when-brighter. The paper stresses that the dataset is heterogeneous, mixing instruments, off-axis pointings, and background levels, so a rigorous regression or correlation coefficient is not reported. Hardness-ratio time series are largely consistent with a constant within individual observations.

A set of representative fits illustrates the range of spectral states:

Epoch	Best-fit summary	$M \approx 2 \times 10^6\, M_\odot$8
XMM 0124710101 (2000-06-21)	$M \approx 2 \times 10^6\, M_\odot$9; $\sim 20$0 keV	$\sim 20$1
XMM 0652310901 (2010-12-05)	$\sim 20$2; $\sim 20$3 keV	$\sim 20$4
XMM 0691610301 (2012-06-04)	$\sim 20$5; $\sim 20$6; $\sim 20$7 keV	$\sim 20$8
Chandra 22648 (2020-03-03)	$\sim 20$9; $7^{\prime\prime}\times 5^{\prime\prime}$0	$7^{\prime\prime}\times 5^{\prime\prime}$1
Chandra 24853 (2020-11-03)	$7^{\prime\prime}\times 5^{\prime\prime}$2; $7^{\prime\prime}\times 5^{\prime\prime}$3	$7^{\prime\prime}\times 5^{\prime\prime}$4

The spectral behavior is one of the main reasons J1257 resists straightforward classification. It shows both a hard X-ray coronal component and only a weak softer-when-brighter signature at the cadence available during flares, which sets it apart from the most archetypal soft QPE sources.

5. Physical interpretations and characteristic scales

The paper considers two principal interpretations: a low-frequency QPO and QPE-like activity in an AGN with unusual spectral properties (Imbrogno et al., 3 Oct 2025).

As a low-frequency QPO candidate, J1257 is supported by the broad Lomb–Scargle peak at $7^{\prime\prime}\times 5^{\prime\prime}$5 Hz, the $7^{\prime\prime}\times 5^{\prime\prime}$6 ks period inferred from sinusoidal fitting, and the persistent modulation seen throughout 2020 in phase-folded profiles. The caveat is that standard mass–frequency scaling makes such a low frequency difficult to reconcile with $7^{\prime\prime}\times 5^{\prime\prime}$7, unless Lense–Thirring precession at relatively large radii is invoked. The observed frequency is therefore unusually low compared with other SMBH QPOs at similar masses.

As a QPE interpretation, the C1 sequence can be described as three close-to-symmetric flares with similar separations, $7^{\prime\prime}\times 5^{\prime\prime}$8 and $7^{\prime\prime}\times 5^{\prime\prime}$9 ks, followed by a shorter, higher-amplitude burst. These can be modeled either as a sequence of Gaussians (model 1) or as a sinusoidal modulation of quiescent emission plus a Gaussian QPE (model 2), motivated by the well-studied GSN 069 behavior where a sinusoidal “QPO-like” modulation accompanies QPEs. The inferred recurrence and duration lie within the known QPE timing plane: $\ge 95\%$0–51.3 ks and $\ge 95\%$1–19.5 ks.

The difficulty for a bona fide QPE classification is spectral. QPEs are characterized by soft thermal spectra with $\ge 95\%$2–200 eV and strong energy dependence, whereas J1257’s flares are comparatively achromatic and show only modestly larger amplitude in the soft band. The clear evidence for a hard X-ray corona, unlike classic QPE hosts, argues against a standard QPE interpretation unless Comptonization by the corona redistributes the flare photons to higher energies. The available data did not permit measurement of soft–hard time lags that could test this scenario.

The paper also connects the observed timescales to accretion-disk scales. For $\ge 95\%$3, the gravitational timescale

$\ge 95\%$4

is $\ge 95\%$5 s, so the $\ge 95\%$6 ks modulation corresponds to $\ge 95\%$7. If the variability is associated with Keplerian motion, the characteristic radius follows from

$\ge 95\%$8

For a baseline timescale near 25.4 ks, and depending on whether one or two impacts occur per orbit in an EMRI-based QPE scenario, the implied semi-major axis is in the range $\ge 95\%$9–88 $\le 4$0, where $\le 4$1. This suggests moderately inner-disc radii for a low-mass SMBH. A plausible implication is that, if the QPE-like interpretation were confirmed, J1257 would occupy the portion of parameter space relevant to proposed EMRI/QPE connections.

The SMBH mass also implies an Eddington luminosity

$\le 4$2

for $\le 4$3. The paper does not provide a bolometric correction or Eddington ratio. It notes only that adopting a typical Seyfert bolometric correction to the X-ray, $\le 4$4–30, would place J1257 at $\le 4$5 few $\le 4$6, but that this is only indicative and not derived in the work.

6. Context, limitations, and observational outlook

J1257 adds to the small but growing set of low-mass SMBHs with extreme short-term X-ray variability (Imbrogno et al., 3 Oct 2025). Its timescales overlap those of QPEs such as GSN 069 and RX J1301.9+2747, but its hard X-ray corona and muted energy dependence during flares distinguish it from bona fide soft QPEs. Relative to the robust AGN QPO in RE J1034+396, with period $\le 4$7 h, and to other reported QPOs, J1257 has a lower frequency for its mass and a much less secure significance because only sparse cycles are available. The source 2XMM J123103.2+110648 has been discussed in the literature as possibly related to long-duty-cycle QPE-like behavior, and J1257’s phenomenology is arguably closer to such intermediate cases than to classic QPEs.

Several important analyses are explicitly absent from the current dataset. The paper does not report fractional variability measurements, rms spectra, or interband time lags. Hardness ratios were examined in multiple epochs and were generally consistent with a constant within individual observations, including during the flare in ObsID 0691610201. The lack of longer uninterrupted coverage prevents robust PSD modeling to account for red noise at low frequencies, epoch-folding significance tests, and Gaussian-process-based period searches.

The main conclusion is therefore classificatory restraint rather than source-type closure. J1257 is a low-mass SMBH showing repeated X-ray variability at $\le 4$8–25 ks, including quasi-coherent modulations and flare-like events, with hints of two flux states and a softer-when-brighter trend at higher flux, but the present data do not decide between a particularly low-frequency QPO and QPE-like activity modified by a persistent corona.

The observational path forward is correspondingly specific. Longer, uninterrupted X-ray observations are needed to secure multiple consecutive cycles, quantify the significance with proper red-noise modeling, and perform time-resolved spectroscopy and lag analysis. \textit{XMM-Newton} and \textit{Chandra} with $\le 4$9 ks continuous exposures are identified as well-suited for this purpose; NICER could provide high-throughput soft-band timing, and future missions such as Athena could refine the spectral–timing diagnostics. If QPE-like behavior is confirmed, the implied EMRI orbital scales of tens of $\simeq 20$00 would make J1257 a source of interest for LISA as a potential electromagnetic counterpart in the EMRI/QPE connection.