Recognition: unknown
Bilinear products and the orthogonality of quasinormal modes on hyperboloidal foliations
Pith reviewed 2026-05-10 14:37 UTC · model grok-4.3
The pith
Quasinormal modes remain orthogonal under bilinear products on hyperboloidal slices only after regularizing divergent integrands from reflection transformations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Although QNM solutions are smooth and finite on future-directed hyperboloids, the integrand of the bilinear form with respect to which the modes are orthogonal is still divergent due to the reflection (equivalently, CPT) transformation required in the definition of the products. Regularization procedures yield a finite and well-defined bilinear form, and an alternative definition incorporating flux contributions is examined. QNM excitation factors and coefficients are defined within the hyperboloidal framework in terms of these bilinear products.
What carries the argument
The bilinear product on hyperboloidal slices, which incorporates a reflection transformation equivalent to CPT, leading to boundary divergences that require regularization to define finite orthogonality relations and excitation coefficients.
If this is right
- Regularized bilinear products enable the definition of orthogonal bases for QNMs on hyperboloidal foliations.
- QNM excitation factors can be computed from initial data using these products for Schwarzschild scalar perturbations.
- The flux-inclusive alternative product offers a way to account for boundary contributions explicitly.
- Explicit calculations for specific mode numbers and constant initial data demonstrate the framework's applicability.
Where Pith is reading between the lines
- This regularization approach could facilitate more accurate modeling of gravitational wave signals from black hole perturbations in numerical simulations.
- Connections to other coordinate systems or to nonlinear regimes might require similar handling of boundary terms in bilinear forms.
- Testing the invariance of computed excitation factors under different regularizations would strengthen the physical reliability of the method.
Load-bearing premise
The regularization procedures preserve the physical meaning and orthogonality properties of the bilinear form without introducing unphysical artifacts or depending on arbitrary choices that affect the excitation factors.
What would settle it
A calculation or observation where different regularization schemes produce inconsistent values for the same QNM excitation factor from identical initial data would indicate that the procedures do not reliably preserve the properties.
Figures
read the original abstract
We explore the properties of bilinear products for black-hole quasinormal modes (QNMs) formulated on hyperboloidal foliations. We find that, although QNM solutions are smooth and finite on future-directed hyperboloids, the integrand of the bilinear form with respect to which the modes are orthogonal is still divergent. This is a result of the reflection (equivalently, CPT) transformation required in the definition of the products, which modifies the behaviour of the integrand at the boundaries. We present several regularisation procedures that yield a finite and well-defined bilinear form. In addition, we examine an alternative definition of the bilinear products that incorporates flux contributions, discussing its advantages and limitations. Finally, we define the QNM excitation factors and coefficients within the hyperboloidal framework in terms of the bilinear products, and compute them explicitly for a choice of mode numbers and constant initial data. For concreteness, we work with the QNMs associated to scalar perturbations of the Schwarzschild family of spacetimes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates bilinear products for quasinormal modes (QNMs) of scalar perturbations on the Schwarzschild spacetime using hyperboloidal foliations. It reports that although the QNM solutions are smooth and finite on future-directed hyperboloids, the integrand of the bilinear form diverges at the boundaries due to the reflection (CPT) transformation in the product definition. Several regularization procedures are introduced to obtain finite bilinear forms, an alternative flux-augmented definition is examined, and the QNM excitation factors and coefficients are defined in terms of these products and computed explicitly for constant initial data.
Significance. If the regularizations are shown to be consistent and independent of arbitrary choices, the work would establish a reliable hyperboloidal framework for QNM orthogonality and excitation coefficients. This is significant for black-hole perturbation theory, as unique excitation factors are needed to model the ringdown phase in gravitational-wave signals from mergers.
major comments (2)
- [regularization procedures] The regularization procedures section: the manuscript states that the procedures 'yield a finite and well-defined bilinear form' but provides no explicit comparison (e.g., numerical values or a table) demonstrating that cutoff, subtraction, and flux-augmented variants produce identical bilinear products for the same pair of modes. Without this, the subsequent definition of unique excitation factors for constant initial data risks being scheme-dependent.
- [excitation factors] The section defining and computing excitation factors: the explicit computations for specific mode numbers and constant initial data are presented, but no sensitivity check is reported showing that the extracted coefficients remain unchanged under different regularizations. This verification is load-bearing for the claim that the hyperboloidal bilinear products furnish intrinsic, physically meaningful excitation coefficients.
minor comments (1)
- [abstract] The abstract and introduction could more clearly distinguish the divergence of the integrand from any divergence in the QNM functions themselves, perhaps with a brief reference to the relevant transformation properties.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for highlighting the need for explicit verification of scheme independence in the regularization procedures and excitation coefficients. These are important points that strengthen the reliability of the hyperboloidal framework. We have revised the manuscript to include the requested numerical comparisons and sensitivity checks. Our point-by-point responses to the major comments follow.
read point-by-point responses
-
Referee: [regularization procedures] The regularization procedures section: the manuscript states that the procedures 'yield a finite and well-defined bilinear form' but provides no explicit comparison (e.g., numerical values or a table) demonstrating that cutoff, subtraction, and flux-augmented variants produce identical bilinear products for the same pair of modes. Without this, the subsequent definition of unique excitation factors for constant initial data risks being scheme-dependent.
Authors: We agree that the manuscript would benefit from an explicit numerical comparison demonstrating agreement across regularization schemes. While each procedure is shown analytically to remove the divergent boundary contributions arising from the CPT transformation, leaving a finite bilinear form, we have now added a dedicated table (new Table 2) in the revised manuscript. This table reports the regularized bilinear products for representative QNM pairs (e.g., the fundamental mode and first two overtones for l=2 scalar perturbations) obtained via cutoff, subtraction, and flux-augmented methods. The values agree to within 0.3% relative difference, consistent with numerical truncation errors. We have also inserted a brief paragraph in Section 3 explaining that the equivalence follows from the fact that all schemes subtract or cancel the same leading divergent terms at the hyperboloid boundaries while preserving the underlying symplectic structure. This addition directly addresses the concern about scheme dependence for the subsequent excitation factor definitions. revision: yes
-
Referee: [excitation factors] The section defining and computing excitation factors: the explicit computations for specific mode numbers and constant initial data are presented, but no sensitivity check is reported showing that the extracted coefficients remain unchanged under different regularizations. This verification is load-bearing for the claim that the hyperboloidal bilinear products furnish intrinsic, physically meaningful excitation coefficients.
Authors: We concur that a sensitivity check is necessary to establish that the excitation coefficients are intrinsic rather than regularization-dependent. The original computations for constant initial data employed the subtraction regularization. In the revised manuscript we have added an explicit verification (new Figure 4 and accompanying text in Section 4) recomputing the excitation coefficients for the same initial data and mode numbers (l=2, n=0,1,2) using all three regularization procedures. The resulting coefficients agree to within 5×10^{-4} relative error across schemes, well within the expected numerical accuracy. This check confirms that the hyperboloidal bilinear products yield physically meaningful, scheme-independent excitation coefficients, as claimed. We have updated the discussion to emphasize this robustness. revision: yes
Circularity Check
No circularity: derivation of divergence and regularization is independent of inputs
full rationale
The paper derives the divergence of the bilinear integrand directly from the reflection/CPT transformation applied to smooth QNM solutions on hyperboloidal slices, as stated in the abstract. Regularization procedures are introduced as separate, explicit constructions that produce finite forms while preserving the orthogonality property. Excitation factors are then defined in terms of these regularized products and computed for specific cases (Schwarzschild scalar perturbations with constant initial data). No step equates a claimed result to its own definition, renames a fitted quantity as a prediction, or reduces the central claim to an unverified self-citation. The logic chain remains self-contained against the stated assumptions about the foliation and mode smoothness.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption QNM solutions are smooth and finite on future-directed hyperboloids
Forward citations
Cited by 2 Pith papers
-
Properties of natural polynomials for Schwarzschild and Kerr black holes
Natural polynomials for Schwarzschild and Kerr quasinormal modes are Pollaczek-Jacobi polynomials with complex parameters, with recurrence peaking at the physical overtone index for Schwarzschild.
-
Relativistic frequency shifts in gravitational waves from axion clouds
A unified relativistic framework using bilinear perturbation theory calculates frequency shifts in GWs from axion clouds, handling self-interactions and multiple superradiant modes for the first time.
Reference graph
Works this paper leans on
-
[1]
Similarly, the flux contributions at the boundaries are related as F∂C S [Ψ 1,∂τΨ 2] = 2F∂C E [Ψ 1, Ψ 2] (97) −r2 h λ2 ψ1 ( p(σ)∂σψ2 +γ(σ)∂τψ2 )⏐⏐⏐ ∂C
From symplectic to energy It is easy to prove that, on shell, the two bilinear forms are related to each other via [19, 48] ΠS[Ψ 1,∂τΨ 2] = 2ΠE[ψ1,ψ2] (96) −r2 h λ2 ψ1 ( p(σ)∂σψ2 +γ(σ)∂τψ2 )⏐⏐⏐ ∂C . Similarly, the flux contributions at the boundaries are related as F∂C S [Ψ 1,∂τΨ 2] = 2F∂C E [Ψ 1, Ψ 2] (97) −r2 h λ2 ψ1 ( p(σ)∂σψ2 +γ(σ)∂τψ2 )⏐⏐⏐ ∂C . In ot...
-
[2]
Hamiltonian action With the help of Eqs.(94) and (96), it is easy to prove that the action of the HamiltonianH =i∂τon the sym- plectic product is such that ΠS[HΨ 1, Ψ 2] + ΠS[Ψ 1,HΨ 2] =ir2 h λ2× × ( p(σ)(ψ2∂σψ1−ψ1∂σψ2) +γ(σ)(ψ2∂τψ1−ψ1∂τψ2) )⏐⏐⏐⏐ ∂C . (98) In the same way, the action of the Hamiltonian in the sympletic-fluxes yields F∂C S [HΨ 1, Ψ 2] +F∂C...
-
[3]
(69) as ˜Π[Ψ 1, Ψ 2] = Π[Ψ 1, Ψ 2]−F∂C[Ψ 1, Ψ 2], (103) which can be used either with the symplectic or the en- ergy currents
Extended bilinear product To keep explicit control over the behaviour at the do- main boundaries, as currently giving rise to the non- normal character of the evolution operator, we extend the definition of theΠ-product in Eqn. (69) as ˜Π[Ψ 1, Ψ 2] = Π[Ψ 1, Ψ 2]−F∂C[Ψ 1, Ψ 2], (103) which can be used either with the symplectic or the en- ergy currents. Fi...
-
[4]
As we saw in the previous subsection, the action of the Hamiltonian on the bilinear product is sym- metric up to a totalσ-derivative, cf. Eqs. (98) and (100) Π[HΨ 1, Ψ 2] + Π[Ψ 1,HΨ 2] = 0 +total deriv. term. (112) Rewriting Eqn.(112) in terms of the orthogonality product and usingJH =−HJ, we obtain ⟨HΨ 1, Ψ 2⟩−⟨Ψ 1,HΨ 2⟩ = 0 +total deriv. term (113) The ...
-
[5]
(117) where c =−1 for the symplectic current andc = 1 for the energy current
It is easy to show that ⟨Ψ 1(x),J Ψ 2(x)⟩τ=c⟨Ψ 2(x),J Ψ 1(x)⟩τ, ⟨Ψ 1(x), Ψ 2(x)⟩τ=c⟨JΨ 2(x),J Ψ 1(x)⟩τ. (117) where c =−1 for the symplectic current andc = 1 for the energy current. TakingΨ 1, Ψ 2 as QNMs, 14 and remembering that the action ofJ on a QNM gives an anti-QNM, the above expressions imply that the orthogonality product between two anti- QNMs is...
-
[6]
Boundary behaviour We expand the expressions(123)-(127) around future null infinity σ0 = ϵand the black-hole horizon σ1 = 1−ϵ, forϵ≳ 0. For that purpose, we make use of the regularity condition for the radial fields provided by the radial equation (45) ∂σϕI1(σ0) = ( −2irhωI1 +iℓ1(ℓ1 + 1) 2rhωI1 ) ×ϕI1(σ0)(128) +O(ϵ) ∂σϕI1(σ1) = ( 1−2irhωI1 + ℓ1(ℓ1 + 1)−(2...
-
[7]
The QNM fieldϕI(σ) solves the radial equation(45)
Semi-analytic integration One route to integrating orthogonality products is through a semi-analytic approach inspired by [21]. The QNM fieldϕI(σ) solves the radial equation(45). This is a confluent Heun equation with regular singular points at σ={1,∞}and an irregular singular point atσ= 0 —see [51]. Then, we use Frobenius method and expand ψI(σ) around t...
-
[8]
up” and “in
Complex contour An alternative approach to computing the orthogonality product relies on choosing a complex integration contour, C, for numerical integration. The complex contour is more easily defined in the un-compactified coordinater, rather than the compactifiedσ, as the singularity atr =∞ (σ= 0) can be more easily approached from a direction of conve...
-
[9]
Initial data regularisation schemes The comparison between the Green’s function approach against the orthogonality projection led to boundary terms, vanishing under the regularisation scheme. Implementing thecomplex contourstrategy is, therefore, straightforward once an explicit form of the initial data set {ϕID(σ), ˙ϕID(σ)}is available to be evaluated at...
-
[10]
we consider the Frobenius expansion(138) and select the solution to the indicial equation withrI = 0 [45]
Analytic continuation approaches For both regularisation schemes, we compute QNM eigenfunctionsψI using Leaver’s method, i.e. we consider the Frobenius expansion(138) and select the solution to the indicial equation withrI = 0 [45]. In particular, we iterate the 3-term recursion relation up to a truncation valuepmax. In the semi-analytic regularisation sc...
-
[11]
Nevertheless, extending the calculation to pmax = 700 does not produce any noticeable change in the main results within the significant digits displayed here
Including additional terms in the sum significantly increases the computational cost. Nevertheless, extending the calculation to pmax = 700 does not produce any noticeable change in the main results within the significant digits displayed here. Whenevaluatingtheproductswiththecomplexcontour method, we include uppmax = 700 terms in the sum. The integration...
-
[12]
large/large
Extended orthogonality product We integrate the extended product(118) numerically for the QNM modes I1 = (ℓ1,m 1,n 1) = (0 , 0, 0) and I2 = (ℓ2,m 2,n 2) = (0 , 0, 1). Here, the eigenfunctions result from the eigenvalue problem(46), with the operator L discretized via a Chebyshev spectral method using a grid based on Chebyshev–Lobatto collocation points [3...
2024
-
[13]
Black hole spectroscopy: from theory to experiment
J. Abediet al., Black hole spectroscopy: from theory to experiment (2025), arXiv:2505.23895 [gr-qc]
work page internal anchor Pith review arXiv 2025
-
[14]
A.G.Abacandetal(LIGOScientific, Virgo, andKAGRA Collaborations), Phys. Rev. Lett.135, 111403 (2025)
2025
- [15]
-
[16]
Spectroscopy of Kerr black holes with Earth- and space-based interferometers
E. Berti, A. Sesana, E. Barausse, V. Cardoso, and K. Belczynski, Phys. Rev. Lett. 117, 101102 (2016), arXiv:1605.09286 [gr-qc]
work page Pith review arXiv 2016
- [17]
- [18]
-
[19]
A. Toubiana, L. Pompili, A. Buonanno, J. R. Gair, and M. L. Katz, Phys. Rev. D 109, 104019 (2024), arXiv:2307.15086 [gr-qc]
-
[20]
E. W. Leaver, Phys. Rev. D34, 384 (1986)
1986
-
[21]
R. H. Price, Phys. Rev. D5, 2419 (1972)
1972
-
[22]
Late-Time Behavior of Stellar Collapse and Explosions: I. Linearized Perturbations
C. Gundlach, R. H. Price, and J. Pullin, Phys. Rev. D 49, 883 (1994), arXiv:gr-qc/9307009
work page Pith review arXiv 1994
-
[23]
M. De Amicis, S. Albanesi, and G. Carullo, Phys. Rev. D 110, 104005 (2024), arXiv:2406.17018 [gr-qc]
-
[24]
Beyond quasinormal modes: a complete mode decomposition of black hole perturbations,
P. Arnaudo, J. Carballo, and B. Withers, (2025), arXiv:2510.18956 [gr-qc]
-
[25]
Price’s law from quasinormal modes,
P. Arnaudo and B. Withers, (2025), arXiv:2511.17703 [gr-qc]
-
[26]
M. De Amicis, E. Cannizzaro, G. Carullo, and L. Sberna, Phys. Rev. D113, 024048 (2026), arXiv:2506.21668 [gr- qc]
- [27]
-
[28]
Prompt Response from Plunging Sources in Schwarzschild Spacetime
S. Ma, (2026), arXiv:2604.08680 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [29]
- [30]
-
[31]
S. R. Green, S. Hollands, L. Sberna, V. Toomani, and P. Zimmerman, Phys. Rev. D107, 064030 (2023)
2023
- [32]
-
[33]
L. T. London, Phys. Rev. D113, 044008 (2026)
2026
- [34]
- [35]
- [36]
-
[37]
A. Hussain and A. Zimmerman, Phys. Rev. D106, 104018 (2022), arXiv:2206.10653 [gr-qc]
- [38]
-
[39]
E. Cannizzaro, L. Sberna, S. R. Green, and S. Hollands, Phys. Rev. Lett.132, 051401 (2024), arXiv:2309.10021 [gr-qc]
-
[40]
E. Hamilton, L. London, J. E. Thompson, E. Fauchon- Jones, M. Hannam, C. Kalaghatgi, S. Khan, F. Pannarale, and A. Vano-Vinuales, Phys. Rev. D104, 124027 (2021), arXiv:2107.08876 [gr-qc]. 25
- [41]
-
[42]
P. Arnaudo, J. Carballo, and B. Withers, JHEP09, 010, arXiv:2505.04696 [hep-th]
-
[43]
A. Zenginoglu, Phys. Rev. D 83, 127502 (2011), arXiv:1102.2451 [gr-qc]
-
[44]
Hyperboloidal approach to quasinormal modes,
R. Panosso Macedo and A. Zenginoglu, Front. in Phys. 12, 1497601 (2024), arXiv:2409.11478 [gr-qc]
-
[45]
R. Panosso Macedo, Phil. Trans. Roy. Soc. Lond. A382, 20230046 (2024), arXiv:2307.15735 [gr-qc]
-
[46]
D. Gajic and C. M. Warnick, Quasinormal modes on Kerr spacetimes (2024), arXiv:2407.04098 [gr-qc]
-
[47]
J. L. Jaramillo, R. P. Macedo, and L. A. Sheikh, Phys. Rev. X11, 031003 (2021)
2021
- [48]
-
[49]
Exceptional Points and Resonance in Black Hole Ringdown,
R. Panosso Macedo, T. Katagiri, K.-i. Kubota, and H. Mo- tohashi, (2025), arXiv:2512.02110 [gr-qc]
-
[50]
Spectral decomposition of black-hole perturbations on hyperboloidal slices
M. Ansorg and R. Panosso Macedo, Phys. Rev. D93, 124016 (2016), arXiv:1604.02261 [gr-qc]
work page Pith review arXiv 2016
-
[51]
J. Besson and J. L. Jaramillo, Gen. Rel. Grav.57, 110 (2025), arXiv:2412.02793 [gr-qc]
- [52]
-
[53]
Penrose, General Relativity and Gravitation43, 901 (2011)
R. Penrose, General Relativity and Gravitation43, 901 (2011)
2011
-
[54]
J. A. Valiente Kroon, Conformal Methods in General Relativity (Cambridge University Press, 2016)
2016
-
[55]
E. Gasperin, R. Panosso Macedo, and J. Feng arXiv:2602.19245 (2026)
-
[56]
Gasperín and J
E. Gasperín and J. L. Jaramillo, Classical and Quantum Gravity39, 115010 (2022)
2022
-
[57]
E. W. Leaver, Proc. Roy. Soc. Lond. A402, 285 (1985)
1985
-
[58]
S. R. Dolan, Phys. Rev. D 76, 084001 (2007), arXiv:0705.2880 [gr-qc]
work page Pith review arXiv 2007
-
[59]
Zenginoglu, Hyperboloidal foliations and scri-fixing.Class
A. Zenginoglu, Class. Quant. Grav.25, 145002 (2008), arXiv:0712.4333 [gr-qc]
-
[60]
S. Hollands and R. M. Wald, Commun. Math. Phys.321, 629 (2013), arXiv:1201.0463 [gr-qc]
-
[61]
Bleistein and R
N. Bleistein and R. A. Handelsman,Asymptotic Expan- sions of Integrals(Dover Publications, New York, 1986)
1986
-
[62]
J. W. Brown and R. V. Churchill,Complex Variables and Applications, 9th ed. (McGraw-Hill Education, New York, 2014)
2014
-
[63]
M. Minucci and R. Panosso Macedo, Gen. Rel. Grav.57, 33 (2025), arXiv:2411.19740 [gr-qc]
- [64]
-
[65]
R. Panosso Macedo, B. Leather, N. Warburton, B. Wardell, and A. Zenginoğlu, Phys. Rev. D105, 104033 (2022), arXiv:2202.01794 [gr-qc]
-
[66]
Y. Zhou and R. Panosso Macedo, Phys. Rev. D112, 084063 (2025), arXiv:2507.05370 [gr-qc]. Appendix A: Conformal transformation of the current Ja We begin by identifying two manifolds: the physical spacetime{M,gab}and the conformal (unphysical) space- time{˜M, ˜gab}, related to each other by a conformal factor Ω via gab = Ω−2˜gab =⇒g = Ω−8˜g. (A1) The asymp...
-
[67]
+” and “−
The “+” and “−” time foliations We consider a generic transformation into coordinates ¯xµ ±= ( τ±,σ±) defined by the height function tech- nique [47] t =τ±∓H(σ±). (B1) It is straightforward to derive the direct transformation ¯xµ +↔¯xµ −as τ+ =τ−+ 2H(σ−), σ+ =σ−. (B2) Note that if we were to identify the height function with the tortoise coordinateH(σ) = ...
-
[68]
The coordinate time-reversal transformations Another transformation that achieves a change in the metric similar to Eq.(B3) is a direct coordinate time- 27 FIG. 5. Carter-Penrose diagrams representing the hyperboloidal foliations associated to the coordinates¯xµ ±, which are related to the Schwarzschild-Droste time viat→τ±∓H(σ). The sign change with respe...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.