Recognition: unknown
Radiation outer boundary conditions and near-to-far field signal transformations for the Bardeen-Press equation
Pith reviewed 2026-05-08 10:29 UTC · model grok-4.3
The pith
Exact radiation outer boundary conditions make the Bardeen-Press equation transparent at any finite radius.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We develop and implement exact radiation outer boundary conditions for the Bardeen-Press equation, making the artificial boundary transparent at any finite radius. We also construct near-to-far field teleportation kernels that map field data recorded at finite radius r1 to the data reaching r2 > r1. The possible choice r2 = infinity corresponds to asymptotic waveform evaluation. Both boundary and teleportation kernels are well approximated by exponential sums, with associated error bounds. Implemented in a time-domain solver, our kernel-based boundary conditions eliminate unphysical late-time growth and give the correct late-time decay rates.
What carries the argument
Exact radiation outer boundary conditions and near-to-far field teleportation kernels for the Bardeen-Press equation, realized through exponential-sum approximations with error bounds.
If this is right
- Long-duration simulations become feasible without corruption by boundary reflections or slowly growing modes.
- Correct late-time decay rates are recovered in the numerical solutions.
- Asymptotic waveforms at future null infinity can be obtained directly from finite-radius data.
- Efficient waveform modeling for extreme mass-ratio inspirals and other black-hole perturbation problems is enabled.
Where Pith is reading between the lines
- The exponential-sum technique for kernels may reduce computational domain size in other linear wave problems with similar outgoing behavior.
- The same mapping approach could be tested for scalar-field or electromagnetic perturbations around Schwarzschild black holes.
- Stable long runs open the possibility of cleanly extracting subtle features such as power-law tails without domain-size artifacts.
Load-bearing premise
The Bardeen-Press equation is assumed to capture the essential dynamics so that the derived boundary conditions and kernels remain exact and stable inside a full time-domain solver over long durations.
What would settle it
A long-duration time-domain evolution that exhibits growing unphysical modes or incorrect late-time power-law decay rates after the proposed boundary conditions and kernels are applied would falsify the exactness and stability claims.
Figures
read the original abstract
Several theoretical and astrophysical problems - including gravitational-wave modeling for extreme mass-ratio inspirals - require accurate time-domain solutions of the spin-weight $s=-2$ Teukolsky equation in Boyer-Lindquist coordinates. Because such simulations are performed on finite computational domains, they typically introduce an artificial outer boundary where nontrivial boundary conditions must be imposed. If these conditions are inaccurate, then spurious reflections and slowly-growing unphysical modes may corrupt long-time evolutions. We develop and implement exact radiation outer boundary conditions for the Bardeen-Press equation (a harmonic moment of the $a=0$ Teukolsky equation), making the artificial boundary transparent at any finite radius. We also construct near-to-far field teleportation kernels that map field data recorded at finite radius $r_1$ to the data reaching $r_2 > r_1$. The possible choice $r_2 = \infty$ corresponds to asymptotic waveform evaluation, that is propagation of the data to future null infinity. We show that both boundary and teleportation kernels are well approximated by exponential sums, with associated error bounds. Implemented in a time-domain solver, our kernel-based boundary conditions eliminate unphysical late-time growth and give the correct late-time decay rates, affording efficient long-duration simulations for waveform modeling and related blackhole perturbation calculations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops and implements exact radiation outer boundary conditions for the Bardeen-Press equation (a harmonic moment of the a=0 Teukolsky equation), rendering artificial outer boundaries transparent at any finite radius in time-domain simulations. It also constructs near-to-far field teleportation kernels that map field data from finite radius r1 to r2 > r1 (including r2 = infinity for asymptotic waveforms at future null infinity). Both the boundary conditions and kernels are shown to be well approximated by exponential sums with associated error bounds. Numerical implementation in a time-domain solver is reported to eliminate unphysical late-time growth and recover the correct late-time decay rates.
Significance. If the central claims hold, this work would be significant for gravitational-wave modeling, especially extreme mass-ratio inspirals, by enabling accurate, efficient long-duration time-domain evolutions of the Teukolsky equation without spurious reflections or growing modes. The exact (within the Bardeen-Press framework) boundary operators and the teleportation kernels, together with the exponential-sum approximations and error bounds, represent a practical advance in boundary handling and far-field extraction for black-hole perturbation calculations.
major comments (2)
- [Numerical results section] Numerical results section: The reported tests confirm correct late-time decay rates for the durations simulated, but the manuscript provides no explicit demonstration or additional long-time runs showing that truncation errors in the exponential-sum approximation of the kernels, when integrated against the solution over arbitrarily long evolutions, remain below the threshold that would produce measurable incoming radiation or alter the known power-law tail (a sensitive diagnostic of exact transparency).
- [Kernel derivation and approximation] Kernel derivation and approximation: The exact integral kernels are derived from the Bardeen-Press equation, but the passage to the exponential-sum surrogate requires a concrete estimate or test of how the stated error bounds behave under repeated convolution over long times; without this, the claim that the approximations 'remain transparent and stable' rests on unverified extrapolation.
minor comments (2)
- [Abstract] Abstract: The statement that the kernels 'are well approximated by exponential sums, with associated error bounds' would be clearer if it referenced the specific section or equation where the bounds are derived and stated.
- [Notation] Notation: The distinction between the exact integral kernels and their exponential-sum surrogates could be made more explicit in the text and figures to avoid potential confusion for readers implementing the method.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment below and have revised the manuscript to incorporate additional numerical evidence and error analysis as suggested.
read point-by-point responses
-
Referee: [Numerical results section] Numerical results section: The reported tests confirm correct late-time decay rates for the durations simulated, but the manuscript provides no explicit demonstration or additional long-time runs showing that truncation errors in the exponential-sum approximation of the kernels, when integrated against the solution over arbitrarily long evolutions, remain below the threshold that would produce measurable incoming radiation or alter the known power-law tail (a sensitive diagnostic of exact transparency).
Authors: We agree that explicit long-time validation is important for confirming that approximation errors do not introduce artifacts. In the revised manuscript we have added results from extended simulations over significantly longer durations. These confirm that the expected power-law tails are preserved without measurable deviations or incoming radiation attributable to the exponential-sum truncation, consistent with the provided error bounds. revision: yes
-
Referee: [Kernel derivation and approximation] Kernel derivation and approximation: The exact integral kernels are derived from the Bardeen-Press equation, but the passage to the exponential-sum surrogate requires a concrete estimate or test of how the stated error bounds behave under repeated convolution over long times; without this, the claim that the approximations 'remain transparent and stable' rests on unverified extrapolation.
Authors: The error bounds on the exponential-sum approximations follow from standard results in approximation theory and are stated explicitly in the manuscript. To address repeated application, the revised version includes a new subsection deriving a bound on the accumulated error under successive convolutions, showing that it remains controlled (growing at most linearly in time but staying below the level that would affect transparency or stability). This is supported by the original bounds together with additional numerical checks over extended evolutions. revision: yes
Circularity Check
Derivation of exact boundary conditions and kernels is self-contained
full rationale
The paper derives radiation outer boundary conditions and near-to-far field teleportation kernels directly from the Bardeen-Press equation by constructing exact integral operators that enforce outgoing-wave behavior at finite radius. These constructions follow from the PDE structure without any parameter fitting, data-driven calibration, or reduction to prior self-citations for the core operators. The subsequent exponential-sum approximation is presented separately with explicit error bounds and is not used to establish the exactness claim. No step in the derivation chain is self-definitional, renames a known result, or imports uniqueness via author-overlapping citations. The central results remain independent of the implementation details and are falsifiable against the analytic properties of the wave equation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The Bardeen-Press equation governs the relevant spin-weight s=-2 perturbations for a=0 black holes.
Reference graph
Works this paper leans on
-
[1]
Black hole spectroscopy: from theory to experiment
J. Abedi et al.,Black hole spectroscopy: from theory to experiment, (2025), https://arxiv. org/abs/2505.23895. 22 Table D.1: ROBC table fors=−2,ℓ= 2,ρ= 30,ε= 10 −13 Beta locations (Re) Beta locations (Im) Gamma residues (Re) Gamma residues (Im) -1.9211957803762693e-01 +0.0e+00 -2.4484092328893475e-09 +0.0e+00 -1.3362488822827323e-01 +0.0e+00 -1.6894339064...
work page internal anchor Pith review arXiv 2025
-
[2]
Abramowitz and I
M. Abramowitz and I. A. Stegun,Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Applied Mathematics Series, vol. 55, 10th printing, December 1972, with corrections
1972
-
[3]
Waveform Modelling for the Laser Interferometer Space Antenna
N. Afshordi et al.,Waveform modelling for the Laser Interferometer Space Antenna, Living Rev. Rel., 28 (2025), 9, https://arxiv.org/abs/2311.01300
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[4]
Alper t, L
B. Alper t, L. Greengard, and T. Hagstrom,Rapid evaluation of nonreflecting boundary kernels for time-domain wave propagation, SIAM Journal on Numerical Analysis, 37 (2000), pp. 1138–1164
2000
-
[5]
Alpert, L
B. Alpert, L. Greengard, and T. Hagstrom,Nonreflecting boundary conditions for the time-dependent wave equation, J. Comput. Phys., 180 (2002), pp. 270–296
2002
-
[6]
Spectral decomposition of black-hole perturbations on hyperboloidal slices
M. Ansorg and R. Panosso Macedo,Spectral decomposition of black-hole perturbations on hyperboloidal slices, Phys. Rev. D, 93 (2016), 124016, https://arxiv.org/abs/1604.02261
work page Pith review arXiv 2016
-
[7]
Appelö and T
D. Appelö and T. Colonius,A high-order super-grid-scale absorbing layer and its application to linear hyperbolic systems, J. Comput. Phys., 228 (2009), pp. 4200–4217
2009
-
[8]
S. Aret akis, G. Khanna, and S. Sabhar w al,An observational signature for extremal black holes, Gen. Rel. Grav., 57 (2025), 160, https://arxiv.org/abs/2307.03963
-
[9]
L. Barack and P. Giudice,Time-domain metric reconstruction for self-force applications, Phys. Rev. D, 95 (2017), 104033, https://arxiv.org/abs/1702.04204
-
[10]
J. M. Bardeen and W. H. Press,Radiation fields in the Schwarzschild background, J. Math. 23 Table D.3: Teleportation table fors=−2,ℓ= 2,ρ 1 = 30,ρ 2 = 60,ε= 10 −13 BetaT locations (Re) BetaT locations (Im) GammaT residues (Re) GammaT residues (Im) -1.9595777859531270e-01 +0.0e+00 -1.2639428338951785e-10 +0.0e+00 -1.3825941343594810e-01 +0.0e+00 -1.0694919...
1973
-
[11]
A. Ba yliss and E. Turkel,Radiation boundary conditions for wave-like equations, Commun. Pure Appl. Math., 33 (1980), pp. 707–725, https://arxiv.org/abs/https://onlinelibrary.wiley. com/doi/pdf/10.1002/cpa.3160330603
- [12]
-
[13]
Quasinormal modes of black holes and black branes
E. Ber ti, V. Cardoso, and A. O. Starinets,Quasinormal modes of black holes and black branes, Class. Quant. Grav., 26 (2009), 163001, https://arxiv.org/abs/0905.2975
work page internal anchor Pith review arXiv 2009
- [14]
-
[15]
Bizzozero, J
D. Bizzozero, J. Ellison, K. Heinemann, and S. Lau,Rapid evaluation of two-dimensional retarded time integrals, J. Comput. Appl. Math., 324 (2017), pp. 118–141
2017
- [16]
-
[17]
L. T. Buchman, M. D. Duez, M. Morales, M. A. Scheel, T. M. Kostersitz, A. M. Ev ans, and K. Mitman,Numerical relativity multimodal waveforms using absorbing boundary conditions, Class. Quant. Grav, 41 (2024), 175011
2024
-
[18]
L. T. Buchman and O. C. A. Sarbach,Towards absorbing outer boundaries in general 24 relativity, Class. Quant. Grav, 23 (2006), pp. 6709–6744
2006
-
[19]
L. T. Buchman and O. C. A. Sarbach,Improved outer boundary conditions for einstein’s field equations, Class. Quant. Grav, 24 (2007), pp. S307–S326
2007
- [20]
-
[21]
V. Cardoso, K. Destounis, F. Duque, R. P. Macedo, and A. Maselli,Black holes in galaxies: Environmental impact on gravitational-wave generation and propagation, Phys. Rev. D, 105 (2022), L061501, https://arxiv.org/abs/2109.00005
-
[22]
Chandrasekhar,On the equations governing the perturbations of the Schwarzschild black hole, Proc
S. Chandrasekhar,On the equations governing the perturbations of the Schwarzschild black hole, Proc. Roy. Soc. Lond. A, 343 (1975), pp. 289–298, https://arxiv.org/abs/https: //royalsocietypublishing.org/rspa/article-pdf/343/1634/289/61331/rspa.1975.0066.pdf. [24]S. Chandrasekhar,The Mathematical Theory of Black Holes, Clarendon Press, 1998
-
[23]
Csukás, I
K. Csukás, I. Rácz, and G. Z. Tóth,Numerical investigation of the dynamics of linear spin s fields on a Kerr background: Late-time tails of spins = ±1,± 2fields, Phys. Rev. D, 100 (2019), 104025
2019
- [24]
-
[25]
Donninger, W
R. Donninger, W. Schlag, and A. Soffer,A proof of Price’s law on Schwarzschild black hole manifolds for all angular momenta, Adv. Math., 226 (2011), pp. 484–540
2011
-
[26]
Erdélyi,Asymptotic Expansions, Dover Books on Mathematics, Dover Publications, New York, 1956
A. Erdélyi,Asymptotic Expansions, Dover Books on Mathematics, Dover Publications, New York, 1956. Unabridged and unaltered republication of Technical report 3, prepared under contract Nonr-220(11) for the Office of Naval Research, 1955
1956
- [27]
- [28]
-
[29]
Z. Gelles and F. Pretorius,Accumulation of charge on an extremal black hole’s event horizon, Phys. Rev. D, 112 (2025), 064003, https://arxiv.org/abs/2503.04881
-
[30]
Greengard, T
L. Greengard, T. Hagstrom, and S. Jiang,The solution of the scalar wave equation in the exterior of a sphere, J. Comput. Phys., 274 (2014), pp. 191–207
2014
-
[31]
Hagstrom and S
T. Hagstrom and S. Hariharan,A formulation of asymptotic and exact boundary conditions using local operators, Appl. Numer. Math., 27 (1998), pp. 403–416. Special Issue on Absorbing Boundary Conditions
1998
-
[32]
S. A. Hughes,Computing radiation from Kerr black holes: Generalization of the Sasaki- Nakamura equation, Phys. Rev. D, 62 (2000), 044029
2000
-
[33]
S. A. Hughes, A. Apte, G. Khanna, and H. Lim,Learning about black hole binaries from their ringdown spectra, Phys. Rev. Lett., 123 (2019), 161101, https://arxiv.org/abs/1901. 05900
2019
- [34]
-
[35]
S. R. Lau,Rapid evaluation of radiation boundary kernels for time-domain wave propagation on black holes: Implementation and numerical tests, Class. Quant. Grav., 21 (2004), pp. 4147–4192
2004
- [36]
- [37]
-
[38]
E. W. Lea ver,An Analytic representation for the quasi-normal modes of Kerr black holes, Proc. Roy. Soc. Lond. A, 402 (1985), pp. 285–298
1985
-
[39]
X. Li, A. Taflove, and V. Backman,Modified fdtd near-to-far-field transformation for improved backscattering calculation of strongly forward-scattering objects, IEEE Antennas Wirel. Propag. Lett., 4 (2005), pp. 35–38
2005
- [40]
-
[41]
O. Long and L. Barack,Time-domain metric reconstruction for hyperbolic scattering, Phys. Rev. D, 104 (2021), 024014, https://arxiv.org/abs/2105.05630
-
[42]
S. Ma et al.,Fully relativistic three-dimensional Cauchy-characteristic matching for physical degrees of freedom, Phys. Rev. D, 109 (2024), 124027, https://arxiv.org/abs/2308.10361. 25
- [43]
-
[44]
Muller, G
J. Muller, G. Parent, G. Jeandel, and D. Lacroix,Finite-difference time-domain and near-field-to-far-field transformation in the spectral domain, J. Opt. Soc. Am. A, 28 (2011), pp. 868–878
2011
-
[45]
Panosso Macedo,Hyperboloidal framework for the Kerr spacetime, Class
R. Panosso Macedo,Hyperboloidal framework for the Kerr spacetime, Class. Quant. Grav., 37 (2020), 065019, https://arxiv.org/abs/1910.13452
-
[46]
R. H. Price,Nonspherical Perturbations of Relativistic Gravitational Collapse. II. Integer-Spin, Zero-Rest-Mass Fields, Phys. Rev. D, 5 (1972), pp. 2439–2454
1972
-
[47]
K. Rink, R. Bachhar, T. Islam, N. E. M. Rif at, K. Gonzalez-Quesada, S. E. Field, G. Khanna, S. A. Hughes, and V. V arma,Gravitational wave surrogate model for spinning, intermediate mass ratio binaries based on perturbation theory and numerical relativity, Phys. Rev. D, 110 (2024), 124069, https://arxiv.org/abs/2407.18319
- [48]
-
[49]
S. Y. Sla vyanov and W. Lay,Special Functions: A Unified Theory Based on Singularities, Oxford Mathematical Monographs, Oxford University Press, Oxford, 2000. [52]L. C. Stein,Probes of strong-field gravity, PhD thesis, MIT, 2012
2000
- [50]
-
[51]
Taflove and S
A. Taflove and S. C. Hagness,Computational Electrodynamics: The Finite-Difference Time-Domain Method, Artech House, Norwood, MA, 3rd ed., 2005
2005
-
[52]
S. A. Teukolsky,Perturbations of a Rotating Black Hole. I. Fundamental Equations for Gravitational, Electromagnetic, and Neutrino-Field Perturbations, Astrophys. J., 185 (1973), pp. 635–648. [56]S. A. Teukolsky and W. H. Press,Perturbations of a rotating black hole. III. Interaction of the hole with gravitational and electromagnetic radiation, Astrophys. ...
1973
-
[53]
K. S. Thorne,Multipole Expansions of Gravitational Radiation, Rev. Mod. Phys., 52 (1980), pp. 299–339
1980
-
[54]
Tokita,Exponential decay of solutions for the wave equation in the exterior domain with spherical boundary, Math
T. Tokita,Exponential decay of solutions for the wave equation in the exterior domain with spherical boundary, Math. Kyoto Univ, 12 (1972), pp. 413–430
1972
-
[55]
M. Vishal, S. E. Field, K. Rink, S. Gottlieb, and G. Khanna,Toward exponentially- convergent simulations of extreme-mass-ratio inspirals: A time-domain solver for the scalar Teukolsky equation with singular source terms, Phys. Rev. D, 110 (2024), 104009, https://arxiv.org/abs/2307.01349
-
[56]
B. W ardell, A. Pound, N. W arburton, J. Miller, L. Durkan, and A. Le Tiec, Gravitational Waveforms for Compact Binaries from Second-Order Self-Force Theory, Phys. Rev. Lett., 130 (2023), 241402, https://arxiv.org/abs/2112.12265
-
[57]
W asow,Asymptotic Expansions for Ordinary Differential Equations, Dover Books on Mathematics, Dover Publications, Inc., New York, 1987
W. W asow,Asymptotic Expansions for Ordinary Differential Equations, Dover Books on Mathematics, Dover Publications, Inc., New York, 1987. Reprint of the 1965 edition
1987
-
[58]
C. H. Wilcox,The initial-boundary value problem for the wave equation in an exterior domain with spherical boundary, Notices of the AMS, 6 (1959), pp. 869–870. Abstract No. 564–20
1959
-
[59]
Xu and S
K. Xu and S. Jiang,A bootstrap method for sum-of-poles approximations, J. Sci. Comput., 55 (2013), pp. 16–39
2013
-
[60]
Zenginoglu,Hyperboloidal layers for hyperbolic equations on unbounded domains, J
A. Zenginoglu,Hyperboloidal layers for hyperbolic equations on unbounded domains, J. Comput. Phys., 230 (2011), pp. 2286–2302, https://arxiv.org/abs/1008.3809
-
[61]
A. Zenginoglu and G. Khanna,Null infinity waveforms from extreme-mass-ratio inspirals in Kerr spacetime, Phys. Rev. X, 1 (2011), 021017, https://arxiv.org/abs/1108.1816
-
[62]
A. Zenginoglu, D. Nunez, and S. Husa,Gravitational perturbations of Schwarzschild spacetime at null infinity and the hyperboloidal initial value problem, Class. Quant. Grav., 26 (2009), 035009, https://arxiv.org/abs/0810.1929. 26
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.