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. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures 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 coefficientmore »
- Authors:
-
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
- Physics Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
- Publication Date:
- Research Org.:
- Georgia Institute of Technology, Atlanta, GA (United States); Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC); USDOE National Nuclear Security Administration (NNSA)
- OSTI Identifier:
- 1867400
- Alternate Identifier(s):
- OSTI ID: 1847017; OSTI ID: 1866178
- Report Number(s):
- LLNL-JRNL-830272
Journal ID: ISSN 0021-9606; TRN: US2306170
- Grant/Contract Number:
- SC0019410; AC52-07NA27344
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Chemical Physics
- Additional Journal Information:
- Journal Volume: 156; Journal Issue: 9; Journal ID: ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; 75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY
Citation Formats
Xu, Qimen, Jing, Xin, Zhang, Boqin, Pask, John E., and Suryanarayana, Phanish. Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature. United States: N. p., 2022.
Web. doi:10.1063/5.0082523.
Xu, Qimen, Jing, Xin, Zhang, Boqin, Pask, John E., & Suryanarayana, Phanish. Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature. United States. https://doi.org/10.1063/5.0082523
Xu, Qimen, Jing, Xin, Zhang, Boqin, Pask, John E., and Suryanarayana, Phanish. Wed .
"Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature". United States. https://doi.org/10.1063/5.0082523. https://www.osti.gov/servlets/purl/1867400.
@article{osti_1867400,
title = {Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature},
author = {Xu, Qimen and Jing, Xin and Zhang, Boqin and Pask, John E. and Suryanarayana, Phanish},
abstractNote = {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. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures 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.},
doi = {10.1063/5.0082523},
journal = {Journal of Chemical Physics},
number = 9,
volume = 156,
place = {United States},
year = {Wed Mar 02 00:00:00 EST 2022},
month = {Wed Mar 02 00:00:00 EST 2022}
}
Works referenced in this record:
SQDFT: Spectral Quadrature method for large-scale parallel Kohn–Sham calculations at high temperature
journal, March 2018
- Suryanarayana, Phanish; Pratapa, Phanisri P.; Sharma, Abhiraj
- Computer Physics Communications, Vol. 224
Frontiers and Challenges in Warm Dense Matter
book, January 2014
- Graziani, Frank; Desjarlais, Michael P.; Redmer, Ronald
- Lecture Notes in Computational Science and Engineering, Vol. 96
On preconditioning the self-consistent field iteration in real-space Density Functional Theory
journal, January 2020
- Kumar, Shashikant; Xu, Qimen; Suryanarayana, Phanish
- Chemical Physics Letters, Vol. 739
SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Isolated clusters
journal, March 2017
- Ghosh, Swarnava; Suryanarayana, Phanish
- Computer Physics Communications, Vol. 212
A mesh-free convex approximation scheme for Kohn–Sham density functional theory
journal, June 2011
- Suryanarayana, Phanish; Bhattacharya, Kaushik; Ortiz, Michael
- Journal of Computational Physics, Vol. 230, Issue 13
\mathcal{O}(N) methods in electronic structure calculations
journal, February 2012
- Bowler, D. R.; Miyazaki, T.
- Reports on Progress in Physics, Vol. 75, Issue 3
Extended application of Kohn-Sham first-principles molecular dynamics method with plane wave approximation at high energy—From cold materials to hot dense plasmas
journal, April 2016
- Zhang, Shen; Wang, Hongwei; Kang, Wei
- Physics of Plasmas, Vol. 23, Issue 4
The Formation of Warm Dense Matter: Experimental Evidence for Electronic Bond Hardening in Gold
journal, February 2009
- Ernstorfer, R.; Harb, M.; Hebeisen, C. T.
- Science, Vol. 323, Issue 5917
Accurate and Scalable Algorithm for First-Principles Molecular-Dynamics Computations on Large Parallel Computers
journal, January 2014
- Osei-Kuffuor, Daniel; Fattebert, Jean-Luc
- Physical Review Letters, Vol. 112, Issue 4
Alternating Anderson–Richardson method: An efficient alternative to preconditioned Krylov methods for large, sparse linear systems
journal, January 2019
- Suryanarayana, Phanish; Pratapa, Phanisri P.; Pask, John E.
- Computer Physics Communications, Vol. 234
Stochastic density functional theory at finite temperatures
journal, March 2018
- Cytter, Yael; Rabani, Eran; Neuhauser, Daniel
- Physical Review B, Vol. 97, Issue 11
Anderson acceleration of the Jacobi iterative method: An efficient alternative to Krylov methods for large, sparse linear systems
journal, February 2016
- Pratapa, Phanisri P.; Suryanarayana, Phanish; Pask, John E.
- Journal of Computational Physics, Vol. 306
On spectral quadrature for linear-scaling Density Functional Theory
journal, October 2013
- Suryanarayana, Phanish
- Chemical Physics Letters, Vol. 584
Introducing ONETEP : Linear-scaling density functional simulations on parallel computers
journal, February 2005
- Skylaris, Chris-Kriton; Haynes, Peter D.; Mostofi, Arash A.
- The Journal of Chemical Physics, Vol. 122, Issue 8
Equation of state of boron nitride combining computation, modeling, and experiment
journal, April 2019
- Zhang, Shuai; Lazicki, Amy; Militzer, Burkhard
- Physical Review B, Vol. 99, Issue 16
Static and dynamic conductivity of warm dense matter within a density-functional approach: Application to aluminum and gold
journal, March 2006
- Dharma-wardana, M. W. C.
- Physical Review E, Vol. 73, Issue 3
Daubechies wavelets for linear scaling density functional theory
journal, May 2014
- Mohr, Stephan; Ratcliff, Laura E.; Boulanger, Paul
- The Journal of Chemical Physics, Vol. 140, Issue 20
Discrete discontinuous basis projection method for large-scale electronic structure calculations
journal, September 2018
- Xu, Qimen; Suryanarayana, Phanish; Pask, John E.
- The Journal of Chemical Physics, Vol. 149, Issue 9
Self-Averaging Stochastic Kohn-Sham Density-Functional Theory
journal, September 2013
- Baer, Roi; Neuhauser, Daniel; Rabani, Eran
- Physical Review Letters, Vol. 111, Issue 10
Development of a Multiphase Beryllium Equation of State and Physics-based Variations
journal, February 2021
- Wu, Christine J.; Myint, Philip C.; Pask, John E.
- The Journal of Physical Chemistry A, Vol. 125, Issue 7
ScaLAPACK Users' Guide
book, January 1997
- Blackford, L. S.; Choi, J.; Cleary, A.
- Society for Industrial and Applied Mathematics
SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Extended systems
journal, July 2017
- Ghosh, Swarnava; Suryanarayana, Phanish
- Computer Physics Communications, Vol. 216
Augmented Lagrangian formulation of orbital-free density functional theory
journal, October 2014
- Suryanarayana, Phanish; Phanish, Deepa
- Journal of Computational Physics, Vol. 275
Non-periodic finite-element formulation of Kohn–Sham density functional theory
journal, February 2010
- Suryanarayana, Phanish; Gavini, Vikram; Blesgen, Thomas
- Journal of the Mechanics and Physics of Solids, Vol. 58, Issue 2
Perspective: Methods for large-scale density functional calculations on metallic systems
journal, December 2016
- Aarons, Jolyon; Sarwar, Misbah; Thompsett, David
- The Journal of Chemical Physics, Vol. 145, Issue 22
Self-consistent-field calculations using Chebyshev-filtered subspace iteration
journal, November 2006
- Zhou, Yunkai; Saad, Yousef; Tiago, Murilo L.
- Journal of Computational Physics, Vol. 219, Issue 1
The SIESTA method for ab initio order- N materials simulation
journal, March 2002
- Soler, José M.; Artacho, Emilio; Gale, Julian D.
- Journal of Physics: Condensed Matter, Vol. 14, Issue 11
Order-N first-principles calculations with the conquest code
journal, July 2007
- Gillan, M. J.; Bowler, D. R.; Torralba, A. S.
- Computer Physics Communications, Vol. 177, Issue 1-2
Extended first-principles molecular dynamics model for high temperature simulations in the Abinit code: Application to warm dense aluminum
journal, February 2022
- Blanchet, A.; Clérouin, J.; Torrent, M.
- Computer Physics Communications, Vol. 271
On the calculation of the stress tensor in real-space Kohn-Sham density functional theory
journal, November 2018
- Sharma, Abhiraj; Suryanarayana, Phanish
- The Journal of Chemical Physics, Vol. 149, Issue 19
Chebyshev-filtered subspace iteration method free of sparse diagonalization for solving the Kohn–Sham equation
journal, October 2014
- Zhou, Yunkai; Chelikowsky, James R.; Saad, Yousef
- Journal of Computational Physics, Vol. 274
Pulay forces from localized orbitals optimized in situ using a psinc basis set
journal, June 2012
- Ruiz-Serrano, Álvaro; Hine, Nicholas D. M.; Skylaris, Chris-Kriton
- The Journal of Chemical Physics, Vol. 136, Issue 23
Large-scale molecular dynamics simulations of dense plasmas: The Cimarron Project
journal, March 2012
- Graziani, Frank R.; Batista, Victor S.; Benedict, Lorin X.
- High Energy Density Physics, Vol. 8, Issue 1
Real-space formulation of the stress tensor for O(N) density functional theory: Application to high temperature calculations
journal, July 2020
- Sharma, Abhiraj; Hamel, Sebastien; Bethkenhagen, Mandy
- The Journal of Chemical Physics, Vol. 153, Issue 3
Periodic Pulay method for robust and efficient convergence acceleration of self-consistent field iterations
journal, March 2016
- Banerjee, Amartya S.; Suryanarayana, Phanish; Pask, John E.
- Chemical Physics Letters, Vol. 647
The ELPA library: scalable parallel eigenvalue solutions for electronic structure theory and computational science
journal, May 2014
- Marek, A.; Blum, V.; Johanni, R.
- Journal of Physics: Condensed Matter, Vol. 26, Issue 21
Spectral Quadrature method for accurate electronic structure calculations of metals and insulators
journal, March 2016
- Pratapa, Phanisri P.; Suryanarayana, Phanish; Pask, John E.
- Computer Physics Communications, Vol. 200
Algorithms and novel applications based on the isokinetic ensemble. I. Biophysical and path integral molecular dynamics
journal, January 2003
- Minary, Peter; Martyna, Glenn J.; Tuckerman, Mark E.
- The Journal of Chemical Physics, Vol. 118, Issue 6
Restarted Pulay mixing for efficient and robust acceleration of fixed-point iterations
journal, August 2015
- Pratapa, Phanisri P.; Suryanarayana, Phanish
- Chemical Physics Letters, Vol. 635
Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma
journal, July 2020
- White, A. J.; Collins, L. A.
- Physical Review Letters, Vol. 125, Issue 5
Higher-order finite-difference formulation of periodic Orbital-free Density Functional Theory
journal, February 2016
- Ghosh, Swarnava; Suryanarayana, Phanish
- Journal of Computational Physics, Vol. 307
SPARC: Simulation Package for Ab-initio Real-space Calculations
journal, July 2021
- Xu, Qimen; Sharma, Abhiraj; Comer, Benjamin
- SoftwareX, Vol. 15