Search Menu
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information
Search Scholarly Publications

Title: Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon

Abstract

We present a new efficient way to perform hybrid density functional theory (DFT)-based electronic structure calculations. The new method uses an interpolative separable density fitting (ISDF) procedure to construct a set of numerical auxiliary basis vectors and a compact approximation of the matrix consisting of products of occupied orbitals represented in a large basis set such as the planewave basis. Such an approximation allows us to reduce the number of Poisson solves from O($$N$$ $$^{2}_{ϵ}$$) to O(Nϵ) when we apply the exchange operator to occupied orbitals in an iterative method for solving the Kohn–Sham equations, where Ne is the number of electrons in the system to be studied. We show that the ISDF procedure can be carried out in O($$N$$ $$^{3}_{ϵ}$$) operations, with a much smaller preconstant compared to methods used in existing approaches. When combined with the recently developed adaptively compressed exchange (ACE) operator formalism, which reduces the number of times the exchange operator needs to be updated, the resulting ACE-ISDF method significantly reduces the computational cost associated with the exchange operator by nearly 2 orders of magnitude compared to existing approaches for a large silicon system with 1000 atoms. We demonstrate that the ACE-ISDF method can produce accurate energies and forces for insulating and metallic systems and that it is possible to obtain converged hybrid functional calculation results for a 1000-atom bulk silicon within 10 min on 2000 computational cores. We also show that ACE-ISDF can scale to 8192 computational cores for a 4096-atom bulk silicon system. We use the ACE-ISDF method to geometrically optimize a 1000-atom silicon system with a vacancy defect using the HSE06 functional and computes its electronic structure. Here, we find that that the computed energy gap from the HSE06 functional is much closer to the experimental value compared to that produced by semilocal functionals in the DFT calculations.

Authors:
ORCiD logo [1]; ORCiD logo [2];  [3]
+ Show Author Affiliations
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
  2. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Univ. of California, Berkeley, CA (United States)
  3. Univ. of California, Berkeley, CA (United States)
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE Office of Science (SC), Basic Energy Sciences (BES)
OSTI Identifier:
1480242
Grant/Contract Number:  
DMS-1652330
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Chemical Theory and Computation
Additional Journal Information:
Journal Volume: 13; Journal Issue: 11; Journal ID: ISSN 1549-9618
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; 97 MATHEMATICS AND COMPUTING

Citation Formats

Hu, Wei, Lin, Lin, and Yang, Chao. Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon. United States: N. p., 2017. Web. doi:10.1021/acs.jctc.7b00807.
Copy to clipboard
Hu, Wei, Lin, Lin, & Yang, Chao. Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon. United States. https://doi.org/10.1021/acs.jctc.7b00807
Copy to clipboard
Hu, Wei, Lin, Lin, and Yang, Chao. Fri . "Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon". United States. https://doi.org/10.1021/acs.jctc.7b00807. https://www.osti.gov/servlets/purl/1480242.
Copy to clipboard
@article{osti_1480242,
title = {Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon},
author = {Hu, Wei and Lin, Lin and Yang, Chao},
abstractNote = {We present a new efficient way to perform hybrid density functional theory (DFT)-based electronic structure calculations. The new method uses an interpolative separable density fitting (ISDF) procedure to construct a set of numerical auxiliary basis vectors and a compact approximation of the matrix consisting of products of occupied orbitals represented in a large basis set such as the planewave basis. Such an approximation allows us to reduce the number of Poisson solves from O($N$ $^{2}_{ϵ}$) to O(Nϵ) when we apply the exchange operator to occupied orbitals in an iterative method for solving the Kohn–Sham equations, where Ne is the number of electrons in the system to be studied. We show that the ISDF procedure can be carried out in O($N$ $^{3}_{ϵ}$) operations, with a much smaller preconstant compared to methods used in existing approaches. When combined with the recently developed adaptively compressed exchange (ACE) operator formalism, which reduces the number of times the exchange operator needs to be updated, the resulting ACE-ISDF method significantly reduces the computational cost associated with the exchange operator by nearly 2 orders of magnitude compared to existing approaches for a large silicon system with 1000 atoms. We demonstrate that the ACE-ISDF method can produce accurate energies and forces for insulating and metallic systems and that it is possible to obtain converged hybrid functional calculation results for a 1000-atom bulk silicon within 10 min on 2000 computational cores. We also show that ACE-ISDF can scale to 8192 computational cores for a 4096-atom bulk silicon system. We use the ACE-ISDF method to geometrically optimize a 1000-atom silicon system with a vacancy defect using the HSE06 functional and computes its electronic structure. Here, we find that that the computed energy gap from the HSE06 functional is much closer to the experimental value compared to that produced by semilocal functionals in the DFT calculations.},
doi = {10.1021/acs.jctc.7b00807},
journal = {Journal of Chemical Theory and Computation},
number = 11,
volume = 13,
place = {United States},
year = {Fri Sep 29 00:00:00 EDT 2017},
month = {Fri Sep 29 00:00:00 EDT 2017}
}
Copy to clipboard
Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
https://doi.org/10.1021/acs.jctc.7b00807
Copyright Statement
Other availability
Search WorldCat Search WorldCat to find libraries that may hold this journal
Citation Metrics:
Cited by: 32 works
Citation information provided by
Web of Science
Save / Share:
Export Metadata
Save to My Library
You must Sign In or Create an Account in order to save documents to your library.

Works referenced in this record:

Inhomogeneous Electron Gas
journal, November 1964

Self-Consistent Equations Including Exchange and Correlation Effects
journal, November 1965

Jacob’s ladder of density functional approximations for the exchange-correlation energy
conference, January 2001

  • Perdew, John P.
  • Density functional theory and its application to materials, AIP Conference Proceedings
  • DOI: 10.1063/1.1390175

A new mixing of Hartree–Fock and local density‐functional theories
journal, January 1993

  • Becke, Axel D.
  • The Journal of Chemical Physics, Vol. 98, Issue 2
  • DOI: 10.1063/1.464304

Rationale for mixing exact exchange with density functional approximations
journal, December 1996

  • Perdew, John P.; Ernzerhof, Matthias; Burke, Kieron
  • The Journal of Chemical Physics, Vol. 105, Issue 22, p. 9982-9985
  • DOI: 10.1063/1.472933

Hybrid functionals based on a screened Coulomb potential
journal, May 2003

  • Heyd, Jochen; Scuseria, Gustavo E.; Ernzerhof, Matthias
  • The Journal of Chemical Physics, Vol. 118, Issue 18
  • DOI: 10.1063/1.1564060

Erratum: “Hybrid functionals based on a screened Coulomb potential” [J. Chem. Phys. 118, 8207 (2003)]
journal, June 2006

  • Heyd, Jochen; Scuseria, Gustavo E.; Ernzerhof, Matthias
  • The Journal of Chemical Physics, Vol. 124, Issue 21
  • DOI: 10.1063/1.2204597

Ground State of the Electron Gas by a Stochastic Method
journal, August 1980

Self-interaction correction to density-functional approximations for many-electron systems
journal, May 1981

Separable dual-space Gaussian pseudopotentials
journal, July 1996

Density-functional exchange-energy approximation with correct asymptotic behavior
journal, September 1988

Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density
journal, January 1988

Generalized Gradient Approximation Made Simple
journal, October 1996

  • Perdew, John P.; Burke, Kieron; Ernzerhof, Matthias
  • Physical Review Letters, Vol. 77, Issue 18, p. 3865-3868
  • DOI: 10.1103/PhysRevLett.77.3865

Climbing the Density Functional Ladder: Nonempirical Meta–Generalized Gradient Approximation Designed for Molecules and Solids
journal, September 2003

Strongly Constrained and Appropriately Normed Semilocal Density Functional
journal, July 2015

Semilocal density functional obeying a strongly tightened bound for exchange
journal, January 2015

  • Sun, Jianwei; Perdew, John P.; Ruzsinszky, Adrienn
  • Proceedings of the National Academy of Sciences, Vol. 112, Issue 3
  • DOI: 10.1073/pnas.1423145112

A scalable and accurate algorithm for the computation of Hartree–Fock exchange
journal, May 2010

Parallel implementation of γ-point pseudopotential plane-wave DFT with exact exchange
journal, November 2010

  • Bylaska, Eric J.; Tsemekhman, Kiril; Baden, Scott B.
  • Journal of Computational Chemistry, Vol. 32, Issue 1
  • DOI: 10.1002/jcc.21598

Improved treatment of exact exchange in Quantum ESPRESSO
journal, May 2017

  • Barnes, Taylor A.; Kurth, Thorsten; Carrier, Pierre
  • Computer Physics Communications, Vol. 214
  • DOI: 10.1016/j.cpc.2017.01.008

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

Auxiliary Density Matrix Methods for Hartree−Fock Exchange Calculations
journal, July 2010

  • Guidon, Manuel; Hutter, Jürg; VandeVondele, Joost
  • Journal of Chemical Theory and Computation, Vol. 6, Issue 8
  • DOI: 10.1021/ct1002225

The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water
journal, August 2014

  • DiStasio, Robert A.; Santra, Biswajit; Li, Zhaofeng
  • The Journal of Chemical Physics, Vol. 141, Issue 8
  • DOI: 10.1063/1.4893377

Performance and Accuracy of Recursive Subspace Bisection for Hybrid DFT Calculations in Inhomogeneous Systems
journal, September 2015

  • Dawson, William; Gygi, François
  • Journal of Chemical Theory and Computation, Vol. 11, Issue 10
  • DOI: 10.1021/acs.jctc.5b00826

Compressed Representation of Kohn–Sham Orbitals via Selected Columns of the Density Matrix
journal, March 2015

  • Damle, Anil; Lin, Lin; Ying, Lexing
  • Journal of Chemical Theory and Computation, Vol. 11, Issue 4
  • DOI: 10.1021/ct500985f

Comparison of Three Efficient Approximate Exact-Exchange Algorithms: The Chain-of-Spheres Algorithm, Pair-Atomic Resolution-of-the-Identity Method, and Auxiliary Density Matrix Method
journal, July 2016

  • Rebolini, Elisa; Izsák, Róbert; Reine, Simen Sommerfelt
  • Journal of Chemical Theory and Computation, Vol. 12, Issue 8
  • DOI: 10.1021/acs.jctc.6b00074

Fast, accurate evaluation of exact exchange: The occ-RI-K algorithm
journal, July 2015

  • Manzer, Samuel; Horn, Paul R.; Mardirossian, Narbe
  • The Journal of Chemical Physics, Vol. 143, Issue 2
  • DOI: 10.1063/1.4923369

Efficient Computation of the Hartree–Fock Exchange in Real-Space with Projection Operators
journal, July 2016

  • Boffi, Nicholas M.; Jain, Manish; Natan, Amir
  • Journal of Chemical Theory and Computation, Vol. 12, Issue 8
  • DOI: 10.1021/acs.jctc.6b00376

Adaptively Compressed Exchange Operator
journal, April 2016

Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene
journal, February 2017

  • Hu, Wei; Lin, Lin; Banerjee, Amartya S.
  • Journal of Chemical Theory and Computation, Vol. 13, Issue 3
  • DOI: 10.1021/acs.jctc.6b01184

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials
journal, September 2009

  • Giannozzi, Paolo; Baroni, Stefano; Bonini, Nicola
  • Journal of Physics: Condensed Matter, Vol. 21, Issue 39, Article No. 395502
  • DOI: 10.1088/0953-8984/21/39/395502

Compression of the electron repulsion integral tensor in tensor hypercontraction format with cubic scaling cost
journal, December 2015

Simplifications in the generation and transformation of two-electron integrals in molecular calculations
journal, October 1977

  • Beebe, Nelson H. F.; Linderberg, Jan
  • International Journal of Quantum Chemistry, Vol. 12, Issue 4
  • DOI: 10.1002/qua.560120408

Reduced scaling in electronic structure calculations using Cholesky decompositions
journal, June 2003

  • Koch, Henrik; Sánchez de Merás, Alfredo; Pedersen, Thomas Bondo
  • The Journal of Chemical Physics, Vol. 118, Issue 21
  • DOI: 10.1063/1.1578621

Low-cost evaluation of the exchange Fock matrix from Cholesky and density fitting representations of the electron repulsion integrals
journal, May 2007

  • Aquilante, Francesco; Pedersen, Thomas Bondo; Lindh, Roland
  • The Journal of Chemical Physics, Vol. 126, Issue 19
  • DOI: 10.1063/1.2736701

A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency
journal, August 2002

  • Weigend, Florian
  • Physical Chemistry Chemical Physics, Vol. 4, Issue 18
  • DOI: 10.1039/b204199p

Tensor hypercontraction. II. Least-squares renormalization
journal, December 2012

  • Parrish, Robert M.; Hohenstein, Edward G.; Martínez, Todd J.
  • The Journal of Chemical Physics, Vol. 137, Issue 22
  • DOI: 10.1063/1.4768233

Discrete variable representation in electronic structure theory: Quadrature grids for least-squares tensor hypercontraction
journal, May 2013

  • Parrish, Robert M.; Hohenstein, Edward G.; Martínez, Todd J.
  • The Journal of Chemical Physics, Vol. 138, Issue 19
  • DOI: 10.1063/1.4802773

Singular value decomposition and least squares solutions
journal, April 1970

  • Golub, G. H.; Reinsch, C.
  • Numerische Mathematik, Vol. 14, Issue 5
  • DOI: 10.1007/BF02163027

On the low-rank approximation by the pivoted Cholesky decomposition
journal, April 2012

  • Harbrecht, Helmut; Peters, Michael; Schneider, Reinhold
  • Applied Numerical Mathematics, Vol. 62, Issue 4
  • DOI: 10.1016/j.apnum.2011.10.001

Some Applications of the Rank Revealing QR Factorization
journal, May 1992

  • Chan, Tony F.; Hansen, Per Christian
  • SIAM Journal on Scientific and Statistical Computing, Vol. 13, Issue 3
  • DOI: 10.1137/0913043

Randomized algorithms for the low-rank approximation of matrices
journal, December 2007

  • Liberty, E.; Woolfe, F.; Martinsson, P. -G.
  • Proceedings of the National Academy of Sciences, Vol. 104, Issue 51
  • DOI: 10.1073/pnas.0709640104

Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions
journal, January 2011

  • Halko, N.; Martinsson, P. G.; Tropp, J. A.
  • SIAM Review, Vol. 53, Issue 2
  • DOI: 10.1137/090771806

Parallel solution of partial symmetric eigenvalue problems from electronic structure calculations
journal, December 2011

Adaptive local basis set for Kohn–Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation
journal, February 2012

DGDFT: A massively parallel method for large scale density functional theory calculations
journal, September 2015

  • Hu, Wei; Lin, Lin; Yang, Chao
  • The Journal of Chemical Physics, Vol. 143, Issue 12
  • DOI: 10.1063/1.4931732

Edge reconstruction in armchair phosphorene nanoribbons revealed by discontinuous Galerkin density functional theory
journal, January 2015

  • Hu, Wei; Lin, Lin; Yang, Chao
  • Physical Chemistry Chemical Physics, Vol. 17, Issue 47
  • DOI: 10.1039/C5CP00333D

Chebyshev polynomial filtered subspace iteration in the discontinuous Galerkin method for large-scale electronic structure calculations
journal, October 2016

  • Banerjee, Amartya S.; Lin, Lin; Hu, Wei
  • The Journal of Chemical Physics, Vol. 145, Issue 15
  • DOI: 10.1063/1.4964861

Relativistic separable dual-space Gaussian pseudopotentials from H to Rn
journal, August 1998

Linear scaling computation of the Hartree–Fock exchange matrix
journal, August 1996

  • Schwegler, Eric; Challacombe, Matt
  • The Journal of Chemical Physics, Vol. 105, Issue 7
  • DOI: 10.1063/1.472135

Linear and sublinear scaling formation of Hartree–Fock-type exchange matrices
journal, August 1998

  • Ochsenfeld, Christian; White, Christopher A.; Head-Gordon, Martin
  • The Journal of Chemical Physics, Vol. 109, Issue 5
  • DOI: 10.1063/1.476741

Implementation of screened hybrid density functional for periodic systems with numerical atomic orbitals: Basis function fitting and integral screening
journal, July 2011

  • Shang, Honghui; Li, Zhenyu; Yang, Jinlong
  • The Journal of Chemical Physics, Vol. 135, Issue 3
  • DOI: 10.1063/1.3610379

Implementation of Exact Exchange with Numerical Atomic Orbitals
journal, November 2009

  • Shang, Honghui; Li, Zhenyu; Yang, Jinlong
  • The Journal of Physical Chemistry A, Vol. 114, Issue 2
  • DOI: 10.1021/jp908836z

Defects in silicon
journal, October 1982

Intrinsic defects in silicon
journal, August 2000

First-principles calculations for point defects in solids
journal, March 2014

  • Freysoldt, Christoph; Grabowski, Blazej; Hickel, Tilmann
  • Reviews of Modern Physics, Vol. 86, Issue 1
  • DOI: 10.1103/RevModPhys.86.253

Energy band gaps and lattice parameters evaluated with the Heyd-Scuseria-Ernzerhof screened hybrid functional
journal, November 2005

  • Heyd, Jochen; Peralta, Juan E.; Scuseria, Gustavo E.
  • The Journal of Chemical Physics, Vol. 123, Issue 17
  • DOI: 10.1063/1.2085170

Photoluminescence of size-separated silicon nanocrystals: Confirmation of quantum confinement
journal, June 2002

  • Ledoux, G.; Gong, J.; Huisken, F.
  • Applied Physics Letters, Vol. 80, Issue 25
  • DOI: 10.1063/1.1485302

Two-Point Step Size Gradient Methods
journal, January 1988

  • Barzilai, Jonathan; Borwein, Jonathan M.
  • IMA Journal of Numerical Analysis, Vol. 8, Issue 1
  • DOI: 10.1093/imanum/8.1.141

Improved semiconductor lattice parameters and band gaps from a middle-range screened hybrid exchange functional
journal, March 2012

    Sort by:
    [ × clear filter / sort ]

    Works referencing / citing this record:

    An auxiliary-Field quantum Monte Carlo perspective on the ground state of the dense uniform electron gas: An investigation with Hartree-Fock trial wavefunctions
    journal, August 2019

    • Lee, Joonho; Malone, Fionn D.; Morales, Miguel A.
    • The Journal of Chemical Physics, Vol. 151, Issue 6
    • DOI: 10.1063/1.5109572
      Sort by:
      [ × clear filter / sort ]

      Similar Records in DOE PAGES and OSTI.GOV collections: