Pith. Sign in

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

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 11, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 8 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDgr-qc/9709033 (version 2)
Pith IDA6F2WLBH, Pith's permanent record of this paper and its review
CitationEric Poisson, William G. Laarakkers. "Quadrupole moments of rotating neutron stars." arXiv:gr-qc/9709033 (1997). https://pith.science/paper/A6F2WLBH
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

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

The paper asks a concrete question: when a neutron star spins, how oblate does its mass distribution become, and how does that oblateness depend on the star's mass, its spin, and the stiffness of the matter inside it? The authors answer by building numerical models of uniformly rotating neutron stars with the publicly available code rns, covering four equations of state (labeled G, FPS, C, and L, in order of increasing stiffness), five masses from $1.0$ to $1.8 M_*$, and spins running from zero all the way to the Keplerian limit where the star sheds mass at the equator. For each model they extract the quadrupole moment $Q$ two ways and report the dimensionless version $q = c^4 Q / G^2 M^3$. Two results come out. At fixed mass and spin, $|q|$ grows with the stiffness of the equation of state, as expected because stiffer matter makes a bigger star and the quadrupole scales with the square of the radius. And at fixed mass and equation of state, the whole run from slow to maximally fast rotation collapses onto $q = -a(M, EOS) chi^2$, with $a$ a fitted number of order unity that shrinks as mass grows. The second result is the one that would matter: a quadratic spin dependence is what slow-rotation theory predicts, but it is not expected to survive to the Keplerian limit, and here it apparently does.

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 4,371 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.02 · states limitation 0.10 · human participants 0.01 · controlled comparison 0.06 · randomized assignment 0.01 · states quantitative result 0.58 · compares with prior results 0.08 · mentions code or data availability 0.04.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

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

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
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 kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1. Plots of
Figure 2. FIG. 2. Same as Fig. 1 but with
Figure 3. FIG. 3. Same as Fig. 1 but with
Figure 4. FIG. 4. Same as Fig. 1 but with
Figure 5. FIG. 5. Same as Fig. 1 but with

Figures reproduced from arXiv:gr-qc/9709033 by the authors.

References
How later papers used it

78 later papers cite it

Evidence: substantial. 78 of the citing papers share no author with it, and 22 of those carry a confident label from Pith's citation classifier.

Pith read 26 of the 103 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

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
How the literature uses it

59 arXiv papers cite this in 78 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 50 sentences without a confident label

Uses its method or tool 3 papers
Uses its data or values 13 papers
Compares with it 1 paper
Cites it as background 7 papers

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.