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.

Third-Generation TB-LMTO

Higher Derivative Muffin Tin Orbitals (HDMTO) and Higher Derivative Koringa Khon and Rostoker (HDKKR) methods

Excited States Calculated by Means of the Linear Muffin-Tin Orbital Method

LAPW vs. LMTO full-potential simulations and anharmonic dynamics of KNbO3

Local Second Order Mo̷ller-Plesset Theory with a Single Threshold Using Orthogonal Virtual Orbitals: A Distributed Memory Implementation

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

Generalized Kohn-Sham Approach for the Electronic Band Structure of Spin-Orbit Coupled Materials

Dual Bethe-Salpeter equation for the multi-orbital lattice susceptibility within dynamical mean-field theory

An efficient solution of the multi-term multi-harmonic electron Boltzmann equation for use in global models

A Programmable K · P Hamiltonian Method and Application to Magnetic Topological Insulator MnBi2Te4

Multiplet ligand-field theory using Wannier orbitals

Theoretical Modeling of Thermodynamic and Mechanical Properties of the Pure Components of Ti and Zr Based Alloys Using the Exact Muffin-Tin Orbitals Method

First-principles calculation of the parameters used by atomistic magnetic simulations

The Orbitally Selective Mott Phase in Electron Doped Twisted TMDs: A Possible Realization of the Kondo Lattice Model

Solving 2D and 3D lattice models of correlated fermions -- combining matrix product states with mean field theory

How HF solutions for the TB interacting electrons in 2D predict SCES properties and suggest phenomenology for attaining RTS

First principles residual resistivity using locally self-consistent multiple scattering method

Open-shell Tensor Hypercontraction

Quadratic to linear magnetoresistance tuning in TmB4

A programmable $k\cdot p$ Hamiltonian method and application to magnetic topological insulator MnBi$_2$Te$_4$

Linear and Nonlinear Optical Properties from TDOMP2 Theory