pith. sign in

arxiv: 2605.20934 · v1 · pith:JWMEXTBInew · submitted 2026-05-20 · 🌌 astro-ph.HE · nucl-th

Superfluid fraction in the crystalline crust of a neutron star: role of quantum zero-point motion of ions

Nicolas Chamel This is my paper

Pith reviewed 2026-05-21 03:27 UTC · model grok-4.3

classification 🌌 astro-ph.HE nucl-th
keywords neutron starsinner crustsuperfluid fractionzero-point motioncrystal latticeband structureeffective mass
0
0 comments X

The pith

The suppression of the neutron superfluid fraction in the inner crust of a neutron star persists despite quantum zero-point motion of the ions.

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

This paper examines how the quantum vibrations of ions around their lattice sites affect the neutron superfluid in the inner crust of cold neutron stars. The author performs detailed three-dimensional calculations for common crystal structures to see if these motions can restore the superfluid fraction that is otherwise reduced by the lattice. Even accounting for this motion, the superfluid fraction stays strongly suppressed in the middle density range, which in turn makes the ions heavier and damps their movements more. This self-consistent picture questions the common assumption that the lattice is perfectly rigid. The findings matter because they influence how neutron stars rotate and cool over time.

Core claim

In both body-centered cubic and face-centered cubic lattices, fully three-dimensional band-structure calculations in the weak-coupling approximation show that the neutron superfluid fraction remains strongly suppressed in the intermediate density region of the inner crust. Consequently, the effective mass of the ions increases substantially, further reducing ion fluctuations. These results are obtained by treating the effects of the superfluid and the ion motion self-consistently.

What carries the argument

Three-dimensional band-structure calculations of the superfluid fraction in the weak-coupling approximation for body-centered and face-centered cubic ion lattices.

If this is right

  • The presence of the neutron superfluid alters the dynamics of the crust.
  • The effective mass of the ions increases dramatically.
  • Ion fluctuations experience additional damping from the higher effective mass.
  • The results bear on the rotational and thermal evolution of neutron stars.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Models of neutron star cooling may need to adjust for the modified ion motion in the inner crust.
  • Rotational evolution calculations could incorporate the self-consistent increase in ion effective mass.

Load-bearing premise

The weak-coupling approximation remains valid for the band-structure calculations of the superfluid fraction across the density range of the inner crust.

What would settle it

A calculation or observation showing that the superfluid fraction recovers significantly in the intermediate density region for either bcc or fcc lattices would falsify the central claim.

read the original abstract

The suppression of the neutron superfluid fraction in the inner crust of a cold neutron star is mitigated by the quantum zero-point motion of ions about their equilibrium position. In turn, the crustal dynamics is altered by the presence of the neutron superfluid. These effects are studied self-consistently to assess the validity of the usual assumption of a perfect rigid lattice. To this end, fully three-dimensional band-structure calculations of the superfluid fraction are carried out in the weak-coupling approximation, considering body- and face-centered cubic lattices. In both cases, the superfluid fraction is still found to be strongly suppressed in the intermediate region of the inner crust. In turn, the effective mass of the ions is dramatically increased, thus further damping the ion fluctuations. These results are of relevance for the rotational and thermal evolutions of neutron stars.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript examines the role of quantum zero-point motion of ions in the inner crust of a neutron star and its effect on the neutron superfluid fraction. Using fully three-dimensional band-structure calculations in the weak-coupling approximation for both body-centered cubic (bcc) and face-centered cubic (fcc) lattices, the authors find that the superfluid fraction remains strongly suppressed in the intermediate density region. This suppression leads to a dramatic increase in the effective mass of the ions, which further damps ion fluctuations, with implications for the rotational and thermal evolution of neutron stars.

Significance. If the central results hold, the work would provide a useful quantitative constraint on crust entrainment and challenge the rigid-lattice assumption in neutron-star models. The self-consistent treatment of zero-point motion together with the use of 3D band-structure methods for two lattice geometries is a clear methodological strength.

major comments (2)
  1. [Weak-coupling section] Weak-coupling section (near Eq. for band-structure superfluid fraction): the central claim of strong suppression in the intermediate inner-crust region rests on the weak-coupling treatment remaining quantitatively accurate when the neutron Fermi energy, pairing gap, and lattice potential vary over the full density window; no independent check against strong-coupling or non-perturbative methods is described, leaving open an O(1) correction precisely where suppression is reported to be strongest.
  2. [Results for bcc and fcc lattices] Results for bcc and fcc lattices: the reported dramatic increase in ion effective mass follows directly from the suppressed superfluid fraction, but the self-consistent damping of fluctuations requires explicit demonstration that the zero-point motion does not invalidate the assumed lattice geometries or the entrainment coefficients used in the band-structure input.
minor comments (1)
  1. [Abstract] Abstract: the statement that calculations are performed 'in the weak-coupling approximation' would benefit from a one-sentence clarification of the precise regime (e.g., ratio of gap to Fermi energy) to help readers assess applicability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We respond to each major comment below and will revise the manuscript accordingly to address the points raised.

read point-by-point responses

  1. Referee: [Weak-coupling section] Weak-coupling section (near Eq. for band-structure superfluid fraction): the central claim of strong suppression in the intermediate inner-crust region rests on the weak-coupling treatment remaining quantitatively accurate when the neutron Fermi energy, pairing gap, and lattice potential vary over the full density window; no independent check against strong-coupling or non-perturbative methods is described, leaving open an O(1) correction precisely where suppression is reported to be strongest.

    Authors: We agree that the quantitative accuracy of the weak-coupling approximation merits explicit discussion, particularly in the intermediate-density region where the suppression is strongest. Our calculations are performed within the weak-coupling framework as stated, which is motivated by the small gap-to-Fermi-energy ratio throughout the inner crust. A full non-perturbative or strong-coupling treatment would require an entirely different computational approach and is beyond the scope of the present study. In the revised manuscript we will add a dedicated paragraph in the discussion section that quantifies the expected regime of validity, cites relevant strong-coupling calculations for uniform neutron matter, and notes the possible size of O(1) corrections. revision: partial

  2. Referee: [Results for bcc and fcc lattices] Results for bcc and fcc lattices: the reported dramatic increase in ion effective mass follows directly from the suppressed superfluid fraction, but the self-consistent damping of fluctuations requires explicit demonstration that the zero-point motion does not invalidate the assumed lattice geometries or the entrainment coefficients used in the band-structure input.

    Authors: The self-consistent loop in our work determines the ion effective mass from the computed superfluid fraction and then uses this mass to evaluate the zero-point fluctuation amplitude. We have verified internally that these amplitudes remain well below the lattice spacing (typically < 10 % of the nearest-neighbor distance) across the density range of interest, thereby preserving the assumed bcc and fcc geometries and the validity of the input entrainment coefficients. To make this demonstration explicit for the reader, we will add a short subsection (or appendix) that reports the root-mean-square displacements and confirms consistency with the rigid-lattice and band-structure assumptions. revision: yes

Circularity Check

0 steps flagged

No circularity: standard band-structure computations yield independent results

full rationale

The paper applies established three-dimensional band-structure methods in the weak-coupling limit to compute the neutron superfluid fraction for bcc and fcc lattices across inner-crust densities. The reported suppression and consequent ion effective-mass increase emerge directly as numerical outputs from these calculations rather than being defined in terms of themselves or obtained by fitting a parameter that is then relabeled as a prediction. No load-bearing self-citations, uniqueness theorems, or smuggled ansatzes are invoked; the derivation chain is self-contained and relies on standard techniques whose validity can be checked externally.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the weak-coupling approximation and the choice of bcc and fcc lattice geometries as representative of the inner crust; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption Weak-coupling approximation is adequate for computing the superfluid fraction throughout the inner crust density range.
    Invoked for the band-structure calculations described in the abstract.
  • domain assumption Body-centered cubic and face-centered cubic lattices adequately represent the crystalline structure in the inner crust.
    Explicitly considered in the calculations.

pith-pipeline@v0.9.0 · 5668 in / 1378 out tokens · 31204 ms · 2026-05-21T03:27:37.968003+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

40 extracted references · 40 canonical work pages

  1. [1]

    The volume of a primitive cell is therefore Ω cell =|a 1a1a1 ·a 2a2a2 ×a 3a3a3|=a 3/2

    Body-centered cubic lattice The primitive basis vectors of a body-centered cubic lattice are given by a1a1a1 = a 2 (−ˆxˆxˆx+ ˆyˆyˆy+ ˆzˆzˆz), a2a2a2 = a 2 (ˆxˆxˆx−ˆyˆyˆy+ ˆzˆzˆz), a3a3a3 = a 2 (ˆxˆxˆx+ ˆyˆyˆy−ˆzˆzˆz),(A1) where ˆxˆxˆx, ˆyˆyˆy, ˆzˆzˆzare the Cartesian unit vectors andais the size of the conventional cubic cell. The volume of a primitive ce...

    work page

  2. [2]

    Face-centered cubic lattice The corresponding primitive basis vectors of a face-centered cubic lattice are given by a1a1a1 = a 2 (ˆyˆyˆy+ ˆzˆzˆz), a2a2a2 = a 2 (ˆxˆxˆx+ ˆzˆzˆz), a3a3a3 = a 2 (ˆxˆxˆx+ ˆyˆyˆy).(A3) The volume of a primitive cell is therefore Ω cell =|a 1a1a1 ·a 2a2a2 ×a 3a3a3|=a 3/4. The reciprocal lattice is a body-centered cubic with basi...

    work page

  3. [3]

    Chamel, Journal of Astrophysics and Astronomy38, 43 (2017)

    N. Chamel, Journal of Astrophysics and Astronomy38, 43 (2017)

    work page 2017

  4. [4]

    Haskell and A

    B. Haskell and A. Sedrakian, inAstrophysics and Space Science Library, edited by L. Rezzolla, P. Pizzochero, D. I. Jones, N. Rea, and I. Vida˜ na (2018), vol. 457 ofAstrophysics and Space Science Library, p. 401

    work page 2018

  5. [5]

    Chamel, Phd thesis, Universit´ e Pierre et Marie Curie - Paris VI, France (2004), URL https://theses.hal.science/tel-00011130

    N. Chamel, Phd thesis, Universit´ e Pierre et Marie Curie - Paris VI, France (2004), URL https://theses.hal.science/tel-00011130

    work page 2004

  6. [6]

    Carter, N

    B. Carter, N. Chamel, and P. Haensel, Nucl. Phys. A748, 675 (2005)

    work page 2005

  7. [7]

    Chamel, Nucl

    N. Chamel, Nucl. Phys. A747, 109 (2005)

    work page 2005

  8. [8]

    Carter, N

    B. Carter, N. Chamel, and P. Haensel, Nuclear Physics A759, 441 (2005)

    work page 2005

  9. [9]

    Chamel, Journal of Low Temperature Physics189, 328 (2017)

    N. Chamel, Journal of Low Temperature Physics189, 328 (2017)

    work page 2017

  10. [10]

    Chauveau, C

    G. Chauveau, C. Maury, F. Rabec, C. Heintze, G. Brochier, S. Nascimbene, J. Dalibard, J. Beugnon, S. M. Roccuzzo, and S. Stringari, Phys. Rev. Lett.130, 226003 (2023), URL https://link.aps.org/doi/10.1103/PhysRevLett.130.226003

    work page doi:10.1103/physrevlett.130.226003 2023

  11. [11]

    J. Tao, M. Zhao, and I. B. Spielman, Phys. Rev. Lett.131, 163401 (2023), URLhttps: //link.aps.org/doi/10.1103/PhysRevLett.131.163401

    work page doi:10.1103/physrevlett.131.163401 2023

  12. [12]

    A. J. Leggett, Phys. Rev. Lett.25, 1543 (1970)

    work page 1970

  13. [13]

    A. J. Leggett, Journal of Statistical Physics93, 927 (1998)

    work page 1998

  14. [14]

    Kashiwaba and T

    Y. Kashiwaba and T. Nakatsukasa, Phys. Rev. C100, 035804 (2019)

    work page 2019

  15. [15]

    Sekizawa, S

    K. Sekizawa, S. Kobayashi, and M. Matsuo, Phys. Rev. C105, 045807 (2022)

    work page 2022

  16. [16]

    Minami and G

    Y. Minami and G. Watanabe, Physical Review Research4, 033141 (2022)

    work page 2022

  17. [17]

    Yoshimura and K

    K. Yoshimura and K. Sekizawa, Phys. Rev. C109, 065804 (2024)

    work page 2024

  18. [18]

    Almirante and M

    G. Almirante and M. Urban, Phys. Rev. C109, 045805 (2024)

    work page 2024

  19. [19]

    Almirante and M

    G. Almirante and M. Urban, Phys. Rev. C110, 065802 (2024), URLhttps://link.aps. org/doi/10.1103/PhysRevC.110.065802. 14

    work page doi:10.1103/physrevc.110.065802 2024

  20. [20]

    Chamel, Phys

    N. Chamel, Phys. Rev. C111, 045803 (2025)

    work page 2025

  21. [21]

    Kittel,Quantum Theory of Solids(John Wiley and Sons, 1966)

    C. Kittel,Quantum Theory of Solids(John Wiley and Sons, 1966)

    work page 1966

  22. [22]

    Kobyakov and C

    D. Kobyakov and C. J. Pethick, Phys. Rev. C87, 055803 (2013)

    work page 2013

  23. [23]

    Andersson, K

    N. Andersson, K. Glampedakis, W. C. G. Ho, and C. M. Espinoza, Phys. Rev. Lett.109, 241103 (2012)

    work page 2012

  24. [24]

    Chamel, Phys

    N. Chamel, Phys. Rev. Lett.110, 011101 (2013)

    work page 2013

  25. [25]

    Chamel, D

    N. Chamel, D. Page, and S. Reddy, Phys. Rev. C87, 035803 (2013)

    work page 2013

  26. [26]

    Delsate, N

    T. Delsate, N. Chamel, N. G¨ urlebeck, A. F. Fantina, J. M. Pearson, and C. Ducoin, Phys. Rev. D94, 023008 (2016)

    work page 2016

  27. [27]

    Allard and N

    V. Allard and N. Chamel, Phys. Rev. C108, 015801 (2023)

    work page 2023

  28. [28]

    Bender, P.-H

    M. Bender, P.-H. Heenen, and P.-G. Reinhard, Reviews of Modern Physics75, 121 (2003)

    work page 2003

  29. [29]

    C. J. Pethick, N. Chamel, and S. Reddy, Progress of Theoretical Physics Supplement186, 9 (2010)

    work page 2010

  30. [30]

    Haensel, A

    P. Haensel, A. Y. Potekhin, and D. G. Yakovlev,Neutron Stars. 1. Equation of State and Structure(Springer, New York, 2007)

    work page 2007

  31. [31]

    D. A. Baiko, Phys. Rev. E66, 056405 (2002), URLhttps://link.aps.org/doi/10.1103/ PhysRevE.66.056405

    work page 2002

  32. [32]

    Chamel and V

    N. Chamel and V. Allard, Phys. Rev. C100, 065801 (2019)

    work page 2019

  33. [33]

    Gilat and L

    G. Gilat and L. J. Raubenheimer, Physical Review144, 390 (1966)

    work page 1966

  34. [34]

    M. Onsi, A. K. Dutta, H. Chatri, S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. C 77, 065805 (2008)

    work page 2008

  35. [35]

    Brack, C

    M. Brack, C. Guet, and H. B. H˚ akansson, Phys. Rep.123, 275 (1985)

    work page 1985

  36. [36]

    Chamel, S

    N. Chamel, S. Naimi, E. Khan, and J. Margueron, Phys. Rev. C75, 055806 (2007)

    work page 2007

  37. [37]

    J. M. Pearson, N. Chamel, S. Goriely, and C. Ducoin, Phys. Rev. C85, 065803 (2012)

    work page 2012

  38. [38]

    D. A. Baiko, A. Y. Potekhin, and D. G. Yakovlev, Phys. Rev. E64, 057402 (2001)

    work page 2001

  39. [39]

    L. G. Cao, U. Lombardo, C. W. Shen, and N. V. Giai, Phys. Rev. C73, 014313 (2006)

    work page 2006

  40. [40]

    A. Y. Potekhin, D. A. Zyuzin, D. G. Yakovlev, M. V. Beznogov, and Y. A. Shibanov, Monthly Notices of the Royal Astronomical Society496, 5052 (2020), ISSN 0035-8711, https://academic.oup.com/mnras/article-pdf/496/4/5052/33507099/staa1871.pdf, URL https://doi.org/10.1093/mnras/staa1871. 15

    work page doi:10.1093/mnras/staa1871 2020