Role of Many-Body Effects in Describing Low-Lying Excited States of π-Conjugated Chromophores: High-Level Equation-of-Motion Coupled-Cluster Studies of Fused Porphyrin Systems.

The unusual photophysical properties of the π-conjugated chromophores make them potential building blocks of various molecular devices. In particular, significant narrowing of the HOMO-LUMO gaps can be observed as an effect of functionalization chromophores with polycyclic aromatic hydrocarbons (PAHs). In this paper we present equation-of-motion coupled cluster (EOMCC) calculations for vertical excitation energies of several functionalized forms of porphyrins. The results for free-base porphyrin (FBP) clearly demonstrate significant differences between functionalization of FBP with one- (anthracene) and two-dimensional (coronene) structures. We also compare the EOMCC results with the experimentally available results for anthracene fused zinc-porphyrin. The impact of various types of correlation effects is illustrated on several benchmark models, where the comparison with the experiment is possible. In particular, we demonstrate that for all excited states considered in this paper, all of them being dominated by single excitations, the inclusion of triply excited configurations is crucial for attaining qualitative agreement with experiment. We also demonstrate the parallel performance of the most computationally intensive part of the completely renormalized EOMCCSD(T) approach (CR-EOMCCSD(T)) across 120 000 cores.

[1]  Robert J. Harrison,et al.  Global arrays: A nonuniform memory access programming model for high-performance computers , 1996, The Journal of Supercomputing.

[2]  Sriram Krishnamoorthy,et al.  Active-space completely-renormalized equation-of-motion coupled-cluster formalism: Excited-state studies of green fluorescent protein, free-base porphyrin, and oligoporphyrin dimer. , 2010, The Journal of chemical physics.

[3]  Rodney J. Bartlett,et al.  The equation-of-motion coupled-cluster method: Excitation energies of Be and CO , 1989 .

[4]  Yukio Kawashima,et al.  Theoretical study of the valence π → π* excited states of polyacenes: anthracene and naphthacene , 1999 .

[5]  S. Hirata.  Tensor Contraction Engine: Abstraction and Automated Parallel Implementation of Configuration-Interaction, Coupled-Cluster, and Many-Body Perturbation Theories , 2003 .

[6]  Akihiko Tsuda,et al.  Fully Conjugated Porphyrin Tapes with Electronic Absorption Bands That Reach into Infrared , 2001, Science.

[7]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[8]  P. Jørgensen,et al.  THE ELECTRONIC SPECTRUM OF FURAN , 1998 .

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

[10]  So Hirata,et al.  Higher-order equation-of-motion coupled-cluster methods. , 2004, The Journal of chemical physics.

[11]  John D. Watts,et al.  Economical triple excitation equation-of-motion coupled-cluster methods for excitation energies , 1995 .

[12]  J. Gauss,et al.  A coupled cluster study of the 1 1A1g and 1 1B2u states of benzene , 1998 .

[13]  Maxwell J. Crossley,et al.  Rigid Fused Oligoporphyrins as Potential Versatile Molecular Wires. 2. B3LYP and SCF Calculated Geometric and Electronic Properties of 98 Oligoporphyrin and Related Molecules , 1999 .

[14]  Poul Jørgensen,et al.  Perturbative triple excitation corrections to coupled cluster singles and doubles excitation energies , 1996 .

[15]  K. Hirao,et al.  Anharmonic vibrational state calculations in the electronic excited states studied by time-dependent density functional theory , 2007 .

[16]  Bryan M. Wong,et al.  Optoelectronic and Excitonic Properties of Oligoacenes: Substantial Improvements from Range-Separated Time-Dependent Density Functional Theory , 2010, Journal of chemical theory and computation.

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

[18]  H. Monkhorst,et al.  Some aspects of the time-dependent coupled-cluster approach to dynamic response functions , 1983 .

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

[20]  K. Kowalski.  Nested variant of the method of moments of coupled cluster equations for vertical excitation energies and excited-state potential energy surfaces. , 2009, The Journal of chemical physics.

[21]  A. Zewail,et al.  Jet spectroscopy of anthracene and deuterated anthracenes , 1984 .

[22]  Kimihiko Hirao,et al.  Multireference Møller—Plesset perturbation theory for high-spin open-shell systems , 1992 .

[23]  Kirk A Peterson,et al.  Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc-Zn. , 2005, The Journal of chemical physics.

[24]  Marat Valiev,et al.  Large-scale parallel calculations with combined coupled cluster and molecular mechanics formalism: Excitation energies of zinc–porphyrin in aqueous solution , 2008 .

[25]  Karina Sendt,et al.  Switchable electronic coupling in model oligoporphyrin molecular wires examined through the measurement and assignment of electronic absorption spectra. , 2002, Journal of the American Chemical Society.

[26]  So Hirata,et al.  Second- and third-order triples and quadruples corrections to coupled-cluster singles and doubles in the ground and excited states. , 2007, The Journal of chemical physics.

[27]  Seigo Ito,et al.  Large pi-aromatic molecules as potential sensitizers for highly efficient dye-sensitized solar cells. , 2009, Accounts of chemical research.

[28]  Ove Christiansen,et al.  Atomic integral driven second order polarization propagator calculations of the excitation spectra of naphthalene and anthracene , 2000 .

[29]  A. Osuka,et al.  Syntheses, structural characterizations, and optical and electrochemical properties of directly fused diporphyrins. , 2001, Journal of the American Chemical Society.

[30]  Klaus Müllen,et al.  Nanographenes as active components of single-molecule electronics and how a scanning tunneling microscope puts them to work. , 2008, Accounts of chemical research.

[31]  P. Jørgensen,et al.  Large-scale calculations of excitation energies in coupled cluster theory: The singlet excited states of benzene , 1996 .

[32]  Piotr Piecuch,et al.  Two new classes of non-iterative coupled-cluster methods derived from the method of moments of coupled-cluster equations , 2006 .

[33]  Niranjan Govind,et al.  Excited states of DNA base pairs using long-range corrected time-dependent density functional theory. , 2009, The journal of physical chemistry. A.

[34]  Marcel Nooijen,et al.  Many‐body similarity transformations generated by normal ordered exponential excitation operators , 1996 .

[35]  K. Kuczera,et al.  Ab initio calculations of S1 excited state vibrational spectra of benzene, naphthalene and anthracene , 1997 .

[36]  P. Piecuch,et al.  New type of noniterative energy corrections for excited electronic states: Extension of the method of moments of coupled-cluster equations to the equation-of-motion coupled-cluster formalism , 2001 .

[37]  R. Bartlett,et al.  Coupled-cluster theory in quantum chemistry , 2007 .

[38]  Piotr Piecuch,et al.  Left‐eigenstate completely renormalized equation‐of‐motion coupled‐cluster methods: Review of key concepts, extension to excited states of open‐shell systems, and comparison with electron‐attached and ionized approaches , 2009 .

[39]  Anna I Krylov,et al.  A noniterative perturbative triples correction for the spin-flipping and spin-conserving equation-of-motion coupled-cluster methods with single and double substitutions. , 2008, The Journal of chemical physics.

[40]  Karol Kowalski,et al.  New coupled-cluster methods with singles, doubles, and noniterative triples for high accuracy calculations of excited electronic states. , 2004, The Journal of chemical physics.

[41]  Kimihiko Hirao,et al.  Multireference Møller-Plesset method , 1992 .

[42]  Hideo Sekino,et al.  A screened potential molecular‐orbital calculation of the π‐electron system of porphyrin , 1981 .

[43]  A. Osuka,et al.  Completely Fused Diporphyrins and Triporphyrin , 2000 .

[44]  Paweł Sałek,et al.  Assessment of a Coulomb-attenuated exchange-correlation energy functional. , 2006, Physical chemistry chemical physics : PCCP.

[45]  Karol Kowalski,et al.  The active-space equation-of-motion coupled-cluster methods for excited electronic states: Full EOMCCSDt , 2001 .

[46]  Ove Christiansen,et al.  Response functions in the CC3 iterative triple excitation model , 1995 .

[47]  John R. Platt,et al.  Classification of Spectra of Cata-Condensed Hydrocarbons , 1949 .

[48]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .

[49]  Harry L Anderson,et al.  Expanding the porphyrin pi-system by fusion with anthracene. , 2008, Organic letters.

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

[51]  Stefan Grimme,et al.  Substantial errors from time-dependent density functional theory for the calculation of excited states of large pi systems. , 2003, Chemphyschem : a European journal of chemical physics and physical chemistry.

[52]  R. Bartlett,et al.  Charge-transfer separability and size-extensivity in the equation-of-motion coupled cluster method: EOM-CCx. , 2011, The Journal of chemical physics.

[53]  V. Kellö,et al.  Medium-size polarized basis sets for high-level-correlated calculations of molecular electric properties , 1991 .

[54]  John E Anthony,et al.  Functionalized acenes and heteroacenes for organic electronics. , 2006, Chemical reviews.

[55]  Stefan Grimme,et al.  A TDDFT study of the lowest excitation energies of polycyclic aromatic hydrocarbons , 2003 .

[56]  Robert J. Harrison,et al.  Global Arrays: a portable "shared-memory" programming model for distributed memory computers , 1994, Proceedings of Supercomputing '94.

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

[58]  Kimihiko Hirao,et al.  Cluster expansion of the wavefunction. Structure of the closed‐shell orbital theory , 1978 .

[59]  John D. Watts,et al.  Iterative and non-iterative triple excitation corrections in coupled-cluster methods for excited electronic states: the EOM-CCSDT-3 and EOM-CCSD(T̃) methods , 1996 .

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

[61]  J. Brédas,et al.  Molecular understanding of organic solar cells: the challenges. , 2009, Accounts of chemical research.

[62]  Anna I. Krylov,et al.  Size-consistent wave functions for bond-breaking: the equation-of-motion spin-flip model , 2001 .

[63]  M. Head‐Gordon,et al.  A fifth-order perturbation comparison of electron correlation theories , 1989 .

[64]  So Hirata,et al.  Perturbative corrections to coupled-cluster and equation-of-motion coupled-cluster energies: A determinantal analysis , 2001 .

[65]  Sriram Krishnamoorthy,et al.  EOMCC, MRPT, and TDDFT studies of charge transfer processes in mixed-valence compounds: application to the spiro molecule. , 2010, The journal of physical chemistry. A.

[66]  Wojciech Pisula,et al.  Graphenes as potential material for electronics. , 2007, Chemical reviews.

[67]  J. Olsen,et al.  Excitation energies of H2O, N2 and C2 in full configuration interaction and coupled cluster theory , 1996 .

[68]  J. Hammond,et al.  Dynamic polarizabilities of polyaromatic hydrocarbons using coupled-cluster linear response theory. , 2007, The Journal of chemical physics.

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