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.

Solving Dirac Equations on A 3d Lattice with Inverse Hamiltonian and Spectral Methods

Efficient Solution for the Dirac Equation in 3D Lattice Space with the Conjugate Gradient Method

Invariant-Related Unitary Transformation Method and Exact Solutions for the Quantum Dirac Field in A Time-Dependent Spatially Homogeneous Electric Field

Hamiltonian Flow Equations for a Dirac Particle in Large Scalar and Vector Potentials

High-accuracy solution of the Dirac equation for three-dimensional hydrogen atom using a variationally optimized diagonalization method

Solving the Dirac equation with nonlocal potential by Imaginary Time Step method

Solving Schrodinger Equation for Three-Electron Quantum Systems by the Use of The Hyperspherical Function Method

Solving the relativistic Hartree-Bogoliubov equation with the finite-difference method

Gaussian Beam Methods for the Dirac Equation in the Semi-classical Regime

On the computation of lattice sums without translational invariance

Exact Boson-Fermion Duality on a 3D Euclidean Lattice

Fast Multiresolution Methods for Density Functional Theory in Nuclear Physics

A Bloch Decomposition-Based Split-Step Pseudospectral Method For Quantum Dynamics With Periodic Potentials

Solution of the 3D logarithmic Schrödinger equation with a central potential

Exponential Collocation Methods Based on Continuous Finite Element Approximations for Efficiently Solving the Cubic Schrödinger Equation

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

A New Study of the Burton and Miller Method for the Solution of A 3d Helmholtz Problem

On Hamiltonian formulations of the Dirac system

Avoid the Tsunami of the Dirac Sea in the Imaginary Time Step Method

Qubit Lattice Algorithms based on the Schrodinger-Dirac representation of Maxwell Equations and their Extensions

Solution of the Dirac Equation with Special Hulthén Potentials