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.

Properties of Molecular Electron Density Functionals

Kinetic-energy functional of the electron density

A new Green's function formalism for kinetic energy density functional for atomic and molecular system: Emergence of $N-$dependence using model potentials

Density functional. Theory and application to atoms and molecules

Generalization of the Kohn-Sham system that can represent arbitary one electron density matrices

The role of electron localization in density functionals

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

A well-tempered density functional theory of electrons in molecules

Multicomponent Density-Functional Theory for Electrons and Nuclei

Kinetic-Energy Density-Functional Theory on a Lattice

Quantifying errors of the electron-proton/muon correlation functionals through the Kohn-Sham inversion of a two-component model system

Molecular Binding Energies from Partition Density Functional Theory

Spectral-partitioned Kohn-Sham density functional theory

Kinetic energy density study of some representative semilocal kinetic energy functionals

Ensemble density functional theory of ground and excited energy levels

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

Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals

Energy Density Functionals From the Strong-Coupling Limit Applied to the Anions of the He Isoelectronic Series

Semilocal Kinetic Energy Density Functionals on Atoms and Diatoms

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

Defining the Contributions of Permanent Electrostatics, Pauli Repulsion, and Dispersion in Density Functional Theory Calculations of Intermolecular Interaction Energies.