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.

FAST ALGORITHM FOR EXTRACTING THE DIAGONAL OF THE INVERSE MATRIX WITH APPLICATION TO THE ELECTRONIC STRUCTURE ANALYSIS OF METALLIC SYSTEMS

A Fast Parallel Algorithm for Selected Inversion of Structured Sparse Matrices with Application to 2D Electronic Structure Calculations.

A Parallel Block Preconditioner-Based VIE-FFT Algorithm for Modeling the Electromagnetic Response from Nanostructures

A Multilevel Method With Novel Correction Strategy for Parallel Finite-Element Analysis of Electromagnetic Problems

The Modified Solutoin for Large Sparse Symmetric Linear Systems in Electromagnetic Field Analysis

SelInv---An Algorithm for Selected Inversion of a Sparse Symmetric Matrix

Fast Simulation of Array Structures Using T-EPA With Hierarchical LU Decomposition

Super Large-Scale Magnetic Data Inversion

Hierarchical Matrices Method and Its Application in Electromagnetic Integral Equations

On-the-fly exact diagonalization solver for quantum electronic models

Diagonalization without Diagonalization: A Direct Optimization Approach for Solid-State Density Functional Theory

FAST ALGORITHMS FOR DETERMINING POSITIVE DEFINITENESS OF LARGE SCALE MATRICES IN SOLID STRUCTURE ANALYSIS1)

Projected Commutator DIIS Method for Accelerating Hybrid Functional Electronic Structure Calculations

A fast low-rank inversion algorithm of dielectric matrix in GW approximation

Fast and spectrally accurate construction of adaptive diagonal basis sets for electronic structure

Efficient Shift-and-Invert Preconditioning for Multi-GPU Accelerated Density Functional Calculations

A Novel Parallel Direct Solver Using Adaptive Cross Approximation for Analysis of Electromagnetic Radiation Problems With Complex Structures

An efficient first-principle approach for electronic structures calculations of nanomaterials.

Combinatorial and Recurrent Approaches for Efficient Matrix Inversion: Sub-cubic algorithms leveraging Fast Matrix products

Fast Phase Retrieval Algrithm for Rotational Symmetric DOE Designing

Superfast direct inversion of the nonuniform discrete Fourier transform via hierarchically semi-separable least squares