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.

Calculation of molecular g-tensors by sampling spin orientations of generalised Hartree-Fock states

Calculation of DFT-GIAO NMR shifts with the inclusion of spin-orbit coupling

Improvements in DFT Calculations of Spin–Spin Coupling Constants

Relativistic magnetic interactions from non-orthogonal basis sets

Four-Component Relativistic Density-Functional Theory Calculations of Nuclear Spin-Rotation Constants: Relativistic Effects in p-Block Hydrides

On the calculation of second-order magnetic properties using subsystem approaches in the relativistic framework

Generalized Kohn-Sham approach for the electronic band structure of spin-orbit coupled materials

Capturing strong correlation in molecules with phaseless auxiliary-field quantum Monte Carlo using generalized Hartree Fock trial wavefunctions

Efficient Four-Component Dirac-Coulomb-Gaunt Hartree--Fock in Pauli Spinor Representation

Efficient Four-Component Dirac-Coulomb-Gaunt Hartree-Fock in the Pauli Spinor Representation.

Spin-pure Stochastic-CASSCF via GUGA-FCIQMC applied to Iron Sulfur Clusters

Efficient Treatment of Relativistic Effects with Periodic Density Functional Methods: Energies, Gradients, and Stress Tensors

Tensor hypercontraction for fully self-consistent imaginary-time GF2 and GWSOX methods: theory, implementation, and role of the Green's function second-order exchange for intermolecular interactions

Projective spin adaptation for the exact diagonalization of isotropic spin clusters

Systematic determination of coupling constants in spin clusters from broken-symmetry mean-field solutions

One-step treatment of spin-orbit coupling and electron correlation in large active spaces

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

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

Spin-flip, Tensor Equation-of-motion Configuration Interaction with a Density-Functional Correction: A Spin-Complete Method for Exploring Excited-State Potential Energy Surfaces.

Perturbative spin-orbit couplings for the simulation of extended framework materials