Positivity Preserving Density Matrix Minimization at Finite Temperatures via Square Root

Abstract

We present a Wave Operator Minimization (WOM) method for calculating the Fermi-Dirac density matrix for electronic structure problems at finite temperature while preserving physicality by construction using the wave operator, i.e., the square root of the density matrix. WOM models cooling a state initially at infinite temperature down to the desired finite temperature. We consider both the grand canonical (constant chemical potential) and canonical (constant number of electrons) ensembles. Additionally, we show that the number of steps required for convergence is independent of the number of atoms in the system. We hope that the discussion and results presented in this article reinvigorates interest in density matrix minimization methods.

10 pages, 5 figures, 1 algorithm, supplementary information 5 pages, 3 figures

References in corpus (11)

• Nearsightedness of Electronic Matter
• SIESTA: recent developments and applications
• Challenges in Large Scale Quantum Mechanical Calculations
• Towards Electronic Structure-Based Ab-Initio Molecular Dynamics Simulations with Hundreds of Millions of Atoms
• Next generation extended Lagrangian first principles molecular dynamics
• Low rank approximation for the numerical simulation of high dimensional Lindblad and Riccati equations
• Beating the House: Fast Simulation of Dissipative Quantum Systems with Ensemble Rank Truncation
• Continuous-time dynamics and error scaling of noisy highly-entangling quantum circuits
• The wave operator representation of quantum and classical dynamics
• Non-stochastic matrix Schrödinger equation for open systems
• Graph-based Quantum Response Theory and Shadow Born-Oppenheimer Molecular Dynamics