Abstract:We present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number $Z = 1$ to $83$.

What problem does this paper attempt to address?

This paper attempts to address the problem of improving the accuracy and efficiency of atomic structure calculations in Density Functional Theory (DFT). Specifically, the authors propose a scheme based on spectral methods for atomic structure calculations in pseudopotential Kohn-Sham DFT. By introducing an exponential transformation of the radial coordinate and using global polynomial interpolation on a Chebyshev grid, this method can efficiently and accurately solve electronic structure problems. ### Main Issues and Solutions 1. **Computational Efficiency and Accuracy**: - **Problem**: Traditional atomic structure calculation methods can achieve relatively accurate results, but the computational cost is high when dealing with large-scale systems. - **Solution**: The authors propose a spectral method that maps the radial coordinate to a finite interval through exponential transformation and performs polynomial interpolation on a Chebyshev grid. This method can achieve high precision with fewer grid points, significantly improving computational efficiency. 2. **Handling Different Types of Exchange-Correlation Functionals**: - **Problem**: Existing atomic structure calculation methods usually can only handle local or semi-local exchange-correlation functionals and cannot effectively deal with hybrid functionals. - **Solution**: The proposed spectral method is not only applicable to local and semi-local exchange-correlation functionals but also capable of handling hybrid functionals, thus enhancing the method's generality and applicability. ### Specific Implementation - **Radial Coordinate Transformation**: The radial coordinate $ r $ is mapped to a finite interval $ y $ through exponential transformation, making the grid points denser near the atomic nucleus and gradually sparser away from the nucleus. - **Chebyshev Grid**: In the transformed coordinate system, a Chebyshev grid is used for discretization, employing the Chebyshev differential matrix to approximate derivative operators and using the Clenshaw-Curtis quadrature method for integration. - **Boundary Condition Handling**: By appropriately modifying the Chebyshev differential matrix, Dirichlet boundary conditions are ensured. ### Validation and Results - **Accuracy Validation**: The proposed method's accuracy and convergence under different exchange-correlation functionals were validated by comparison with reference codes ONCVPSP and SPARC. - **Efficiency Evaluation**: Results show that for optimized van der Waals pseudopotentials (ONCV), only about 200 grid points are needed to achieve an accuracy of 1 microhartree. In summary, this paper proposes an efficient spectral method that can significantly reduce the computational cost of atomic structure calculations while maintaining high accuracy, applicable to various exchange-correlation functionals.

Spectral scheme for atomic structure calculations in density functional theory

A spectral scheme for Kohn-Sham density functional theory of clusters

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

Spectral-partitioned Kohn-Sham density functional theory

Spectral approximation scheme for a hybrid, spin-density Kohn–Sham density-functional theory in an external (nonuniform) magnetic field and a collinear exchange-correlation energy

Time-Dependent Density Functional Theory Excited-State Forces with the Gaussian and Augmented Plane Wave Method

Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory

Adaptive Local Basis Set for Kohn-Sham Density Functional Theory in a Discontinuous Galerkin Framework I: Total Energy Calculation

Quasiparticle spectra from a non-empirical optimally-tuned range-separated hybrid density functional

Functional Renormalization Group and Kohn-Sham Scheme in Density Functional Theory

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations

Efficient periodic resolution-of-the-identity Hartree–Fock exchange method with k -point sampling and Gaussian basis sets

Plane-Wave-Based Stochastic-Deterministic Density Functional Theory for Extended Systems

Pseudospectral implementations of long‐range corrected density functional theory

Density Functional Prediction of Quasiparticle, Excitation, and Resonance Energies of Molecules with a Global Scaling Correction Approach

Pseudopotential density functional treament of atoms and molecules in cartesian coordinate grid

A Mixed Basis Density Functional Approach for One-Dimensional Systems with B-splines

Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations

Convergence of density functional iterative procedures with a Newton-Raphson algorithm