An arXiv version-one preprint dated 20 August 2026 offers a mathematical explanation for why neutron-star universal relations can be relatively insensitive to the equation of state—the model used to describe the star’s matter—while showing that sharp first-order phase transitions can be associated with larger departures from the relations in the model. The authors interpret the pattern as memory loss: only a limited amount of information from the star’s interior is transmitted through the outermost segment to surface observables, and the surviving contribution is reduced as compactness, a measure of how concentrated the star is, increases.
The paper examines two families of relations. I–Love–Q links a star’s dimensionless moment of inertia, its tidal deformability—how readily it is distorted by tides—and its spin-induced quadrupole moment; Love–C links tidal deformability with compactness. The asymptotic construction is not tied to a particular observable, so the same differential-equation system can be used to study both.
Why the relations look universal
The researchers derive asymptotic expansions—controlled series for the behavior of the stellar equations in a defined limit—directly from the differential equations and boundary conditions for slowly rotating, tidally deformed stars described by a general piecewise-polytropic equation of state. For the single-polytrope analysis, the coefficient functions are obtained by numerical integration through the sixth post-Newtonian order.
After imposing the required surface behavior, the main expansion keeps only basis I, where the surface values σp(C) and σr(C) are both small compared with one. In that regime, the observable expansions are controlled by three outer-segment parameters: the integration constants Σp and Σr, transmitted outward by each segment, and the polytrope index γ, a parameter describing the outermost segment.
One part of the expansion is especially restrictive. At Σr = 0, the relation is independent of Σp to all orders. Terms involving Σp that move the result across the relation must also contain Σr, and they begin at quadratic order.
For moderate Σr, the contribution carrying inner-equation-of-state information is suppressed with increasing compactness as C^-3/2. The authors call this memory loss: in the model, greater compactness means less influence from the inner structure on the surface relation.
Small shifts in the modeled realistic range
The size of the modeled variation depends on the range of outer-segment models included. For γ from 0.2 to 0.6 and the logarithm of tidal deformability from 5.8 to 7.4, the half-spread in the logarithm of dimensionless moment of inertia is 1.0% to 1.2%. Across the full stable range, 0 ≤ γ < 3/4, the half-spread is 9% to 14%; the stability analysis treats γ = 3/4 as marginal and γ = 0 separately as the uniform-density endpoint.
Those percentages are model-based, not observational confidence intervals. For its realism check, the paper uses an equation-of-state ensemble represented by five logarithmically spaced segments from ε0 = 150 MeV/fm³ to 8ε0. The ensemble was inferred from 20 observed sources, and the paper reports 3 × 10^5 equations of state.
Where the relation can move
The paper’s sharp first-order phase-transition test links a sudden density discontinuity with a large Σr. The leading contribution is set by the density jump relative to the mean density enclosed by the transition, and the effect is less suppressed when the transition lies closer to the surface.
Numerical sequences show larger displacements for larger density discontinuities. At larger central pressures, the sequences turn back toward the single-polytrope band, the reference band from models without a transition.
The expansion has a reported breakdown condition, C ≲ |Σr|^(2/3). The paper presents this as a leading-order estimate, not a general quantitative threshold for the transition pressure and density jump associated with significant breaking.
Within this framework, a departure does not identify one explanation by itself. The paper classifies possible violations as growth of the ordinary Σr channel, modifications to the equations or sources, additional free constants, or an independent radial function.
A specific model, not a universal guarantee
The boundaries of the analysis matter. Its main universal-relation expansion is restricted to basis I, where surface σp and σr are small, and the underlying model is for slowly rotating, tidally deformed stars with piecewise-polytropic equations of state. The phase-transition result relies on the leading-order estimate above, so it should be read as a model condition rather than a universal cutoff.
Open questions include how the memory-loss mechanism behaves in more general stellar models, what transition parameters predict significant breaking in realistic models, and whether an added radial function can be reduced to finite modes or must be retained in full. The authors propose applying the same differential-equation and boundary-condition strategy to other systems.
Paper data and sources
Original title: Universal Relations for Neutron Stars from Asymptotic Analysis
Authors: Syo Kamata, Josuke Minamiguchi, Shuhei Minato
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-20
DOI: Not available
Original paper · Full text