pith. sign in

arxiv: 2605.26028 · v1 · pith:ISIAG7MOnew · submitted 2026-05-25 · 🌀 gr-qc · cs.NA· math.NA

Hyperboloidal evolution for scalar scattering in Minkowski space

Ekrem S Demirbo\u{g}a , An{\i}l Zengino\u{g}lu This is my paper

Pith reviewed 2026-06-29 20:45 UTC · model grok-4.3

classification 🌀 gr-qc cs.NAmath.NA
keywords hyperboloidal evolutionconformal compactificationscalar wave scatteringMinkowski spacetimenumerical relativitynull infinityPenrose domainglobal evolution
0
0 comments X

The pith

Exact conformal matching of hyperboloidal and Penrose domains enables global scalar wave evolution from past to future null infinity without artificial boundaries.

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

The paper presents a time-domain numerical method that glues three compactified regions together to evolve scalar waves across all of Minkowski spacetime. A past hyperboloidal domain reaches back to past null infinity, a central Penrose patch covers the neighborhood of spatial infinity, and a future hyperboloidal domain reaches forward to future null infinity. The regions join along identical conformal hypersurfaces so the evolution crosses the interfaces directly. Linear tests achieve stable fourth-order convergence with radiation extracted at future null infinity. Nonlinear tests mostly retain this accuracy except when the rescaled source fails to vanish at the boundaries.

Core claim

The central claim is that an exact conformal matching of a past hyperboloidal domain attached to past null infinity, a Penrose domain covering a neighborhood of spatial infinity, and a future hyperboloidal domain attached to future null infinity yields a global evolution scheme for scalar waves without artificial timelike outer boundaries and without interpolation between scri-fixing gauges.

What carries the argument

The exact conformal matching of the three compactified regions along identical conformal hypersurfaces.

If this is right

  • Stable propagation occurs across the matching interfaces for free waves and linear scattering potentials.
  • Fourth-order convergence holds in free and linear-potential tests with direct extraction of radiation at future null infinity.
  • Quintic and septic nonlinearities recover approximately fourth-order convergence and the expected late-time tail rates.
  • The cubic nonlinearity yields only first-order convergence because the rescaled source remains non-vanishing near the boundaries.

Where Pith is reading between the lines

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

  • The same matching surfaces could support longer evolutions if the boundary regularity issue for nonlinear sources is resolved.
  • Avoiding artificial outer boundaries removes a common source of spurious reflections in long-time wave simulations.
  • The construction supplies a template for linking different foliations across spatial infinity in other asymptotically flat settings.
  • Improving the treatment of non-vanishing sources at compactified boundaries would extend the method to a wider class of nonlinear problems.

Load-bearing premise

The conformally rescaled nonlinear source term must vanish near the compactified boundaries or be treatable without loss of regularity for the matching to stay stable and high-order accurate.

What would settle it

A simulation of the cubic nonlinearity case that produces only first-order convergence instead of fourth-order, or any instability detected when crossing the matching surfaces.

Figures

Figures reproduced from arXiv: 2605.26028 by An{\i}l Zengino\u{g}lu, Ekrem S Demirbo\u{g}a.

Figure 1
Figure 1. Figure 1: Coordinate systems used in the numerical construction: (a) Penrose coordinates, (b) past hyperboloidal coordinates, (c) future hyperboloidal coordinates. The hypersurfaces T = ±π/2 in the Penrose diagram coincide with the hyperboloidal slices τ± = 0 in the respective hyperboloidal diagrams (emphasized as thick, horizontal black lines). Instead, we switch to time-translation invariant (stationary) hyperbolo… view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of strategies for bridging past null infinity I − and future null infinity I + in the scattering problem. Left: The method used in [22], which relies exclusively on hyperboloidal foliations. This approach yields overlapping slices and leaves an uncovered causal gap near spatial infinity i0 (red line). Right: The hybrid approach proposed in this work. By combining hyperboloidal foliations near th… view at source ↗
Figure 3
Figure 3. Figure 3: Snapshots of the potentials used in this study in different coordinate systems. Left: P¨oschl–Teller. Right: shaking barrier [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Global evolution of a free scalar wave superimposed on Penrose diagrams for two distinct incident profiles: a standard Gaussian pulse (left) and a modulated Gaussian packet G(s) = A sin(ω(s − τ0)) exp(−(s − τ0) 2/(2σ 2 )) (right). The plotted asymptotic profiles are the oriented radiation fields R− on I − and R+ on I +, defined in (64). With this normalization the free scattering map is the identity. 2 1 0… view at source ↗
Figure 5
Figure 5. Figure 5: Convergence analysis against the exact analytical solution. Left Panel: Time evolution of the L2 error between the numerical and exact solutions of the free conformal wave equation. Right Panel: The estimated convergence order, showing the expected fourth-order behavior across the coordinate patches. 5.2. Wave equation with potentials We consider the propagation of scalar waves in the presence of localized… view at source ↗
Figure 6
Figure 6. Figure 6: Global scalar field evolution across Penrose diagrams for the four localized potential models described above. The incident Gaussian pulse at past null infinity (I −) and the corresponding extracted radiation at future null infinity (I +) are plotted along the null boundaries. The resulting deformations in the asymptotic outgoing profiles show the scattering response of each potential. Because no exact sol… view at source ↗
Figure 7
Figure 7. Figure 7: Left Panel: Self-convergence factor Q for the four localized potentials. Right Panel: Pointwise difference between successive resolutions after restriction to a common grid. Norms in the Penrose patch are computed over the active physical interval. numerically [45, 46]. The P¨oschl–Teller spectrum is ωn = ± r V0 − α2 4 − iα  n + 1 2  , n = 0, 1, 2, . . . . (74) In our spherically symmetric 3D framework, … view at source ↗
Figure 8
Figure 8. Figure 8: Characteristic quasinormal ringdown of the P¨oschl–Teller potential for two different barrier configurations. Left panel: Configuration with V0 = 10 and α = 3. The numerical fit gives ω1 = 2.7837 − 4.5004i, differing from the exact analytical value ω1 = 2.7839 − 4.5000i by 0.005% in the real part and 0.008% in the imaginary part. Right panel: Configuration with V0 = 10 and α = 1. The numerical fit gives ω1… view at source ↗
Figure 9
Figure 9. Figure 9: Global scalar field solutions on Penrose diagrams for the cubic (p = 3), quintic (p = 5), and septic (p = 7) semilinear wave equations. The proper resolution of these nonlinearities places stringent demands on the numerical scheme near the conformal boundaries. As shown in Section 3, the conformally rescaled source term scales as Ωp−3 [PITH_FULL_IMAGE:figures/full_fig_p022_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Convergence factor Q for the cubic (p = 3), quintic (p = 5), and septic (p = 7) semilinear wave equations. The quintic and septic cases show approximately fourth-order behavior. The cubic case shows only about first-order convergence, reflecting a limitation of the present finite-difference treatment when the conformally rescaled nonlinear source remains nonzero at the conformal boundary [PITH_FULL_IMAGE… view at source ↗
Figure 11
Figure 11. Figure 11: Late-time radiation tails extracted at I + for the cubic (p = 3), quintic (p = 5), and septic (p = 7) semilinear wave equations. The agreement with the expected tail rates at late retarded times validates the extraction of nonlinear radiation tails at I +. 6. Discussion We have presented a global numerical framework for wave scattering from past null infinity to future null infinity in Minkowski spacetime… view at source ↗
read the original abstract

We develop a time-domain numerical framework for global scalar wave scattering in Minkowski spacetime. The main contribution is an exact conformal matching of three compactified regions: a past hyperboloidal domain attached to $\mathscr I^-$, a Penrose domain covering a neighborhood of spatial infinity $i^0$, and a future hyperboloidal domain attached to $\mathscr I^+$. The matching surfaces are identical conformal hypersurfaces in the adjacent charts. This yields a global evolution scheme connecting $\mathscr I^-$, the neighborhood of $i^0$, and $\mathscr I^+$ without artificial timelike outer boundaries and without interpolation between scri-fixing gauges. We implement the construction for spherically symmetric scalar waves, including free propagation, localized linear scattering potentials such as the P\"oschl--Teller potential, and semilinear wave equations with cubic, quintic, and septic nonlinearities. The numerical experiments demonstrate stable propagation across the matching interfaces, direct extraction of radiation at $\mathscr I^+$, and fourth-order convergence for the free and linear-potential tests. The quintic and septic nonlinear tests exhibit approximately fourth-order convergence and recover the expected late-time tail rates. The cubic case, by contrast, shows only first-order convergence, revealing a limitation of our treatment near compactified boundaries when the conformally rescaled nonlinear source remains non-vanishing. These results validate the conformal matching strategy for long-time simulations, while identifying the boundary regularity issues that must be addressed using a more robust treatment of spatial infinity.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 3 minor

Summary. The paper claims to develop a global time-domain numerical framework for scalar wave scattering in Minkowski spacetime via an exact conformal matching of three compactified regions: a past hyperboloidal domain attached to I^-, a Penrose domain around i^0, and a future hyperboloidal domain attached to I^+. This matching on identical conformal hypersurfaces enables stable evolution connecting I^-, the neighborhood of i^0, and I^+ without artificial timelike boundaries or gauge interpolation. Tests for spherically symmetric scalars (free, linear Pöschl-Teller potential, and semilinear cubic/quintic/septic nonlinearities) demonstrate stable interface propagation, direct radiation extraction at I^+, fourth-order convergence in free/linear/quintic/septic cases, and recovery of expected late-time tails, while the cubic case yields only first-order convergence due to non-vanishing conformally rescaled source.

Significance. If the central matching construction holds, the work provides a technically novel route to boundary-free global evolutions in asymptotically flat settings, with direct access to radiation at I^+ and no need for artificial outer boundaries. The recovery of expected tail rates in the quintic and septic cases, together with fourth-order convergence where the source-regularity assumption is satisfied, constitutes a concrete strength. The explicit identification of the cubic limitation is also a positive feature, as it clarifies the scope of the current implementation.

major comments (1)
  1. [Numerical experiments / abstract] Numerical experiments (cubic nonlinearity test): only first-order convergence is reported, which directly follows from the non-vanishing of the conformally rescaled nonlinear source near the compactified boundaries. This is load-bearing for the claim of a high-order global scheme applicable to general semilinear equations, as the abstract itself flags the assumption failure.
minor comments (3)
  1. [Abstract] Abstract: verify spelling of 'P"oschl--Teller' (standard form is Pöschl-Teller) and ensure consistent math-mode rendering of symbols such as ℐ^-, ℐ^+, i^0.
  2. [Matching construction] The description of the three-domain matching would benefit from an explicit statement (even if brief) of the continuity conditions imposed on the rescaled field and its first derivatives at the interfaces.
  3. [Numerical results] Minor: the abstract states 'approximately fourth-order convergence' for quintic/septic; a table or plot caption should quantify the observed rates more precisely.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment and for highlighting this point. We address the major comment below.

read point-by-point responses

  1. Referee: Numerical experiments (cubic nonlinearity test): only first-order convergence is reported, which directly follows from the non-vanishing of the conformally rescaled nonlinear source near the compactified boundaries. This is load-bearing for the claim of a high-order global scheme applicable to general semilinear equations, as the abstract itself flags the assumption failure.

    Authors: We agree that the first-order convergence observed in the cubic case follows directly from the non-vanishing of the conformally rescaled nonlinear source near the compactified boundaries. The manuscript already states this explicitly in both the abstract and the numerical experiments section, and flags the assumption failure for the cubic nonlinearity. The framework is shown to achieve fourth-order convergence precisely when the source-regularity condition is satisfied (free evolution, linear potential, quintic, and septic cases). We do not claim a high-order scheme for arbitrary semilinear equations without this regularity; the cubic result is presented to illustrate the limitation. No revision is required. revision: no

Circularity Check

0 steps flagged

No circularity in conformal matching construction

full rationale

The paper presents a numerical framework whose central step is the explicit geometric construction of matching surfaces as identical conformal hypersurfaces across three compactified domains (past hyperboloidal, Penrose neighborhood of i^0, future hyperboloidal). This matching is defined directly from the coordinate charts and is not derived from or reduced to any fitted parameter, self-referential equation, or prior self-citation. Validation proceeds via independent numerical experiments that measure convergence orders and late-time tails against known analytic expectations; the cubic nonlinearity failure is reported as a direct consequence of the stated regularity assumption on the rescaled source term. No load-bearing step reduces by construction to its own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The construction relies on established properties of conformal compactifications and hyperboloidal foliations in Minkowski space; no new free parameters, invented entities, or ad-hoc axioms are introduced beyond the standard mathematical background.

axioms (1)
  • standard math Minkowski spacetime admits a conformal compactification allowing hyperboloidal slices attached to I^- and I^+ together with a Penrose chart near i^0.
    The matching strategy presupposes the existence and compatibility of these standard compactified coordinate systems.

pith-pipeline@v0.9.1-grok · 5812 in / 1446 out tokens · 33597 ms · 2026-06-29T20:45:03.941816+00:00 · methodology

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

47 extracted references · 11 canonical work pages · 4 internal anchors

  1. [1]

    Penrose R 1963Physical Review Letters1066–68

  2. [2]

    Series A

    Penrose R 1965Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 284159–203

  3. [3]

    Friedlander F G 1980Mathematical Proceedings of the Cambridge Philosophical Society88483–515

  4. [4]

    Friedrich H 1981Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences378 401–421

  5. [5]

    Baez J C, Segal I E and Zhou Z F 1990Journal of Functional Analysis93239–269

  6. [6]

    Mason L J and Nicolas J P 2004Journal of Hyperbolic Differential Equations1197–233

  7. [7]

    Bishop N T and Rezzolla L 2016Living Reviews in Relativity192 (PreprintarXiv:1606.02532)

    work page internal anchor Pith review Pith/arXiv arXiv

  8. [8]

    Bishop N T, G´ omez R, Lehner L, Maharaj M and Winicour J 1996Physical Review D546153–6165 (PreprintarXiv:gr-qc/9705033)

    work page internal anchor Pith review Pith/arXiv arXiv

  9. [9]

    Moxon J, Scheel M A, Teukolsky S A, Deppe N, Vu N, H´ ebert F, Kidder L E and Throwe W 2023Physical Review D107064013 Hyperboloidal evolution for scalar scattering in Minkowski space25

  10. [10]

    Bernuzzi S, Fontbut´ e J, Albanesi S and Zengino˘ glu A 2025Physical Review D112084036

  11. [11]

    Lovelace G 2009Classical and Quantum Gravity26114002

  12. [12]

    Pfeiffer H P and York J W 2003Physical Review D67

  13. [13]

    Comp` ere G, Gralla S E and Wei H 2023 An asymptotic framework for gravitational scattering (Preprint arXiv:2303.17124)

    work page arXiv 2023

  14. [14]

    Valiente Kroon J A and Gomes Da Silva L J 2024General Relativity and Gravitation5685

  15. [15]

    Friedrich H 1983Communications in Mathematical Physics91445–472

  16. [16]

    Frauendiener J 2004Living Reviews in Relativity71

  17. [17]

    Lindblad H and Schlue V 2023Journal of Hyperbolic Differential Equations20155–218

  18. [18]

    Lindblad H and Schlue V 2025Pure and Applied Analysis7865–926

  19. [19]

    Kadar I and Kehrberger L 2025 Scattering, polyhomogeneity and asymptotics for quasilinear wave equations from past to future null infinity (PreprintarXiv:2501.09814)

    work page arXiv 2025

  20. [20]

    Marajh J, Taujanskas G and Kroon J A V 2025 Controlled regularity at future null infinity from past asymptotic initial data: the wave equation (PreprintarXiv:2508.04690)

    work page arXiv 2025

  21. [21]

    Zengino˘ glu A 2026 From Penrose to Melrose: computing scattering amplitudes at infinity for unbounded media (PreprintarXiv:2601.04167)

    work page arXiv 2026

  22. [22]

    Duarte-Baptista P, Va˜ n´ o-Vi˜ nuales A and del R´ ıo A 2026Classical and Quantum Gravity43045009

  23. [23]

    Frauendiener J, Stevens C and Thwala S 2025Physical Review Letters134161401

  24. [24]

    Friedrich H 1998Journal of Geometry and Physics2483–163

  25. [25]

    Zengino˘ glu A 2008Classical and Quantum Gravity25145002

  26. [26]

    Zenginoglu A 2008Class. Quant. Grav.25195025 (PreprintarXiv:0808.0810)

    work page internal anchor Pith review Pith/arXiv arXiv

  27. [27]

    Va˜ n´ o-Vi˜ nuales A, Husa S and Hilditch D 2015Class. Quant. Grav.32175010 (PreprintarXiv:1412.3827)

    work page internal anchor Pith review Pith/arXiv arXiv

  28. [28]

    Gasperin E, Gautam S, Hilditch D and Va˜ n´ o-Vi˜ nuales A 2020Classical and Quantum Gravity37035006

  29. [29]

    Gautam S, Va˜ n´ o-Vi˜ nuales A, Hilditch D and Bose S 2021Physical Review D103084045

  30. [30]

    Peterson C, Gautam S, Va˜ n´ o-Vi˜ nuales A and Hilditch D 2024Physical Review D110124033

  31. [31]

    Va˜ n´ o-Vi˜ nuales A and Valente T 2024Gen. Rel. Grav.56135 (PreprintarXiv:2408.08952)

    work page arXiv

  32. [32]

    ´Alvares J D and Va˜ n´ o-Vi˜ nuales A 2025Phys. Rev. D112104053 [Erratum: Phys.Rev.D 113, 049902 (2026)] (PreprintarXiv:2506.15311)

    work page arXiv 2026

  33. [33]

    Kroon J A V 2023Conformal Methods in General Relativity(Cambridge University Press)

  34. [34]

    Carroll S M 2019Spacetime and geometry(Cambridge University Press)

  35. [35]

    Zengino˘ glu A 2024Am. J. Phys.92965–974 (PreprintarXiv:2404.01528)

    work page arXiv

  36. [36]

    Zengino˘ glu A 2008Classical and Quantum Gravity25175013

  37. [37]

    Bizo´ n P, Rostworowski A and Zengino˘ glu A 2010Classical and Quantum Gravity27175003

  38. [38]

    Zengino˘ glu A 2011Journal of Computational Physics2302286–2302

  39. [39]

    Ripley J L, Loutrel N, Giorgi E and Pretorius F 2021Physical Review D103104018

  40. [40]

    Panosso Macedo R 2024Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences382

  41. [41]

    Bourg P, Macedo R P, Spiers A, Leather B, Bonga B and Pound A 2025Physical Review D112044049

  42. [42]

    Bizo´ n P and Zengino˘ glu A 2009Nonlinearity222473

  43. [43]

    Donninger R and Zengino˘ glu A 2014Analysis & PDE7461–495

  44. [44]

    Rinne O 2025Nonlinearity38105026

  45. [45]

    Ferrari V and Mashhoon B 1984Physical Review D30295

  46. [46]

    Cardona A F and Molina C 2017Classical and Quantum Gravity34245002

  47. [47]

    Kroon J A V 2003Classical and Quantum Gravity20L53 Appendix A. The free solution and its conformal representations In this section we construct an explicit non-trivial spherically symmetric solution of the standard flat wave equation □ηϕ= 0,(A.1) and rewrite it as a solution of the corresponding conformal wave equations in (i) hyperboloidal coordinates (τ...