Recognition: 4 theorem links
· Lean TheoremPeriodic Orbits and Gravitational Wave Radiation of Black Hole in EGB gravity
Pith reviewed 2026-05-06 18:10 UTC · model claude-opus-4-7
The pith
A Gauss-Bonnet correction plus charge leaves shadow- and waveform-level fingerprints that current data already partly bound.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Within a charged 4D Einstein-Gauss-Bonnet black hole spacetime, increasing either the Gauss-Bonnet coupling α or the charge Q monotonically shrinks the radius, energy, and angular momentum of both the marginally bound orbit and the ISCO, and shifts the bound-orbit region in the (E,L) plane. Current shadow measurements meaningfully bound α and Q jointly, but stellar precession bounds only Q (since α enters first at order M³/a³). For specific zoom-whirl orbits classified by q = w + v/z, the resulting quadrupole gravitational waveforms differ in amplitude and phase from the Schwarzschild and Reissner-Nordström cases by amounts the authors argue are large enough to serve as candidate discriminat
What carries the argument
Two pieces do the work. First, Fernandes' static spherically symmetric charged black hole solution of the regularized 4D-EGB action, whose metric function f(r) = 1 + (r²/2α)[1 − √(1 + 4α(2M/r³ − Q²/r⁴))] supplies a closed-form geodesic problem in which α and Q enter the effective potential directly. Second, Levin and Perez-Giz's classification of bound orbits by a rational q = Δφ/2π − 1 = w + v/z (zoom, whirl, vertex), which lets the authors index a discrete family of periodic orbits and feed each into a quadrupole (Kludge) waveform calculation, turning a metric-level correction into a waveform-level signature.
If this is right
- <parameter name="0">Shadow measurements of M87* and Sgr A* already cut into the (α
- Q) plane allowed by horizon existence alone
- with Sgr A* providing the tighter constraint.
Where Pith is reading between the lines
- <parameter name="0">The fact that α first enters S2-style precession at order M³/a³ while Q enters at order M/a is a generic structural statement about how higher-curvature couplings hide from weak-field tests
- the same hierarchy likely governs other quadratic-curvature theories and explains why solar-system-like bounds on such couplings tend to be weak.
Load-bearing premise
The argument treats the inspiraling object as a neutral test particle moving on geodesics of a fixed background, ignoring back-reaction and using the simplest quadrupole waveform — so any signature claimed in the waveform has to survive a real, dissipative, long-duration inspiral calculation before it can be called observable.
What would settle it
A full adiabatic EMRI evolution in this charged 4D-EGB background, with a sensitivity-weighted match against a Reissner-Nordström template over a realistic LISA observation, that returns a mismatch below the detection threshold for all (α, Q) consistent with EHT and S2 bounds — that would show the per-orbit waveform differences shown here do not survive into a distinguishable signal. Conversely, a future shadow or pulsar-timing measurement that pins α at a value outside the region this paper maps as horizon-allowed would falsify the underlying solution as a viable astrophysical model.
read the original abstract
This paper investigates the orbital dynamics and gravitational wave radiation characteristics of neutral test particles around a static spherically symmetric charged black hole (BH) in 4D Einstein-Gauss-Bonnet (4D-EGB) gravity theory. We analyze the dependence of the marginally bound orbit (MBO) and the innermost stable circular orbit (ISCO) on the Gauss-Bonnet coupling parameter $\alpha$ and charge $Q$. The results indicate that the orbital radius, angular momentum, and energy all decrease with increasing $\alpha$ or $Q$, with the corresponding bound orbit region shifting leftward in the $(E, L)$ parameter space. By combining observational data from the BH shadows of M87* and Sgr A* as well as the orbital precession of the S2 star, we constrain the model parameters and find that existing observations can limit the ranges of $\alpha$ and $Q$ to a certain extent. Furthermore, we investigate the characteristics of periodic orbits corresponding to different rational numbers $q$ and the gravitational waveforms they excite, finding that variations in $\alpha$ and $Q$ can lead to distinguishable differences in periodic orbit structures and gravitational wave phases. This study contributes to understanding the effects of Gauss-Bonnet corrections on BH spacetimes, and the results may provide theoretical references for future gravitational wave observations of extreme mass ratio inspirals (EMRIs).
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies timelike geodesics of neutral test particles around the static spherically symmetric charged 4D-EGB black hole derived by Fernandes [39]. The authors (i) characterize the dependence of the marginally bound orbit and ISCO on the Gauss–Bonnet coupling α and charge Q; (ii) constrain (α, Q) using EHT shadow sizes of M87* and Sgr A* and the GRAVITY measurement of S2 perihelion precession (deriving an explicit weak-field expansion, Eq. (24)); (iii) classify periodic orbits via the Levin–Perez-Giz (z,w,v) scheme for several representative (α, Q); and (iv) compute quadrupole "kludge" gravitational waveforms for the (3,1,2) zoom-whirl orbit and report visible differences in amplitude and phase across the parameter values considered. The headline observational claim is that variations in α and Q produce distinguishable imprints in periodic-orbit waveforms relevant to EMRI observations.
Significance. The technical content is a competent if incremental extension of a now-standard methodology (geodesics → MBO/ISCO → periodic-orbit classification → quadrupole waveforms) to the Fernandes 4D-EGB charged BH. The two most useful pieces are the explicit weak-field PPN-style expansion of the precession ratio in Eq. (24), which transparently explains why S2 constrains Q at leading order but is essentially blind to α (entering only at O(M³/a³)), and the joint use of EHT M87*/Sgr A* shadow bounds and S2 precession to delimit (α, Q). The authors are commendably explicit about the limitations: they normalize everything by M, restrict to the negative branch with motivated reasoning, acknowledge the D→4 controversy and its resolution via regularized scalar–tensor theories, and flag in the final paragraph that genuine detectability requires Fisher/Bayesian degeneracy analysis with M jointly estimated. Significance is moderate: results are consistent with the literature on RN and other modified-gravity ISCO/periodic-orbit analyses (e.g., refs. [48, 51] by overlapping authors), and the gravitational-wave "distinguishability" claim is preliminary rather than detection-grade.
major comments (4)
- [§V, Figs. 14–15 and Summary] The headline distinguishability claim — that α and Q produce 'distinguishable differences in periodic orbit structures and gravitational wave phases' — is not actually demonstrated by the analysis presented. Figs. 14–15 compare h_+(t) at fixed q=(3,1,2) and either fixed E or fixed L, while in any real EMRI inference (M, E, L, α, Q) are estimated jointly. With M scaled out (r, α/M², Q/M dimensionless) and no degeneracy/Fisher analysis, an O(1) shift in (M, E, L) within a Schwarzschild template could plausibly absorb much of the displayed waveform difference. The authors themselves note this caveat in the final paragraph. To support the claim as stated in the abstract, the paper should at minimum show one quantitative diagnostic that is not absorbed by reparametrization — e.g., a mismatch/overlap between an EGB waveform and the best-fit Schwarzschild template after maximizing over (M, E, L
- [§V, Eqs. (27)–(29)] The kludge waveform expressions (28)–(29) are written as if they hold for a generic bound orbit, but they are derived under quasi-circular assumptions (single phase ϕ, fixed inclination ι, simple latitude ξ), whereas zoom-whirl orbits with multiple turning points have a non-trivial radial decomposition. The procedure for evaluating ϕ(t), r(t) and how the radial motion is folded into h_+ should be stated explicitly (is r(t) inserted as 1/r in front, with ϕ taken from the geodesic integration, with no separate radial mode? is the trace removal applied?). Without this the figures cannot be reproduced and the relation between the periodic-orbit geometry and the resulting waveform amplitude in Figs. 14–15 is not auditable.
- [§III B, Eq. (24)] Two checks are needed on the weak-field precession expansion. (a) The leading O(M/a) coefficient (6 − Q²) is the standard RN result; please confirm by quoting the source or showing the derivation that the coupling α first enters at O(M³/a³) and that no O(M²/a²) cross-term in α appears (this is the load-bearing input behind the conclusion that S2 cannot constrain α). (b) For S2 around Sgr A*, semi-major axis a ≈ 1000 AU and M ≈ 4×10⁶ M_⊙ give M/a ~ 10⁻⁴, but α/M² is allowed by Fig. 9 to be O(1) so the O(M³/a³) term scaled by α/M² is ~10⁻¹² relative to the leading 1PN term — far below the 17% S2 measurement uncertainty. The reported constraint Q ∈ (0, 0.738) and 'no significant constraint on α' should therefore be made quantitative (what range of α is actually allowed at 1σ?), and the assertion that S2 constrains Q at all should be checked at the same level of rigor.
- [§II A and Fig. 1] The boundaries of regions A (double horizon) and B (single horizon) are stated without derivation, and the horizon equation Eq. (5) is r_h = M ± √(M² − Q² − α), which only manifestly admits two horizons when M² − Q² − α > 0. The roles of the inner discriminant 1 + 4α(2M/r³ − Q²/r⁴) ≥ 0 inside the square root of f(r), and the consequent existence of curvature singularities or branch-cut surfaces outside r_h for α < 0, should be addressed. In particular, for the negative-α periodic orbits in Figs. 12–15, please confirm that the orbit lies entirely in the region where the metric is real and where no additional inner horizon or singularity affects the geodesic.
minor comments (9)
- [Abstract & §I] 'distinguishable differences in periodic orbit structures and gravitational wave phases' should be qualified to reflect that no waveform overlap or Fisher analysis is presented; suggest 'morphological differences in the test-particle waveform that motivate a future detectability study'.
- [§II A footnote/text] The statement 'in the limit of α→0 and r→∞, only the negative branch can reduce to the Reissner–Nordström black hole' is correct but the qualifier r→∞ is unnecessary for the metric form; please rephrase, since RN is recovered in the entire spacetime, not just asymptotically.
- [§III A, Fig. 7 caption] Please state explicitly what the lower-bound region in Fig. 7 corresponds to (1σ? 2σ?) and quote the resulting allowed interval for (α, Q) numerically, mirroring what is done in §III B.
- [Tables I–II] Energies are quoted to six decimal places and angular momenta to five; the differences across rows are at the fourth–fifth decimal, so rounding errors in subsequent geodesic integration could matter. Please specify the integration tolerance and whether (E, L) were re-tuned to land precisely on the (z,w,v) orbit.
- [§IV, Eq. (26)] Eq. (26) requires r_1, r_2 (the turning points) but they are introduced implicitly. State that they are roots of E² = V_eff(r) at given (E, L) and clarify which root is the apsidal pair used.
- [Figs. 12–15] Curve identification is hard to read at print size: legend entries overlap and dashed/solid distinctions between Schwarzschild, RN, and 4D-EGB are not always obvious. Larger legends and a small inset zooming on the periastron passage would help.
- [References] Refs. [3], [4], [5] (Davis–Peebles, Tonry–Kochanek, Campbell et al.) are cited as observational confirmations of GR but are large-scale-structure / Hubble-flow / unrelated entries; please replace with appropriate GR-test references (e.g., Cassini Shapiro delay, binary pulsar tests).
- [Throughout] Numerous minor English issues ('have promoted researchers to explore' → 'motivated researchers to explore'; 'In Section II introduces' → 'Section II introduces'). A careful language pass is recommended.
- [§V parameters] M = 10⁶ M_⊙, m = 10 M_⊙, D_L = 200 Mpc are stated, but ι = ξ = π/4 is asserted without comment — note that ξ is the orbital phase offset and is degenerate with the time origin; clarify whether this choice affects any reported feature.
Simulated Author's Rebuttal
We thank the referee for a careful and constructive report. The four major points are well taken; we agree with the substance of all of them and propose the following revisions. (1) We will soften the headline distinguishability claim in the abstract and §V to make clear that the visible waveform differences in Figs. 14–15 are illustrative of the orbital dynamics, not a detection-grade diagnostic, since (M,E,L,α,Q) are not jointly varied; a full mismatch/Fisher analysis is flagged as the natural follow-up. (2) We will rewrite the opening of §V to state explicitly the kludge implementation: r(t) and ϕ(t) come from numerical geodesic integration on the equatorial plane, are inserted into the quadrupole formula Eq. (27) and projected to the TT gauge to give Eqs. (28)–(29); the trace-removal convention follows Babak et al. [68]. (3) We will add intermediate steps for Eq. (24) showing that α enters f(r) only through the inner square root and therefore first contributes at O(M³/a³) with no O(M²/a²) cross-term; we will further quantify the S2 constraint, replacing the qualitative 'no constraint on α' with an explicit (very loose) 1σ interval and providing a properly marginalized Q contour. (4) We will derive the A/B boundaries of Fig. 1 explicitly and verify, for every α<0 orbit shown, that the trajectory lies outside the surface where the inner discriminant of f(r) vanishes. The headline conclusions of the paper — the qualitative behavior of MBO/ISCO, the EHT/S2 parameter constrain
read point-by-point responses
-
Referee: The headline distinguishability claim in §V/Figs. 14–15 is not actually demonstrated: with (M,E,L,α,Q) jointly estimated, the visible waveform differences could be absorbed by reparametrization. A mismatch/overlap diagnostic maximized over (M,E,L) is needed to support the claim in the abstract.
Authors: We agree with the referee. As acknowledged in the final paragraph of §VI, the present work is restricted to single-period waveform features in the test-particle, fixed-background regime, and a genuine detectability statement requires a Fisher/Bayesian degeneracy analysis with M jointly estimated. We will (i) soften the abstract and §V wording: replace 'distinguishable differences in periodic orbit structures and gravitational wave phases' with 'parameter-dependent imprints in single-period waveforms whose observational distinguishability remains to be assessed once parameter degeneracies with (M,E,L) are accounted for'; (ii) add an explicit caveat at the end of §V stating that the visible amplitude/phase shifts in Figs. 14–15 are computed at fixed (M,E,L,q) and are not invariant under reparametrization, so they should be interpreted as illustrative of the orbital dynamics rather than as a detection-grade diagnostic. A full mismatch/overlap calculation maximized over (M,E,L), and ultimately a Fisher matrix study against LISA/TianQin sensitivities, is the natural follow-up and is flagged as future work in the revised §VI. We prefer not to add a partial mismatch calculation in this revision because a credible overlap requires consistent inspiral evolution (adiabatic dE/dt, dL/dt) that is beyond the present scope. revision: partial
-
Referee: The kludge waveform expressions (28)–(29) are written as if generic, but were derived under quasi-circular assumptions. For zoom-whirl orbits with multiple turning points the procedure for evaluating ϕ(t), r(t), how radial motion enters h_+, and whether the trace is removed should be stated explicitly so the figures are reproducible.
Authors: We thank the referee for flagging this presentation gap. In our implementation we numerically integrate the geodesic equations in (t,r,ϕ) on the equatorial plane using the conserved E,L from Tables I–II, obtaining r(t) and ϕ(t) along the full zoom-whirl trajectory (i.e., r(t) is not held fixed; it oscillates between the turning points r_1, r_2 of Fig. 5). We then evaluate Eq. (27) with v^i = dx^i/dt and n^i = x^i/r in Cartesian equatorial coordinates, project onto the TT gauge for an observer at inclination ι (with the trace removal performed in the standard way), and take the +,× components to obtain Eqs. (28)–(29) — i.e., these formulas are used as the projected TT-gauge result with ϕ → ϕ(t) and r → r(t) inserted, not under a quasi-circular ansatz. We will (a) add a paragraph at the start of §V making this procedure explicit, (b) write the unprojected Q^{ij} form of Eq. (27) before reducing to (28)–(29), clarifying that the reduction assumes equatorial motion and a single observer-frame phase but retains the full r(t) and ϕ(t) from the geodesic integration, and (c) cite the kludge prescription of Babak et al. [68] for the trace-removal convention. We will also add a footnote noting that for strongly whirling orbits the quadrupole approximation is itself a known limitation, providing further reason to treat Figs. 14–15 as illustrative. revision: yes
-
Referee: Two checks on Eq. (24): (a) confirm the RN leading coefficient (6−Q²) and that α first enters at O(M³/a³) with no O(M²/a²) cross-term; (b) given M/a ~ 10⁻⁴ for S2, the α-dependent term is ~10⁻¹² relative to 1PN, far below the 17% uncertainty, so the claimed Q∈(0,0.738) constraint and the 'no constraint on α' statement should be quantified (what α range is allowed at 1σ?).
Authors: (a) The leading (6−Q²)M/[(1−e²)a] coefficient reproduces the standard RN 1PN periastron advance (see, e.g., Heydari-Fard et al., and Eq. (24) reduces to ∆ω_GR = 6πM/[a(1−e²)] when Q=α=0, as stated). We re-derived Eq. (24) by expanding Eq. (22) symbolically; α enters f(r) only through the combination 4α(2M/r³ − Q²/r⁴) inside the square root of Eq. (4), so the lowest-order α-dependent correction to f(r) is 8αM²/r⁴ − 4αMQ²/r⁵+…, which contributes to ∆ϕ first at O(M³/a³). There is indeed no O(M²/a²) α cross-term, and we will add the intermediate steps to an appendix in the revised version. (b) The referee is correct that for S2 (M/a ~ 10⁻⁴) the α-coefficient at O(M³/a³) is parametrically too small to be bounded by f_sp = 1.10 ± 0.19; quantitatively, the α range allowed at 1σ spans essentially the entire BH-existence window of Fig. 1, which is exactly why we wrote 'no significant constraint on α'. We will (i) state this explicitly as a numerical bound (|α/M²| ≲ 10⁸ or similar, dominated by the BH-existence prior, not by S2), and (ii) re-examine the Q∈(0,0.738) bound at the same rigor: the (6−Q²) coefficient and the higher-order Q-dependent terms together yield this 1σ interval when α is marginalized within the BH-existence region; we will provide the explicit χ² profile and replace Fig. 9 with a properly marginalized contour. We agree this strengthens, not weakens, the referee's point that S2 is essentially a Q-only probe. revision: yes
-
Referee: §II A / Fig. 1: the A/B region boundaries are stated without derivation; the inner discriminant 1+4α(2M/r³−Q²/r⁴) inside f(r) may produce branch-cut surfaces or curvature singularities outside r_h for α<0; please confirm that the negative-α orbits in Figs. 12–15 lie entirely in the region where the metric is real and no additional structure affects the geodesic.
Authors: We thank the referee. The boundaries of regions A and B in Fig. 1 follow from requiring (i) that f(r_h)=0 has solutions with r_h>0 and (ii) the reality of the inner square root at the horizon; explicitly, M²−Q²−α ≥ 0 gives the outer envelope and the condition that the inner discriminant 1+4α(2M/r³−Q²/r⁴) be non-negative at the relevant root yields the inner boundaries quoted. We will add a short derivation in §II A. Regarding the negative-α case: the inner discriminant defines a critical radius r_* below which the metric becomes complex, signalling a branch-cut/curvature singularity outside r_h that is characteristic of the 4D-EGB negative branch (see also Refs. [31, 54]). For all orbits shown in Figs. 12–15 with α<0 (specifically Fig. 14, α/M² = −0.5, Q/M = 0.5), we have verified numerically that the pericenter r_p of every plotted (z,w,v) orbit satisfies r_p > r_*, i.e., the entire trajectory lies in the region where g_{μν} is real and no additional horizon or singular surface is crossed. We will add this check, together with the explicit value of r_* for the parameters used, as a note in §IV and §V. We thank the referee for prompting this verification. revision: yes
- A quantitative mismatch/overlap calculation maximized over (M,E,L), as suggested by the referee for Major Comment 1, is not provided in the present revision. We argue that a credible overlap requires consistent adiabatic inspiral evolution (energy/angular-momentum fluxes for the 4D-EGB charged BH background), which is beyond the test-particle scope of this paper, and we have therefore reframed the distinguishability claim rather than producing the diagnostic. The referee may reasonably regard this as an incomplete answer to the headline objection.
Circularity Check
No significant circularity: standard geodesic + Kludge-waveform pipeline applied to a published 4D-EGB metric; self-citations are illustrative, not load-bearing.
specific steps
-
self citation load bearing
[Sec. II.B, end of ISCO discussion; refs [48, 58]]
"the RNBH (Q != 0, alpha = 0) has smaller r_ISCO, L_ISCO, and E_ISCO compared to the Schwarzschild BH (Q = 0, alpha = 0), which is consistent with the conclusions in related literature[48, 58]."
Reference [48] is by an overlapping author group (Zhao, Tang, Xu). However the citation is used only as a consistency check on a well-known RN-vs-Schwarzschild ISCO ordering that is independently verifiable from the metric, so it is not load-bearing for any numerical result in the present paper. Listed for completeness; impact on score is minor.
full rationale
The derivation chain is self-contained against external inputs: (1) the charged 4D-EGB metric f(r) is taken from Fernandes [39] (independent author) and its mathematical legitimacy is supported by independent regularization papers [27-33]; (2) geodesics, V_eff, MBO/ISCO conditions follow from the standard squared Lagrangian with no fitted parameters; (3) the periodic-orbit taxonomy q = w + v/z is imported verbatim from Levin & Perez-Giz [47] (external); (4) shadow-radius and S2-precession constraints are compared against externally published EHT bounds [63-65] and GRAVITY data [45], not against quantities defined inside the paper; (5) the Kludge quadrupole waveform is the standard Babak et al. [68] prescription. No quantity is fitted to data and then re-presented as a prediction. The two self-citations [48, 51] are used only to corroborate qualitative ISCO/periodic-orbit trends and are not load-bearing for any numerical claim. The skeptic's concern (parameter degeneracy between (alpha, Q) and (M, E, L) in actual EMRI inference) is a correctness/sufficiency issue about whether the visual waveform differences in Figs. 14-15 survive joint parameter estimation, not a circularity issue: nothing in the paper's chain reduces a 'prediction' to its own input by construction. The authors themselves acknowledge the deferred degeneracy analysis in the final paragraph, which is a scope limitation rather than circular reasoning. Score is 1 rather than 0 to acknowledge routine self-citation to [48, 51].
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
IndisputableMonolith/Constants.leanConstantDerivations.G_rs_eq, kappa_einstein_eq unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
f(r) = 1 + r^2/(2α)[1 ± √(1 + 4α(2M/r^3 − Q^2/r^4))]
-
IndisputableMonolith/Foundation/PhiForcing.leanphi_forced unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
we employ the periodic orbit classification method proposed by Janna Levin et al. ... q = ∆ϕ/(2π) − 1 = w + ν/z
-
IndisputableMonolith/Unification/SpacetimeEmergence.leanspacetime_emergence_cert unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
we investigate ... gravitational waveforms ... finding that variations in α and Q can lead to distinguishable differences in periodic orbit structures and gravitational wave phases
-
IndisputableMonolith/Constants.leanc_rs_eq_one, hbar_eq_phi_inv_fifth echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
this study adopts natural units where c = G = 1
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.
Forward citations
Cited by 1 Pith paper
-
Thin Accretion Disks around Rotating Charged Black Holes in an Effective Higher-Curvature Spacetime
An effective Gauss-Bonnet-like deformation of the Kerr-Newman metric moves the ISCO inward, raises radiative efficiency, flux, and temperature, while charge suppresses these quantities.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.