Coordinate-space solver for finite-temperature Hartree-Fock-Bogoliubov calculation using the shifted Krylov method

Abstract

In order to study structure of proto-neutron stars and those in subsequent cooling stages, it is of great interest to calculate inhomogeneous hot and cold nuclear matter in a variety of phases. The finite-temperature Hartree-Fock-Bogoliubov (FT-HFB) theory is a primary choice for this purpose, however, its numerical calculation for superfluid (superconducting) many-fermion systems in three dimensions requires enormous computational costs. To study a variety of phases in the crust of hot and cold neutron stars, we propose an efficient method to perform the FT-HFB calculation with the three-dimensional (3D) coordinate-space representation. Recently, an efficient method based on the contour integral of Green's function with the shifted conjugate-orthogonal conjugate-gradient method has been proposed [Phys. Rev. C 95, 044302 (2017)]. We extend the method to the finite temperature, using the shifted conjugate-orthogonal conjugate-residual method. We benchmark the 3D coordinate-space solver of the FT-HFB calculation for hot isolated nuclei and fcc phase in the inner crust of neutron stars at finite temperature. The computational performance of the present method is demonstrated. Different critical temperatures of the quadrupole and the octupole deformations are confirmed for ${}^{146}$Ba. The robustness of the shape coexistence feature in ${}^{184}$Hg is examined. For the neutron-star crust, the deformed neutron-rich Se nuclei embedded in the sea of superfluid low-density neutrons appear in the fcc phase at the nucleon density of 0.045 fm${}^{-3}$ and the temperature of $k_{B}T=200$ keV. The efficiency of the developed solver is demonstrated for nuclei and inhomogeneous nuclear matter at finite temperature. It may provide a standard tool for nuclear physics, especially for the structure of the hot and cold neutron-star matters.

arXiv comment13 pages, 9 figures

References in corpus (33)

• Pairing in nuclear systems: from neutron stars to finite nuclei
• Nuclear Quantum Many-Body Dynamics: From Collective Vibrations to Heavy-Ion Collisions
• Time-dependent density-functional description of nuclear dynamics
• Coordinate-space solution of the Skyrme-Hartree-Fock-Bogolyubov equations within spherical symmetry. The program HFBRAD (v1.0)
• Renormalization of the Hartree-Fock-Bogoliubov Equations in the Case of a Zero Range Pairing Interaction
• Spatial structure of neutron Cooper pair in low density uniform matter
• Canonical-basis time-dependent Hartree-Fock-Bogoliubov theory and linear-response calculations
• Fission Barriers of Compound Superheavy Nuclei
• Isovector Giant Dipole Resonance from the 3D Time-Dependent Density Functional Theory for Superfluid Nuclei
• Broyden's Method in Nuclear Structure Calculations
• Time-Dependent Density Functional Theory and the Real-Time Dynamics of Fermi Superfluids
• Axially deformed solution of the Skyrme-Hartree-Fock-Bogolyubov equations using the transformed harmonic oscillator basis (III) hfbtho (v3.00): a new version of the program
• Application of the gradient method to Hartree-Fock-Bogoliubov theory
• Vortex pinning and dynamics in the neutron star crust
• Novel Role of Superfluidity in Low-Energy Nuclear Reactions
• Solution of the Skyrme-Hartree-Fock-Bogolyubov equations in the Cartesian deformed harmonic-oscillator basis. (VIII) hfodd (v2.73y): a new version of the program
• Solution of the Skyrme HF+BCS equation on a 3D mesh. II. A new version of the Ev8 code
• Energy Density Functional Approach to Superfluid Nuclei
• Microscopic Structure of High-Spin Vibrational Excitations in Superdeformed 190,192,194Hg
• Systematic Study of Fission Barriers of Excited Superheavy Nuclei
• Systematic investigation of low-lying dipole modes using the canonical-basis time-dependent Hartree-Fock-Bogoliubov theory
• Compilation of Giant Electric Dipole Resonances Built on Excited States
• Octupole deformation in the nuclear chart based on the 3D Skyrme Hartree-Fock plus BCS model
• Density functional approaches to collective phenomena in nuclei: Time-dependent density-functional theory for perturbative and non-perturbative nuclear dynamics
• Systematic of isovector and isoscalar giant quadrupole resonances in normal and superfluid deformed nuclei
• Linear Responses in Time-dependent Hartree-Fock-Bogoliubov Method with Gogny Interaction
• Gauge angle dependence in TDHFB calculations of ${}^{20}$O + ${}^{20}$O head-on collisions with the Gogny interaction
• Relativistic Coulomb excitation within Time Dependent Superfluid Local Density Approximation
• Coordinate-Space Solver for Superfluid Many-Fermion Systems with Shifted Conjugate Orthogonal Conjugate Gradient Method
• Transfer probabilities for the reactions ${}^{14,20}$O+${}^{20}$O in terms of multiple time-dependent Hartree-Fock-Bogoliubov trajectories
• Self-consistent band calculation of slab phase in neutron-star crust
• Density-constraint Time-dependent Hartree-Fock-Bogoliubov method
• Reduced-Shifted Conjugate-Gradient Method for a Green's Function: Efficient Numerical Approach in a Nano-structured Superconductor

Cited by in corpus (16)

• Quantal diffusion approach for multinucleon transfer processes in the ${}^{58,64}$Ni+${}^{208}$Pb reactions: Toward the production of unknown neutron-rich nuclei
• Unified neutron star EOSs and neutron star structures in RMF models
• Nuclear pasta structures and symmetry energy
• Unified nuclear matter EOSs constrained by the in-medium balance in density-dependent covariant density functionals
• Superfluid extension of the self-consistent time-dependent band theory for neutron star matter: Anti-entrainment versus superfluid effects in the slab phase
• Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation
• Microscopic derivation of density functional theory for superfluid systems based on effective action formalism
• Nuclear pasta structures at high temperatures
• Cranked Skyrme-Hartree-Fock-Bogoliubov approach for a mean-field description of nuclear rotations near the drip line
• Properties and microscopic structures of dense stellar matter in RMF models
• Enhanced moments of inertia for rotation in neutron-rich nuclei
• Fermi operator expansion method for nuclei and inhomogeneous matter with nuclear energy density functional
• Self-consistent energy density functional approaches to the crust of neutron stars
• The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions
• Framework for Polarized Superfluid Fermion Systems
• Fermi Operator Expansion for the Hartree-Fock-Bogoliubov Theory