Real-space formulation of the stress tensor for $\mathcal{O}(N)$ density functional theory: application to high temperature calculations

Abstract

We present an accurate and efficient real-space formulation of the Hellmann-Feynman stress tensor for $\mathcal{O}(N)$ Kohn-Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi-Dirac distribution becomes smoother and density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to $\mathcal{O}(N)$ methods. We then describe its evaluation within the $\mathcal{O}(N)$ infinite-cell Clenshaw-Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic as well as insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results, and show convergence with respect to mesh size to established planewave results. We employ the new formulation to compute the viscosity of hydrogen at a million kelvin from Kohn-Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

arXiv comment13 pages, 7 figures

References in corpus (27)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• Optimized norm-conserving Vanderbilt pseudopotentials
• Nearsightedness of Electronic Matter
• O(N) methods in electronic structure calculations
• The first principles calculation of transport coefficients
• The Decay Properties of the Finite Temperature Density Matrix in Metals
• Daubechies Wavelets for Linear Scaling Density Functional Theory
• SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Extended systems
• SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Isolated clusters
• Periodic Pulay method for robust and efficient convergence acceleration of self-consistent field iterations
• The static and dynamic conductivity of warm dense Aluminum and Gold calculated within a density functional approach
• SQDFT: Spectral Quadrature method for large-scale parallel $\mathcal{O}(N)$ Kohn-Sham calculations at high temperature
• Coarse-graining Kohn-Sham Density Functional Theory
• Augmented Lagrangian formulation of Orbital-Free Density Functional Theory
• Ionic and electronic transport properties in dense plasmas by orbital-free density functional theory
• Equation of state of warm-dense boron nitride combining computation, modeling, and experiment
• Implicit Purification for Temperature-Dependent Density Matrices
• Linear scaling DFT calculations for large Tungsten systems using an optimized local basis
• Augmented orbital minimization method for linear scaling electronic structure calculations
• Configurational forces in electronic structure calculations using Kohn-Sham density functional theory
• On nearsightedness in metallic systems for $\mathcal{O}(N)$ Density Functional Theory calculations: A case study on Aluminum
• Spatial decay of the single-particle density matrix in tight-binding metals: +AFwAXA- analytic results in two dimensions
• On preconditioning the self-consistent field iteration in real-space Density Functional Theory
• Canonical density matrix perturbation theory
• On the calculation of the stress tensor in real-space Kohn-Sham Density Functional Theory
• On real-space Density Functional Theory for non-orthogonal crystal systems: Kronecker product formulation of the kinetic energy operator

Cited by in corpus (4)

• Analytical representation of the Local Field Correction of the Uniform Electron Gas within the Effective Static Approximation
• Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics
• Real-space density kernel method for Kohn-Sham density functional theory calculations at high temperature
• Virial Expansion of the Electrical Conductivity of Hydrogen Plasmas