arxiv: 2605.05634 · v1 · submitted 2026-05-07 · ⚛️ nucl-th

Recognition: unknown

Relativistic mean-field study of the neutron star inner crust using the asymmetric finite difference method

Jinzhe Zhang , Hong Shen , Ying Zhang , Jinniu Hu

Authors on Pith no claims yet

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

classification ⚛️ nucl-th

keywords neutron star inner crustrelativistic mean-fieldWigner-Seitz approximationsymmetry energy slope Leffective nucleon massDirac equationsquantum shell effects

0 comments

The pith

Neutron star inner crust neutron radii and chemical potentials depend sensitively on the symmetry-energy slope L and effective nucleon mass.

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

This paper examines the ground-state properties of neutron-rich nuclear clusters in the inner crust of neutron stars by applying a relativistic mean-field model inside the Wigner-Seitz approximation. The radial Dirac equations are discretized with an asymmetric finite-difference scheme that maintains hermiticity and removes spurious states. Representative cells are computed with TM1-based forces that vary the symmetry-energy slope parameter L and with one parametrization that has a larger nucleon effective mass. Binding energy per nucleon falls steadily as L rises, and the larger effective mass produces an additional drop that grows stronger at higher densities. Quantum shell effects, missing from Thomas-Fermi treatments, create oscillatory density profiles that alter the neutron root-mean-square radius and chemical potential inside each cell.

Core claim

Within the relativistic mean-field model applied to representative Wigner-Seitz cells, the binding energy per nucleon decreases with larger symmetry-energy slope L, and a parametrization with higher nucleon effective mass leads to further decrease particularly at higher densities. The quantum shell effects absent in Thomas-Fermi lead to oscillatory density profiles that modify neutron properties. Consequently, the neutron root-mean-square radius and chemical potential in the cell are sensitive to L and the effective mass, highlighting their roles in the microscopic structure of the inner crust.

What carries the argument

The asymmetric finite-difference discretization of the radial Dirac equations inside the Wigner-Seitz cell, combined with TM1-based relativistic mean-field interactions that differ in the symmetry-energy slope L and in the effective nucleon mass.

If this is right

Binding energy per nucleon decreases systematically as the symmetry-energy slope L increases.
A larger effective nucleon mass produces an additional reduction in binding energy that becomes more pronounced at higher densities.
Quantum shell effects produce oscillatory neutron density distributions inside each cell.
The neutron root-mean-square radius and chemical potential inside the cell respond clearly to changes in both L and the effective mass.

Where Pith is reading between the lines

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

These parameter dependences could propagate into the overall equation of state of the crust and alter predictions for neutron-star radius or thermal evolution.
The method's removal of spurious states may improve reliability when shell effects dominate at densities near the crust-core transition.
Repeating the survey with other relativistic mean-field parametrizations would test whether the reported sensitivities are generic or specific to the TM1 family.

Load-bearing premise

The Wigner-Seitz approximation that replaces the full lattice with independent spherical cells remains valid for the chosen densities and interactions.

What would settle it

A full three-dimensional lattice calculation at the same average densities and with the same interactions would show whether the neutron root-mean-square radii and chemical potentials differ substantially from the spherical-cell results.

Figures

Figures reproduced from arXiv: 2605.05634 by Hong Shen, Jinniu Hu, Jinzhe Zhang, Ying Zhang.

**Figure 1.** Figure 1: FIG. 1 view at source ↗

**Figure 2.** Figure 2: FIG. 2 view at source ↗

read the original abstract

The ground-state properties of neutron-rich nuclear clusters in the inner crust of neutron stars are investigated within the Wigner-Seitz approximation using a relativistic mean-field framework. The radial Dirac equations are solved with an asymmetric finite-difference scheme, by which the hermiticity is preserved and spurious states are eliminated. Calculations are performed for representative Wigner-Seitz cells employing TM1-based interactions with different symmetry-energy slope parameters $L$, as well as a parametrization with a larger nucleon effective mass. It is found that the binding energy per nucleon decreases systematically with increasing $L$, while a larger effective mass leads to further reduction, particularly at higher densities. Quantum shell effects, which are absent in the Thomas-Fermi approximation, give rise to oscillatory density distributions and modify neutron properties. Within the Wigner-Seitz cell, the resulting neutron root-mean-square radius and chemical potential are shown to be sensitive to both $L$ and the effective nucleon mass, underscoring their important roles in determining the microscopic structure of the neutron-star inner crust.

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.

Desk Editor's Note private letter to a colleague

The asymmetric finite-difference solver is a practical technical step forward, and the paper maps clear trends in binding energy and neutron properties with L and effective mass, but the spherical Wigner-Seitz approximation is used without any error estimate or cross-check.

read the letter

The main takeaway is that this work gives a new asymmetric finite-difference scheme for the radial Dirac equations in relativistic mean-field calculations. The method keeps hermiticity and removes spurious states, which is a useful fix for this kind of problem. They apply it to Wigner-Seitz cells with TM1-based forces that vary the symmetry energy slope L and include one parametrization with larger nucleon effective mass. The results show binding energy per nucleon dropping steadily with higher L, and dropping more with the larger effective mass at higher densities. Shell effects produce oscillatory densities and shift the neutron root-mean-square radii and chemical potentials, with both quantities sensitive to L and effective mass. This is solid incremental work because it adds quantum shell corrections that Thomas-Fermi calculations miss and gives systematic trends that could matter for crust inputs in cooling or glitch models. The numerical implementation looks carefully done and the trends are presented directly. The softer part is the Wigner-Seitz spherical-cell assumption. The paper does all its calculations under isolated spherical cells with uniform background but does not quantify the error from this averaging, does not compare to cylindrical or slab geometries, and does not check how close the densities sit to the pasta transition where the approximation is known to be less reliable. If the geometry bias is comparable in size to the L-induced changes, the claimed sensitivities could be affected. This paper is for nuclear astrophysicists who build or use RMF crust models and want to see the effect of L and effective mass with shell corrections included. Readers working on numerical methods for Dirac equations or on neutron-star interior structure will find the concrete results useful. It shows clear, honest engagement with the equations and the literature, so it deserves a serious referee, though the review should ask for checks on the approximation's accuracy and some convergence or error analysis.

Referee Report

2 major / 2 minor

Summary. The paper investigates the ground-state properties of neutron-rich nuclear clusters in the neutron-star inner crust using relativistic mean-field theory within the Wigner-Seitz approximation. The radial Dirac equations are solved via an asymmetric finite-difference scheme that preserves hermiticity and eliminates spurious states. Employing TM1-based parametrizations with varying symmetry-energy slope L and a variant with larger nucleon effective mass, the work reports that binding energy per nucleon decreases systematically with increasing L and with larger effective mass (especially at higher densities). Quantum shell effects produce oscillatory density profiles absent in Thomas-Fermi treatments and modify neutron properties; specifically, the neutron root-mean-square radius and chemical potential inside the Wigner-Seitz cell are found to be sensitive to both L and the effective mass.

Significance. If the central results hold, the study usefully isolates the roles of the symmetry-energy slope L and the nucleon effective mass in shaping the microscopic structure of the inner crust, while demonstrating a numerically stable method for solving the Dirac equation in this setting. The explicit inclusion of quantum shell effects provides a clear advance over Thomas-Fermi calculations. The work is technically sound in its direct numerical approach and does not rely on circular fitting of parameters to the observables it predicts.

major comments (2)

[Wigner-Seitz cell calculations and results] The Wigner-Seitz approximation (employed throughout the calculations and results): all reported sensitivities of the neutron RMS radius and chemical potential to L and effective mass are obtained under the assumption of isolated spherical cells with uniform background. No error estimate, comparison to cylindrical or slab geometries, or assessment of multipole corrections is provided, even though the chosen densities approach the pasta transition where non-spherical structures become relevant. This unquantified systematic uncertainty directly affects the load-bearing claim that L and effective mass play important roles in determining the microscopic structure.
[Numerical method and results sections] Convergence and numerical validation: the abstract and results assert that the asymmetric finite-difference scheme eliminates spurious states and yields reliable neutron properties, yet no grid-size convergence tests, comparison of binding energies or radii against Thomas-Fermi solutions in the identical RMF framework, or error bars on the reported sensitivities are shown. Without these, the magnitude of the L- and m*-induced variations cannot be assessed relative to numerical or approximation errors.

minor comments (2)

[Abstract] The abstract introduces the 'asymmetric finite difference method' without a one-sentence description or citation to the original formulation; a brief parenthetical explanation would improve accessibility.
[Introduction] Notation for the effective mass and symmetry-energy slope L is used consistently but could be defined explicitly on first appearance in the main text for readers outside the immediate subfield.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the positive assessment of our work and the valuable suggestions for improvement. We address the major comments below and will revise the manuscript accordingly.

read point-by-point responses

Referee: The Wigner-Seitz approximation (employed throughout the calculations and results): all reported sensitivities of the neutron RMS radius and chemical potential to L and effective mass are obtained under the assumption of isolated spherical cells with uniform background. No error estimate, comparison to cylindrical or slab geometries, or assessment of multipole corrections is provided, even though the chosen densities approach the pasta transition where non-spherical structures become relevant. This unquantified systematic uncertainty directly affects the load-bearing claim that L and effective mass play important roles in determining the microscopic structure.

Authors: We agree that the spherical Wigner-Seitz approximation has limitations, especially as densities approach the transition to pasta phases. Our calculations are performed within this standard framework to focus on the effects of varying L and the effective mass in a controlled manner. We will revise the manuscript to include a more detailed discussion of the applicability of the Wigner-Seitz approximation, including references to studies on non-spherical geometries, and specify the density range where spherical cells are appropriate. However, performing calculations in cylindrical or slab geometries would require a substantial extension of our numerical method, which is beyond the scope of the present study. revision: partial
Referee: Convergence and numerical validation: the abstract and results assert that the asymmetric finite-difference scheme eliminates spurious states and yields reliable neutron properties, yet no grid-size convergence tests, comparison of binding energies or radii against Thomas-Fermi solutions in the identical RMF framework, or error bars on the reported sensitivities are shown. Without these, the magnitude of the L- and m*-induced variations cannot be assessed relative to numerical or approximation errors.

Authors: We appreciate this point. Although the asymmetric finite-difference method is designed to preserve hermiticity and eliminate spurious states, we did not present explicit convergence tests or comparisons in the original manuscript. In the revised version, we will add an appendix or subsection detailing the grid-size convergence, demonstrating stability of the results for the adopted mesh. We will also provide comparisons of binding energies and neutron radii with Thomas-Fermi calculations using the same RMF interactions to highlight the role of shell effects. These additions will allow readers to assess the numerical reliability and the magnitude of the reported sensitivities relative to approximation errors. revision: yes

standing simulated objections not resolved

Direct comparisons to cylindrical or slab geometries, multipole corrections, or quantitative error estimates from non-spherical structures.

Circularity Check

0 steps flagged

No circularity: direct numerical solution of RMF equations

full rationale

The paper solves the radial Dirac equations numerically inside fixed Wigner-Seitz cells using the asymmetric finite-difference method for given TM1-based RMF parametrizations (varied only in L and effective mass). The neutron RMS radius and chemical potential are computed outputs of this standard boundary-value problem; they are not fitted to the same observables, not defined in terms of themselves, and not obtained by renaming a prior result. The WS approximation is an input modeling choice whose validity is external to the derivation chain, not a self-referential step.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the relativistic mean-field Lagrangian, the Wigner-Seitz spherical-cell approximation, and the chosen TM1 parametrizations. No new particles or forces are introduced.

free parameters (2)

symmetry-energy slope L
Varied across several TM1-based interactions; its value is taken from prior fits to nuclear data and is not derived inside the paper.
nucleon effective mass
A larger value is tested as an additional parametrization; again taken from existing RMF sets.

axioms (2)

domain assumption Wigner-Seitz approximation replaces the periodic lattice with independent spherical cells
Invoked in the first sentence of the abstract; no error estimate is provided.
domain assumption Relativistic mean-field Lagrangian with TM1 parameters accurately describes neutron-rich matter at sub-saturation densities
Underlying framework; the paper tests variations but does not re-derive the Lagrangian.

pith-pipeline@v0.9.0 · 5482 in / 1660 out tokens · 38978 ms · 2026-05-08T04:23:26.696261+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

40 extracted references

[1]

Oertel, M

M. Oertel, M. Hempel, T. Klähn, and S. Typel, Rev. Mod. Phys. 89, 015007 (2017)

2017
[2]

Haensel, A

P. Haensel, A. Y. Potekhin, and D. G. Yakovlev, Neutron Stars 1: Equation of State and Structure, Astrophysics and Space Science Library, Vol. 326 (Springer New York, 2007)

2007
[3]

J. M. Lattimer and M. Prakash, Science 304, 536 (2004)

2004
[4]

Chamel and P

N. Chamel and P. Haensel, Living Reviews in Relativity 11, 10 (2008)

2008
[5]

Chamel, J

N. Chamel, J. M. Pearson, A. F. Fantina, C. Ducoin, S. Goriely, and A. Pastore, Acta Phys. Pol. B 46, 349 (2015)

2015
[6]

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

2011
[7]

Kreim, M

S. Kreim, M. Hempel, D. Lunney, and J. Schaffner-Bielich, International Journal of Mass Spectrometry 349, 63 (2013)

2013
[8]

Hashimoto, H

M. Hashimoto, H. Seki, and M. Yamada, Prog. Theor. Phys. 71, 320 (1984)

1984
[9]

Chamel, S

N. Chamel, S. Naimi, E. Khan, and J. Margueron, Physical Review C—Nuclear Physics 75, 055806 (2007)

2007
[10]

Chamel, Physical Review C 85, 035801 (2012)

N. Chamel, Physical Review C 85, 035801 (2012)

2012
[11]

Yoshimura and K

K. Yoshimura and K. Sekizawa, Physical Review C 109, 065804 (2024)

2024
[12]

Almirante, Particles 9, 22 (2026)

G. Almirante, Particles 9, 22 (2026)

2026
[13]

J. W. Negele and D. Vautherin, Nucl. Phys. A 207, 298 (1973)

1973
[14]

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

2012
[15]

Mondal, X

C. Mondal, X. Viñas, M. Centelles, and J. N. De, Phys. Rev. C 102, 015802 (2020)

2020
[16]

Ring, Prog

P. Ring, Prog. Part. Nucl. Phys. 37, 193 (1996)

1996
[17]

Nikšić, D

T. Nikšić, D. Vretenar, and P. Ring, Prog. Part. Nucl. Phys. 66, 519 (2011)

2011
[18]

J. Meng, H. Toki, S. G. Zhou, S. Q. Zhang, W. H. Long, and L. S. Geng, Prog. Part. Nucl. Phys. 57, 470 (2006)

2006
[19]

J. Cao, D. Yang, Z. Ma, and N. V. Giai, Chin. Phys. Lett. 25, 73 (2008)

2008
[20]

Boguta and J

J. Boguta and J. Rafelski, Phys. Lett. B 71, 22 (1977)

1977
[21]

H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi, Nucl. Phys. A 637, 435 (1998)

1998
[22]

N. N. Shchechilin, N. Chamel, J. M. Pearson, A. I. Chugunov, and A. Y. Potekhin, Phys. Rev. C 109, 055802 (2024) . 19

2024
[23]

C. J. Horowitz and B. D. Serot, Nuclear Physics A 368, 503 (1981)

1981
[24]

Y. K. Gambhir, P. Ring, and A. Thimet, Annals of Physics 198, 132 (1990)

1990
[25]

S.-G. Zhou, J. Meng, and P. Ring, Physical Review C 68, 034323 (2003)

2003
[26]

L.-S. Geng, J. Meng, and H. Toki, Chinese Physics Letters 24, 1865 (2007)

2007
[27]

B.-N. Lu, J. Zhao, E.-G. Zhao, and S.-G. Zhou, Physical Review C 89, 014323 (2014)

2014
[28]

Zhang, Y

Y. Zhang, Y. Bao, H. Shen, and J. Hu, Phys. Rev. C 106, L051303 (2022)

2022
[29]

Salomonson and P

S. Salomonson and P. Öster, Physical Review A 40, 5548 (1989)

1989
[30]

J. Y. Fang, S. W. Chen, and T. H. Heng, Nuclear Science and Techniques 31, 15 (2020)

2020
[31]

Zhao, Acta Physica Sinica 65, 052401 (2016) , (in Chinese)

B. Zhao, Acta Physica Sinica 65, 052401 (2016) , (in Chinese)

2016
[32]

Susskind, Physical Review D 16, 3031 (1977)

L. Susskind, Physical Review D 16, 3031 (1977)

1977
[33]

Stacey, Physical Review D 26, 468 (1982)

R. Stacey, Physical Review D 26, 468 (1982)

1982
[34]

S. S. Bao, J. N. Hu, Z. W. Zhang, and H. Shen, Phys. Rev. C 90, 045802 (2014)

2014
[35]

F. J. Fattoyev, C. J. Horowitz, J. Piekarewicz, and G. Shen, Physical Review C 82, 055803 (2010)

2010
[36]

Li, The European Physical Journal Special Topics , 1 (2026)

B.-A. Li, The European Physical Journal Special Topics , 1 (2026)

2026
[37]

S. Li, H. Shen, and J. Hu, Universe 11, 255 (2025)

2025
[38]

Sugahara and H

Y. Sugahara and H. Toki, Nuclear Physics A 579, 557 (1994)

1994
[39]

Sandulescu, Physical Review C 70, 025801 (2004)

N. Sandulescu, Physical Review C 70, 025801 (2004)

2004
[40]

Grill, J

F. Grill, J. Margueron, and N. Sandulescu, Physical Review C—Nuclear Physics 84, 065801 (2011). 20

2011