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.

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

Higher orbital integrals, rho numbers and index theory

Efficient simulation of complex Ginzburg--Landau equations using high-order exponential-type methods

Computationally Efficient Molecular Integrals of Solid Harmonic Gaussian Orbitals Using Quantum Entanglement of Angular Momentum

On a class of oscillatory integrals and their application to the time dependent Schrödinger equation

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

Efficient integration of the variational equations of multi-dimensional Hamiltonian systems: Application to the Fermi-Pasta-Ulam lattice

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

Continuous approximation for linear combination coefficients: Exploring a neglected concept

Efficient evaluation of two-center Gaussian integrals in periodic systems

Analytic Evaluation of some 2-, 3- and 4- Electron Atomic Integrals Containing Exponentially Correlated Functions of $r_{ij}$

Influence of gauges in the numerical simulation of the time-dependent Ginzburg-Landau model

On a power series involving classical orthogonal polynomials

Application of complex energy integration to selfconsistent electronic structure calculations

Analytic evaluation of Coulomb integrals for one, two and three-electron distance operators, $R_{C1}^{-n}R_{D1}^{-m}$, $R_{C1}^{-n}r_{12}^{-m}$ and $r_{12}^{-n}r_{13}^{-m}$ with $n, m=0,1,2$

On The Time Splitting Spectral Method For The Complex Ginzburg-Landau Equation In The Large Time And Space Scale Limit

Excitations and spectra from equilibrium real-time Green's functions

One-Range Addition Theorems in Terms of -ETOs for STOs and Coulomb-Yukawa Like Correlated Interaction Potentials of Integer and Noninteger Indices

High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics

The Operator Product Expansion for Radial Lattice Quantization of 3D $ϕ^4$ Theory

Multiplicity of critical orbits to nonlinear, strongly indefinite functionals with sign-changing nonlinear part