Large-scale parallel calculations with combined coupled cluster and molecular mechanics formalism: Excitation energies of zinc–porphyrin in aqueous solution

Abstract The need for combined QM/MM approaches with high-level coupled cluster (CC) methods in the QM core is unquestionable. In an effort to make these methods widely applicable, the scalability of CC approaches across hundreds/thousands of CPUs will play a central role. In this Letter we report the performance of our CC codes in calculations of vertical excitation energies of low-lying excited states of the zinc–porphyrin molecule in aqueous solution. Coupled cluster description of excited states is based on equation-of-motion approach with singles and doubles (EOMCCSD) as well as its non-iterative extension for triply excited configurations. These results are compared with those obtained with time-dependent density functional theory (TD-DFT).

[1]  Marat Valiev,et al.  Hybrid coupled cluster and molecular dynamics approach: application to the excitation spectrum of cytosine in the native DNA environment. , 2006, The Journal of chemical physics.

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

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

[4]  Rodney J. Bartlett,et al.  Similarity transformed equation-of-motion coupled-cluster theory: Details, examples, and comparisons , 1997 .

[5]  J. Cizek.  On the Correlation Problem in Atomic and Molecular Systems. Calculation of Wavefunction Components in Ursell-Type Expansion Using Quantum-Field Theoretical Methods , 1966 .

[6]  K. Kowalski.  Excitation energies through the locally renormalized equation-of-motion formalism: singles and doubles model. , 2006, The Journal of chemical physics.

[7]  A. Ellis,et al.  Second-order overtone and combination modes in the LOLA region of acid treated double-walled carbon nanotubes. , 2006, The Journal of chemical physics.

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

[9]  N. Nakatsuji,et al.  Cluster expansion of the wavefunction. Excited states , 1978 .

[10]  Marat Valiev,et al.  Phosphorylation reaction in cAPK protein kinase-free energy quantum mechanical/molecular mechanics simulations. , 2007 .

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

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

[13]  F. Coester,et al.  Bound states of a many-particle system , 1958 .

[14]  A. Osuka,et al.  Doubly meso‐β‐Linked Diporphyrins from Oxidation of 5,10,15‐Triaryl‐Substituted NiII– and PdII – Porphyrins , 2000 .

[15]  Marat Valiev,et al.  Asymptotic extrapolation scheme for large-scale calculations with hybrid coupled cluster and molecular dynamics simulations. , 2006, The journal of physical chemistry. A.

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

[17]  Anna I. Krylov,et al.  Spin-flip configuration interaction: an electronic structure model that is both variational and size-consistent , 2001 .

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

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

[20]  T. H. Dunning.  Gaussian Basis Functions for Use in Molecular Calculations. III. Contraction of (10s6p) Atomic Basis Sets for the First‐Row Atoms , 1970 .

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

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

[23]  Tomoo Miyahara,et al.  Ground and excited states of linked and fused zinc porphyrin dimers: Symmetry adapted cluster (SAC)—configuration interaction (CI) study , 2002 .

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

[25]  A. Burrell,et al.  Synthetic routes to multiporphyrin arrays. , 2001, Chemical reviews.

[26]  T. Ogawa,et al.  Completely Regioselective Synthesis of Directly Linkedmeso,meso andmeso,β Porphyrin Dimers by One-Pot Electrochemical Oxidation of Metalloporphyrins , 1999 .

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

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

[29]  Karol Kowalski,et al.  The method of moments of coupled-cluster equations and the renormalized CCSD[T], CCSD(T), CCSD(TQ), and CCSDT(Q) approaches , 2000 .

[30]  Marat Valiev,et al.  Hybrid approach for free energy calculations with high-level methods: application to the SN2 reaction of CHCl3 and OH- in water. , 2007, The Journal of chemical physics.

[31]  F. Coester,et al.  Short-range correlations in nuclear wave functions , 1960 .

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

[33]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[34]  Josef Paldus,et al.  Correlation problems in atomic and molecular systems III. Rederivation of the coupled-pair many-electron theory using the traditional quantum chemical methodst†‡§ , 1971 .

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

[36]  Michael R. Wasielewski,et al.  Photoinduced electron transfer in supramolecular systems for artificial photosynthesis , 1992 .

[37]  Rodney J. Bartlett,et al.  A new method for excited states: Similarity transformed equation-of-motion coupled-cluster theory , 1997 .

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

[39]  R. Bartlett,et al.  A coupled‐cluster based effective Hamiltonian method for dynamic electric polarizabilities , 1993 .

[40]  Michael J. McGuire,et al.  Recent advances in electronic structure theory: Method of moments of coupled-cluster equations and renormalized coupled-cluster approaches , 2002 .

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

[42]  K. Kowalski.  Completely Renormalized EOMCCSD(T) Method Employing Independent Optimization of the Cluster Product Terms , 2005 .

[43]  Kimihiko Hirao,et al.  Cluster expansion of the wavefunction. Symmetry-adapted-cluster expansion, its variational determination, and extension of open-shell orbital theory , 1978 .

[44]  Hiroshi Nakatsuji,et al.  Cluster expansion of the wavefunction. Calculation of electron correlations in ground and excited states by SAC and SAC CI theories , 1979 .

[45]  Hiroshi Nakatsuji,et al.  Cluster expansion of the wavefunction. Electron correlations in ground and excited states by SAC (symmetry-adapted-cluster) and SAC CI theories , 1979 .