DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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 » 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.« less

Authors:
 [1]; ORCiD logo [1];  [1];  [2]; ORCiD logo [1]
  1. College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
  2. 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 O ( N ) Kohn–Sham calculations at high temperature
journal, March 2018

  • Suryanarayana, Phanish; Pratapa, Phanisri P.; Sharma, Abhiraj
  • Computer Physics Communications, Vol. 224
  • DOI: 10.1016/j.cpc.2017.12.003

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
  • DOI: 10.1007/978-3-319-04912-0

On preconditioning the self-consistent field iteration in real-space Density Functional Theory
journal, January 2020


SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Isolated clusters
journal, March 2017


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
  • DOI: 10.1016/j.jcp.2011.03.018

\mathcal{O}(N) methods in electronic structure calculations
journal, February 2012


The Formation of Warm Dense Matter: Experimental Evidence for Electronic Bond Hardening in Gold
journal, February 2009


Accurate and Scalable O ( N ) Algorithm for First-Principles Molecular-Dynamics Computations on Large Parallel Computers
journal, January 2014


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
  • DOI: 10.1016/j.cpc.2018.07.007

Stochastic density functional theory at finite temperatures
journal, March 2018


Electronic Structure
book, January 2004


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
  • DOI: 10.1016/j.jcp.2015.11.018

On spectral quadrature for linear-scaling Density Functional Theory
journal, October 2013


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
  • DOI: 10.1063/1.1839852

Equation of state of boron nitride combining computation, modeling, and experiment
journal, April 2019


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
  • DOI: 10.1063/1.4871876

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
  • DOI: 10.1063/1.5037794

Self-Averaging Stochastic Kohn-Sham Density-Functional Theory
journal, September 2013


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
  • DOI: 10.1021/acs.jpca.0c09809

ScaLAPACK Users' Guide
book, January 1997

  • Blackford, L. S.; Choi, J.; Cleary, A.
  • Society for Industrial and Applied Mathematics
  • DOI: 10.1137/1.9780898719642

SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Extended systems
journal, July 2017


Augmented Lagrangian formulation of orbital-free density functional theory
journal, October 2014


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
  • DOI: 10.1016/j.jmps.2009.10.002

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
  • DOI: 10.1063/1.4972007

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
  • DOI: 10.1016/j.jcp.2006.03.017

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
  • DOI: 10.1088/0953-8984/14/11/302

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
  • DOI: 10.1016/j.cpc.2007.02.075

Extended first-principles molecular dynamics model for high temperature simulations in the Abinit code: Application to warm dense aluminum
journal, February 2022


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
  • DOI: 10.1063/1.5057355

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
  • DOI: 10.1016/j.jcp.2014.06.056

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
  • DOI: 10.1063/1.4728026

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
  • DOI: 10.1016/j.hedp.2011.06.010

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
  • DOI: 10.1063/5.0016783

Periodic Pulay method for robust and efficient convergence acceleration of self-consistent field iterations
journal, March 2016


The ELPA library: scalable parallel eigenvalue solutions for electronic structure theory and computational science
journal, May 2014


Spectral Quadrature method for accurate O ( N ) electronic structure calculations of metals and insulators
journal, March 2016

  • Pratapa, Phanisri P.; Suryanarayana, Phanish; Pask, John E.
  • Computer Physics Communications, Vol. 200
  • DOI: 10.1016/j.cpc.2015.11.005

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
  • DOI: 10.1063/1.1534582

Restarted Pulay mixing for efficient and robust acceleration of fixed-point iterations
journal, August 2015


Ab Initio Molecular Dynamics
book, January 2009


Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma
journal, July 2020


Higher-order finite-difference formulation of periodic Orbital-free Density Functional Theory
journal, February 2016


SPARC: Simulation Package for Ab-initio Real-space Calculations
journal, July 2021