Calculation of two-photon absorption strengths with the approximate coupled cluster singles and doubles model CC2 using the resolution-of-identity approximation.

An implementation of two-photon absorption matrix elements using the approximate second-order coupled-cluster singles and doubles model CC2 is presented. In this implementation we use the resolution-of-the-identity approximation for the two-electron repulsion integrals to reduce the computational cost. To avoid storage of large arrays we introduce in addition a numerical Laplace transformation of orbital energy denominators for the response of the doubles amplitudes. The error due to the numerical Laplace transformation is found to be negligible. Using this new implementation, we performed a series of benchmark calculations on substituted benzene and azobenzene derivatives to get reference values for TD-DFT results. We show that results obtained with the Coulomb-attenuated B3LYP functional are in reasonable agreement with the coupled-cluster results, whereas other density functionals which do not have a long-range correction give considerably less accurate results. Applications to the AF240 dye molecule and a weakly bound molecular tweezer complex demonstrate that this new RI-CC2 implementation allows for the first time to compute two-photon absorption cross sections with a correlated wave function method for molecules with more than 70 atoms and to apply this method for benchmarking TD-DFT calculations on molecules which are of particular relevance for experimental studies of two-photon absorption.

[1]  R. Birge,et al.  A theoretical analysis of the two‐photon properties of linear polyenes and the visual chromophores , 1979 .

[2]  J. C. Slater.  A Simplification of the Hartree-Fock Method , 1951 .

[3]  K. Ruud,et al.  Large two-photon absorption cross section: molecular tweezer as a new promising class of compounds for nonlinear optics. , 2009, Physical chemistry chemical physics : PCCP.

[4]  Trygve Helgaker,et al.  Analytical Calculation of Geometrical Derivatives in Molecular Electronic Structure Theory , 1988 .

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

[6]  L. Misoguti,et al.  Degenerate two-photon absorption spectra in azoaromatic compounds. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.

[7]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[8]  Christof Hättig,et al.  Analytic gradients for excited states in the coupled-cluster model CC2 employing the resolution-of-the-identity approximation , 2003 .

[9]  C. Garrett,et al.  Two-Photon Excitation in CaF 2 : Eu 2+ , 1961 .

[10]  Luis Serrano-Andrés,et al.  Does density functional theory contribute to the understanding of excited states of unsaturated organic compounds , 1999 .

[11]  Paweł Sałek,et al.  Benchmarking two-photon absorption with CC3 quadratic response theory, and comparison with density-functional response theory. , 2006, The Journal of chemical physics.

[12]  P. Dirac.  Quantum Mechanics of Many-Electron Systems , 1929 .

[13]  M. Head‐Gordon,et al.  Long-range charge-transfer excited states in time-dependent density functional theory require non-local exchange , 2003 .

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

[15]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

[16]  Thomas Bondo Pedersen,et al.  Polarizability and optical rotation calculated from the approximate coupled cluster singles and doubles CC2 linear response theory using Cholesky decompositions. , 2004, The Journal of chemical physics.

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

[18]  Christof Hättig,et al.  Transition moments and excited-state first-order properties in the coupled-cluster model CC2 using the resolution-of-the-identity approximation , 2002 .

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

[20]  Maria Goeppert-Mayer.  Über Elementarakte mit zwei Quantensprüngen , 1931 .

[21]  N. Handy,et al.  A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP) , 2004 .

[22]  Timothy Clark,et al.  Efficient diffuse function‐augmented basis sets for anion calculations. III. The 3‐21+G basis set for first‐row elements, Li–F , 1983 .

[23]  P. Krawczyk,et al.  Theoretical study of one- and two-photon absorption spectra of azoaromatic compounds. , 2009, The Journal of chemical physics.

[24]  Poul Jørgensen,et al.  The second-order approximate coupled cluster singles and doubles model CC2 , 1995 .

[25]  Christof Hättig,et al.  CC2 excitation energy calculations on large molecules using the resolution of the identity approximation , 2000 .

[26]  C. Hättig,et al.  Scaled opposite-spin CC2 for ground and excited states with fourth order scaling computational costs. , 2011, The Journal of chemical physics.

[27]  Paras N. Prasad,et al.  ERRATUM: Nonlinear multiphoton processes in organic and polymeric materials , 1996 .

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

[29]  M. Head‐Gordon,et al.  Failure of time-dependent density functional theory for long-range charge-transfer excited states: the zincbacteriochlorin-bacteriochlorin and bacteriochlorophyll-spheroidene complexes. , 2004, Journal of the American Chemical Society.

[30]  H. Ågren,et al.  Calculations of two-photon absorption cross sections by means of density-functional theory , 2003 .

[31]  Mark S. Gordon,et al.  Self‐consistent molecular orbital methods. XXIII. A polarization‐type basis set for second‐row elements , 1982 .

[32]  Ove Christiansen,et al.  Multiphoton transition moments and absorption cross sections in coupled cluster response theory employing variational transition moment functionals , 1998 .

[33]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[34]  Paweł Sałek,et al.  Calculations of two-photon charge-transfer excitations using Coulomb-attenuated density-functional theory. , 2005, The Journal of chemical physics.

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

[36]  F. Weigend,et al.  Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations , 2002 .

[37]  M. Head‐Gordon,et al.  Chlorophyll fluorescence quenching by xanthophylls , 2003 .

[38]  Ann G. Dombroskie,et al.  Degenerate two-photon-absorption spectral studies of highly two-photon active organic chromophores. , 2004, The Journal of chemical physics.

[39]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

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