An efficient and accurate decomposition of the Fermi operator

Abstract

We present a method to compute the Fermi function of the Hamiltonian for a system of independent fermions, based on an exact decomposition of the grand-canonical potential. This scheme does not rely on the localization of the orbitals and is insensitive to ill-conditioned Hamiltonians. It lends itself naturally to linear scaling, as soon as the sparsity of the system's density matrix is exploited. By using a combination of polynomial expansion and Newton-like iterative techniques, an arbitrarily large number of terms can be employed in the expansion, overcoming some of the difficulties encountered in previous papers. Moreover, this hybrid approach allows us to obtain a very favorable scaling of the computational cost with increasing inverse temperature, which makes the method competitive with other Fermi operator expansion techniques. After performing an in-depth theoretical analysis of computational cost and accuracy, we test our approach on the DFT Hamiltonian for the metallic phase of the LiAl alloy.

8 pages, 7 figures

References in corpus (6)

• Separable Dual Space Gaussian Pseudo-potentials
• Relativistic separable dual-space Gaussian Pseudopotentials from H to Rn
• Nearsightedness of Electronic Matter
• An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics
• Linear scaling electronic structure calculations and accurate sampling with noisy forces
• Stochastic linear scaling for metals and non metals

Cited by in corpus (20)

• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• O(N) methods in electronic structure calculations
• A Transferable Machine-Learning Model of the Electron Density
• Ab-Initio Molecular Dynamics
• Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion
• Recent achievements in ab initio modelling of liquid water
• Pole-based approximation of Fermi-Dirac function
• A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
• Gradient-based stochastic estimation of the density matrix
• Self-consistent field theory based molecular dynamics with linear system-size scaling
• Multipole Representation of the Fermi Operator with Application to the Electronic Structure Analysis of Metallic Systems
• Extended Lagrangian free energy molecular dynamics
• On the optimal calculation of the pair correlation function for an orthorombic system
• Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations
• Disordered Crystals from First Principles II: Transport Coefficients
• Large-scale DFT methods for calculations of materials with complex structures
• Liquid methane at extreme temperature and pressure: Implications for models of Uranus and Neptune
• Evaluating functions of positive-definite matrices using colored noise thermostats
• A General Algorithm to Calculate the Inverse Principal $p$-th Root of Symmetric Positive Definite Matrices
• Using Approximate Computing for the Calculation of Inverse Matrix p-th Roots