REVIEW 4 major objections 4 minor 1 cited by
Low-Mass Neutron Stars and Effective Phase Transitions from a Hybrid Van der Waals-Polytropic Equation of State
T0 review · 4 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read A two-branch analytic equation of state can produce neutron stars as light as 0.99 solar masses.
desk verdict A load-bearing singularity in the core EoS at the matching density invalidates the TOV results as written. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the two-branch pressure–energy-density relation known as the piecewise hybrid van der Waals–polytropic EoS. In the core it combines an excluded-volume-type repulsive term with a power-law interaction term; in the crust it is a two-term polytrope. The two core exponents τ1 and σ1 are the control knobs: they fix the number and location of critical densities (Eq. 14) via a discriminant, and they determine the effective conductivity A(η) in the generalized Lane–Emden equation that governs the smooth transition layer. A(η) multiplies the highest derivative in the interface equation, so its sign and magnitude decide mechanical stability and the width of the core–crust boundar
What would settle it
Integrate the TOV equations for each parameter set in Tables III–V while recording dp/dε and v_s^2 at every radial step; a single negative pressure slope or a sound speed above the speed of light at any density would invalidate that configuration. Alternatively, find a genuinely observed neutron star below 0.99 M_sun with well-constrained mass and radius that no parameter set in the model can reproduce.
Extended reading notes
Core claim
The paper's central claim is that a piecewise EoS built from a modified van der Waals core, p(ε)=α1/(ε0/ε−1)(ε/ε0)^τ1 + β1(ε/ε0)^σ1 for ε>ε0, and a two-term polytropic crust can generate the full range of phase-transition-like richness: monotonic EoS, weak curvature changes, and spinodal-like regions with two positive critical densities where p(ε) has two inflection points. Within the physically allowed regions of the (τ1,σ1) plane, TOV integration yields stable neutron stars spanning 0.99–2.05 M_sun. The paper finds no plateau in μ(ε), consistent with weak or smooth transitions, and the strongest curvature signatures occur at densities significantly above the matching point ε0, showing that
Load-bearing premise
The results rest on the assumption that the listed parameter sets keep dp/dε≥0 and sound speed ≤1 throughout the whole star, and that the approximate matching condition with δ≪1 faithfully represents the core–crust interface; the paper states this but does not separately demonstrate the full audit.
Editorial extensions
If this is right
- If correct, the model shows that stable neutron stars near 1.0–1.1 M_sun are structurally allowed, so observational scarcity of such stars must come from formation physics or microphysics, not hydrostatic equilibrium alone.
- Weak phase-transition signatures in μ(ε) occur deep in the core, implying that matching observations of smooth transitions to EoS features should target core curvature, not the assumed crust–core density.
- The two-parameter map gives a fast screening tool for analytic EoS before detailed nuclear calculations are attempted.
- The two-critical-density regime produces an S-shaped p(ε) resembling spinodal instability without an actual first-order transition, serving as a caution against over-interpreting curvature features as phase coexistence.
- Parameter sets that reproduce low-mass candidates can be cross-checked against X-ray timing radius measurements and gravitational-wave tidal-deformability constraints, narrowing the viable (τ1,σ1) regions.
Reading between the lines
- A direct numerical audit of every tabulated parameter set, checking dp/dε≥0 and v_s^2≤1 at all densities from center to surface, would confirm whether the claimed mass–radius families are all causally consistent; the paper asserts but does not display this check.
- The same piecewise construction could be extended to include temperature-dependence or multiple components (hyperons, quarks); whether the curvature-based transition signatures persist would test how generic the effect is.
- If future surveys find no neutron stars below ~1.2 M_sun, the parameter choices that produce the 0.99 M_sun configurations would be strongly disfavored, giving a concrete observational falsifier.
- The model's interface width scaling ℓ²∼A(η)/(4πη) suggests a testable prediction: the smoothness of the density/composition gradient in real neutron stars encodes the local EoS curvature, potentially constraining A(η) via oscillation modes.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a piecewise, analytically tractable equation of state for neutron-star matter: a modified van der Waals term for the core (ε>ε0) matched approximately to a two-term polytropic crust (ε≤ε0). It derives critical densities, a non-relativistic generalized Lane–Emden description, and then integrates the TOV equations for selected parameter sets in the (τ1,σ1) plane. The main claims are that this simple EoS can produce stable low-mass neutron stars with masses in the range (0.99–2.05) M⊙ and that the chemical potential exhibits phase-transition-like curvature changes at densities well above the matching point.
Significance. If correct, the paper would provide a useful pedagogical and heuristic testbed: a closed-form EoS that generates families of TOV solutions, including low-mass configurations, with explicit parameter sets that could be audited and reused. The authors are transparent that the model is not microphysical and explicitly note that the observed low-mass population is sparse. The paper also has the virtue of making its parameter choices explicit. However, the central construction has a singular matching-point defect and the tabulated configurations are not demonstrably solutions of the stated EoS. As submitted, the claimed phenomenology is not supported by a well-defined stellar model.
major comments (4)
- [Eqs. (8), (12); Section IV.B] The core branch (8) has a pole at ε=ε0: for ε→ε0+, p_core ≈ −α1 ε0/(ε−ε0), so it diverges for every parameter set in Tables III–V. Condition (12) avoids the pole by evaluating the core only at ε+ with ε0/ε+ ≪ 1, i.e., at ε+ far above ε0. Thus no pressure is defined for the interval (ε0, ε+) and, because the TOV equations require p(ε) over the whole density range from center to surface, the configurations in Tables III–V cannot be solutions of the stated piecewise EoS. The later statement in Section IV.B that η is continuous across r1 and that A(η)>0 describes a smooth transition layer is therefore not derived from the EoS; it assumes an unstated regularization of the singular branch.
- [Eqs. (10)–(11), near ε0] Independently of the matching problem, the core branch violates the paper's own thermodynamic condition (1) in any neighborhood of ε0. Writing x=ε/ε0−1, one obtains p≈−α1(1+x)^{τ1+1}/x and dp/dε≈α1 ε0/(ε−ε0)^2. For α1<0 the pressure is positive but dp/dε<0; for α1>0 the pressure is negative but dp/dε>0. Since all tabulated sets have α1≠0, no set satisfies both positivity and dp/dε≥0 across the matching region. High-density behavior also fails for some sets: for example, Table IV, 0<τ1, σ1<0, row 3 (α1=0.70, β1=2.21, τ1=0.01, σ1=−0.21) gives p→−∞ as ε→∞ because the −α1(ε/ε0)^{τ1} term dominates. The results therefore are not physically admissible TOV solutions of the stated EoS.
- [Eq. (13)] I could not reproduce inequality (13) from the stated EoS. Direct differentiation of the core term (10) with ε0=1 gives dp/dε = α1(1+x)^{τ1}(1−τ1 x)/x² + β1σ1(1+x)^{σ1−1}, which does not reduce to (1+x)^{τ1−σ1+2}(τ1 x−1) ≤ β1σ1/α1 ε². The powers of (1+x) and the role of the factor ε² are inconsistent, and the claimed condition σ1≠τ1 as a divergence requirement is not evident. Because the critical-density formula (14) and the regime classification in Table I rest on this calculation, the algebraic basis for those results is currently unsupported.
- [Table III and abstract] The text above Table III states that the listed sets 'produce stable configurations with masses ranging from 0.99 to 2.05 M⊙', and the abstract repeats the range (0.99−2.05) M⊙. However, the first row of Table III gives M=0.79 M⊙. This is not a harmless typo: the low-end mass is a central advertised outcome, and the table is the only place the reader can check which configurations actually exist. The discrepancy must be resolved, and the status of the 0.79 M⊙ row must be clarified.
minor comments (4)
- [Tables III–V, notation] The symbol εc/ε0 is used in the tables for the central density, while εc in Eq. (14) denotes a critical density. This collision makes the tables and the text hard to read; please use different symbols, e.g., ε_center and ε_crit.
- [Eq. (14)] The denominator 2τ1(σ1−τ1) vanishes for τ1=0 or σ1=τ1. These cases are excluded only implicitly; they should be discussed explicitly, since some tabulated sets have τ1=0.01 or τ1=0.02 and the condition σ1≠τ1 is not always stated.
- [Section IV.B, Eq. (42)] The baryon density n_B is determined only up to a multiplicative constant by (42), but the chemical potential µ=(ε+p)/n_B is subsequently treated as having a well-defined magnitude. The chosen normalization should be stated, because an arbitrary rescaling of n_B changes µ by a constant factor.
- [Section II.A, Eq. (7)] The approximation in Eq. (7) is written as a dimensionless ratio that is said to be '≈0'. As written it has no clear expansion parameter; please state explicitly that it requires |τ1|≪1 and specify the typical values used.
Circularity Check
Low-mass 'prediction' is a parameter-search target; chemical-potential signatures restate the input EoS's curvature.
-
fitted input called prediction
[Sec. IV.A (Parameter regimes and mass–radius relations); Abstract; Tables III–V]
"For each regime specified in Table I, we adjust parameters to produce NS configurations with realistic densities and masses, including low-mass configurations. ... The resulting mass–radius sequences yield low neutron star masses (0.99−2.05)M⊙"
The abstract presents the low-mass range as a derived result, but the parameter selection is explicitly an inverse search: Tables III–V are chosen by adjusting α1,β1,τ1,σ1,α2,β2,τ2,σ2 to yield realistic and low-mass configurations. The TOV integration then returns masses near the targeted range. No fixed, parameter-free EoS is used, and no independent benchmark fixes the parameters before integration. Thus the reported 0.99–2.05 M⊙ output is equivalent to the fitting target by construction, not an independent prediction of the model.
-
self definitional
[Sec. IV.B (Phase-transition signatures and chemical potential), Eqs. (41)–(43), Fig. 7; Sec. II.A Eq. (14)]
"µ(ε) = (ε+p(ε))/n_B(ε) ... The resulting curves show that the most significant deviations from monotonicity occur in the parameter regions where two positive critical densities exist, in agreement with the classification of Fig. 2."
The chemical potential is not an independent observable: Eq. (42), dn_B/n_B = dε/(ε+p), and Eq. (43), µ=(ε+p)/n_B, give dµ/dε = (dp/dε)/n_B. Therefore the sign and curvature features of µ(ε) are exactly the features of the input EoS p(ε), whose zeros of dp/dε define the critical densities (14) used to build Fig. 2 and Table I. Claiming that the µ curves 'agree with' the critical-density classification is a restatement of the same derivative of the same function, not an independent confirmation of phase-transition-like behavior.
full rationale
No self-citation chain or imported uniqueness theorem is load-bearing here; the paper's references to Alford et al. are external and not used to force the central construction. The main circularity is in the low-mass claim: Section IV.A explicitly states that parameters are adjusted to produce realistic and low-mass configurations, and the abstract then presents the resulting 0.99–2.05 M⊙ range as an outcome. Because Tables III–V are the product of that inverse search, the TOV masses are selected-for outputs rather than predictions from a fixed EoS. The phase-transition-like interpretation is also internal to the ansatz: critical densities are roots of dp/dε=0 for the modified van der Waals form, and the chemical-potential 'signatures' are computed from that same p(ε), so the agreement with the critical-density classification is definitional (dµ/dε = p'/n_B). The paper is transparent that the model is a simplified testbed 'not intended to reproduce detailed microphysics,' which mitigates but does not remove the circularity of presenting fitted targets and built-in curvature as findings. A separate mathematical difficulty—the pole of (ε0/ε−1)^{-1} at ε0 and the undefined interval between ε0 and the ε+ used in matching—is a correctness concern rather than a circularity concern and is not scored here.
Assumptions & free parameters
free parameters (8)
- ε0 (saturation/reference density) =
set to 1 (natural units)
- α1 (core repulsion amplitude) =
e.g., 1.75, -0.37, 0.70, -1.37 across Tables III–V
- β1 (core interaction amplitude) =
e.g., 0.25, 2.01, 2.40, 2.50
- τ1 (core repulsion exponent) =
e.g., -1.89, 3.25, 0.10, 0.20
- σ1 (core interaction exponent) =
e.g., 2.50, -2.50, -0.01, 2.50
- α2, β2, τ2, σ2 (crust parameters) =
listed per row in Tables III–V
- δ (core–crust pressure mismatch) =
δ≪1, unquantified
- n_B integration constant C =
arbitrary (sets absolute scale of μ)
assumptions (5)
- standard math Tolman–Oppenheimer–Volkoff equations describe hydrostatic equilibrium of the star
- domain assumption Matter is cold, β-equilibrated, and approximated as a neutral free Fermi gas
- ad hoc to paper The star can be split at ε0 into a core EoS and a crust EoS with approximate matching (12), δ≪1
- ad hoc to paper A density-dependent temperature T=T0(ε/ε0)^{τ1} with |τ1|≪1 and Eq. (7) approximated as ≈0
- ad hoc to paper All tabulated parameter sets satisfy dp/dε≥0 and vs²≤1 over the star
Cite this review
Pith. "Pith review of Low-Mass Neutron Stars and Effective Phase Transitions from a Hybrid Van der Waals-Polytropic Equation of State." pith.science (2026). https://pith.science/paper/O5KX24XC
@misc{pith2026251208672,
author = {Pith},
title = {Pith review of: Low-Mass Neutron Stars and Effective Phase Transitions from a Hybrid Van der Waals-Polytropic Equation of State},
year = {2026},
howpublished = {\url{https://pith.science/paper/O5KX24XC}},
note = {Machine review of arXiv:2512.08672}
}
abstract
We study phase-transition-like behavior in neutron stars using a simplified, piecewise equation of state that couples a modified van der Waals-type core to a polytropic crust. The model remains analytically tractable while allowing for nonlinear density dependence. We impose thermodynamic and causal consistency conditions and determine the critical densities at which the curvature of the pressure-energy density relation changes. In the non-relativistic limit, the generalized Lane-Emden equations describe a smooth core-crust transition layer. We integrate the Tolman-Oppenheimer-Volkoff equations across different $(\tau_1,\sigma_1)$ regimes, where these parameters encode thermal and interaction effects in the core. The resulting mass-radius sequences yield low neutron star masses $(0.99-2.05)M_{\odot}$, and the chemical potential exhibits the characteristic signatures of phase-transition behavior at densities well above the matching point. Our results show that analytic EOS models can reproduce the key phenomenology of phase transitions and provide a controlled framework for exploring low-mass neutron star configurations.
Figures
Figures from the paper (5 more)
Forward citations
Cited by 1 Pith paper
-
Strongly interacting matter with criticality induced by modified excluded volume in core-collapse supernova simulations
Modified-excluded-volume EOS with continuous van der Waals phase transition yields CCSN explosions and a several-millisecond neutrino burst, while the same EOS under Gibbs construction fails to explode.
Reference graph
Works this paper leans on
-
[1]
R. C. Tolman,Effect of inhomogeneity on cosmological models,Proc. Nat. Acad. Sci.20(1934) pg. 3
1934
-
[2]
R. C. Tolman,Static solutions of Einstein ’s field equations for spheres of fluid,Phys. Rev.55 (1939) pg. 364
1939
-
[3]
J. R. Oppenheimer and G. M. Volkoff,On massive neutron cores,Phys. Rev.55(1939) pg. 334
1939
-
[4]
Fermi,Zur Quantelung des idealen einatomigen Gases,Zeitschrift f¨ ur Physik36(1926) pg
E. Fermi,Zur Quantelung des idealen einatomigen Gases,Zeitschrift f¨ ur Physik36(1926) pg. 902
1926
-
[5]
C. G. Bassaet al.,LOF AR discovery of the fastest-spinning millisecond pulsar in the galactic field,Astrophys. J. Lett.846(2017) pg. L20 27
2017
-
[6]
Antoniadiset al.,A massive pulsar in a compact relativistic binary,Science340(2013) pg
J. Antoniadiset al.,A massive pulsar in a compact relativistic binary,Science340(2013) pg. 1233232
2013
-
[7]
Linares, T
M. Linares, T. Shahbaz and J. Casares,Peering into the dark side: magnesium lines establish a massive neutron star in PSR J2215+5135,Astrophys. J.859(2018) pg. 1
2018
-
[8]
M. C. Milleret al.,PSR J0030+0451 mass and radius from NICER data and implications for the properties of neutron star matter,Astrophys. J. Lett.887(2019) pg. L24
2019
Show all 31 references
-
[9]
T. E. Riley et al.,A NICER view of PSR J0030+0451: millisecond pulsar parameter estima- tion,Astrophys. J. Lett.887(2019) pg. L21
2019
-
[10]
M. C. Milleret al.,The radius of PSR J0740+6620 from NICER and XMM-Newton data, Astrophys. J. Lett.918(2021) pg. L28
2021
-
[11]
T. E. Rileyet al.,A NICER view of the massive pulsar PSR J0740+6620 informed by radio timing and XMM-Newton spectroscopy,Astrophys. J. Lett.918(2021) pg. L27
2021
-
[12]
Doroshenkoet al.,A strangely light neutron star within a supernova remnant,Nature Astronomy6(2022) pg
V. Doroshenkoet al.,A strangely light neutron star within a supernova remnant,Nature Astronomy6(2022) pg. 1444
2022
-
[13]
J. E. Horvathet al.,A light strange star in the remnant HESS J1731-347: Minimal consistency checks,A&A672(2023) pg. L11
2023
-
[14]
isolated
J. Linet al.,An X-Ray-dim “isolated” neutron star in a binary?,Astrophys. J. Lett.944 (2023) pg. L4
2023
-
[15]
M¨ uller, A
B. M¨ uller, A. Heger and J. Powell,Minimum neutron star mass in neutrino-driven supernova explosions,Phys. Rev. Lett.134(2025) pg. 071403
2025
-
[16]
Wang and A
T. Wang and A. Burrows,Supernova explosions of the lowest-mass massive star progenitors, Astrophys. J.969(2024) pg. 74
2024
-
[17]
B. P. Abbottet al.,GW170817: measurements of neutron star radii and equation of state, Phys. Rev. Lett.121(2018) pg. 161101
2018
-
[18]
B. P. Abbottet al.,Binary black hole population properties inferred from the first and second observing runs of advanced LIGO and advanced Virgo,Astrophys. J. Lett.882(2019) pg. L24
2019
-
[19]
A. G. W. Cameron,Neutron star models,Astrophys. J.130(1959) pg. 884
1959
-
[20]
G. Baym, C. Pethick, and D. Pines,Superfluidity in neutron stars,Nature224(1969) pg. 673
1969
-
[21]
G. Baym, C. Pethick, and P. Sutherland,The ground state of matter at high densities: equation of state and stellar models,Astrophys. J.170(1971) pg. 299 28
1971
-
[22]
Akmal, V
A. Akmal, V. R. Pandharipande, and D. G. Ravenhall,Equation of state of nucleon matter and neutron star structure,Phys. Rev. C58(1998) pg. 1804
1998
-
[23]
A. R. Raduta, F. Nacu, and M. Oertel,EoS for hot neutron stars,Eur. Phys. J. A57(2021) pg. 329
2021
-
[24]
Ji and J
Z. Ji and J. Chen,The equation of state of neutron stars: theoretical models, observational constraints, and future perspectives, arXiv:2502.05513v1
-
[25]
Baacke,Thermodynamics of a gas of MIT bags,Acta Phys
J. Baacke,Thermodynamics of a gas of MIT bags,Acta Phys. Pol.88(1977) pg. 625
1977
-
[26]
Chodos, R
A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn and V. F. Weisskopf,A new extended model of hadrons,Phys. Rev. D9(1974) pg. 3471
1974
-
[27]
Alvarez-Castillo, S
D. Alvarez-Castillo, S. Benic, D. Blaschke, S. Han and S. Type,Neutron star mass limit at 2 M⊙ supports the existence of a CEP,Eur. Phys. J. A52(2016) pg. 232
2016
-
[28]
Vovchenko,Hadron resonance gas with van der Waals interactions,Int
V. Vovchenko,Hadron resonance gas with van der Waals interactions,Int. J. Mod. Phys. E 29(2020) pg. 2040002
2020
-
[29]
Verma, A
A. Verma, A. K. Saha and R. Mallick,Comparison of equations of state for neutron stars with first-order phase transitions: a qualitative study,Astrophys. J.985(2025) pg. 1
2025
-
[30]
M. G. Alford, S. Han and M. Prakash,Generic conditions for stable hybrid stars,Phys. Rev. D88(2013) pg. 083013
2013
-
[31]
M. G. Alford, G. F. Burgio, S. Han, G. Taranto and D. Zappal` a,Constraining and applying a generic high-density equation of state,Phys. Rev. D92(2015) pg. 083002
2015
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.