Robust determination of the chemical potential in the pole expansion and selected inversion method for solving KohnSham density functional theory
Abstract
Fermi operator expansion (FOE) methods are powerful alternatives to diagonalization type methods for solving KohnSham density functional theory (KSDFT). One example is the pole expansion and selected inversion (PEXSI) method, which approximates the Fermi operator by rational matrix functions and reduces the computational complexity to at most quadratic scaling for solving KSDFT. Unlike diagonalization type methods, the chemical potential often cannot be directly read off from the result of a single step of evaluation of the Fermi operator. Hence multiple evaluations are needed to be sequentially performed to compute the chemical potential to ensure the correct number of electrons within a given tolerance. This hinders the performance of FOE methods in practice. In this paper, we develop an efficient and robust strategy to determine the chemical potential in the context of the PEXSI method. The main idea of the new method is not to find the exact chemical potential at each selfconsistentfield (SCF) iteration but to dynamically and rigorously update the upper and lower bounds for the true chemical potential, so that the chemical potential reaches its convergence along the SCF iteration. Instead of evaluating the Fermi operator for multiple times sequentially, our method uses a twolevel strategy that evaluatesmore »
 Authors:

 University of California, Berkeley, CA (United States). Dept. of Mathematics
 University of California, Berkeley, CA (United States). Dept. of Mathematics; Lawrence Berkeley National Laboratory, Berkeley, CA (United States). Computational Research Division
 Publication Date:
 Research Org.:
 Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center
 Sponsoring Org.:
 USDOE Office of Science (SC)
 OSTI Identifier:
 1497884
 Grant/Contract Number:
 AC0205CH11231
 Resource Type:
 Accepted Manuscript
 Journal Name:
 Journal of Chemical Physics
 Additional Journal Information:
 Journal Volume: 147; Journal Issue: 14; Journal ID: ISSN 00219606
 Publisher:
 American Institute of Physics (AIP)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY
Citation Formats
Jia, Weile, and Lin, Lin. Robust determination of the chemical potential in the pole expansion and selected inversion method for solving KohnSham density functional theory. United States: N. p., 2017.
Web. doi:10.1063/1.5000255.
Jia, Weile, & Lin, Lin. Robust determination of the chemical potential in the pole expansion and selected inversion method for solving KohnSham density functional theory. United States. doi:https://doi.org/10.1063/1.5000255
Jia, Weile, and Lin, Lin. Wed .
"Robust determination of the chemical potential in the pole expansion and selected inversion method for solving KohnSham density functional theory". United States. doi:https://doi.org/10.1063/1.5000255. https://www.osti.gov/servlets/purl/1497884.
@article{osti_1497884,
title = {Robust determination of the chemical potential in the pole expansion and selected inversion method for solving KohnSham density functional theory},
author = {Jia, Weile and Lin, Lin},
abstractNote = {Fermi operator expansion (FOE) methods are powerful alternatives to diagonalization type methods for solving KohnSham density functional theory (KSDFT). One example is the pole expansion and selected inversion (PEXSI) method, which approximates the Fermi operator by rational matrix functions and reduces the computational complexity to at most quadratic scaling for solving KSDFT. Unlike diagonalization type methods, the chemical potential often cannot be directly read off from the result of a single step of evaluation of the Fermi operator. Hence multiple evaluations are needed to be sequentially performed to compute the chemical potential to ensure the correct number of electrons within a given tolerance. This hinders the performance of FOE methods in practice. In this paper, we develop an efficient and robust strategy to determine the chemical potential in the context of the PEXSI method. The main idea of the new method is not to find the exact chemical potential at each selfconsistentfield (SCF) iteration but to dynamically and rigorously update the upper and lower bounds for the true chemical potential, so that the chemical potential reaches its convergence along the SCF iteration. Instead of evaluating the Fermi operator for multiple times sequentially, our method uses a twolevel strategy that evaluates the Fermi operators in parallel. In the regime of full parallelization, the wall clock time of each SCF iteration is always close to the time for one single evaluation of the Fermi operator, even when the initial guess is far away from the converged solution. We demonstrate the effectiveness of the new method using examples with metallic and insulating characters, as well as results from ab initio molecular dynamics.},
doi = {10.1063/1.5000255},
journal = {Journal of Chemical Physics},
number = 14,
volume = 147,
place = {United States},
year = {2017},
month = {10}
}
Web of Science
Figures / Tables:
Works referenced in this record:
XIX. A demonstration of the theorem that every homogeneous quadratic polynomial is reducible by real orthogonal substitutions to the form of a sum of positive and negative squares
journal, August 1852
 Sylvester, J. J.
 The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, Vol. 4, Issue 23
SIESTAPEXSI: massively parallel method for efficient and accurate ab initio materials simulation without matrix diagonalization
journal, July 2014
 Lin, Lin; García, Alberto; Huhs, Georg
 Journal of Physics: Condensed Matter, Vol. 26, Issue 30
Accelerating atomic orbitalbased electronic structure calculation via pole expansion and selected inversion
journal, June 2013
 Lin, Lin; Chen, Mohan; Yang, Chao
 Journal of Physics: Condensed Matter, Vol. 25, Issue 29
SelInvAn Algorithm for Selected Inversion of a Sparse Symmetric Matrix
journal, February 2011
 Lin, Lin; Yang, Chao; Meza, Juan C.
 ACM Transactions on Mathematical Software, Vol. 37, Issue 4
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
Linearscaling densityfunctional theory with tens of thousands of atoms: Expanding the scope and scale of calculations with ONETEP
journal, July 2009
 Hine, N. D. M.; Haynes, P. D.; Mostofi, A. A.
 Computer Physics Communications, Vol. 180, Issue 7
SuperLU_DIST: A scalable distributedmemory sparse direct solver for unsymmetric linear systems
journal, June 2003
 Li, Xiaoye S.; Demmel, James W.
 ACM Transactions on Mathematical Software, Vol. 29, Issue 2
PoleBased approximation of the FermiDirac function
journal, August 2009
 Lin, Lin; Lu, Jianfeng; Ying, Lexing
 Chinese Annals of Mathematics, Series B, Vol. 30, Issue 6
Monotonic cubic spline interpolation
conference, January 1999
 Wolber,
 Proceedings Computer Graphics International CGI99
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
Recent progress in linear scaling ab initio electronic structure techniques
journal, March 2002
 Bowler, D. R.; Miyazaki, T.; Gillan, M. J.
 Journal of Physics: Condensed Matter, Vol. 14, Issue 11
Minimax rational approximation of the FermiDirac distribution
journal, October 2016
 Moussa, Jonathan E.
 The Journal of Chemical Physics, Vol. 145, Issue 16
Integral representation of the Fermi distribution and its applications in electronicstructure calculations
journal, December 1993
 Goedecker, S.
 Physical Review B, Vol. 48, Issue 23
\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
Enhancing Scalability and Load Balancing of Parallel Selected Inversion via TreeBased Asynchronous Communication
conference, May 2016
 Jacquelin, Mathias; Lin, Lin; Wichmann, Nathan
 2016 IEEE International Parallel and Distributed Processing Symposium (IPDPS)
Ab initio molecular simulations with numeric atomcentered orbitals
journal, November 2009
 Blum, Volker; Gehrke, Ralf; Hanke, Felix
 Computer Physics Communications, Vol. 180, Issue 11
Nearsightedness of electronic matter
journal, August 2005
 Prodan, E.; Kohn, W.
 Proceedings of the National Academy of Sciences, Vol. 102, Issue 33
Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach
journal, April 2005
 VandeVondele, Joost; Krack, Matthias; Mohamed, Fawzi
 Computer Physics Communications, Vol. 167, Issue 2
EdgeModified Phosphorene Nanoflake Heterojunctions as Highly Efficient Solar Cells
journal, February 2016
 Hu, Wei; Lin, Lin; Yang, Chao
 Nano Letters, Vol. 16, Issue 3
Electronic structure and aromaticity of largescale hexagonal graphene nanoflakes
journal, December 2014
 Hu, Wei; Lin, Lin; Yang, Chao
 The Journal of Chemical Physics, Vol. 141, Issue 21
Adaptive local basis set for Kohn–Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation
journal, February 2012
 Lin, Lin; Lu, Jianfeng; Ying, Lexing
 Journal of Computational Physics, Vol. 231, Issue 4
PSelInv—A Distributed Memory Parallel Algorithm for Selected Inversion: The Symmetric Case
journal, December 2016
 Jacquelin, Mathias; Lin, Lin; Yang, Chao
 ACM Transactions on Mathematical Software, Vol. 43, 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
Figures / Tables found in this record: