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.

Single-Point Grand Potential and Nuclear Derivatives of Finite-Temperature Kohn-Sham Density-Functional Theory

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

An Efficient Implementation of Analytical Nuclear Gradients for Linear-Response Time-Dependent Density Functional Theory in the Plane Wave Basis

On the Analysis of the Discretized Kohn-Sham Density Functional Theory

D4FT: A Deep Learning Approach to Kohn-Sham Density Functional Theory

Gradient Type Optimization Methods for Electronic Structure Calculations

Very-Large-Scale GPU-Accelerated Nuclear Gradient of Time-Dependent Density Functional Theory with Tamm-Dancoff Approximation and Range-Separated Hybrid Functionals

A Proximal Gradient Method for Ensemble Density Functional Theory

Plane-Wave-Based Stochastic-Deterministic Density Functional Theory for Extended Systems

Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces

Computation of the Kohn-Sham Orbital Kinetic Energy Density in the Full-Potential Linearized Augmented Plane-Wave Method

Sparse tensor based nuclear gradients for periodic Hartree-Fock and low-scaling correlated wave function methods in the CP2K software package: A massively parallel and GPU accelerated implementation

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

Time-dependent Density Functional Theory Quantum Transport Simulation in Non-Orthogonal Basis

Time-Dependent Density Functional Theory Excited-State Forces with the Gaussian and Augmented Plane Wave Method

Efficient periodic resolution-of-the-identity Hartree–Fock exchange method with k -point sampling and Gaussian basis sets

Beyond single-reference fixed-node approximation in ab initio Diffusion Monte Carlo using antisymmetrized geminal power applied to systems with hundreds of electrons

Structure Optimization with Stochastic Density Functional Theory

An Asymptotic Green's Function Method for Time-Dependent Schrödinger Equations with Application to Kohn-Sham Equations

DFT+ within the framework of linear combination of numerical atomic orbitals