Abstract:The ab initio description of the spectral interior of the absorption spectrum poses both a theoretical and computational challenge for modern electronic structure theory. Due to the often spectrally dense character of this domain in the quantum propagator's eigenspectrum for medium-to-large sized systems, traditional approaches based on the partial diagonalization of the propagator often encounter oscillatory and stagnating convergence. Electronic structure methods which solve the molecular response problem through the solution of spectrally shifted linear systems, such as the complex polarization propagator, offer an alternative approach which is agnostic to the underlying spectral density or domain location. This generality comes at a seemingly high computational cost associated with solving a large linear system for each spectral shift in some discretization of the spectral domain of interest. In this work, we present a novel, adaptive solution to this high computational overhead based on model order reduction techniques via interpolation. Model order reduction reduces the computational complexity of mathematical models and is ubiquitous in the simulation of dynamical systems and control theory. The efficiency and effectiveness of the proposed algorithm in the ab initio prediction of X-ray absorption spectra is demonstrated using a test set of challenging water clusters which are spectrally dense in the neighborhood of the oxygen K-edge. On the basis of a single, user defined tolerance we automatically determine the order of the reduced models and approximate the absorption spectrum up to the given tolerance. We also illustrate that, for the systems studied, the automatically determined model order increases logarithmically with the problem dimension, compared to a linear increase of the number of eigenvalues within the energy window. Furthermore, we observed that the computational cost of the proposed algorithm only scales quadratically with respect to the problem dimension.

A Structure Preserving Lanczos Algorithm for Computing the Optical Absorption Spectrum

Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT.

Accelerating Optical Absorption Spectra and Exciton Energy Computation via Interpolative Separable Density Fitting

Structure Preserving Parallel Algorithms for Solving the Bethe-Salpeter Eigenvalue Problem

Fast Optical Absorption Spectra Calculations for Periodic Solid State Systems

Computing linear optical spectra in the presence of nonadiabatic effects on Graphics Processing Units using molecular dynamics and tensor-network approaches

Structure-preserving ΓQR and Γ-Lanczos Algorithms for Bethe-Salpeter Eigenvalue Problems.

A new and efficient approach to time-dependent density-functional perturbation theory for optical spectroscopy

Computing the density of states for optical spectra by low-rank and QTT tensor approximation

Linear Scaling Approach for Optical Excitations Using Maximally Localized Wannier Functions

Fast Algorithms for Estimating the Absorption Spectrum within the Linear Response Time-dependent Density Functional Theory

Linear Scaling Time-Dependent Density-Functional Tight-Binding Method for Absorption Spectra of Large Systems

Energy-Specific Bethe-Salpeter Equation Implementation for Efficient Optical Spectrum Calculations

Block Lanczos for lattice QCD spectroscopy and matrix elements

Lanczos, the transfer matrix, and the signal-to-noise problem

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid TDDFT Calculations

Lawrence Berkeley National Laboratory Recent Work Title Properties of definite Bethe – Salpeter eigenvalue problems

Using symmetry-adapted optimized sum-of-products basis functions to calculate vibrational spectra

Model Order Reduction Algorithm for Estimating the Absorption Spectrum

An efficient quantum algorithm for generation of ab initio n-th order susceptibilities for non-linear spectroscopies

A continuous approach for computing the pseudospectra of linear operators