Pith. sign in

REVIEW 2 major objections 6 minor 83 references

Inner crust of neutron stars: Polymorphism and superconductivity in the liquid drop model

T0 review · 2 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read Self-consistent curvature removes bubble pasta from the neutron-star crust's ground state, making lasagna the stable phase over a wide density range.

desk verdict Self-consistent LDM with curvature kills bubble phases and gives a plausible superconductivity crossover, but the polymorphism claim needs a proper Gibbs construction. read the letter →

arxiv 2411.17303 v2 pith:2H34ENNN submitted 2024-11-26 nucl-th astro-ph.HEastro-ph.SRhep-thphysics.plasm-ph

classification nucl-thastro-ph.HEastro-ph.SRhep-thphysics.plasm-ph
keywords neutronstarinnercrustnuclearpastaliquiddropmodelcurvatureenergycrystalpolymorphismprotonsuperconductivitylasagnaphaseSkyrmeinteractionSkχ450
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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}$.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 6 minor

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.

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 (2)
  1. [Section VI A, Fig. 6] 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.
  2. [Table I, Case 3; Section IV C] 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.
minor comments (6)
  1. [Throughout] The word 'discreet' should be 'discrete' in all occurrences (e.g., Fig. 8 caption, Section VI B, Conclusions).
  2. [Section VII A] The text contains 'Newtown's second law'; this should be 'Newton's second law'.
  3. [Fig. 6 caption] The phrase 'thermodynamically allowed pasta configurations' is misleading because the selection criterion is only internal energy proximity; suggest 'energetically close configurations.'
  4. [Section VI A] 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.
  5. [References] Ref. [43] is an unpublished self-citation and Ref. [55] is a private communication; both should be updated or substantiated with publicly available work.
  6. [Eq. (11)] 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.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the central structure calculation is self-consistent rather than fitted, and the cited prior-work inputs are not used to force the paper's main conclusions.

full rationale

The paper's main structural results are obtained by solving the liquid-drop variational equations with bulk, planar surface, and curvature terms all derived from the same Skχ450 interaction. This is an internal-consistency feature, not a circular reduction: the surface tension and curvature are computed from extended Thomas-Fermi solutions of the same energy functional (Appendix D), and no pasta-phase transition density or coexistence interval is used as input to fix those quantities. The authors explicitly avoid fitting the surface functions in the main calculation, stating that they 'represent these functions numerically with the help of spline interpolation – in this way we avoid choosing any fitting parameters such as in Eq. (9)' (Section IV C). The verification against Lim and Holt [24] (Case 2 in Table I) is a benchmark check, not a prediction built from fitted pasta data. The pairing gap is imported from a separate uniform-matter calculation [49], and the entrainment and stress expressions are taken from prior published work [41, 61]; these are external inputs and do not constitute a self-citation chain that forces the polymorphism or superconductivity-crossover conclusions. The polymorphism claim is based on comparing internal energy per baryon at fixed average density, and the paper explicitly acknowledges that a full thermodynamic treatment would require minimizing a thermodynamic potential including entropy, stating: 'In this work we focus solely on the internal energy, while leaving the thermal effects for the future work' (Section VI A). That is an honest limitation of the coexistence interpretation, not a circular derivation. The self-citations to [41], [43], and [63] are present but are not load-bearing for the paper's central claims: [43] is an unpublished motivational estimate, and [41]/[63] supply a stress-tensor formula that is then evaluated with new structural inputs. No step was found in which a 'prediction' reduces by construction to a fitted parameter or to a self-citation.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The central calculation introduces no new fitted parameters and no invented physical entities. The main assumptions are the liquid drop idealizations (sharp interfaces, uniform densities, no neutron skin), the validity of the Skχ450 interaction for both bulk and surface terms, and the extension of uniform-matter pairing gaps to the lasagna geometry. The stress estimate uses illustrative scenario parameters D and delta_v rather than fitted values.

free parameters (3)
  • D (lasagna column height) = 1 cm
    Chosen for the stress estimate in Eq. (46). The authors note D is likely much smaller in real neutron stars, so the quoted force is illustrative, not a fitted parameter.
  • delta_v (neutron-proton velocity lag) = 1 cm/s
    Assumed from pulsar spin lag observations (max lag ~1e-6 s^-1); used to scale the magnetic stress formula.
  • Thermal energy threshold temperatures = 10^8, 5x10^8, 10^9 K
    Three illustrative temperatures used to define the polymorphism window in Fig. 6; not fitted to data.
assumptions (6)
  • domain assumption Liquid drop model with sharp interfaces and uniform densities inside and outside the cluster
    Introduced in Section III A; neglects neutron skin and density smearing, relying on B.K. Sharma private communication [55] that these do not significantly affect transitions.
  • domain assumption Electrons form a uniform, degenerate ultrarelativistic background within the Wigner-Seitz cell
    Section III A and Eq. (B19); standard approximation for cold neutron star crust matter.
  • domain assumption Skχ450 Skyrme parametrization from Lim-Holt [24] describes bulk and surface nuclear energies
    The interaction is constrained by chiral EFT and nuclear data in prior work; this paper does not derive it.
  • domain assumption Extended Thomas-Fermi expansion truncated at second order in hbar provides accurate surface tension and curvature correction
    Appendix D; only the first two terms of the hbar expansion are kept, and the curvature correction is added as a linear term.
  • ad hoc to paper Uniform-matter proton pairing gaps from [49] apply locally inside lasagna slabs with a sharp superconducting density profile
    Section V; authors explicitly call this a crude approximation because the gap may be anisotropic and matter inside clusters is inhomogeneous.
  • ad hoc to paper Thermodynamic coexistence can be assessed by energy-per-baryon differences at fixed average baryon density, without computing Gibbs free energy or entropy
    Section VI A states the authors focus solely on internal energy and leave thermal effects for future work, yet the conclusions invoke polymorphic coexistence.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Inner crust of neutron stars: Polymorphism and superconductivity in the liquid drop model." pith.science (2026). https://pith.science/paper/2H34ENNN

@misc{pith2026241117303,
  author       = {Pith},
  title        = {Pith review of: Inner crust of neutron stars: Polymorphism and superconductivity in the liquid drop model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2H34ENNN}},
  note         = {Machine review of arXiv:2411.17303}
}
abstract

Within the liquid drop model built up with the nuclear interaction parametrization Sk$\chi$450, which is based on the chiral effective field theory, we calculate numerically the internal energy density for each of nuclear pasta phases and for the uniform nuclear matter. We provide quantitative arguments in favor of coexistence of various nuclear matter phases at a significant range of total pressure within the inner crust of neutron stars, a concept known as crystal polymorphism. Specifically, we find that differences of the internal energy per baryon for various phases are typically less than the thermal energy per a freedom degree at temperature about $10^8$--$10^9$ K, which sets the energetic scale for thermal fluctuations of state of Fermi liquid from the ground state. The nuclear energy contributions are described using the same parametrization Sk$\chi$450 for the bulk, plain surface and curvature terms. We find that the introduction of the curvature correction changes the ground state in a relevant way. This may be understood as a consequence of the corresponding change in size of the nucleus, which significantly modifies the phase transition densities. Using the calculated structural parameters from liquid drop model, we explore the physical consequences of the expected Cooper pairing of protons in lasagna phase. In this case, we find a crossover between the discreet layered and the three-dimensional anisotropic regimes of superconductivity. Additionally, we study the magnetic stress in lasagna accounting for a rotational lag between superfluid neutrons and the crystal lattice, which is believed to develop naturally in pulsars and magnetars. Our results offer a preliminary insight into rich magnetic properties of the inner crust of neutron stars.

Figures

Figures reproduced from arXiv: 2411.17303 by the authors.

Figure 1
Figure 1. FIG. 1. Solution to the basic equations of equilibrium, Eqs. (4)-(7): [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Solution to the basic equations of equilibrium, Eqs. (4)- [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 5
Figure 5. FIG. 5. The workfunction to transfer a proton from the nucleus to [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Solution to the basic equations of equilibrium, Eqs. (4)-(7): [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 7
Figure 7. Figure 7: of [49], where the kF p-dependence has been converted into the nb-dependence with the help of the functions shown in Figs. 1–4 and Eq. (20). where the quantity Yp (the proton fraction averaged over the Wigner-Seitz cell) utilized in [49] is given in variables of the pr…
Figure 8
Figure 8. Figure 8: FIG. 8. Schematic representation of a unit length of ideally ordered [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]
Figure 10
Figure 10. Figure 10: FIG. 10. (Color online). Panel (a): A schematic representation (not [PITH_FULL_IMAGE:figures/full_fig_p011_10.png]
Figure 11
Figure 11. Figure 11: FIG. 11. Numerical evaluation of Eq. (46) in lasagna with the vari [PITH_FULL_IMAGE:figures/full_fig_p013_11.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

83 extracted references · 33 canonical work pages

  1. [55]

    D. N. Kobyakov, Phys. Rev. C98, 045803 (2018). Application of Superconducting-Superfluid Magneto- hydrodynamics to Nuclear “Pasta” in Neutron Stars. https://doi.org/10.1103/physrevc.98.045803

  2. [1]

    Energy density The protons are uniformly distributed within the nucleus with volumeV N with the number density np = Np VN .(B3) Here, VN =    (4/3)πr 3 N,(3N) πr 2 N ×L,(2N) 2rN ×L 2,(1N) (B4) wherer N is the radius of the nucleus,L→∞is the length of the rod-like nucleus andL 2 →∞is the area of the slab- like nucleus, see Fig. 11 (b). The limitL→∞sho...

  3. [2]

    In order to reveal the physical signifi- cance of the variational equations we shall discuss in detail each of the equations

    Variational equations The variational equations have the form ∂w tot ∂y =0,(B22) whereyis either of the variables from the set {n,x,n no,r N,u}. In order to reveal the physical signifi- cance of the variational equations we shall discuss in detail each of the equations. 17

  4. [3]

    suggests that, the inclusion ofσ c(x)in this paper induces a qualitative difference for the ground state of the inner crust matter, as compared with the ground state predicted earlier in [24]. B. Energies of pairing, neutron skin and density smearing In the main numerical calculations of this paper, we have ignored energies of pairing, neutron skin and de...

  5. [4]

    Notice that, again, there are no contributions from neutron skin even if the corresponding energy density were included intow tot

    Chemical equilibrium The extremum ∂(w tot) ∂x n,nno,rN ,u =0 (B31) implies that the isospin of nucleus is in equilibrium with re- spect to weak interactions: µe + (mp −m n)c2 =− ∂ ε ∂x n (B32) − 1 un ∂(w p.surf +w C+L +w curv) ∂x n,nno,rN ,u . Notice that, again, there are no contributions from neutron skin even if the corresponding energy density were in...

  6. [5]

    Application of the derivative to Eq

    Optimum size of nucleus The extremum ∂(w tot) ∂r N n,x,nno,u =0 (B23) with respect to the nucleus sizer N at fixed{n,x,n no,u}im- plies that the energy per baryon is minimized by nucleus size while the number densities are fixed. Application of the derivative to Eq. (B12) yields the result known as the nuclear virial theorem with a correction due to the p...

  7. [6]

    As a result, the pressure in the nucleus is equal to the pressure in the dripped neutrons,P i =P o

    Continuity of pressure across the nucleus interface The extremum ∂(w tot) ∂u Nn,Nno,x,Vc =0 (B47) implies that the energy density is minimized by choosing the optimal size of the nucleus while keeping fixed the cell size, the number of neutrons, both inside and outside the nucleus, and the proton fraction. As a result, the pressure in the nucleus is equal...

  8. [7]

    As a result, the energy to add a neutron to the nucleus is equal to the energy to add a neutron to the dripped neutrons, µni =µ no

    Continuity of neutron chemical potential across the nucleus interface The extremum ∂(w tot) ∂n ni Nn,npi,rN ,u =0 (B35) implies that the total energy density is minimized by choosing the optimal number of neutrons inside and outside of the nu- cleus while keeping their total number fixed as well as fixing the proton number, the size of the nucleus and tha...

Show all 83 references
  1. [9]

    The bubble volume fraction is ubub = VB Vc = rB rc d ,(C5) wherer c is the radius of the unit cell,d=3 for 3B andd=2 for 2B

    Energy density The protons are uniformly distributed outside the bubble with volumeV B with the number density np = Np Vc −VB .(C3) Here, VB =    (4/3)πr 3 B,(3B) πr 2 B ×L,(2B) 2rB ×L 2,(1B) (C4) wherer B is the radius of the bubble,L→∞is the length of the rod-like bubbl...

  2. [10]

    Variational equations The variational equations have the form ∂w bub tot ∂y =0,(C20) whereyis either of the variables from the set {n,x,n no,r B,u bub}

  3. [11]

    Optimum size of bubble The extremum ∂(w bub tot ) ∂r B n,x,nno,ubub =0 (C21) leads to the nuclear virial theorem: wbub p.surf +2w bub curv =2w bub C+L.(C22) From Eq. (C22) we find the equation that determinesr N: r4 B −4q bubrB +3r bub =0,(C23) where qbub = 1 4 σs(x)d 4π(enx)2...

  4. [12]

    (C39) and (C41), respectively

    Chemical equilibrium The extremum ∂(w bub tot ) ∂x n,nno,rB,ubub =0 (C29) provides the beta-equilibrium condition: µbub e + (mp −m n)c2 =(C30) − ∂ ε ∂x n − 1 (1−u bub)n ∂(w bub s +w bub curv +w bub C+L) ∂x n,nno,rB,ubub , where µbub e = ¯hc h 3π2(1−u bub)xn i 1 3 .(C31) Again,...

  5. [13]

    Continuity of neutron chemical potential across the bubble interface The extremum ∂(w bub tot ) ∂n ni Nn,npi,rB,ubub =0 (C33) implies that the energy to add a neutron to the nuclear matter outside the bubble is equal to the energy to add a neutron to the dripped neutrons insid...

  6. [14]

    As a result, the pressure inside the bubble is equal to that outside the bubble,P bub i =P bub o

    Continuity of pressure across the bubble interface The extremum ∂(w bub tot ) ∂u bub Nn,Nno,x,Vc =0 (C42) implies that the energy density is minimized by choosing the optimal size of the bubble while keeping fixed the cell size, the number of neutrons, both inside and outside ...

  7. [15]

    G. E. Brown, H. A. Bethe and Gordon Baym, Nucl. Phys. A375, 481 (1982). Supernova theory. https://doi.org/10.1016/0375-9474(82)90025-2

  8. [16]

    Phys.: Conf

    Kazuhiro Oyamatsu, Kei Iida and Hajime Sotani, J. Phys.: Conf. Ser.1643, 012059 (2020). Systematic study of pasta nuclei in neutron stars with families of the empiri- cal nuclear equations of state. https://doi.org/10.1088/1742- 6596/1643/1/012059

  9. [17]

    Turolla, S

    R. Turolla, S. Zane and A. L. Watts, Rep. Progr. Phys.78, 116901 (2015). Magnetars: the physics behind observations. A review. https://doi.org/10.1088/0034-4885/78/11/116901

  10. [19]

    D. G. Ravenhall, C. J. Pethick and J. R. Wilson, Phys. Rev. Lett. 50, 2066 (1983). Structure of Matter below Nuclear Saturation Density. https://doi.org/10.1103/PhysRevLett.50.2066

  11. [20]

    Hashimoto, H

    M. Hashimoto, H. Seki and M. Yamada, Progr. Theor. Phys. 71, 320 (1984). Shape of Nuclei in the Crust of Neutron Star. https://doi.org/10.1143/PTP.71.320

  12. [21]

    M. E. Caplan and C. J. Horowitz, Rev. Mod. Phys.89, 041002 (2017). Colloquium: Astromaterial science and nuclear pasta. https://doi.org/10.1103/RevModPhys.89.041002

  13. [22]

    Ken’ichiro Nakazato, Kei Iida and Kazuhiro Oyamatsu, Phys. Rev. C83, 065811 (2011). Curvature effect on nuclear “pasta”: Is it helpful for gyroid appearance? https://doi.org/10.1103/PhysRevC.83.065811

  14. [24]

    Magierski and P.-H

    P. Magierski and P.-H. Heenen, Phys. Rev. C65, 045804 (2002). Structure of the inner crust of neu- tron stars: Crystal lattice or disordered phase? https://doi.org/10.1103/PhysRevC.65.045804

  15. [25]

    J. W. Negele and D. Vautherin, Nucl. Phys. A207, 298 (1973). Neutron Star Matter at Sub-Nuclear Densities. https://doi.org/10.1016/0375-9474(73)90349-7

  16. [26]

    Grill, J

    F. Grill, J. Margueron and N. Sandulescu, Phys. Rev. C 84, 065801 (2011). Cluster structure of the inner crust of neutron stars in the Hartree-Fock-Bogoliubov approach. https://doi.org/10.1103/PhysRevC.84.065801

  17. [27]

    Comparison between the Thomas–Fermi and Hartree–Fock–Bogoliubov Methods in the Inner Crust of a Neutron Star: The Role of Pairing Correlations

    Matthew Shelley and Alessandro Pastore, Universe6, 206 (2020). Comparison between the Thomas–Fermi and Hartree–Fock–Bogoliubov Methods in the Inner Crust of a Neutron Star: The Role of Pairing Correlations. https://doi.org/10.3390/universe6110206

  18. [28]

    Mondal, X

    C. Mondal, X. Vi ˜nas, M. Centelles and J.N. De, Phys. Rev. C102, 015802 (2020). Structure and Composition of the Inner Crust of Neutron Stars from Gogny Interactions. https://doi.org/10.1103/PhysRevC.102.015802

  19. [29]

    W. G. Newton, S. Cantu, S. Wang, A. Stinson, M. A. Kaltenborn and J. R. Stone, Phys. Rev. C105, 025806 (2022). Glassy Quantum Nuclear Pasta in Neutron Star Crusts. https://doi.org/10.1103/PhysRevC.105.025806

  20. [31]

    Kazuyuki Sekizawa, Sorataka Kobayashi and Masayuki Mat- suo, Phys. Rev. C105, 045807 (2022). Time-dependent extension of the self-consistent band theory for neutron star matter: Anti-entrainment effects in the slab phase. https://doi.org/10.1103/PhysRevC.109.065804

  21. [32]

    J. M. Pearson and N. Chamel, Phys. Rev. C105, 015803 (2022). Unified equations of state for cold nonac- creting neutron stars with Brussels-Montreal functionals. III. Inclusion of microscopic corrections to pasta phases. https://doi.org/10.1103/PhysRevC.105.015803

  22. [33]

    B. K. Sharma, M. Centelles, X. Vi ˜nas, M. Baaldo and F. Burgio, Astron. Astrophys.584, A103 (2015). Unified Equa- tion of State for Neutron Stars on a Microscopic Basis. https://doi.org/10.1051/0004-6361/201526642

  23. [34]

    W. D. Myers and W. J. Swiatecki, Ann. Phys.55, 395 (1969). Average nuclear properties. https://doi.org/10.1016/0003- 4916(69)90202-4

  24. [35]

    W. D. Myers and W. J. Swiatecki, Ann. Phys.84, 186 (1974). The nuclear droplet model for arbitrary shapes. https://doi.org/10.1016/0003-4916(74)90299-1

  25. [36]

    P. B. Jones, MNRAS243, 257 (1990). Rotation of the neutron-drip superfluid in pulsars: the resistive force. https://ui.adsabs.harvard.edu/abs/1990MNRAS.243..257J

  26. [37]

    W. G. Newton, M. Gearheart and Bao-An Li, Astrophys. J. Supp. Ser.204, 9 (2013). A Survey of the Parameter Space of the Compressible Liquid Drop Model as Applied to the Neutron Star Inner Crust. https://doi.org/10.1088/0067-0049/204/1/9

  27. [38]

    Lim and J

    Y . Lim and J. W. Holt, Phys. Rev. C95, 065805 (2017). Structure of Neutron Star Crusts from New Skyrme Effec- tive Interactions Constrained by Chiral Effective Field Theory. https://doi.org/10.1103/PhysRevC.95.065805

  28. [39]

    Carreau, F

    T. Carreau, F. Gulminelli, N. Chamel, A.F. Fantina and J.M. Pearson, Astron. Astrophys.635, A84 (2020). Crystallization of the Inner Crust of a Neutron Star and the Influence of Shell Effects. https://doi.org/10.1051/0004-6361/201937236

  29. [40]

    Carreau, A.F

    T. Carreau, A.F. Fantina and F. Gulminelli, Astron. As- tropys.640, A77 (2020). Inner Crust of a Neutron Star at the Point of Crystallization in a Multicomponent Approach. https://doi.org/10.1051/0004-6361/202038347

  30. [41]

    Unit of Excellence Mar´ıa de Maeztu 2020-2023

    (see also a corrected expression in [63]). As compared with the work done in [41], here we will specify the quantita- tive evaluation of the force based on the microscopic param- eters calculated in Sec. IV and generalize the expression for the force for the case whenA 0 and p...

  31. [42]

    Dinh Thi, A

    H. Dinh Thi, A. F. Fantina and F. Gulmonelli, Eur. Phys. J. A57, 296 (2021). The Effect of the Energy Functional 24 on the Pasta-Phase Properties of Catalysed Neutron Stars. https://doi.org/10.1140/epja/s10050-021-00605-6

  32. [43]

    Dinh Thi, T

    H. Dinh Thi, T. Carreau, A. F. Fantina and F. Gul- minelli, Astron. Astrophys.654, A77 (2020). Uncertainties in the Pasta-Phase Properties of Catalysed Neutron Stars. https://doi.org/10.1051/0004-6361/202141192

  33. [44]

    L. E. Balliet, W. G. Newton, S. Cantu and S. Budimir, As- trophys. J.918, 79 (2021). Prior Probability Distributions of Neutron Star Crust Models. https://doi.org/10.3847/1538- 4357/ac06a4

  34. [45]

    N. N. Shchechilin, N. A. Zemlyakov, A. I. Chugunov and M.E. Gusakov, Symmetry8, 582 (2022). Pasta Phases in Neutron Star Mantle: Extended Thomas–Fermi vs. Compressible Liquid Drop Approaches. https://doi.org/10.3390/universe8110582

  35. [46]

    D. G. Ravenhall, C. J. Pethick, J. M. Lattimer, Nucl. Phys. A, 407, 571 (1983). Nuclear Interface Energy at Finite Tempera- tures. https://doi.org/10.1016/0375-9474(83)90667-X

  36. [47]

    C. P. Lorentz, D. G. Ranvenhall and C. J. Pethick, Phys. Rev. Lett.70, 379 (1993). Neutron Star Crusts. https://doi.org/10.1103/PhysRevLett.70.379

  37. [48]

    Brack, C

    M. Brack, C. Guet and H.-B. H ˚akanson, Phys. Rep.123, 275 (1985). Selfconsistent Semiclassical Description of Average Nuclear Properties - a Link Between Microscopic and Macro- scopic Models. https://doi.org/10.1016/0370-1573(86)90078-5

  38. [49]

    Centeles, M

    M. Centeles, M. Del Estal and M. Vinas, Nucl. Phys. A635, 193 (1998). Semiclassical Treatment of Asymmetric Semi-infinite Nuclear Matter: Surface and Curvature Properties in Relativis- tic and Non-relativistic Models. https://doi.org/10.1016/S0375- 9474(98)00167-5

  39. [50]

    Douchin, P

    F. Douchin, P. Haensel and J. Meyer, Nucl. Phys. A,665, 419 (2000). Nuclear Surface and Curvature Properties for SLy Skyrme Forces and Nuclei in the Inner Neutron-Star Crust. https://doi.org/10.1016/S0375-9474(99)00397-8

  40. [51]

    Kobyakov and C

    D. Kobyakov and C. J. Pethick, MNRAS Lett.449, L110— L112 (2015). Elastic Properties of Polycrystalline Dense Mat- ter. https://doi.org/10.1093/mnrasl/slv027

  41. [52]

    C. J. Pethick and A. Y . Potekhin, Phys. Lett. B427, 7 (1998). Liquid Crystals in the Mantles of Neutron Stars. https://doi.org/10.1016/s0370-2693(98)00341-4

  42. [53]

    D. N. Kobyakov and C. J. Pethick, Sov. Phys. JETP127, 851 (2018). Superfluid Liquid Crystals: Pasta Phases in Neutron Star Crusts. https://doi.org/10.1134/s1063776118110067

  43. [54]

    G. Baym, C. Pethick and D. Pines, Nature (Lon- don),224, 673 (1969). Superfluidity in Neutron Stars. https://doi.org/10.1038/224673a0

  44. [56]

    Zhang and C

    Z.-W. Zhang and C. J. Pethick, Phys. Rev. C103, 055807 (2021). Proton Superconductivity in Pasta Phases in Neutron Star Crusts. https://doi.org/10.1103/physrevc.103.055807

  45. [57]

    Kobyakov, (2024).Unpublished

    D. Kobyakov, (2024).Unpublished

  46. [58]

    J. B. Ketterson and S. N. Song,Superconductivity, (Cambridge University Press, 1999)

  47. [59]

    J. R. Clem, Supercond. Sci. Technol.5, S33 (1992). Funda- mentals of vortices in the high-temperature superconductors. https://doi.org/10.1088/0953-2048/5/1S/006

  48. [60]

    R. A. Klemm, A. Luther and M. R. Beasley, Phys. Rev. B12, 877 (1975). Theory of the Upper Critical Field in Layered Su- perconductors. https://doi.org/10.1103/PhysRevB.12.877

  49. [61]

    Deutscher and O

    G. Deutscher and O. Entin-Wohlman, Phys. Rev. B17, 1249 (1978). Critical Fields of Weakly Coupled Superconductors. https://doi.org/10.1103/PhysRevB.17.1249

  50. [62]

    L. Wang, H. S. Lim and C. K. Ong, Supercond. Sci. Tech- nol.14, 252 (2001). Continuous Ginzburg-Landau Descrip- tion of Layered Superconductors. https://doi.org/10.1088/0953- 2048/14/5/305

  51. [63]

    Lim and J

    Y . Lim and J. W. Holt, Phys. Rev. C103, 025807 (2021). Proton Pairing in Neutron Stars from Chiral Effective Field Theory. https://doi.org/10.1103/physrevc.103.025807

  52. [64]

    Del Estal, M

    M. Del Estal, M. Centelles, X. Vinas, and S. K. Patra, Phys. Rev. C63, 044321 (2001). Pairing properties in relativis- tic mean field models obtained from effective field theory. https://doi.org/10.1103/PhysRevC.63.044321

  53. [65]

    Nucl. Phys. A811, 127 (2008). The microscopic pairing gap in a slab of nuclear matter for the Argonne v18 NN-potential. https://doi.org/10.1016/j.nuclphysa.2008.07.002

  54. [66]

    Schuck and X

    P. Schuck and X. Vinas, Phys. Rev. Lett.107, 205301 (2011). Suppression of Superfluidity upon Overflow of Trapped Fermions: Quantal and Thomas-Fermi Studies. https://doi.org/10.1103/PhysRevLett.107.205301

  55. [67]

    Minoru Okamoto, Toshiki Maruyama, Kazuhiro Ya- bana, and Toshitaka Tatsumi, Phys. Rev. C88, 025801 (2013). Nuclear “pasta” structures in low-density nu- clear matter and properties of the neutron-star crust. https://doi.org/10.1103/PhysRevC.88.025801

  56. [68]

    U. J. Furtado, S. S. Avancini and J. R. Marinelli, J. Phys. G: Nucl. Part. Phys.49, 025202 (2022). Pairing effects in nuclear pasta phase within the relativistic Thomas–Fermi formalism. https://doi.org/10.1088/1361-6471/ac3c33

  57. [69]

    B. K. Sharma,private communication(2024)

  58. [70]

    C. J. Pethick, D. G. Ravenhall, and C. P. Lorenz, Nucl. Phys. A584, 675 (1995). The inner boundary of a neutron star crust. https://doi.org/10.1016/0375-9474(94)00506-I

  59. [71]

    Keller, K

    J. Keller, K. Hebeler, C. J. Pethick, and A. Schwenk, Phys. Rev. Lett.132, 232701 (2024). Neutron Star Matter as a Dilute Solution of Protons in Neutrons. https://doi.org/10.1103/PhysRevLett.132.232701

  60. [72]

    P. M. Pizzochero, L. Viverit and R. A. Broglia, Phys. Rev. Lett79, 3347 (1997). V ortex-Nucleus Interaction and Pinning Forces in Neutron Stars. https://doi.org/10.1103/PhysRevLett.79.3347

  61. [73]

    Ian Easson and C. J. Pethick, Phys. Rev. D16, 275 (1977). Stress Tensor of Cosmic and Laboratory Type-II Superconduc- tors. https://doi.org/10.1103/PhysRevD.16.275

  62. [74]

    A. D. Sedrakian and D. M. Sedrakian, Astrophys. J.447, 305 (1995). Superfluid Core Rotation in Pulsars. I. V ortex Cluster Dynamics. https://doi.org/10.1086/175876

  63. [75]

    D. N. Kobyakov, C. J. Pethick, S. Reddy, A. Schwenk, Phys. Rev. C96, 025805 (2017). Dispersion and Decay of Collective Modes in Neutron Star Cores. https://doi.org/10.1103/physrevc.96.025805

  64. [76]

    Chabanat, P

    E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Schaeffer, Nucl. Phys. A627, 710 (1997). A Skyrme parametrization from subnuclear to neutron star densities. https://doi.org/10.1016/S0375-9474(97)00596-4

  65. [77]

    D. N. Kobyakov, Phys. Rev. C108, 049901 (2023). Erratum: Application of superconducting-superfluid magnetohydrodynamics to nuclear “pasta” in neu- tron stars [Phys. Rev. C98, 045803 (2018)]. https://doi.org/10.1103/PhysRevC.108.049901

  66. [78]

    C. A. van Eysden and A. Melatos, MNRAS409, 1253 (2010). Pulsar glitch recovery and the superfluidity coeffi- cients of bulk nuclear matter. https://doi.org/10.1111/j.1365- 25 2966.2010.17387.x

  67. [79]

    Dutra, O

    M. Dutra, O. Lourenco, J. S. Sa Martins, A. Delfino, J. R. Stone and P. D. Stevenson, Phys. Rev. C85, 035201 (2012). Skyrme Interaction and Nuclear Matter Constraints. https://doi.org/10.1103/PhysRevC.85.035201

  68. [80]

    Gonzalez-Boquera, Ph.D

    C. Gonzalez-Boquera, Ph.D. thesis, Universitat de Barcelona, September 2019; arXiv:2003.00490. Neutron- rich Matter in Atomic Nuclei and Neutron Stars. https://doi.org/10.48550/arXiv.2003.00490

  69. [81]

    Grams, R

    G. Grams, R. Somasundaram, J. Margueron and S. Reddy, Phys. Rev. C105, 035806 (2022). Properties of the Neutron Star Crust: Quantifying and Corre- lating Uncertainties with Improved Nuclear Physics. https://doi.org/10.1103/PhysRevC.105.035806

  70. [82]

    Grams, J

    G. Grams, J. Margueron, R. Somasundaram and S. Reddy, Eur. Phys. J. A58, 56 (2022). Confronting a Set of Skyrme andχ EFT Predictions for the Crust of Neutron Stars. https://doi.org/10.1140/epja/s10050-022-00706-w

  71. [83]

    Vinas, M

    X. Vinas, M. Centelles, X. Roca-Maza, and M. Warda, Eur. Phys. J. A50, 27 (2014). Density Dependence of the Sym- metry Energy from Neutron Skin Thickness in Finite Nuclei. https://doi.org/10.1140/epja/i2014-14027-8

  72. [84]

    C. J. Pethick and D. G. Ravenhall, Annu. Rev. Nucl. Part. Sci.45, 429 (1995). Matter at Large Neu- tron Excess and the Physics of Neutron-Star Crusts. https://doi.org/10.1146/annurev.ns.45.120195.002241

  73. [85]

    Richard Witmer, Ars magna or the rules of algebra, (Constable Dover, 1993)

    Girolamo Cardano, translated and edited by T. Richard Witmer, Ars magna or the rules of algebra, (Constable Dover, 1993)

  74. [86]

    Bonche, S

    P. Bonche, S. Levit and D. Vautherin, Nucl. Phys. A,436, 265 (1985). Statistical Properties and Stability of Hot Nuclei. https://doi.org/10.1016/0375-9474(85)90199-X

  75. [87]

    Suraud, Nucl

    E. Suraud, Nucl. Phys. A462, 109 (1987). Semi-Classical Calculations of Hot Nuclei. https://doi.org/10.1016/0375- 9474(87)90382-4

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.