Projected Commutator DIIS Method for Accelerating Hybrid Functional Electronic Structure Calculations.

The commutator direct inversion of the iterative subspace (commutator DIIS or C-DIIS) method developed by Pulay is an efficient and the most widely used scheme in quantum chemistry to accelerate the convergence of self-consistent field (SCF) iterations in Hartree-Fock theory and Kohn-Sham density functional theory. The C-DIIS method requires the explicit storage of the density matrix, the Fock matrix, and the commutator matrix. Hence, the method can only be used for systems with a relatively small basis set, such as the Gaussian basis set. We develop a new method that enables the C-DIIS method to be efficiently employed in electronic structure calculations with a large basis set such as planewaves for the first time. The key ingredient is the projection of both the density matrix and the commutator matrix to an auxiliary matrix called the gauge-fixing matrix. The resulting projected commutator-DIIS method (PC-DIIS) only operates on matrices of the same dimension as that consists of Kohn-Sham orbitals. The cost of the method is comparable to that of standard charge mixing schemes used in large basis set calculations. The PC-DIIS method is gauge-invariant, which guarantees that its performance is invariant with respect to any unitary transformation of the Kohn-Sham orbitals. We demonstrate that the PC-DIIS method can be viewed as an extension of an iterative eigensolver for nonlinear problems. We use the PC-DIIS method for accelerating Kohn-Sham density functional theory calculations with hybrid exchange-correlation functionals, and demonstrate its superior performance compared to the commonly used nested two-level SCF iteration procedure. Furthermore, we demonstrate that in the context of ab initio molecular dynamics (MD) simulation with hybrid functionals one can extrapolate the gauge-fixing matrix to achieve the goal of extrapolating the entire density matrix implicitly along the MD trajectory. Numerical results indicate that the new method significantly reduces the number of SCF iterations per MD step, compared to the commonly used strategy of extrapolating the electron density.

[1]  Matthias Krack,et al.  Efficient and accurate Car-Parrinello-like approach to Born-Oppenheimer molecular dynamics. , 2007, Physical review letters.

[2]  K. Burke,et al.  Rationale for mixing exact exchange with density functional approximations , 1996 .

[3]  G. Scuseria,et al.  Climbing the density functional ladder: nonempirical meta-generalized gradient approximation designed for molecules and solids. , 2003, Physical review letters.

[4]  Donald G. M. Anderson.  Iterative Procedures for Nonlinear Integral Equations , 1965, JACM.

[5]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[6]  G. Kerker.  Efficient iteration scheme for self-consistent pseudopotential calculations , 1981 .

[7]  Gustavo E. Scuseria,et al.  Erratum: “Hybrid functionals based on a screened Coulomb potential” [J. Chem. Phys. 118, 8207 (2003)] , 2006 .

[8]  G. Scuseria,et al.  On the equivalence of LIST and DIIS methods for convergence acceleration. , 2014, The Journal of chemical physics.

[9]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[10]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[11]  Chao Yang,et al.  DGDFT: A massively parallel method for large scale density functional theory calculations. , 2015, The Journal of chemical physics.

[12]  Lexing Ying,et al.  Compressed Representation of Kohn-Sham Orbitals via Selected Columns of the Density Matrix. , 2014, Journal of chemical theory and computation.

[13]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[14]  Lin Lin,et al.  Adaptively Compressed Exchange Operator. , 2016, Journal of chemical theory and computation.

[15]  Gustavo E. Scuseria,et al.  Converging self-consistent field equations in quantum chemistry – recent achievements and remaining challenges , 2007 .

[16]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[17]  Lexing Ying,et al.  Compression of the electron repulsion integral tensor in tensor hypercontraction format with cubic scaling cost , 2014, J. Comput. Phys..

[18]  Anders Nilsson,et al.  Benchmark oxygen-oxygen pair-distribution function of ambient water from x-ray diffraction measurements with a wide Q-range. , 2013, The Journal of chemical physics.

[19]  Chao Yang,et al.  Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations , 2012, J. Comput. Phys..

[20]  J. Zaanen,et al.  Density-functional theory and strong interactions: Orbital ordering in Mott-Hubbard insulators. , 1995, Physical review. B, Condensed matter.

[21]  Chao Yang,et al.  Chebyshev polynomial filtered subspace iteration in the discontinuous Galerkin method for large-scale electronic structure calculations. , 2016, The Journal of chemical physics.

[22]  Chao Yang,et al.  Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene. , 2017, Journal of chemical theory and computation.

[23]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[24]  Joost VandeVondele,et al.  Auxiliary Density Matrix Methods for Hartree-Fock Exchange Calculations. , 2010, Journal of chemical theory and computation.

[25]  Yanli Wang,et al.  Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials , 2009 .

[26]  E Weinan,et al.  Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation , 2011, J. Comput. Phys..

[27]  Arias,et al.  Ab initio molecular dynamics: Analytically continued energy functionals and insights into iterative solutions. , 1992, Physical review letters.

[28]  Stefan Grimme,et al.  Semiempirical GGA‐type density functional constructed with a long‐range dispersion correction , 2006, J. Comput. Chem..

[29]  Chao Yang,et al.  Edge reconstruction in armchair phosphorene nanoribbons revealed by discontinuous Galerkin density functional theory. , 2015, Physical chemistry chemical physics : PCCP.

[30]  S. Goedecker,et al.  Relativistic separable dual-space Gaussian pseudopotentials from H to Rn , 1998, cond-mat/9803286.

[31]  G. Scuseria,et al.  A black-box self-consistent field convergence algorithm: One step closer , 2002 .

[32]  Trygve Helgaker,et al.  The augmented Roothaan-Hall method for optimizing Hartree-Fock and Kohn-Sham density matrices. , 2008, The Journal of chemical physics.

[33]  Juan C. Meza,et al.  A Trust Region Direct Constrained Minimization Algorithm for the Kohn-Sham Equation , 2007, SIAM J. Sci. Comput..

[34]  Adrienn Ruzsinszky,et al.  Strongly Constrained and Appropriately Normed Semilocal Density Functional. , 2015, Physical review letters.

[35]  François Gygi,et al.  Performance and accuracy of recursive subspace bisection for hybrid DFT calculations in inhomogeneous systems. , 2015, Journal of chemical theory and computation.

[36]  Reinhold Schneider,et al.  An analysis for the DIIS acceleration method used in quantum chemistry calculations , 2011 .

[37]  Andrew V. Knyazev,et al.  Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method , 2001, SIAM J. Sci. Comput..

[38]  G. Scuseria,et al.  Hybrid functionals based on a screened Coulomb potential , 2003 .

[39]  Teter,et al.  Separable dual-space Gaussian pseudopotentials. , 1996, Physical review. B, Condensed matter.

[40]  Jií Kolafa,et al.  Time‐reversible always stable predictor–corrector method for molecular dynamics of polarizable molecules , 2004, J. Comput. Chem..

[41]  Gustavo E Scuseria,et al.  Comparison of self-consistent field convergence acceleration techniques. , 2012, The Journal of chemical physics.

[42]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[43]  D. Bowler,et al.  O(N) methods in electronic structure calculations. , 2011, Reports on progress in physics. Physical Society.

[44]  Adrienn Ruzsinszky,et al.  Semilocal density functional obeying a strongly tightened bound for exchange , 2015, Proceedings of the National Academy of Sciences.

[45]  B. Alder,et al.  THE GROUND STATE OF THE ELECTRON GAS BY A STOCHASTIC METHOD , 2010 .

[46]  Chao Yang,et al.  Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon. , 2017, Journal of chemical theory and computation.

[47]  Roberto Car,et al.  The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water. , 2014, The Journal of chemical physics.

[48]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[49]  Chao Yang,et al.  A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix , 2014, J. Comput. Phys..

[50]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[51]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[52]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[53]  Roberto Car,et al.  Nuclear quantum effects in water. , 2008, Physical review letters.

[54]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .

[55]  Anders M N Niklasson,et al.  Extended Born-Oppenheimer molecular dynamics. , 2008, Physical review letters.