Phase-Plane Analysis Solves Weak Gravitational Deflection Across Static Spherical Spacetimes

Theoretical physicists Reggie C. Pantig and Ali Övgün have formulated gravitational light deflection in static spherical spacetimes through the amplitude and angle of a null orbit in its phase plane, offering an alternative mathematical framework that sidesteps traditional perturbative trajectory tracking. Gravitational deflection of light holds an established role in general relativity by connecting spacetime geometry directly with observable propagation effects, tracing back to Albert Einstein’s relativistic predictions and early empirical verification through solar eclipse observations. While conventional weak-deflection derivations determine bending by expanding a nonlinear geodesic equation and tracking how an outgoing trajectory asymptote shifts away from flat-space coordinates, the newly developed phase-plane approach treats the null path as an oscillatory system whose total deflection angle emerges as an accumulated difference between physical azimuth and an intrinsic phase-plane advance across a fixed half-cycle.

To establish the dynamical system, the researchers examined equatorial null scattering in the exterior geometry of a static, spherically symmetric spacetime. Normalizing timelike and rotational Killing vectors in standard Schwarzschild coordinates yields the conserved energy and angular momentum of a passing photon, whose ratio defines an invariant asymptotic impact parameter that prevents coordinate-dependent ambiguities. Differentiating the exact radial energy equation along the equatorial path yields a second-order nonlinear differential equation for the dimensionless inverse areal radius as a function of the physical azimuthal angle. By pairing this inverse radius with its azimuthal derivative, the authors mapped the photon scattering trajectory into a two-dimensional phase plane characterized by a positive phase-plane amplitude and an intrinsic angular phase variable.

In flat Minkowski spacetime, the unperturbed null orbit traces a unit semicircle in this phase plane, where the incoming asymptote corresponds to an intrinsic phase of zero, closest approach occurs at pi over two, and the outgoing asymptote terminates precisely at pi. Under the influence of Schwarzschild gravity, the phase-plane amplitude deforms outward as the light ray nears the massive body, reaching a maximum at closest approach before symmetrically returning to its asymptotic unit radius. By proving that the intrinsic phase angle increases monotonically along every returning null orbit outside the critical photon sphere, the authors established that the intrinsic scattering endpoints remain strictly fixed between zero and pi. Gravitational deflection is therefore expressed entirely as an accumulated phase lag, where the bending angle equals the excess physical azimuth swept out over this invariant intrinsic interval.

A critical step in the investigation was removing the differential evolution of the amplitude through an exact algebraic reduction. For Schwarzschild geometry, the coupled first-order differential equations governing the phase-plane amplitude and phase depend on an identical composite variable composed of the dimensionless weak-scattering parameter and the cube of the sine of the phase. Dividing the amplitude derivative by the phase rate produces an autonomous first-order equation that becomes linear when the dependent and independent variables are inverted. Integrating this relation using an integrating factor reduces the determination of the reciprocal amplitude along the trajectory to the physical branch of a cubic algebraic map, matching the first integral of the geodesic equation and eliminating the need to integrate differential orbit equations.

Exploiting this algebraic structure, the authors applied the Lagrange inversion formula to invert the cubic relation and extract an explicit all-order power series for the reciprocal amplitude displacement. Differentiating this series yields the exact algebraic phase factor, which is represented either as a finite combinatorial sum over binomial coefficients or via a closed algebraic generating equation. When inserted into the fixed-interval bending integral, the calculation separates the nonlinear gravitational dynamics from the global orbital geometry. The local Schwarzschild interaction contributes purely rational algebraic coefficients, while the integration across the fixed intrinsic phase interval contributes universal trigonometric moments expressed through Euler beta and gamma functions.

This clean mathematical factorization explains the alternating analytic structure of the Schwarzschild weak-deflection series, clarifying why even powers of the gravitational parameter carry an explicit factor of pi while odd powers remain purely rational numbers. Comparing the resulting series against established literature, the formulation successfully reproduces the invariant Schwarzschild deflection coefficients through eighth order, matching independent results obtained from inhomogeneous Picard-Fuchs equations, Padé approximants, and post-post-Newtonian perturbation schemes. The authors verified these findings through direct comparisons with the exact radial scattering integral, second-order conventional orbit perturbations, and high-precision numerical evaluations of successive series truncations, noting that their work contributes an amplitude-phase representation rather than newly discovered asymptotic numerical coefficients.

The investigation subsequently established that the phase-plane methodology extends beyond Schwarzschild geometry to general asymptotically flat, static spherical metrics described by arbitrary temporal and radial metric functions. By adopting the invariant areal radius, the general null geodesic equation generates an implicit relation connecting the phase-plane amplitude to the radial metric function. Whenever the radial function exhibits a single nonlinear monomial, the amplitude satisfies an exact algebraic polynomial relation. Applying this general formalism to the charged Reissner-Nordström spacetime produces an exact quartic amplitude equation that correctly reproduces known weak-deflection coefficients through third order in the mass and charge parameters, while recovering the leading post-Newtonian and second-order invariant coefficients for general static spherical compact objects.

The study also adapted the exact phase formulation to realistic astrophysical situations involving static sources and receivers located at finite radial coordinates. Because an observer at a finite distance measures photon angles relative to local orthonormal frames rather than asymptotic directions, truncating the intrinsic phase limits alone cannot produce the correct finite-distance deflection. By incorporating local directional angles alongside the adjusted phase integration boundaries, the authors derived an exact expression for finite-distance bending and recovered the leading-order Schwarzschild finite-distance deflection formula. Independent numerical tests with unequal finite source and receiver radii confirmed that the analytical finite-distance phase relation matches direct numerical integration of the radial geodesic equations.

Analyzing the analytical boundaries of the formulation revealed a direct mathematical connection between the convergence of the weak-field expansion and the onset of strong gravitational deflection. The convergence radius of the weak power series is strictly governed by the nearest square-root branch point of the inverse cubic amplitude map, which occurs at the critical threshold where the dimensionless impact parameter equals three times the square root of three times the Schwarzschild mass. Because this algebraic singularity coincides precisely with the unstable photon sphere where the intrinsic phase rate vanishes, the large-order asymptotic behavior of the weak series coefficients exhibits a harmonic form that matches the known logarithmic divergence of the bending angle as the impact parameter approaches its critical value.

By isolating the regular component of the complete angular integral near the critical photon orbit, the researchers analytically recovered both the unit logarithmic coefficient and the established finite strong-deflection constant for Schwarzschild spacetime, verifying the limiting approach through numerical evaluations. The authors noted that the analysis remains subject to specific physical constraints, including the requirement of static asymptotic geometry, a monotonic areal radius across the scattering domain, returning trajectories with an ordinary radial minimum, and a positive phase-advance rate. In spacetimes characterized by non-polynomial metric functions, multiple photon spheres, or non-monotonic radial coordinates, the amplitude relation remains an implicit relation rather than a single solvable cubic, marking a clear path for future geometric classifications.

Subject of Research: Theoretical physics and gravitational lensing in static spherical spacetimes

Article Title: Phase-plane formulation of weak gravitational deflection in static spherical spacetimes

Article References: Pantig, R. C., & Övgün, A. (2026). Phase-plane formulation of weak gravitational deflection in static spherical spacetimes. The European Physical Journal C, 86(11), Article 1265. https://doi.org/10.1140/epjc/s10052-026-16566-3

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16566-3

Keywords: gravitational lensing, general relativity, Schwarzschild spacetime, light deflection, theoretical physics, phase plane, asymptotic scattering, photon orbit, Phase-plane, formulation, weak, gravitational