Neural-Network-Based Variational Method in Nuclear Density Functional Theory: Application to the Kohn--Sham method
Abstract
We extend the neural-network-based variational method for nuclear density functional theory to the Kohn--Sham scheme, representing the complex spinor components of the single-particle orbitals by multi-layer perceptrons. We show that neural-network optimization of a given energy density functional, combined with an orthonormalization post-processing step, is mathematically equivalent to the variational condition projected onto the tangent space of the wave-function manifold spanned by the network parameters, and that the training optimizes not only the expansion coefficients but also the basis functions themselves. We assess the method from three points of view. In the first place, we examine how the results depend on the number of units, the number of layers, and the arithmetic precision, and find that quantitative accuracy requires both a sufficient width and a sufficient depth, while single-precision arithmetic is sufficient to represent the nuclear density distribution. In the second place, the binding energies and charge radii of several closed-shell nuclei agree with conventional Skyrme--Hartree--Fock results, and the quadrupole deformations of open-shell nuclei are consistent with reference calculations that include pairing and with experiment. In the third place, we confirm that a neural-network single-particle basis can represent the three-dimensional configurations of the fundamental pasta phases: spheres, rods, and slabs. The framework offers a new perspective on computational nuclear theory, well suited to the forthcoming generation of GPU- and AI-oriented high-throughput supercomputers.
Related
- Microscopic Statistical Calculation of Nuclear Level Density Based on Relativistic Density Functional Theory
- Skyrme-quark-meson-coupling energy density functional predictions to nuclear ground state properties
- Nuclear landscape based on point-coupling density functional with localized exchange terms
- Revisiting nuclear chirality in 128Cs with relativistic configuration-interaction density functional theory
- Configuration-interaction time-dependent density functional theory for nuclear dynamics
BibTeX
@misc{yoshimura2026neuralnetworkbased,
title = {{Neural-Network-Based Variational Method in Nuclear Density Functional Theory: Application to the Kohn--Sham method}},
author = {Kenta Yoshimura and Kazuyuki Sekizawa},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.00836},
eprint = {2609.00836},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.00836},
}
RIS
TY - UNPB
TI - Neural-Network-Based Variational Method in Nuclear Density Functional Theory: Application to the Kohn--Sham method
AU - Kenta Yoshimura
AU - Kazuyuki Sekizawa
PY - 2026
DA - 2026/09/01
PB - arXiv
DO - 10.48550/arxiv.2609.00836
UR - https://arxiv.org/abs/2609.00836
AB - We extend the neural-network-based variational method for nuclear density functional theory to the Kohn--Sham scheme, representing the complex spinor components of the single-particle orbitals by multi-layer perceptrons. We show that neural-network optimization of a given energy density functional, combined with an orthonormalization post-processing step, is mathematically equivalent to the variational condition projected onto the tangent space of the wave-function manifold spanned by the network parameters, and that the training optimizes not only the expansion coefficients but also the basis functions themselves. We assess the method from three points of view. In the first place, we examine how the results depend on the number of units, the number of layers, and the arithmetic precision, and find that quantitative accuracy requires both a sufficient width and a sufficient depth, while single-precision arithmetic is sufficient to represent the nuclear density distribution. In the second place, the binding energies and charge radii of several closed-shell nuclei agree with conventional Skyrme--Hartree--Fock results, and the quadrupole deformations of open-shell nuclei are consistent with reference calculations that include pairing and with experiment. In the third place, we confirm that a neural-network single-particle basis can represent the three-dimensional configurations of the fundamental pasta phases: spheres, rods, and slabs. The framework offers a new perspective on computational nuclear theory, well suited to the forthcoming generation of GPU- and AI-oriented high-throughput supercomputers.
ER -