{"id":"4fb972e0-8cb2-4b1f-aefb-75a290c38e15","arxiv_id":"2411.17303","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"With a self-consistent Skχ450 liquid drop model, curvature corrections remove bubble pasta phases from the zero-temperature ground state, and lasagna shows a density-dependent crossover between layered and 3D proton superconductivity.","lead":"This paper calculates the shapes of nuclear 'pasta' in neutron star crusts with a liquid drop model, finding that several shapes have nearly equal energy and can coexist. It also estimates how superconducting slabs respond to magnetic fields and pulsar spin lag.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Polymorphism claim rests on fixed-density energy comparisons; Gibbs/Maxwell construction at fixed pressure is needed to establish coexistence.","rationale":"I read the paper in good faith. The liquid-drop machinery appears sound: the reproduction of Lim-Holt's transition densities (Case 1 vs Case 2 in Table I) with sigma_c=0 is a genuine independent check, and computing sigma_s and sigma_c self-consistently with the same SkChi450 interaction is a real improvement over earlier mixed-parametrization treatments. The zero-temperature ground-state statement that curvature removes bubble phases is a well-defined model result, though it would merit uncertainty estimates on sigma_c. The load-bearing weakness is the paper's own central interpretive step: converting small internal-energy differences at fixed n_b into a statement of crystal polymorphism. The energy criterion k_B T does not by itself establish coexistence; genuine phase coexistence requires equality of temperature, pressure, and chemical potentials. The paper is explicit that it does not compute the thermodynamic potential or entropy (Section VI A), so the abstract and conclusions overstate what the calculation shows. This is not an internal inconsistency or a numerics issue; it is a gap between the computed quantity and the claimed physical conclusion. The reader's weakest assumption identifies the same issue, so I agree. The proposed Maxwell-construction test is directly implemented from the paper's own energy curves and would settle whether the polymorphism claim survives; if it does, the paper's conclusions can stand, but the current manuscript should be conditional on that check or on a careful rewording.","tokens_in":33379,"tokens_out":5916,"duration_ms":58266,"concrete_test":"Take the internal energy per baryon curves e(n_b) for the six phases shown in Fig. 6, compute P(n_b) = n_b^2 d(e/n_b)/dn_b, and perform a Maxwell construction (find the convex hull of e/n_b as a function of 1/n_b, or equivalently minimize the Gibbs free energy per baryon at fixed pressure). Determine the equilibrium phase sequence and coexistence density ranges at T=0. Then, at T=10^8 K, add a simple estimate of the leading entropy contributions (e.g., free electrons and phonons) and repeat the free-energy minimization. If the resulting coexistence ranges differ materially from the k_B T proximity bands in Fig. 6, or if bubble phases reappear under the Maxwell construction, the polymorphism claim needs to be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central polymorphism claim (abstract; Section VI A; Fig. 6) is inferred by comparing internal energy per baryon at fixed average baryon density n_b and declaring phases within k_B T of the minimum as coexisting. This is not a thermodynamic coexistence condition. At fixed n_b the phases generally have different pressures, so the equilibrium at fixed pressure requires minimizing the Gibbs free energy per baryon g = e/n_b + P/n_b (or a Maxwell construction). The paper explicitly states (Section VI A) that it focuses solely on the internal energy, leaving thermal effects for future work; it does not compute entropy, the Gibbs free energy, or pressure equality between phases. Because energy differences between pasta phases are of order 0.01-0.1 MeV, comparable to k_B T at T=10^8-10^9 K (0.0086-0.086 MeV), the neglected PV and entropic terms can plausibly change which phases lie within the coexistence window. The 'within k_B T' criterion is an energetic proximity condition, not a phase-equilibrium condition; it indicates possible metastability, not established polymorphism. Thus the main qualitative conclusion of the paper is not supported by the calculation as presented.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript develops a liquid drop model for neutron star inner crust pasta phases using the Skχ450 interaction, with the planar surface tension and curvature correction computed self-consistently from the same interaction via extended Thomas-Fermi calculations. The authors solve the full set of variational equations (Eqs. 4–7) for spherical, cylindrical, and planar nuclei and bubbles. They report that including the curvature correction removes bubble phases from the zero-temperature ground state and makes lasagna the ground state over a significant density range (Table I, Case 3). They then identify phases within k_B T of the minimum internal energy per baryon and interpret this set as coexisting polymorphs (Fig. 6). The paper also introduces a model of proton superconductivity in lasagna, predicts a crossover between discrete layered and anisotropic three-dimensional superconductivity, and estimates the magnetic stress caused by a rotational lag between superfluid neutrons and the lattice.","tokens_in":33632,"tokens_out":6374,"duration_ms":57370,"significance":"The self-consistent treatment of surface and curvature energies is a methodological advance, and the numerical machinery is convincingly verified by reproducing the Lim–Holt results when curvature is omitted. The computed structural parameters and the estimates for the superconducting coherence length, penetration depth, and magnetic stress are useful quantitative inputs for crust physics. However, the central polymorphism claim rests on an energetic proximity criterion rather than a thermodynamic phase-equilibrium calculation, which limits the support the paper provides for coexistence of multiple phases.","major_comments":[{"comment":"The polymorphism conclusion is drawn by comparing the internal energy per baryon at fixed average baryon density n_b and declaring phases within k_B T of the minimum as coexisting. This is not a thermodynamic coexistence condition: at fixed n_b, the phases generally have different pressures, so equilibrium at fixed pressure requires minimizing the Gibbs free energy per baryon g = e/n_b + P/n_b or performing a Maxwell construction. The paper itself states in Section VI A that it focuses solely on the internal energy and leaves thermal effects for future work, and the neglected PV and entropic terms are of the same order as the energy differences (0.01–0.1 MeV versus k_B T ≈ 0.0086–0.086 MeV at T = 10^8–10^9 K). Consequently, the 'within k_B T' criterion can only indicate energetic proximity, not coexistence, and the main qualitative claim of polymorphism is not supported by the calculation as presented.","section":"Section VI A, Fig. 6"},{"comment":"The disappearance of bubble phases from the ground state is a strong claim that depends sensitively on the curvature correction, whose energy scale (σ_c ≲ 0.6 MeV/fm, Eq. (C19) and Table III) is much smaller than the planar surface energy. The authors do not quantify the sensitivity of the phase boundaries to uncertainties in σ_c(x) or to the neglected neutron skin and density smearing corrections (the latter is dismissed by a private communication, ref. [55]). A robustness check, for example varying σ_c by its numerical uncertainty or including the next-order term, would be needed to establish that the disappearance of bubbles is not an artifact of the particular ETF evaluation.","section":"Table I, Case 3; Section IV C"}],"minor_comments":[{"comment":"The word 'discreet' should be 'discrete' in all occurrences (e.g., Fig. 8 caption, Section VI B, Conclusions).","section":"Throughout"},{"comment":"The text contains 'Newtown's second law'; this should be 'Newton's second law'.","section":"Section VII A"},{"comment":"The phrase 'thermodynamically allowed pasta configurations' is misleading because the selection criterion is only internal energy proximity; suggest 'energetically close configurations.'","section":"Fig. 6 caption"},{"comment":"The statement that T = 10^8–10^9 K is 'very small compared to the typical nuclear energy scale of 1 MeV' is misleading, since the relevant comparison is with the inter-phase energy differences, which are comparable to k_B T. Please rephrase.","section":"Section VI A"},{"comment":"Ref. [43] is an unpublished self-citation and Ref. [55] is a private communication; both should be updated or substantiated with publicly available work.","section":"References"},{"comment":"The analytic fit to σ_s(x) is presented, but the main calculation uses spline interpolation of the tabulated values; please clarify whether the fit is used anywhere in the reported results.","section":"Eq. (11)"}],"recommendation":"major_revision","confidential_remarks":"The paper's central polymorphism claim is likely to be controversial because it does not use a thermodynamic potential. Even if the energy differences are correct, the interpretation as coexistence is not established. The authors should either compute the Gibbs free energy or weaken the claim. The heavy reliance on unpublished references (private communication and unpublished self-citation) should be addressed. The manuscript is otherwise well-structured and the numerical work is transparent."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe genuine content here is the self-consistent liquid-drop calculation with surface and curvature terms from the same Skχ450 interaction, and the concrete superconductivity crossover estimate in lasagna. The numerical verification against Lim-Holt (curvature off) is clean and gives real confidence in the machinery. The main new outcome—curvature killing the bubble phases in the ground state and shifting transition densities—is a solid, useful result for the crust subfield.\n\nThe soft spot is the polymorphism claim. Comparing internal energies per baryon at fixed n_b and calling phases within k_B T of the minimum 'coexisting' is not a thermodynamic coexistence condition. The phases generally have different pressures at fixed n_b, so the equilibrium at fixed pressure requires Gibbs free energy or a Maxwell construction. The authors explicitly say they only consider internal energy and leave thermal effects to future work, but the abstract and conclusions still assert the crust is likely polymorphic. Given the energy differences are ~0.01–0.1 MeV, comparable to k_B T at 10^8–10^9 K, the PV and entropic terms could plausibly change the 'coexisting' set. This makes the polymorphism an intriguing conjecture, not a demonstrated result. It should be either upgraded with a Gibbs construction or clearly labeled as an energetic proximity criterion.\n\nThe stress estimate is honest—they state D=1 cm is illustrative and probably an overestimate—and the superconductivity model is admittedly crude but still gives a reasonable first orientation. The main numerical results lack propagated uncertainties, though that's typical for LDM studies.\n\nWho benefits: anyone working on inner crust structure, pasta phases, or superconducting pasta. It deserves a serious referee; the self-consistent curvature treatment and the superconductivity crossover are worth referee time. I'd expect major revision on the thermodynamics before acceptance, but this is not a desk reject.","headline":"Self-consistent LDM with curvature kills bubble phases and gives a plausible superconductivity crossover, but the polymorphism claim needs a proper Gibbs construction.","tokens_in":34089,"tokens_out":3423,"would_cite":true,"duration_ms":33360,"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":"Self-consistent curvature removes bubble pasta from the neutron-star crust's ground state, making lasagna the stable phase over a wide density range.","keywords":["neutron star inner crust","nuclear pasta","liquid drop model","curvature energy","crystal polymorphism","proton superconductivity","lasagna phase","Skyrme interaction Skχ450"],"falsifier":"Compute the Gibbs free energy per baryon, including phonon and pairing entropies, at fixed total pressure and temperature over $0.16 \\lesssim n_b/n_0 \\lesssim 0.55$ for the same Skχ450 liquid-drop model. If the free-energy differences between the 3N, 2N, 1N, and uniform phases are everywhere larger than $k_B T$ at $T = 10^8$ K, the polymorphic-coexistence claim fails; alternatively, a full Hartree-Fock-Bogoliubov calculation with the same interaction that finds stable 2B or 3B phases would refute the curvature-driven disappearance of bubbles.","tokens_in":33162,"feed_emoji":"⭐","tokens_out":12293,"duration_ms":103533,"temperature":0.7,"pith_summary":"Using the liquid drop model with one chiral-effective-field-theory interaction (Skχ450) for the bulk, the planar surface, and the curvature terms, the paper calculates the internal energy of every pasta phase in the inner crust of a neutron star. Its central claim is that the curvature correction, computed self-consistently from the same interaction as the bulk energy, changes the ground state: the bubble phases 2B and 3B disappear and lasagna (slab-shaped nuclei) is stable over a wide density range, with the sequence 3N → 2N → 1N → uniform matter. The paper also shows that at temperatures $10^8$–$10^9$ K the energy per baryon of several phases lies within one thermal energy of the minimum, so it argues the crust is polymorphic rather than a single crystal. From the same solution it finds that proton Cooper pairing in lasagna crosses from discrete layered to anisotropic three-dimensional superconductivity, and it estimates a magnetic stress from the neutron–proton rotational lag that can reach the crust's breaking scale. These results matter because they change the mechanical, electrical, and magnetic properties of the inner crust, the layer that governs pulsar glitches and part of magnetar behaviour.","feed_headline":"Self-consistent curvature erases bubble phases in neutron-star crust","feed_subtitle":"Same chiral interaction for bulk and surface leaves lasagna stable, with several pasta shapes coexisting.","key_machinery":"The central object is the compressible liquid drop model of the Wigner-Seitz cell, whose total energy density is split into bulk, planar surface, Coulomb plus lattice, dripped-neutron, electron, and curvature terms. The load-bearing piece is the curvature term $w_{\\rm curv} = u d(d-1)\\sigma_c(x)/r_N^2$, with $\\sigma_c(x)$ taken from the same extended Thomas-Fermi calculation, using Skχ450, that supplies the planar surface tension $\\sigma_s(x)$. It enters through the nuclear virial theorem $w_{\\rm p.surf} + 2w_{\\rm curv} = 2w_{\\rm C+L}$, which fixes the nucleus radius $r_N$ and hence the densities at which one pasta symmetry gives way to another; the other variational equations enforce $\\beta$-equilibrium, equal neutron chemical potentials, and pressure balance across the interface. For the superconductivity part, the deciding ratios are the proton coherence length $\\xi$ against the slab width $2r_N$ and the inter-slab distance $d_L$, with the London depth $\\lambda$ setting the type-II condition $\\lambda/\\xi > \\sqrt{2}$.","core_discovery":"The paper's central discovery is that the curvature contribution $\\sigma_c(x)$ to the nuclear surface tension, extracted from the same extended Thomas-Fermi calculation with the Skχ450 interaction that gives the planar tension $\\sigma_s(x)$, is not a small correction but a phase-diagram-changing one. Through the virial balance $w_{\\rm p.surf} + 2 w_{\\rm curv} = 2 w_{\\rm C+L}$ it enlarges the equilibrium nucleus size and shifts the pasta transition densities upward, so that the bubble phases 2B and 3B, which appear in earlier liquid-drop calculations without curvature, drop out of the zero-temperature ground state; the ground-state sequence becomes 3N → 2N → 1N → uniform, with lasagna occupying a substantial density interval below the uniform-matter instability at $n_b/n_0 \\approx 0.5508$. At $T = 10^8$–$10^9$ K, the internal energy per baryon of non-minimal phases (3N, 2N, 1N, sometimes 2B or uniform) lies within $k_B T$ of the minimum, which the authors read as evidence for a polymorphic crust. Using the same structure solution, proton pairing in lasagna is type-II ($\\lambda/\\xi > \\sqrt{2}$), with a crossover from discrete layered superconductivity to anisotropic three-dimensional superconductivity when the coherence length exceeds the inter-slab spacing near $n_b \\approx 0.48 n_0$; the associated magnetic stress, evaluated for $B_0 = 5\\times 10^{14}$ G and a neutron–proton velocity lag of 1 cm s$^{-1}$, can exceed the crust's yielding stress.","pith_inferences":["Extending beyond the paper: the disappearance of bubble phases is likely specific to the Skχ450 interaction; repeating the identical liquid-drop construction with other chiral-EFT parameter sets would show whether 2B/3B survive and how much the lasagna window shrinks or grows.","The polymorphism band defined by internal-energy differences at fixed average density is not a Maxwell construction; a proper Gibbs treatment at fixed pressure and temperature with phonon entropy could widen, narrow, or split the coexistence windows, and it would decide which phase fractions actually populate the crust.","The superconductivity crossover near $n_b \\approx 0.48 n_0$ implies a thin density shell where the magnetic response changes character; in a polymorphic crust this makes vortex pinning spatially patchy, a feature that could be tested against magnetar quasi-periodic oscillation spectra and glitch recovery statistics.","The magnetic-stress estimate assumes defect-free, periodic lasagna; real crusts with grain boundaries and dislocations might yield at lower stress, or the disorder could decouple slabs and suppress the coherent force, so the shattering threshold is an upper bound rather than a precise prediction."],"forward_implications":["The zero-temperature ground state of the inner crust is 3N → 2N → 1N → uniform matter, with no bubble phases; the lasagna (1N) layer is much wider than in curvature-free models, so crust transport, elasticity, and cooling calculations should be re-done with slab geometry.","At $T \\sim 10^8$–$10^9$ K several pasta symmetries have energy within $k_B T$ of the minimum, so the crust is expected to solidify as a polymorphic mixture of 3N, 2N, 1N, and sometimes 2B or uniform matter rather than a single crystal; the local microstructure then depends on thermal history.","Proton superconductivity in lasagna is type-II, and it switches from discrete layered to continuous anisotropic three-dimensional behaviour near $n_b \\approx 0.48 n_0$; magnetic vortex, pinning, and critical-field calculations must use different models on either side of that density.","A rotational lag between superfluid neutrons and the lattice in superconducting lasagna produces a magnetic force that, for $B_0 = 5\\times10^{14}$ G and $\\delta v \\sim 1$ cm s$^{-1}$, is of order $10^{33}$ dyn per $1\\,\\mathrm{cm}^2 \\times 1\\,\\mathrm{cm}$ column, exceeding the crust's yielding stress; spin-down and glitches could therefore crack the inner crust."],"supporting_citations":[{"why":"Provides the Skχ450 interaction, the bulk energy functional, the tabulated surface tension and curvature data, and the Case-1 critical densities reproduced for verification.","marker":"[24]"},{"why":"Supplies the extended Thomas-Fermi method used to compute the planar surface tension and curvature correction self-consistently.","marker":"[34]"},{"why":"Extends the same surface-curvature formalism to neutron-star crust conditions and fixes the sign and form of the curvature term.","marker":"[35]"},{"why":"Gives the proton pairing gap as a function of Fermi momentum from chiral EFT, the input used for the superconductivity analysis in lasagna.","marker":"[49]"},{"why":"Establishes the superconducting-superfluid magnetohydrodynamic stress model for lasagna that the paper quantifies with its microscopic parameters.","marker":"[41]"},{"why":"Supplies the continuous anisotropic superconductivity model for pasta; the paper shows it applies only above the crossover density.","marker":"[42]"},{"why":"Classic liquid-drop treatment showing small energy differences between pasta phases and providing the surface-tension parametrization used in earlier cases.","marker":"[8]"},{"why":"Foundational liquid-drop description of neutron-star matter that defines the unit-cell energy decomposition this paper refines.","marker":"[4]"}],"fun_headline_variants":["Curvature kills bubble pasta in neutron-star crust","Neutron-star crust: curvature erases bubbles, keeps lasagna","Curvature reshapes neutron-star crust, drops bubble phases","Superconductivity crossover in neutron-star lasagna","Neutron-star pasta: bubbles vanish, lasagna stable"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The polymorphism conclusion rests on comparing internal energies per baryon at fixed average density and calling phases within one thermal energy of the minimum 'coexisting'; the paper does not compute the free energy at fixed pressure and temperature or include entropy, so the coexistence claim is not thermodynamically established.","fun_headline_variants_meta":{"raw":{"variants":["Curvature kills bubble pasta in neutron-star crust","Neutron-star crust: curvature erases bubbles, keeps lasagna","Curvature reshapes neutron-star crust, drops bubble phases","Superconductivity crossover in neutron-star lasagna","Neutron-star pasta: bubbles vanish, lasagna stable"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000531,"raw_usage":{"total_tokens":2683,"prompt_tokens":1197,"completion_tokens":1486,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":813,"completion_tokens_details":{"reasoning_tokens":1403}},"tokens_in":813,"tokens_out":1486,"duration_ms":11171,"temperature":1.0,"reasoning_tokens":1403,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:15:13.044421+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Gibbs free energy per baryon, including phonon and pairing entropies, at fixed total pressure and temperature over $0.16 \\lesssim n_b/n_0 \\lesssim 0.55$ for the same Skχ450 liquid-drop model. If the free-energy differences between the 3N, 2N, 1N, and uniform phases are everywhere larger than $k_B T$ at $T = 10^8$ K, the polymorphic-coexistence claim fails; alternatively, a full Hartree-Fock-Bogoliubov calculation with the same interaction that finds stable 2B or 3B phases would refute the curvature-driven disappearance of bubbles.","supporting_citations":[{"cited_title":"Unit of Excellence Mar´ıa de Maeztu 2020-2023","cited_arxiv_id":null,"evidence_quote":"Establishes the superconducting-superfluid magnetohydrodynamic stress model for lasagna that the paper quantifies with its microscopic parameters."},{"cited_title":"Notice that, again, there are no contributions from neutron skin even if the corresponding energy density were included intow tot","cited_arxiv_id":null,"evidence_quote":"Foundational liquid-drop description of neutron-star matter that defines the unit-cell energy decomposition this paper refines."}],"review_version":1}