Orbit classification in a pseudo-Newtonian Copenhagen problem with Schwarzschild-like primaries
Pith reviewed 2026-05-24 18:12 UTC · model grok-4.3
The pith
In a model of two black holes, the energy of a test particle and the Schwarzschild radius together set the final states of orbits and the degree of fractality in their boundaries.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors classify orbits in the planar pseudo-Newtonian Copenhagen problem with Schwarzschild-like primaries and show that the Jacobi constant and the Schwarzschild radius highly influence the character of the orbits as well as the degree of fractality of the dynamical system, demonstrated by basin diagrams that reveal the distribution of final states across multiple parameter planes.
What carries the argument
Basin diagrams on multiple two-dimensional planes that partition initial conditions according to the final state of each orbit (escape, collision, or bounded motion) while the Jacobi constant and Schwarzschild radius are varied.
If this is right
- Increasing the Jacobi constant changes the relative sizes of escape and collision basins.
- Raising the Schwarzschild radius alters both the locations of regular regions and the fractal dimension of the boundaries.
- The same diagrams show systematic differences from the corresponding Newtonian Copenhagen problem.
- The fractality measure itself becomes a function of both energy and the relativistic parameter.
Where Pith is reading between the lines
- The model implies that modest changes in orbital energy could move a real black-hole binary from mostly regular to highly mixed escape regions.
- The observed dependence on Schwarzschild radius suggests that full general-relativity simulations of comparable systems would need to resolve similar fractal boundaries at small separations.
- The basin-diagram technique could be applied directly to other pseudo-Newtonian or post-Newtonian binary potentials to test robustness of the fractality result.
Load-bearing premise
The chosen pseudo-Newtonian potential gives a sufficiently accurate picture of Schwarzschild gravity for determining the final states of test-particle orbits.
What would settle it
Numerical orbit integrations performed in full general relativity for identical initial conditions and the same range of energies would produce different basin structures or different degrees of fractality.
read the original abstract
We examine the orbital dynamics of the planar pseudo-Newtonian Copenhagen problem, in the case of a binary system of Schwarzschild-like primaries, such as super-massive black holes. In particular, we investigate how the Jacobi constant (which is directly connected with the energy of the orbits) influences several aspects of the orbital dynamics, such as the final state of the orbits. We also determine how the relativistic effects (i.e., the Schwarzschild radius) affect the character of the orbits, by comparing our results with the classical Newtonian problem. Basin diagrams are deployed for presenting all the different basin types, using multiple types of planes with two dimensions. We demonstrate that both the Jacobi constant as well as the Schwarzschild radius highly influence the character of the orbits, as well as the degree of fractality of the dynamical system.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the planar pseudo-Newtonian Copenhagen problem for equal-mass Schwarzschild-like primaries. Using numerical forward integration of test-particle trajectories, it classifies orbits by final state (collision, escape, or bounded) and constructs basin diagrams in several two-dimensional planes. The central claim is that both the Jacobi constant and the Schwarzschild-radius parameter strongly affect the character of the orbits and the degree of fractality of the basins, with explicit comparisons to the Newtonian limit.
Significance. If the pseudo-Newtonian model is shown to be sufficiently faithful, the work would provide concrete evidence that relativistic corrections alter both the topology of escape/collision basins and their fractal properties in equal-mass binary systems. The deployment of basin diagrams across multiple planes is a useful visualization choice that strengthens the presentation of the numerical results.
major comments (3)
- [§3] §3 (numerical setup): the manuscript provides no information on the integrator, tolerance settings, step-size control, or any convergence/accuracy tests used to classify orbits. Because the reported basin structures and fractality measures rest entirely on these classifications, the absence of such verification leaves the central numerical claims without verifiable support.
- [§2] §2 (model definition): the chosen pseudo-Newtonian potential is introduced without any error estimate or cross-check against geodesic motion in the corresponding binary Schwarzschild spacetime (or even against post-Newtonian expansions). The claim that the Schwarzschild radius “highly influences” orbit character and fractality therefore lacks a demonstrated link to full general-relativistic dynamics.
- [Results] Results section: the phrase “degree of fractality” is used repeatedly but never accompanied by a quantitative measure (box-counting dimension, uncertainty exponent, etc.). Without an explicit metric, it is impossible to assess how the reported changes with Jacobi constant or Schwarzschild radius are evaluated.
minor comments (2)
- Figure captions should state the precise ranges of the plotted variables and the number of initial conditions used to generate each basin diagram.
- Notation for the pseudo-Newtonian potential and the effective Schwarzschild radius should be introduced once in §2 and used consistently thereafter.
Simulated Author's Rebuttal
We thank the referee for the constructive comments, which help improve the clarity and rigor of the manuscript. We address each major comment below and indicate the revisions we will make.
read point-by-point responses
-
Referee: §3 (numerical setup): the manuscript provides no information on the integrator, tolerance settings, step-size control, or any convergence/accuracy tests used to classify orbits. Because the reported basin structures and fractality measures rest entirely on these classifications, the absence of such verification leaves the central numerical claims without verifiable support.
Authors: We agree that the numerical methods were insufficiently documented. In the revised manuscript we will add a dedicated subsection in §3 describing the integrator (a 7-8th order Runge-Kutta-Fehlberg scheme with adaptive step-size control), the absolute and relative tolerance settings (both set to 10^{-12}), the criteria used to classify final states, and the results of convergence tests in which tolerances were varied by two orders of magnitude and basin membership was verified to remain unchanged for the great majority of initial conditions. revision: yes
-
Referee: §2 (model definition): the chosen pseudo-Newtonian potential is introduced without any error estimate or cross-check against geodesic motion in the corresponding binary Schwarzschild spacetime (or even against post-Newtonian expansions). The claim that the Schwarzschild radius “highly influences” orbit character and fractality therefore lacks a demonstrated link to full general-relativistic dynamics.
Authors: The pseudo-Newtonian potential employed is a standard approximation in the literature for incorporating leading relativistic corrections into the restricted three-body problem. We will expand §2 to include a brief discussion of its derivation, cite existing single-body validations against post-Newtonian expansions, and explicitly qualify the claim by noting that the reported influence is within the pseudo-Newtonian framework rather than a direct demonstration of full geodesic dynamics. A quantitative error estimate against the binary Schwarzschild spacetime is not feasible within the scope of this work. revision: partial
-
Referee: Results section: the phrase “degree of fractality” is used repeatedly but never accompanied by a quantitative measure (box-counting dimension, uncertainty exponent, etc.). Without an explicit metric, it is impossible to assess how the reported changes with Jacobi constant or Schwarzschild radius are evaluated.
Authors: We accept that the fractality discussion remained qualitative. In the revised Results section we will introduce the uncertainty exponent as a quantitative measure of basin fractality, compute its value for representative slices in the (x,y), (C,J) and (C, r_s) planes, and report how it varies with Jacobi constant and Schwarzschild radius, thereby replacing the previous descriptive statements with explicit numerical values. revision: yes
- A direct numerical cross-check of the pseudo-Newtonian trajectories against geodesic motion in the full binary Schwarzschild spacetime.
Circularity Check
No circularity: numerical orbit classification uses independent forward integration
full rationale
The paper performs direct numerical integration of the equations of motion in a fixed pseudo-Newtonian potential (with Jacobi constant and Schwarzschild radius as independent inputs) and classifies final states via basin diagrams on multiple planes. No parameters are fitted to the classified data, no equations reduce reported basin structures or fractality measures to quantities defined from the same outputs, and no self-citations are invoked as load-bearing uniqueness theorems or ansatzes. The comparison to the Newtonian case is a straightforward parameter variation. The derivation chain is therefore self-contained and does not reduce to its inputs by construction.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.