Abstract: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 the 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.

General technique for analytical derivatives of post-projected Hartree-Fock

Analytic gradient and derivative couplings for the spin-flip extended configuration interaction singles method: Theory, implementation, and application to proton transfer

Analytical Gradients for Projection-Based Wavefunction-in-DFT Embedding

Projected Commutator DIIS Method for Accelerating Hybrid Functional Electronic Structure Calculations

Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces

Analytical Hartree-Fock gradients for periodic systems

Analytic Derivative Couplings For Spin-Flip Configuration Interaction Singles And Spin-Flip Time-Dependent Density Functional Theory

Analytical Hartree-Fock gradients with respect to the cell parameter for systems periodic in three dimensions

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

The power of Hellmann-Feynman theorem: Kohn-Sham DFT energy derivatives with respect to the parameters of exchange-correlation functional at linear cost

Efficient simulations of Hartree--Fock equations by an accelerated gradient descent method

Analytic gradients for compressed multistate pair-density functional theory

Interpolative separable density fitting decomposition for accelerating Hartree-Fock exchange calculations within numerical atomic orbitals

Analytical Energy Gradient Evaluation in Relativistic and Nonrelativistic Density Functional Calculations

Analytical Gradient Theory for Spin-Free State-Averaged Second-Order Driven Similarity Renormalization Group Perturbation Theory (SA-DSRG-MRPT2) and Its Applications for Conical Intersection Optimizations

Efficient Treatment of Relativistic Effects with Periodic Density Functional Methods: Energies, Gradients, and Stress Tensors

Ground State Energy Functional with Hartree-Fock Efficiency and Chemical Accuracy.

Analytic Gradients for Density Fitting MP2 Using Natural Auxiliary Functions

Structural Relaxation of Materials with Spin-Orbit Coupling: Analytical Forces in Spin-Current DFT

Sparse tensor based nuclear gradients for periodic Hartree-Fock and low-scaling correlated wave function methods in the CP2K software package: A massively parallel and GPU accelerated implementation