pith. sign in

arxiv: 2501.08075 · v1 · pith:BY6PEYEBnew · submitted 2025-01-14 · ⚛️ nucl-th · astro-ph.HE

Pressure and chemical potentials in the inner crust of a cold neutron star within Hartree-Fock and extended Thomas-Fermi methods

Nicolas Chamel , Nikolai N. Shchechilin , Andrey I. Chugunov This is my paper

Pith reviewed 2026-05-23 05:17 UTC · model grok-4.3

classification ⚛️ nucl-th astro-ph.HE
keywords neutron star inner crustequation of statechemical potentialspressureHartree-Fockextended Thomas-FermiSkyrme interactionsaccreted crust
0
0 comments X

The pith

Exact formulas for chemical potentials and pressure allow consistent full equation of state calculation for the inner crust of neutron stars in mean-field methods.

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

The paper derives exact expressions for the chemical potentials and the pressure directly from the energy functional in self-consistent Hartree-Fock and extended Thomas-Fermi calculations with Skyrme interactions. These formulas replace numerical interpolation of energies computed at a few selected densities, thereby removing systematic errors when building the global equation of state. The same expressions apply to both catalyzed and accreted crusts once the appropriate equilibrium conditions for the composition are imposed. Numerical illustrations are given for the BSk24 model together with the resulting adiabatic index.

Core claim

Exact formulas for the chemical potentials and for the pressure can be derived from the energy functional within the Hartree-Fock and extended Thomas-Fermi approximations, so that the full equation of state of the inner crust is obtained consistently inside the same framework for both catalyzed and accreted matter.

What carries the argument

Exact analytic expressions for chemical potentials and pressure obtained directly from the energy functional of Skyrme-type interactions in the Hartree-Fock plus extended Thomas-Fermi treatment of nuclear clusters coexisting with free neutrons and electrons.

If this is right

  • The full equation of state follows directly from the energy functional without post-hoc numerical interpolation.
  • The same formulas apply to both catalyzed and accreted crusts once the correct composition conditions are chosen.
  • Refined tabulations of the BSk24 equation of state and the associated adiabatic index become available for the inner crust.
  • Global neutron-star structure and dynamical evolution calculations can be performed with the same mean-field framework used for the composition.

Where Pith is reading between the lines

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

  • The method could be applied to other density functionals or semiclassical approximations used for neutron-star matter.
  • Improved crust equations of state may affect integrated quantities such as the total crust mass or the moment of inertia.
  • The approach supplies a consistent starting point for later calculations of transport or superfluid properties that also rely on the same energy functional.

Load-bearing premise

The Skyrme effective interactions together with the Hartree-Fock and extended Thomas-Fermi approximations remain thermodynamically consistent when applied to the inhomogeneous inner-crust matter.

What would settle it

Implementation of the derived formulas yields pressure or chemical potentials that violate the relation between the energy density, pressure, and chemical potentials expected from the variational principle or from the Gibbs-Duhem relation.

Figures

Figures reproduced from arXiv: 2501.08075 by Andrey I. Chugunov, Nicolas Chamel, Nikolai N. Shchechilin.

Figure 1
Figure 1. Figure 1: FIG. 1: Panel (a): relative deviations between the following formulas for the pressure and the numerical evaluation of Eq. (9) [PITH_FULL_IMAGE:figures/full_fig_p017_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Adiabatic index in the inner crust (black and red curves) and core (blue curve) of a nonaccreted neutron star as a [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Relative deviation in the adiabatic index in the inner crust of a nonaccreted neutron star as a function of the average [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
read the original abstract

Self-consistent mean-field methods with Skyrme-type effective interactions and semiclassical approximations, such as the Thomas-Fermi approach and its extensions are particularly well-suited for describing in a thermodynamically consistent way the various phases of the dense matter present in the interior of neutron stars. These methods have been applied to predict the composition of the different regions, including the inner crust constituted by nuclear clusters coexisting with free neutrons and electrons. Because of the computational cost, the energy is typically calculated for a few selected average baryon number densities, and the results are interpolated to obtain the pressure numerically. However, this may introduce systematic errors in the calculations of the global structure of a neutron star and its dynamical evolution. In this paper, we show how the full equation of state can be consistently calculated within the same framework by deriving exact formulas for the chemical potentials and for the pressure that can be easily implemented in existing computer codes. These formulas are applicable to both catalyzed and accreted crusts. We discuss in each case the suitable conditions to impose to determine the composition. Numerical examples are also presented and discussed. Results from refined calculations of the BSk24 equation of state for the inner crust of nonaccreted neutron stars and the corresponding adiabatic index are provided.

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

1 major / 2 minor

Summary. The paper claims to derive exact formulas for the chemical potentials and pressure from the energy functional in the Hartree-Fock plus extended Thomas-Fermi framework with Skyrme interactions. These expressions are presented as directly implementable in existing codes to obtain the full equation of state for the inner crust without numerical interpolation, and are stated to apply to both catalyzed and accreted matter; numerical results are shown for the BSk24 interaction including the adiabatic index.

Significance. If the central derivations hold exactly, the work would remove a source of systematic error in neutron-star structure calculations by ensuring thermodynamic consistency between the minimized energy and the derived P and μ within the same mean-field approximation. The explicit treatment of both catalyzed and accreted cases and the provision of concrete numerical examples for BSk24 add practical value.

major comments (1)
  1. [§3] §3 (formulas for chemical potentials and pressure): the claim that the expressions are exact and free of additional surface corrections rests on the assumption that the ETF variational equations plus Skyrme effective-mass and gradient terms preserve the thermodynamic identities inside a finite Wigner-Seitz cell with the chosen boundary conditions. No explicit numerical verification is shown that the derived pressure equals −dE/dV (or the equivalent grand-potential expression) to machine precision for the same functional; this verification is load-bearing for the central claim of exact consistency.
minor comments (2)
  1. [Abstract] The abstract refers to 'refined calculations' of the BSk24 EOS but does not state what refinements (e.g., mesh size, convergence criteria, or additional terms) distinguish the present results from earlier BSk24 tabulations.
  2. [Formalism section] Notation for the cell volume V and the averaging procedure over the Wigner-Seitz cell should be introduced once and used consistently when the pressure formula is written in integrated form.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We respond to the single major comment below, agreeing that additional verification would strengthen the central claim while defending the analytical basis of the derivations.

read point-by-point responses

  1. Referee: §3 (formulas for chemical potentials and pressure): the claim that the expressions are exact and free of additional surface corrections rests on the assumption that the ETF variational equations plus Skyrme effective-mass and gradient terms preserve the thermodynamic identities inside a finite Wigner-Seitz cell with the chosen boundary conditions. No explicit numerical verification is shown that the derived pressure equals −dE/dV (or the equivalent grand-potential expression) to machine precision for the same functional; this verification is load-bearing for the central claim of exact consistency.

    Authors: The expressions in §3 are obtained by differentiating the total energy functional while enforcing the Euler-Lagrange equations from the ETF minimization (including effective-mass and gradient contributions) together with the chosen Wigner-Seitz boundary conditions. This procedure guarantees thermodynamic consistency by construction within the mean-field approximation, without extra surface terms, because any variation that would produce such corrections is already set to zero by the variational solution. The same logic applies to both catalyzed and accreted matter. We nevertheless acknowledge that an explicit numerical check is not presented in the current version. In the revised manuscript we will add a short verification (e.g., in an appendix) showing that the analytically derived pressure agrees with −dE/dV to machine precision for representative BSk24 densities. revision: yes

Circularity Check

0 steps flagged

Derivation of P and μ from energy functional is self-contained

full rationale

The paper derives exact formulas for chemical potentials and pressure directly from the energy functional using standard thermodynamic relations within the Hartree-Fock and extended Thomas-Fermi approximations. These formulas are presented as generally applicable to the inhomogeneous matter in the inner crust for both catalyzed and accreted cases, with numerical examples provided. No quoted step reduces a prediction to a fitted input by construction, no self-citation is load-bearing for the central claim, and the derivation does not rely on uniqueness theorems or ansatzes imported from prior self-work. The result is independent of the specific Skyrme parametrization chosen and remains falsifiable via direct comparison to numerical derivatives of the energy.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard thermodynamic identities applied to mean-field energy functionals; Skyrme parameters are taken from prior literature (BSk24). No new free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption Thermodynamic relations (pressure and chemical potentials derivable from the energy density functional) hold exactly under the Hartree-Fock and extended Thomas-Fermi approximations for the inner crust.
    Invoked to justify the exact formulas replacing numerical interpolation.

pith-pipeline@v0.9.0 · 5771 in / 1269 out tokens · 34678 ms · 2026-05-23T05:17:28.986005+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

60 extracted references · 60 canonical work pages

  1. [1]

    dR(q) α dr d dr ∂R(q) α ∂R + ℓα(ℓα + 1) r2 R(q) α (r) ∂R(q) α ∂R # = 1 2π X α(q) g(q) α Z R 0 dr

    Hartree-Fock theory Adopting the spherical Wigner-Seitz approximation, the pressure formula (83) in the HF theory reads PHF = −EHF(R) − 1 R2 Z R 0 dr r2X q Uq(r) ∂nq(r) ∂R + ¯h2 2m⊕q (r) ∂τq(r) ∂R + Wq(r) ∂Jq(r) ∂R − 1 R2 Z R 0 dr r2Ue(r) ∂ne(r) ∂R . (D1) Let us examine the first term in the integral over nucleon fields: Z R 0 dr r2Uq(r) ∂nq(r) ∂R = 1 2π ...

    work page

  2. [2]

    ∂nq(r) ∂R − nΛq ∂fq(r) ∂R xqxqxq # = X q ˜µq

    Extended Thomas-Fermi approach with parametrized profiles Quite generally, expressing the profiles as in Eq. (69), we have ∂nq(r) ∂R = ∂nq(r) ∂nBq ∂nBq ∂R + ∂nq(r) ∂xqxqxq · ∂xqxqxq ∂R + ∂nq(r) ∂nΛq ∂nΛq ∂R + nΛq ∂fq(r) ∂R xqxqxq , (D10) where the last term arises from the explicit dependence of the parametrized distributions on R. Let us analyze separate...

    work page

  3. [3]

    Blaschke and N

    D. Blaschke and N. Chamel, in Astrophysics 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. 337

    work page 2018

  4. [4]

    Chamel and A

    N. Chamel and A. F. Fantina, Phys. Rev. C 94, 065802 (2016)

    work page 2016

  5. [5]

    B. K. Harrison, M. Wakano, and J. A. Wheeler, in Onzi` eme Conseil de Physique Solvay, Stoops, Brussels (1958), p. 124

    work page 1958

  6. [6]

    R. N. Wolf, D. Beck, K. Blaum, C. B¨ ohm, C. Borgmann, M. Breitenfeldt, N. Chamel, S. Goriely, F. Herfurth, M. Kowalska, et al., Phys. Rev. Lett. 110, 041101 (2013)

    work page 2013

  7. [7]

    Chamel, Phys

    N. Chamel, Phys. Rev. C 101, 032801 (2020)

    work page 2020

  8. [8]

    Chamel and A

    N. Chamel and A. F. Fantina, Phys. Rev. D 93, 063001 (2016)

    work page 2016

  9. [9]

    J. M. Pearson, S. Goriely, and N. Chamel, Phys. Rev. C 83, 065810 (2011)

    work page 2011

  10. [10]

    Chamel, A

    N. Chamel, A. F. Fantina, J. L. Zdunik, and P. Haensel, Phys. Rev. C 91, 055803 (2015)

    work page 2015

  11. [11]

    Bender, P.-H

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

    work page 2003

  12. [12]

    J. W. Negele and D. Vautherin, Nuclear Physics A 207, 298 (1973)

    work page 1973

  13. [13]

    Baldo, E

    M. Baldo, E. E. Saperstein, and S. V. Tolokonnikov, Nucl. Phys. A 775, 235 (2006)

    work page 2006

  14. [14]

    G¨ ogelein and H

    P. G¨ ogelein and H. M¨ uther, Phys. Rev. C76, 024312 (2007)

    work page 2007

  15. [15]

    Grill, J

    F. Grill, J. Margueron, and N. Sandulescu, Phys. Rev. C 84, 065801 (2011)

    work page 2011

  16. [16]

    P. N. McDermott, H. M. van Horn, and C. J. Hansen, Astrophys. J. 325, 725 (1988)

    work page 1988

  17. [17]

    B. K. Harrison, K. S. Thorne, M. Wakano, and J. A. Wheeler, Gravitation Theory and Gravitational Collapse (Chicago: University of Chicago Press, 1965)

    work page 1965

  18. [18]

    D. A. Kirzhnits, Field theoretical methods in many body systems (Pergamon, Oxford, 1967)

    work page 1967

  19. [19]

    Grammaticos and A

    B. Grammaticos and A. Voros, Annals of Physics 123, 359 (1979)

    work page 1979

  20. [20]

    Grammaticos and A

    B. Grammaticos and A. Voros, Annals of Physics 129, 153 (1980)

    work page 1980

  21. [21]

    Brack, C

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

    work page 1985

  22. [22]

    Shelley and A

    M. Shelley and A. Pastore, Phys. Rev. C 103, 035807 (2021)

    work page 2021

  23. [23]

    J. M. Pearson, N. Chamel, A. Pastore, and S. Goriely, Phys. Rev. C 91, 018801 (2015)

    work page 2015

  24. [24]

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

    work page 2012

  25. [25]

    J. M. Pearson, N. Chamel, A. Y. Potekhin, A. F. Fantina, C. Ducoin, A. K. Dutta, and S. Goriely, Mon. Not. R. Astron. Soc. 481, 2994 (2018). 25

    work page 2018

  26. [26]

    Ferreira and C

    M. Ferreira and C. Providˆ encia, Universe6, 220 (2020)

    work page 2020

  27. [27]

    Suleiman, M

    L. Suleiman, M. Fortin, J. L. Zdunik, and P. Haensel, Phys. Rev. C 104, 015801 (2021)

    work page 2021

  28. [28]

    Chamel, J

    N. Chamel, J. M. Pearson, and N. N. Shchechilin, Phys. Rev. C 110, 045808 (2024), URL https://link.aps.org/doi/ 10.1103/PhysRevC.110.045808

    work page doi:10.1103/physrevc.110.045808 2024

  29. [29]

    Okamoto, T

    M. Okamoto, T. Maruyama, K. Yabana, and T. Tatsumi, Phys. Rev. C 88, 025801 (2013)

    work page 2013

  30. [30]

    Kittel, Introduction to solid state physics (John Wiley and Sons, 1996)

    C. Kittel, Introduction to solid state physics (John Wiley and Sons, 1996)

    work page 1996

  31. [31]

    M. E. Gusakov and A. I. Chugunov, Phys. Rev. Lett. 124, 191101 (2020)

    work page 2020

  32. [32]

    M. E. Gusakov, E. M. Kantor, and A. I. Chugunov, Phys. Rev. D 104, L081301 (2021)

    work page 2021

  33. [33]

    M. E. Gusakov and A. I. Chugunov, Phys. Rev. D 109, 123032 (2024), 2406.05680

    work page arXiv 2024

  34. [34]

    Maruyama, T

    T. Maruyama, T. Tatsumi, D. N. Voskresensky, T. Tanigawa, and S. Chiba, Phys. Rev. C 72, 015802 (2005), URL https://link.aps.org/doi/10.1103/PhysRevC.72.015802

    work page doi:10.1103/physrevc.72.015802 2005

  35. [35]

    E. E. Salpeter, Astrophys. J. 134, 669 (1961)

    work page 1961

  36. [36]

    Kohn and L

    W. Kohn and L. J. Sham, Physical Review 140, 1133 (1965)

    work page 1965

  37. [37]

    Wigner and F

    E. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933)

    work page 1933

  38. [38]

    Vautherin and D

    D. Vautherin and D. M. Brink, Phys. Rev. C 5, 626 (1972)

    work page 1972

  39. [39]

    P. M. Pizzochero, F. Barranco, E. Vigezzi, and R. A. Broglia, Astrophys. J. 569, 381 (2002)

    work page 2002

  40. [40]

    Barranco, R

    F. Barranco, R. A. Broglia, and E. Vigezzi, Journal of Physics G Nuclear Physics 37, 064023 (2010)

    work page 2010

  41. [41]

    Baroni, A

    S. Baroni, A. Pastore, F. Raimondi, F. Barranco, R. A. Broglia, and E. Vigezzi, Phys. Rev. C 82, 015807 (2010)

    work page 2010

  42. [42]

    Chamel, S

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

    work page 2007

  43. [43]

    Montani, C

    F. Montani, C. May, and H. M¨ uther, Phys. Rev. C69, 065801 (2004)

    work page 2004

  44. [44]

    Pastore, S

    A. Pastore, S. Baroni, and C. Losa, Phys. Rev. C84, 065807 (2011), URL https://link.aps.org/doi/10.1103/PhysRevC. 84.065807

    work page doi:10.1103/physrevc 2011

  45. [45]

    Pastore, M

    A. Pastore, M. Shelley, S. Baroni, and C. A. Diget, Journal of Physics G Nuclear Physics 44, 094003 (2017)

    work page 2017

  46. [46]

    E. H. Lieb and B. Simon, Phys. Rev. Lett. 31, 681 (1973)

    work page 1973

  47. [47]

    S. L. Shapiro and S. A. Teukolsky, Black holes, white dwarfs, and neutron stars: The physics of compact objects (Cornell University, Ithaca, New York, 1983), ISBN 9780471873167

    work page 1983

  48. [48]

    A. M. Abrahams and S. L. Shapiro, Phys. Rev. A 42, 2530 (1990)

    work page 1990

  49. [49]

    D. G. Yakovlev, L. R. Gasques, A. V. Afanasjev, M. Beard, and M. Wiescher, Phys. Rev. C 74, 035803 (2006)

    work page 2006

  50. [50]

    N. N. Shchechilin, M. E. Gusakov, and A. I. Chugunov, arXiv e-prints arXiv:2202.07486 (2022)

    work page arXiv 2022

  51. [51]

    N. N. Shchechilin, M. E. Gusakov, and A. I. Chugunov, Monthly Notices of the Royal Astronomical Society 523, 4643 (2023), ISSN 0035-8711, https://academic.oup.com/mnras/article-pdf/523/3/4643/50667204/stad1731.pdf, URL https: //doi.org/10.1093/mnras/stad1731

    work page doi:10.1093/mnras/stad1731 2023

  52. [52]

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

    work page 2008

  53. [53]

    N. N. Shchechilin, N. Chamel, J. M. Pearson, A. I. Chugunov, and A. Y. Potekhin, Phys. Rev. C 109, 055802 (2024), URL https://link.aps.org/doi/10.1103/PhysRevC.109.055802

    work page doi:10.1103/physrevc.109.055802 2024

  54. [54]

    Goriely, N

    S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. C 88, 024308 (2013)

    work page 2013

  55. [55]

    J. M. Pearson and N. Chamel, Phys. Rev. C 105, 015803 (2022)

    work page 2022

  56. [56]

    N. N. Shchechilin, N. Chamel, and J. M. Pearson, Phys. Rev. C 108, 025805 (2023)

    work page 2023

  57. [57]

    R. C. Tolman, Physical Review 55, 364 (1939)

    work page 1939

  58. [58]

    J. R. Oppenheimer and G. M. Volkoff, Physical Review 55, 374 (1939)

    work page 1939

  59. [59]

    Huxford, R

    R. Huxford, R. Kashyap, S. Borhanian, A. Dhani, I. Gupta, and B. S. Sathyaprakash, Phys. Rev. D 109, 103035 (2024), URL https://link.aps.org/doi/10.1103/PhysRevD.109.103035

    work page doi:10.1103/physrevd.109.103035 2024

  60. [60]

    G. Baym, H. A. Bethe, and C. J. Pethick, Nucl. Phys. A 175, 225 (1971)

    work page 1971