REVIEW
Quadrupole moments of rotating neutron stars
General Relativity and Quantum Cosmology · Submitted 12 Sep 1997 · 6 references
Pith's verdict: A solid, useful numerical reference whose tables are the real product — the quadratic fit is genuine but only good to about 10%, and the paper never says so out loud.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
Rotating neutron stars have a quadrupole moment that follows a single quadratic law in spin, with a coefficient set by mass and by how stiff the matter is.
The reading
What carries the argument
Two pieces of machinery carry the argument. The first is the invariance proof in Sec. II: for a stationary, axially symmetric, equatorially reflection-symmetric metric in the Bardeen-Wagoner form, the authors compute the angular velocity $Omega$ and the orbital energy per unit mass $tilde{E}$ of a circular equatorial test-particle orbit, expand both in powers of $v = (M Omega)^{1/3}$, and form Ryan's combination $Delta tilde{E} = -(Omega/3), dtilde{E}/dOmega$. Matching the resulting series to the general expansion $Delta tilde{E} = frac{1}{3}v^2 - frac{1}{2}v^4 + frac{20}{9}chi v^5 + (q - frac{27}{8})v^6 + O(v^7)$ confirms $J = chi M^2$ and $Q = q M^3$. That licenses reading $q$ directly off the asymptotic potential, which is what the numerics then do. The second is the extraction itself, done two ways: the asymptotic integral $q = -frac{5}{2} lim_{r to infty} (r/M)^3 int_{-1}^{1} nu(r,theta) P_2(cos theta), dcos theta$, which is clean in principle but numerically delicate because the quadrupole piece is tiny next to the spherical piece; and the volume integral $Q = frac{1}{8pi} int S_rho(r',theta'), r'^2 P_2(cos theta'), d^3x'$, where $S_rho$ is the Komatsu-Eriguchi-Hachisu source function that reduces to $8pi epsilon$ only in the Newtonian limit. The second method is the workhorse, the first a check.
Extended reading notes
Core claim
The central claim is that the quadrupole moment of a rotating neutron star is well described, across the entire allowed range of rotation, by $q = -a(M, EOS) chi^2$, where $chi = cJ/GM^2$ is the dimensionless spin and $a$ is a single number fixed by the mass and the equation of state, tabulated in Table VII and ranging from about $2.0$ (soft equation of state G at $1.4 M_*$) up to about $12.1$ (stiff equation of state L at $1.0 M_*$). Equivalently, $Q = -a c^2 J^2 / M$. The sign is negative because rotation makes the star oblate. Alongside this, and controlling for mass and spin, the magnitude of $q$ increases monotonically with the stiffness of the equation of state: $|q(G)| < |q(FPS)| < |q(C)| < |q(L)|. The authors are explicit that the quadratic law's survival to fast rotation is not something their theory anticipated; they compare against Maclaurin spheroids in Newtonian gravity, where the relation is quadratic only for slow rotation and departs strongly by $chi ~ 0.4$, and note that $q = -$chi^{2}$$ holds exactly for rotating black holes. Before any of this, the paper establishes that the parameter $q$ read off the $P_2(cos $\theta$)/$r^{3}$$ part of the metric potential really is the coordinate-invariant quadrupole moment, by computing the orbital quantity $$\Delta$ tilde{E}$ for a test particle on circular equatorial orbits and matching its $$v^{6}$$ coefficient to the general formula.
Load-bearing premise
The claim that the quadratic fit holds to one part in a thousand across the whole rotation range rests on a numerical accuracy the paper asserts rather than demonstrates: it states that the computations are accurate to at least one part in $10^3$ but reports no resolution scan, no convergence study, and no side-by-side table of the two extraction methods' results.
Editorial extensions
If this is right
- The gravitational-wave phasing of inspiraling binary neutron stars can be corrected using a one-parameter formula for each star's quadrupole moment, rather than requiring a full numerical model for every system.
- Because the effect is comparable in size to the general-relativistic spin-spin interaction, the quadrupole correction belongs in the same order of the waveform expansion and cannot be dropped while keeping spin-spin.
- A measurement of a neutron star's quadrupole moment, combined with its mass and spin, would constrain the stiffness of the equation of state through the fitted coefficient $a$.
- The tabulated values of $a$ give a compact interface between numerical relativity and analytic binary inspiral calculations: four equations of state and five masses cover a usable chunk of the plausible parameter space.
Reading between the lines
- If the quadratic law really holds to the Keplerian limit for realistic equations of state, it suggests the oblateness of a rotating star is governed mainly by how the centrifugal support redistributes mass, and only weakly by the detailed radial structure, which is why a single coefficient suffices.
- The fitted $a$ falls with increasing mass at fixed equation of state; a reader could test whether $a$ scales roughly with the star's radius squared over $M^2$, which is the naive dimensional guess the authors invoke but do not fit.
- The paper's two extraction routes are a genuine independent check only if their difference is reported; publishing that difference would turn a plausibility argument into a measurable error bar and would let outsiders judge whether the tightness of the quadratic fit is real or an artifact of the fit absorbing numerical scatter.
- Because the relation is exact with $a = 1$ for black holes, the neutron-star coefficient $a$ can be read as a single number interpolating between the Newtonian and black-hole limits, which makes it a natural target for a first-principles calculation rather than a purely numerical output.
Original abstract
Numerical models of rotating neutron stars are constructed for four equations of state using the computer code RNS written by Stergioulas. For five selected values of the star's gravitational mass (in the interval between 1.0 and 1.8 solar masses) and for each equation of state, the star's angular momentum is varied from J=0 to the Keplerian limit J=J_{max}. For each neutron-star configuration we compute Q, the quadrupole moment of the mass distribution. We show that for given values of M and J, |Q| increases with the stiffness of the equation of state. For fixed mass and equation of state, the dependence on J is well reproduced with a simple quadratic fit, Q \simeq - aJ^2/M c^2, where c is the speed of light, and a is a parameter of order unity depending on the mass and the equation of state.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
Numerical models of rotating neutron stars are constructed for four equations of state using the computer code RNS written by Stergioulas. For five selected values of the star's gravitational mass (in the interval between 1.0 and 1.8 solar masses) and for each equation of state, the star's angular momentum is varied from J=0 to the Keplerian limit J=J_{max}. For each neutron-star configuration we compute Q, the quadrupole moment of the mass distribution. We show that for given values of M and J, |Q| increases with the stiffness of the equation of state. For fixed mass and equation of state, the dependence on J is well reproduced with a simple quadratic fit, Q \simeq - aJ^2/M c^2, where c is the speed of light, and a is a parameter of order unity depending on the mass and the equation of state.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | Numerical models of rotating neutron stars are constructed for four equations of state using the computer code RNS written by Stergioulas. | method 0.99 | 0.83 | 0.02 | asserted 0.98 |
| s1 | For five selected values of the star's gravitational mass (in the interval between 1.0 and 1.8 solar masses) and for each equation of state, the star's angular momentum is varied from J=0 to the Keplerian limit J=J_{max}. | method 0.99 | 0.87 | 0.03 | asserted 0.79 |
| s2 | For each neutron-star configuration we compute Q, the quadrupole moment of the mass distribution. | method 0.94 | 0.96 | 0.05 | asserted 0.94 |
| s3 | We show that for given values of M and J, |Q| increases with the stiffness of the equation of state. | result 1.00 | 0.94 | 0.55 | asserted 1.00 |
| s4 | For fixed mass and equation of state, the dependence on J is well reproduced with a simple quadratic fit, Q \simeq - aJ^2/M c^2, where c is the speed of light, and a is a parameter of order unity depending on the mass and the equation of state. | result 0.99 | 0.93 | 0.32 | asserted 0.96 |
What was examined, and what was not
- arXiv listing abstract v2
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:gr-qc/9709033 by the authors.
References
How later papers used it
78 later papers cite it
4 later papers by other authors are listed below: 3 used its method or tool and 1 used its data or values.
Uses its method or tool 3 of 4 papers with this label
-
The physical quadrupole moment $Q$ is then obtained through [equation] where $b$ is a dimensionless coefficient determined from the asymptotic behavior of the metric potentials (see Eq. (9) in [ref] ), and $M_2$ corresponds to the coefficient of the $l=2$ (quadrupole) term in the expansion (see Section 2 of [ref] ).
Rotational effects in quark stars: comparing different models (2025) · Rotating Neutron Stars · 1.00
-
The latter is extracted from the asymptotic expansion of the $g_tt$ component, following the method described in [ref] .
Neutron stars and the cosmological constant problem (2024) · Fast Rotating Neutron Stars · 1.00
-
The physical quadrupole moment \(Q\) is then obtained as [equation] where \(b\) is a dimensionless coefficient determined from the asymptotic behavior of the metric potentials (see Eq. (9) of Ref. [ref] ), and \(M_2\) corresponds to the \(l=2\) multipole coefficient of the expansion [ref] .
Rotational enhancement and stability of protoquark stars during thermal evolution (2026) · Rotational Properties · 0.97
Uses its data or values 1 of 17 papers with this label
-
Averaged this paper's values of alpha for each equation of state for a 1.4 solar mass system, yielding a quadrupole moment.
Classical general relativity effects by magnetars with massive quadrupole, angular momentum and a magnetic dipole (2024) · 0.99 · paraphrased: the licence does not allow quoting
How the literature uses it
59 arXiv papers cite this in 78 sentences
Papers citing it · 50 sentences without a confident label
Uses its method or tool 3 papers
-
The latter is extracted from the asymptotic expansion of the $g_tt$ component, following the method described in [ref] .
Neutron stars and the cosmological constant problem · Fast Rotating Neutron Stars · 1.00
-
Within the framework, the coefficient \(M_2\) is extracted from the \(r^-3 P_2(cosθ)\) term in the expansion of the \(ρ\)-potential [ref] .
Rotational enhancement and stability of protoquark stars during thermal evolution · Rotational Properties · 1.00
-
In addition, the inclusion of quadrupole moments is due to Morsink based on the method by [ref] .
Astrophysical aspects of general relativistic mass twin stars · Full solution for uniform rotational bodies · 0.98
Uses its data or values 13 papers
-
According to Table VII. of Ref. [ref] the respective quadrupole parameters of the neutron stars with a FPS (Friedman-Pandharipande + Skyrme) equation of state are $w_NS_1=2.55$ and $w_NS_2=4.3$.
Stability analysis of the spin evolution fix points in inspiraling compact binaries with black hole, neutron star, gravastar, or boson star components · Numerical investigation for high spins · 1.00
-
For BHs, $C_ = 1$ [ref] , while for NSs, its value depends on the equation of state [ref] .
Spin-Precessing Black Hole Binaries in Dynamical Chern-Simons Gravity · Quadrupole Approximation · 0.84
-
For NSs, $λ$ takes value between $∼4-30$ [ref] , and $λ∼10-200$ for boson stars consisting of a self-interacting massive scalar field [ref] .
Constraining the nature of dark compact objects with spin-induced octupole moment measurement · Spin-induced octupole moment · 0.83
Compares with it 1 paper
-
It should be noted that an analogous result was obtained by Laarakkers and Poisson [ref] .
Gravitational field of compact objects in general relativity · The exterior Fock solution · 0.98
Cites it as background 7 papers
-
With PN theory, the spin precession equations up to 2PN order contain the spin-orbit, spin-spin, and monopole-quadrupole interactions, the latter of these being dependent on the structure of the component objects of the binary [ref] .
Parity Violation in Spin-Precessing Binaries: Gravitational Waves from the Inspiral of Black Holes in Dynamical Chern-Simons Gravity · Introduction · 0.87
-
Due to the no-hair conjecture [ref] , $κ$ is unity for BHs in GR, while it may take a range of values for neutron stars or BH mimickers [ref] .
Tests of General Relativity with Binary Black Holes from the second LIGO-Virgo Gravitational-Wave Transient Catalog · Spin-induced quadrupole moment · 0.87
-
These include a broad spectrum of oscillation modes [ref] , non-constant spin induced quadrupole moments [ref] , higher order tides beyond the quadrupole [ref] and current multipolar effects [ref] .
Investigating the impact of quasi-universal relations on neutron star constraints in third-generation detectors · Introduction · 0.86
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.