Fermi operator expansion method for nuclei and inhomogeneous matter with nuclear energy density functional
arXiv:2211.09448 · doi:10.1103/PhysRevC.107.015802
Abstract
The nuclear energy density functional method at finite temperature is a useful tool for studies of nuclear structure at high excitation, and also for researches of nuclear matter involved in explosive stellar phenomena and neutron stars. However, its unrestricted calculation requires large computational costs for the three-dimensional coordinate-space solvers, especially for the Hamiltonian matrix diagonalization and (or) the Gram-Schmidt orthonormalization of the single-particle wave functions. We test numerical performance of the Fermi operator expansion method, that requires neither the diagonalization nor the Gram-Schmidt orthonormalization, for finite nuclei and inhomogeneous nuclear matter. The method is applied to isolated finite N=Z nuclei and to non-uniform symmetric nuclear matter at finite temperature, which turns out be very effective with the three-dimensional coordinate-space representation, especially at high temperature. The Fermi operator expansion method is a useful tool for studies of various nuclear phases at finite temperature with the energy density functional calculations. The method is suitable for massively parallel computing with distributed memory. Furthermore, when the space size is large, the calculation may benefit from its order-N scaling property.
arXiv comment14 pages, 16 figures
References in corpus (18)
- Relativistic Equation of State of Nuclear Matter for Supernova and Neutron Star
- Physics of Neutron Star Crusts
- The TDHF Code Sky3D
- Disordered nuclear pasta, magnetic field decay, and crust cooling in neutron stars
- Nuclear equation of state for core-collapse supernova simulations with realistic nuclear forces
- Linear response theory in the continuum for deformed nuclei: Green's function vs. time-dependent Hartree-Fock with the absorbing-boundary condition
- The Decay Properties of the Finite Temperature Density Matrix in Metals
- Band structure effects for dripped neutrons in neutron star crust
- Quantum corrections to the ground state energy of inhomogeneous neutron matter
- Astromaterial Science and Nuclear Pasta
- Broyden's Method in Nuclear Structure Calculations
- Nuclear "pasta" structures in low-density nuclear matter and neutron star crust
- Phases of hot nuclear matter at subnuclear densities
- Role of the boundary conditions in the Wigner-Seitz approximation applied to the neutron star inner crust
- Coordinate-Space Solver for Superfluid Many-Fermion Systems with Shifted Conjugate Orthogonal Conjugate Gradient Method
- Coordinate-space solver for finite-temperature Hartree-Fock-Bogoliubov calculation using the shifted Krylov method
- Nuclear pasta structures at high temperatures
- A finite-temperature Hartree-Fock code for shell-model Hamiltonians