Vibrational absorption spectra from vibrational coupled cluster damped linear response functions calculated using an asymmetric Lanczos algorithm.

We report the theory and implementation of vibrational coupled cluster (VCC) damped response functions. From the imaginary part of the damped VCC response function the absorption as function of frequency can be obtained, requiring formally the solution of the now complex VCC response equations. The absorption spectrum can in this formulation be seen as a matrix function of the characteristic VCC Jacobian response matrix. The asymmetric matrix version of the Lanczos method is used to generate a tridiagonal representation of the VCC response Jacobian. Solving the complex response equations in the relevant Lanczos space provides a method for calculating the VCC damped response functions and thereby subsequently the absorption spectra. The convergence behaviour of the algorithm is discussed theoretically and tested for different levels of completeness of the VCC expansion. Comparison is made with results from the recently reported [P. Seidler, M. B. Hansen, W. Györffy, D. Toffoli, and O. Christiansen, J. Chem. Phys. 132, 164105 (2010)] vibrational configuration interaction damped response function calculated using a symmetric Lanczos algorithm. Calculations of IR spectra of oxazole, cyclopropene, and uracil illustrate the usefulness of the new VCC based method.

[1]  A. Császár,et al.  On the variational computation of a large number of vibrational energy levels and wave functions for medium-sized molecules. , 2009, The Journal of chemical physics.

[2]  O. Christiansen.  Response theory for vibrational wave functions. , 2005, The Journal of chemical physics.

[3]  Barbara H. Stuart,et al.  Infrared Spectroscopy: Fundamentals and Applications: Stuart/Infrared Spectroscopy: Fundamentals and Applications , 2005 .

[4]  Ove Christiansen,et al.  A Lanczos-chain driven approach for calculating damped vibrational configuration interaction response functions. , 2010, The Journal of chemical physics.

[5]  O. Christiansen,et al.  Vibrational coupled cluster response theory: a general implementation. , 2011, The Journal of chemical physics.

[6]  Lorenz S. Cederbaum,et al.  Multimode Molecular Dynamics Beyond the Born‐Oppenheimer Approximation , 2007 .

[7]  Claude Pouchan,et al.  The VCI-P code: an iterative variation–perturbation scheme for efficient computations of anharmonic vibrational levels and IR intensities of polyatomic molecules , 2010 .

[8]  O. Christiansen,et al.  A hierarchy of potential energy surfaces constructed from energies and energy derivatives calculated on grids. , 2009, The Journal of chemical physics.

[9]  Hans-Dieter Meyer,et al.  A band‐Lanczos method for computing matrix elements of a resolvent , 1989 .

[10]  Ove Christiansen,et al.  Towards fast computations of correlated vibrational wave functions: vibrational coupled cluster response excitation energies at the two-mode coupling level. , 2008, The Journal of chemical physics.

[11]  Ove Christiansen,et al.  An adaptive density-guided approach for the generation of potential energy surfaces of polyatomic molecules , 2009 .

[12]  R. Haydock.  The recursive solution of the Schrödinger equation , 1980 .

[13]  T. Carrington,et al.  A comparison of filter diagonalisation methods with the Lanczos method for calculating vibrational energy levels , 1999 .

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

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

[16]  T. Carrington,et al.  A new iterative method for calculating energy levels and wave functions , 2000 .

[17]  Tucker Carrington,et al.  Accelerating the calculation of energy levels and wave functions using an efficient preconditioner with the inexact spectral transform method , 2001 .

[18]  P. Cassam-Chenaï,et al.  The VMFCI method: A flexible tool for solving the molecular vibration problem , 2006, J. Comput. Chem..

[19]  Ove Christiansen,et al.  Vibrational coupled cluster theory. , 2004, The Journal of chemical physics.

[20]  Hua Guo,et al.  A single Lanczos propagation method for calculating transition amplitudes , 1999 .

[21]  M. Palmer.  Comparison of theoretical and experimental studies of infrared spectral data for the 5-membered ring heterocycles , 2007 .

[22]  R. W. Larsen,et al.  High-resolution infrared and theoretical study of gaseous oxazole in the 600–1400 cm−1 region , 2007 .

[23]  O. Christiansen,et al.  Potential energy surfaces for vibrational structure calculations from a multiresolution adaptive density-guided approach: implementation and test calculations. , 2009, The journal of physical chemistry. A.

[24]  N. Handy,et al.  Extensions and tests of “multimode”: a code to obtain accurate vibration/rotation energies of many-mode molecules , 1998 .

[25]  J. Kongsted,et al.  Automatic generation of potential energy and property surfaces of polyatomic molecules in normal coordinates. , 2007, The Journal of chemical physics.

[26]  O. Christiansen.  First-order nonadiabatic coupling matrix elements using coupled cluster methods. I. Theory , 1999 .

[27]  C. Leforestier,et al.  Calculation of highly excited vibrational levels: a prediagonalized Davidson scheme , 2002 .

[28]  J. A. R. Fris,et al.  Time evolution of dynamic propensity in a model glass former: the interplay between structure and dynamics. , 2008, The Journal of chemical physics.

[29]  H. O. Karlsson.  Reduced-order modeling, error estimation, and the role of the start-vector: the recursive residue generation method revisited. , 2007, The journal of physical chemistry. A.

[30]  Peter F. Bernath,et al.  The infrared spectra of uracil, thymine, and adenine in the gas phase , 1997 .

[31]  Poul Jørgensen,et al.  Quasienergy formulation of damped response theory. , 2009, The Journal of chemical physics.

[32]  Ove Christiansen,et al.  Vibrational absorption spectra calculated from vibrational configuration interaction response theory using the Lanczos method. , 2010, The Journal of chemical physics.

[33]  Ove Christiansen,et al.  Vibrational excitation energies from vibrational coupled cluster response theory. , 2007, The Journal of chemical physics.

[34]  G. Moro,et al.  Calculation of ESR spectra and related Fokker–Planck forms by the use of the Lanczos algorithm , 1981 .

[35]  A. Zecchina,et al.  Infra-red spectra of oxazole and its alkyl derivatives—I , 1966 .

[36]  J. Olsen,et al.  Linear and nonlinear response functions for an exact state and for an MCSCF state , 1985 .

[37]  Hans-Dieter Meyer,et al.  Calculation and selective population of vibrational levels with the Multiconfiguration Time-Dependent Hartree (MCTDH) algorithm , 2006 .

[38]  H. Meyer,et al.  Block Lanczos and many-body theory: Application to the one-particle Green's function , 1996 .

[39]  T. Carrington,et al.  Vibrational energy levels of CH5(+). , 2008, The Journal of chemical physics.

[40]  Yohann Scribano,et al.  Iterative active-space selection for vibrational configuration interaction calculations using a reduced-coupling VSCF basis , 2008 .

[41]  R. Whitehead.  Moment Methods and Lanczos Methods , 1980 .

[42]  So Hirata,et al.  Vibrational quasi-degenerate perturbation theory: applications to fermi resonance in CO2, H2CO, and C6H6. , 2008, Physical chemistry chemical physics : PCCP.

[43]  R. Wyatt.  The Recursive Residue Generation Method , 2007 .

[44]  T. Carrington,et al.  A symmetry-adapted Lanczos method for calculating energy levels with different symmetries from a single set of iterations , 2001 .

[45]  W. Haxton,et al.  Piecewise moments method: Generalized Lanczos technique for nuclear response surfaces , 2005, nucl-th/0508034.

[46]  T. Crawford,et al.  An Introduction to Coupled Cluster Theory for Computational Chemists , 2007 .

[47]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

[48]  Daniel Neuhauser,et al.  Extraction, through filter‐diagonalization, of general quantum eigenvalues or classical normal mode frequencies from a small number of residues or a short‐time segment of a signal. I. Theory and application to a quantum‐dynamics model , 1995 .

[49]  R. L. Erskine,et al.  Vibrational Spectrum of Cyclopropene , 1967 .

[50]  C. Lanczos.  An iteration method for the solution of the eigenvalue problem of linear differential and integral operators , 1950 .

[51]  Patrick Norman,et al.  A perspective on nonresonant and resonant electronic response theory for time-dependent molecular properties. , 2011, Physical chemistry chemical physics : PCCP.

[52]  N. Matsunaga,et al.  Degenerate perturbation theory corrections for the vibrational self-consistent field approximation: Method and applications , 2002 .

[53]  L. Adamowicz,et al.  The infrared spectra of matrix isolated uracil and thymine: An assignment based on new theoretical calculations , 1992 .

[54]  Gábor Czakó,et al.  Vibrational energy levels with arbitrary potentials using the Eckart-Watson Hamiltonians and the discrete variable representation. , 2007, The Journal of chemical physics.

[55]  Eduard Matito,et al.  Using Electronic Energy Derivative Information in Automated Potential Energy Surface Construction for Vibrational Calculations. , 2010, Journal of chemical theory and computation.

[56]  Guntram Rauhut,et al.  Configuration selection within vibrational multiconfiguration self-consistent field theory: application to bridged lithium compounds. , 2011, The Journal of chemical physics.

[57]  D. M. Bishop,et al.  Nonlinear response theory with relaxation: the first-order hyperpolarizability. , 2005, The Journal of chemical physics.

[58]  V. Vaida,et al.  Photolysis of Sulfuric Acid Vapor by Visible Solar Radiation , 2003, Science.

[59]  Ove Christiansen,et al.  Automatic derivation and evaluation of vibrational coupled cluster theory equations. , 2009, The Journal of chemical physics.

[60]  E. A. Dorko,et al.  Vibrational spectra of cyclopropene and cyclopropene-1,2-d2 , 1968 .

[61]  Hua Guo,et al.  Effect of spectral range on convergence in Lanczos algorithm, a numerical study , 2003 .