Analytical approximation of the exterior gravitational field of rotating neutron stars

It is known that B"acklund transformations can be used to generate stationary axisymmetric solutions of Einstein’s vacuum field equations with any number of constants. We will use this class of exact solutions to describe the exterior vacuum region of numerically calculated neutron stars. Therefore we study how an Ernst potential given on the rotation axis and containing an arbitrary number of constants can be used to determine the metric everywhere. Then we review two methods to determine those constants from a numerically calculated solution. Finally, we compare the metric and physical properties of our analytic solution with the numerical data and find excellent agreement even for a small number of parameters.

💡 Research Summary

The paper presents a novel analytical approximation scheme for the exterior vacuum field of rotating neutron stars, based on multi‑soliton solutions generated by Bäcklund transformations applied to the flat‑space (Minkowski) seed. The authors start from the Ernst formulation of stationary, axisymmetric vacuum spacetimes, where the complex Ernst potential f(ρ,ζ)=e^{2U}+i b satisfies a nonlinear partial differential equation. Crucially, the Ernst potential on the symmetry axis, f(0,ζ), uniquely determines the full solution. By assuming that the axis data can be written as a rational function Z(ζ)/N(ζ), the authors show that the entire spacetime metric can be reconstructed from a finite set of Bäcklund parameters {K_i, α_i}.

The construction proceeds as follows. First, the linear problem associated with Ernst’s equation is expressed as a Riemann‑Hilbert problem for a 2×2 matrix Φ(λ, z, \bar z) depending on a spectral parameter λ. A special ansatz forces Φ to be a finite polynomial in λ, which corresponds to applying n successive Bäcklund transformations. The resulting Ernst potential is given by a quotient of two polynomials in λ, and setting λ=1 yields f(ρ,ζ). The parameters K_i are the zeros of a 2n‑degree polynomial built from Z and N, obtained via the algebraic relation (14). The corresponding α_i are then fixed by (15), which involves evaluating Z and N at the K_i. This algorithm works irrespective of whether the K_i are real or complex and automatically enforces the reality condition of the metric.

Two practical ways to obtain the axis potential are discussed. The first uses the Geroch‑Hansen multipole moments: expanding the function X(ζ)= (1−f)/(1+f) at infinity yields coefficients m_j that are directly related to the physical multipoles (mass, angular momentum, quadrupole, etc.). Substituting these coefficients into the rational ansatz determines Z and N, guaranteeing the correct asymptotic behavior. The second method fits a rational function directly to numerically computed axis data g(ζ) at selected points ζ_i. By choosing enough points, especially near the stellar surface, the fit reproduces the interior‑exterior matching region accurately; however, additional constraints are required to preserve the correct mass and angular momentum at infinity.

The authors implement the scheme for n=1, 2, 3, 4 (i.e., 2, 4, 6, 8 free parameters) and compare the resulting metric functions U, a, and k with high‑precision numerical solutions obtained from the AKM code (Ansorg et al., 2003). The test case is a uniformly rotating neutron star with a constant‑density equation of state and a polar‑to‑equatorial radius ratio of 0.7. Relative errors in e^{2U} and the other potentials drop below 10^{-4} already for n=2, and the agreement improves systematically as n increases. The results demonstrate that even a modest number of Bäcklund parameters yields an analytically tractable metric that reproduces the full numerical spacetime to very high accuracy.

In the discussion, the authors emphasize the advantages of their approach: the analytical form is simple, requires no integration of differential equations, and the parameter extraction is purely algebraic. This makes the method attractive for applications where an explicit metric is needed, such as computing geodesics, modeling accretion disks, studying oscillation modes, or coupling to external electromagnetic fields. Moreover, the framework naturally incorporates reflection symmetry (f(ζ) f(−ζ)=1) which reduces the number of independent parameters by half, and it respects asymptotic flatness by construction.

Overall, the paper provides a clear, systematic pathway from axis data to a full exterior solution via Bäcklund transformations, validates the method against state‑of‑the‑art numerical models, and opens the door to a wide range of analytic investigations of rapidly rotating compact objects.