PaperPanorama
FAQ

Nuclear Theory·nucl-th

Friday·January 3, 2020

11 papers—2 primary·9 cross-listed

All11Primary2Cross-lists9
  1. 01

    Non-relativistic expansion of Dirac equation with spherical scalar and vector potentials by reconstituted Foldy-Wouthuysen transformation

    Yixin Guo · Haozhao Liang

    Inspired by the reconstituted similarity renormalization group method, the reconstituted Foldy-Wouthuysen (FW) transformation is proposed. Applied to the Dirac equation in the covariant density functional theory, the reconstituted FW transformation shows a fast convergence of the spectrum of the single-particle energy. The single-particle densities and the single-particle scalar densities obtained by this new method are also investigated. In particular, the relativistic corrections to the densities from the picture-change error between the Schrödinger and Dirac pictures are discussed in detail. Taking these relativistic corrections into account, both the single-particle densities and the single-particle scalar densities are almost identical to their exact values.

    nucl-thmath-phmath.MPquant-phPRC(2020)·10 citations
  2. 02

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

    Yu Kashiwaba · Takashi Nakatsukasa

    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 Ba. The robustness of the shape coexistence feature in 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 and the temperature of 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.

    nucl-thPRC(2020)·20 citations

Affiliations

first authorsco-authorsvia INSPIRE