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.

Genuine converging solution of self-consistent field equations for extended many-electron systems

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

A Unified Hamiltonian Solution to Maxwell–Schrödinger Equations for Modeling Electromagnetic Field–particle Interaction

Untangling Strongly and Weakly Interacting Configurations in Many-electron Wave Functions

Convergence criteria for the Hartree iterative solution of the general self-consistent field equations

On the Convergence of the Self-Consistent Field Iteration in Kohn--Sham Density Functional Theory

Solving the Bethe-Salpeter Equation in Real Frequencies at Finite Temperature

Converging self-consistent field equations in quantum chemistry – recent achievements and remaining challenges

Eigen-Functional Method And Eikonal-Type Equations In One-Dimensional Strongly Correlated Electron Systems

Solutions of some Schrodinger equations

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

Converging Perturbative Solutions of the Schroedinger Equation for a Two-Level System with a Hamiltonian Depending Periodically on Time

N ov 2 00 0 Eigen-functional bosonization and Eikonal-type equations in one-dimensional strongly correlated electron systems

Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods

A unified quantum framework for electrons and ions: The self-consistent harmonic approximation on a neural network curved manifold

A combined R-matrix eigenstate basis set and finite-differences propagation method for the time-dependent Schr\"{od}dinger equation: the one-electron case

Self-Consistent Approximations to Non-Equilibrium Many-Body Theory

A General Second Order Complete Active Space Self-Consistent-field Solver for Large-Scale Systems

Dielectric screening in extended systems using the self-consistent Sternheimer equation and localized basis sets

Convergence study of the h-adaptive PUM and the hp-adaptive FEM applied to eigenvalue problems in quantum mechanics

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