{"id":"9ae09b88-9a04-424f-be36-062a9bc8dfb8","arxiv_id":"2505.18309","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Nuclear pasta onset in neutron star crusts occurs at a quasi-universal filling fraction u_sp ≈ 0.13-0.15, explained by curvature corrections in the liquid-drop model.","lead":"This paper computes where non-spherical nuclear 'pasta' phases appear in neutron star crusts, using two different models. It finds the onset is nearly universal at a filling fraction of about 0.13 to 0.15, and explains this with curvature corrections.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The CLDM explanation hinges on σc/σs extrapolated to y_p≈0.03–0.15 via Eq. (17); fit uncertainty is unquantified and no first-principles check is given, so the u_sp≈0.14 match with ETF is not yet established as causal.","rationale":"The central claim has two parts: the ETF quasi-universality (Sec. 2) and the CLDM curvature explanation (Sec. 3). The ETF part is supported by three BSk functionals plus SLy4 and is consistent with earlier TF studies that found u_sp≈1/8 [16–19], so I do not see a compelling reason to doubt it. The CLDM part is where the causal assertion lives: the paper claims the agreement with ETF occurs because of the surface curvature correction. That assertion is entirely controlled by σc/σs at y_p≈0.03–0.15, which is not measured or first-principles computed in this work. The authors' own caveat in the Conclusions is an explicit admission of missing support, and the review rules require me to flag it. The internal cross-checks (reproducing Ref. [4] without curvature, matching Ref. [38] for BSk24/exp masses) show the criterion is correctly implemented, but they do not validate the extrapolation of Eq. (17) to the pasta region. The lack of code/data is a reproducibility issue, not a correctness issue, and the reader already conditioned on it; it is secondary to the scientific claim. The proposed direct computation of σs and σc in semi-infinite asymmetric matter, followed by recomputation of u_sp, would settle whether the apparent agreement in Fig. 6 is causal or coincidental. The paper should remain CONDITIONAL pending this check and/or a quantified uncertainty budget for u_sp derived from the fit covariance.","tokens_in":12516,"tokens_out":15361,"duration_ms":135568,"concrete_test":"Compute σs and σc for BSk24 (and, if feasible, one additional functional such as SLy4) directly from the ETF functional in semi-infinite nuclear matter at the neutron chemical potential and proton fraction of the inner crust (y_p≈0.03–0.15), following Refs. [13,28,31,32], instead of using the finite-nucleus fits of Eq. (17). Then re-evaluate the instability threshold u_sp from Eqs. (12)–(13) with this first-principles σc/σs and the same ETF inputs a(u) and y_p(u). If the recomputed u_sp remains within roughly 0.13–0.15, the curvature explanation survives; if it moves outside this band by more than 0.02, the agreement in Fig. 6 is an artifact of the fitted σc/σs extrapolation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 3.3 the CLDM instability criterion (Eqs. 12–13) is used to explain the ETF onset u_sp≈0.13–0.14. The only term that distinguishes the curvature-corrected prediction from the no-curvature one is Xσ/a = (g_c/g_s)(σc/σs)/a (Eq. 11). The ratio σc/σs is taken from Eq. (17), an ansatz fitted either to ETF mass tables (Ref. [36]) or to experimental atomic masses (Ref. [33]). Pasta-relevant proton fractions are y_p≈0.03–0.15 (yellow band, Fig. 5). The fitted nuclei, even in the ETF mass table, have proton fractions well above this band; Eq. (17) is thus an uncontrolled extrapolation into a regime where a neutron gas surrounds the cluster and the surface/curvature thermodynamics may change qualitatively. The authors themselves concede in the Conclusions that σs and σc 'remain very uncertain at the extreme isospin asymmetries prevailing in the deepest layers of neutron-star crusts.' No uncertainty on σc/σs is propagated into the u_sp predictions of Fig. 6, and no comparison with directly computed semi-infinite-matter surface tensions is provided. Since the no-curvature prediction is u_sp≈0.215 (dotted lines) and the curvature shift to ≈0.14 scales linearly with σc/σs, a modest error in the fitted ratio—for example, a factor of two at y_p≈0.05—could move u_sp by ~0.05 and erase the agreement. The ETF quasi-universality itself would survive, but the causal claim that the curvature correction is responsible for u_sp≈1/8 would not be established. The fact that the two fits (ETF vs experimental masses) give similar u_sp does not bound the extrapolation error, because both share the same ansatz and the same scarcity of very neutron-rich data.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates the conditions for nuclear pasta formation in neutron-star inner crusts. Using the second-order extended Thomas-Fermi (ETF) method with three Brussels-Montreal Skyrme functionals (BSk22, BSk24, BSk25) and a soft-damping density-profile parametrization, it computes the equilibrium pasta sequence and finds that the filling fractions at shape transitions are quasi-universal across functionals, with the sphere-to-pasta onset at u_sp≈0.13–0.14. The paper then derives a compressible liquid-drop model (CLDM) shape-instability criterion that includes curvature and neutron-skin corrections, and shows that with curvature tensions fitted either to ETF mass tables or to experimental atomic masses, the CLDM criterion predicts u_sp≈0.14–0.15, close to the ETF result. The paper proposes this criterion as a fast estimator of pasta abundance and as a guide for ETF searches.","tokens_in":12962,"tokens_out":11283,"duration_ms":96627,"significance":"If the claims hold, the paper provides a simple, falsifiable prediction that pasta appears near a filling fraction of 1/8 in neutron-star crusts, and an analytic criterion to estimate pasta abundance without costly ETF calculations. The ETF results are transparent and cover several functionals, and the use of experimental-mass-based surface tensions for BSk24 and SLy4 provides a partial check against circularity. However, the explanatory power of the CLDM analysis depends on the extrapolation of surface and curvature tensions to extreme isospin asymmetries, which the authors themselves flag as uncertain, and on a consistent definition of the filling fraction between the ETF and CLDM frameworks. The paper does not ship code or deposited data, but the numerical results are described in enough detail to be reproduced in principle.","major_comments":[{"comment":"The central explanatory claim that the curvature correction shifts u_sp from about 0.215 to about 0.14–0.15 relies entirely on the ratio sigma_c/sigma_s evaluated from Eq. (17) at y_p ≈ 0.03–0.15. This is an extrapolation from fits to finite nuclei, either from ETF mass tables or experimental atomic masses, into a regime that the authors themselves describe in the Conclusions as 'very uncertain.' The paper neither propagates fit uncertainties into Fig. 6 nor compares with direct semi-infinite-matter calculations of sigma_s and sigma_c for the same functionals at the relevant isospin asymmetries. Because the curvature shift scales linearly with sigma_c/sigma_s (Eq. 11), a factor-of-two error in the extrapolated ratio at y_p ≈ 0.05 could move u_sp by roughly 0.05 and potentially erase the agreement. Please quantify this sensitivity and provide an independent check.","section":"Sec. 3.3, Eq. (17), and Conclusions"},{"comment":"The ETF filling fraction u is defined from the proton density profile in Eq. (1), whereas the CLDM instability criterion used in Sec. 3.2 assumes u is the geometric volume fraction of the dense total-nucleon phase, as implied by the bulk term E_bulk(ni, yp, non, nop, u) and the geometry functions in Table 1. In neutron-rich crustal matter with y_p ≈ 0.03–0.15, the neutron skin makes these two definitions appreciably different. If the ETF proton-based u is inserted into Eqs. (12)–(13), the comparison in Fig. 6 may mix two different quantities. Please justify that the proton-based u is the appropriate variable for the CLDM geometry functions, or repeat the comparison using a total-baryon-density filling fraction.","section":"Sec. 2, Eq. (1), versus Sec. 3.1, Table 1"}],"minor_comments":[{"comment":"The row for lasagna reads '2 0 πu^2(1−u)^2/6', which is ambiguous; the entries for g_s, g_c, and w should be clearly separated, presumably g_s = 2, g_c = 0, and w = πu^2(1−u)^2/6.","section":"Table 1"},{"comment":"The yellow band labeled 'pasta' does not state the corresponding y_p interval; please give the numerical range, e.g., y_p ≈ 0.03–0.15, in the caption.","section":"Fig. 5"},{"comment":"The abstract quotes u_sp ≈ 0.13–0.15 while Sec. 2 reports 0.13–0.14 for the ETF calculations and Sec. 3.3 reports 0.14–0.15 for the CLDM criterion; stating the separate ranges explicitly in the abstract would avoid confusion.","section":"Abstract and Secs. 2, 3.3"},{"comment":"The data availability statement first says 'This manuscript has no associated data' and then says the data are available from N.N.S. upon reasonable request; please harmonize these two statements.","section":"Declarations, Data availability"}],"recommendation":"major_revision","confidential_remarks":"The proton-based versus total-density filling fraction issue is the point I would most want resolved before publication; if the authors can show that the two definitions give similar u at the transitions, the causal claim about the curvature correction becomes much stronger. The manuscript is well within the scope of the journal and the ETF results are likely to be of lasting value even if the CLDM explanation needs refinement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real result here is the ETF calculation: across BSk22/24/25, the onset filling fraction for pasta is quasi-universal at u_sp ≈ 0.13–0.14, and the density ranges vary with symmetry energy in a way that is clearly explained by the u(¯n) formula. That is a clean, internally consistent numerical finding, and it settles the earlier tension between the TF results near 1/8 and the CLDM spread of 0.08–0.28. The paper also gives a neat generalized CLDM instability criterion (Eqs. 12–13), which reduces to the known results without curvature and is numerically checked against full CLDM calculations from Ref. [38]. The demonstration that the curvature correction is what shifts u_sp from ~0.215 to ~0.14 is plausible and well argued.\n\nThe soft spot is the causal claim in Sec. 3.3. The CLDM agreement depends on the ratio σc/σs at proton fractions 0.03–0.15, but Eq. (17) is fitted to nuclei far less neutron-rich. The authors acknowledge this in the Conclusions, and the stress-test note is right that no uncertainty is propagated. However, the concern should be kept in proportion: the CLDM is used as an explanatory tool, not as the primary evidence. The ETF quasi-universality stands on its own, and the fact that two independent fits (ETF masses and experimental masses) give similar σc/σs in the pasta band is at least suggestive. The remaining transitions (lasagna, bucatini) degrade, and the paper says so. That is honest. The bigger practical weakness is the absence of code or data; the paper says data are available on request, but that limits reproducibility.\n\nThis is a useful paper for people working on neutron-star crust structure and pasta transport/cooling. It deserves a serious referee. I would recommend acceptance after the authors provide the ETF filling-fraction data or code and add a quantitative caveat—even a rough one—on how much σc/σs would have to change to erase the CLDM match. The central ETF result holds up.","headline":"Solid ETF result on quasi-universal pasta filling fractions, but the CLDM curvature explanation is not yet causal because of unquantified extrapolation of surface tensions to extreme isospin.","tokens_in":13531,"tokens_out":558,"would_cite":true,"duration_ms":6449,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Nuclear pasta in neutron stars appears when clusters fill about one-eighth of their cell, and the liquid-drop curvature correction is why.","keywords":["neutron star inner crust","nuclear pasta","filling fraction","extended Thomas-Fermi","compressible liquid drop model","curvature correction","generalized Skyrme functionals","symmetry energy"],"falsifier":"A decisive check would be to compute the sphere-to-pasta onset for the same three functionals with a method that does not rely on the liquid-drop geometry or the fitted density-profile ansatz—for instance a full quantum band-structure or molecular-dynamics calculation at inner-crust densities—and see whether the onset still sits at $u\\approx0.13$–$0.15$; alternatively, computing $\\sigma_c/\\sigma_s$ from first principles for semi-infinite neutron-rich matter and feeding it into Eqs. (12)–(13) would test whether the curvature explanation, rather than the calibration, is responsible.","tokens_in":12349,"feed_emoji":"🍝","tokens_out":11514,"duration_ms":108066,"temperature":0.7,"pith_summary":"Deep in a neutron star's inner crust, nuclei are packed so tightly that theory predicts they deform into rods, slabs, tubes, and bubbles known collectively as nuclear pasta. This paper tries to establish that the filling fraction at which spheres first give way to pasta is quasi-universal, around $u_\\mathrm{sp}\\approx0.13$–$0.15$, across three precision-fitted generalized Skyrme functionals, and that the low value is caused by the curvature correction to the nuclear surface energy in a liquid-drop description. The result matters because pasta formation becomes a simple geometric condition: once the local filling fraction crosses about one-eighth, the spherical shape should be abandoned, and a cheap algebraic criterion can estimate how much pasta a neutron star contains. The ETF calculations place the onset at $u_\\mathrm{sp}\\approx0.13$–$0.14$, while the curvature-corrected liquid-drop stability analysis gives $u_\\mathrm{sp}\\approx0.14$–$0.15$, with later transitions near 0.26–0.27, 0.48–0.51, and 0.66–0.69.","feed_headline":"Pasta forms in neutron stars at a near-universal filling fraction","feed_subtitle":"Curvature, not the nuclear interaction, sets when spheres morph into spaghetti, lasagna, and bubbles.","key_machinery":"The carrying machinery is the shape-instability comparison inside the compressible liquid-drop model: for a reference shape A and a competing shape B at the same mean density, the paper forms the ratio $\\lambda(\\zeta)$ of the two energy densities after each cell size is optimized, and declares A unstable when $\\lambda(\\zeta)<1$. The surface, curvature, and Coulomb terms enter through the geometric functions $g_s(u)$, $g_c(u)$, and $w(u)$, with the curvature entering through $X_\\sigma=g_c\\sigma_c/(g_s\\sigma_s a)$. This generalizes earlier criteria by solving the cell-size equilibrium condition, Eq. (7), without a perturbative expansion in the curvature. On the ETF side, the mechanism is the extended Thomas-Fermi energy functional with soft-damping nucleon density profiles and the proton-density definition of $u$ given by Eqs. (1)–(2). Using the same functionals in both calculations is what allows the low quasi-universal $u_\\mathrm{sp}$ to be attributed to the curvature correction rather than to the nuclear interaction.","core_discovery":"The central claim is that the shape transitions of nuclear pasta occur at nearly the same filling fractions for the BSk22, BSk24, and BSk25 generalized Skyrme functionals, even though the densities at which those transitions occur are strongly functional-dependent because they are set by the symmetry energy. Spheres lose to spaghetti at $u\\approx0.13$–$0.14$; spaghetti give way to lasagna near 0.26–0.27, lasagna to bucatini near 0.48–0.51, and bucatini to Swiss cheese near 0.66–0.69. The paper then shows that a compressible liquid-drop stability criterion that includes both curvature and neutron-skin corrections reproduces the low onset, $u_\\mathrm{sp}\\approx0.14$–$0.15$, provided the surface and curvature tensions are calibrated to very neutron-rich systems, as in fits to ETF mass tables. Without curvature, the same criterion yields the older universal thresholds near 0.19–0.215 and misses the pasta onset seen in the ETF calculations; with curvature, the thresholds shift down. The paper concludes that the curvature correction is the reason pasta appears near the long-quoted filling fraction of 1/8, and proposes the algebraic criterion as a fast estimator of pasta abundance and as a guide for future ETF searches. For the SLy4 functional, which in older liquid-drop models produced no pasta, the same criterion gives $u_\\mathrm{sp}\\approx0.14$, consistent with the ETF result.","pith_inferences":["[Editorial inference] The quasi-universal ladder suggests the pasta-phase sequence is primarily a geometric phenomenon; one could precompute the four threshold filling fractions once and apply them to any equation of state that provides the proton filling fraction, rather than re-running shape searches for every model.","[Editorial inference] The criterion could be extended to finite temperature for supernova and merger remnants, where surface and curvature tensions vary; the ratio $\\sigma_c/\\sigma_s$ rather than their absolute values would be the controlling input.","[Editorial inference] The explanation implies a testable ordering: models with larger $\\sigma_c/\\sigma_s$ at the relevant proton fractions should show pasta at smaller $u$, and this correlation is sharp enough to distinguish curvature physics from other surface effects.","[Editorial inference] Because the CLDM thresholds depend on whether the tensions are fitted to experimental masses or to ETF mass tables, improved data or many-body calculations for very neutron-rich nuclei would directly tighten the pasta-onset prediction."],"forward_implications":["Pasta onset can be located from geometry alone: once the local filling fraction reaches about 0.13–0.15, a neutron-star crust model should switch from spheres to non-spherical shapes regardless of which of the tested functionals supplies the equation of state.","The density thickness of each pasta layer is controlled by the symmetry energy, so different functionals will disagree on where pasta sits but should agree on the filling-fraction ladder 0.13–0.15, 0.26–0.27, 0.48–0.51, 0.66–0.69.","Equations (12)–(13) provide an inexpensive way to estimate pasta abundance without full ETF minimization, making large surveys of nuclear models tractable.","The same criterion can direct an ETF search: candidate non-spherical shapes need only be evaluated near the predicted filling-fraction thresholds, reducing the cost of the realistic calculation.","For the SLy4 functional, the older liquid-drop finding of no pasta is not reproduced; with curvature included the onset appears near 0.14, consistent with the ETF result."],"supporting_citations":[{"why":"Proved the traditional pasta sequence and its universal filling-fraction thresholds in the curvatureless liquid-drop picture, the baseline the paper extends.","marker":"[3]"},{"why":"Supplies the refined shape-instability criterion without curvature, which the paper generalizes by solving cell-size equilibrium with curvature included.","marker":"[4]"},{"why":"Provides the thermodynamically consistent curvatureless CLDM result with onset at about 0.21, the reference value curvature must lower.","marker":"[6]"},{"why":"Introduced the curvature correction to pasta thresholds; the paper replaces its perturbative treatment with an unapproximated criterion.","marker":"[9]"},{"why":"Shows curvature corrections in a liquid-drop survey can bring the pasta onset into the 0.08–0.14 range, motivating the curvature explanation.","marker":"[12]"},{"why":"Earlier Thomas-Fermi calculations that found the pasta onset near filling fraction 1/8, which the ETF results confirm.","marker":"[16–18]"},{"why":"Prior ETF study of symmetry energy's role in pasta abundance whose framework the present ETF calculation follows.","marker":"[22]"},{"why":"Introduces the soft-damping parametrization of nucleon density profiles used in the ETF calculations.","marker":"[23]"},{"why":"Fit of surface and curvature tensions to experimental atomic masses, used for comparison CLDM predictions.","marker":"[33]"},{"why":"Fit of the surface and curvature tensions to ETF mass tables for BSk22/24/25, the calibration that makes CLDM thresholds agree with ETF.","marker":"[36]"}],"fun_headline_variants":["Curvature, not interaction, sets neutron-star pasta onset","Neutron-star pasta appears at universal fill: curvature explains","Spheres to spaghetti at 13% fill: curvature sets the threshold","Pasta in neutron stars: universal filling fraction tied to curvature","Neutron-star pasta: why shape-shift happens at the same fill"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the plane-surface and curvature tensions taken from mass-table fits remain valid at the extreme neutron excesses and high densities of the deepest crustal layers; the authors themselves state that these tensions are very uncertain there, and if the calibration fails, the curvature-based explanation of $u_\\mathrm{sp}$ would break down even though the direct ETF result would be less affected.","fun_headline_variants_meta":{"raw":{"variants":["Curvature, not interaction, sets neutron-star pasta onset","Neutron-star pasta appears at universal fill: curvature explains","Spheres to spaghetti at 13% fill: curvature sets the threshold","Pasta in neutron stars: universal filling fraction tied to curvature","Neutron-star pasta: why shape-shift happens at the same fill"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000483,"raw_usage":{"total_tokens":2442,"prompt_tokens":1059,"completion_tokens":1383,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":675,"completion_tokens_details":{"reasoning_tokens":1292}},"tokens_in":675,"tokens_out":1383,"duration_ms":12745,"temperature":1.0,"reasoning_tokens":1292,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:32:48.396765+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check would be to compute the sphere-to-pasta onset for the same three functionals with a method that does not rely on the liquid-drop geometry or the fitted density-profile ansatz—for instance a full quantum band-structure or molecular-dynamics calculation at inner-crust densities—and see whether the onset still sits at $u\\approx0.13$–$0.15$; alternatively, computing $\\sigma_c/\\sigma_s$ from first principles for semi-infinite neutron-rich matter and feeding it into Eqs. (12)–(13) would test whether the curvature explanation, rather than the calibration, is responsible.","supporting_citations":[{"cited_title":"Progress of Theoretical Physics 71(2), 320–326 (1984) https://doi.org/10","cited_arxiv_id":null,"evidence_quote":"Proved the traditional pasta sequence and its universal filling-fraction thresholds in the curvatureless liquid-drop picture, the baseline the paper extends."},{"cited_title":"Progress of Theoretical Physics 72(2), 373–375 (1984) https://doi.org/10.1143/PTP.72.373","cited_arxiv_id":null,"evidence_quote":"Supplies the refined shape-instability criterion without curvature, which the paper generalizes by solving cell-size equilibrium with curvature included."},{"cited_title":"Particles 5(3), 225–234 (2022) https://doi.org/10.3390/ particles5030020","cited_arxiv_id":null,"evidence_quote":"Provides the thermodynamically consistent curvatureless CLDM result with onset at about 0.21, the reference value curvature must lower."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced the curvature correction to pasta thresholds; the paper replaces its perturbative treatment with an unapproximated criterion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior ETF study of symmetry energy's role in pasta abundance whose framework the present ETF calculation follows."}],"review_version":1}