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

We present a special symmetric Lanczos algorithm and a kernel polynomial method (KPM) for approximating the absorption spectrum of molecules within the linear response time-dependent density functional theory (TDDFT) framework in the product form. In contrast to existing algorithms, the new algorithms are based on reformulating the original non-Hermitian eigenvalue problem as a product eigenvalue problem and the observation that the product eigenvalue problem is self-adjoint with respect to an appropriately chosen inner product. This allows a simple symmetric Lanczos algorithm to be used to compute the desired absorption spectrum. The use of a symmetric Lanczos algorithm only requires half of the memory compared with the nonsymmetric variant of the Lanczos algorithm. The symmetric Lanczos algorithm is also numerically more stable than the nonsymmetric version. The KPM algorithm is also presented as a low-memory alternative to the Lanczos approach, but the algorithm may require more matrix-vector multiplications in practice. We discuss the pros and cons of these methods in terms of their accuracy as well as their computational and storage cost. Applications to a set of small and medium-sized molecules are also presented.

[1]  Christopher C. Paige,et al.  The computation of eigenvalues and eigenvectors of very large sparse matrices , 1971 .

[2]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[3]  A. Voter,et al.  Kernel Polynomial Approximations for Densities of States and Spectral Functions , 1996 .

[4]  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.

[5]  Maurice Cohen,et al.  HARTREE–FOCK WAVE FUNCTIONS FOR EXCITED STATES: II. SIMPLIFICATION OF THE ORBITAL EQUATIONS , 1966 .

[6]  Mykhaylo Krykunov,et al.  The implementation of a self-consistent constricted variational density functional theory for the description of excited states. , 2012, The Journal of chemical physics.

[7]  L Jensen,et al.  Finite lifetime effects on the polarizability within time-dependent density-functional theory. , 2005, The Journal of chemical physics.

[8]  Stefan Grimme,et al.  A simplified Tamm-Dancoff density functional approach for the electronic excitation spectra of very large molecules. , 2013, The Journal of chemical physics.

[9]  Patrick Norman,et al.  Efficient Calculations of Molecular Linear Response Properties for Spectral Regions. , 2014, Journal of chemical theory and computation.

[10]  Niranjan Govind,et al.  Modeling Fast Electron Dynamics with Real-Time Time-Dependent Density Functional Theory: Application to Small Molecules and Chromophores. , 2011, Journal of chemical theory and computation.

[11]  A. B. Gordienko,et al.  Kernel polynomial method as an efficient O ( N2 ) scheme for optical spectra calculations including electron–hole interaction , 2014 .

[12]  Benjamin T. Miller,et al.  A parallel implementation of the analytic nuclear gradient for time-dependent density functional theory within the Tamm–Dancoff approximation , 1999 .

[13]  Y. Saad,et al.  Turbo charging time-dependent density-functional theory with Lanczos chains. , 2006, The Journal of chemical physics.

[14]  Jiaya Jia,et al.  To appear in , 2004 .

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

[16]  Tjerk P. Straatsma,et al.  NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations , 2010, Comput. Phys. Commun..

[17]  L. Reining,et al.  Electronic excitations: density-functional versus many-body Green's-function approaches , 2002 .

[18]  Alexander Gaenko,et al.  Two-component relativistic density functional method for computing nonsingular complex linear response of molecules based on the zeroth order regular approximation. , 2009, The Journal of chemical physics.

[19]  P. Joergensen,et al.  Second Quantization-based Methods in Quantum Chemistry , 1981 .

[20]  Sergei Tretiak,et al.  Representation independent algorithms for molecular response calculations in time-dependent self-consistent field theories. , 2009, The Journal of chemical physics.

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

[22]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[23]  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.

[24]  Á. Rubio,et al.  Time-dependent density-functional theory. , 2009, Physical chemistry chemical physics : PCCP.

[25]  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.

[26]  D. Jackson,et al.  The theory of approximation , 1982 .

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

[28]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[29]  Yong Wang,et al.  Optical Absorption and Band Gap Reduction in (Fe1–xCrx)2O3 Solid Solutions: A First-Principles Study , 2013 .

[30]  D. Chong.  Recent Advances in Density Functional Methods Part III , 2002 .

[31]  Jack H. Freed,et al.  Classical time‐correlation functions and the Lanczos algorithm , 1981 .

[32]  Matt Challacombe,et al.  Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations , 2010, Comput..

[33]  E. Gross,et al.  Density-Functional Theory for Time-Dependent Systems , 1984 .

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

[35]  Stefan Grimme,et al.  A simplified time-dependent density functional theory approach for electronic ultraviolet and circular dichroism spectra of very large molecules , 2014 .

[36]  J. Chelikowsky,et al.  Ab initio absorption spectra and optical gaps in nanocrystalline silicon. , 2001, Physical review letters.

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

[38]  Zhaojun Bai,et al.  Minimization Principles for the Linear Response Eigenvalue Problem II: Computation , 2013, SIAM J. Matrix Anal. Appl..

[39]  Roi Baer,et al.  Efficient linear-response method circumventing the exchange-correlation kernel: theory for molecular conductance under finite bias. , 2005, The Journal of chemical physics.

[40]  Niranjan Govind,et al.  Comparison of Real-Time and Linear-Response Time-Dependent Density Functional Theories for Molecular Chromophores Ranging from Sparse to High Densities of States. , 2015, Journal of chemical theory and computation.

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

[42]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[43]  E. Gross,et al.  Time-dependent density functional theory. , 2004, Annual review of physical chemistry.

[44]  Feliciano Giustino,et al.  Linear optical response of finite systems using multishift linear system solvers. , 2014, The Journal of chemical physics.