A block variational procedure for the iterative diagonalization of non-Hermitian random-phase approximation matrices.

We present a technique for the iterative diagonalization of random-phase approximation (RPA) matrices, which are encountered in the framework of time-dependent density-functional theory (TDDFT) and the Bethe-Salpeter equation. The non-Hermitian character of these matrices does not permit a straightforward application of standard iterative techniques used, i.e., for the diagonalization of ground state Hamiltonians. We first introduce a new block variational principle for RPA matrices. We then develop an algorithm for the simultaneous calculation of multiple eigenvalues and eigenvectors, with convergence and stability properties similar to techniques used to iteratively diagonalize Hermitian matrices. The algorithm is validated for simple systems (Na(2) and Na(4)) and then used to compute multiple low-lying TDDFT excitation energies of the benzene molecule.

[1]  S. Dancoff.  Non-Adiabatic Meson Theory of Nuclear Forces , 1950 .

[2]  D. Thouless.  Vibrational states of nuclei in the random phase approximation , 1961 .

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

[4]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[5]  Testa,et al.  Green's-function approach to linear response in solids. , 1987, Physical review letters.

[6]  Ed Anderson,et al.  LAPACK Users' Guide , 1995 .

[7]  M. E. Casida.  Time-Dependent Density Functional Response Theory for Molecules , 1995 .

[8]  R. Ahlrichs,et al.  Treatment of electronic excitations within the adiabatic approximation of time dependent density functional theory , 1996 .

[9]  L. Benedict,et al.  Theory of optical absorption in diamond, Si, Ge, and GaAs , 1998 .

[10]  Stefan Albrecht Lucia Reining Rodolfo Del Sole Giovanni Onida.  Ab Initio Calculation of Excitonic Effects in the Optical Spectra of Semiconductors , 1998 .

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

[12]  Jack Dongarra,et al.  LAPACK Users' Guide, 3rd ed. , 1999 .

[13]  E. Tsiper,et al.  LETTER TO THE EDITOR: A classical mechanics technique for quantum linear response , 2000, cond-mat/0003164.

[14]  S. Louie,et al.  Electron-hole excitations and optical spectra from first principles , 2000 .

[15]  N. Rösch,et al.  An efficient method for calculating molecular excitation energies by time-dependent density-functional theory , 2000 .

[16]  Stefano de Gironcoli,et al.  Phonons and related crystal properties from density-functional perturbation theory , 2000, cond-mat/0012092.

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

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

[19]  K. Yabana,et al.  Solving the RPA Eigenvalue Equation in Real-Space , 2002 .

[20]  Friedhelm Bechstedt,et al.  Molecular electronic excitations calculated from a solid-state approach: Methodology and numerics , 2005 .

[21]  R. Gebauer,et al.  Efficient approach to time-dependent density-functional perturbation theory for optical spectroscopy. , 2005, Physical review letters.

[22]  M. L. Tiago,et al.  Optical excitations in organic molecules, clusters, and defects studied by first-principles Green's function methods , 2006, cond-mat/0605248.

[23]  M. Challacombe,et al.  Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory. , 2008, The Journal of chemical physics.

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

[25]  A. Marini,et al.  Exciton-plasmon States in nanoscale materials: breakdown of the Tamm-Dancoff approximation. , 2008, Nano letters.

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

[27]  D. Lu,et al.  Ab initio calculations of optical absorption spectra: solution of the Bethe-Salpeter equation within density matrix perturbation theory. , 2010, The Journal of chemical physics.

[28]  Ralph Gebauer,et al.  turboTDDFT - A code for the simulation of molecular spectra using the Liouville-Lanczos approach to time-dependent density-functional perturbation theory , 2011, Comput. Phys. Commun..