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 Formalism and Matrix Elements of a Complete Potential-Harmonic Scheme for Directly Solving the Schrödinger Equation of the Helium Atom

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

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

Analytical Structure Matching and Very Precise Approach to the Coulombic Quantum Three-Body Problem

Analytically exact solution of the Schrodinger equation for neutral helium in the ground state

The energy level structure of the modified Schrodinger equation can be consistent with Lamb shift

Practical phase-space electronic Hamiltonians for ab initio dynamics

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

Solution of the Schrodinger equation for quasi-one-dimensional materials using helical waves

Supersymmetric Expansion Algorithm and complete analytical solution for the Hulthén and anharmonic potentials

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

A Unified Scheme of Central Symmetric Shape-Invariant Potentials

Accurate and Efficient Calculations of Hellmann-Feynman Forces for Quantum Computation.

Molecular structure calculations: a unified quantum mechanical description of electrons and nuclei using explicitly correlated Gaussian functions and the global vector representation

A phase-space view of vibrational energies without the Born-Oppenheimer framework

Analytical Expression for Continuum-continuum Transition Amplitude of Hydrogen-like Atoms with Angular-momentum Dependence

Anharmonic Vibrational States of Double-Well Potentials in the Solid State from DFT Calculations

Perturbation Theory Made Efficient and Effective for Predictions of Ionization Potential and Electron Affinity.

Quantum Hydrodynamic Model of Density Functional Theory

Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules