pith. sign in

arxiv: 2604.22643 · v1 · submitted 2026-04-24 · 🌌 astro-ph.HE

Radiative feedbacks as drivers for quasi-periodic-oscillation activity in black-hole X-ray binaries

Apostolos Mastichiadis This is my paper

Pith reviewed 2026-05-08 10:14 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords black hole X-ray binariestype-C QPOsquasi-periodic oscillationsradiative feedbackaccretion diskcoronalimit cycle oscillationsinverse Compton
0
0 comments X

The pith

Radiative feedback from reprocessed disk photons drives limit-cycle oscillations that match type-C QPOs in black-hole X-ray binaries.

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

The paper investigates whether radiative coupling in the disk-corona system of black-hole X-ray binaries can explain type-C quasi-periodic oscillations through limit-cycle behavior. It uses a one-zone model to solve the time-dependent equations for electrons and photons, finding that limit cycles arise when cooling of coronal electrons is dominated by soft photons that the disk reprocesses from the corona's hard X-rays. These cycles produce variations whose frequencies depend on the coronal radius and energization timescale, while the spectra are power laws up to 10-100 keV. For a range of parameters, the oscillations reproduce important observed properties of type-C QPOs, confirming and extending earlier simplified models, and the scale-free nature allows application to active galactic nuclei.

Core claim

When electron cooling is dominated by soft photons reprocessed in the accretion disk, the disk-corona system undergoes limit-cycle oscillations. The model solves the time-dependent kinetic equations for electrons and photons in a one-zone radiation framework, with electrons energized by an unspecified process and cooling through inverse Compton scattering of soft photons from the disk and from disk reprocessing of coronal hard radiation. For a subset of model parameters, these oscillations reproduce key properties of type-C QPOs observed in BHXRBs, with frequencies depending on coronal radius and energization timescale, and X-ray spectra described by power laws up to ~10-100 keV.

What carries the argument

The radiative feedback mechanism where hard X-rays produced in the corona are reprocessed by the accretion disk into soft photons that cool the coronal electrons; this feedback is the driver for the limit-cycle oscillations when it dominates the cooling.

Load-bearing premise

The dominance of disk-reprocessed soft photons in the cooling of coronal electrons, along with treating the energization as an unspecified constant timescale, is required for the limit cycles to develop.

What would settle it

If type-C QPOs are observed in spectral states where the fraction of soft photons from disk reprocessing is low compared to direct disk emission, or if the QPO frequency does not vary with estimated coronal size as predicted.

Figures

Figures reproduced from arXiv: 2604.22643 by Apostolos Mastichiadis.

Figure 1
Figure 1. Figure 1: Schematic representation of the feedback loop between electrons and photons in the coronal environment. For artistic reasons the corona is depicted as a cloud, while in the paper its shape is assumed to be spherical. (The figure is AI generated.) Electrons are assumed to be energized inside the corona by some (unspecified) mechanism and, at the same time, lose energy on soft photons. These are not independ… view at source ↗
Figure 2
Figure 2. Figure 2: Light curves of ICS compactness in the 2-200 keV band for various values of the parameter ℓds. All other parameters have been kept constant and have the values Rc = 109 cm, M˙ in j = 3.2 × 10−7M⊙/yr, t˜ en = 2.5, tesc/ten = 103 , Te f f = 4 × 106 K, α = 0.1. The colors of the lines correspond to ℓds = 10−6 , 10−5 , 10−4 , 10−3 , 10−2 and 10−1 from red to violet. For the sake of clarity each curve is displa… view at source ↗
Figure 3
Figure 3. Figure 3: Average photon spectra for the same cases shown in the previous figure. The color codes have been kept the same. x dℓ dx = x 2n˙γ(x, t), where x = ϵ/mec 2 . This is equivalent to ϵ dLγ dϵ where Lγ is the spectral luminosity. As one can observe, low values of ℓds produce harder ICS spectra; however, as ℓds increases, the disk photons hinder elec￾tron acceleration, and as a result both electrons is lower and… view at source ↗
Figure 4
Figure 4. Figure 4: Plot of a set of light curves all sharing the same value of Qe,in j (Min j = 16 × 10−8M⊙/yr) but having different energization timescales. Here t˜ en = 1.25 (violet), 2.5 (dark blue), 5 (light blue), 10 (green), 20 (yellow), 40 (orange) and 80 (red), all given in units of tcr. The other parameters are Rc = 109 cm, tesc/ten = 103 , α = 0.1, Te f f = 4 × 106 K, ℓds = 10−6 . For the sake of clarity each curve… view at source ↗
Figure 5
Figure 5. Figure 5: Plot of the average multiwavelength spectrum for the case shown in the previous figure keeping the color codes the same. – The fractional rms frms is of the order of tens of percent. However, these values must be considered strictly as upper limits, especially in those cases that are characterized by strong to intermediate damping. Generally, frms and T˜ 10 have an erratic behavior that will be discussed f… view at source ↗
Figure 6
Figure 6. Figure 6: Heat map of the parameter T˜ 10 as a function of t˜ en and M˙ in j for Rc = 109 cm. Other parameters are α = 0.1 and Te f f = 4 × 106 K. Dark regions correspond to light curves with strong decay view at source ↗
Figure 7
Figure 7. Figure 7: Plot of the QPO frequency νQPO versus the energization timescale t˜ en for various values of M˙ in j. Values of M˙ in j, in units of 10−8M⊙/yr, are 64 (violet), 32 (dark blue), 16 (light blue), 8 (green), 4 (yellow), 2 (orange) and 1 (red). Other parameters are Rc = 109 cm, α = 0.1 and Te f f = 4 × 106 K. The black line indicates the relation νQPO ∝ t −1/2 en given by Eq. 21. to the X-ray luminosity but it… view at source ↗
Figure 10
Figure 10. Figure 10: Heat map of the time lag ∆Tsh (in units of tcr) as a function of t˜ en and M˙ in j for Rc = 109 cm. Other parameters are α = 0.1 and Te f f = 4 × 106K. Violet color indicates negative lags, i.e. hard photons come before the soft ones in the prescribed energy bands. ceed values of around 3 irrespective of M˙ in j. Note that the high M˙ in j cases cannot reach low νQPO values due to damping. As an example o… view at source ↗
Figure 11
Figure 11. Figure 11: Plot of the phase lag ∆ϕ versus the QPO frequency νQPO. The violet line corresponds to M˙ in j = 8×10−8M⊙/yr for Rc between 109.3 cm and Rc = 107.95 cm, the dark blue line corresponds to M˙ in j = 16 × 10−8M⊙/yr for Rc between 109.6 cm and 108.25 cm, while the light blue line corresponds to M˙ in j = 32 × 10−8M⊙/yr for Rc between 109.9 cm to 108.85 cm. In all cases the radius Rc is decreasing from left to… view at source ↗
Figure 12
Figure 12. Figure 12: Plot of the fractional RMS versus the QPO frequency νQPO for the cases depicted in view at source ↗
read the original abstract

Black-hole X-ray binaries (BHXRBs) in the hard and hard-intermediate spectral states commonly exhibit prominent type-C quasi-periodic oscillations (QPOs) in their X-ray power spectra. Despite extensive observational and theoretical efforts, the physical mechanism responsible for these oscillations has not yet been firmly established. The disk-corona system in BHXRBs is radiatively coupled, as hard X-ray emission from the corona can be reprocessed by the accretion disk and re-emitted as soft photons that contribute to cooling the coronal electrons. Aim of the present study is to examine whether this feedback can give rise to limit cycles having the spectro-temporal properties of QPOs. We model the coronal emission using a one-zone radiation framework and solve the time-dependent kinetic equations for electrons and photons. Electrons are energized by some unspecified process and cool via inverse Compton scattering of soft photons originating from (i) the accretion disk and (ii) disk reprocessing of the hard radiation produced in the corona. When electron cooling is dominated by soft photons reprocessed in the accretion disk, the disk-corona system undergoes limit-cycle oscillations. For a subset of the model parameters, these oscillations reproduce key properties of type-C QPOs observed in BHXRBs. The oscillation frequency depends on the coronal radius and on the energization timescale, while the resulting X-ray spectra are well described by power laws extending up to energies of ~ 10-100 keV. These calculations confirm and extend earlier semi-analytical results obtained with simplified treatments. Owing to the scale-invariant nature of the model, the results can be readily extrapolated to other accreting systems, such as Active Galactic Nuclei.

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

2 major / 1 minor

Summary. The manuscript develops a one-zone time-dependent kinetic model for the coupled disk-corona system in black-hole X-ray binaries. It demonstrates that limit-cycle oscillations arise in the electron and photon distributions when cooling of coronal electrons is dominated by inverse Compton scattering of soft photons reprocessed by the accretion disk. For suitable choices of coronal radius and energization timescale, the resulting oscillation frequencies and X-ray spectra are shown to be consistent with observed type-C QPOs, extending previous semi-analytical treatments.

Significance. If the oscillatory solutions prove robust under variations in geometry and additional cooling channels, this work would establish radiative feedback as a plausible driver for QPO activity, providing a physically motivated mechanism that is scale-invariant and thus applicable to AGN as well. The confirmation of earlier results with a more complete kinetic treatment is a strength.

major comments (2)
  1. Abstract: The central result that limit-cycle oscillations appear 'when electron cooling is dominated by soft photons reprocessed in the accretion disk' relies on selecting parameters such that reprocessed photons dominate over direct disk photons; the manuscript does not show that this dominance emerges naturally from the model equations or realistic coronal geometry.
  2. Abstract: The oscillation frequency is explicitly stated to depend on the coronal radius and the energization timescale, both of which are free parameters tuned to match the observed QPO frequency range; this makes the reproduction of QPO properties more of a consistency check than a parameter-free prediction.
minor comments (1)
  1. The main text should clarify how the reprocessing fraction is calculated and whether it is self-consistently determined from the geometry or imposed externally.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments, which help clarify the scope and limitations of our one-zone model. We address each major comment below and have revised the manuscript accordingly where appropriate.

read point-by-point responses

  1. Referee: Abstract: The central result that limit-cycle oscillations appear 'when electron cooling is dominated by soft photons reprocessed in the accretion disk' relies on selecting parameters such that reprocessed photons dominate over direct disk photons; the manuscript does not show that this dominance emerges naturally from the model equations or realistic coronal geometry.

    Authors: We agree that the one-zone framework parameterizes the relative importance of reprocessed versus direct photons through the assumed coronal geometry and covering factor rather than deriving it self-consistently from the equations. The central result is therefore conditional on the regime in which reprocessed photons dominate cooling. We have revised the abstract to state this condition more explicitly and added a paragraph in the discussion section that outlines the geometric requirements (e.g., coronal height and solid angle) needed for this dominance in more realistic setups. This remains a limitation of the simplified model, but the kinetic equations themselves demonstrate that limit-cycle behavior arises robustly once the condition is satisfied. revision: partial

  2. Referee: Abstract: The oscillation frequency is explicitly stated to depend on the coronal radius and the energization timescale, both of which are free parameters tuned to match the observed QPO frequency range; this makes the reproduction of QPO properties more of a consistency check than a parameter-free prediction.

    Authors: The referee correctly notes that the frequency depends on these two parameters, which are adjusted within physically plausible ranges to reproduce observed QPO frequencies. The model therefore provides a consistency check rather than a parameter-free prediction of the frequency itself. Once the parameters are fixed, however, the resulting spectra, rms amplitudes, and phase lags are obtained without further tuning and match type-C QPO properties. We have updated the abstract and conclusions to emphasize that the work demonstrates a physically motivated mechanism capable of producing the observed spectro-temporal characteristics for reasonable parameter choices, extending earlier semi-analytical results. This is an inherent feature of exploratory models of this type. revision: yes

Circularity Check

1 steps flagged

QPO reproduction and oscillation frequency reduce to tuning of free parameters (coronal radius, energization timescale) that enforce reprocessed-photon dominance

specific steps
  1. fitted input called prediction [Abstract]
    "When electron cooling is dominated by soft photons reprocessed in the accretion disk, the disk-corona system undergoes limit-cycle oscillations. For a subset of the model parameters, these oscillations reproduce key properties of type-C QPOs observed in BHXRBs. The oscillation frequency depends on the coronal radius and on the energization timescale"

    The dominance condition required for oscillations is enforced by choosing values of coronal radius and energization timescale; the frequency itself is a direct function of those same parameters. Selecting the subset that matches observed QPO frequencies and spectra therefore reproduces the target properties by construction rather than predicting them from the kinetic equations without prior tuning.

full rationale

The paper solves time-dependent kinetic equations in a one-zone model but explicitly conditions the emergence of limit-cycle oscillations on a dominance regime that is controlled by adjustable parameters rather than arising self-consistently. The frequency is stated to depend directly on those same parameters, which are then chosen so that the solutions fall inside the observed QPO range. This makes the claimed reproduction of type-C QPO properties a consequence of parameter selection inside the model rather than an independent derivation. No self-citation chain or uniqueness theorem is invoked as load-bearing; the circularity is limited to the fitted-input pattern.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on a one-zone approximation, an unspecified energization process treated as a fixed timescale, and the requirement that reprocessed photons dominate cooling; these are introduced without independent derivation and directly control the existence of the reported oscillations.

free parameters (2)
  • coronal radius
    Sets the light-travel and cooling timescales that determine oscillation frequency
  • energization timescale
    Controls the rate at which electrons gain energy and therefore the cycle period
axioms (2)
  • domain assumption one-zone radiation framework with uniform electron and photon distributions
    Invoked to reduce the problem to ordinary differential equations for the corona
  • ad hoc to paper electrons are energized by an unspecified process whose net effect is a constant timescale
    Required to close the energy equation without specifying the microphysical heating mechanism

pith-pipeline@v0.9.0 · 5604 in / 1608 out tokens · 44315 ms · 2026-05-08T10:14:20.567550+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

    work page

  3. [3]

    2022, , 515, 2099

    Bellavita , C., Garc \' a , F., M \'e ndez , M., & Karpouzas , K. 2022, , 515, 2099

    work page 2022

  4. [4]

    Belloni , T. M. 2010, in Lecture Notes in Physics, Berlin Springer Verlag, ed. T. Belloni , Vol. 794, 53

    work page 2010

  5. [5]

    Blumenthal , G. R. & Gould , R. J. 1970, Reviews of Modern Physics, 42, 237

    work page 1970

  6. [6]

    2026, , 282, 22

    Cerruti , M., Rudolph , A., Petropoulou , M., et al. 2026, , 282, 22

    work page 2026

  7. [7]

    J., & Reimer , A

    Dimitrakoudis , S., Mastichiadis , A., Protheroe , R. J., & Reimer , A. 2012, , 546, A120

    work page 2012

  8. [8]

    2007, , 15, 1

    Done , C., Gierli \'n ski , M., & Kubota , A. 2007, , 15, 1

    work page 2007

  9. [9]

    2021, , 501, 3173

    Garc \' a , F., M \'e ndez , M., Karpouzas , K., et al. 2021, , 501, 3173

    work page 2021

  10. [10]

    & Spruit , H

    Giannios , D. & Spruit , H. C. 2004, , 427, 251

    work page 2004

  11. [11]

    2025, Large Language Model, [AI image generation tool]

    Google . 2025, Large Language Model, [AI image generation tool]

    work page 2025

  12. [12]

    & Maraschi , L

    Haardt , F. & Maraschi , L. 1991, , 380, L51

    work page 1991

  13. [13]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357

    work page 2020

  14. [14]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90

    work page 2007

  15. [15]

    Ingram , A., Done , C., & Fragile , P. C. 2009, , 397, L101

    work page 2009

  16. [16]

    Ingram , A. R. & Motta , S. E. 2019, , 85, 101524

    work page 2019

  17. [17]

    2021, , 503, 5522

    Karpouzas , K., M \'e ndez , M., Garc \' a , F., et al. 2021, , 503, 5522

    work page 2021

  18. [18]

    M., et al

    Karpouzas , K., M \'e ndez , M., Ribeiro , E. M., et al. 2020, , 492, 1399

    work page 2020

  19. [19]

    G., Rieger , F

    Kirk , J. G., Rieger , F. M., & Mastichiadis , A. 1998, , 333, 452

    work page 1998

  20. [20]

    & Misra , R

    Kumar , N. & Misra , R. 2014, , 445, 2818

    work page 2014

  21. [21]

    Lee , H. C. & Miller , G. S. 1998, , 299, 479

    work page 1998

  22. [22]

    Li , H., Kusunose , M., & Liang , E. P. 1996, , 460, L29

    work page 1996

  23. [23]

    & Kirk , J

    Mastichiadis , A. & Kirk , J. G. 1995, , 295, 613

    work page 1995

  24. [24]

    & Kirk , J

    Mastichiadis , A. & Kirk , J. G. 1997, , 320, 19

    work page 1997

  25. [25]

    Mastichiadis , A., Petropoulou , M., & Kylafis , N. D. 2022, , 662, A118

    work page 2022

  26. [26]

    H., Remillard , R

    Morgan , E. H., Remillard , R. A., & Greiner , J. 1997, , 482, 993

    work page 1997

  27. [27]

    2024, , 691, A314

    Petropoulou , M., Ponti , G., Stel , G., & Mastichiadis , A. 2024, , 691, A314

    work page 2024

  28. [28]

    & Coppi , P

    Poutanen , J. & Coppi , P. S. 1998, Physica Scripta Volume T, 77, 57

    work page 1998

  29. [29]

    1999, , 520, 262

    Psaltis , D., Belloni , T., & van der Klis , M. 1999, , 520, 262

    work page 1999

  30. [30]

    L., Lu , F

    Qu , J. L., Lu , F. J., Lu , Y., et al. 2010, , 710, 836

    work page 2010

  31. [31]

    J., McClintock , J

    Reid , M. J., McClintock , J. E., Steiner , J. F., et al. 2014, , 796, 2

    work page 2014

  32. [32]

    2000, , 541, 883

    Reig , P., Belloni , T., van der Klis , M., et al. 2000, , 541, 883

    work page 2000

  33. [33]

    & Kylafis , N

    Reig , P. & Kylafis , N. D. 2021, , 646, A112

    work page 2021

  34. [34]

    Remillard , R. A. & McClintock , J. E. 2006, , 44, 49

    work page 2006

  35. [35]

    & Petrosian , V

    Stawarz , . & Petrosian , V. 2008, , 681, 1725

    work page 2008

  36. [36]

    & Vietri , M

    Stella , L. & Vietri , M. 1998, , 492, L59

    work page 1998

  37. [37]

    Strogatz, S. H. 1994, Nonlinear Dynamics and Chaos (Reading, MA: Addison Wesley)

    work page 1994

  38. [38]

    & Pellat , R

    Tagger , M. & Pellat , R. 1999, , 349, 1003

    work page 1999

  39. [39]

    & Fiorito , R

    Titarchuk , L. & Fiorito , R. 2004, , 612, 988

    work page 2004

  40. [40]

    P., Churazov , E

    Trudolyubov , S. P., Churazov , E. M., & Gilfanov , M. R. 1999, Astronomy Letters, 25, 718

    work page 1999

  41. [41]

    1989, , 27, 517

    van der Klis , M. 1989, , 27, 517

    work page 1989

  42. [42]

    E., et al

    Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261

    work page 2020

  43. [43]

    2020, , 494, 1375

    Zhang , L., M \'e ndez , M., Altamirano , D., et al. 2020, , 494, 1375

    work page 2020