Real-space density kernel method for Kohn-Sham density functional theory calculations at high temperature

Abstract

Kohn-Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. We present a density matrix based method for Kohn-Sham calculations at high temperature that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann-Feynman forces, and Hellmann-Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw-Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials; and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116,045 K from Kohn-Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

arXiv comment13 pages, 2 figures, 1 table

References in corpus (24)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• Optimized norm-conserving Vanderbilt pseudopotentials
• O(N) methods in electronic structure calculations
• The first principles calculation of transport coefficients
• Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration
• Daubechies Wavelets for Linear Scaling Density Functional Theory
• SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Extended systems
• Self-averaging stochastic Kohn-Sham density functional theory
• 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
• Extended First-Principles Molecular Dynamics Method From Cold Materials to Hot Dense Plasmas
• Stochastic Density Functional Theory at Finite Temperatures
• The static and dynamic conductivity of warm dense Aluminum and Gold calculated within a density functional approach
• Fast and Universal Kohn Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma
• SQDFT: Spectral Quadrature method for large-scale parallel $\mathcal{O}(N)$ Kohn-Sham calculations at high temperature
• 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
• Augmented orbital minimization method for linear scaling electronic structure calculations
• On preconditioning the self-consistent field iteration in real-space Density Functional Theory
• On the calculation of the stress tensor in real-space Kohn-Sham Density Functional Theory
• Real-space formulation of the stress tensor for $\mathcal{O}(N)$ density functional theory: application to high temperature calculations
• Bypassing the malfunction junction in warm dense matter simulations