Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

Nature of phase transitions and metastability in scalar-tensor theories

T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash

Pith's one-line read Spontaneous scalarization of neutron stars usually proceeds as a first-order phase transition, not the commonly assumed second-order one.

desk verdict The core claim—first-order scalarization is the rule in the canonical DEF model—is numerically well supported; the 'norm' generalization to all EOSs is thinner than the headline. read the letter →

arxiv 2502.01781 v2 pith:UVCHZ2LK submitted 2025-02-03 gr-qc

classification gr-qc
keywords spontaneousscalarizationscalar-tensortheoryfirst-orderphasetransitionneutronstarmetastabilityLandauscalarchargestrong-fieldgravity
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Spontaneous scalarization, the sudden growth of a scalar field around neutron stars above a critical mass, is usually described as a second-order phase transition in which the field grows continuously from zero. This paper argues that, for the standard massive scalar-tensor theory action, that smooth picture holds only in a narrow part of parameter space. In most cases where scalarization occurs, the onset is a first-order phase transition: the star jumps discontinuously from an unscalarized to a scalarized configuration. As a result, there is an interval of baryon masses in which both a scalarized and an unscalarized neutron star exist as locally stable equilibria, one of them metastable. A reader should care because metastable stars and sudden transitions are far more promising observational targets than the smooth growth previously assumed.

What carries the argument

The load-bearing object is the Landau free-energy ansatz, Eq. (3): $M_{\rm ADM}=M_0(M_b)+a(M_b)Q^2+\frac{1}{2}b(M_b)Q^4+\frac{1}{3}c(M_b)Q^6$, where the scalar charge $Q$, defined by $\phi(r\to\infty)=\phi_\infty+Q/r+\cdots$, is the order parameter and only even powers appear because of the $\phi\to-\phi$ symmetry. The sign of $a(M_b)$ controls whether the unscalarized $Q=0$ solution is stable, while the sign of $b(M_b)$ at the branching mass $M_{\rm crit}$ decides the order: $b>0$ gives a continuous, second-order onset, and $b<0$ gives a first-order transition with a coexistence region, metastable states, and a discontinuous jump. To determine $b_0=b(M_{\rm crit})$, the paper fits the near-critical relation $2(M_0-M_{\rm ADM})=b_0Q^4$ to numerically constructed neutron-star solutions, increasing numerical precision to control catastrophic cancellation.

What would settle it

Compute the same neutron-star families with quadruple-precision arithmetic or an independent solver and check whether three equilibrium solutions with identical baryon mass still exist in the interval $M_{\rm bottom}<M_b<M_{\rm crit}$ for a representative case such as $\beta=-40$, $m_\phi=4\times10^{-11}$ eV. If no such coexistence interval survives, or if the fitted slope $b_0$ in the relation $2(M_0-M_{\rm ADM})=b_0Q^4$ changes sign under improved precision, the first-order classification collapses.

Watch

Extended reading notes

Core claim

The paper's central claim is that the onset of spontaneous scalarization in the scalar-tensor theory with coupling $A(\phi)=e^{\beta\phi^2/2}$ and a scalar mass $m_\phi$ is generically a first-order phase transition, not the second-order transition invoked since the earliest treatments. For a given baryon mass in the coexistence region, there are three equilibrium configurations: an unscalarized star, an unstable scalarized star with lower scalar charge, and a stable scalarized star. The energy of a configuration is its ADM mass, and the globally stable solution switches from the unscalarized to the scalarized one at a specific baryon mass, with the scalar charge jumping discontinuously from zero to a finite value. The paper supports this by computing equilibrium neutron-star families across a broad grid of $(\beta, m_\phi)$ values, extracting from them the sign of the Landau coefficient $b$ at the branching point, and showing that $b$ is negative, hence first order, for almost the entire scalarizing parameter space, including the massless $m_\phi=0$ theory.

Load-bearing premise

The classification rests on the assumption that the truncated energy expansion in the scalar charge is accurate near the branching point and that the fitted coefficient $b_0$ has the correct sign despite severe loss of precision when subtracting nearly equal masses; if that fit is wrong, the transition could actually be second order or of a shape the expansion misses.

Editorial extensions

If this is right

  • For essentially all scalarizing parameters in this theory, scalarized and unscalarized neutron stars of the same baryon mass coexist over a finite mass interval, with one configuration metastable.
  • A star that crosses the critical mass can suddenly jump into or out of a scalarized state, releasing energy comparable to the ADM mass difference between the two branches.
  • Such discontinuous transitions would emit gravitational waves carrying extra polarization modes beyond general relativity, with potentially prolonged signals and accompanying electromagnetic and neutrino signatures.
  • Observational searches should target the intermediate parameter region where both the ADM mass difference and the branch-off ADM mass are moderately large, since that region gives the most detectable events.
  • The theory with $m_\phi=0$, strongly constrained by binary-pulsar observations, already exhibits first-order scalarization, and adding a scalar mass leaves almost the whole parameter space unconstrained.

Reading between the lines

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

  • Editorial inference: if first-order scalarization is generic, the final state of a neutron star can depend on its formation and accretion history, not just on its baryon mass, because a metastable branch can persist until a sufficiently large perturbation knocks it over.
  • Editorial inference: the same Landau logic with a different order parameter should apply to curvature-induced scalarization of black holes in scalar-Gauss-Bonnet theory; the paper's explanation via higher-order coupling terms already suggests why those transitions are preferentially first order, a prediction that could be tested by extracting $b_0$ for those theories.
  • Editorial inference: a direct dynamical test would be to evolve a neutron star slowly accreting across the coexistence region and look for hysteresis, meaning scalarization happening at a different mass than descalarization.
  • Editorial inference: the tricritical point at $\beta\approx-10.6$ for $m_\phi=0$ is a sharp prediction of the paper's framework, and measuring how $b_0(\beta)$ changes sign with an independent equation of state would check whether the boundary between first- and second-order scalarization is robust.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The manuscript studies spontaneous scalarization in the scalar-tensor theory of Eq. (1), with A(phi)=exp(beta phi^2/2) and a scalar mass m_phi, by solving the static, spherically symmetric neutron-star equations with the HB piecewise-polytropic equation of state. For a domain of the (beta, m_phi) parameter space the authors find that the baryon-mass versus central-density diagram contains a region with three equilibrium solutions at fixed baryon mass, and that the fractional binding energy crosses between the scalarized and unscalarized branches. They interpret this as first-order scalarization, describe it with the Landau free-energy ansatz of Eq. (3) using the scalar charge Q as order parameter, and extract the Landau coefficient b0=b(M_crit) from fits to near-critical numerical data via End Matter Eq. (E10). The central claim is that first-order scalarization is the most common mechanism for the onset of scalarization in this theory, so that metastable configurations and discontinuous jumps are more likely than previously believed.

Significance. If the central claim holds, the paper reframes the onset of spontaneous scalarization: continuous second-order scalarization would occur only in a narrow part of the parameter space, while first-order transitions with metastable branches would dominate. This is observationally relevant because discontinuous scalarization/descalarization would produce distinctive gravitational-wave, electromagnetic, and neutrino signatures. The paper's concrete strengths are the explicit solution families in Fig. 1, the clean binding-energy crossing in Fig. 2, the parameter-space map in Fig. 3, and the economical explanation of branch ordering through Landau theory. The main caveats are that the prevalence claim is demonstrated with only one EOS, and that the sign of b0, which sets the transition order, is obtained from low-signal fits in a catastrophic-cancellation regime without reported uncertainties.

major comments (3)
  1. [Sec. II and Fig. 3; Discussion] The central claim that first-order scalarization is “the norm” for action (1) rests on a parameter scan computed with a single equation of state, the HB piecewise polytrope of Ref. [30]. The Discussion correctly states that the Landau parameters a, b, c, M0, and M_crit depend on the EOS, but the only support for EOS independence is the sentence in Sec. II that results are “qualitatively similar for other choices.” No second EOS is tested, no reference is given, and no argument is supplied. Because the sign of b(Mb; beta, m_phi, EOS) controls the transition order, the observed dominance of first-order scalarization could be an artifact of the HB EOS. Please repeat the b0(beta) extraction and at least the Delta-M map for a softer and a stiffer EOS, or else soften the prevalence claim to the HB EOS; a concrete test such as using APR, SLy, or a family of polytropes would resolve this concern.
  2. [End Matter, Eq. (E10) and Fig. 5] The classification of the transition order hinges on the sign of b0, which is obtained by fitting 2(M0 - M_ADM) = b0 Q^4 to equilibrium solutions near M_crit. This fit assumes that the quartic leading term dominates at the fitted Q values and that the near-critical points lie in the asymptotic regime. The text acknowledges the loss-of-precision problem and states that the points closest to the origin are excluded from the fit, but it does not report uncertainties on b0, fit-range sensitivity, or convergence checks with respect to numerical precision and tolerance. Since the first-to-second-order boundary occurs where b0 crosses zero near beta = -10.6, a small systematic error in the fitted slope can move the boundary substantially. Please provide error bars and robustness tests for b0(beta), for example by varying the number of fitted points, the maximum Q^4 included, and the numerical precision, and state how the boundary in Fig. 3 depends on this uncertainty.
  3. [End Matter, Eqs. (3) and (E8)–(E10)] The explanatory content of the Landau analysis is weaker than it appears. Eq. (3) is a phenomenological ansatz, not derived from action (1), and the authors state that no first-principles relation between b0 and the theory parameters is known. The coefficients are fitted to the same numerical solutions whose branch structure they are used to explain, so Fig. 4 should be read as a qualitative analogy rather than an independent confirmation of first-order behavior. This does not invalidate the direct numerical observation of three solutions and metastability, which stands on its own. A more decisive test would be to compare the predicted scaling Q_-^2 proportional to (M_crit - M_b) from Eq. (E8) with the numerical values over a range of masses, which would validate the quartic ansatz independently of the single b0 fit.
minor comments (4)
  1. [Fig. 5 and surrounding text] Fig. 5 shows b0(beta) only for m_phi = 0, while the text says that “the case being similar for all scalar masses.” Because the scalar charge Q is not well defined for massive scalars, the analogous extraction is not shown; please clarify how this statement is supported, or mark it as an expectation rather than a computed result.
  2. [Fig. 3 caption] The threshold Delta M < 10^-5 M_sun, below which first-order scalarization is indistinguishable from numerical noise, determines the location of the white regions in Fig. 3. A brief statement of the numerical resolution at that scale, or a convergence test showing that the boundary is stable, would make the prevalence claim easier to assess.
  3. [Sec. II, Fig. 1] The text refers to “the dotted curve” and “the solid red” in Fig. 1, but the figure legend labels the scalarized branches by color and the unstable segment is not explicitly identified in the legend. Adding a legend entry such as “unstable scalarized branch (dotted)” would improve readability.
  4. [End Matter, Eq. (E8)] Equation (E8) is written with a sign that gives Q^2 negative for M_b > M_crit; the text correctly explains that the relevant branch is M_b < M_crit, but the sign convention could be stated explicitly for clarity.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the first-order prevalence claim is read directly from computed star sequences; the Landau fit is a post-hoc phenomenological description, not a prediction.

full rationale

The paper's central claim, that first-order scalarization dominates the parameter space of action (1), is established by directly solving the TOV equations and reading the structure of the Mb(rho_c) and binding-energy curves (Figs. 1-3 and E1), not by the Landau free-energy fit. The Landau ansatz Eq. (3) is introduced as a later phenomenological re-description, and its parameters a, b, c, and M0 are explicitly extracted from the same numerical solutions via Eq. (E10); the paper even acknowledges that "there is no understanding for the behavior of b0 in Fig. 5 from first principles." Thus the fit is not used to predict a quantity independent of its input; it is an interpretation of already-computed equilibria. The first/second-order classification in Fig. 3 is based on the direct observable Delta M being nonzero, i.e. on the presence of two stable solutions with the same baryon mass, not on the fitted b0. Self-citations to Tuna et al. [29] and Ramazanoğlu & Pretorius [24] provide the numerical method and the massive-scalar model, respectively; these are method/model citations and the present solutions are computed, not imported as the conclusion. The unverified statement that results are "qualitatively similar for other choices" of EOS is a missing-support caveat, not a circular step. Overall, the derivation chain from the field equations to the numerical solutions to the phase-transition language is self-contained, with the Landau model serving as an explanatory overlay rather than as the evidential basis for the headline claim.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The paper's central claim is a numerical survey plus a phenomenological interpretation. The Landau coefficients are fitted to the same numerical solutions they are used to explain, so the report lists them as free parameters. No new particles or forces are introduced, and the model action is taken as input from prior literature. The most consequential assumed inputs are the HB equation of state, the turning-point stability criterion, and the validity of the Landau expansion.

free parameters (3)
  • b0 = b(Mcrit) = Functional of beta and mphi; zero near beta = -10.6 for mphi = 0
    Determines the order of scalarization; extracted by fitting 2(M0-MADM)=b0 Q^4 to numerical solutions near the branch-off point (Eq. E10, Fig. E3).
  • Landau coefficients a(Mb) and c(Mb) = Not reported explicitly
    The ansatz Eq. (3) uses a, b, and c; a1 is used to relate Q^2 to (Mcrit-Mb), and c(Mb)>0 is assumed. Only b0 is quantified in the paper.
  • Numerical threshold for Delta M = 1e-5 solar masses
    Used to distinguish first-order scalarization from numerical noise in Fig. 3; points below the threshold are excluded, which affects the prevalence map.
assumptions (6)
  • domain assumption The scalar-tensor action (1) with A(phi)=exp(beta phi^2/2) and scalar mass mphi is the model under study.
    Defines the theoretical framework; not derived in this paper.
  • domain assumption Static, spherically symmetric TOV solutions obtained by the relaxation method of Tuna et al. represent all relevant equilibrium neutron star configurations.
    Invoked throughout; no dynamical or non-spherical perturbations are considered.
  • standard math The turning-point criterion: along a one-parameter family of equilibria, stability changes at dMb/drho_c = 0, and the dMb/drho_c < 0 section is unstable.
    Used to label branches in Fig. 1 and Fig. 2; cites Sorkin and Shapiro-Teukolsky.
  • ad hoc to paper Landau expansion Eq. (3) with scalar charge Q as order parameter and c(Mb)>0 is a valid free energy ansatz.
    Phenomenological; coefficients are fitted, not derived from the action.
  • domain assumption The unscalarized solution is locally stable for Mb < Mcrit because no tachyonic scalar mode exists and the GR branch has no hydrodynamical instability.
    Needed for the metastability region; asserted in the main text with references.
  • domain assumption The HB equation of state from Read et al. is representative for neutron stars.
    Used for all plots; the paper claims qualitative similarity for other EOS but does not show it.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Nature of phase transitions and metastability in scalar-tensor theories." pith.science (2026). https://pith.science/paper/UVCHZ2LK

@misc{pith2026250201781,
  author       = {Pith},
  title        = {Pith review of: Nature of phase transitions and metastability in scalar-tensor theories},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UVCHZ2LK}},
  note         = {Machine review of arXiv:2502.01781}
}
read the original abstract

Compact stars above a critical stellar mass develop large scalar fields in some scalar-tensor theories. This scenario, called spontaneous scalarization, has been an intense topic of study since it passes weak-field gravity tests naturally while providing clear observables in the strong-field regime. The underlying mechanism for the onset of scalarization is often depicted as a second-order phase transition. Here, we show that a first-order phase transition is in fact the most common mechanism. This means metastability and transitions between locally stable compact object configurations are much more likely than previously believed, opening vast new avenues for observational prospects.

Figures

Figures reproduced from arXiv: 2502.01781 by the authors.

Figure 2
Figure 2. FIG. 2. Fractional binding energy ( [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Dependence of scalarization characteristics on the [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. FIG. 4. The Landau ansatz in Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Underlying mechanisms of phase transitions in scalar-tensor theories

    gr-qc 2026-04 unverdicted novelty 8.0 of 10

    Landau coefficients for scalarization phase transitions are calculated from first principles via reduction of the theory's energy functional to an effective energy function.

Reference graph

Works this paper leans on

48 extracted references · 12 canonical work pages · cited by 1 Pith paper

  1. [30]

    S. Tuna, K. I. ¨Unl¨ ut¨ urk, and F. M. Ramazano˘ glu, Con- straining scalar-tensor theories using neutron star mass and radius measurements, Phys. Rev. D 105, 124070 (2022), arXiv:2204.02138 [gr-qc]

  2. [1]

    B. P. Abbott et al. (LIGO Scientific, Virgo), Tests of general relativity with GW150914, Phys. Rev. Lett. 116, 221101 (2016), [Erratum: Phys.Rev.Lett. 121, 129902 (2018)], arXiv:1602.03841 [gr-qc]

  3. [2]

    B. P. Abbott et al. (LIGO Scientific, Virgo), Tests of Gen- eral Relativity with GW170817, Phys. Rev. Lett. 123, 011102 (2019), arXiv:1811.00364 [gr-qc]

  4. [3]

    Abbott et al

    R. Abbott et al. (KAGRA, VIRGO, LIGO Scientific), Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3, Phys. Rev. X13, 011048 (2023), arXiv:2111.03634 [astro-ph.HE]

  5. [4]

    Abbott et al

    R. Abbott et al. (KAGRA, VIRGO, LIGO Scien- tific), GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run, Phys. Rev. X 13, 041039 (2023), arXiv:2111.03606 [gr-qc]

  6. [5]

    K. G. Arun et al. (LISA), New horizons for fundamen- tal physics with LISA, Living Rev. Rel. 25, 4 (2022), arXiv:2205.01597 [gr-qc]

  7. [6]

    Barack et al

    L. Barack et al. , Black holes, gravitational waves and fundamental physics: a roadmap, Class. Quant. Grav. 36, 143001 (2019), arXiv:1806.05195 [gr-qc]

  8. [7]

    Gendreau and Z

    K. Gendreau and Z. Arzoumanian, Searching for a pulse, Nat. Astron. 1, 895 (2017)

Show all 48 references
  1. [8]

    Damour and G

    T. Damour and G. Esposito-Far` ese, Nonperturbative strong-field effects in tensor-scalar theories of gravitation, Phys. Rev. Lett. 70, 2220 (1993)

  2. [9]

    D. D. Doneva, F. M. Ramazano˘ glu, H. O. Silva, T. P. Sotiriou, and S. S. Yazadjiev, Spontaneous scalarization, 6 Rev. Mod. Phys. 96, 015004 (2024), arXiv:2211.01766 [gr-qc]

  3. [10]

    Damour and G

    T. Damour and G. Esposito-Farese, Tensor-scalar gravity and binary-pulsar experiments, Phys. Rev. D 54, 1474 (1996), arXiv:gr-qc/9602056

  4. [11]

    Sperhake, C

    U. Sperhake, C. J. Moore, R. Rosca, M. Agathos, D. Gerosa, and C. D. Ott, Long-lived inverse chirp sig- nals from core collapse in massive scalar-tensor gravity, Phys. Rev. Lett. 119, 201103 (2017), arXiv:1708.03651 [gr-qc]

  5. [12]

    D. D. Doneva and S. S. Yazadjiev, New Gauss-Bonnet Black Holes with Curvature-Induced Scalarization in Ex- tended Scalar-Tensor Theories, Phys. Rev. Lett. 120, 131103 (2018), arXiv:1711.01187 [gr-qc]

  6. [13]

    H. O. Silva, J. Sakstein, L. Gualtieri, T. P. Sotiriou, and E. Berti, Spontaneous scalarization of black holes and compact stars from a Gauss-Bonnet coupling, Phys. Rev. Lett. 120, 131104 (2018), arXiv:1711.02080 [gr-qc]

  7. [14]

    R. F. P. Mendes and N. Ortiz, New class of quasinormal modes of neutron stars in scalar-tensor gravity, Phys. Rev. Lett. 120, 201104 (2018), arXiv:1802.07847 [gr-qc]

  8. [15]

    C. A. R. Herdeiro, E. Radu, N. Sanchis-Gual, and J. A. Font, Spontaneous Scalarization of Charged Black Holes, Phys. Rev. Lett. 121, 101102 (2018), arXiv:1806.05190 [gr-qc]

  9. [16]

    C. A. R. Herdeiro, E. Radu, H. O. Silva, T. P. Sotiriou, and N. Yunes, Spin-induced scalarized black holes, Phys. Rev. Lett. 126, 011103 (2021), arXiv:2009.03904 [gr-qc]

  10. [17]

    Berti, L

    E. Berti, L. G. Collodel, B. Kleihaus, and J. Kunz, Spin-induced black-hole scalarization in Einstein-scalar- Gauss-Bonnet theory, Phys. Rev. Lett. 126, 011104 (2021), arXiv:2009.03905 [gr-qc]

  11. [18]

    H.-J. Kuan, A. G. Suvorov, D. D. Doneva, and S. S. Yazadjiev, Gravitational Waves from Accretion-Induced Descalarization in Massive Scalar-Tensor Theory, Phys. Rev. Lett. 129, 121104 (2022), arXiv:2203.03672 [gr-qc]

  12. [19]

    D. D. Doneva and S. S. Yazadjiev, Beyond the sponta- neous scalarization: New fully nonlinear mechanism for the formation of scalarized black holes and its dynam- ical development, Phys. Rev. D 105, L041502 (2022), arXiv:2107.01738 [gr-qc]

  13. [20]

    D. D. Doneva, L. G. Collodel, and S. S. Yazadjiev, Spon- taneous nonlinear scalarization of Kerr black holes, Phys. Rev. D 106, 104027 (2022), arXiv:2208.02077 [gr-qc]

  14. [21]

    D. D. Doneva, C. J. Kr¨ uger, K. V. Staykov, and P. Y. Yordanov, Neutron stars in Gauss-Bonnet gravity: Non- linear scalarization and gravitational phase transitions, Phys. Rev. D 108, 044054 (2023), arXiv:2306.16988 [gr- qc]

  15. [22]

    A. M. Pombo and D. D. Doneva, Effects of mass and self- interaction on nonlinear scalarization of scalar-Gauss- Bonnet black holes, Phys. Rev. D 108, 124068 (2023), arXiv:2310.08638 [gr-qc]

  16. [23]

    K. V. Staykov and D. D. Doneva, Nonlinear black hole scalarization in multi-scalar Gauss-Bonnet gravity, J. Phys. Conf. Ser. 2719, 012007 (2024)

  17. [24]

    F. M. Ramazano˘ glu and F. Pretorius, Spontaneous Scalarization with Massive Fields, Phys. Rev. D93, 064005 (2016), arXiv:1601.07475 [gr-qc]

  18. [25]

    Alsing, E

    J. Alsing, E. Berti, C. M. Will, and H. Zaglauer, Grav- itational radiation from compact binary systems in the massive Brans-Dicke theory of gravity, Phys. Rev. D 85, 064041 (2012), arXiv:1112.4903 [gr-qc]

  19. [26]

    Antoniadis, P

    J. Antoniadis, P. C. C. Freire, N. Wex, T. M. Tauris, R. S. Lynch, M. H. van Kerkwijk, M. Kramer, C. Bassa, V. S. Dhillon, T. Driebe, J. W. T. Hessels, V. M. Kaspi, V. I. Kondratiev, N. Langer, T. R. Marsh, M. A. McLaugh- lin, T. T. Pennucci, S. M. Ransom, I. H. Stairs, J. van...

  20. [27]

    Kramer et al

    M. Kramer et al. , Strong-Field Gravity Tests with the Double Pulsar, Phys. Rev. X 11, 041050 (2021), arXiv:2112.06795 [astro-ph.HE]

  21. [28]

    J. Zhao, P. C. C. Freire, M. Kramer, L. Shao, and N. Wex, Closing a spontaneous-scalarization window with binary pulsars, Class. Quant. Grav. 39, 11LT01 (2022), arXiv:2201.03771 [astro-ph.HE]

  22. [29]

    End Matter

    to obtain the static, spherically symmetric NS solu- tions for the action (1), also see the “End Matter.” Both scalarized and unscalarized solutions are plotted in Fig. 1 as one-parameter families in terms of the central density of the star ˜ρc. We use the HB equation of state...

  23. [31]

    J. S. Read, C. Markakis, M. Shibata, K. Uryu, J. D. E. Creighton, and J. L. Friedman, Measuring the neutron star equation of state with gravitational wave observa- tions, Phys. Rev. D 79, 124033 (2009), arXiv:0901.3258 [gr-qc]

  24. [32]

    We will be talking about points on the scalarized branches of the Mb(˜ρc) curve as if they are single solu- tions to avoid confusion, since they look as such on the Mb(˜ρc) curves

    There is a ϕ → −ϕ symmetry for our choice of A(ϕ), which means each scalarized solution in the figure rep- resents two configurations with opposite scalar field par- ity. We will be talking about points on the scalarized branches of the Mb(˜ρc) curve as if they are single solu...

  25. [33]

    Sorkin, A Criterion for the onset of instability at a turning point, Astrophys

    R. Sorkin, A Criterion for the onset of instability at a turning point, Astrophys. J. 249, 254 (1981)

  26. [34]

    R. D. Sorkin, A Stability criterion for many-parameter equilibrium families, Astrophys. J. 257, 847 (1982)

  27. [35]

    S. L. Shapiro and S. A. Teukolsky, Black Holes, White Dwarfs, and Neutron Stars: The Physics of Compact Ob- jects (WILEY-VCH, 1983)

  28. [36]

    Harada, Neutron stars in scalar tensor theories of grav- ity and catastrophe theory, Phys

    T. Harada, Neutron stars in scalar tensor theories of grav- ity and catastrophe theory, Phys. Rev. D57, 4802 (1998), arXiv:gr-qc/9801049

  29. [37]

    Doroshenko, V

    V. Doroshenko, V. Suleimanov, G. P¨ uhlhofer, and A. Santangelo, A strangely light neutron star within a supernova remnant, Nature Astron. 6, 1444 (2022)

  30. [38]

    Khalil, N

    M. Khalil, N. Sennett, J. Steinhoff, and A. Buonanno, Theory-agnostic framework for dynamical scalarization of compact binaries, Phys. Rev. D 100, 124013 (2019), arXiv:1906.08161 [gr-qc]

  31. [39]

    The end, rather than the beginning, of scalarization was called a first- order phase transition in Kuan et al

    Muniz and Mendes [47] is an exception, which was posted after the submission of our manuscript. The end, rather than the beginning, of scalarization was called a first- order phase transition in Kuan et al. [18], though this study was not specifically concerned with a phenomen...

  32. [40]

    Plischke and B

    M. Plischke and B. Bergersen, Equilibrium Statistical Physics (World Scientific, Singapore, 2006)

  33. [41]

    Goldenfeld, Lectures on Phase Transitions and the Renormalization Group (CRC Press, Boca Raton, 2018)

    N. Goldenfeld, Lectures on Phase Transitions and the Renormalization Group (CRC Press, Boca Raton, 2018)

  34. [42]

    However, in general, ϕ∞ has an effect on scalarization analogous to the effect of an external magnetic field on magnetization

    We are only interested in solutions with ϕ∞ = 0, it therefore plays no role in our phase transition picture. However, in general, ϕ∞ has an effect on scalarization analogous to the effect of an external magnetic field on magnetization

  35. [43]

    Andreou, N

    N. Andreou, N. Franchini, G. Ventagli, and T. P. 7 Sotiriou, Spontaneous scalarization in generalized scalar- tensor theory, Phys. Rev. D99, 124022 (2019), arXiv:1904.06365 [gr-qc]

  36. [44]

    E. S. Demirbo˘ ga, Y. E. S ¸ahin, and F. M. Ra- mazano˘ glu, Subtleties in constraining gravity theories with mass-radius data, Phys. Rev. D 108, 024028 (2023), arXiv:2303.01910 [gr-qc]

  37. [45]

    B. D. Metzger, L. Hui, and M. Cantiello, Fragmentation in Gravitationally Unstable Collapsar Disks and Subso- lar Neutron Star Mergers, Astrophys. J. Lett. 971, L34 (2024), arXiv:2407.07955 [astro-ph.HE]

  38. [46]

    B. C. Stephens, W. E. East, and F. Pretorius, Eccen- tric Black Hole-Neutron Star Mergers, Astrophys. J. Lett. 737, L5 (2011), arXiv:1105.3175 [astro-ph.HE]

  39. [47]

    W. E. East, F. Pretorius, and B. C. Stephens, Eccentric black hole-neutron star mergers: effects of black hole spin and equation of state, Phys. Rev. D 85, 124009 (2012), arXiv:1111.3055 [astro-ph.HE]

  40. [48]

    J. a. V. M. Muniz and R. F. P. Mendes, On the phase transition mechanism of spontaneous scalarization (2025), arXiv:2503.11385 [gr-qc]. End Matter Computing the equilibrium NS solutions: A static, spherically-symmetric metric can be put in the form gµνdxµdxν = −eν(r)dt2 + dr2 ...

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.