Linear Scaling Density Matrix Real Time TDDFT: Propagator Unitarity \& Matrix Truncation

Abstract

Real time, density matrix based, time dependent density functional theory proceeds through the propagation of the density matrix, as opposed to the Kohn-Sham orbitals. It is possible to reduce the computational workload by imposing spatial cut-off radii on sparse matrices, and the propagation of the density matrix in this manner provides direct access to the optical response of very large systems, which would be otherwise impractical to obtain using the standard formulations of TDDFT. Following a brief summary of our implementation, along with several benchmark tests illustrating the validity of the method, we present an exploration of the factors affecting the accuracy of the approach. In particular we investigate the effect of basis set size and matrix truncation, the key approximation used in achieving linear scaling, on the propagator unitarity and optical spectra. Finally we illustrate that, with an appropriate density matrix truncation range applied, the computational load scales linearly with the system size and discuss the limitations of the approach.

arXiv commentAccepted for publication in J. Chem. Phys

References in corpus (8)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• O(N) methods in electronic structure calculations
• Turbo charging time-dependent density-functional theory with Lanczos chains
• Calculation of the optical response of C60 and Na8 using time-dependent density functional theory and local orbitals
• Calculations on millions of atoms with DFT: Linear scaling shows its potential
• Linear-scaling time-dependent density-functional theory in the linear response formalism
• Calculating optical absorption spectra for large systems using linear-scaling density-functional theory

Cited by in corpus (6)

• Linear-scaling time-dependent density-functional theory (TDDFT) beyond the Tamm-Dancoff approximation: obtaining efficiency and accuracy with in situ optimised local orbitals
• Quantum simulation of exact electron dynamics can be more efficient than classical mean-field methods
• Quantum mechanics in an evolving Hilbert space
• Cost-Efficient High-Resolution Linear Absorption Spectra Through Extrapolating the Dipole Moment from Real-Time Time-Dependent Electronic-Structure Theory
• Linear and Nonlinear Optical Properties of Molecules from Real-Time Propagation Based on the Bethe-Salpeter Equation
• Absolutely Localized Projection-Based Embedding for Excited States