Abstract:We present the first full-potential method that solves the fully relativistic 4-component Dirac-Kohn-Sham equation for materials in the solid state within the framework of atom-centered Gaussian-type orbitals (GTOs). Our GTO-based method treats one-, two-, and three-dimensional periodic systems on an equal footing, and allows for a seamless transition to the methodology commonly used in studies of molecules with heavy elements. The scalar relativistic effects and spin-orbit interaction are handled variationally. The full description of the electron-nuclear potential in the core region of heavy nuclei is straightforward due to the local nature of the GTOs and does not pose any computational difficulties. We show how the time-reversal symmetry and a quaternion algebra-based formalism can be exploited to significantly reduce the increased methodological complexity and computational cost associated with multiple wave-function components coupled by the spin-orbit interaction. We provide a description of how to employ the matrix form of the multipole expansion and an iterative renormalization procedure to evaluate the conditionally convergent lattice sums arising for periodic systems. We investigate the problem of inverse variational collapse that arises if the Dirac operator containing a repulsive periodic potential is expressed in a basis containing diffuse functions, and suggest a possible solution. We demonstrate the validity of the method on silver halide crystals with large relativistic effects, and two-dimensional honeycomb structures (silicene and germanene) exhibiting the spin-orbit-driven quantum spin Hall effect. Our results are well-converged with respect to the basis limit using standard bases developed for molecular calculations, and indicate that the rule of removing basis functions with small exponents should not be applied when transferring the molecular basis to solids.

Computation of the Kohn-Sham Orbital Kinetic Energy Density in the Full-Potential Linearized Augmented Plane-Wave Method

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

Partially deorbitalized meta-GGA

All-electron fully relativistic Kohn-Sham theory for solids based on the Dirac-Coulomb Hamiltonian and Gaussian-type functions

Implementation of the SCAN Exchange-Correlation Functional with Numerical Atomic Orbitals

Implementation of the meta-GGA exchange-correlation functional in numerical atomic orbital basis: With systematic testing on SCAN, rSCAN, and r2SCAN functionals

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

Developing Orbital-Dependent Corrections for the Non-Additive Kinetic Energy in Subsystem Density Functional Theory

Achieving accurate ionization potential of semiconductors by the efficient Kohn-Sham scheme of density functional theory

Efficient iterative method for solving the Dirac–Kohn–Sham density functional theory

Kinetic-Energy Density-Functional Theory on a Lattice

Kohn-Sham Potentials From Electron Densities Using A Matrix Representation Within Finite Atomic Orbital Basis Sets

Self-averaging stochastic Kohn-Sham density functional theory

Effective one-particle energies from generalized Kohn-Sham random phase approximation: A direct approach for computing and analyzing core ionization energies

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

Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities

Importance of the Kinetic Energy Density for Band Gap Calculations in Solids with Density Functional Theory

Accurate Ionization Potential of Semiconductors from Efficient Density Functional Calculations

Effect of orbital-overlap dependence in density functionals

Efficient treatment of local meta-generalized gradient density functionals via auxiliary density expansion: the density fitting (DF) J+X approximation

An Asymptotic Green's Function Method for Time-Dependent Schrödinger Equations with Application to Kohn-Sham Equations