{"id":"a323cd4f-4583-4837-ac5c-4b46ac2d5003","arxiv_id":"2412.12817","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Oscillating standing shocks in low-angular-momentum relativistic accretion flows produce 0.1-10 Hz luminosity oscillations for stellar-mass black holes and 10^-6 to 10^-5 Hz oscillations for Sgr A*, matching observed X-ray variability.","lead":"Simulations of low-angular-momentum gas spiraling into a black hole show that oscillating shocks make the emitted light flicker at frequencies resembling observed X-ray variability in black hole X-ray binaries and Sgr A*. The result suggests a shared hydrodynamic origin for quasi-periodic oscillations across stellar-mass and supermassive black holes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The observational claim rests on Eq. (13) free-free emission from an unmagnetized flow, while BHXRB and Sgr A* X-rays are dominated by nonthermal or Comptonized processes; a radiative post-processing check is needed before the frequency match is accepted.","rationale":"I read the paper as a numerical study of low-angular-momentum relativistic accretion flows whose observational significance is carried by the claim that shock oscillations produce luminosity oscillations matching LFQPOs in BHXRBs and day-scale variability in Sgr A*. The hydrodynamic machinery is mostly standard: PLUTO with SRHD, a pseudo-Newtonian potential, and a parameter scan in lambda0. The identification of shocks for lambda0 >= 1.75 is interesting and the authors are appropriately cautious about the contrast with the 1D transonic threshold near 1.854 and with Ryu et al. (1997). However, the paper does not include a convergence study, the frequency extraction uses flexible multi-component fits on short simulated durations, and the observational mapping assumes that free-free emission alone tracks the observed X-rays. Of these, the emission-mechanism assumption is the most load-bearing, because even a perfectly converged hydrodynamic simulation would not validate the QPO match if the observed X-ray band is dominated by nonthermal or magnetized emission. The reader identified exactly this weakest assumption, and the requested release of simulation inputs and a resolution test are appropriate. I do not see grounds to reject the paper: the numerical phenomenology may stand on its own, and the observational claim is explicitly qualified by the authors. Therefore the conditional verdict is appropriate and I recommend no change.","tokens_in":17856,"tokens_out":8829,"duration_ms":93201,"concrete_test":"Recompute the light curves and power spectra from the saved lambda0 = 1.75 and 1.80 simulation snapshots with a multi-frequency post-processing code that adds thermal synchrotron and Compton scattering for a plausible magnetic-field prescription (for example, plasma beta in the range 1-100, with field strength normalized to the gas pressure), and compare the 2-10 keV X-ray light curve and its power spectrum with the free-free result from Eq. (13). If the dominant X-ray variability shifts frequency, disappears, or is overwhelmed by nonthermal flare emission, the claimed agreement with observed QPOs and Sgr A* day-scale variability does not hold. This test requires releasing the time-dependent density and temperature fields from the simulations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The decisive step connecting the simulation to observations is the mapping in Fig. 7 through Eq. (13): optically thin thermal bremsstrahlung with a constant Gaunt factor, computed from an adiabatic, unmagnetized, axisymmetric, pseudo-Newtonian SRHD flow. For the central claim to hold, the 0.1-10 Hz oscillations for a 10 solar-mass black hole and the 10^-6 to 10^-5 Hz oscillations for Sgr A* must correspond to what Chandra, Swift, and XMM-Newton actually detect. But the X-ray emission from black hole X-ray binaries in the hard/LFQPO state is generally power-law and Comptonized, not free-free, and Sgr A* X-ray emission is typically attributed to synchrotron self-Compton and nonthermal processes in a magnetized accretion flow. The authors themselves acknowledge in the final discussion paragraph that the model omits magnetic fields and full general relativity. Therefore, even if the hydrodynamic shock oscillations are real and correctly timed, the paper has not established that these oscillations appear in the observed X-ray band. The free-free proxy is the single most load-bearing assumption because it is the only bridge between the simulated flow and the observational signatures that motivate the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 2D special-relativistic hydrodynamic simulations of zero-energy, low-angular-momentum accretion flows onto a Schwarzschild black hole, using the PLUTO code with a pseudo-Newtonian potential. For specific angular momentum values in the range 1.70 <= lambda0 <= 1.80, the authors report oscillating standing shocks, and they compute free-free luminosities from the simulation data. They claim that the resulting luminosity oscillations fall in the 0.1-10 Hz band for a 10 solar-mass black hole, matching low-frequency QPOs in sources like GX 339-4, and in the 10^-6 to 10^-5 Hz band for Sgr A*, matching day-scale X-ray variability observed by Chandra, Swift, and XMM-Newton. The paper also discusses multiple shock merger events, advection/convection properties, and compares the simulation results with the earlier work of Ryu et al. (1997).","tokens_in":18068,"tokens_out":3544,"duration_ms":34936,"significance":"If the central result stands, the paper offers a simple, mass-scalable hydrodynamic mechanism for low-frequency quasi-periodic oscillations and for Sgr A*'s day-scale X-ray variability, linking shock oscillations to observable timing features. The simulations are reproducible in principle, using the public PLUTO code and an open-source ShockFinder tool, and the parameter study across lambda0 is systematic. However, the observational significance is currently limited by the reliance on optically thin thermal free-free emission as the only radiation proxy, which is not representative of the X-ray emission processes believed to operate in the target sources. The manuscript would be a useful contribution to the study of shock instabilities in low-angular-momentum accretion flows, but the observational claims require substantial additional justification or a clearly stated narrower scope.","major_comments":[{"comment":"The claimed match to observed X-ray QPOs and Sgr A* variability rests entirely on Eq. (13), which computes optically thin thermal bremsstrahlung with a constant Gaunt factor from an adiabatic, unmagnetized, axisymmetric flow. However, X-ray emission in the hard/LFQPO state of black hole X-ray binaries and in Sgr A* flares is generally attributed to Comptonization, synchrotron, or synchrotron self-Compton processes, not free-free emission. The authors themselves acknowledge in the final paragraph of Section 5 that the model does not account for magnetic fields or full general relativity. Therefore, the paper has not established that the simulated luminosity oscillations are what Chandra, Swift, and XMM-Newton actually detect; this is a load-bearing gap in the observational claim. The authors should either add a radiative post-processing step with a more realistic emission model or explicitly limit the claim to the hydrodynamic oscillation mechanism without asserting direct observational correspondence.","section":"Section 3, Fig. 3 caption and Section 2.4"},{"comment":"The funnel-region cuts for v > c, described as a way to circumvent numerical errors for lambda0 = 1.75 and 1.85, are post-hoc and suggest that parts of the computational domain violate the special-relativistic constraint. No convergence tests or resolution studies are reported; the manuscript uses a single grid setup (512 logarithmic radial zones, N_theta determined by Eq. 11). The shock positions and the oscillation frequencies extracted from these simulations could be affected by numerical artifacts in the regions where v > c. A resolution study and a quantitative statement about the influence of the v > c regions on the shock radii and power spectra are needed before the claimed frequencies can be considered robust.","section":"Section 4.1, Eq. (14)"},{"comment":"The 'confirmation' of the oscillation peaks in Figure 8 uses Eq. (14) with the shock location R_s and compression ratio R taken from the same simulation. Since the predicted frequency and the measured frequency are derived from the same simulated quantities, this is an internal consistency check rather than an independent confirmation. The wording 'this is in agreement with the oscillation peaks found in Figure 8' overstates the evidential value. The authors should clarify that Eq. (14) is being used as a scaling relation, not as a validation, or compare against an independent analytical prediction based on the inflow parameters alone.","section":"Section 3, Fig. 7 caption"},{"comment":"The luminosity normalization is arbitrary and inconsistently presented. The caption states that luminosity levels span 10^30 to 10^33 erg/s with an inflow density of 10^12 m_p, but then adds a range of 10^33 to 10^37 erg/s for typical gas densities around BHXRBs of 10^10 to 10^11 g/cm^3. Since the simulation density is dimensionless and the physical luminosity depends entirely on the assumed density scale and Gaunt factor, these numbers are not predictions. The authors should state the assumed conversion explicitly and provide a sensitivity analysis for both the luminosity and, if relevant, the oscillation frequencies.","section":"Section 4.1, discussion of shocks below the Rankine-Hugoniot threshold"}],"minor_comments":[{"comment":"There are minor language issues, including 'we shows' in the abstract and some awkward constructions in Section 5; these should be corrected during copyediting.","section":"Abstract"},{"comment":"The notation for the Bernoulli constant is unclear: the text says 'specific energy of the flow or Bernoulli constant, epsilon = ...' but Eq. (4) is written with the right-hand side in code units and does not explicitly define the sign convention. Please clarify the units and the meaning of setting epsilon = 0.","section":"Section 2.2"},{"comment":"The column heading 'T_min' is used but the text refers to 'minimal temperature' and 'initial temperature'; also the units of P_infinity are written as rho_u c^2, which is fine, but the definition of rho_u should be stated. Please reconcile these notations.","section":"Table 1"},{"comment":"The derivation of the velocity distribution at the outer boundary is not fully specified. The sentence 'Opting for the smallest value in this distribution, specifically theta = pi/2' is ambiguous: does this mean the largest equatorial velocity? Please clarify how theta is chosen for the initial conditions.","section":"Section 2.3"},{"comment":"The luminosity ranges quoted in the caption appear to mix two different normalizations without a clear connection. Please revise the caption to state a single, well-defined normalization or explicitly show how the two ranges are related.","section":"Fig. 7"},{"comment":"A few references have formatting issues, such as 'kwan Chan et al. 2009' where the first author's name capitalization is nonstandard, and some arXiv identifiers are absent (e.g., Liska et al. 2021, Kaaz et al. 2022). Please standardize the bibliography.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper's hydrodynamic shock-oscillation mechanism is interesting and plausibly correct as a numerical finding, but the observational mapping via free-free emission is the weakest link and is central to the paper's stated claims. The lack of convergence tests and the post-hoc v>c cuts add to the concern. I would advise the editor that a major revision is required, with a clear request either to implement a more realistic radiative post-processing or to reframe the paper as a purely hydrodynamic simulation study with observational implications discussed only qualitatively. The paper is within the scope of the journal, and the authors' comparison with Ryu et al. (1997) is a useful contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know about this paper: it probably does find new shock behavior in low-angular-momentum accretion flows, but the observational connection it advertises is built on a free-free emission proxy from an unmagnetized flow. The shock physics deserves a serious look; the light-curve mapping does not.\n\nWhat is genuinely new: standing shocks at λ0=1.75 where Ryu et al. 1997 and 1D transonic theory found none; a two-shock merger at λ0=1.80 that produces a luminosity burst; and the mass-scaling that puts oscillation frequencies at 0.1–10 Hz for a 10 Msun BH and 10^-6 to 10^-5 Hz for Sgr A*. The simulations use PLUTO, a public code, and the setup is described in enough detail that a competent group could reproduce the runs. That is real work.\n\nThe soft spots are concentrated in the last third of the paper. No convergence tests or resolution studies are shown. The funnel-region cuts for v>c in Fig. 3 are post-hoc. The luminosity normalization is arbitrary because it depends on a dimensionless injection density. The power spectra are extracted with a flexible multi-component fit (Lorentz + power law + Gaussian), so the peak frequencies should not be over-interpreted. The more serious issue is Eq. (13): thermal bremsstrahlung from an adiabatic, unmagnetized, pseudo-Newtonian flow. The X-ray band in black hole X-ray binaries in the LFQPO state is dominated by Comptonized power-law emission, and Sgr A* X-rays are usually attributed to synchrotron self-Compton from a magnetized flow. The authors themselves acknowledge in the final discussion paragraph that magnetic fields and full GR are omitted. So the claimed matches to GX 339-4 and Sgr A* are not established; they are an interesting scaling exercise. Also, using Eq. (14) with simulated shock locations and compression ratios to \"confirm\" the simulated oscillation frequency is a consistency check, not an independent confirmation.\n\nThe core simulation phenomenology—shock formation, oscillation, merger—is plausible and internally consistent, and the paper engages honestly with the prior literature. The advertised observational alignment is the load-bearing weakness, but it is a weakness in interpretation, not in the hydrodynamics itself.\n\nThis paper is for accretion theorists who care about low-angular-momentum flow dynamics and shock instabilities. It will not convince observers. I would send it to peer review, but with a strong request for a resolution test, release of the simulation setup and outputs, and a major rewrite of the observational section—either add radiative post-processing or explicitly downgrade the QPO/Sgr A* comparison to a speculative scaling.\n\nYes, it deserves a serious referee.","headline":"A competent 2D RHD shock study with genuinely new results at lower angular momentum, but the advertised match to observed QPOs and Sgr A* variability rests on a free-free emission proxy that likely does not represent the actual X-ray bands.","tokens_in":18681,"tokens_out":2620,"would_cite":true,"duration_ms":25395,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Oscillating standing shocks in low-angular-momentum accretion flows can explain both few-Hz X-ray QPOs in stellar-mass black holes and multi-day X-ray variability in Sgr A*.","keywords":["accretion flows","standing shocks","low angular momentum","quasi-periodic oscillations","Sgr A*","black hole X-ray binaries","relativistic hydrodynamics","bremsstrahlung"],"falsifier":"Check the X-ray spectrum of GX 339-4 during a type-C low-frequency QPO: if the few-Hz oscillation is produced by free-free emission, the oscillating component should be thermal and follow the emissivity in Eq. 13; detection of a synchrotron- or Compton-dominated oscillating component, or of QPOs persisting when free-free is negligible, would falsify the mapping. Similarly, if Sgr A* day-scale variability persists in states where the predicted free-free luminosity is negligible, the shock-oscillation explanation fails.","tokens_in":17606,"feed_emoji":"🕳️","tokens_out":8448,"duration_ms":73321,"temperature":0.7,"pith_summary":"This paper tries to show that oscillating standing shocks in zero-energy, low-angular-momentum advective accretion flows can account for two apparently different X-ray variability phenomena with one mechanism. In two-dimensional special-relativistic hydrodynamic simulations, shocks appear for specific angular momentum $\\lambda_0 \\ge 1.75$ at radii $10$–$50\\,R_g$, and for $\\lambda_0$ between $1.70$ and $1.80$ the shocks oscillate strongly enough to modulate the free-free luminosity. Because the flow is dimensionless, the same oscillation frequencies rescale with black hole mass: roughly $0.1$–$10$ Hz for a $10\\,M_\\odot$ black hole, matching low-frequency quasi-periodic oscillations in sources like GX 339-4, and $10^{-6}$–$10^{-5}$ Hz for Sgr A*, matching multi-day X-ray variability seen by long-term monitors. If correct, the paper provides a mass-scalable, shock-driven explanation for low-frequency QPOs in black hole X-ray binaries and for Sgr A*'s day-scale X-ray swings, with shock mergers producing luminosity bursts.","feed_headline":"Oscillating shocks tie X-ray QPOs to Sgr A*'s day cycles","feed_subtitle":"Mass-scaled shock oscillations give 0.1-10 Hz for stellar black holes and multi-day swings for Sgr A*, matching X-ray monitors.","key_machinery":"The load-bearing objects are the standing shock and its oscillatory instability in a zero-energy transonic advective flow. The flow is initialized from a Bernoulli-type energy equation with a pseudo-Newtonian potential, and the shock is located where the radial entropy divergence peaks with the Mach number crossing below unity; its radius $R_s$ is tracked in time. The emission proxy is thermal free-free (bremsstrahlung) luminosity with a constant Gaunt factor, computed from the local density and temperature. The frequency mapping uses the infall-time scaling $\\nu_{\\rm QPO} = \\nu_{s0}\\sqrt{R_s-1}/(\\mathcal{R}\\,R_s)$, where $\\mathcal{R}$ is the compression ratio, which is what rescales the simulated oscillations from stellar-mass to supermassive black holes. The authors suggest the advective-acoustic cycle as the physical driver of the shock oscillations.","core_discovery":"The paper's central claim is that standing shocks are a generic feature of zero-energy low-angular-momentum relativistic accretion flows onto a non-rotating black hole, with shock properties controlled by the conserved specific angular momentum $\\lambda_0$. It finds discernible standing shocks within $10$–$50\\,R_g$ for $\\lambda_0 \\ge 1.75$, whereas the classical transonic one-dimensional analysis would require $\\lambda_0 > 1.854$; the authors attribute the difference to additional sonic points appearing in two dimensions. For $\\lambda_0 \\in [1.70, 1.80]$, the shocks are unstable and oscillate, and the oscillations show up in the thermal free-free light curve computed from the simulated density and temperature fields. At $\\lambda_0 = 1.80$, two shocks merge and produce a dramatic luminosity increase, attributed to convection and the expanding-shock cycle. The recovered oscillation frequencies, $0.1$–$10$ Hz for $10\\,M_\\odot$ and $10^{-6}$–$10^{-5}$ Hz for Sgr A*, are presented as an explanation of type-C low-frequency QPOs in black hole X-ray binaries and of the multi-day X-ray variability of Sgr A*.","pith_inferences":["If this mechanism is right, the same oscillating shock should modulate not only free-free luminosity but also any emission process that traces density and temperature, so multi-wavelength monitoring of Sgr A* could look for correlated day-scale oscillations in infrared and X-ray light curves.","The threshold at $\\lambda_0 \\ge 1.75$ is derived in an unmagnetized flow; adding magnetic fields could shift it, since prior magnetized simulations of similar flows show shocks at lower angular momentum, which would broaden the predicted QPO frequency range.","The paper's use of a constant Gaunt factor and optically thin free-free emission could be replaced with a fuller radiative transfer treatment, and differences between the two light curves would reveal how much of the predicted variability is an artifact of the emission proxy."],"forward_implications":["A $10\\,M_\\odot$ black hole accreting with $\\lambda_0 \\approx 1.75$ should show type-C low-frequency QPOs at roughly $0.4$–$30$ Hz, with the strongest peaks near the few-Hz range reported for sources like GX 339-4.","For Sgr A*, the same shock oscillations translate to periods of one to several days, which is the variability timescale seen in long-term X-ray monitoring.","Shocks with $\\lambda_0 \\ge 1.85$ are stable or outflow-dominated and should not produce detectable low-frequency QPOs, so systems with strong disk winds should be QPO-quiet.","Shock mergers, such as the one seen at $\\lambda_0 = 1.80$, should appear as quasi-periodic luminosity bursts correlated with inward and outward shock motion.","The oscillation frequency is inversely proportional to the infall time, so for fixed $\\lambda_0$ the QPO frequency should scale roughly inversely with black hole mass."],"supporting_citations":[{"why":"Baseline two-dimensional simulations of low-angular-momentum advective flows; the paper's parameter range and injection setup extend this work.","marker":"Ryu et al. (1997)"},{"why":"Transonic flow theory and Rankine-Hugoniot analysis that set the expected shock-forming angular momentum and the supersonic-point framework used to initialize the flow.","marker":"Chakrabarti (1990)"},{"why":"Supplies the pseudo-Newtonian potential used to mimic the Schwarzschild spacetime of a non-rotating black hole.","marker":"Paczyński & Wiita (1980)"},{"why":"Provides the thermal free-free emissivity formula (Eq. 13) that converts simulation density and temperature into luminosity.","marker":"Rybicki & Lightman (1979)"},{"why":"Earlier simulation-based explanation of Sgr A* multi-day luminosity variability against which the present results are compared.","marker":"Okuda et al. (2019)"},{"why":"Reports multiple shocks with inner and outer shocks evolving oppositely; used to interpret the shock merger and luminosity bursts.","marker":"Okuda et al. (2022)"},{"why":"Observational review that anchors the type-C low-frequency QPO frequencies and their association with accretion-rate fluctuations.","marker":"Ingram & Motta (2019)"},{"why":"Shock-oscillation model relating QPO frequency to infall time, providing the scaling used in Eq. (14).","marker":"Molteni et al. (1996)"},{"why":"The shock-capturing relativistic hydrodynamics code used for all simulations.","marker":"Mignone et al. (2007)"}],"fun_headline_variants":["Oscillating shocks explain QPOs and Sgr A* day cycles","Accretion shocks tie X-ray QPOs to Sgr A* variability","Relativistic shocks produce observed black hole flickers","Standing shocks oscillate: link to QPOs and Sgr A*","Shock oscillations mimic GX 339-4 and Sgr A* X-ray signals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the X-ray light curve is faithfully represented by thermal free-free emission from an unmagnetized, adiabatic, axisymmetric flow with a constant Gaunt factor, so if the observed variability is dominated by nonthermal or magnetized emission, the simulated oscillations would not correspond to what the telescopes see.","fun_headline_variants_meta":{"raw":{"variants":["Oscillating shocks explain QPOs and Sgr A* day cycles","Accretion shocks tie X-ray QPOs to Sgr A* variability","Relativistic shocks produce observed black hole flickers","Standing shocks oscillate: link to QPOs and Sgr A*","Shock oscillations mimic GX 339-4 and Sgr A* X-ray signals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00058,"raw_usage":{"total_tokens":2815,"prompt_tokens":1108,"completion_tokens":1707,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":724,"completion_tokens_details":{"reasoning_tokens":1619}},"tokens_in":724,"tokens_out":1707,"duration_ms":12394,"temperature":1.0,"reasoning_tokens":1619,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:42:12.886747+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Check the X-ray spectrum of GX 339-4 during a type-C low-frequency QPO: if the few-Hz oscillation is produced by free-free emission, the oscillating component should be thermal and follow the emissivity in Eq. 13; detection of a synchrotron- or Compton-dominated oscillating component, or of QPOs persisting when free-free is negligible, would falsify the mapping. Similarly, if Sgr A* day-scale variability persists in states where the predicted free-free luminosity is negligible, the shock-oscillation explanation fails.","supporting_citations":[{"cited_title":"K., & Molteni, D","cited_arxiv_id":null,"evidence_quote":"Baseline two-dimensional simulations of low-angular-momentum advective flows; the paper's parameter range and injection setup extend this work."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Transonic flow theory and Rankine-Hugoniot analysis that set the expected shock-forming angular momentum and the supersonic-point framework used to initialize the flow."},{"cited_title":", & Lightman, A","cited_arxiv_id":null,"evidence_quote":"Provides the thermal free-free emissivity formula (Eq. 13) that converts simulation density and temperature into luminosity."},{"cited_title":"B., Das, S., et al","cited_arxiv_id":null,"evidence_quote":"Earlier simulation-based explanation of Sgr A* multi-day luminosity variability against which the present results are compared."},{"cited_title":"B., & Aktar, R","cited_arxiv_id":null,"evidence_quote":"Reports multiple shocks with inner and outer shocks evolving oppositely; used to interpret the shock merger and luminosity bursts."},{"cited_title":"R., & Motta, S","cited_arxiv_id":null,"evidence_quote":"Observational review that anchors the type-C low-frequency QPO frequencies and their association with accretion-rate fluctuations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shock-oscillation model relating QPO frequency to infall time, providing the scaling used in Eq. (14)."},{"cited_title":"Our Non-Stable Universe","cited_arxiv_id":null,"evidence_quote":"The shock-capturing relativistic hydrodynamics code used for all simulations."}],"review_version":1}