A hybrid approach to Fermi operator expansion

Abstract

In a recent paper we have suggested that the finite temperature density matrix can be computed efficiently by a combination of polynomial expansion and iterative inversion techniques. We present here significant improvements over this scheme. The original complex-valued formalism is turned into a purely real one. In addition, we use Chebyshev polynomials expansion and fast summation techniques. This drastically reduces the scaling of the algorithm with the width of the Hamiltonian spectrum, which is now of the order of the cubic root of such parameter. This makes our method very competitive for applications to ab-initio simulations, when high energy resolution is required.

preprint of ICCMSE08 proceedings

Cited by in corpus (9)

• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Ab-Initio Molecular Dynamics
• Recent achievements in ab initio modelling of liquid water
• Pole-based approximation of Fermi-Dirac function
• Self-consistent field theory based molecular dynamics with linear system-size scaling
• On the optimal calculation of the pair correlation function for an orthorombic system
• Disordered Crystals from First Principles II: Transport Coefficients
• 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