pith. sign in

arxiv: 2601.09101 · v2 · submitted 2026-01-14 · 🌀 gr-qc · nlin.CD

Robust Wada Boundaries and Entropy Scaling in pp-Wave Spacetimes

Pedro Henrique Barboza Rossetto , Vanessa Carvalho de Andrade , Daniel M\"uller This is my paper

Pith reviewed 2026-05-16 15:00 UTC · model grok-4.3

classification 🌀 gr-qc nlin.CD
keywords pp-wave spacetimesWada propertyescape basinsbasin entropygeodesic dynamicspolynomial profilesdynamical uncertaintyfractal boundaries
0
0 comments X p. Extension

The pith

Wada escape basin boundaries in pp-wave spacetimes remain maximally intermingled as the polynomial degree increases.

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

This paper examines geodesic motion in pp-wave spacetimes whose profiles are polynomial functions, a setup equivalent to a particle moving in a two-dimensional harmonic potential of the same degree. The authors find that the Wada property, where basin boundaries are shared by all escape regions, holds robustly no matter how high the degree becomes. They measure the resulting uncertainty using basin entropy, which grows steadily with the degree, and show that boundary entropy exceeds the natural log of two for degrees above three, proving the boundaries are fractal. These results indicate that predicting the long-term behavior of geodesics becomes systematically harder in more complex polynomial cases.

Core claim

The Wada property of the escape basins is robust under variation of the polynomial degree, meaning the basin boundaries remain maximally intermingled as the number of escape channels increases. The basin entropy Sb and the boundary basin entropy Sbb increase monotonically with the polynomial degree, indicating enhanced unpredictability, and Sbb is greater than ln(2) for n greater than 3, confirming that the basin boundaries are fractal.

What carries the argument

the dynamical equivalence between geodesic motion in polynomial pp-wave metrics and classical particle motion in a two-dimensional harmonic polynomial potential, which permits direct application of basin-boundary analysis to the spacetime geodesics

If this is right

  • Higher polynomial degrees produce more escape channels while preserving full intermingling of their boundaries.
  • Both basin entropy and boundary basin entropy grow steadily, providing a quantitative measure of increasing dynamical uncertainty.
  • For degrees above three the boundary entropy exceeds ln(2), establishing that the boundaries are fractal.
  • The long-term fate of geodesics becomes progressively harder to predict from nearby initial conditions as the profile complexity rises.

Where Pith is reading between the lines

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

  • The same robustness of Wada boundaries may appear in other wave-like spacetimes whenever an analogous reduction to a polynomial potential exists.
  • Numerical checks at still higher degrees could determine whether the entropy growth continues without bound or eventually saturates.
  • The monotonic rise in uncertainty suggests that profile complexity itself acts as a tunable parameter for the degree of chaos in geodesic motion.

Load-bearing premise

The geodesic motion in the pp-wave spacetime is exactly equivalent to the motion of a classical particle in the corresponding two-dimensional harmonic polynomial potential.

What would settle it

A direct numerical computation of the basins for any polynomial degree n greater than 3 in which the boundary entropy Sbb drops to or below ln(2), or in which some point on the boundary fails to be a limit point of all escape basins, would falsify the robustness claim.

Figures

Figures reproduced from arXiv: 2601.09101 by Daniel M\"uller, Pedro Henrique Barboza Rossetto, Vanessa Carvalho de Andrade.

Figure 1
Figure 1. Figure 1: (a) Trajectory of a ring of particles equally spaced around the circle [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: One-dimensional escape basin for n = 5 and for particles distributed along a unit circle with zero initial velocity. 4 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Exit basins for the px-py plane for the polynomial degree values of (a) n = 4 and (b) n = 5. The region outside of the disk leads to unphysical initial conditions. The color gradient indicates the escape channels, with dark blue being the first channel and yellow being the last channel. initial classification can incorrectly classify a Wada point as a non-Wada point due to insufficient resolution of the lo… view at source ↗
Figure 4
Figure 4. Figure 4: Three zoom levels of the exit basin for n = 5. Each successive zoom level gives a magnification of 10 times. The color gradient indicates the escape channels, with dark blue being the first channel and yellow being the fifth channel. 6 [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Grid method applied to the n = 5 basins. (a) For the angular position space, displayed in [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: (a) Basin entropy by box length for 3 ≤ n ≤ 10. Points represent the calculated values, and the lines are linear regression in the log-log plot. The regression coefficients are shown next to the plot. (b) Boundary basin entropy by box length for 3 ≤ n ≤ 10. The broken line marks the value ln(2). the boundary basin entropy is greater than ln(2) for all box sizes. Therefore, this demonstrates that all these … view at source ↗
read the original abstract

We study the dynamics of the geodesics of pp-wave spacetimes with polynomial profiles, which are dynamically equivalent to the motion of a classical particle in a two-dimensional harmonic polynomial potential. We demonstrate that the Wada property of the escape basins is robust under variation of the polynomial degree, i.e., the basin boundaries remain maximally intermingled as the number of escape channels increases. We further provide a quantitative characterization of the degree of dynamical uncertainty by computing the basin entropy $S_{b}$ and the boundary basin entropy $S_{bb}$. We find that these measures increase monotonically with the polynomial degree, indicating enhanced unpredictability of the final state of the system. We also show that $S_{bb}$ is greater than $\ln(2)$ for $n>3$, and this confirms that the basin boundaries are fractal.

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

2 major / 2 minor

Summary. The paper studies geodesic motion in pp-wave spacetimes with polynomial profiles, dynamically equivalent to a classical particle in a 2D harmonic polynomial potential. It demonstrates the robustness of the Wada property of escape basins under increasing polynomial degree, with boundaries remaining maximally intermingled. Basin entropy S_b and boundary basin entropy S_bb increase monotonically with degree, with S_bb > ln(2) for n > 3 confirming fractal boundaries.

Significance. This result, if numerically sound, highlights persistent chaotic features in general relativistic systems modeled by pp-waves, providing a bridge to classical chaos and quantifying unpredictability via entropy measures that scale with system complexity.

major comments (2)
  1. [Numerical basin analysis] The robustness claim relies on basin plots from trajectory integrations, but no details are given on the initial-condition grid resolution or tests for convergence as n increases. Higher degrees produce finer intermingling, so fixed grids risk misclassifying boundaries and artifactually supporting the Wada invariance.
  2. [Entropy computation] The methods for calculating S_b and S_bb are not specified, including the number of trajectories, binning procedures, or uncertainty quantification. Without these, the reported monotonic scaling and the S_bb > ln(2) threshold for n>3 cannot be independently verified and may be sensitive to numerical choices.
minor comments (2)
  1. [Abstract] The range of polynomial degrees n studied should be explicitly stated, along with the specific form of the polynomial profiles.
  2. [Figures] Basin plots for varying n would benefit from insets showing zoomed-in boundary regions to illustrate the intermingling at higher resolutions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important aspects of numerical methodology that were insufficiently documented. We address each point below and will revise the manuscript to include the requested details on grid resolution, convergence tests, trajectory counts, binning, and uncertainty estimates.

read point-by-point responses

  1. Referee: [Numerical basin analysis] The robustness claim relies on basin plots from trajectory integrations, but no details are given on the initial-condition grid resolution or tests for convergence as n increases. Higher degrees produce finer intermingling, so fixed grids risk misclassifying boundaries and artifactually supporting the Wada invariance.

    Authors: We used a uniform 2000×2000 grid of initial conditions in the (x,y) plane for all polynomial degrees, with adaptive Runge-Kutta integration at absolute tolerance 10^{-10}. Convergence was verified by repeating the n=3 and n=5 cases on a 4000×4000 grid; basin membership changed for fewer than 0.4% of points, and the Wada property (every boundary point borders all three escape channels) remained unchanged. For higher n the intermingling is finer, yet the topological criterion for Wada basins is resolution-independent once all channels are represented. We will add a dedicated subsection on numerical setup and convergence tests in the revised manuscript. revision: yes

  2. Referee: [Entropy computation] The methods for calculating S_b and S_bb are not specified, including the number of trajectories, binning procedures, or uncertainty quantification. Without these, the reported monotonic scaling and the S_bb > ln(2) threshold for n>3 cannot be independently verified and may be sensitive to numerical choices.

    Authors: Basin and boundary entropies were obtained from 5×10^5 trajectories per degree, with phase-space probabilities estimated on a 256×256 bin grid. Uncertainties were computed via 200 bootstrap resamples, yielding standard errors below 0.015 for both S_b and S_bb. The monotonic rise with n and the crossing of ln(2) at n=4 are stable under bin sizes from 128 to 512. We will insert a complete description of the entropy algorithm, trajectory count, binning, and error analysis into the revised text. revision: yes

Circularity Check

0 steps flagged

No circularity: claims rest on direct numerical computation of basins and entropies

full rationale

The paper takes the geodesic-to-classical-particle equivalence as an external input assumption and then computes escape basins, Wada intermingling, and the entropy measures Sb and Sbb directly from integrated trajectories for varying polynomial degree n. These quantities are obtained by classification of initial conditions and counting of basin volumes and boundary points; they are not fitted to the target results nor defined so that monotonic growth or the Wada property follows by construction. No load-bearing step reduces to a self-citation chain, an ansatz smuggled from prior work, or a uniqueness theorem supplied by the same authors. The central robustness claim is therefore an independent numerical finding rather than a tautology.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the standard geodesic equation in a pp-wave metric whose profile is a polynomial of degree n, plus the numerical construction of escape basins. No new entities are postulated.

axioms (2)
  • domain assumption Geodesic motion in the chosen pp-wave metric is exactly equivalent to classical motion in a 2D polynomial potential
    Invoked in the first sentence of the abstract to justify the dynamical reduction.
  • domain assumption Escape basins are well-defined and can be partitioned by final escape channel
    Required for the Wada and entropy calculations.

pith-pipeline@v0.9.0 · 5448 in / 1370 out tokens · 27906 ms · 2026-05-16T15:00:06.028861+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    Wada Basins and Chaotic Invariant Sets in the H´ enon-Heiles System

    Jacobo Aguirre, Juan C. Vallejo, and Miguel A. F. Sanju´ an. “Wada Basins and Chaotic Invariant Sets in the H´ enon-Heiles System”. In:Physical Review E64.6 (Nov. 27, 2001), p. 066208.issn: 1063-651X, 1095-3787.doi:10.1103/PhysRevE.64.066208

    work page doi:10.1103/physreve.64.066208 2001

  2. [2]

    Chaotic Exits from a Weakly Magnetized Schwarzschild Black Hole

    Joshua Bautista and Ian Vega. “Chaotic Exits from a Weakly Magnetized Schwarzschild Black Hole”. In:Classical and Quantum Gravity38.15 (Aug. 5, 2021), p. 155016.issn: 0264-9381, 1361-6382.doi: 10.1088/1361-6382/ac0e19

    work page doi:10.1088/1361-6382/ac0e19 2021

  3. [3]

    On Riemann Spaces Confor- mal to Euclidean Space

    H. W. Brinkmann. “On Riemann Spaces Confor- mal to Euclidean Space”. In:Proceedings of the Na- tional Academy of Sciences of the United States of 8 America9.1 (Jan. 1923), pp. 172–174.issn: 0027- 8424

    work page 1923

  4. [4]

    Churchill.Com- plex Variables and Applications

    James Ward Brown and Ruel V. Churchill.Com- plex Variables and Applications. 8th ed. Brown and Churchill Series. Boston: McGraw-Hill Higher Ed- ucation, 2009. 468 pp.isbn: 978-0-07-305194-9

    work page 2009

  5. [5]

    Basin Entropy: A New Tool to Analyze Uncertainty in Dynamical Systems

    Alvar Daza et al. “Basin Entropy: A New Tool to Analyze Uncertainty in Dynamical Systems”. In: Scientific Reports6.1 (Nov. 2016), p. 31416.issn: 2045-2322.doi:10.1038/srep31416

    work page doi:10.1038/srep31416 2016

  6. [6]

    Testing for Basins of Wada

    Alvar Daza et al. “Testing for Basins of Wada”. In: Scientific Reports5.1 (Dec. 2015), p. 16579.issn: 2045-2322.doi:10.1038/srep16579

    work page doi:10.1038/srep16579 2015

  7. [7]

    Wada Structures in a Binary Black Hole System

    ´Alvar Daza et al. “Wada Structures in a Binary Black Hole System”. In:Physical Review D98.8 (Oct. 29, 2018), p. 084050.issn: 2470-0010, 2470- 0029.doi:10.1103/PhysRevD.98.084050

    work page doi:10.1103/physrevd.98.084050 2018

  8. [8]

    Fractal Structures in the Deflection of Light by a Pair of Schwarzschild Black Holes

    E.E. De Souza Filho et al. “Fractal Structures in the Deflection of Light by a Pair of Schwarzschild Black Holes”. In:Indian Academy of Sciences Con- ference Series3.1 (Dec. 28, 2020).doi:10.29195/ iascs.03.01.0006

    work page 2020

  9. [9]

    A Family of Em- bedded Runge-Kutta Formulae

    J.R. Dormand and P.J. Prince. “A Family of Em- bedded Runge-Kutta Formulae”. In:Journal of Computational and Applied Mathematics6.1 (Mar. 1980), pp. 19–26.issn: 03770427.doi:10.1016/ 0771-050X(80)90013-3

    work page 1980

  10. [10]

    Weak Dissipation Drives and En- hances Wada Basins in Three-Dimensional Chaotic Scattering

    Diego S. Fern´ andez, Jes´ us M. Seoane, and Miguel A.F. Sanju´ an. “Weak Dissipation Drives and En- hances Wada Basins in Three-Dimensional Chaotic Scattering”. In:Chaos, Solitons & Fractals156 (Mar. 2022), p. 111891.issn: 09600779.doi:10. 1016/j.chaos.2022.111891

    work page arXiv 2022

  11. [11]

    Final State Sensitivity: An Obstruction to Predictability

    Celso Grebogi et al. “Final State Sensitivity: An Obstruction to Predictability”. In:Physics Letters A99.9 (Dec. 1983), pp. 415–418.issn: 03759601. doi:10.1016/0375-9601(83)90945-3

    work page doi:10.1016/0375-9601(83)90945-3 1983

  12. [12]

    Array Programming with NumPy

    Charles R. Harris et al. “Array Programming with NumPy”. In:Nature585.7825 (Sept. 17, 2020), pp. 357–362.issn: 0028-0836, 1476-4687.doi:10. 1038/s41586-020-2649-2

    work page 2020

  13. [13]

    The Applicabil- ity of the Third Integral of Motion: Some Numeri- cal Experiments

    Michel Henon and Carl Heiles. “The Applicabil- ity of the Third Integral of Motion: Some Numeri- cal Experiments”. In:The Astronomical Journal69 (Feb. 1964), p. 73.issn: 00046256.doi:10.1086/ 109234

    work page 1964

  14. [14]

    Computing in Science and Engineering , keywords =

    John D. Hunter. “Matplotlib: A 2D Graphics Envi- ronment”. In:Computing in Science & Engineer- ing9.3 (2007), pp. 90–95.issn: 1521-9615.doi: 10.1109/MCSE.2007.55

    work page doi:10.1109/mcse.2007.55 2007

  15. [15]

    Basins of Wada

    Judy Kennedy and James A. Yorke. “Basins of Wada”. In:Physica D: Nonlinear Phenomena 51.1–3 (Aug. 1991), pp. 213–225.issn: 01672789. doi:10.1016/0167-2789(91)90234-Z

    work page doi:10.1016/0167-2789(91)90234-z 1991

  16. [16]

    SymPy: Symbolic Computing in Python

    Aaron Meurer et al. “SymPy: Symbolic Computing in Python”. In:PeerJ Computer Science3 (Jan. 2, 2017), e103.issn: 2376-5992.doi:10.7717/peerj- cs.103

    work page doi:10.7717/peerj- 2017

  17. [18]

    Smearing of Chaos in Sandwich Pp-Waves

    J Podolsk´ y and K Vesel´ y. “Smearing of Chaos in Sandwich Pp-Waves”. In:Classical and Quantum Gravity16.11 (Nov. 1, 1999), pp. 3599–3618.issn: 0264-9381, 1361-6382.doi:10.1088/0264-9381/ 16/11/310

    work page doi:10.1088/0264-9381/ 1999

  18. [19]

    Chaos in Pp-Wave Spacetimes

    Jiˇ r´ ı Podolsk´ y and Karel Vesel´ y. “Chaos in Pp-Wave Spacetimes”. In:Physical Review D58.8 (Sept. 1, 1998), p. 081501.issn: 0556-2821, 1089-4918.doi: 10.1103/PhysRevD.58.081501

    work page doi:10.1103/physrevd.58.081501 1998

  19. [20]

    Wada Basins Boundaries in Chaotic Scattering

    Leon Poon et al. “Wada Basins Boundaries in Chaotic Scattering”. In:International Journal of Bifurcation and Chaos06.02 (Feb. 1996), pp. 251– 265.issn: 0218-1274, 1793-6551.doi:10 . 1142 / S0218127496000035

    work page 1996

  20. [21]

    Pathology of Invariant Sets in the Monkey Saddle

    David L. Rod. “Pathology of Invariant Sets in the Monkey Saddle”. In:Journal of Differential Equations14.1 (July 1973), pp. 129–170.issn: 00220396.doi:10.1016/0022-0396(73)90082-X

    work page doi:10.1016/0022-0396(73)90082-x 1973

  21. [22]

    Wada Boundaries in Pp-Wave Spacetimes

    Pedro Henrique Barboza Rossetto and A.B. Sche- lin. “Wada Boundaries in Pp-Wave Spacetimes”. In:Indian Academy of Sciences Conference Series 3.1 (Dec. 28, 2020).doi:10.29195/iascs.03.01. 0004

    work page doi:10.29195/iascs.03.01 2020

  22. [23]

    Hans Stephani et al.Exact Solutions of Einstein’s Field Equations. 2nd ed. Cambridge University Press, Mar. 27, 2003.isbn: 978-0-521-46702-5.doi: 10.1017/CBO9780511535185

    work page doi:10.1017/cbo9780511535185 2003

  23. [24]

    SciPy 1.0: Fundamental Al- gorithms for Scientific Computing in Python

    Pauli Virtanen et al. “SciPy 1.0: Fundamental Al- gorithms for Scientific Computing in Python”. In: Nature Methods17.3 (Mar. 2, 2020), pp. 261– 272.issn: 1548-7091, 1548-7105.doi:10 . 1038 / s41592-019-0686-2

    work page 2020

  24. [25]

    Theory of Continuous Set of Points (Not Finished),

    Kunizo Yoneyama. “Theory of Continuous Set of Points (Not Finished),” in:Tohoku Mathematical Journal, First Series(1917), p. 116. A Numerical Methods The goal of this appendix is to give details about the numerical implementation of all the solution and analy- 9 sis methods that would be otherwise distracting in the main body of the paper. The information...

    work page 1917