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.

Numerical solution of the radial Dirac equation in pseudopotential construction

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

Acceptable solutions of the Schrodinger radial equation for a particle in a two-dimensional central potential

Non-variational computation of the eigenstates of Dirac operators with radially symmetric potentials

On the calculation of irregular solutions of the Schrödinger equation for non-spherical potentials

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

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

On numerical solutions of the time-dependent Schrödinger equation

Calculating the many-potential vacuum polarization density of the Dirac equation in the finite-basis approximation

Energy spectra with the Dirac equation of the q-deformed generalized Pöschl-Teller potential via the Feynman approach for 39 K 2 a 3 ∑ u +

Exact Solutions to the Dirac Equation for A Coulomb Potential in D+1 Dimensions

A iterative finite method approach to Kohn-Sham equations of light atoms: Density Functional Theory tutorial for undergraduate students

In Situ Pseudopotentials for Electronic Structure Theory

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

Accurate Quantum States for a 2D-Dipole

Remarks on the Solution of the Position Dependent Mass (PDM) Schrödinger Equation

Calculation of Relativistic Single-Particle States

LHC and VLHC Based ep Colliders: E-Linac vs e-Ring

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

Numerical Solution to Poisson's Equation for Estimating Electrostatic Properties Resulting from an Axially Symmetric Gaussian Charge Density Distribution: Charge in Free Space

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