Pith. sign in

REVIEW 3 major objections 5 minor 74 references

The paper constructs leading-order spacetimes for rotating black holes carrying massive scalar fields in three modified-gravity theories, with scalar residuals near 10^-5 and horizon spin and temperature shifts computed from them.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 01:05 UTC pith:5B6HU6AO

load-bearing objection First massive-scalar-hair Kerr spacetimes in three modified-gravity theories, built on a coherent spectral extension; the main risk is an unproven metric-ansatz completeness, and the accuracy claims need restating. the 3 major comments →

arxiv 2602.10462 v3 pith:5B6HU6AO submitted 2026-02-11 gr-qc

Spacetime of rotating black holes surrounded by massive scalar charges

Adrian Ka-Wai Chung This is my paper
classification gr-qc PACS 04.70.Bw04.50.Kd04.30.-w
keywords rotating black holesmassive scalar fieldsscalar Gauss-Bonnet gravitydynamical Chern-Simons gravityaxi-dilaton gravityspectral methodsblack hole hairgravitational waves
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Working at leading order in a small coupling that measures the strength of the curvature-scalar interaction, the paper solves the massive Klein-Gordon equation and the linearized Einstein equations on a fixed Kerr background. The key move is to factor out the exponential decay of the massive field analytically and treat the remaining smooth function with a spectral expansion; this yields solutions with residuals near 10^-5 for scalar fields with masses up to 0.2 in units of the black-hole mass, and near 10^-3 for the metric, for dimensionless spins up to 0.8. The resulting spacetimes are used to compute the leading-order changes in the horizon angular velocity and surface gravity—the quantities needed to predict how black-hole ringdown and inspiral signals would be altered by such scalar hair. A sympathetic reader would care because massive scalar fields are well-motivated dark-matter candidates, and these spacetimes connect them to observable gravitational-wave and electromagnetic signatures.

Core claim

The central claim is that the spacetime of a rotating black hole surrounded by a massive scalar charge can be constructed at leading order in the coupling without expanding in spin, for scalar Gauss-Bonnet, dynamical Chern-Simons, and axi-dilaton gravity. The author reports that the spectral scheme resolves scalar fields with Compton wavelengths as short as five times the black-hole mass (μM ≤ 0.2) with residual errors ≲10^-5, and yields the corresponding leading-order metric deformations with residual errors ≲10^-3, for dimensionless spins a ≤ 0.8. The constructed spacetimes are then used to compute the leading-order shifts in the horizon angular velocity and surface gravity, which control

What carries the argument

The central object is the auxiliary field φ = e^{μr} ϑ, which removes the exponential Yukawa tail from the unknown; the exponential factor is kept attached to the differential operator, and the weighted integrals that appear in the spectral projection, I(i,j,k,l|ξ) = ∫ z^i T_k^{(j)}(z) T_l(z)/√(1-z²) exp(−ξ/(1+z)) dz, are evaluated numerically at high precision. Together with a compactified radial coordinate z = 2r+/r − 1 and a Chebyshev-Legendre basis, this isolates the stiff radial behavior of the massive field and turns the Klein-Gordon equation and the six independent linearized Einstein equations into augmented linear systems that are solved by least squares. The metric is captured by t

Load-bearing premise

The load-bearing assumption is that the four-function metric ansatz of Eq. (24) spans the entire leading-order deformation of Kerr induced by the massive scalar field; any sourced stationary, axisymmetric component outside this ansatz, such as a g_{rχ} term or a different parity combination, would be missed, and then the constructed spacetime would not be the true solution.

What would settle it

Solve the same leading-order equations with a fully general stationary, axisymmetric metric ansatz and check whether the components outside Eq. (24) are sourced by the massive scalar field; if any of them develops a nonzero source, the four-function ansatz is incomplete. An independent check at a single point (e.g., a=0.8, μ=0.2) using a different numerical scheme, such as finite differences or time-domain evolution, should reproduce the reported shifts in Ω_H and κ to within the claimed residuals.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The computed horizon quantities enter standard quasinormal-mode frameworks, so the leading-order effect of massive scalar hair on ringdown frequencies can now be calculated for spins up to 0.8.
  • The metric perturbations provide a starting point for post-Newtonian and post-Minkowskian inspiral models, extending existing searches for scalar charges to massive fields in dynamical Chern-Simons and axi-dilaton gravity.
  • The result that mass affects magnitude more than multipole structure suggests observational templates built for massless scalars can be rescaled for a range of masses, a simplification for data analysis.
  • The spacetimes can serve as backgrounds for superradiance studies, since stationary massive scalar configurations are natural seeds for superradiant instabilities.
  • Axi-dilaton spacetimes are obtained by superposition of the dynamical Chern-Simons and scalar Gauss-Bonnet pieces, so a single computation covers three theories.

Where Pith is reading between the lines

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

  • A direct test of completeness would be to free up the metric ansatz—for instance adding an independent g_{rχ} function—and check whether its equation is identically sourced at leading order by the massive scalar; the paper verifies residuals only within the four-function ansatz.
  • The spectral instability at larger masses (N_min scaling like 1/(μ r_+)) suggests that basis functions with an exponential weight, such as generalized Laguerre polynomials, may extend the method beyond μM = 0.2, which the paper itself suggests.
  • The same machinery could be recycled for vector or tensor ultralight fields, whose asymptotic behavior is also Yukawa-like, although the parity and tensor structure of the sources would differ.
  • The metric residual grows with spin to about 10^-3 at a = 0.8, implying that practical forecasts for gravitational-wave constraints should marginalize over this systematic error when translating residuals into physical bounds.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper extends spectral methods previously developed for massless scalar fields to construct the leading-order (in the coupling parameter) massive scalar field configurations and metric deformations around Kerr black holes in dynamical Chern-Simons (dCS), scalar-Gauss-Bonnet (sGB), and axi-dilaton gravity. The scalar field is factored as e^{-μr}φ and solved with Chebyshev/Legendre spectral expansions on a compactified radial coordinate; the metric perturbations are restricted to the four-function ansatz (24), and the resulting linearized Einstein equations are solved for H1..H4. The paper reports residual errors of order 10^{-5} for the scalar and 10^{-3} for the metric, computes the leading-order shifts in horizon angular velocity Ω_H^{(1)} and surface gravity κ^{(1)}, and argues that adding a scalar-field mass leaves the multipolar structure qualitatively unchanged while modifying the overall magnitude. All numerical results are available as Supplementary Material.

Significance. If correct, this is the first systematic numerical construction of leading-order scalar hair and metric deformations for rotating Kerr black holes with massive scalar fields in these three beyond-GR theories. The exponential-factor treatment in the spectral scheme is a useful technical contribution, and the analytic evaluation of residual integrals (Eqs. 50-52, 69-70) is commendable. The resulting Ω_H^{(1)} and κ^{(1)} are potentially important inputs for quasinormal-mode and gravitational-wave searches. The paper is clear about its parameter range (a≤0.8, μ≤0.2/M) and provides convergence diagnostics, which are strengths. However, the accuracy claims rest on normalized residuals and the assumed completeness of the four-function metric ansatz, neither of which is fully established.

major comments (3)
  1. [§II.C / §IV.A, Eq. (24)] The four-function ansatz (24) is assumed to span the entire leading-order stationary, axisymmetric metric deformation. It sets g_{tr}, g_{tχ}, g_{rφ}, g_{χφ}, and g_{rχ} to zero and identifies δg_rr with δg_χχ through one function H3. The paper does not demonstrate that the trace-reversed stress tensor (9) and A^{(0)} in Eq. (14) satisfy the circularity conditions that make this form general, nor that the omitted components are pure gauge in the massive case. The rχ component is solved within the ansatz, and the text in §IV.A that this component is odd and needs odd Legendre polynomials shows it is nontrivial. The residual (69) is evaluated only for metrics of the form (24), so it cannot detect a missing g_{rχ}^{(1)}. This is load-bearing: if a nonzero g_{rχ}^{(1)} (or another omitted component) is sourced, the constructed spacetime and the derived Ω_H and κ shifts are not the physical l
  2. [§III.C / §IV.C, Eqs. (50), (69); abstract/conclusion] The residual errors quoted in the abstract are normalized so that E(N=1)=1 for each (a,μ), as stated in §III.C.2 and §IV.C.2. Thus the numbers '10^{-5}' and '10^{-3}' are relative convergence measures, not absolute residuals. The abstract's wording 'achieving residual errors ≲10^{-5}' and similar statements are misleading. Moreover, the body reports E(Nopt) ∼10^{-3} at a=0.8 for metric modifications (§IV.C.2, Fig. 7), while the conclusion says 'error ≲10^{-4}'. Please report an unnormalized or gauge-invariant error, or rephrase the abstract/conclusion, and reconcile the 10^{-3}/10^{-4} discrepancy.
  3. [§V, Eqs. (73)-(76), Figs. 11 and 13] The advertised physical outputs are Ω_H^{(1)} and κ^{(1)}. The horizon-constancy diagnostics show ||dκ^{(1)}/dχ||_2 ∼10^{-2} for sGB at a=0.8 (Fig. 13), which is roughly 10% of κ^{(1)}∼0.1. The paper states this, but the abstract and conclusion present the construction as accurate without noting this 10% uncertainty in the key observables at high spin. Please quantify the error in Ω_H^{(1)} and κ^{(1)} (e.g., variation with N around Nopt, or comparison with an extended ansatz) and state the limitation in the abstract.
minor comments (5)
  1. [§II.B, Eq. (23)] Since Y0 diverges as the argument vanishes, physical solutions demand C2=0, not C1=0 as written.
  2. [§IV.A, Eq. (58)] The text says the functions Hi are even in χ but then retains both even and odd Legendre polynomials in the expansion. Clarify whether the odd polynomials are used only as projection/test functions, or whether odd modes are included in the solution and how parity is enforced or checked.
  3. [Throughout] Grammar/typos: 'we computes' in the abstract; 'ingradients' in §V; 'diamons' in Fig. 11 caption; 'Identification to Fig.' in several figure captions; 'to future work to future work' in §VI.
  4. [§III.A / §IV.B] The subtraction of 'bare terms' (Eqs. 45-46, 66-67) is an ad hoc part of the algorithm; please justify that it does not absorb physical r^{-1} components and that the resulting solution still converges to the same continuum solution.
  5. [Data availability] The statement that numerical results are 'available from the author upon reasonable request' is less reproducible than depositing the data/code in a public repository; consider making the Supplementary Material self-contained.

Circularity Check

0 steps flagged

No circularity found: horizon observables are read off an independently solved linearized system; self-cited ansatz is an assumption, not a reduction.

full rationale

The derivation chain is self-contained in the relevant sense. The paper (i) solves the linearized Klein-Gordon equation (13) on the fixed Kerr background for the rescaled field phi = e^{mu r} bar-vartheta via spectral projection (39)-(44); (ii) substitutes that field into the modified Einstein equations (14)/(55) and solves the linear algebraic system (59)-(65) for the metric functions H_i; and (iii) evaluates the leading-order horizon angular velocity and surface gravity from the solved H_i through Eqs. (73) and (75). No parameter is fitted to Omega_H^(1) or kappa^(1); those quantities are computed as functionals of the constructed metric, so they are not forced by construction. The residual measures (50) and (69) and the horizon-constancy checks (74) and (76) are independent diagnostics of the numerical solutions. The paper does cite the author's prior work [31,32] for the spectral framework and for the four-function metric ansatz (24), but the massive-scalar construction itself is new and is solved in this paper; no uniqueness theorem from the author's prior work is invoked to forbid alternative metric components. The residual in Eq. (69) is computed over the field equations within the retained ansatz, and the completeness of the four-function ansatz is an unproved modeling assumption — Sec. IV A even notes the need to retain odd Legendre polynomials for the (r,chi) equation. If that ansatz were incomplete, the constructed spacetime would not be the true leading-order solution; that is a correctness risk, not circularity, because the derivation does not assume the value of the predicted observables. The self-citations are present but not used to derive the target result, so the score is 1 rather than 0.

Axiom & Free-Parameter Ledger

4 free parameters · 6 axioms · 0 invented entities

No new particles, forces, or fields are introduced; ϑ1 and ϑ2 are from existing theories. The free parameters are exploratory choices and numerical convergence parameters, not fitted physical constants. The key load-bearing assumptions are the small-coupling expansion and the completeness of the four-function metric ansatz.

free parameters (4)
  • scalar field mass µ = 0.01, 0.1, 0.2 (M^-1)
    Explored values chosen by hand; not fitted, but determines the Compton-wavelength regime the method must resolve.
  • dimensionless spin a = 0 to 0.8
    Scan range; a=0.8 is the claimed upper limit of accuracy.
  • optimal spectral order Nopt = varies with (a, µ)
    Selected as argmin of the normalized residual Eϑ or E; a numerical convergence parameter, not a physical one.
  • residual normalization = E(N=1)=1
    Convention chosen to make relative convergence meaningful; absolute residual scale is arbitrary.
axioms (6)
  • domain assumption Kerr metric is the exact zeroth-order background
    Perturbative expansion in ζ around Kerr; justifies solving linearized equations on fixed Kerr (Eqs. (12)-(14)).
  • domain assumption Small-coupling truncation at O(ζ) captures the leading-order scalar hair and metric effect
    ζ from Eq. (5) is small; no O(ζ²) terms considered (Eqs. (10)-(14)).
  • ad hoc to paper Metric ansatz Eq. (24) is complete for the leading-order stationary axisymmetric deformation
    Four H_i functions assumed to span all sourced components; inherited from [21,31,32]; residual is only checked within this ansatz.
  • domain assumption The massive scalar field is strictly stationary and axisymmetric (∂t=∂ϕ=0)
    Klein-Gordon solved as elliptic problem without time dependence; superradiant growth neglected (Eq. (16)).
  • standard math Regularity at the horizon selects the physical solution
    Bessel-function analysis near r=r+; the physical branch must be finite. The text swaps C1/C2 in Eq. (23), but the intended regularity condition is standard.
  • domain assumption Asymptotic falloff ϑ ~ e^{-µr} φ with φ polynomial in 1/r
    Guides the spectral ansatz; removes exponential stiffness (Eqs. (18)-(19), (27)).

pith-pipeline@v1.3.0-alltime-deepseek · 28304 in / 16585 out tokens · 161401 ms · 2026-08-03T01:05:05.989212+00:00 · methodology

0 comments
read the original abstract

Massive scalar charges are ubiquitous in extensions to General Relativity and the Standard Model in particle physics. We describe spectral methods which can accurately construct the spacetime of rotating black holes with dimensionless spin up to $a \leq 0.8$ surrounded by massive scalar fields nonminimally coupled to spacetime curvature. We consider axi dilaton, dynamical Chern Simons, and scalar Gauss Bonnet couplings, and obtain leading order solutions for both the scalar field and the associated metric modifications. Our method accurately resolves massive scalar fields with Compton wavelengths as short as 5 times the black hole mass, achieving residual errors $\lesssim 10^{-5}$, and yields the corresponding leading order spacetime modifications with residual errors $\lesssim 10^{-3}$. Using the constructed spacetimes, we computes the leading-order shifts in the surface gravity and the angular velocity of the event horizon, important information for computing the quasinormal modes. These results pave the way to incorporate massive scalar charges into electromagnetic observations and gravitational-wave detections of black holes, potentially enabling new probes of fundamental scalar degrees of freedom.

Figures

Figures reproduced from arXiv: 2602.10462 by Adrian Ka-Wai Chung.

Figure 1
Figure 1. Figure 1: FIG. 1. The backward-modulus difference [see Eq. [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The error [see Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The least error [see the text around Eq. [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Meridional cross-sections of the massive scalar field around a rotating black hole with dimensionless spin [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The backward-modulus difference [see Eq. [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The error [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. The least residual of the modified Einstein equations [see Eq. [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Meridional cross-sections of [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Identical to Fig [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. The leading-order modification to the angular velocity of the event horizon, Ω [PITH_FULL_IMAGE:figures/full_fig_p019_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11 [PITH_FULL_IMAGE:figures/full_fig_p019_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Identification to Fig [PITH_FULL_IMAGE:figures/full_fig_p020_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Identification to Fig [PITH_FULL_IMAGE:figures/full_fig_p021_13.png] view at source ↗

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

74 extracted references · 55 linked inside Pith

  1. [1]

    distance

    Backward modulus difference To assess the convergence of the scalar field ϑ with increasing spectral order N , we define the backward mod- ulus difference (BMD) as Bϑ(N) = "Z ∞ r+ Z +1 −1 [φ(N)−φ(N−1)] 2 drdχ #1/2 . (47) The exponential factor e−µr is deliberately excluded from this definition, since it deforms φ(N ) and φ(N− 1) in the same manner and the...

  2. [2]

    [31, 32]: Eϑ ∝ "Z ∞ r+ Z +1 −1 E2 ϑ h −g(0) i6 dr r26 dχ # 1 2 .(50) If ϑ were an exact solution of the Klein–Gordon equation, then Eϑ would vanish identically

    Absolute error To quantify how well the scalar fieldϑ constructed using the spectral scheme satisfies the Klein–Gordon equation, Eϑ = 0, we define the following absolute error, based on the definition in Refs. [31, 32]: Eϑ ∝ "Z ∞ r+ Z +1 −1 E2 ϑ h −g(0) i6 dr r26 dχ # 1 2 .(50) If ϑ were an exact solution of the Klein–Gordon equation, then Eϑ would vanish...

  3. [3]

    3 shows the least error, Eϑ(Nopt), as a function of the dimensionless spin a for µ = 0.01, 0.1, and 0 .2 in dCS (left panel) and sGB (right panel) gravity

    Results fora≤0.8 Fig. 3 shows the least error, Eϑ(Nopt), as a function of the dimensionless spin a for µ = 0.01, 0.1, and 0 .2 in dCS (left panel) and sGB (right panel) gravity. We find that the least error exhibits no significant dependence on the spin parameter a, but increases noticeably as the scalar-field mass is raised from µ = 0.1 to µ = 0.2. This ...

  4. [4]

    The Sloan Lens ACS Survey. 4. The mass density profile of early-type galaxies out to 100 effective radii,

    R. Gavazzi, T. Treu, J. D. Rhodes, L. V. Koopmans, A. S. Bolton, S. Burles, R. Massey, and L. A. Moustakas, “The Sloan Lens ACS Survey. 4. The mass density profile of early-type galaxies out to 100 effective radii,” Astrophys. J.667, 176–190 (2007), arXiv:astro-ph/0701589

    Pith/arXiv arXiv 2007

  5. [5]

    (68) We include a factor of r−2 as it is the simplest regulariza- tion factor such that B(N ) is finite and well defined, as Hi can approach to a constant asr→+∞

    Backward modulus difference To gauge the convergence of the metric modifications, we define the backward modulus difference B(N ) that the changes of Hi(r, χ) computed at a given spectral N from that computed at the previous spectral order N− 1 as follows B(N) = "Z +∞ r+ Z +1 −1 4X i=1 [Hi(N)−H i(N−1)] 2 dr r2 dχ # 1 2 . (68) We include a factor of r−2 as...

  6. [6]

    Absolute error To quantify how well the metric perturbations Hi and the scalar field ¯ϑ satisfy the modified Einstein equations, we define the following error, E ∝ "Z ∞ r+ Z +1 −1 Eβ νEν β ∆4 r32 (1−χ 2)2 −g(0) 6 drdχ #1/2 (69) The factor ∆ 4(1 −χ 2)2 −g(0) 6 /r32 is the simplest reg- ularization weight we have identified that renders the integral finite ...

  7. [7]

    Rotation curves of spiral galax- ies,

    Y. Sofue and V. Rubin, “Rotation curves of spiral galax- ies,” Ann. Rev. Astron. Astrophys.39, 137–174 (2001), arXiv:astro-ph/0010594

    Pith/arXiv arXiv 2001

  8. [8]

    History of dark matter,

    G. Bertone and D. Hooper, “History of dark matter,” Rev. Mod. Phys.90, 045002 (2018), arXiv:1605.04909 [astro-ph.CO]

    Pith/arXiv arXiv 2018

  9. [9]

    The dark matter of gravitational lensing,

    R. Massey, T. Kitching, and J. Richard, “The dark matter of gravitational lensing,” Rept. Prog. Phys.73, 086901 (2010), arXiv:1001.1739 [astro-ph.CO]

    Pith/arXiv arXiv 2010

  10. [10]

    Dilaton-axion inflation with PBHs and GWs,

    R. Kallosh and A. Linde, “Dilaton-axion inflation with PBHs and GWs,” JCAP08, 037 (2022), arXiv:2203.10437 [hep-th]

    Pith/arXiv arXiv 2022

  11. [11]

    Virial masses of galactic halos from galaxy-galaxy lensing: Theoretical modeling and application to SDSS,

    J. Guzik and U. Seljak, “Virial masses of galactic halos from galaxy-galaxy lensing: Theoretical modeling and application to SDSS,” Mon. Not. Roy. Astron. Soc.335, 311 (2002), arXiv:astro-ph/0201448

    Pith/arXiv arXiv 2002

  12. [12]

    Galaxy halo masses and satellite fractions from galaxy-galaxy lensing in the sdss: stellar mass, luminosity, morphology, and environment depen- dencies,

    R. Mandelbaum, U. Seljak, G. Kauffmann, C. M. Hirata, and J. Brinkmann, “Galaxy halo masses and satellite fractions from galaxy-galaxy lensing in the sdss: stellar mass, luminosity, morphology, and environment depen- dencies,” Mon. Not. Roy. Astron. Soc.368, 715 (2006), arXiv:astro-ph/0511164

    Pith/arXiv arXiv 2006

  13. [13]

    Measurements of Ω and Λ from 42 high redshift super- novae,

    S. Perlmutteret al.(Supernova Cosmology Project), “Measurements of Ω and Λ from 42 high redshift super- novae,” Astrophys. J.517, 565–586 (1999), arXiv:astro- ph/9812133

    arXiv 1999

  14. [14]

    Obser- vational evidence from supernovae for an accelerating universe and a cosmological constant,

    A. G. Riesset al.(Supernova Search Team), “Obser- vational evidence from supernovae for an accelerating universe and a cosmological constant,” Astron. J.116, 1009–1038 (1998), arXiv:astro-ph/9805201

    Pith/arXiv arXiv 1998

  15. [15]

    Baryon asym- metry of the universe in the minimal Standard Model,

    G. R. Farrar and M. E. Shaposhnikov, “Baryon asym- metry of the universe in the minimal Standard Model,” Phys. Rev. Lett.70, 2833–2836 (1993), [Erratum: Phys.Rev.Lett. 71, 210 (1993)], arXiv:hep-ph/9305274

    Pith/arXiv arXiv 1993

  16. [16]

    Im- proved gravitational-wave constraints on higher-order cur- vature theories of gravity,

    S. E. Perkins, R. Nair, H. O. Silva, and N. Yunes, “Im- proved gravitational-wave constraints on higher-order cur- vature theories of gravity,” Phys. Rev. D104, 024060 (2021), arXiv:2104.11189 [gr-qc]

    Pith/arXiv arXiv 2021

  17. [17]

    String gravity in D=4,

    P. A. Cano and A. Ruip´ erez, “String gravity in D=4,” Phys. Rev. D105, 044022 (2022), arXiv:2111.04750 [hep- th]

    Pith/arXiv arXiv 2022

  18. [18]

    Chern-Simons modification of general relativity,

    R. Jackiw and S. Y. Pi, “Chern-Simons modification of general relativity,” Phys. Rev. D68, 104012 (2003), arXiv:gr-qc/0308071

    Pith/arXiv arXiv 2003

  19. [19]

    Chern-Simons Modified General Relativity,

    S. Alexander and N. Yunes, “Chern-Simons Modified General Relativity,” Phys. Rept.480, 1–55 (2009), arXiv:0907.2562 [hep-th]

    Pith/arXiv arXiv 2009

  20. [20]

    Challenging the Presence of Scalar Charge and Dipolar Radiation in Binary Pulsars,

    K. Yagi, L. C. Stein, and N. Yunes, “Challenging the Presence of Scalar Charge and Dipolar Radiation in Binary Pulsars,” Phys. Rev. D93, 024010 (2016), arXiv:1510.02152 [gr-qc]

    Pith/arXiv arXiv 2016

  21. [21]

    Funda- mental Physics Implications for Higher-Curvature Theo- ries from Binary Black Hole Signals in the LIGO-Virgo Catalog GWTC-1,

    R. Nair, S. Perkins, H. O. Silva, and N. Yunes, “Funda- mental Physics Implications for Higher-Curvature Theo- ries from Binary Black Hole Signals in the LIGO-Virgo Catalog GWTC-1,” Phys. Rev. Lett.123, 191101 (2019), arXiv:1905.00870 [gr-qc]

    Pith/arXiv arXiv 2019

  22. [22]

    Superstring modifications of einstein’s equations,

    D. J. Gross and E. Witten, “Superstring modifications of einstein’s equations,” Nuclear Physics B277, 1–10 (1986)

    1986

  23. [23]

    Gravitational collapse in Einstein dilaton-Gauss–Bonnet gravity,

    J. L. Ripley and F. Pretorius, “Gravitational collapse in Einstein dilaton-Gauss–Bonnet gravity,” Class. Quant. Grav.36, 134001 (2019), arXiv:1903.07543 [gr-qc]

    Pith/arXiv arXiv 2019

  24. [24]

    Inspiral- merger-ringdown waveforms in Einstein-scalar-Gauss- Bonnet gravity within the effective-one-body formalism,

    F.-L. Juli´ e, L. Pompili, and A. Buonanno, “Inspiral- merger-ringdown waveforms in Einstein-scalar-Gauss- Bonnet gravity within the effective-one-body formalism,” Phys. Rev. D111, 024016 (2025), arXiv:2406.13654 [gr- qc]

    Pith/arXiv arXiv 2025

  25. [25]

    Post-Newtonian dynamics and black hole thermodynamics in Einstein-scalar-Gauss- Bonnet gravity,

    F.-L. Juli´ e and E. Berti, “Post-Newtonian dynamics and black hole thermodynamics in Einstein-scalar-Gauss- Bonnet gravity,” Phys. Rev. D100, 104061 (2019), arXiv:1909.05258 [gr-qc]

    Pith/arXiv arXiv 2019

  26. [26]

    Tensor multiscalar theories of gravitation,

    T. Damour and G. Esposito-Farese, “Tensor multiscalar theories of gravitation,” Class. Quant. Grav.9, 2093–2176 (1992)

    2093

  27. [27]

    Leading higher-derivative corrections to kerr geometry,

    P. A. Cano and A. Ruip´ erez, “Leading higher-derivative corrections to kerr geometry,” Journal of High Energy Physics2019, 189 (2019), arXiv: 1901.01315

    Pith/arXiv arXiv 2019

  28. [28]

    Are black holes in alternative theories serious astrophysical candidates? The Case for Einstein-Dilaton-Gauss-Bonnet black holes,

    P. Pani and V. Cardoso, “Are black holes in alternative theories serious astrophysical candidates? The Case for Einstein-Dilaton-Gauss-Bonnet black holes,” Phys. Rev. D79, 084031 (2009), arXiv:0902.1569 [gr-qc]

    Pith/arXiv arXiv 2009

  29. [29]

    Heterotic String Theory. 2. The Interacting Heterotic String,

    D. J. Gross, J. A. Harvey, E. J. Martinec, and R. Rohm, “Heterotic String Theory. 2. The Interacting Heterotic String,” Nucl. Phys. B267, 75–124 (1986)

    1986

  30. [30]

    Can the string scale be related to the cosmic baryon asymmetry?

    S. H. S. Alexander and S. J. Gates, Jr., “Can the string scale be related to the cosmic baryon asymmetry?” JCAP 06, 018 (2006), arXiv:hep-th/0409014

    Pith/arXiv arXiv 2006

  31. [31]

    Leptogenesis from gravity waves in models of inflation,

    S. H. S. Alexander, M. E. Peskin, and M. M. Sheikh- Jabbari, “Leptogenesis from gravity waves in models of inflation,” Phys. Rev. Lett.96, 081301 (2006)

    2006

  32. [32]

    Dilatonic black holes in higher curva- ture string gravity,

    P. Kanti, N. E. Mavromatos, J. Rizos, K. Tamvakis, and E. Winstanley, “Dilatonic black holes in higher curva- ture string gravity,” Phys. Rev. D54, 5049–5058 (1996), arXiv:hep-th/9511071

    Pith/arXiv arXiv 1996

  33. [33]

    Dynamical Chern-Simons Modified Gravity. I. Spinning Black Holes in the Slow- Rotation Approximation,

    N. Yunes and F. Pretorius, “Dynamical Chern-Simons Modified Gravity. I. Spinning Black Holes in the Slow- Rotation Approximation,” Phys. Rev. D79, 084043 (2009), arXiv:0902.4669 [gr-qc]

    Pith/arXiv arXiv 2009

  34. [34]

    Axion Cosmology,

    D. J. E. Marsh, “Axion Cosmology,” Phys. Rept.643, 1–79 (2016), arXiv:1510.07633 [astro-ph.CO]

    Pith/arXiv arXiv 2016

  35. [35]

    Numerical black hole solutions in modified gravity theories: Spher- ical symmetry case,

    A. Sullivan, N. Yunes, and T. P. Sotiriou, “Numerical black hole solutions in modified gravity theories: Spher- ical symmetry case,” Phys. Rev. D101, 044024 (2020), arXiv:1903.02624 [gr-qc]

    Pith/arXiv arXiv 2020

  36. [36]

    Numerical black hole solutions in modified gravity theories: Ax- ial symmetry case,

    A. Sullivan, N. Yunes, and T. P. Sotiriou, “Numerical black hole solutions in modified gravity theories: Ax- ial symmetry case,” Phys. Rev. D103, 124058 (2021), arXiv:2009.10614 [gr-qc]

    Pith/arXiv arXiv 2021

  37. [37]

    Near- Extremal Black Holes in Modified Gravity via Spectral Methods,

    K. K.-H. Lam, A. K.-W. Chung, and N. Yunes, “Near- Extremal Black Holes in Modified Gravity via Spectral Methods,” (2025), arXiv:2509.07061 [gr-qc]

    arXiv 2025

  38. [38]

    Ana- lytic and accurate approximate metrics for black holes with arbitrary rotation in beyond-Einstein gravity using spectral methods,

    K. K.-H. Lam, A. K.-W. Chung, and N. Yunes, “Ana- lytic and accurate approximate metrics for black holes with arbitrary rotation in beyond-Einstein gravity using spectral methods,” (2025), arXiv:2510.05208 [gr-qc]

    arXiv 2025

  39. [39]

    Cold and fuzzy dark matter,

    W. Hu, R. Barkana, and A. Gruzinov, “Cold and fuzzy dark matter,” Phys. Rev. Lett.85, 1158–1161 (2000), arXiv:astro-ph/0003365. 23

    Pith/arXiv arXiv 2000

  40. [40]

    Black holes as particle detectors: evolution of superradiant instabilities,

    R. Brito, V. Cardoso, and P. Pani, “Black holes as particle detectors: evolution of superradiant instabilities,” Class. Quant. Grav.32, 134001 (2015), arXiv:1411.0686 [gr-qc]

    Pith/arXiv arXiv 2015

  41. [41]

    Ultralight scalars as cosmological dark matter,

    L. Hui, J. P. Ostriker, S. Tremaine, and E. Witten, “Ultralight scalars as cosmological dark matter,” Phys. Rev. D95, 043541 (2017), arXiv:1610.08297 [astro-ph.CO]

    Pith/arXiv arXiv 2017

  42. [42]

    Amplification of waves during reflec- tion from a rotating

    A. A. Starobinskii, “Amplification of waves during reflec- tion from a rotating ”black hole”,” Sov. Phys. JETP37, 28–32 (1973)

    1973

  43. [43]

    Superradiance: New Frontiers in Black Hole Physics,

    R. Brito, V. Cardoso, and P. Pani, “Superradiance: New Frontiers in Black Hole Physics,” Lect. Notes Phys.906, pp.1–237 (2015), arXiv:1501.06570 [gr-qc]

    Pith/arXiv arXiv 2015

  44. [44]

    Floating and sinking: The Imprint of mas- sive scalars around rotating black holes,

    V. Cardoso, S. Chakrabarti, P. Pani, E. Berti, and L. Gualtieri, “Floating and sinking: The Imprint of mas- sive scalars around rotating black holes,” Phys. Rev. Lett. 107, 241101 (2011), arXiv:1109.6021 [gr-qc]

    Pith/arXiv arXiv 2011

  45. [45]

    Gravitational Waves from Quasicircular Extreme Mass-Ratio Inspirals as Probes of Scalar-Tensor Theories,

    N. Yunes, P. Pani, and V. Cardoso, “Gravitational Waves from Quasicircular Extreme Mass-Ratio Inspirals as Probes of Scalar-Tensor Theories,” Phys. Rev. D85, 102003 (2012), arXiv:1112.3351 [gr-qc]

    Pith/arXiv arXiv 2012

  46. [46]

    Or- bital fingerprints of ultralight scalar fields around black holes,

    M. C. Ferreira, C. F. B. Macedo, and V. Cardoso, “Or- bital fingerprints of ultralight scalar fields around black holes,” Phys. Rev. D96, 083017 (2017), arXiv:1710.00830 [gr-qc]

    Pith/arXiv arXiv 2017

  47. [47]

    Gravitational wave searches for ultralight bosons with LIGO and LISA,

    R. Brito, S. Ghosh, E. Barausse, E. Berti, V. Cardoso, I. Dvorkin, A. Klein, and P. Pani, “Gravitational wave searches for ultralight bosons with LIGO and LISA,” Phys. Rev. D96, 064050 (2017), arXiv:1706.06311 [gr-qc]

    Pith/arXiv arXiv 2017

  48. [48]

    Detecting scalar fields with Extreme Mass Ratio Inspirals,

    A. Maselli, N. Franchini, L. Gualtieri, and T. P. Sotiriou, “Detecting scalar fields with Extreme Mass Ratio Inspirals,” Phys. Rev. Lett.125, 141101 (2020), arXiv:2004.11895 [gr-qc]

    Pith/arXiv arXiv 2020

  49. [49]

    Detecting fundamental fields with LISA observations of gravitational waves from ex- treme mass-ratio inspirals,

    A. Maselli, N. Franchini, L. Gualtieri, T. P. Sotiriou, S. Barsanti, and P. Pani, “Detecting fundamental fields with LISA observations of gravitational waves from ex- treme mass-ratio inspirals,” Nature Astron.6, 464–470 (2022), arXiv:2106.11325 [gr-qc]

    Pith/arXiv arXiv 2022

  50. [50]

    Ultralight boson cloud depletion in binary systems,

    E. Berti, R. Brito, C. F. B. Macedo, G. Raposo, and J. L. Rosa, “Ultralight boson cloud depletion in binary systems,” Phys. Rev. D99, 104039 (2019), arXiv:1904.03131 [gr-qc]

    Pith/arXiv arXiv 2019

  51. [51]

    Gravitational floating orbits around hairy black holes,

    J. Zhang and H. Yang, “Gravitational floating orbits around hairy black holes,” Phys. Rev. D99, 064018 (2019), arXiv:1808.02905 [gr-qc]

    Pith/arXiv arXiv 2019

  52. [52]

    Quasinormal mode fre- quencies and gravitational perturbations of black holes with any subextremal spin in modified gravity through METRICS: The scalar-Gauss-Bonnet gravity case,

    A. K.-W. Chung and N. Yunes, “Quasinormal mode fre- quencies and gravitational perturbations of black holes with any subextremal spin in modified gravity through METRICS: The scalar-Gauss-Bonnet gravity case,” Phys. Rev. D110, 064019 (2024), arXiv:2406.11986 [gr-qc]

    Pith/arXiv arXiv 2024

  53. [53]

    Testing General Relativity with the Event Horizon Telescope,

    D. Psaltis, “Testing General Relativity with the Event Horizon Telescope,” Gen. Rel. Grav.51, 137 (2019), arXiv:1806.09740 [astro-ph.HE]

    Pith/arXiv arXiv 2019

  54. [54]

    Probing Ultralight Bosons with Binary Black Holes,

    D. Baumann, H. S. Chia, and R. A. Porto, “Probing Ultralight Bosons with Binary Black Holes,” Phys. Rev. D99, 044001 (2019), arXiv:1804.03208 [gr-qc]

    Pith/arXiv arXiv 2019

  55. [55]

    First Sagit- tarius A* Event Horizon Telescope Results. VI. Testing the Black Hole Metric,

    K. Akiyamaet al.(Event Horizon Telescope), “First Sagit- tarius A* Event Horizon Telescope Results. VI. Testing the Black Hole Metric,” Astrophys. J. Lett.930, L17 (2022), arXiv:2311.09484 [astro-ph.HE]

    Pith/arXiv arXiv 2022

  56. [56]

    The future ability to test theories of gravity with black-hole shadows,

    A. Uniyal, I. K. Dihingia, Y. Mizuno, and L. Rez- zolla, “The future ability to test theories of gravity with black-hole shadows,” Nature Astron.1, 8 (2025), arXiv:2511.03789 [gr-qc]

    arXiv 2025

  57. [57]

    Black hole superradiance in dynamical Chern-Simons grav- ity,

    S. Alexander, G. Gabadadze, L. Jenks, and N. Yunes, “Black hole superradiance in dynamical Chern-Simons grav- ity,” Phys. Rev. D107, 084016 (2023), arXiv:2201.02220 [gr-qc]

    Pith/arXiv arXiv 2023

  58. [58]

    R. P. Feynman, R. B. Leighton, and M. Sands,The Feynman Lectures on Physics(1963)

    1963

  59. [59]

    Quasi- normal mode frequencies and gravitational perturbations of spinning black holes in modified gravity through MET- RICS: The dynamical Chern-Simons gravity case,

    A. K.-W. Chung, K. K.-H. Lam, and N. Yunes, “Quasi- normal mode frequencies and gravitational perturbations of spinning black holes in modified gravity through MET- RICS: The dynamical Chern-Simons gravity case,” (2025), arXiv:2503.11759 [gr-qc]

    Pith/arXiv arXiv 2025

  60. [60]

    Evolution of Einstein-scalar- Gauss-Bonnet gravity using a modified harmonic formula- tion,

    W. E. East and J. L. Ripley, “Evolution of Einstein-scalar- Gauss-Bonnet gravity using a modified harmonic formula- tion,” Phys. Rev. D103, 044040 (2021), arXiv:2011.03547 [gr-qc]

    Pith/arXiv arXiv 2021

  61. [61]

    Probing quadratic gravity with black-hole ringdown gravitational waves measured by LIGO-Virgo-KAGRA detectors,

    A. K.-W. Chung and N. Yunes, “Probing quadratic gravity with black-hole ringdown gravitational waves measured by LIGO-Virgo-KAGRA detectors,” (2025), arXiv:2506.14695 [gr-qc]

    Pith/arXiv arXiv 2025

  62. [62]

    Bayesian Search of Massive Scalar Fields from LIGO- Virgo-KAGRA Binaries,

    Y. Xie, A. K.-W. Chung, T. P. Sotiriou, and N. Yunes, “Bayesian Search of Massive Scalar Fields from LIGO- Virgo-KAGRA Binaries,” Phys. Rev. Lett.134, 191402 (2025), arXiv:2410.14801 [gr-qc]

    Pith/arXiv arXiv 2025

  63. [63]

    Instability of the massive Klein-Gordon field on the Kerr spacetime,

    S. R. Dolan, “Instability of the massive Klein-Gordon field on the Kerr spacetime,” Phys. Rev. D76, 084001 (2007), arXiv:0705.2880 [gr-qc]

    Pith/arXiv arXiv 2007

  64. [64]

    Detecting Massive Scalar Fields with Extreme Mass- Ratio Inspirals,

    S. Barsanti, A. Maselli, T. P. Sotiriou, and L. Gualtieri, “Detecting Massive Scalar Fields with Extreme Mass- Ratio Inspirals,” Phys. Rev. Lett.131, 051401 (2023), arXiv:2212.03888 [gr-qc]

    Pith/arXiv arXiv 2023

  65. [65]

    On massive neu- tron cores,

    J. R. Oppenheimer and G. M. Volkoff, “On massive neu- tron cores,” Phys. Rev.55, 374–381 (1939)

    1939

  66. [66]

    Spectral method for metric perturbations of black holes: Kerr background case in general relativity,

    A. K.-W. Chung, P. Wagle, and N. Yunes, “Spectral method for metric perturbations of black holes: Kerr background case in general relativity,” Phys. Rev. D109, 044072 (2024), arXiv:2312.08435 [gr-qc]

    Pith/arXiv arXiv 2024

  67. [67]

    Poisson,A Relativist’s Toolkit: The Mathematics of Black-Hole Mechanics(Cambridge University Press, 2009)

    E. Poisson,A Relativist’s Toolkit: The Mathematics of Black-Hole Mechanics(Cambridge University Press, 2009)

    2009

  68. [68]

    The Four laws of black hole mechanics,

    J. M. Bardeen, B. Carter, and S. W. Hawking, “The Four laws of black hole mechanics,” Commun. Math. Phys.31, 161–170 (1973)

    1973

  69. [69]

    Leading effective field theory cor- rections to the Kerr metric at all spins,

    P. G. S. Fernandes, “Leading effective field theory cor- rections to the Kerr metric at all spins,” (2025), arXiv:2512.02338 [gr-qc]

    Pith/arXiv arXiv 2025

  70. [71]

    Slowly Rotating Black Holes in Dynamical Chern-Simons Gravity: De- formation Quadratic in the Spin,

    K. Yagi, N. Yunes, and T. Tanaka, “Slowly Rotating Black Holes in Dynamical Chern-Simons Gravity: De- formation Quadratic in the Spin,” Phys. Rev. D86, 044037 (2012), [Erratum: Phys.Rev.D 89, 049902 (2014)], arXiv:1206.6130 [gr-qc]

    Pith/arXiv arXiv 2012

  71. [72]

    Ringing out General Rel- ativity: Quasinormal Mode Frequencies for Black Holes of Any Spin in Modified Gravity,

    A. K.-W. Chung and N. Yunes, “Ringing out General Rel- ativity: Quasinormal Mode Frequencies for Black Holes of Any Spin in Modified Gravity,” Phys. Rev. Lett.133, 181401 (2024), arXiv:2405.12280 [gr-qc]

    Pith/arXiv arXiv 2024

  72. [73]

    Extremal black holes in dynamical Chern–Simons gravity,

    R. McNees, L. C. Stein, and N. Yunes, “Extremal black holes in dynamical Chern–Simons gravity,” Class. Quant. Grav.33, 235013 (2016), arXiv:1512.05453 [gr-qc]

    Pith/arXiv arXiv 2016

  73. [74]

    Note on Legendre decomposition of the Pontryagin density in Kerr,

    L. C. Stein, “Note on Legendre decomposition of the Pontryagin density in Kerr,” (2014), arXiv:1407.0744 [gr-qc]

    Pith/arXiv arXiv 2014

  74. [350]

    Recall that H1(r, χ) is required to satisfy a vanishing boundary condition at spatial infinity

    The inverse of ˜DT ˜D is computed using the built-in Mathematica function Inverse, with a working precision of 300. Recall that H1(r, χ) is required to satisfy a vanishing boundary condition at spatial infinity. In practice, how- ever, due to numerical truncation errors, the spectral solution for H1(r, χ) does not generally vanish in the limit r→ ∞. As in...