Improving the inference of the stellar quantities using the extended $I$-Love-$Q$-$δM$ relations

In relativistic Astrophysics the $I$-Love-$Q$ relations refer to approximately EoS-independent relations involving the moment of inertia, Love number, and quadrupole moment through some quantities that are normalised by the mass $M_0$ of the background configuration of the perturbative scheme. Since $M_0$ is not an observable quantity, this normalisation hinders the direct applicability of the relations. A common remedy assumes that $M_0$ coincides with the actual mass of the star $M_S$; however, this approximation is only adequate for very slow rotation (when the dimensionless spin parameter is $χ_S<0.1$). The more accurate alternative approach, based on the $I$-Love-$Q$-$δM$ set of relations, circumvents this limitation by enabling the inference of $M_0$. Here we review both approaches and provide numerical comparisons.

💡 Research Summary

The paper addresses a practical limitation of the widely used I‑Love‑Q universal relations for rotating compact stars. These relations connect the moment of inertia (I), the tidal Love number (k₂, often expressed through the dimensionless tidal deformability λ), and the quadrupole moment (Q) via dimensionless quantities that are normalized by the background mass M₀ of the non‑rotating spherical configuration. While M₀ is a convenient theoretical parameter, it is not directly observable; astrophysical measurements provide the actual stellar mass Mₛ and the dimensionless spin χₛ = Iₛ Ωₛ / Mₛ². In most applications the standard approach simply assumes M₀ ≈ Mₛ, an approximation that is only justified for very slowly rotating stars (χₛ ≲ 0.1). For faster rotators, the discrepancy between M₀ and Mₛ introduces systematic errors that break the apparent universality of the I‑Love‑Q relations.

To overcome this problem the authors introduce an extended set of universal relations that also involve the mass‑correction term δM, defined through the second‑order Hartle‑Thorne perturbative expansion as δM = (Mₛ − M₀)/(Ωₛ² M₀³ I²). The key idea is that, given a measurement of the tidal deformability λ (or equivalently k₂), one can use a δM‑Love relation to infer δM. Once δM is known, the true background mass follows from the algebraic relation M₀ = Mₛ / (1 + δM χₛ²). With M₀ in hand, the standard I‑Love‑Q relations can be applied to obtain the dimensionless I and Q, and the physical quantities are then reconstructed as Iₛ = M₀³ I and Qₛ = χₛ² M₀³ Q. This “extended approach” eliminates the need for the M₀ ≈ Mₛ assumption and therefore remains accurate even for moderately fast rotation.

The authors test both the standard and extended methods using a simple polytropic equation of state (P = 100 ρ²). They compute the exact values of M₀, Iₛ, and Qₛ from the full Hartle‑Thorne perturbative solution for a grid of central pressures P_c and spin parameters χₛ. Relative errors are defined as the absolute difference between the inferred and exact values, normalized by the exact value. The results, displayed in Figure 1, show that the extended approach consistently yields smaller errors across the entire parameter space. The advantage becomes more pronounced as χₛ increases: for χₛ ≈ 0.4 the standard method can produce errors exceeding 10 % in M₀, Iₛ, and Qₛ, whereas the extended method keeps errors below a few percent. The authors note that the residual scatter in the extended method’s errors likely reflects intrinsic noise in the universal δM‑Love relation itself, but even this scatter is smaller than the systematic bias introduced by the M₀ ≈ Mₛ approximation.

In the discussion, the paper emphasizes that the extended I‑Love‑Q‑δM framework enables a more reliable inference of stellar properties from gravitational‑wave observations, where λ and χ can be extracted from the inspiral waveform. By providing a way to recover the unobservable background mass, the method opens the possibility of applying universal relations to a broader class of neutron stars, including millisecond pulsars and other rapidly rotating objects that are now being observed in gravitational‑wave data. The authors suggest that future work should explore a wider variety of realistic equations of state, incorporate higher‑order rotational corrections, and refine the δM‑Love fit to further reduce the remaining uncertainties.

Overall, the paper makes a clear methodological contribution: it identifies the source of error in the traditional use of I‑Love‑Q relations, proposes a mathematically consistent extension that includes the mass‑correction term, and demonstrates through numerical experiments that this extension markedly improves the accuracy of inferred stellar quantities, especially for stars with non‑negligible spin. This work therefore represents a significant step toward more precise astrophysical inference from multimessenger observations of compact objects.