Universal Relations for Neutron Stars from Asymptotic Analysis

Abstract: Dimensionless observables of neutron stars, such as the moment of inertia, the tidal deformability, the spin-induced quadrupole moment, and the compactness, satisfy the I--Love--Q and Love--$C$ universal relations to percent-level accuracy over a wide range of equations of state. We investigate analytically the origin of this insensitivity in the stellar structure equations. We derive asymptotic expansions of these relations directly from the differential equations and boundary conditions for slowly rotating, tidally deformed stars described by a general piecewise-polytropic equation of state. Although the differential equations allow several forms of the asymptotic expansion, we find that only one class is consistent with the observed universality, and we use this class to analyze the universal relations. Because the observables are determined at the stellar surface, information about the deep-interior equation of state can enter them only through the parameters that survive in the asymptotic expansion at the surface. We find that the universal relations depend only on three parameters: two integration constants and one polytrope index of the outermost segment. By determining how these parameters deform the relations, we show that, throughout the region of parameter space occupied by realistic equations of state, the resulting deviations remain at the percent level, consistent with the observed accuracy of the universal relations. Within the same formalism, we classify violations of universality according to the additional degrees of freedom or input data responsible for them. Since the construction relies only on the underlying differential equations and boundary conditions, the same procedure can be applied to other systems once the corresponding equations and boundary conditions are specified.

• The paper derives I–Love–Q and Love–C relations from asymptotic analysis of relativistic stellar equations, showing that surface observables depend on a small set of segment parameters and that realistic polytropic indices yield approximately 1% accuracy.
• The analysis identifies memory loss as the mechanism suppressing deep-interior equation-of-state effects, with the parameter contribution scaling as the compactness factor $(c_t/C)^{3/2}$ and diminishing at higher compactness.
• The paper classifies universality violations from strong phase transitions, anisotropy, magnetic fields, elasticity, and additional matter or gravitational degrees of freedom, providing diagnostic clues for interpreting departures from standard relations.

The I–Love–Q relations connect the dimensionless moment of inertia $\bar I$, tidal deformability $\bar\lambda$, and spin-induced quadrupole moment $\bar Q$ of neutron stars with equation-of-state (EoS) insensitivity at the percent level, together with the companion Love–$C$ relation tying $\bar\lambda$ to the compactness $C$. While these relations have been verified numerically across broad EoS ensembles, a first-principles explanation of why EoS dependence is suppressed has remained incomplete, with earlier analytic treatments restricted to single polytropes, phenomenological density profiles, or Newtonian expansions. The paper by Kamata, Minamiguchi, and Minato addresses this gap through an asymptotic analysis of the full relativistic stellar structure equations for slowly rotating, tidally deformed stars described by piecewise-polytropic EoS (2608.19939).

Formulation

The authors work within the standard perturbative framework of Hartle slow rotation and static tides about Tolman–Oppenheimer–Volkoff backgrounds. They adopt $c=m/r$ as a flow variable, which vanishes at the center, equals the compactness $C=M/R$ at the surface, and is bounded by the black-hole limit $c<1/2$, making it a natural organizing parameter. The observables $\bar I$, $\bar\lambda$0, and $\bar\lambda$1 are shown to be determined entirely by surface quantities and exterior matching data; in particular, the tidal Love number follows from the single interior combination $\bar\lambda$2 through the closed-form expression of Hinderer and Damour–Nagar.

Asymptotic structure and the three segment parameters

For each polytrope segment, the TOV equations carry two integration constants beyond those fixed by central regularity. These are packaged as $\bar\lambda$3, where $\bar\lambda$4 measures the mass offset relative to the regular solution of that segment and accumulates outward across interfaces — it is the only channel through which inner-segment information reaches the surface. A leading-order shift $\bar\lambda$5 can be absorbed into a redefinition of the central pressure.

A key methodological step is the classification of possible multivariable expansions in the three scales $\bar\lambda$6, $\bar\lambda$7, and $\bar\lambda$8. Six orderings collapse into four inequivalent expansion bases. Only basis I — where both $\bar\lambda$9 and $\bar Q$0 are small at the surface — yields relations whose log–log slope is fixed by the observable definitions alone ($\bar Q$1, $\bar Q$2, $\bar Q$3 for the three I–Love–Q pairs), matching the observed near-constant slopes of the universal relations. Bases II–IV instead produce $\bar Q$4-independent intercepts set by the segment constants, with no universal logarithmic structure. This provides a selection criterion grounded purely in the differential equations: the known universality singles out basis I without any microscopic restriction on the EoS.

Within basis I, Lagrange inversion of the surface condition produces a master series relating any pair of observables, with all EoS dependence entering through three parameters: the outermost segment's polytrope index $\bar Q$5 and the two integration constants. A tangential/transverse decomposition shows that $\bar Q$6 displaces stars exactly along the universal curve (it is a reparametrization of central pressure), so transverse deviations are controlled only by $\bar Q$7 and $\bar Q$8.

Polytrope-index dependence

Reducing to single-polytrope stars, the coefficients $\bar Q$9 of the I–Love relation are computed by quadrature through sixth post-Newtonian order over the stable range $C$0. Over the realistic core band $C$1 — motivated by the posterior ensemble of Fujimoto et al., whose median effective indices remain below 0.46 up to densities $C$2 — the half-spread of $C$3 is only 1.0–1.2% for $C$4–7.4, whereas it grows to 9–14% over the full stable range. Higher-order coefficients steepen sharply as $C$5 (e.g., $C$6 changes by 82% at $C$7), but these enter at suppressed powers of $C$8 and contribute little. The analytic relation agrees with numerical sequences for SLy, APR4, and QHC21 and reproduces the Yagi–Yunes fit to better than 1%. The weak $C$9 sensitivity for realistic matter is thus tied to the empirical fact that sampled indices stay well below the stability bound — a point the paper states plainly rather than deriving from deeper principles.

Memory loss and its failure at strong phase transitions

The $\bar\lambda$0 contribution enters as $\bar\lambda$1, suppressed at high compactness — a mechanism the authors call memory loss: deep-interior information progressively decouples from surface observables. A sharp Maxwell first-order transition generates $\bar\lambda$2 proportional to the density jump relative to the mean enclosed energy density, transmitted outward with suppression factor $\bar\lambda$3. Sizable breaking requires both a large discontinuity and a transition occurring at $\bar\lambda$4. Numerical two-segment tests confirm that sequences with $\bar\lambda$5 of order unity depart from the single-polytrope band, then return toward it at higher compactness, consistent with memory loss. For negative $\bar\lambda$6, universality also breaks but via continuation into bases III/IV, where the definition-fixed slope is absent. Within this framework, an observed deviation from universality could indicate either a near-boundary polytrope index or a strong phase transition, though a quantitative criterion in terms of $\bar\lambda$7 is left open.

Extensions and taxonomy of violations

The paper demonstrates robustness against generalizing the piecewise-polytrope form to power-series corrections within segments: higher coefficients enter with additional factors $\bar\lambda$8 and do not alter the memory counting. For static perfect-fluid stars, composition stratification introduces no new independent information because the integrability condition $\bar\lambda$9 forces the static response to depend only on the equilibrium sound speed — though this argument fails at finite frequency (where $C$0 modes appear) and for rotating stars, where barotropy is required for a first integral of the Euler equation.

Anisotropic stress, magnetic fields, and elasticity are treated schematically, yielding a three-type classification of universality violations: Type I, growth of the ordinary $C$1 channel (strong first-order transitions, twin stars); Type IIa, additional parameters or sources modifying the equations (anisotropy parameters, field strength, spin, modified-gravity couplings); Type IIb, additional free constants after equations and sources are fixed (one-mode stress channels, dark-matter admixture); and Type III, additional radial functions (magnetic-field configurations, differential rotation, frozen composition profiles). Fluid envelopes screen internal elastic stress analogously to how outer EoS segments screen deep-interior structure, consistent with prior results on solid-core/fluid-envelope configurations.

Limitations

Several restrictions bear directly on the quantitative claims. The analysis is confined to slowly rotating, statically tidally deformed, barotropic-response perfect fluids; rapid rotation, dynamical tides, and superfluidity enter only through the taxonomy. The neglect of $C$2 is justified only because its transverse effect starts at quadratic order $C$3. The realistic range $C$4 is inferred from one particular posterior ensemble and could shift with future data. The treatment of Gibbs transitions, mixed phases, and elastic constitutive models remains schematic, with quantitative corrections deferred to solving the full coupled perturbation equations. Finally, the claim that sizable breaking requires $C$5 rests on a leading-order estimate valid when $C$6.

Conclusion

This work supplies an analytic derivation of the I–Love–Q and Love–$C$7 universality from the stellar structure equations themselves, identifying the precise parametric channel — two integration constants and one polytrope index per segment — through which EoS information reaches surface observables, and showing that percent-level accuracy follows from the smallness of realistic $C$8 and $C$9. The classification of violation mechanisms by how extra degrees of freedom enter the radial problem offers a transferable criterion, applicable in principle to other self-gravitating systems once their equations and boundary conditions are specified.