Optimized Dirac Woods-Saxon basis for covariant density functional theory

TL;DR

This paper introduces an optimized Dirac Woods-Saxon basis that simplifies covariant density functional calculations by reducing basis size and computational resources, while maintaining convergence accuracy.

Contribution

The paper proposes an optimized Dirac Woods-Saxon basis with a potential close to the nuclear mean field, significantly reducing basis size and computational cost in covariant density functional theory.

Findings

Reduced basis space needed for convergence

Elimination of continuum basis in Dirac sea

Significant decrease in computational resources

Abstract

The Woods-Saxon basis has achieved great success in both nonrelativistic and covariant density functional theories in recent years. Due to its nonanalytical nature, however, applications of the Woods-Saxon basis are numerically complicated and computationally time consuming. In this paper, based on the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc), we check in detail the convergence with respect to the basis space in the Dirac sea. An optimized Dirac Woods-Saxon basis is proposed, whose corresponding potential is close to the nuclear mean field. It is shown that the basis space of the optimized Dirac Woods-Saxon basis required for convergence is substantially reduced compared with the original one. In particular, it does not need to contain the bases from continuum in the Dirac sea. The application of the optimized Woods-Saxon basis would greatly reduce computing resource for large-scale density functional calculations.