{"id":"29110989-8d97-47f5-bf7f-e2b1a44b7166","arxiv_id":"2512.08672","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A hand-tuned van der Waals-plus-polytrope equation of state can generate low-mass neutron-star configurations and curvature-only 'phase-transition' fingerprints in the chemical potential.","lead":"This study builds a two-part toy equation of state—a van der Waals-like core matched to a polytropic crust—and integrates the relativistic structure equations to obtain neutron-star mass-radius sequences down to roughly one solar mass. The authors argue that simplified analytic models can reproduce the qualitative thermodynamic fingerprints of weak phase transitions, but the key parameters are hand-adjusted rather than derived from observations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pressure pole at ε0 invalidates the claimed smooth core–crust transition and any TOV solutions using the stated piecewise EoS.","rationale":"The central claim—that a simple analytic EoS can produce low-mass neutron-star sequences and phase-transition-like diagnostics—depends entirely on the TOV integrations reported in Tables III–V. Those integrations use a piecewise EoS whose core branch has a pressure pole at the nominal matching density ε0. This is not a matter of tuning or missing numerical audits; it makes the stated model ill-defined at the core–crust boundary. The paper's own interface analysis (Section IV.B) claims a smooth, continuous transition, which is impossible given the divergent core pressure. The approximate matching condition (12) avoids the pole only by evaluating the core at densities far above ε0, leaving a density interval with no EoS at all. This invalidates the reported mass–radius relations, chemical-potential curves, and phase-transition interpretations as consequences of the written model. I agree with the reader that the matching condition and stability/causality audit are the weak spot, but the failure is sharper than 'not demonstrated': the core branch provably violates Eq. (1) near ε0 for the published parameters. A fix may be possible by regularizing the vdW term or implementing a genuine Maxwell construction, but as written the central claim is unsupported. The lack of released code or data makes it impossible to check what the authors actually integrated, which further weakens confidence. This is an internal consistency objection, not a disagreement with the field's consensus, and it is checkable by direct evaluation of Eqs. (8) and (9).","tokens_in":1132,"tokens_out":1546,"duration_ms":120155,"concrete_test":"Reproduce the TOV integration exactly as stated, switching from Eq. (8) to Eq. (9) at ε=ε0, for Table V row 1 parameters (α1=-0.22, β1=2.43, τ1=0.10, σ1=-0.01; α2=β2=0.19, τ2=0.25, σ2=2.21). Show that no finite central density yields a solution reaching ε=0: as r increases and ε decreases toward ε0, hydrostatic equilibrium requires dP/dr<0, but P→∞ at ε0, so the integration cannot pass the interface. Then repeat with a minimal regularization, replacing the core denominator by (ε0/ε−1)+λ with λ>0, and check whether the published 0.99–2.05M⊙ sequences survive in the λ→0 limit. If they require λ≈0, the reported mass–radius curves and phase-transition diagnostics are artifacts of the unmodeled pole.","verdict_should_be":"REJECT","load_bearing_attack":"The core branch of the piecewise EoS, Eq. (8), contains the factor α1/(ε0/ε−1)·(ε/ε0)^{τ1}. At the nominal matching density ε=ε0, the denominator vanishes. For α1<0, p_core→+∞ as ε→ε0+; for α1>0, p_core→−∞. All parameter sets in Tables III–V have α1≠0, and the crust branch is finite at ε0 (α2+β2 after ε0=1). Thus the approximate matching condition (12) with δ≪1 cannot hold on any interval containing ε0; the paper sidesteps this by requiring ε0/ε+≪1, i.e., evaluating the core only at ε+≫ε0 and leaving the interval (ε0, ε+) with no defined EoS. Section IV.B nevertheless claims that η is continuous across r1 and that A(η)>0, describing a smooth diffusive transition rather than a discontinuity. A continuous density profile crossing ε0 would require finite p_core there. Moreover, near ε0 the core branch violates the thermodynamic stability condition (1): for Table V row 1, p_core(1.001)≈222 while p_core(2)≈2.9, so dp/dε<0 over a broad density interval. Therefore the TOV solutions in Tables III–V are not solutions of the stated EoS; they require an implicit cutoff or an unstated density jump that the paper does not define. This is a stronger version of the reader's concern: the model is not merely unaudited, it is internally inconsistent at the matching point.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":15376,"tokens_out":10282,"duration_ms":103394,"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":[{"comment":"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.","section":"Eqs. (8), (12); Section IV.B"},{"comment":"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.","section":"Eqs. (10)–(11), near ε0"},{"comment":"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.","section":"Eq. (13)"},{"comment":"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.","section":"Table III and abstract"}],"minor_comments":[{"comment":"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.","section":"Tables III–V, notation"},{"comment":"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":"Eq. (14)"},{"comment":"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":"Section IV.B, Eq. (42)"},{"comment":"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.","section":"Section II.A, Eq. (7)"}],"recommendation":"reject","confidential_remarks":"The singular matching-point defect is decisive: the piecewise EoS as written does not define a pressure for the density interval connecting core and crust, and the tabulated TOV solutions therefore are not solutions of the stated equations. The instability/negative-pressure behavior near ε0 worsens the problem. These are load-bearing issues rather than presentation problems; a repair would require redefining the functional form of the core EoS and recomputing all numerical results, which goes beyond a revision of the present manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline: this paper has a load-bearing singularity. The core EoS (Eq. 8) diverges at the matching density ε0, yet the authors present smooth core–crust transitions and TOV solutions built from it. As written, the numerical results in Tables III–V are not solutions of the stated piecewise EoS.\n\nThat said, the genuinely new piece is the classification of the (τ1, σ1) parameter plane into regimes with zero, one, or two critical densities. That mapping is transparent and, as far as I can tell, not in the cited literature. The Lane-Emden treatment of the transition layer is a reasonable exercise, and the chemical-potential analysis is a sensible diagnostic even if it yields no plateaus. If the construction were regularized, it could be a useful toy testbed for low-mass NS configurations.\n\nThe pole is the main problem. With α1 ≠ 0 in every listed set, the core pressure goes to ±∞ as ε → ε0+. The paper only evaluates the core at densities ε+ with ε0/ε+ ≪ 1, leaving the interval (ε0, ε+) with no defined EoS. That's not a smooth transition; it's an undefined gap. The TOV curves cannot be trusted.\n\nSmaller issues: Eq. (13) is algebraically garbled and doesn't match a direct derivative. Table III's first row gives 0.79 M☉, contradicting the abstract's 0.99–2.05 M☉ range. The claimed dp/dε ≥ 0 and causality conditions are not audited for all parameter sets—at least one row in Table V shows a wide interval with dp/dε < 0. The low-mass configurations are partly a fitted outcome, since the authors say they adjust parameters to produce realistic densities and masses. No code or data are provided, so the curves can't be independently checked.\n\nOn the positive side, the paper is honest about its scope: it repeatedly says this is not microphysics, just a controlled toy. It also acknowledges that μ(ε) shows no true phase-transition plateau. The singularity is fixable in principle—replace the (ε0/ε − 1) denominator with a regularized form and define a matching density—but that is substantial revision, not copyediting.\n\nWho should read this? Someone cataloging analytic toy EoS might note the classification as a cautionary example. I wouldn't cite it in its current form, and I wouldn't bring it to reading group without the objection attached. The internal inconsistency is severe enough that I'd desk reject rather than spend referee time, though I'd allow resubmission if the EoS is regularized and the TOV runs redone with a clear matching and a full stability/causality audit.\n\nMy recommendation: reject for now, with a path to resubmit.\n\nRegards,","headline":"A load-bearing singularity in the core EoS at the matching density invalidates the TOV results as written.","tokens_in":15742,"tokens_out":8340,"would_cite":false,"duration_ms":70405,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["97.60.Jd","26.60.+c","26.60.-c"],"model":"deepseek-v4-flash","headline":"A two-branch analytic equation of state can produce neutron stars as light as 0.99 solar masses.","keywords":["neutron star","equation of state","phase transition","van der Waals","Tolman-Oppenheimer-Volkoff","low-mass neutron star","polytrope","chemical potential"],"falsifier":"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.","tokens_in":14738,"feed_emoji":"⭐","tokens_out":5095,"duration_ms":45525,"temperature":0.7,"pith_summary":"This paper argues that a deliberately simple, analytically tractable equation of state — a modified van der Waals-style core matched to a polytropic crust — can reproduce the qualitative signatures of phase transitions in neutron stars and yield stellar masses between 0.99 and 2.05 solar masses. The authors show that just two exponents governing the core pressure control whether the pressure–energy-density curve has one inflection, two inflection points, or none, and they interpret these regimes as effective phase-transition diagnostics. They then integrate the Tolman–Oppenheimer–Volkoff equations and find families of stable configurations, some with mass near 1.0–1.1 solar masses and radius near 11 km, which line up with recent low-mass neutron star candidates. The baryonic chemical potential stays smooth but changes slope most sharply at densities far above the artificial core–crust junction, indicating that any real transition would lie deep inside the star. The upshot: qualitative neutron-star phenomenology does not require complex microphysics, and this minimal model offers a controlled setting to explore it.","feed_headline":"Simple equation of state yields neutron stars as light as 0.99 suns","feed_subtitle":"Two core parameters control stiffness and mimic phase transitions in mass-radius models.","key_machinery":"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","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Tiny neutron stars from a simple two-piece equation","Phase-transition mimicry in neutron stars via a van der Waals core","Low-mass neutron stars emerge from an analytic EoS model","A simple EoS reproduces neutron star phase-transition signs","Neutron stars down to 0.99 solar masses from a hybrid EoS"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Tiny neutron stars from a simple two-piece equation","Phase-transition mimicry in neutron stars via a van der Waals core","Low-mass neutron stars emerge from an analytic EoS model","A simple EoS reproduces neutron star phase-transition signs","Neutron stars down to 0.99 solar masses from a hybrid EoS"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000781,"raw_usage":{"total_tokens":3294,"prompt_tokens":756,"completion_tokens":2538,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":500,"completion_tokens_details":{"reasoning_tokens":2447}},"tokens_in":500,"tokens_out":2538,"duration_ms":16004,"temperature":1.0,"reasoning_tokens":2447,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T17:36:04.601154+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}