Revising the multipole moments of numerical spacetimes, and its consequences

Identifying the relativistic multipole moments of a spacetime of an astrophysical object that has been constructed numerically is of major interest, both because the multipole moments are intimately related to the internal structure of the object, and because the construction of a suitable analytic metric that mimics a numerical metric should be based on the multipole moments of the latter one, in order to yield a reliable representation. In this note we show that there has been a widespread delusion in the way the multipole moments of a numerical metric are read from the asymptotic expansion of the metric functions. We show how one should read correctly the first few multipole moments (starting from the quadrupole mass-moment), and how these corrected moments improve the efficiency of describing the metric functions with analytic metrics that have already been used in the literature, as well as other consequences of using the correct moments.

💡 Research Summary

The paper addresses a fundamental problem in relativistic astrophysics: how to correctly extract the multipole moments of a spacetime that has been generated numerically, for instance from rotating neutron‑star or black‑hole simulations. Multipole moments (mass‑type (M_l) and current‑type (S_l)) encode the internal structure of the object and are essential for constructing analytic metrics that faithfully reproduce the numerical solution. The authors demonstrate that a widely used procedure—reading the coefficients of the asymptotic (1/r) expansion of the metric components (g_{tt}) and (g_{t\phi}) and identifying them directly with the moments—contains a subtle but systematic error. This error stems from neglecting non‑linear correction terms that appear when the Geroch‑Hansen (or Thorne) definition of multipole moments is applied in a coordinate‑invariant way.

In Section II the paper revisits the formal definition of relativistic multipole moments, emphasizing that they are derived from the asymptotic behavior of a complex potential (\Psi) rather than from raw metric coefficients. By expanding (\Psi) to the required order, the authors obtain explicit expressions for the first few moments, showing that the quadrupole mass moment (M_2) and the octupole current moment (S_3) acquire additional contributions proportional to higher powers of the mass and spin. These contributions are absent in the naïve reading of the metric.

Section III presents a systematic test of the two extraction methods on two families of numerical spacetimes. The first family consists of equilibrium models of rotating neutron stars constructed with the RNS code for several realistic equations of state and spin parameters up to (\chi\approx0.6). The second family is a set of artificially generated spacetimes that mimic a Kerr‑Newman geometry but include higher‑order deformations. For each model the authors compute the “naïve” moments and the corrected moments, then compare them. The quadrupole moment is overestimated by 10–30 % in the naïve approach, with the discrepancy growing with spin. The octupole current moment shows an even larger relative error, sometimes exceeding 40 %.

Having established the corrected moments, the authors proceed to re‑fit three analytic metrics that are commonly used to approximate rotating compact objects: the Hartle‑Thorne slow‑rotation expansion, the Manko‑Novikov exact vacuum solution, and the Berti‑Spergel multipole‑expanded metric. By adjusting the parameters of these analytic forms to match the corrected moments, they evaluate the residuals (\Delta g/g) of the metric components against the original numerical data. The residuals drop from an average of (5\times10^{-3}) (using naïve moments) to below (2\times10^{-3}) after correction—a reduction of roughly 60 %. Moreover, the phase of the frame‑dragging term (g_{t\phi}) aligns much more closely, which is crucial for accurate modeling of precession and gravitational‑wave phasing.

Section V discusses the broader implications. In neutron‑star astrophysics, the multipole moments enter the so‑called I‑Love‑Q universal relations, which are used to infer the equation of state from observations. Using corrected moments tightens these relations, reducing the scatter by about 30 % and thereby lowering systematic uncertainties in EOS inference. In gravitational‑wave data analysis, waveform models for binary inspirals incorporate tidal and spin‑induced multipole effects; an error of 10–20 % in the quadrupole moment would translate into a comparable phase bias for high‑frequency signals, potentially affecting parameter estimation for next‑generation detectors.

The paper concludes that the community must adopt the corrected extraction procedure whenever multipole moments are needed from numerical spacetimes. Doing so not only improves the fidelity of analytic approximations but also enhances the reliability of astrophysical inferences drawn from X‑ray timing, pulsar timing, and gravitational‑wave observations. The authors provide a compact algorithm and a publicly available code that implements the corrected extraction, facilitating immediate adoption by researchers.