Three-dimensional orbital-free density functional theory description of nuclear pasta in the inner crust of neutron stars
Abstract
Background: In the bottom layer of the inner crust of neutron stars, various crystalline structures are expected to emerge that are collectively called ``nuclear pasta.'' It is desirable to know properties of nuclear pasta in a wide variety of conditions for astrophysical applications. However, three-dimensional fully-microscopic calculations require huge computational effort that makes it still challenging to carry out systematic calculations. Purpose: In this paper, we propose an efficient method to calculate various nuclear pasta configurations in a non-empirical manner, based on three-dimensional orbital-free density functional theory (OF-DFT). We demonstrate the feasibility of the proposed approach by applying it to densities across the inner crust of neutron stars. Methods: As a first application of OF-DFT for nuclear pasta, we employ the second-order extended Thomas-Fermi (ETF) expansion of Skyrme energy density functional (EDF) to construct an EDF that depends only on neutron and proton number densities. Based on the variational principle, we derive Euler-Lagrange equations to determine optimal neutron and proton density distributions and solve them self-consistently. In this work, we call this approach the self-consistent ETF (SC-ETF) method. Results: We perform three-dimensional SC-ETF calculations with various box sizes. We successfully obtain various pasta structures, depending on given average nucleon number densities, consistent with earlier studies. Moreover, we find other exotic structures, such as bending and/or connected rods, slabs with a hole, etc., underlining the advantage of the self-consistent formalism. Conclusions: We demonstrate that the SC-ETF method proposed in this study, which can be regarded as a realization of OF-DFT, is a promising tool that can efficiently describe complex pasta structures without empirical assumptions on geometric shapes.
Related
- Impact of geomagnetic storm intensity on thermospheric density and orbital decay of LEO satellites: a G1–G5 per-event analysis
- Nuclear equation-of-state at high density and multi-messenger astronomy: contribution of heavy-ion collisions
- Nuclear equation-of-state effects on the two-dimensional post-outburst thermal evolution of magnetized neutron-star crusts
- Three-dimensional spatial modelling of nuclear de-excitation lines
- Bayesian Constraints on Bosonic Dark Matter in Neutron Stars within a Covariant Density-Functional Equation-of-State Ensemble
BibTeX
@misc{nakamura2026threedimensional,
title = {{Three-dimensional orbital-free density functional theory description of nuclear pasta in the inner crust of neutron stars}},
author = {Yo Nakamura and Kazuyuki Sekizawa},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2605.28783},
eprint = {2605.28783},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2605.28783},
}
RIS
TY - UNPB
TI - Three-dimensional orbital-free density functional theory description of nuclear pasta in the inner crust of neutron stars
AU - Yo Nakamura
AU - Kazuyuki Sekizawa
PY - 2026
DA - 2026/09/18
PB - arXiv
DO - 10.48550/arxiv.2605.28783
UR - https://arxiv.org/abs/2605.28783
AB - Background: In the bottom layer of the inner crust of neutron stars, various crystalline structures are expected to emerge that are collectively called ``nuclear pasta.'' It is desirable to know properties of nuclear pasta in a wide variety of conditions for astrophysical applications. However, three-dimensional fully-microscopic calculations require huge computational effort that makes it still challenging to carry out systematic calculations. Purpose: In this paper, we propose an efficient method to calculate various nuclear pasta configurations in a non-empirical manner, based on three-dimensional orbital-free density functional theory (OF-DFT). We demonstrate the feasibility of the proposed approach by applying it to densities across the inner crust of neutron stars. Methods: As a first application of OF-DFT for nuclear pasta, we employ the second-order extended Thomas-Fermi (ETF) expansion of Skyrme energy density functional (EDF) to construct an EDF that depends only on neutron and proton number densities. Based on the variational principle, we derive Euler-Lagrange equations to determine optimal neutron and proton density distributions and solve them self-consistently. In this work, we call this approach the self-consistent ETF (SC-ETF) method. Results: We perform three-dimensional SC-ETF calculations with various box sizes. We successfully obtain various pasta structures, depending on given average nucleon number densities, consistent with earlier studies. Moreover, we find other exotic structures, such as bending and/or connected rods, slabs with a hole, etc., underlining the advantage of the self-consistent formalism. Conclusions: We demonstrate that the SC-ETF method proposed in this study, which can be regarded as a realization of OF-DFT, is a promising tool that can efficiently describe complex pasta structures without empirical assumptions on geometric shapes.
ER -