Expansion algorithm for the density matrix

Abstract

A purification algorithm for expanding the single-particle density matrix in terms of the Hamiltonian operator is proposed. The scheme works with a predefined occupation and requires less than half the number of matrix-matrix multiplications compared to existing methods at low (<10%) and high (>90%) occupancy. The expansion can be used with a fixed chemical potential in which case it is an asymmetric generalization of and a substantial improvement over grand canonical McWeeny purification. It is shown that the computational complexity, measured as number of matrix multiplications, essentially is independent of system size even for metallic materials with a vanishing band gap.

5 pages, 4 figures, to appear in Phys. Rev. B

Cited by in corpus (43)

• O(N) methods in electronic structure calculations
• Locality in Quantum and Markov Dynamics on Lattices and Networks
• Calculations on millions of atoms with DFT: Linear scaling shows its potential
• An assessment of initial guesses for self-consistent field calculations. Superposition of Atomic Potentials: simple yet efficient
• Density Matrix Perturbation Theory
• ELSI -- An Open Infrastructure for Electronic Structure Solvers
• Ab initio linear scaling response theory: Electric polarizability by perturbed projection
• Generalized Extended Lagrangian Born-Oppenheimer Molecular Dynamics
• Gradient-based stochastic estimation of the density matrix
• Implicit Purification for Temperature-Dependent Density Matrices
• Linear scaling computation of the Fock matrix. VIII. Periodic boundaries for exact exchange at the $\Gamma$-point
• Graph-based linear scaling electronic structure theory
• Efficient Computation of Sparse Matrix Functions for Large-Scale Electronic Structure Calculations: The CheSS Library
• Non-monotonic recursive polynomial expansions for linear scaling calculation of the density matrix
• Fast method for quantum mechanical molecular dynamics
• Spin unrestricted linear scaling electronic structure theory and its application to magnetic carbon doped BN nanotubes
• Generalized canonical purification for density matrix minimization
• Higher Order Response in $\mathcal{O}(N)$ by Perturbed Projection
• Non-Orthogonal Density Matrix Perturbation Theory
• Extended Lagrangian free energy molecular dynamics
• Linear scaling computation of the Fock matrix. IX. Parallel computation of the Coulomb matrix
• Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory
• Mixed Precision Fermi-Operator Expansion on Tensor Cores From a Machine Learning Perspective
• Linear scaling computation of the Fock matrix VII. Periodic Density Functional Theory at the $\Gamma$-point
• Parameterless stopping criteria for recursive density matrix expansions
• Linear scaling calculation of band edge states and doped semiconductors
• Using Graph Partitioning for Scalable Distributed Quantum Molecular Dynamics
• Linear scaling calculation of maximally-localized Wannier functions with atomic basis set
• Large-scale DFT methods for calculations of materials with complex structures
• Efficient Parallel Linear Scaling Method to get the Response Density Matrix in All-Electron Real-Space Density-Functional Perturbation Theory
• Graph-based Quantum Response Theory and Shadow Born-Oppenheimer Molecular Dynamics
• Localized density matrix minimization and linear scaling algorithms
• Notes on density matrix perturbation theory
• On-the-fly computation of frontal orbitals in density matrix expansions
• Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures
• Interior Eigensolver for Sparse Hermitian Definite Matrices Based on Zolotarev's Functions
• Approximate multiplication of nearly sparse matrices with decay in a fully recursive distributed task-based parallel framework
• Efficient computation of the density matrix with error control on distributed computer systems
• Sparse approximate matrix-matrix multiplication for density matrix purification with error control
• Multiple eigenvectors around the homo-lumo gap as a cheap by-product in linear scaling electronic structure calculations
• Accelerated density matrix expansions for Born-Oppenheimer molecular dynamics
• Performance Optimizations of Recursive Electronic Structure Solvers targeting Multi-Core Architectures (LA-UR-20-26665)
• Solvers for $\mathcal{O}(N)$ Electronic Structure in the Strong Scaling Limit