pith. sign in

arxiv: 2604.13199 · v1 · submitted 2026-04-14 · 🌌 astro-ph.HE · astro-ph.GA

Radiation hydrodynamic simulations for the origin of quasi-periodic oscillations for accretion onto supermassive black holes

Yiyang Lin , Erlin Qiao , Jifeng Liu , Meng Guo , Zikun Lin This is my paper

Pith reviewed 2026-05-10 13:56 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA
keywords quasi-periodic oscillationssupermassive black holesaccretion flowsradiation hydrodynamicsradial epicyclic frequencyactive galactic nucleitidal disruption eventsblack hole mass scaling
0
0 comments X

The pith

Radiation hydrodynamic simulations show quasi-periodic oscillations arise from the maximum radial epicyclic frequency around supermassive black holes.

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

The paper runs radiation hydrodynamic simulations of accretion onto a 10 million solar mass black hole by injecting gas at a fixed radius of 10 Schwarzschild radii. Time series of the mass inflow rates at different locations inside the flow display clear periodic signals. These periods match the radial epicyclic frequencies computed from analytic orbital mechanics for all radii beyond 3.8 Schwarzschild radii. The match occurs precisely where the epicyclic frequency reaches its highest value, leading the authors to identify that peak frequency as a reliable proxy for the QPOs seen in real active galactic nuclei and tidal disruption events. The same relation, when scaled by black hole mass, reproduces the observed frequency-mass trend across both classes of objects.

Core claim

In radiation hydrodynamic simulations of accretion flows with mass injected at 10 Schwarzschild radii around a 10^7 solar mass black hole, the mass inflow rates exhibit quasi-periodic oscillations whose frequencies agree with the radial epicyclic frequencies obtained from analytic calculations for radii larger than 3.8 Schwarzschild radii. This radius marks the location of the maximum radial epicyclic frequency ν_r,max. The authors therefore propose that ν_r,max serves as a direct proxy for the observed QPO frequency ν_QPO, and that the resulting frequency-mass relation reproduces the trends measured in samples of AGN and TDE systems.

What carries the argument

The radial epicyclic frequency of orbital motion in the accretion flow, specifically its maximum value ν_r,max at 3.8 Schwarzschild radii, which the inflow-rate time series in the simulations reproduce as the QPO signal.

If this is right

  • QPO frequency scales inversely with black hole mass in a manner that already matches the compiled AGN and TDE data.
  • The inner accretion flow at radii greater than 3.8 Schwarzschild radii dominates the observed periodicity.
  • Changes in mass-injection radius, viscosity, or magnetic field strength can be explored in follow-up runs while preserving the core frequency match.
  • General-relativistic corrections to the epicyclic frequency remain small enough that the Newtonian analytic curve still provides a usable proxy.

Where Pith is reading between the lines

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

  • If the identification is correct, QPO detections could be inverted to estimate the radial location of the emitting gas in the disk.
  • The same mechanism could be tested in stellar-mass black-hole binaries by running otherwise identical simulations at lower masses.
  • Full general-relativistic radiation hydrodynamics would provide a direct check on whether the 3.8 Schwarzschild radius location shifts appreciably.

Load-bearing premise

The simulation results obtained at one black hole mass can be scaled directly to other masses to predict observed QPO frequencies without further changes in disk structure or radiative physics.

What would settle it

A measured QPO frequency in a new AGN or TDE whose black hole mass places the point significantly off the predicted ν_r,max versus M_BH curve would show the proposed proxy does not hold.

Figures

Figures reproduced from arXiv: 2604.13199 by Erlin Qiao, Jifeng Liu, Meng Guo, Yiyang Lin, Zikun Lin.

Figure 1
Figure 1. Figure 1: Snapshots of gas density (in 𝑟 − 𝜃 plane) at 𝑡 = 10, 50, 300 day since the injection of matter at 𝑅inject with 𝑀¤ inject = 0.01𝑀¤ Edd [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Snapshots of gas density (in 𝑟 − 𝜃 plane) at 𝑡 = 10, 50, 300 day since the injection of matter at 𝑅inject with 𝑀¤ inject = 0.1𝑀¤ Edd [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Snapshots of gas density (in 𝑟 − 𝜃 plane) at 𝑡 = 10, 20, 40, 300 day since the injection of matter at 𝑅inject with 𝑀¤ inject = 2𝑀¤ Edd. that the matter in the accretion flow is very well constrained in the equatorial plane for 𝑀¤ inject = 2𝑀¤ Edd as for 𝑀¤ inject = 0.01𝑀¤ Edd and 𝑀¤ inject = 0.1𝑀¤ Edd, in which the radial oscillation naturally can occur. However, if the accretion flow is in the high state,… view at source ↗
Figure 4
Figure 4. Figure 4: Evolution of the scaled mass inflow rate 𝑀¤ in/𝑀¤ Edd, the corresponding power spectrum and the dynamical power spectrum with mass injection rate 𝑀¤ inject = 0.01𝑀¤ Edd. Left panel: 𝑀¤ in/𝑀¤ Edd as a function of 𝑡 at 𝑟 = 2𝑅S (top) and 𝑟 = 3.8𝑅S (bottom) respectively. The purple region and the red region refer to phase 1 (0 − 41) day and phase 2 (41 − 400) day respectively. Middle panel: Power spectrum calc… view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of the scaled mass inflow rate 𝑀¤ in/𝑀¤ Edd, the corresponding power spectrum and the dynamical power spectrum with mass injection rate 𝑀¤ inject = 0.1𝑀¤ Edd. Left panel: 𝑀¤ in/𝑀¤ Edd as a function of 𝑡 at 𝑟 = 2𝑅S (top) and 𝑟 = 3.8𝑅S (bottom) respectively. The purple region and the red region refer to phase 1 (0 − 25) day and phase 2 (25 − 400) day respectively. Middle panel: Power spectrum calcu… view at source ↗
Figure 6
Figure 6. Figure 6: Evolution of the scaled mass inflow rate 𝑀¤ in/𝑀¤ Edd, the corresponding power spectrum and the dynamical power spectrum with mass injection rate 𝑀¤ inject = 2𝑀¤ Edd. Left panel: 𝑀¤ in/𝑀¤ Edd as a function of 𝑡 at 𝑟 = 2𝑅S (top) and 𝑟 = 3.8𝑅S (bottom) respectively. The purple region and the red region refer to phase 1 (20 − 30) day and phase 2 (33 − 47) day respectively. Middle panel: Power spectrum calcula… view at source ↗
Figure 8
Figure 8. Figure 8: Relation between the observed QPO frequency 𝜈QPO and the BH mass 𝑀BH for five sources, including two AGNs, i.e, RE J1034+396, 1ES 1927+654, as well as three TDEs, i.e., Swift J164449.3+573451, ASASSN￾14li and 2XMM J123103.2+110648 (TDE candidate). The black solid line is the maximum radial epicyclic frequency 𝜈r,max as a function of 𝑀BH, i.e equation (5) in the paper. and 2 times lower than that of 𝑎 = 0.9… view at source ↗
Figure 7
Figure 7. Figure 7: The theoretical Keplerian orbital frequency, the epicyclic fre￾quency of the accretion disk, as well as the QPO frequency derived from our simulations as a function of radius 𝑟/𝑅S for 𝑀¤ inject = 0.01𝑀¤ Edd (top), 𝑀¤ inject = 0.1𝑀¤ Edd (middle) and 𝑀¤ inject = 2𝑀¤ Edd (bottom) respectively. The red solid lines and the blue solid line represent the theoretical Keplerian or￾bital frequency (equation (3)) and… view at source ↗
read the original abstract

Quasi-periodic oscillation (QPO) has been detected in several accreting supermassive black hole (SMBH) systems, including active galactic nuclei (AGNs) and tidal disruption events (TDEs). However, despite that several models have been proposed, the physical origin of QPO is still unclear. In this paper, we performed radiation hydrodynamic simulations of accretion flow by injecting mass at a fixed radius, i.e. 10 Schwarzschild radius with different mass accretion rates, and setting the black hole (BH) mass to $10^7M_{\odot}$. We find that there are QPO signals by analyzing the mass inflow rates as a function of time from the simulations for different radii. The QPO frequencies from our simulations are well consistent with the radial epicyclic frequencies from analytic calculations for radius greater than a critical radius 3.8 Schwarzschild radius. This critical radius corresponds to the maximum epicyclic frequency, i.e. $\nu_{\rm r,max}$, in the radial direction. We proposed that $\nu_{\rm r,max}$ can be a good proxy for the observed QPO $\nu_{\rm QPO}$. Furthermore, assuming that our simulation results can be scaled to different BH masses $M_{\rm BH}$, we find that the theoretical relation of $\nu_{\rm r,max}$ as a function of $M_{\rm BH}$ can well match $\nu_{\rm QPO}$ as a function of $M_{\rm BH}$ for a sample of AGN and TDE. Finally, we discuss the effects of the BH mass, general relativity (GR), and other possible factors including the size of the mass injecting radius, viscosity and magnetic field on the simulation results.

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 / 2 minor

Summary. The paper performs radiation-hydrodynamic simulations of accretion onto a 10^7 M_⊙ black hole with mass injection at a fixed radius of 10 R_s and varying accretion rates. It reports QPO signals in the mass-inflow-rate time series whose frequencies match analytic radial epicyclic frequencies for radii > 3.8 R_s, proposes ν_r,max as a proxy for observed ν_QPO, and shows that scaling the simulation results to other M_BH produces a relation that matches a sample of AGN and TDE QPOs. Effects of BH mass, GR, injection radius, viscosity, and magnetic fields are discussed.

Significance. If the reported frequency match is robust and the 1/M scaling preserves the QPO–epicyclic link under radiation-hydrodynamics, the work would provide a concrete physical origin for QPOs in SMBH accretion flows and a predictive relation across mass scales. The radiation-hydrodynamic treatment is a strength relative to purely hydrodynamic or test-particle models, and the direct comparison to analytic epicyclic frequencies is a clear, falsifiable test.

major comments (3)
  1. [§4] §4 (scaling to different M_BH and comparison to observations): The central claim that ν_r,max serves as a proxy for observed ν_QPO across AGN/TDE samples rests on rescaling the fixed-M_BH=10^7 M_⊙ radiation-hydro runs by simple 1/M frequency scaling. The radiation-hydrodynamic equations contain dimensionful terms (radiative cooling, opacity, diffusion) whose dimensionless ratios change with M_BH unless density, temperature, and optical depth are explicitly rescaled to hold all timescales fixed; no such rescaling test or demonstration that the inflow-rate QPO remains locked to the analytic ν_r is provided.
  2. [§3] §3 (QPO extraction and frequency comparison): The abstract states that QPO frequencies are extracted from mass-inflow-rate time series and found consistent with radial epicyclic frequencies beyond 3.8 R_s, yet no information is given on the periodogram method, windowing, statistical significance thresholds, error bars on the reported frequencies, or robustness against variations in viscosity parameter or grid resolution.
  3. [§3.1] §3.1 (critical radius and ν_r,max): The identification of 3.8 R_s as the radius of maximum epicyclic frequency and the assertion that simulated QPO frequencies track ν_r only outside this radius are load-bearing for the proposed proxy; the manuscript does not show the radial profile of simulated QPO power or demonstrate that the match degrades inside 3.8 R_s in a manner consistent with the analytic ν_r curve.
minor comments (2)
  1. [Abstract and §5] The abstract and §5 mention that GR effects are discussed, but it is unclear whether the simulations themselves are performed in Newtonian gravity or include any GR corrections to the metric or radiative transfer.
  2. [Figures] Figure captions and axis labels for the time-series and frequency spectra should explicitly state the time baseline, sampling cadence, and units of the power spectra to allow independent verification of the reported periods.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough and constructive review. We address each major comment point by point below, indicating where revisions will be made to strengthen the manuscript.

read point-by-point responses

  1. Referee: [§4] §4 (scaling to different M_BH and comparison to observations): The central claim that ν_r,max serves as a proxy for observed ν_QPO across AGN/TDE samples rests on rescaling the fixed-M_BH=10^7 M_⊙ radiation-hydro runs by simple 1/M frequency scaling. The radiation-hydrodynamic equations contain dimensionful terms (radiative cooling, opacity, diffusion) whose dimensionless ratios change with M_BH unless density, temperature, and optical depth are explicitly rescaled to hold all timescales fixed; no such rescaling test or demonstration that the inflow-rate QPO remains locked to the analytic ν_r is provided.

    Authors: We thank the referee for highlighting this important caveat in the scaling argument. Our simulations are performed at a fixed black hole mass of 10^7 M_⊙, and the QPO signals arise primarily from the dynamical epicyclic motions in the accretion flow, which are governed by the gravitational potential. While radiation-hydrodynamic effects are self-consistently included, the frequency matching is to the analytic radial epicyclic frequency, which scales purely as 1/M_BH. We agree that a full verification of the scaling under radiation hydrodynamics would require additional simulations with appropriately rescaled parameters for different masses, which is computationally intensive. In the revised manuscript, we will expand the discussion in §4 to explicitly state the assumptions underlying the 1/M scaling and acknowledge that the radiation terms may introduce mass-dependent effects not captured here. We will also suggest this as a direction for future work. revision: partial

  2. Referee: [§3] §3 (QPO extraction and frequency comparison): The abstract states that QPO frequencies are extracted from mass-inflow-rate time series and found consistent with radial epicyclic frequencies beyond 3.8 R_s, yet no information is given on the periodogram method, windowing, statistical significance thresholds, error bars on the reported frequencies, or robustness against variations in viscosity parameter or grid resolution.

    Authors: We acknowledge that the manuscript lacks sufficient detail on the QPO analysis procedure. In the revised version, we will add a dedicated subsection in §3 describing the time series analysis, including the periodogram method employed, any window functions applied, the criteria for identifying significant peaks (e.g., false alarm probability thresholds), and how uncertainties on the frequencies were estimated. Additionally, we will include a brief discussion or appendix on the robustness of the results with respect to changes in the viscosity parameter α and grid resolution, based on the simulations we have performed. revision: yes

  3. Referee: [§3.1] §3.1 (critical radius and ν_r,max): The identification of 3.8 R_s as the radius of maximum epicyclic frequency and the assertion that simulated QPO frequencies track ν_r only outside this radius are load-bearing for the proposed proxy; the manuscript does not show the radial profile of simulated QPO power or demonstrate that the match degrades inside 3.8 R_s in a manner consistent with the analytic ν_r curve.

    Authors: The critical radius is determined from the analytic calculation of the radial epicyclic frequency ν_r(r), which peaks at approximately 3.8 R_s for the Schwarzschild metric. In our simulations, we find that the extracted QPO frequencies from the mass inflow rate time series at different radii align with the analytic ν_r for r > 3.8 R_s. To strengthen this, in the revised manuscript we will include a new figure in §3.1 showing the power spectrum or QPO frequency as a function of radius from the simulations, overlaid with the analytic ν_r profile, to illustrate where the match holds and where it deviates inside the critical radius. revision: yes

Circularity Check

0 steps flagged

No significant circularity; simulations and comparisons are independent

full rationale

The paper's central chain consists of independent radiation-hydrodynamic simulations (mass injection at fixed 10 Rs, M_BH=10^7 M_⊙, varying accretion rates) that generate time series of mass inflow rates, from which QPO frequencies are extracted via analysis. These frequencies are then compared to standard analytic radial epicyclic frequencies computed from general-relativistic orbital mechanics, which are external formulas not derived or fitted within the paper. The identification of a critical radius at 3.8 Rs (corresponding to ν_r,max) and the proposal that it proxies observed ν_QPO follow directly from this comparison. Scaling the results to other M_BH values is stated explicitly as an assumption to enable observational matching, but introduces no self-definitional loop, fitted-parameter renaming, or load-bearing self-citation. No step reduces a claimed prediction to an input by construction.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard radiation-hydrodynamic assumptions, Newtonian orbital mechanics for epicyclic frequencies, and the untested premise that simulated inflow-rate variability directly corresponds to observed QPO signals.

free parameters (2)
  • mass injection radius = 10 r_s
    Fixed at 10 Schwarzschild radii for all runs; choice directly sets the inner boundary condition.
  • black hole mass in base runs = 10^7 M_sun
    Set to 10^7 solar masses before any scaling is applied.
axioms (2)
  • domain assumption Radiation hydrodynamics without full general relativity or magnetic fields suffices to capture the relevant variability
    Invoked by the choice of simulation method and the later discussion of GR effects as secondary.
  • domain assumption Quasi-periodic signals extracted from mass inflow rate time series at fixed radii correspond to the physical QPO mechanism
    Central step linking simulation output to the proposed ν_r,max proxy.

pith-pipeline@v0.9.0 · 5631 in / 1452 out tokens · 35377 ms · 2026-05-10T13:56:12.974568+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

1 extracted references · 1 canonical work pages

  1. [1]

    A., Fragile P

    Abramowicz M. A., Fragile P. C., 2013, Living Reviews in Relativity, 16, 1 AbramowiczM.A.,CzernyB.,LasotaJ.P.,SzuszkiewiczE.,1988,ApJ,332, 646 Alston W. N., Markeviciute J., Kara E., Fabian A. C., Middleton M., 2014, MNRAS, 445, L16 Blandford R. D., Payne D. G., 1982, MNRAS, 199, 883 Bloom J. S., et al., 2011, Science, 333, 203 Brown J. S., Holoien T. W.-...

    work page doi:10.1007/978-981-16-4544-0_127-1 2013