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.

The use of Slater-type spinor orbitals in algebraic solution of two-center Dirac equation

Relativistic exponential-type spinor orbitals and their use in many-electron Dirac equation solution

Advantages of Slater-type spinor orbitals in the Dirac–Hartree–Fock method. Results for hydrogen-like atoms with super-critical nuclear charge

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

B-Spline Solution of the Two-Center Dirac Equation in the Electronic Continuum for Relativistic Molecular Photoionization

Analytical solutions of the Dirac equation using the Tridiagonal Representation Approach: General study, limitations, and possible applications

Solution of Schrödinger Equation for Simple Diatomic Molecules Using One‐Parameter 1s Slater‐Type Orbitals Wave Function

Direct numerical solution of the two-particle Lippmann-Schwinger equation in coordinate space using the multi-variable Nystrom method

Towards the explicit computation of Bohm velocities associated to N-electron wave-functions with arbitrary spin-orientations

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

Self-Spin-Controlled Rotation of Spatial States of a Dirac Electron in a Cylindrical Potential via Spin-Orbit Interaction

Solving the Time-Dependent Klein-Gordon and Square-Root Klein-Gordon Equations with Krylov-subspace Methods

Methods for accurate calculations of multi-center integrals of squared Coulomb potentials for lower bounds to energy levels of molecular systems

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

An algorithm for exact analytical solutions for tilted anisotropic Dirac materials

Numerical solution of the radial Dirac equation in pseudopotential construction

B-spline finite element studies of the non-linear Schrödinger equation

On the use of complex GTOs for the evaluation of radial integrals involving oscillating functions

Exact solution of two dimensional Dunkl harmonic oscillator in Non-Commutative phase-space

Solution of the Dirac–Coulomb equation using the Rayleigh–Ritz method: Results for He-like atoms