A Model Order Reduction Algorithm for Estimating the Absorption Spectrum

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.

[1]  Patrick Norman,et al.  Linear complex polarization propagator in a four-component Kohn-Sham framework. , 2010, The Journal of chemical physics.

[2]  D. M. Bishop,et al.  Near-resonant absorption in the time-dependent self-consistent field and multiconfigurational self-consistent field approximations , 2001 .

[3]  John F. Stanton,et al.  The equation of motion coupled‐cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties , 1993 .

[4]  M. Linares,et al.  Predicting near-UV electronic circular dichroism in nucleosomal DNA by means of DFT response theory. , 2015, Physical chemistry chemical physics : PCCP.

[5]  Poul Jørgensen,et al.  Response functions from Fourier component variational perturbation theory applied to a time-averaged quasienergy , 1998 .

[6]  Xiaosong Li,et al.  Efficient first-principles electronic dynamics. , 2011, The Journal of chemical physics.

[7]  A. D. McLACHLAN,et al.  Time-Dependent Hartree—Fock Theory for Molecules , 1964 .

[8]  Chao Yang,et al.  New algorithms for iterative matrix‐free eigensolvers in quantum chemistry , 2015, J. Comput. Chem..

[9]  P. Jørgensen,et al.  Polarization propagator methods in atomic and molecular calculations , 1984 .

[10]  G. Scuseria,et al.  An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules , 1998 .

[11]  O. Christiansen,et al.  Communication: A reduced-space algorithm for the solution of the complex linear response equations used in coupled cluster damped response theory. , 2013, The Journal of chemical physics.

[12]  H. Walker.  Implementation of the GMRES method using householder transformations , 1988 .

[13]  Michael J Frisch,et al.  Energy-Specific Linear Response TDHF/TDDFT for Calculating High-Energy Excited States. , 2011, Journal of chemical theory and computation.

[14]  H. Monkhorst,et al.  Calculation of properties with the coupled-cluster method , 2009 .

[15]  K. Lopata,et al.  Accelerated Broadband Spectra Using Transition Dipole Decomposition and Padé Approximants. , 2016, Journal of chemical theory and computation.

[16]  A. Munkholm,et al.  Capturing a Photoexcited Molecular Structure Through Time-Domain X-ray Absorption Fine Structure , 2001, Science.

[17]  O. Christiansen,et al.  Coupled-cluster response theory for near-edge x-ray-absorption fine structure of atoms and molecules , 2012 .

[18]  E. Gross,et al.  Fundamentals of time-dependent density functional theory , 2012 .

[19]  K. Hodgson,et al.  Ligand K‐Edge X‐Ray Absorption Spectroscopic Studies: Metal‐Ligand Covalency in a Series of Transition Metal Tetrachlorides. , 1995 .

[20]  Mauro Stener,et al.  Time dependent density functional theory of core electrons excitations , 2003 .

[21]  N. Govind,et al.  Linear-Response and Real-Time Time-Dependent Density Functional Theory Studies of Core-Level Near-Edge X-Ray Absorption. , 2012, Journal of chemical theory and computation.

[22]  Y. Joly.  Calculating X-ray absorption near-edge structure at very low energy. , 2003, Journal of synchrotron radiation.

[23]  J. Tully,et al.  Ab initio time-resolved density functional theory for lifetimes of excited adsorbate states at metal surfaces , 2007 .

[24]  Frank Neese,et al.  Description of the ground-state covalencies of the bis(dithiolato) transition-metal complexes from X-ray absorption spectroscopy and time-dependent density-functional calculations. , 2007, Chemistry.

[25]  Xiaosong Li,et al.  Energy-Specific Equation-of-Motion Coupled-Cluster Methods for High-Energy Excited States: Application to K-edge X-ray Absorption Spectroscopy. , 2015, Journal of chemical theory and computation.

[26]  Henrik Koch,et al.  Coupled cluster response functions , 1990 .

[27]  R. Morgan.  Generalizations of davidson's method for computing eigenvalues of large nonsymmetric matrices , 1992 .

[28]  Donald C. Comeau,et al.  The equation-of-motion coupled-cluster method. Applications to open- and closed-shell reference states , 1993 .

[29]  Andreas Dreuw,et al.  Calculating X-ray Absorption Spectra of Open-Shell Molecules with the Unrestricted Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator. , 2014, Journal of chemical theory and computation.

[30]  Andreas Dreuw,et al.  Calculating core‐level excitations and x‐ray absorption spectra of medium‐sized closed‐shell molecules with the algebraic‐diagrammatic construction scheme for the polarization propagator , 2014, J. Comput. Chem..

[31]  Diling Zhu,et al.  Ultrafast Excited State Relaxation of a Metalloporphyrin Revealed by Femtosecond X-ray Absorption Spectroscopy. , 2016, Journal of the American Chemical Society.

[32]  V. McKoy,et al.  Some applications of excited-state-excited-state transition densities , 1975 .

[33]  V. Mehrmann,et al.  A numerically stable, structure preserving method for computing the eigenvalues of real Hamiltonian or symplectic pencils , 1998 .

[34]  Thomas Fransson,et al.  Asymmetric-Lanczos-Chain-Driven Implementation of Electronic Resonance Convergent Coupled-Cluster Linear Response Theory. , 2012, Journal of chemical theory and computation.

[35]  Chao Yang,et al.  Accelerating nuclear configuration interaction calculations through a preconditioned block iterative eigensolver , 2016, Comput. Phys. Commun..

[36]  Diling Zhu,et al.  Imaging ultrafast excited state pathways in transition metal complexes by X-ray transient absorption and scattering using X-ray free electron laser source. , 2016, Faraday discussions.

[37]  K. Hodgson,et al.  A Systematic K-edge X-ray Absorption Spectroscopic Study of Cu(III) Sites , 2000 .

[38]  S. Allard,et al.  Molecular Orientation in Polymer Films for Organic Solar Cells Studied by NEXAFS , 2012 .

[39]  David J Rowe,et al.  EQUATIONS-OF-MOTION METHOD AND THE EXTENDED SHELL MODEL. , 1968 .

[40]  P. Norman,et al.  K- and L-edge X-ray absorption spectrum calculations of closed-shell carbon, silicon, germanium, and sulfur compounds using damped four-component density functional response theory. , 2016, Physical chemistry chemical physics : PCCP.

[41]  Poul Jørgensen,et al.  On the Efficiency of Algorithms for Solving Hartree-Fock and Kohn-Sham Response Equations. , 2011, Journal of chemical theory and computation.

[42]  Chao Yang,et al.  BSEPACK User's Guide , 2016, ArXiv.

[43]  Isaiah Shavitt,et al.  Many-Body Methods in Chemistry and Physics: MBPT and Coupled-Cluster Theory , 2009 .

[44]  Zhaojun Bai,et al.  Minimization Principles for the Linear Response Eigenvalue Problem I: Theory , 2012, SIAM J. Matrix Anal. Appl..

[45]  Trygve Helgaker,et al.  Recent advances in wave function-based methods of molecular-property calculations. , 2012, Chemical reviews.

[46]  Mark Earl Casida,et al.  In Recent Advances in Density-Functional Methods , 1995 .

[47]  Chao Yang,et al.  Structure preserving parallel algorithms for solving the Bethe-Salpeter eigenvalue problem , 2015, 1501.03830.

[48]  H. Ågren,et al.  Complex polarization propagator approach in the restricted open-shell, self-consistent field approximation: the near K-edge X-ray absorption fine structure spectra of allyl and copper phthalocyanine. , 2011, The journal of physical chemistry. B.

[49]  Robert A. Harris,et al.  Oscillator Strengths and Rotational Strengths in Hartree–Fock Theory , 1969 .

[50]  Carbon X-ray absorption spectra of fluoroethenes and acetone: a study at the coupled cluster, density functional, and static-exchange levels of theory. , 2013, The Journal of chemical physics.

[51]  Yousef Saad,et al.  Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT. , 2015, Journal of chemical theory and computation.

[52]  R. Morgan,et al.  Generalizations of Davidson's method for computing eigenvalues of sparse symmetric matrices , 1986 .

[53]  H. Ågren,et al.  Polarization propagator for x-ray spectra. , 2006, Physical review letters.

[54]  Vincenzo Barone,et al.  Recent Advances in Density Functional Methods Part III , 2002 .

[55]  J. Olsen,et al.  Solution of the large matrix equations which occur in response theory , 1988 .

[56]  Nicholas A Besley,et al.  Time-dependent density functional theory calculations of the spectroscopy of core electrons. , 2010, Physical chemistry chemical physics : PCCP.

[57]  M. Strayer,et al.  The Nuclear Many-Body Problem , 2004 .

[58]  Patrick Norman,et al.  A Polarization Propagator for Nonlinear X-ray Spectroscopies. , 2016, The journal of physical chemistry letters.

[59]  Dmitri A Romanov,et al.  A time-dependent Hartree-Fock approach for studying the electronic optical response of molecules in intense fields. , 2005, Physical chemistry chemical physics : PCCP.

[60]  Horst Weiss,et al.  A direct algorithm for self‐consistent‐field linear response theory and application to C60: Excitation energies, oscillator strengths, and frequency‐dependent polarizabilities , 1993 .

[61]  Paul,et al.  SECOND QUANTIZATION- BASED METIJODS IN QUANTUM CHEMISTRY , 2003 .

[62]  E. Koch,et al.  Orientation of aromatic hydrocarbons on metal surfaces as determined by NEXAFS , 1987 .

[63]  Andreas Dreuw,et al.  Single-reference ab initio methods for the calculation of excited states of large molecules. , 2005, Chemical reviews.

[64]  P. D. Nguyen,et al.  Calibration of Energy-Specific TDDFT for Modeling K-edge XAS Spectra of Light Elements. , 2015, Journal of chemical theory and computation.

[65]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .