{"id":"2b8a9951-deb9-4b35-bc1c-5e5a78210f66","arxiv_id":"2507.06837","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using the Gibbs construction for the phase transition in the CCT model, the authors find a mixed-phase region extending down to the inner crust, giving neutron stars a thick, massive crust.","lead":"A nuclear matter model already used to explain the small, light neutron star HESS J1731-347 is reanalyzed with a more correct treatment of a phase transition. The result is a predicted thick and massive crust that changes how this object's mass, radius, and spin evolution would be interpreted.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that the Gibbs mixed phase is a thick crust rests on neglecting surface and Coulomb effects that are not small away from the splitting point, so the quantitative crustal properties are not yet substantiated.","rationale":"","tokens_in":8697,"tokens_out":7732,"duration_ms":77944,"concrete_test":"Recompute the mixed-phase EOS for the C_sigma^2=12 model including standard surface and Coulomb terms in a Wigner-Seitz approximation: add epsilon_surf = (3 sigma / r_c) * w (or the appropriate geometric factor) and epsilon_Coul to Eq. (18), minimize the cell size r_c and volume fraction w at each average density, with surface tension sigma ranging from 10 to 70 MeV/fm^2. Then solve the TOV equations with this EOS and recompute the crustal thickness, mass fraction, and moment of inertia. If these quantities change by more than 20% relative to Fig. 5, or if the two-phase region no longer extends to the low njoin values in Table II, the central claim of a thick crust is not robust and the paper's stated justification is incomplete.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the CCT model with Gibbs construction yields a thick, massive crust depends on identifying the entire two-phase region (Eqs. 17-18, Fig. 2) as the inner crust. The paper itself states in Section III that finite-size (surface and Coulomb) effects are not included, arguing that they do not change the density reached by the two-phase system because the density difference between phases vanishes at the splitting point. That argument only protects the upper onset density n2ph, not the lower-density part of the mixed phase where the phase densities differ substantially. Away from the splitting point, surface and Coulomb energies are not negligible; they alter the energy and pressure of the mixed phase, which directly affects the TOV solution and hence the crustal thickness, mass, and moment of inertia shown in Fig. 5. Moreover, without finite-size effects the Gibbs mixed phase is a homogeneous mixture, not a lattice of quasi-nuclei; the identification of the proton-rich phase with 'quasi-nuclei' and the inner crust structure (Section III) is an interpretive leap that is not derived from the calculation. If finite-size effects are included, the EOS may change enough that the crustal properties quoted in the abstract and Fig. 5 are not robust, and the phase boundary itself could shift, undermining the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the CCT relativistic mean-field model, previously used to explain the compact low-mass object HESS J1731-347, and replaces the Maxwell construction with a Gibbs construction for the phase transition in beta-equilibrated nucleon-lepton matter. The authors find that the system separates into a low-density neutron-lepton phase and a high-density neutron-proton-lepton phase over densities from about 0.0025 to 0.0656 fm^-3 up to about 0.26 to 0.30 fm^-3. Interpreting the entire two-phase region as the inner crust, they solve the Tolman-Oppenheimer-Volkoff equations and report that for HESS J1731-347-like masses the crust is about 50% thicker, several times more massive, and has a much larger moment of inertia than in standard crust models. Mass-radius relations for three parametrizations are compared with NICER, X-ray burst, and gravitational-wave constraints.","tokens_in":9048,"tokens_out":8641,"duration_ms":101434,"significance":"If the central claim were fully substantiated, the result would be interesting: it would connect a non-exotic, soft low-density equation of state with a qualitatively different prediction for crustal properties, and it would give a concrete interpretation of HESS J1731-347's compactness in terms of a mixed-phase region rather than exotic matter. The manuscript has clear strengths: the Gibbs construction is applied in a thermodynamically consistent way for a system with two conserved charges, the phase diagram and EOS are presented quantitatively in tables and figures, and the mass-radius curves are confronted with a wide set of observations. However, the central claim that this object has a thick crust is not yet established because the calculation neglects surface and Coulomb effects and because the identification of the homogeneous mixed phase with a crustal microstructure is an interpretive step rather than a derived result.","major_comments":[{"comment":"The central quantitative claims—crust thickness, crustal mass fraction, and moment-of-inertia fraction—are computed from the Gibbs mixed-phase EOS without surface or Coulomb terms. The paper's justification that finite-size effects do not change \"the density reached by the two-phase system\" protects only the upper endpoint n2ph, where the densities of the two phases merge; in the lower-density part of the mixed phase, down to njoin ~ 0.0025-0.0656 fm^-3, the phase densities differ substantially, so surface and Coulomb energies are not negligible and would alter the pressure-density relation used in the TOV integration. The authors should either include a finite-size (pasta) treatment or demonstrate quantitatively that the neglected terms leave the crustal properties in Fig. 5 unchanged.","section":"Section III, Eqs. (17)-(18) and Fig. 5"},{"comment":"The statement that the residual proton-neutron phase \"may be interpreted as quasi-nuclei sparsely placed in the region dominated by neutron fluid; such a system corresponds to the inner crust structure\" is not derived from the calculation. In the absence of finite-size effects, the two-phase system is a homogeneous mixture of a positively charged npl phase and a negatively charged nl phase; it is not a lattice of quasi-nuclei embedded in a neutron gas. The negatively charged phase is not the standard inner-crust neutron gas, and the charge-separated microstructure (droplets, rods, slabs) that defines a crust is precisely what the neglected surface and Coulomb terms would produce. The title's claim that HESS J1731-347 has a thick crust therefore requires either a finite-size calculation or a more careful restatement that the model predicts a mixed-phase region, not a crust in the usual sense.","section":"Section III, inner-crust interpretation"},{"comment":"The model parameters, in particular the symmetry-energy slope L = 40 MeV, are chosen to reproduce the compactness of HESS J1731-347, and the three C_sigma^2 values vary only the stiffness at fixed L. No sensitivity study over L or the other couplings is presented, and no uncertainties from the HESS J1731-347 measurement are propagated. Since the existence and density range of the two-phase region depend on the low-density symmetry energy, the answer to the title's question is model-dependent. The authors should show how the crustal properties change when L is varied within its experimentally allowed range, or at least identify the range of L for which the thick-crust phenomenon persists.","section":"Section II, Table I, and Section IV"}],"minor_comments":[{"comment":"The text gives \"ncc = 0.056 fm2\"; the units should be fm^-3.","section":"Section IV"},{"comment":"The SLy4 EOS is spelled \"Sly4\" in one place; please use a consistent spelling.","section":"Section III"},{"comment":"The njoin values for C_sigma^2 = 13 and 14 are identical to four significant digits (0.0656 fm^-3); please clarify whether this is a coincidence or a typographical error.","section":"Table II"},{"comment":"The caption says \"energy difference\" while the text says \"energy of the two-phase system and that of homogeneous, beta-equilibrated matter\"; please clarify the quantity and sign convention shown in the lower panel.","section":"Figure 2"},{"comment":"The speed-of-sound panel would benefit from a legend distinguishing the Gibbs and Maxwell curves, and from an explicit check that the plotted values satisfy the inequality in Eq. (19) over the entire density range.","section":"Figure 3"},{"comment":"The curves for the three C_sigma^2 values are not labelled directly in the figure; adding labels or a legend would improve readability.","section":"Figure 5"},{"comment":"The conclusion that a thick crust \"should significantly alter the thermal evolution and rotational properties\" is not accompanied by any quantitative estimate; a benchmark calculation (for example, of the crustal moment of inertia or a cooling indicator) would make the observational consequences more concrete.","section":"Section V"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the Gibbs construction is handled cleanly, but the abstract and title claim more than the calculation supports. The finite-size neglect is acknowledged in Section III, yet it directly affects the computed crustal mass, thickness, and moment of inertia, and the identification of the mixed phase with an inner crust is an interpretive leap. I would ask the authors either to add a finite-size treatment or to reframe the claims as properties of a mixed-phase region rather than a crust. A sensitivity scan over L would also strengthen the paper, since the model is tuned to reproduce HESS J1731-347. The manuscript is promising but needs revision before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is a legitimate Gibbs-versus-Maxwell point: for npl matter with two conserved charges, the Maxwell construction breaks chemical equilibrium, and the authors correctly apply the Gibbs conditions. That is not new by itself, but applying it to their CCT model gives a genuinely new result: the two-phase region persists all the way down to densities around 0.0025–0.066 fm^-3, far below the ~0.1 fm^-3 lower bound of their earlier paper [17]. So the claim that this model has a crust-core transition near 0.26–0.30 fm^-3 is new and worth taking seriously.\n\nSecond thing: treat the quantitative crustal properties as provisional. The paper states in Section III that finite-size (surface and Coulomb) effects are not included, arguing they don't change the density reached by the two-phase system because the density difference between phases vanishes at the splitting point. That argument only protects the upper onset density n2ph. Away from that point, the two phases have substantially different densities, so surface and Coulomb energies are not negligible. They would change the energy and pressure of the mixed phase, which feeds directly into the TOV solution and hence the crust thickness, mass, and moment of inertia shown in Fig. 5. Without those terms, the mixed phase is a homogeneous mixture, not a lattice of quasi-nuclei—so calling the whole region an \"inner crust\" is an interpretation, not a derivation.\n\nThere is also a mild circularity: L=40 was chosen to reproduce HESS J1731-347's compactness, and then the same model is used to infer the crust properties of that same object. That is not fatal, but it means the quantitative predictions are tuned to the target. No error bars are propagated either, though the Cσ^2 variation gives some sense of robustness.\n\nWhat the paper does well: the calculation is internally consistent, the Gibbs construction is applied in a standard way, the mass-radius curves still respect the major constraints (with the usual caveat about PSR J1231-1411), and the authors are transparent about what they left out. The distinction between [17] and this work is clear.\n\nWho is this for: nuclear astrophysicists working on EOS, crust structure, and NICER/gravitational-wave interpretations. A serious referee should engage with it. The main thing to ask for is either a finite-size calculation for the mixed phase or at least a quantitative estimate of its effect on the EOS and crustal observables. As it stands, the qualitative idea is plausible and the quantitative claim is unproven. I'd send it out, and I'd be skeptical until the finite-size issue is resolved.","headline":"The Gibbs treatment is right and the two-phase region reaching the crust is genuinely new, but the thick-crust numbers are not yet substantiated because finite-size effects are neglected away from the splitting point.","tokens_in":26,"tokens_out":2623,"would_cite":false,"duration_ms":38076,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Using thermodynamically consistent Gibbs conditions instead of the Maxwell construction in the CCT nuclear model pushes the crust-core boundary of HESS J1731-347 to roughly twice saturation density, making the crust about 50 percent…","keywords":["HESS J1731-347","neutron star crust","Gibbs construction","phase transition","relativistic mean-field model","equation of state","meson crossing terms","symmetry energy"],"falsifier":"A calculation of the finite-size structured mixed phase (nuclear pasta) for the CCT equation of state over the density range 0.0025 to 0.30 fm$^{-3}$ would settle the quantitative claim: if inclusion of surface and Coulomb energies shifts the crust-core boundary below roughly 0.1 fm$^{-3}$, or changes the crustal moment-of-inertia fraction outside the predicted several-times-larger range, the thick-crust conclusion fails. Alternatively, a precise glitch measurement on a star with the mass of HESS J1731-347 that implies a crustal moment of inertia comparable to thin-crust models would contradict the prediction.","tokens_in":8550,"feed_emoji":"⭐","tokens_out":9157,"duration_ms":79596,"temperature":0.7,"pith_summary":"This paper argues that the CCT model of nuclear matter, a relativistic mean-field model with meson crossing terms that can explain the unusually small and light neutron star HESS J1731-347, has been solved with the wrong phase-transition rule. The Maxwell construction used earlier cannot enforce chemical equilibrium when matter carries two conserved charges, so the paper redoes the phase transition with the thermodynamically consistent Gibbs conditions. Doing so splits the dense matter into a proton-free neutron-lepton phase and a neutron-proton-lepton phase over a wide density range, from the deep inner crust up to roughly twice saturation density. The paper interprets this two-phase region as the star's crust, which makes the crust of HESS J1731-347 about 50 percent thicker and several times more massive, with a much larger moment of inertia, than standard crust models predict. That matters because crustal mass and moment of inertia control thermal evolution, rotation, and glitch behavior.","feed_headline":"Thermodynamically correct phase rules thicken HESS J1731-347's crust 50%","feed_subtitle":"The crust-core boundary shifts to about twice saturation density, multiplying crustal mass and moment of inertia.","key_machinery":"The load-bearing object is the Gibbs construction for two-phase matter with two conserved charges (baryon number and electric charge), applied to the CCT relativistic mean-field Lagrangian that includes $\\sigma$-$\\delta$ and $\\omega$-$\\delta$ meson crossing terms. The paper requires $P^{(I)}=P^{(II)}$ and equality of chemical potentials $\\mu^{(I)}_i = \\mu^{(II)}_i$ for every species present in both phases; because protons cannot satisfy this, they are confined to the proton-rich phase (inequality (15)), while neutrons and leptons equilibrate across the phases. Global charge neutrality then fixes the volume fraction $w = V^{(I)}/V$ through Eq. (16), and the instability criterion $K_\\mu = (\\partial P/\\partial n)_\\mu < 0$ marks where the homogeneous phase must split. This machinery converts a single first-order transition into a broad mixed-phase band, and it is the reason the crust-core boundary moves up to about twice saturation density.","core_discovery":"The central claim is that applying the proper Gibbs phase-equilibrium conditions to the CCT model changes where the neutron star crust ends. The authors find that $\\beta$-equilibrated matter is unstable for charge fluctuations over a much broader density interval than previously thought: the two-phase system forms at $n_{2\\mathrm{ph}} \\simeq 0.26$ to $0.30\\,\\mathrm{fm}^{-3}$ and persists down to $n_{\\mathrm{join}} \\simeq 0.0025$ to $0.0656\\,\\mathrm{fm}^{-3}$, where it is joined to the SLy4 crust equation of state. The coexisting phases are a low-density neutron-lepton phase with no protons and a higher-density phase containing protons and leptons, with volume fraction fixed by global charge neutrality; the proton-free phase is described as quasi-nuclei embedded in a neutron fluid, the usual picture of inner-crust matter. Treating the entire two-phase region as the crust places the crust-core transition at roughly $2n_0$ instead of the standard $0.03$ to $0.08\\,\\mathrm{fm}^{-3}$, and for HESS J1731-347-like stars this makes the crust about 50 percent thicker and raises its mass and moment-of-inertia contributions by a factor of several.","pith_inferences":["If the same Gibbs construction is applied to other relativistic mean-field models with a low symmetry-energy slope, the charge-fluctuation instability may produce comparably thick crusts, so the result may be generic rather than specific to the CCT model.","A concrete next step would be computing the finite-size structured mixed phase (nuclear pasta) over the two-phase density range; the paper argues the onset density is unaffected, but crustal mass and moment of inertia would shift if surface and Coulomb energies are significant.","Observational discrimination could come from pulsar glitches: the predicted several-fold larger crustal moment of inertia would change the glitch-activity pattern, so a precise measurement on a compact low-mass pulsar could support or exclude the thick-crust identification."],"forward_implications":["For HESS J1731-347-like stars, the crust is about 50 percent thicker, and its mass and moment of inertia are several times larger than with the Maxwell construction or standard crust models.","The crust-core transition density rises to about 0.26 to 0.30 fm$^{-3}$, roughly twice saturation density, instead of the usual 0.03 to 0.08 fm$^{-3}$.","The mass-radius relation remains compatible with HESS J1731-347, NICER measurements of PSR J0740+6620, PSR J0030+0451 and PSR J0437-4715, GW170817 constraints, and the maximum-mass lower bound, with the $C_\\sigma^2 = 14$ model in best agreement with the most massive pulsar constraint.","In the mixed-phase region the speed of sound stays finite and positive, avoiding the density discontinuity of the Maxwell-construction EOS, and it satisfies the general relativistic kinetic-theory bound even though it exceeds $c/\\sqrt{3}$.","The enlarged crustal mass and moment of inertia would alter the star's thermal evolution and its rotational and glitch properties."],"supporting_citations":[{"why":"The CCT model and its Maxwell-construction EOS; the paper re-solves this same model with Gibbs conditions and compares the results.","marker":"[11]"},{"why":"Provides the thermodynamic argument that two conserved charges require full Gibbs conditions rather than a Maxwell construction, and is the source of the finite-size effect caveat.","marker":"[14]"},{"why":"The SLy4 equation of state used to model the crust below the joining density.","marker":"[15]"},{"why":"Establishes the mechanism by which the low-density phase drops protons when the neutron-lepton chemical-potential difference is too small.","marker":"[16]"},{"why":"Earlier analysis of a two-phase system slightly above saturation density; the present work extends it down to inner-crust densities.","marker":"[17]"},{"why":"Original discovery and mass-radius measurement of HESS J1731-347 that the models must reproduce.","marker":"[1]"},{"why":"PSR J0740+6620 mass measurement that sets the lower bound on maximum mass the EOS must satisfy.","marker":"[3]"},{"why":"General analysis of phase-transition softening and speed-of-sound limits used to justify the mixed-phase EOS.","marker":"[5]"}],"fun_headline_variants":["Gibbs conditions push HESS J1731-347 crust boundary to 2n0","Mixed phase in HESS J1731-347 core thickens crust 50%","Correct phase rules multiply HESS J1731-347 crust mass","HESS J1731-347's crust-core transition at double saturation density","Thermodynamic consistency gives HESS J1731-347 a thick crust"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole thick-crust picture rests on neglecting the surface tension and electric (Coulomb) energy of the tiny mixed-phase structures; if those effects are sizable, the crust's mass, thickness, and moment of inertia would change even though the density where the mixed phase first appears would not.","fun_headline_variants_meta":{"raw":{"variants":["Gibbs conditions push HESS J1731-347 crust boundary to 2n0","Mixed phase in HESS J1731-347 core thickens crust 50%","Correct phase rules multiply HESS J1731-347 crust mass","HESS J1731-347's crust-core transition at double saturation density","Thermodynamic consistency gives HESS J1731-347 a thick crust"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001349,"raw_usage":{"total_tokens":5463,"prompt_tokens":912,"completion_tokens":4551,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":4447}},"tokens_in":528,"tokens_out":4551,"duration_ms":37758,"temperature":1.0,"reasoning_tokens":4447,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:52:29.511201+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A calculation of the finite-size structured mixed phase (nuclear pasta) for the CCT equation of state over the density range 0.0025 to 0.30 fm$^{-3}$ would settle the quantitative claim: if inclusion of surface and Coulomb energies shifts the crust-core boundary below roughly 0.1 fm$^{-3}$, or changes the crustal moment-of-inertia fraction outside the predicted several-times-larger range, the thick-crust conclusion fails. Alternatively, a precise glitch measurement on a star with the mass of HESS J1731-347 that implies a crustal moment of inertia comparable to thin-crust models would contradict the prediction.","supporting_citations":[{"cited_title":"Miyatsu, M","cited_arxiv_id":null,"evidence_quote":"The CCT model and its Maxwell-construction EOS; the paper re-solves this same model with Gibbs conditions and compares the results."},{"cited_title":"Lope-Oter and A","cited_arxiv_id":null,"evidence_quote":"Provides the thermodynamic argument that two conserved charges require full Gibbs conditions rather than a Maxwell construction, and is the source of the finite-size effect caveat."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The SLy4 equation of state used to model the crust below the joining density."},{"cited_title":"Kubis, Acta Phys","cited_arxiv_id":null,"evidence_quote":"Earlier analysis of a two-phase system slightly above saturation density; the present work extends it down to inner-crust densities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Original discovery and mass-radius measurement of HESS J1731-347 that the models must reproduce."},{"cited_title":"Veselsky, P","cited_arxiv_id":null,"evidence_quote":"PSR J0740+6620 mass measurement that sets the lower bound on maximum mass the EOS must satisfy."},{"cited_title":"Laskos-Patkos, P","cited_arxiv_id":null,"evidence_quote":"General analysis of phase-transition softening and speed-of-sound limits used to justify the mixed-phase EOS."}],"review_version":1}