Pith. Sign in

REVIEW

Neutron stars in scalar-tensor theories of gravity and catastrophe theory

General Relativity and Quantum Cosmology · Submitted 14 Jan 1998 · 29 references

Pith's verdict: A careful, correctly executed catastrophe-theory classification of scalarized neutron star equilibria; the cusp claim is sound in structure but rests on numerical envelopes shown without resolution or convergence data.

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 →

· 3 major objections · 7 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDgr-qc/9801049 (version 2)
Pith IDWYZPVR34, Pith's permanent record of this paper and its review
CitationTomohiro Harada. "Neutron stars in scalar-tensor theories of gravity and catastrophe theory." arXiv:gr-qc/9801049 (1998). https://pith.science/paper/WYZPVR34
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

Neutron star equilibrium sequences in scalar-tensor gravity undergo fold catastrophes, and for sufficiently negative coupling strength a cusp catastrophe at zero asymptotic scalar field, which is the same structure as spontaneous scalariz

The reading

The paper asks what happens to the set of stable neutron-star solutions in scalar-tensor gravity when two knobs are turned independently: the star's total baryon number $N$ and the value the scalar field takes far away, $phi_0$. Stability is settled with the turning point method, which says a solution's stability can flip only where the map from the star's central parameters to those two knobs becomes degenerate, which is exactly the bifurcation set of a catastrophe. For the quadratic coupling $A(phi)=exp(frac{1}{2}beta phi^2)$ the answer depends on $beta$: for $beta gtrsim -4.35$ only fold catastrophes occur, while for $beta lesssim -4.35$ a cusp catastrophe appears at $phi_0=0$. The payoff is a classification: spontaneous scalarization of neutron stars is not a numerical curiosity of one tuned coupling function but a structurally stable cusp, so it survives small perturbations of the theory.

What carries the argument

The carrying device is the catastrophe map $chi_m$ from the equilibrium space $M_m={(phi_0,N,omega): partial m/partialomega=0}$ down to the control plane $(phi_0,N)$. Its singular points, where $partial m/partialomega=partial^2m/partialomega^2=0$, project onto a bifurcation set that coincides with the envelope of the numerically drawn curves of constant central density $tilde{rho}_c$. The turning point method (criteria (i) and (ii)) then fixes which branch beyond a crossing is unstable, and Thom's theorem restricts structurally stable catastrophes with two controls to the fold and the cusp, which is what licenses the classification. A Landau-style expansion of the ADM energy in the scalar charge supplies the local normal forms.

Extended reading notes

Core claim

Tracking secularly stable neutron-star equilibria in the two-dimensional control space $(phi_0,N)$, the paper shows that their number and identity reorganize discontinuously as the controls vary. For $beta simeq -4$ the picture is a single fold: one stable solution exists on one side of an envelope in the $(phi_0,N)$ plane and none on the other, and at $phi_0=0$ the surviving branch is exactly the general-relativistic one. For $beta=-6$ the envelope acquires a cusp at $(phi_0,N)=(0,N_a)$: below the critical baryon number the unique stable solution has $omega=0$ and is general-relativistic, while above it the general-relativistic branch becomes unstable and two new stable branches with $omega=pm(D/C)^{1/2}(N-N_a)^{1/2}$ appear, degenerate at $phi_0=0$ and splitting for $phi_0 neq 0$. The local energy $m=frac{C}{4}omega^4-frac{D}{2}(N-N_a)omega^2-phi_0omega+m_a$ yields the cubic equation $partial m/partialomega=0$ whose one-versus-three root structure is the cusp, and gives a scalar susceptibility $chi_phi$ diverging as $|N-N_a|^{-1}$ with a factor-of-two asymmetry between the two sides of the transition.

Load-bearing premise

The cusp identification at the critical point depends on the local energy expansion having no cubic term in the scalar charge, and the paper rules that term out from the visual shape of the numerically computed bifurcation set rather than from a derivation; the passage in question even states that a cubic term would change the number of roots at the critical point, which is the opposite of what the surrounding argument needs, so the symmetric envelope in the figure is doing load-bearing work.

Editorial extensions

If this is right

  • For $beta gtrsim -4.35$, a neutron star that accretes past the critical baryon number has no stable equilibrium left and must collapse, with the critical mass shifting as $phi_0$ changes.
  • For $beta lesssim -4.35$ at $phi_0=0$, the general-relativistic branch turns unstable at $N_a$ and the star settles onto one of two mirror branches with opposite scalar charge; the mass changes continuously but its slope in $N$ does not.
  • Because fold and cusp are the only structurally stable catastrophes with two controls, the qualitative picture is expected to survive small higher-order corrections to the exponent of the coupling function.
  • The scalar charge and scalar susceptibility near the critical point obey fixed power laws, $omega propto (N-N_a)^{1/2}$ and $chi_phi propto |N-N_a|^{-1}$, giving a template for reading off critical behaviour from numerical scans of other couplings.
  • If $phi_0$ drifts cosmologically or $N$ grows by accretion, the star can jump between branches and radiate scalar gravitational waves, making branch switching a candidate source for scalar-wave observations.

Reading between the lines

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

  • The fold/cusp dichotomy should apply to any stellar family with two independent controls, such as rotating or charged configurations, so it can be used as a cheap consistency check on numerical stability scans rather than as a fact about neutron stars alone.
  • Because the cusp is generic rather than tuned, observational bounds on $beta$ from binary-pulsar timing are bounds on a structurally stable phenomenon, which makes them harder to evade by small changes to the coupling function.
  • The factor-of-two asymmetry in $chi_phi$ across $N_a$ is a sharp and inexpensive numerical target: measuring it directly would confirm the quartic Landau form without needing to resolve the whole bifurcation set.
  • If for some coupling a cubic term in $omega$ is genuinely present, the transition degrades to a fold and the jump between branches becomes discontinuous rather than continuous; mapping which couplings do this would locate where scalarization can be strong.
Original abstract

We investigate neutron stars in scalar-tensor theories. We examine their secular stability against spherically symmetric perturbations by use of a turning point method. For some choices of the coupling function contained in the theories, the number of the stable equilibrium solutions changes and the realized equilibrium may change discontinuously as the asymptotic value of the scalar field or total baryon number is changed continuously. The behaviour of the stable equilibrium solutions is explained by fold and cusp catastrophes. Whether the cusp catastrophe appears or not depends on the choices of the coupling function. These types of the catastrophes are structurally stable. Recently discovered spontaneous scalarization, which is non-perturbative strong-field phenomenon due to the presence of the gravitational scalar field, is well described in terms of the cusp catastrophe.

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, 5,463 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.15 · human participants 0.01 · controlled comparison 0.03 · randomized assignment 0.01 · states quantitative result 0.05 · compares with prior results 0.21 · mentions code or data availability 0.02.

The abstract, sentence by sentence

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

We investigate neutron stars in scalar-tensor theories. We examine their secular stability against spherically symmetric perturbations by use of a turning point method. For some choices of the coupling function contained in the theories, the number of the stable equilibrium solutions changes and the realized equilibrium may change discontinuously as the asymptotic value of the scalar field or total baryon number is changed continuously. The behaviour of the stable equilibrium solutions is explained by fold and cusp catastrophes. Whether the cusp catastrophe appears or not depends on the choices of the coupling function. These types of the catastrophes are structurally stable. Recently discovered spontaneous scalarization, which is non-perturbative strong-field phenomenon due to the presence of the gravitational scalar field, is well described in terms of the cusp catastrophe.

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 We investigate neutron stars in scalar-tensor theories. contribution 0.89 0.87 0.03 asserted 0.73
s1 We examine their secular stability against spherically symmetric perturbations by use of a turning point method. method 0.97 0.93 0.03 asserted 0.96
s2 For some choices of the coupling function contained in the theories, the number of the stable equilibrium solutions changes and the realized equilibrium may change discontinuously as the asymptotic value of the scalar field or total baryon number is changed continuously. result 1.00 0.82 0.53 hedged possibility 0.58
s3 The behaviour of the stable equilibrium solutions is explained by fold and cusp catastrophes. interpretation 0.98 0.75 0.35 asserted 1.00
s4 Whether the cusp catastrophe appears or not depends on the choices of the coupling function. result 0.91 0.74 0.13 asserted 0.99
s5 These types of the catastrophes are structurally stable. interpretation 0.53 0.50 0.15 asserted 1.00
s6 Recently discovered spontaneous scalarization, which is non-perturbative strong-field phenomenon due to the presence of the gravitational scalar field, is well described in terms of the cusp catastrophe. interpretation 0.62 0.25 0.25 asserted 1.00

What kind of paper the classifier reports this is: theoretical result 0.95 · evidence from proof or derivation 0.52.

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

References
How later papers used it

57 later papers cite it

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

Pith read 8 of the 90 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 1 of 2 papers with this label
Cites it as support 1 of 2 papers with this label
How the literature uses it

47 arXiv papers cite this in 75 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 · 53 sentences without a confident label

Uses its method or tool 1 paper
  • Finally, we also ensure that our neutron star solutions belong to a stable branch following the prescription described in [ref] , i.e. we consider solutions where the following stability criterion is satisfied: [equation]

    Angular momentum sensitivities in scalar-tensor theories · Results · 0.97

Uses its data or values 1 paper
Compares with it 2 papers
  • This is in stark contrast to the standard case of scalarized stars in the DEF model and other scalarization theories inspired by it, where the tachyonic instability is eventually quenched, leading to stability [ref] .

    Instability of vectorized stars · Conclusions · 0.97

  • [ref] reached similar conclusions, but worked in the context of catastrophe theory.

    Spontaneous scalarization · Linearized dynamics · 1.00

Cites it as support 2 papers
  • Moreover, using the catastrophe theory, it is shown that the onset of spontaneous scalarization in the quadratic coupling function corresponds to cusp catastrophe [ref] which is structurally stable, i.e. the theory is stable against adding higher order terms in the coupling function.

    Spontaneous Scalarization in Scalar-Tensor Theories with Conformal Symmetry as an Attractor · Choice of Coupling Function · 0.98

  • Detailed investigations of stellar structure [ref] , numerical simulations of collapse [ref] and stability studies [ref] confirmed that spontaneously scalarized configurations would indeed be the end-state of stellar collapse in these theories.

    Numerical Relativity and High Energy Physics: Recent Developments · Spontaneous scalarization in compact stars · 0.97

Cites it as background 14 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.