{"id":"1c33803e-b70d-47a8-9dcb-c3fe611fe5bb","arxiv_id":"2505.19984","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For the same Skyrme functional BSk24, the simplified liquid-drop model matches the more detailed Thomas-Fermi method for finite-temperature thermodynamics, but misses the neutron skin and overestimates proton radii.","lead":"This paper compares a simple liquid-drop model with a more detailed semiclassical calculation of the hot inner crust of neutron stars, and finds they agree on thermodynamics but differ on density profiles. The result helps simulations of neutron star mergers and supernovae decide which approximation is good enough.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Composition agreement between CLDM and TETF may be partly an artifact of the common single-radius surface-energy extraction: the same spurious bulk contamination enters the fit and the validation.","rationale":"The paper is a careful, well-qualified benchmark: it uses one functional consistently, compares three surface-fitting protocols, discloses the single-radius limitation in Appendix B, and does not overstate its conclusions. The central claim most in need of scrutiny is not the thermodynamic agreement, which appears robust because it is not very sensitive to finite-size details, but the composition statement. For that statement, the surface energy is the controlling input. The weakest point is the definition of the cluster radius and density in Appendix A: n_i is equated with the central total density and n_Np is fixed by charge conservation through Eq. (A5). For a neutron-rich cluster with a skin, the sharp-interface bulk term V_N F_B(n_Np,n_Nn) is not the actual bulk energy of the diffuse TETF solution, and the missing piece is silently assigned to E_surf+curv. The authors explicitly flag this ('may include some spurious bulk contribution', Appendix B), but no estimate of its magnitude is given. Because the same extraction enters both the CLDM fit (i) and the TETF-side surface energy plotted in Fig. A1, the apparent agreement of surface energies is not independent evidence; it can reflect that the CLDM was fitted to the same contaminated decomposition. The composition comparison in Fig. 3 is a nontrivial finite-temperature test, and the similarity of fit (i) and fit (ii) results, plus footnote 8's claim of a worse skin-aware fit, independently supports the authors' position. Still, without a quantitative bound on the spurious contribution, the 'reasonable reproduction of composition' claim is conditional: it rests on the assumption that the spurious bulk/surface separation is small everywhere in the crust. This is exactly the reader's weakest_assumption, and the proposed check would settle it. It does not overturn the paper; it sharpens the condition under which the central composition claim should be accepted. The current CONDITIONAL verdict is therefore appropriate.","tokens_in":20194,"tokens_out":6424,"duration_ms":80000,"concrete_test":"Take a converged TETF cell (e.g., nB=0.04 fm^-3, T=1 MeV) and compute the bulk part of Eq. (A6) two ways: with the sharp-interface densities of Eqs. (A3)-(A5) versus the actual profile integral of F_B(n_n(r), n_p(r)) over the same r_N. The difference is the spurious bulk contribution attributed to E_surf. If this difference exceeds about 10% of the extracted E_surf at any crust density, the fit (i) parameters are contaminated; one should then repeat the Fig. 3 comparison with a two-radius (skin-aware) extraction to see whether the reported Z and Ntot agreement survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that a CLDM can reproduce the TETF composition relies on surface parameters fitted to 'ETF calculations in the medium' (fit i). The extraction in Appendix A, Eqs. (A1)-(A5), collapses the diffuse TETF cluster to one sharp radius r_N with n_Nn = n_i - n_Np, where n_i is the central total density. Because the ETF cluster has a neutron skin, Eqs. (A3)-(A5) do not give the true bulk contribution, so the residual E_surf+curv 'may include some spurious bulk contribution,' as the authors concede in Appendix B. The same contaminated decomposition is then used in Eq. (A6) for the TETF surface energy compared in Fig. A1. A fit to the contaminated quantity is expected to reproduce it, so the excellent T=0 agreement in Fig. A1 and part of the composition agreement in Fig. 3 could be an artifact of the common single-radius mapping rather than evidence that the fitted sharp-interface free energy captures the TETF finite-size effects. The authors report in footnote 8 that including a skin in the fit yields a higher chi2, and fit (ii) (ETF mass table) is not exposed to the medium extraction and gives broadly similar composition trends; these points soften, but do not quantify away, the concern. What is missing is a bound on the spurious bulk/surface separation, and a demonstration that a skin-aware fit changes Z and Ntot by less than the visual spread in Fig. 3.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares two semiclassical descriptions of the inner crust of a non-accreting neutron star at finite temperature: the temperature-dependent extended Thomas-Fermi (TETF) method and the compressible liquid-drop model (CLDM), both implemented with the same BSk24 functional. The CLDM surface and curvature parameters are fitted in three ways: to zero-temperature ETF calculations of clusters in the medium (protocol i), to a full ETF mass table (protocol ii), and to experimental AME2020 masses (protocol iii). The authors report that the CLDM reproduces the TETF thermodynamic quantities closely, that the TETF composition is best approximated when the CLDM surface energy is fitted to ETF-based data (protocols i and ii), and that the CLDM's neglect of the neutron skin leads to an overestimation of proton radii. They also provide the first second-order TETF equation of state for the inner crust and release the numerical data as supplementary material.","tokens_in":20538,"tokens_out":10449,"duration_ms":107046,"significance":"The paper addresses a question of practical importance: whether the fast CLDM can be substituted for the more microscopic TETF when constructing finite-temperature equations of state for supernova and merger simulations. Its design is clean: using one functional (BSk24) in both frameworks isolates the effect of the finite-size treatment, and the three fitting protocols bracket the sensitivity to the surface-energy input. The main deliverables—a first TETF-based hot inner-crust EoS, benchmark comparisons, and public data—are useful to the community. If the composition agreement survives a quantitative and skin-aware robustness analysis, the conclusions would justify the use of CLDM-based tables for thermodynamic quantities while reserving TETF for density-profile-sensitive microphysics. The present version, however, does not yet provide the quantitative support needed to turn the visual agreement into a firm benchmark statement.","major_comments":[{"comment":"The central claim of agreement between CLDM and TETF for the thermodynamic properties and composition rests entirely on visual inspection of overlaid curves. No quantitative deviations are reported for F/A_tot, P, μ_n, Y_e, r_WS, Z, or N_tot. Please add explicit comparison metrics (e.g., mean and maximum absolute or relative differences over a stated density range, for each temperature and for each of the three fit protocols). This is particularly important for the composition plots in Fig. 3, where the CLDM underestimates nucleon numbers for roughly 0.01 fm^-3 < n_B < 0.05 fm^-3 and overestimates them near the crust-core transition; without numbers the statement that the agreement is 'reasonable' cannot be assessed.","section":"Sec. 3.1, Figs 1-3"},{"comment":"The favorable T=0 comparison in Fig. A1 is partly built into the fitting procedure. Protocol (i) extracts the surface energy from ETF calculations by mapping the diffuse cluster onto a single sharp radius r_N with n_i taken as the central total density, n_Np from charge conservation, and n_Nn = n_i - n_Np (Eqs. A3-A5). The same mapping is then used in Eq. (A6) to define the TETF surface energy that is compared with the CLDM in Fig. A1. Because the ETF/TETF clusters possess a neutron skin, the bulk subtraction is not the true bulk contribution, and the residual 'may include some spurious bulk contribution,' as acknowledged in Appendix B. A fit to this contaminated residual is expected to reproduce it, so the agreement does not by itself demonstrate that the fitted sharp-interface free energy captures the finite-size effects. Please provide a quantitative estimate of the spurious bulk contribution—for example, by repeating the extraction with a two-radius, skin-aware decomposition and reporting how the fitted parameters in Table 1 and the resulting Z and N_tot change. Footnote 8 reports only that a skin-including fit yields a higher χ2, with no value, no degrees of freedom, and no parameter variation; as it stands it does not rule out a sizable bias.","section":"Appendix A, Eqs. (A1)-(A5); Appendix B, Eq. (A6)"},{"comment":"The best-fit surface parameters are quoted to six significant figures without uncertainties or fit-quality indicators. Moreover, for protocols (ii) and (iii) the parameter p is fixed to 3 rather than fitted, whereas for protocol (i) it is fitted; the table and text do not make this distinction explicit. Please state which parameters were floated in each fit, report their uncertainties (or covariance) and the χ2 or rms residual, and justify the choice p=3 with a sensitivity analysis. This is necessary because the composition curves in Fig. 3 depend directly on these values.","section":"Table 1"}],"minor_comments":[{"comment":"The sentence 'The third term on the right-hand side of Equation (15), -uFg, accounts for the excluded volume' is imprecise: Eq. (15) has three terms, and the excluded-volume correction enters through the factor (1-u) in the second term, not as a separate third term. Please rephrase.","section":"Sec. 2.2, after Eq. (15)"},{"comment":"Several figure labels and captions contain typographical errors or garbled text, e.g., 'CLDM(fit ETF in the medi m)' in Fig. 1, 'CLDM(fit ETF ass table)' and 'CLDM(fit ETF in the ediu )' in Fig. A1, and 'E_su f' in Fig. A1. Please correct these labels.","section":"Figs 1, 3, A1"},{"comment":"The phrase 'all data used to produce the figures of this work are avaible as supplementary material' contains a typo ('avaible' should be 'available').","section":"Sec. 3.1"},{"comment":"The fitting routine is referred to as 'scipy.curvfit'; the correct name is scipy.optimize.curve_fit.","section":"Appendix A, last paragraph"},{"comment":"The sentence 'The nuclear energetics is described employing a CLDM model approach' is grammatically awkward and redundantly says 'CLDM model'; please rephrase.","section":"Sec. 2.2, opening sentence"},{"comment":"In the discussion of the neutron chemical potential, the sentence 'due to the degeneracy of the neutron gas, which is, however, lifted by the finite-temperature effect at very low densities leading to the noticeable smaller µn at T = 2 MeV' is a run-on; consider splitting it into two sentences for clarity.","section":"Sec. 3.1, bottom-panel discussion"}],"recommendation":"major_revision","confidential_remarks":"The main technical concern is the possible contamination of the protocol-(i) surface energy by the single-radius decomposition. The independent protocol (ii), fit to the ETF mass table, gives similar composition trends, which is evidence that the qualitative conclusion is robust; however, the manuscript should make this quantitative rather than rely on the reader's eye. The paper is within scope for Universe and the data release is a plus."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about arXiv:2505.19984. First, it delivers the first systematic finite-temperature comparison of CLDM and TETF for the neutron star inner crust using the same BSk24 functional, and the first second-order TETF inner-crust EoS at finite T. Second, the agreement it reports between CLDM and TETF thermodynamics is credible; the composition agreement is real but partly dependent on how the CLDM surface parameters are fitted.\n\nThe setup is clean: same functional in both methods, three different surface-fitting protocols (ETF in medium, ETF mass table, AME2020), and the authors openly discuss the neutron-skin limitation. The thermodynamic quantities—free energy, pressure, chemical potential, electron fraction—match well between CLDM and TETF, which is the practically important result for people building finite-T EoS tables. The composition comparison shows CLDM can reproduce TETF Z and Ntot trends, with the ETF-based fits doing better than the AME fit. That is a useful, honest benchmark.\n\nThe soft spots are real but disclosed. The strongest concern, which the stress-test note picks up, is that the surface-energy extraction in Appendix A uses a single-radius prescription that collapses the diffuse TETF profile with a neutron skin into one sharp radius. As the authors admit in Appendix B, the extracted surface energy 'may include some spurious bulk contribution.' Since the same extraction is used both to fit the CLDM surface parameters (fit i) and to define the TETF surface energy compared in Figure A1, that comparison is partly circular. The paper softens this by showing fit (ii) (ETF mass table, not exposed to the medium extraction) gives broadly similar composition trends, and footnote 8 says a skin-aware fit gives worse chi-squared. But there is no quantitative bound on the spurious bulk/surface separation; a skeptic would want to know how much the composition results shift if the extraction allowed a skin. Another minor issue: no numerical measures of agreement, only visual inspection, and the Table 1 surface parameters have no uncertainties.\n\nOverall, the central claims hold up. The paper is careful, the limitations are in the text, and the first-TETF-EoS plus systematic CLDM comparison is a contribution people in the EoS community will want. It deserves a serious referee, though the referee should push for a quantitative estimate of the surface-extraction bias. I would cite it for benchmark purposes.","headline":"Useful, honest CLDM vs TETF benchmark for the hot inner crust; the composition agreement depends partly on a circular surface-energy fit, but the paper discloses it and the thermodynamic results are solid.","tokens_in":21024,"tokens_out":2038,"would_cite":true,"duration_ms":21083,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["26.60.-c","21.65.Mn"],"model":"deepseek-v4-flash","headline":"At finite temperature, the fast liquid-drop model matches the extended Thomas-Fermi crust calculation—except for proton radii.","keywords":["neutron stars","inner crust","finite temperature","equation of state","compressible liquid drop model","extended Thomas-Fermi","surface energy","BSk24 functional"],"falsifier":"A finite-temperature Hartree-Fock-Bogoliubov calculation with the BSk24 functional, or a TETF run whose density profiles allow a separate neutron radius (a two-radius parametrization), would settle the claimed cause of the proton-radius overestimation: if the overestimation persists when a neutron skin is allowed, the single-radius geometry is not the explanation. A second, cheaper check is to redo the Appendix A surface-energy extraction with a two-radius prescription and see whether the fitted surface parameters move by more than their fit uncertainties, which would confirm that the spurious bulk contribution the authors warn about is materially affecting the CLDM fit.","tokens_in":20018,"feed_emoji":"🌟","tokens_out":12799,"duration_ms":115981,"temperature":0.7,"pith_summary":"This paper tests whether the computationally cheap compressible liquid-drop model (CLDM) can stand in for the more microscopic, temperature-dependent extended Thomas-Fermi (TETF) method when describing the inner crust of a neutron star at finite temperatures up to about 2 MeV. Building both approaches on the same nuclear functional (BSk24), it finds that the CLDM closely tracks the TETF predictions for the thermodynamic quantities—free energy, pressure, and neutron chemical potential—and reproduces the composition trends (proton and neutron numbers per Wigner-Seitz cell), especially when the CLDM surface energy is fitted to ETF calculations rather than to experimental masses. The systematic failure is geometric: the CLDM represents each cluster by a single radius with no neutron skin, so it places protons too far out and overestimates their radii, a defect that matters for transport and elastic properties that depend on the charge distribution. If correct, the benchmark tells builders of supernova and merger equations of state where fast CLDM-based tables can be trusted and where a TETF-level treatment is needed.","feed_headline":"Fast liquid-drop model matches Thomas-Fermi, except proton radii","feed_subtitle":"The liquid-drop equation of state stays reliable up to 2 MeV, but neutron-skin microphysics needs the fuller model.","key_machinery":"The argument is carried by the contrast between two finite-size treatments sharing one functional. The TETF method is a semi-classical approximation built from the second-order Wigner-Kirkwood (gradient) expansion of the free energy, with neutron and proton density profiles parametrized by smooth soft-damping shapes and minimized at each density and temperature. The CLDM condenses each cluster into a uniform sphere of radius $r_N$ and internal density $n_i$, splitting the free energy into bulk, Coulomb, and surface-plus-curvature parts, with surface and curvature tensions $\\sigma_s$ and $\\sigma_c$ depending on cluster asymmetry; the load-bearing step is how the surface parameters are obtained, and the paper compares three fits: to ETF calculations in the medium at $T = 0$, to an ETF mass table from the proton to the neutron drip line, and to the AME2020 experimental masses. The single-radius closure relation—baryon conservation fixing $r_N$ and charge conservation fixing the cluster proton density—is what makes the CLDM fast, and also what removes the neutron skin; all three fits ignore the skin, and a fit that included it raised the $\\chi^2$ on the zero-temperature fits, which is why the paper keeps the one-radius prescription.","core_discovery":"On the paper's own terms, the central discovery is that the one-component-plasma, Wigner-Seitz treatment of the inner crust in beta equilibrium is largely independent of which of the two finite-size descriptions is used. At temperatures of 1–2 MeV and across inner-crust baryon densities, the CLDM and TETF results for free energy per nucleon, pressure, and neutron chemical potential nearly overlap, and temperature moves them the same way; the electron fraction and Wigner-Seitz radius show only small, low-density deviations. For the composition, both methods predict that the proton number $Z$ grows with temperature at high densities and falls with temperature at low densities, and the CLDM reproduces the TETF values best when its surface and curvature parameters are fitted to zero-temperature ETF calculations in the medium or to a full ETF mass table—fits that include highly isospin-asymmetric matter—rather than to the AME2020 experimental masses alone, though even that fit is acceptable. The qualitatively different outcome is the nucleon density profiles: CLDM-deduced profiles are step-like with a single radius, and because the neutron skin is absent the proton radius is overestimated relative to TETF, with the gap growing toward the crust-core transition. The paper also delivers the first inner-crust equation of state and composition from the second-order TETF method with the BSk24 functional, which is the benchmark the CLDM is measured against.","pith_inferences":["A minimal upgrade of the CLDM to two radii (cluster radius plus a fitted neutron-skin thickness), fitted on the same ETF benchmarks at finite temperature, might recover TETF proton radii while keeping the model fast; the authors rejected a neutron-skin fit only because it raised the chi-square on the zero-temperature fits they used.","The systematic tendency of the AME2020-fitted CLDM to predict lower proton and neutron numbers than TETF suggests that equation-of-state tables calibrated to terrestrial masses alone carry a bias toward smaller, less neutron-rich clusters in the inner crust; repeating this comparison with functionals of different symmetry-energy slope would show whether the bias scales with the skin size.","Since the authors set the CLDM surface tension to its zero-temperature value at all temperatures studied, the temperature dependence of the crust composition they report is carried by bulk and Coulomb terms; a dedicated fit of the surface tension's temperature dependence from the TETF profiles would test whether surface melting matters below 2 MeV."],"forward_implications":["Fast CLDM-based finite-temperature equation-of-state tables remain trustworthy for crust thermodynamics, since pressure, energy, and chemical potentials track the TETF benchmark closely up to about 2 MeV.","CLDM composition accuracy improves when the surface fit includes extremely isospin-asymmetric matter (ETF calculations in the medium or ETF mass tables), information absent from terrestrial mass data.","Microphysics that depends on the proton density distribution—electron conductivity, neutrino transport, elastic properties—should be taken from TETF-class calculations, because the one-radius CLDM overestimates proton radii.","A new second-order TETF inner-crust equation of state and composition table for BSk24 is now available as a benchmark, including first finite-temperature TETF results of this kind."],"supporting_citations":[{"why":"Defines the BSk24 functional used identically in both TETF and CLDM, so the comparison isolates the finite-size treatment.","marker":"[28]"},{"why":"Supplies the temperature-dependent extended Thomas-Fermi formalism (second-order Wigner-Kirkwood expansion) that produces the TETF results.","marker":"[20]"},{"why":"Provides the finite-temperature CLDM treatment (bulk free energy via the meta-model, excluded volume) that this work follows.","marker":"[7]"},{"why":"The earlier finite-temperature comparison and source of the ETF mass-table surface fit; this paper extends it systematically.","marker":"[16]"},{"why":"Provides the zero-temperature ETF calculations in the medium used to extract the surface energy for fit protocol (i).","marker":"[29]"},{"why":"Supplies the procedure for parametrizing the CLDM surface energy from ETF calculations, adapted for the medium fit.","marker":"[64]"},{"why":"The AME2020 experimental mass table used for the third surface-parameter fit.","marker":"[65]"},{"why":"The zero-temperature CLDM-versus-ETF comparison whose finite-size conclusions this work extends to finite temperature.","marker":"[27]"}],"fun_headline_variants":["Liquid-drop and Thomas-Fermi inner crusts agree, but proton radii differ","Neutron star crust models match until neutron skin skews proton radii","Temperature-dependent crust: CLDM and TETF agree, except for proton size","Inner crust thermodynamics robust, but proton radius needs full model","Hot neutron star crust: liquid drop matches Thomas-Fermi except proton size"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison stands or falls on the assumption that a one-radius cluster with no neutron skin is an adequate representation of the ETF cluster, with the surface energy extracted under that assumption; the authors state this extraction may include a spurious bulk contribution, and if that spurious term is significant the fitted surface parameters—and hence the CLDM composition agreement—would be biased.","fun_headline_variants_meta":{"raw":{"variants":["Liquid-drop and Thomas-Fermi inner crusts agree, but proton radii differ","Neutron star crust models match until neutron skin skews proton radii","Temperature-dependent crust: CLDM and TETF agree, except for proton size","Inner crust thermodynamics robust, but proton radius needs full model","Hot neutron star crust: liquid drop matches Thomas-Fermi except proton size"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000239,"raw_usage":{"total_tokens":1518,"prompt_tokens":949,"completion_tokens":569,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":471}},"tokens_in":565,"tokens_out":569,"duration_ms":6758,"temperature":1.0,"reasoning_tokens":471,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:02:52.253582+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A finite-temperature Hartree-Fock-Bogoliubov calculation with the BSk24 functional, or a TETF run whose density profiles allow a separate neutron radius (a two-radius parametrization), would settle the claimed cause of the proton-radius overestimation: if the overestimation persists when a neutron skin is allowed, the single-radius geometry is not the explanation. A second, cheaper check is to redo the Appendix A surface-energy extraction with a two-radius prescription and see whether the fitted surface parameters move by more than their fit uncertainties, which would confirm that the spurious bulk contribution the authors warn about is materially affecting the CLDM fit.","supporting_citations":[{"cited_title":"Further explorations of Skyrme-Hartree-Fock-Bogoliubov mass formulas","cited_arxiv_id":null,"evidence_quote":"Defines the BSk24 functional used identically in both TETF and CLDM, so the comparison isolates the finite-size treatment."},{"cited_title":"Equation of state of stellar nuclear matter in the temperature-dependent extended Thomas–Fermi formalism","cited_arxiv_id":null,"evidence_quote":"Supplies the temperature-dependent extended Thomas-Fermi formalism (second-order Wigner-Kirkwood expansion) that produces the TETF results."},{"cited_title":"Parametrization of the surface energy in the ETF approximation","cited_arxiv_id":null,"evidence_quote":"Supplies the procedure for parametrizing the CLDM surface energy from ETF calculations, adapted for the medium fit."}],"review_version":1}