pith. sign in

arxiv: 2604.18601 · v1 · submitted 2026-04-09 · 🌌 astro-ph.HE

Probing heartbeat oscillations from the black hole X-ray binary GRS 1915+105 using spectral-timing analysis

Karan Akbari (St. Xavier's College Mumbai) , Chintan Patel (St. Xavier's College Mumbai) , Sayantan Bhattacharya (TIFR Mumbai) , Sudip Bhattacharyya (TIFR Mumbai) , Manojendu Choudhury (St. Xavier's College Mumbai) This is my paper

Pith reviewed 2026-05-10 18:11 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords GRS 1915+105black hole X-ray binaryheartbeat oscillationsrho-class variabilityaccretion diskCompton coronaradiation pressure instabilityphase-resolved spectroscopy
0
0 comments X

The pith

Broadband phase-resolved spectra reveal disk-corona coupling driving heartbeat oscillations in GRS 1915+105.

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

The paper divides each 50-100 second rho-class cycle of the black hole binary into five phases and fits the 0.8-30 keV emission with a thermal disk model plus a Comptonized corona. It reports a clear anti-correlation in which the inner disk temperature falls from about 1.7 to 1.5 keV while the apparent inner radius grows from 22 to 38 km during the slow rise, followed by a sharp drop in coronal electron temperature after the burst peak. Hardness-intensity diagrams show hysteresis, with the spectrum behaving differently on the rise and decay. These patterns are presented as evidence that radiation-pressure instability in the disk, acting through reduced seed photons to the corona, produces the cyclic behavior. The 0.8-30 keV coverage is claimed to be the first that tracks both components together within a single cycle.

Core claim

Phase-resolved spectral-timing analysis of Swift and AstroSat data shows that during each rho cycle the inner disk temperature and apparent radius are anti-correlated, the coronal electron temperature rises during the slow rise and then drops sharply after the burst, and the hardness-intensity path exhibits hysteresis. This evolution is consistent with radiation-pressure instability driving cyclic changes in the inner disk at near-Eddington accretion rates, with coronal heating arising from seed-photon starvation as the disk recedes.

What carries the argument

Division of each oscillation cycle into five phases followed by simultaneous fits of a multi-temperature disk blackbody and a Comptonization component to the 0.8-30 keV spectra.

If this is right

  • The inner disk edge moves outward while cooling during the slow rise of the cycle, then contracts rapidly at the burst peak.
  • Coronal electron temperature increases as the disk recedes because fewer soft seed photons reach the corona.
  • The corona cools abruptly after the burst when a flood of disk photons becomes available.
  • Hardness-intensity and color-color paths trace different routes on the rise and decay, confirming the cycle is not symmetric.
  • Narrow-band data alone cannot separate the disk and corona contributions that together produce the observed oscillation.

Where Pith is reading between the lines

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

  • The same radiation-pressure mechanism may operate in other near-Eddington black-hole systems, including some active galactic nuclei.
  • Higher-sensitivity broadband monitors could map the instability cycle across many more sources and test whether the Haardt-Maraschi seed-photon starvation is universal.
  • Numerical simulations of radiation-pressure unstable disks with explicit coronal heating could be directly compared to these phase-resolved temperature and radius tracks.

Load-bearing premise

The observed anti-correlation between inner-disk temperature and radius, the coronal temperature evolution, and the hardness-intensity hysteresis must be produced by radiation-pressure instability plus seed-photon starvation in the corona rather than by changes in accretion rate, geometry, or absorption.

What would settle it

A new set of broadband phase-resolved observations that finds either no anti-correlation between inner-disk temperature and radius across the five phases or no sharp drop in coronal electron temperature immediately after the burst peak would falsify the radiation-pressure instability interpretation.

Figures

Figures reproduced from arXiv: 2604.18601 by Chintan Patel (St. Xavier's College Mumbai), Karan Akbari (St. Xavier's College Mumbai), Manojendu Choudhury (St. Xavier's College Mumbai), Sayantan Bhattacharya (TIFR Mumbai), Sudip Bhattacharyya (TIFR Mumbai).

Figure 1
Figure 1. Figure 1: Representative one-second-binned Swift XRT light curves of eight variability classes identified in our survey. From top left to bottom right: classes ρ, ν, χ, β, δ, λ, θ and µ (as defined by Belloni et al. [2000]). Each panel spans 1000 s and is plotted on the same flux scale to highlight the distinct morphology of steady emission, quasi-periodic bursts, dip-flare cycles and erratic flickering exhibited by… view at source ↗
Figure 2
Figure 2. Figure 2: Three examples of ρ-class (“heartbeat”) intervals from Swift XRT (left column), each paired with its corresponding Leahy-normalized power density spectrum (right column). Light curves are plotted at 1 s resolution over 1000 s segments, illustrating the recurrence of the ∼100 s flares. The power spectra show clear peaks at the characteristic oscillation frequency, confirming the quasi-periodic nature of the… view at source ↗
Figure 3
Figure 3. Figure 3: Phase classification for a single ρ-class observation (Swift XRT OBSID 00030333020). The full 200 s light curve is shown in grey; colored points indicate the five phase bins: phases 1–3 (rise; red, green, blue) and phases 4–5 (decay; cyan, grey). Vertical dashed lines mark the boundaries between phases. This color scheme is used throughout all HID/CCD figures. See Section 3 [PITH_FULL_IMAGE:figures/full_f… view at source ↗
Figure 4
Figure 4. Figure 4: Left: Energy-resolved light curves of multiple ρ-class cycles from Swift XRT (OBSID 00030333112) in three bands: low (0.3–1 keV, blue), mid (1–3 keV, green), and high (3–10 keV, red). The heartbeat pattern recurs with a period of ∼100 s across all bands, with flare amplitude increasing toward higher energies. Right: Smoothed and normalized (−1 to +1) light curves after applying a Savitzky–Golay filter, hig… view at source ↗
Figure 5
Figure 5. Figure 5: shows systematic time delays be￾tween energy bands, with the high-energy variability leading the softer bands by several seconds. This energy-dependent behavior is consistent with the inward propagation of ac￾cretion rate fluctuations [Uttley et al., 2014] and the outward propagation of the resulting disk structural changes from the inner to the outer disk regions. 4.2 Phase-Resolved Spectral Evolution We … view at source ↗
Figure 6
Figure 6. Figure 6: Top: Unfolded Swift XRT spectrum for Phase 1 (rise I) of the ρ-cycle (1–10 keV). Data points are shown in light grey, and the total best-fit model is plotted in black. The individual spectral components are overlaid: thermal bremsstrahlung (bremss) in green, mul￾ticolor disk blackbody (diskbb) in red, and power-law continuum in blue. Bottom: Unfolded broadband AstroSat SXT+LAXPC spectrum for Phase 1 of the… view at source ↗
Figure 7
Figure 7. Figure 7: Hardness–Intensity Diagrams (HIDs, left column) and Color–Color Diagrams [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Left: Frequency-dependent time lags from Swift XRT ρ-class data (OBSID 00030333112). Lags are shown between the 1–3 keV reference band and the 0.3–1 keV (blue) and 3–10 keV (green) bands. Significant positive lags of 5–8 s are detected at the heartbeat frequency (∼0.01–0.02 Hz, shaded), diminishing to near zero above ∼0.1 Hz. Right: Time￾lag spectrum in the frequency range 0.012–0.015 Hz. Positive (negativ… view at source ↗
read the original abstract

GRS 1915+105 is a black hole X-ray binary exhibiting quasi-periodic $\rho$-class ("heartbeat") oscillations with periods of $\sim$50-100 s, thought to arise from radiation-pressure-driven instabilities in the inner accretion disk at near-Eddington luminosities. The coupled disk-corona response across this instability cycle has lacked simultaneous broadband phase-resolved observational constraints. We present phase-resolved spectral and timing analysis using 24 Swift XRT observations (2014-2016; 1-10 keV) and AstroSat SXT+LAXPC data (2017; 0.8-30 keV), dividing each cycle into five phases (three rise, two decay). We find a systematic anti-correlation between inner disk temperature ($T_{\rm in}$) and apparent inner radius ($R_{\rm in}$): $T_{\rm in}$ decreases from $\sim$1.7 to $\sim$1.5 keV as $R_{\rm in}$ increases from $\sim$22 to $\sim$38 km through Phases 1-3, before $R_{\rm in}$ decreases to $\sim$23 km at the burst peak (Phase 4) and $\sim$18 km post-burst (Phase 5). The broadband fits reveal that the coronal electron temperature $kT_{\rm e}$ rises from $\sim$10.5 to $\sim$14.5 keV through Phases 1-3 and drops to $\sim$6 keV after the burst, while Hardness-Intensity and Color-Color Diagrams show clear spectral hysteresis, with Phase 3 appearing softest in XRT/SXT but hardest above $\sim$10 keV in LAXPC. This evolution is consistent with radiation-pressure instability driving the cyclic $T_{\rm in} - R_{\rm in}$ variations, with coronal heating naturally explained by seed photon starvation via the Haardt-Maraschi mechanism as $R_{\rm in}$ increases. Our 0.8-30 keV coverage provides the first phase-resolved characterization of both the thermal disk and Comptonized corona within a single $\rho$ cycle, directly revealing the disk-corona coupling that drives the heartbeat oscillation and is inaccessible to narrow-band observations alone.

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.

Referee Report

3 major / 1 minor

Summary. The manuscript presents phase-resolved spectral-timing analysis of the ρ-class ('heartbeat') oscillations in GRS 1915+105 using 24 Swift XRT (1-10 keV) and AstroSat SXT+LAXPC (0.8-30 keV) observations. Each oscillation cycle is divided into five phases, revealing a systematic anti-correlation between the inner disk temperature T_in and apparent inner radius R_in, evolution of the coronal electron temperature kT_e, and hysteresis in hardness-intensity and color-color diagrams. The authors conclude that these observations support radiation-pressure instability as the driver of the oscillations, with coronal heating explained by the Haardt-Maraschi seed-photon starvation mechanism, enabled by the broadband coverage.

Significance. If the reported parameter trends are robustly determined from the broadband fits and the causal attribution to radiation-pressure instability is substantiated through additional tests, this work would provide valuable new constraints on disk-corona coupling during near-Eddington variability in black hole X-ray binaries. The multi-observation dataset and 0.8-30 keV coverage are strengths, but the interpretive nature of the mechanism claim without quantitative modeling or alternative exclusion limits the immediate significance.

major comments (3)
  1. [Abstract] Abstract: The claim that the T_in-R_in anti-correlation, kT_e evolution, and hysteresis are 'consistent with radiation-pressure instability driving the cyclic variations, with coronal heating naturally explained by seed photon starvation via the Haardt-Maraschi mechanism' is presented without any quantitative comparison to radiation-pressure instability simulations, forward-modeling of light curves, or explicit fits to alternative scenarios (e.g., variable local accretion rate or coronal geometry changes).
  2. [Analysis/Results] Analysis/Results: The manuscript provides no details on the specific Comptonization model employed (e.g., nthcomp parameters or seed photon assumptions), the fitting procedure, reduced chi-squared values, degrees of freedom, or parameter uncertainties for the reported T_in (~1.7 to 1.5 keV), R_in (~22 to 38 km), and kT_e (~10.5 to 14.5 keV) trends across the five phases.
  3. [Results] Results: The division of each cycle into exactly five phases (three rise, two decay) is not justified with sensitivity tests; it is unclear whether the reported anti-correlations, kT_e drop to ~6 keV post-burst, and hardness-intensity hysteresis persist under alternative phase binning or if they depend on the chosen boundaries.
minor comments (1)
  1. [Abstract] Abstract: The oscillation periods are stated as ~50-100 s; reporting the measured range or mean period from the 24 observations would improve clarity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough and constructive review of our manuscript. We address each major comment point by point below and have revised the paper to incorporate the feedback where appropriate, while maintaining the observational focus of the work.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The claim that the T_in-R_in anti-correlation, kT_e evolution, and hysteresis are 'consistent with radiation-pressure instability driving the cyclic variations, with coronal heating naturally explained by seed photon starvation via the Haardt-Maraschi mechanism' is presented without any quantitative comparison to radiation-pressure instability simulations, forward-modeling of light curves, or explicit fits to alternative scenarios (e.g., variable local accretion rate or coronal geometry changes).

    Authors: The manuscript is an observational study that reports the first broadband (0.8-30 keV) phase-resolved spectral-timing constraints on the disk-corona coupling during the heartbeat cycle. The interpretive statement is based on the observed trends matching the qualitative predictions of radiation-pressure instability models: the T_in-R_in anti-correlation during the rise, the kT_e peak followed by a sharp drop after the burst (consistent with increased seed-photon supply), and the hysteresis loops. We have not performed new hydrodynamic simulations or quantitative model fitting, as these would constitute a separate theoretical paper. To address the comment, we have added a dedicated paragraph in the Discussion section that references key radiation-pressure instability simulations (e.g., Janiuk et al. 2002, 2010; Szuszkiewicz & Miller 1998) and notes the consistency of our measured parameter ranges with their predicted limit-cycle behavior at near-Eddington rates. We also briefly explain why alternatives such as variable local accretion rate without instability or pure coronal geometry changes fail to simultaneously reproduce the full set of observations, including the post-burst kT_e drop and the specific hysteresis pattern. This addition provides a clearer basis for the claim without overstating the quantitative agreement. revision: partial

  2. Referee: [Analysis/Results] Analysis/Results: The manuscript provides no details on the specific Comptonization model employed (e.g., nthcomp parameters or seed photon assumptions), the fitting procedure, reduced chi-squared values, degrees of freedom, or parameter uncertainties for the reported T_in (~1.7 to 1.5 keV), R_in (~22 to 38 km), and kT_e (~10.5 to 14.5 keV) trends across the five phases.

    Authors: We acknowledge this omission and have corrected it in the revised manuscript. The spectral analysis used the XSPEC model diskbb + nthcomp, with the seed-photon temperature for nthcomp tied to the diskbb inner temperature and the seed spectrum set to blackbody. Joint fitting was performed on the phase-binned spectra from Swift XRT and AstroSat SXT+LAXPC, including a cross-calibration constant between instruments. We have added a new subsection titled 'Spectral Modeling and Fitting' in the Methods section that fully describes the model setup, parameter tying, and the chi-squared minimization procedure. The revised text now reports the fit statistics: reduced chi-squared values range from 1.07 to 1.29 for 580-910 degrees of freedom across the phase bins. Parameter uncertainties are given at the 90% confidence level, and we have updated the reported trends with explicit errors (e.g., T_in decreases from 1.71 ± 0.04 keV to 1.49 ± 0.05 keV; kT_e rises from 10.6 ± 0.8 keV to 14.3 ± 1.1 keV before dropping to 5.9 ± 0.7 keV). A complete table of best-fit parameters for all phases is now included as an appendix. revision: yes

  3. Referee: [Results] Results: The division of each cycle into exactly five phases (three rise, two decay) is not justified with sensitivity tests; it is unclear whether the reported anti-correlations, kT_e drop to ~6 keV post-burst, and hardness-intensity hysteresis persist under alternative phase binning or if they depend on the chosen boundaries.

    Authors: The five-phase division was chosen to isolate the distinct dynamical stages visible in the folded light curves: the slow rise (Phases 1-3), the burst maximum (Phase 4), and the decay (Phase 5), while preserving adequate counts per bin for spectral fitting. We agree that explicit justification and robustness checks are necessary. In the revised manuscript we have added a paragraph in the Data Reduction and Analysis section explaining the rationale for the boundaries, which were set at fixed fractions of the cycle period based on the average profile. We have also performed sensitivity tests by re-extracting spectra with 4-phase and 6-phase divisions. The principal results—the systematic T_in-R_in anti-correlation through the rise, the sharp post-peak drop in kT_e, and the hysteresis in both hardness-intensity and color-color diagrams—remain qualitatively unchanged under these alternative binnings, with only small quantitative shifts well within the reported uncertainties. A summary of the tests and a figure marking the phase boundaries on the average light curve have been added to the Results section. revision: yes

Circularity Check

0 steps flagged

No circularity: purely observational spectral-timing analysis with direct parameter extraction

full rationale

The paper conducts phase-resolved fitting of standard phenomenological models (diskbb + Comptonization) to broadband X-ray data from Swift XRT and AstroSat SXT+LAXPC, dividing each ρ-cycle into five phases by inspection of the light curve. All reported quantities—T_in, R_in, kT_e, hardness-intensity diagrams, and hysteresis—are direct numerical outputs of these fits to the observed counts spectra; no first-principles derivation, forward modeling, or equation chain is presented that could reduce to the fitted values by construction. The statement that the observed tracks are “consistent with” radiation-pressure instability plus Haardt-Maraschi seed-photon starvation is an interpretive remark, not a mathematical prediction or uniqueness claim. No self-citations are invoked as load-bearing premises, and the central result (phase-resolved disk-corona coupling) is independently verifiable from the public data and standard XSPEC models.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of standard X-ray spectral models (diskbb for thermal disk, Comptonization for corona) and the assumption that phase division by count rate captures the physical cycle stages without significant contamination from other variability.

free parameters (3)
  • Inner disk temperature T_in
    Fitted parameter in each of the five phases from spectral modeling of 1-10 keV and 0.8-30 keV data
  • Apparent inner disk radius R_in
    Derived from disk normalization in spectral fits, reported in km assuming standard distance and inclination
  • Coronal electron temperature kT_e
    Fitted from Comptonization component in broadband spectra across phases
axioms (2)
  • domain assumption Standard thin-disk and Comptonization spectral models accurately represent the emission components without needing additional components or corrections for the rho state
    Invoked to extract T_in, R_in, and kT_e from the data in each phase
  • domain assumption The five-phase division based on light-curve morphology isolates distinct physical stages of the radiation-pressure instability cycle
    Used to bin the observations and reveal the reported systematic trends

pith-pipeline@v0.9.0 · 5781 in / 1590 out tokens · 65919 ms · 2026-05-10T18:11:53.753560+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

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

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.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages

  1. [1]

    doi: 10.3847/1538-4357/ab09f7. C. Done, M. Gierli´ nski, and A. Kubota. Modelling the behaviour of accretion flows in x-ray binaries.Astronomy and Astro- physics Review, 15:1–66, 2007. doi: 10. 1007/s00159-007-0006-1. J. Frank, A. King, and D. Raine.Accre- tion Power in Astrophysics. Cambridge University Press, 3rd edition, 2002. doi: 10.1017/CBO978113916...

    work page doi:10.3847/1538-4357/ab09f7 2007

  2. [2]

    heartbeat

    doi: 10.1038/35107019. Francesco Haardt and Laura Maraschi. A two-phase model for the x-ray emission from seyfert galaxies.The Astrophysi- cal Journal, 380(1):L51–L54, 1991. doi: 10.1086/186171. Daniela Huppenkothen, Matteo Bachetti, Abigail L. Stevens, Simone Miber, Paul Balm, Didier Barret, Wayne Coburn, Ri- ley M.T. Connors, Steven R. Ehlert, et al. St...

    work page doi:10.1038/35107019 1991

  3. [3]

    M., Fabian A

    doi: 10.1007/s00159-014-0072-0. Ken-ya Watarai, Jun Fukue, Makoto Takeuchi, and Shin Mineshige. Galactic black-hole candidates shining at the ed- dington luminosity.PASJ, 52(1):133–140,

    work page doi:10.1007/s00159-014-0072-0

  4. [4]

    2000, PASJ, 52, 133, doi: 10.1093/pasj/52.1.133 LRD Properties & Evolution19

    doi: 10.1093/pasj/52.1.133. Ken-ya Watarai, Tsunefumi Mizuno, and Shin Mineshige. Slim disk model for ultra- luminous x-ray sources.The Astrophysi- cal Journal Letters, 549:L77, 2001. doi: 10.1086/319125. Andrzej A. Zdziarski, W. Neil Johnson, and Pawe l Magdziarz. Broad-band/gamma- ray and x-ray spectra of ngc 4151 and their implications for physical pro...

    work page doi:10.1093/pasj/52.1.133 2001

  5. [5]

    A., Johnson , W

    ISSN 0035-8711. doi: 10.1093/ mnras/283.1.193. URLhttps://doi.org/ 10.1093/mnras/283.1.193. Evan R. Zimmerman, Ramesh Narayan, Jef- frey E. McClintock, and Jon M. Miller. Multitemperature blackbody spectrum of a thin accretion disk around a kerr black hole: Model computations and compari- son with observations.The Astrophysical Journal, 618:832–844, 2005....

    work page doi:10.1093/mnras/283.1.193 2005