Imaginary-Time Time-Dependent Density Functional Theory and Its Application for Robust Convergence of Electronic States

Positivity preserving density matrix minimization at finite temperatures via square root

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 reinvigorate interest in density matrix minimization methods.

The Structure of Density-Potential Mapping. Part I: Standard Density-Functional Theory

The Hohenberg-Kohn theorem of density-functional theory (DFT) is broadly considered the conceptual basis for a full characterization of an electronic system in its ground state by just the one-body particle density. Part I of this review aims at clarifying the status of the Hohenberg-Kohn theorem within DFT and Part II at different extensions of the theory that include magnetic fields. We collect evidence that the Hohenberg-Kohn theorem does not so much form the basis of DFT, but is rather the consequence of a more comprehensive mathematical framework. Such results are especially useful when it comes to the construction of generalized DFTs.

All-electron real-time and imaginary-time time-dependent density functional theory within a numeric atom-centered basis function framework

Real-time time-dependent density functional theory (RT-TDDFT) is an attractive tool to model quantum dynamics by real-time propagation without the linear response approximation. Sharing the same technical framework of RT-TDDFT, imaginary-time time-dependent density functional theory (it-TDDFT) is a recently developed robust-convergence ground state method. Presented here are high-precision all-electron RT-TDDFT and it-TDDFT implementations within a numerical atom-centered orbital (NAO) basis function framework in the FHI-aims code. We discuss the theoretical background and technical choices in our implementation. First, RT-TDDFT results are validated against linear-response TDDFT results. Specifically, we analyze the NAO basis sets' convergence for Thiel's test set of small molecules and confirm the importance of the augmentation basis functions for adequate convergence. Adopting a velocity-gauge formalism, we next demonstrate applications for systems with periodic boundary conditions. Taking advantage of the all-electron full-potential implementation, we present applications for core level spectra. For it-TDDFT, we confirm that within the all-electron NAO formalism, it-TDDFT can successfully converge systems that are difficult to converge in the standard self-consistent field method. We finally benchmark our implementation for systems up to ∼500 atoms. The implementation exhibits almost linear weak and strong scaling behavior.

Imaginary-time time-dependent density functional theory for periodic systems

Imaginary-time time-dependent density functional theory (it-TDDFT) has been proposed as an alternative method for obtaining the ground state within density functional theory (DFT) which avoids some of the difficulties with convergence encountered by the self-consistent-field (SCF) iterative method. It-TDDFT was previously applied to clusters of atoms where it was demonstrated to converge in select cases where SCF had difficulty with convergence. In the present work we implement it-TDDFT propagation forperiodic systemsby modifying the Quantum ESPRESSO (QE) package, which uses a plane-wave basis with multiplekpoints, and has the options of non-collinear and DFT + U calculations using ultra-soft or norm-conserving pseudo potentials. We demonstrate that our implementation of it-TDDFT propagation with multiplekpoints is correct for DFT + U non-collinear calculations and for DFT + U calculations with ultra-soft pseudo potentials. Our implementation of it-TDDFT propagation converges to the exact SCF energy (up to the decimal guaranteed by double precision) in all but one case where it converged to a slightly lower value than SCF, suggesting a useful alternative for systems where SCF has difficulty reaching the Kohn-Sham (KS) ground state. In addition, we demonstrate that more rapid convergence can be achieved if we use adaptive-size imaginary-time-steps for different kinetic-energy plane-waves.