Turbo charging time-dependent density-functional theory with Lanczos chains

Abstract

We introduce a new implementation of time-dependent density-functional theory which allows the \emph{entire} spectrum of a molecule or extended system to be computed with a numerical effort comparable to that of a \emph{single} standard ground-state calculation. This method is particularly well suited for large systems and/or large basis sets, such as plane waves or real-space grids. By using a super-operator formulation of linearized time-dependent density-functional theory, we first represent the dynamical polarizability of an interacting-electron system as an off-diagonal matrix element of the resolvent of the Liouvillian super-operator. One-electron operators and density matrices are treated using a representation borrowed from time-independent density-functional perturbation theory, which permits to avoid the calculation of unoccupied Kohn-Sham orbitals. The resolvent of the Liouvillian is evaluated through a newly developed algorithm based on the non-symmetric Lanczos method. Each step of the Lanczos recursion essentially requires twice as many operations as a single step of the iterative diagonalization of the unperturbed Kohn-Sham Hamiltonian. Suitable extrapolation of the Lanczos coefficients allows for a dramatic reduction of the number of Lanczos steps necessary to obtain well converged spectra, bringing such number down to hundreds (or a few thousands, at worst) in typical plane-wave pseudopotential applications. The resulting numerical workload is only a few times larger than that needed by a ground-state Kohn-Sham calculation for a same system. Our method is demonstrated with the calculation of the spectra of benzene, C${}_{60}$ fullerene, and of chlorofyll a.

15 pages, 7 figures, to be pdflatex + bibtex

References in corpus (5)

• Phonons and related properties of extended systems from density-functional perturbation theory
• A new and efficient approach to time-dependent density-functional perturbation theory for optical spectroscopy
• Calculation of the optical response of C60 and Na8 using time-dependent density functional theory and local orbitals
• Time-Dependent Density Functional Theory with Ultrasoft Pseudopotential: Real-Time Electron Propagation across Molecular Junction
• Classical mechanics technique for quantum linear response

Cited by in corpus (56)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Advanced capabilities for materials modelling with Quantum ESPRESSO
• Quantum ESPRESSO toward the exascale
• Many-body perturbation theory calculations using the yambo code
• GW quasi-particle spectra from occupied states only
• turboEELS -- A code for the simulation of the electron energy loss and inelastic X-ray scattering spectra using the Liouville-Lanczos approach to time-dependent density-functional perturbation theory
• turboTDDFT 2.0 - Hybrid functionals and new algorithms within time-dependent density-functional perturbation theory
• Linear-scaling time-dependent density-functional theory in the linear response formalism
• Revealing the Mechanism of Low-Energy Electron Yield Enhancement from Sensitizing Nanoparticles
• Time-Dependent Self Consistent Harmonic Approximation: Anharmonic nuclear quantum dynamics and time correlation functions
• A finite field approach to solving the Bethe Salpeter equation
• Electron energy-loss and inelastic X-ray scattering cross sections from time-dependent density-functional perturbation theory
• Cubic-scaling iterative solution of the Bethe-Salpeter equation for finite systems
• First-principles modeling of plasmons in aluminum under ambient and extreme conditions
• Electron production by sensitizing gold nanoparticles irradiated by fast ions
• Linear-scaling time-dependent density-functional theory (TDDFT) beyond the Tamm-Dancoff approximation: obtaining efficiency and accuracy with in situ optimised local orbitals
• In situ Characterization of Nanoparticles Using Rayleigh Scattering
• Magnon-phonon interactions enhance the gap at the Dirac point in the spin-wave spectra of CrI${}_3$ two-dimensional magnets
• Spin dynamics from time-dependent density functional perturbation theory
• The Wannier Function Software Ecosystem for Materials Simulations
• Linear Scaling Density Matrix Real Time TDDFT: Propagator Unitarity \& Matrix Truncation
• Full spectrum of the Liouvillian of open dissipative quantum systems in the Zeno limit
• A structure preserving Lanczos algorithm for computing the optical absorption spectrum
• Fast iterative solution of the Bethe-Salpeter eigenvalue problem using low-rank and QTT tensor approximation
• Wigner Gaussian dynamics: simulating the anharmonic and quantum ionic motion
• Efficient dielectric matrix calculations using the Lanczos algorithm for fast many-body $G_0{W}_0$ implementations
• Quantum and classical features of the photoionization spectrum of C${}_{60}$
• Solid-State Effects on the Optical Excitation of Push-Pull Molecular J-Aggregates by First-Principles Simulations
• Electron energy loss spectroscopy of bulk gold with ultrasoft pseudopotentials and the Liouville-Lanczos method
• Multimodel Approach to the Optical Properties of Molecular Dyes in Solution
• turboMagnon -- A code for the simulation of spin-wave spectra using the Liouville-Lanczos approach to time-dependent density-functional perturbation theory
• Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory
• Molecular Polarizability of Water from the Local Dielectric Response Theory
• First-principles DFT + GW study of the Te antisite in CdTe
• On the Kohn--Sham density response in a localized basis set
• Computational protocol to evaluate electron-phonon interactions within density matrix perturbation theory
• Photoionization of multishell fullerenes studied by ab initio and model approaches
• Photoabsorption spectrum of the Xe@C${}_{60}$ endohedral fullerene
• First-principles study of the gap in the spin excitation spectrum of the CrI${}_3$ honeycomb ferromagnet
• Applying the Liouville-Lanczos Method of Time-Dependent Density-Functional Theory to Warm Dense Matter
• Optical spectrum of proflavine and its ions
• Fast construction of the Kohn--Sham response function for molecules
• Effect of C${}_{60}$ giant resonance on the photoabsorption of encaged atoms
• GPU-Accelerated Solution of the Bethe-Salpeter Equation for Large and Heterogeneous Systems
• Formation of C${}_{18}$H and C${}_{18}$H${}_2$ molecules by low energy irradiation with atomic and molecular hydrogen
• Magnons from time-dependent density-functional perturbation theory and nonempirical Hubbard functionals
• Locality and Computational Reliability of Linear Response Calculations for Molecular Systems
• Wavelets for Density-Functional Theory and Post-Density-Functional-Theory Calculations
• A study of the density functional methods on the photoabsorption of Bodipy dyes
• Analytical Excited-State Gradients and Derivative Couplings in TDDFT with Minimal Auxiliary Basis Set Approximation and GPU Acceleration
• Formation of Micrometer-Sized Textured Hexagonal Silicon Crystals via Nanoindentation
• Comparative study of magnetic exchange parameters and magnon dispersions in NiO and MnO from first principles
• Variational formulation of dynamical electronic response functions in presence of nonlocal exchange interactions
• Influence of C${}_6$H${}_4$(OH)${}_2$ isomers on water disinfection by photocatalysis: a computational study
• Computing the density of states for optical spectra by low-rank and QTT tensor approximation
• Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange