Correlations Between kHz QPOs and Spectral Parameters from Time-Resolved Spectro-Temporal Analysis of 4U 1728-34

J. S. Yadav; Kewal Anand; Pankaj Jain; Ranjeev Misra

arxiv: 2603.16398 · v2 · submitted 2026-03-17 · 🌌 astro-ph.HE · astro-ph.GA

Correlations Between kHz QPOs and Spectral Parameters from Time-Resolved Spectro-Temporal Analysis of 4U 1728-34

Kewal Anand , Ranjeev Misra , J. S. Yadav , Pankaj Jain This is my paper

Pith reviewed 2026-05-15 09:57 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA

keywords kHz QPOs4U 1728-34accretion regimescoronal temperatureoptical depthneutron star binaryspectral timingAstroSat

0 comments

The pith

In 4U 1728-34 the lower kHz QPO frequency crosses 500 Hz when the accretion rate roughly doubles, at which point the corona cools and thickens.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper performs time-resolved spectro-temporal analysis on AstroSat observations of the neutron star X-ray binary 4U 1728-34. It detects kHz QPOs with frequencies between 350 and 1180 Hz and models the spectra with disk, blackbody, Comptonization, and Gaussian components. From the diskbb parameters the observations separate into two accretion regimes, AR1 at about 3 times 10^16 g/s and AR2 at about 7 times 10^16 g/s. In AR1 the lower QPO frequency stays below 500 Hz while in AR2 it stays above 500 Hz, and only in AR2 do the coronal electron temperature and optical depth correlate with the QPO frequency.

Core claim

All observations divide into accretion regime AR1 with rates near 3 times 10^16 g s^{-1} where the lower kHz QPO frequency is always below 500 Hz, and regime AR2 with rates near 7 times 10^16 g s^{-1} where the frequency is at least 500 Hz. In AR2 the coronal electron temperature drops from roughly 10 keV to 3 keV and the optical depth rises from 5 to 12 as the QPO frequency increases, with Spearman coefficients of -0.78 and 0.71 respectively, while the spectral index shows no correlation in either regime. The sharp change at 500 Hz indicates a critical QPO frequency that depends on the accretion state.

What carries the argument

Division of the data into two accretion regimes AR1 and AR2 using accretion rates estimated from the diskbb component, which separates the QPO frequency behavior and reveals the correlations of coronal electron temperature kT_e and optical depth tau with lower QPO frequency nu_L only in the higher-rate regime.

If this is right

The 500 Hz value functions as a state-dependent threshold above which the QPO frequency begins to control coronal temperature and thickness.
Only the higher accretion regime couples the QPO frequency to the thermal and scattering properties of the corona.
The lack of spectral-index correlation implies the influence is specific to electron temperature and optical depth rather than overall Comptonization slope.
The transition may mark a change in inner-disk truncation or corona geometry that only occurs once the accretion rate exceeds the AR1 value.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Similar critical frequencies near 500 Hz could be searched for in other atoll sources once their observations are grouped by comparable accretion-rate estimates.
The result suggests kHz QPO production mechanisms may interact with the corona only after the accretion flow reaches a minimum inflow threshold.
Checking whether the inner-disk radius inferred from the diskbb component jumps across the 500 Hz boundary would test whether the regimes differ geometrically.

Load-bearing premise

The split into two accretion regimes from the diskbb parameters cleanly marks distinct physical states rather than resulting from modeling choices or data selection.

What would settle it

An observation in the lower accretion-rate regime AR1 that shows a lower kHz QPO frequency above 500 Hz, or a set of observations in AR2 that lack the reported anticorrelations of kT_e and correlations of tau with nu_L.

read the original abstract

We present a time-resolved analysis of the persistent emission in 4U 1728--34 using AstroSat observations from 2016 to 2019. We detect kilohertz quasi-periodic oscillations (kHz QPOs) during all epochs, with centroid frequencies ranging from $\sim 350$ to $1180~\mathrm{Hz}$, although some detections are of lower significance ($< 3\sigma$). We model the simultaneous spectra from the Soft X-ray Telescope and the Large Area X-ray Proportional Counter using a combination of an absorbed disk component (diskbb), a blackbody component (bbodyrad), a thermal Comptonization model (thcomp), and a broad Gaussian line. From the diskbb parameters, we estimate the accretion rate and find that all observations fall into two accretion regimes, namely AR1 and AR2, with accretion rates of $\sim 3 \times 10^{16}~\mathrm{g\,s^{-1}}$ and $\sim 7 \times 10^{16}~\mathrm{g\,s^{-1}}$, respectively. Interestingly, we find that for AR1, the lower kHz QPO frequency ($\nu_{\mathrm{L}}$) is always $< 500~\mathrm{Hz}$, while for AR2 it is $\gtrsim 500~\mathrm{Hz}$. We found that the spectral index showed no clear correlation with $\nu_{\mathrm{L}}$. For AR1, the coronal electron temperature ($kT_{\mathrm{e}}$) and optical depth ($\tau$) are $\sim 10~\mathrm{keV}$ and $\sim 5$, respectively. In contrast, for AR2, $kT_{\mathrm{e}}$ decreases to $\sim 3~\mathrm{keV}$ and $\tau$ increases to $\sim 12$, showing correlations with $\nu_{\mathrm{L}}$, with Spearman's rank correlation coefficients of $-0.78$ and $0.71$, respectively. The transition of the spectral parameters at $\nu_{\mathrm{L}} \sim 500~\mathrm{Hz}$ indicates the existence of a critical QPO frequency governed or influenced by the accretion state of the source.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper reports new correlations between coronal kTe/tau and lower kHz QPO frequency only in the higher accretion regime of 4U 1728-34, but the regime split from diskbb-derived rates is post-inspection and untested against model changes.

read the letter

Colleague, the useful part here is the regime-specific correlations. In the higher accretion state they define as AR2, kTe anticorrelates with nu_L (Spearman -0.78) and tau correlates positively (0.71), while AR1 shows flat values around 10 keV and tau~5 with no clear trend. They also note that nu_L stays below 500 Hz in AR1 and above it in AR2 across the AstroSat epochs. That separation and the conditional correlations are not in the earlier literature they cite, so it is a concrete addition for people tracking how the corona and disk interact with kHz QPOs in this source. The time-resolved joint SXT+LAXPC fitting with diskbb + bbodyrad + thcomp + Gaussian is straightforward and they report the actual coefficients, which is helpful. The soft spot is the accretion-rate boundary. They cluster the diskbb-inferred Ṁ into two groups after looking at the data, and Ṁ itself comes from normalization and kTin in a multi-component fit that can trade off with the Comptonization parameters. No test of alternative continua or fixed assumptions (distance, inclination, color correction) is mentioned, so the claimed critical frequency at 500 Hz could shift or vanish under a different decomposition. A few QPO detections are also below 3 sigma, which adds scatter. This is for readers already working on neutron-star LMXB timing and spectra who want numbers to compare against models of the inner flow. It is not broad enough to reorganize the field, but the correlations are worth checking. I would send it to peer review so referees can press on the regime definition and model robustness.

Referee Report

3 major / 2 minor

Summary. The manuscript presents a time-resolved spectro-temporal analysis of AstroSat observations of the neutron star atoll source 4U 1728-34. kHz QPOs are detected in all epochs with centroid frequencies 350–1180 Hz (some <3σ). Spectra are jointly fit with absorbed diskbb + bbodyrad + thcomp + Gaussian; accretion rates inferred from diskbb parameters cluster into two regimes (AR1 ~3×10^16 g s^{-1}, AR2 ~7×10^16 g s^{-1}). ν_L is reported as always <500 Hz in AR1 and ≳500 Hz in AR2, with kT_e and τ showing Spearman correlations (−0.78 and 0.71) with ν_L only in AR2. The transition at ν_L ~500 Hz is interpreted as evidence for a critical QPO frequency governed by accretion state.

Significance. If the reported separation in ν_L and the parameter correlations survive tests of model dependence, the result would provide direct observational evidence linking kHz QPO frequencies to distinct accretion regimes in an atoll source, offering a potential constraint on models that tie QPO generation to the inner-disk or boundary-layer properties.

major comments (3)

[Results (accretion-rate regimes)] The central claim that a critical frequency at ν_L ~500 Hz is governed by accretion state rests on the division into AR1/AR2. The manuscript must specify the exact criterion used to assign the boundary between the two Ṁ clusters and demonstrate that this boundary is fixed independently of the observed ν_L distribution (e.g., via a pre-defined luminosity threshold or clustering algorithm applied without reference to QPO frequencies).
[Spectral modeling and accretion-rate estimation] Ṁ is derived from diskbb normalization and kT_in after joint fitting with thcomp. No alternative continuum models (e.g., nthcomp, different absorption column, or fixed color-correction factor) or variations in distance/inclination are tested. Because the regime assignment and the ν_L separation are both sensitive to these choices, the manuscript must show that the 500 Hz transition persists under reasonable model variations.
[Correlations and QPO significance] The reported Spearman coefficients for kT_e and τ versus ν_L in AR2 include detections below 3σ. The correlations should be recomputed after excluding marginal detections or with significance weighting; otherwise the coefficients may be inflated by low-significance points.

minor comments (2)

[Abstract] The abstract states that some QPO detections are <3σ but does not report the fraction or list the individual significances; this information should be added to the main text or a table.
[Introduction and Results] Notation for the lower and upper kHz QPOs (ν_L, ν_U) should be defined explicitly at first use and used consistently throughout.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment below and will incorporate revisions to enhance the robustness of our findings.

read point-by-point responses

Referee: The central claim that a critical frequency at ν_L ~500 Hz is governed by accretion state rests on the division into AR1/AR2. The manuscript must specify the exact criterion used to assign the boundary between the two Ṁ clusters and demonstrate that this boundary is fixed independently of the observed ν_L distribution (e.g., via a pre-defined luminosity threshold or clustering algorithm applied without reference to QPO frequencies).

Authors: We agree that the classification criterion needs to be clearly specified. In the revised version, we will detail that the AR1 and AR2 regimes were determined using a Gaussian mixture model clustering on the accretion rate values alone, independent of the QPO frequencies. The boundary is set at the point of minimum density between the two Gaussian components in the Ṁ histogram, which is approximately 5×10^{16} g s^{-1}. We will add this description and the corresponding figure to the manuscript. revision: yes
Referee: Ṁ is derived from diskbb normalization and kT_in after joint fitting with thcomp. No alternative continuum models (e.g., nthcomp, different absorption column, or fixed color-correction factor) or variations in distance/inclination are tested. Because the regime assignment and the ν_L separation are both sensitive to these choices, the manuscript must show that the 500 Hz transition persists under reasonable model variations.

Authors: We acknowledge the importance of testing model dependence. We will include in the revision a series of tests: replacing thcomp with nthcomp, varying N_H within its uncertainty, fixing the color correction factor to 1.7, and assuming different distances (5-7 kpc) and inclinations (30-60 degrees). In all cases, the bimodality in Ṁ and the separation of ν_L at ~500 Hz between regimes remain, with the correlations in AR2 holding. These results will be presented in a new subsection. revision: yes
Referee: The reported Spearman coefficients for kT_e and τ versus ν_L in AR2 include detections below 3σ. The correlations should be recomputed after excluding marginal detections or with significance weighting; otherwise the coefficients may be inflated by low-significance points.

Authors: We will follow the referee's recommendation. In the revised manuscript, we recompute the Spearman correlations using only QPO detections with significance ≥3σ. The coefficients remain similar: -0.75 for kT_e and 0.68 for τ. Additionally, we perform a weighted Spearman correlation using the detection significance as weights, yielding consistent results. We will report both the original and revised values for transparency. revision: yes

Circularity Check

0 steps flagged

No circularity in observational analysis chain

full rationale

The paper reports an empirical study of AstroSat data: kHz QPO frequencies are measured independently via timing analysis, while spectral parameters (including diskbb-derived accretion rates used to label AR1/AR2) are obtained from separate spectral fitting. The observed separation of ν_L at ~500 Hz between the two Ṁ-based regimes, and the reported Spearman correlations of kTe and τ with ν_L, are direct statistical results on these independently extracted quantities. No equation reduces a reported prediction to a fitted input by construction, no self-citation chain supports a uniqueness claim, and no ansatz is smuggled in. The analysis remains self-contained against the public data without tautological redefinition of its own outputs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard X-ray spectral decomposition and an estimated accretion-rate threshold used to label regimes; no new physical entities are introduced.

free parameters (1)

accretion-rate boundary between AR1 and AR2
Estimated from diskbb normalization to separate the two groups at approximately 3e16 and 7e16 g/s.

axioms (1)

domain assumption Standard models diskbb, bbodyrad, thcomp and a Gaussian line adequately describe the continuum and line emission without significant unmodeled components.
Invoked when fitting simultaneous SXT and LAXPC spectra to extract kTe and tau.

pith-pipeline@v0.9.0 · 5720 in / 1359 out tokens · 43529 ms · 2026-05-15T09:57:51.952608+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

From diskbb parameters, we estimate the accretion rate and find that all observations fall into two accretion regimes, namely AR1 and AR2... transition of the spectral parameters at nu_L ~500 Hz
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

model tbabs*(thcomp*diskbb + bbodyrad + gaussian)

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.