Equation of motion coupled cluster methods for electron attachment and ionization potential in fullerenes C60 and C70.

In both molecular and periodic solid-state systems there is a need for the accurate determination of the ionization potential and the electron affinity for systems ranging from light harvesting polymers and photocatalytic compounds to semiconductors. The development of a Green's function approach based on the coupled cluster (CC) formalism would be a valuable tool for addressing many properties involving many-body interactions along with their associated correlation functions. As a first step in this direction, we have developed an accurate and parallel efficient approach based on the equation of motion-CC technique. To demonstrate the high degree of accuracy and numerical efficiency of our approach we calculate the ionization potential and electron affinity for C60 and C70. Accurate predictions for these molecules are well beyond traditional molecular scale studies. We compare our results with experiments and both quantum Monte Carlo and GW calculations.

[1]  H. Nakatsuji,et al.  Symmetry‐adapted cluster–configuration interaction method applied to high‐spin multiplicity , 1993 .

[2]  C. Deibel,et al.  Fullerene dimers (C60/C70) for energy harvesting. , 2009, Chemistry.

[3]  Anna I Krylov,et al.  Equation-of-motion coupled-cluster methods for open-shell and electronically excited species: the Hitchhiker's guide to Fock space. , 2008, Annual review of physical chemistry.

[4]  R. Bartlett,et al.  Equation-of-motion coupled cluster method with full inclusion of connected triple excitations for electron-attached states: EA-EOM-CCSDT , 2003 .

[5]  Rodney J. Bartlett,et al.  Equation-of-motion coupled cluster method with full inclusion of the connected triple excitations for ionized states: IP-EOM-CCSDT , 2003 .

[6]  Dmitry I. Lyakh,et al.  Potential energy curves via double electron-attachment calculations: dissociation of alkali metal dimers. , 2013, The Journal of chemical physics.

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

[8]  Karol Kowalski,et al.  Coupled cluster calculations for static and dynamic polarizabilities of C60. , 2008, The Journal of chemical physics.

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

[10]  S. Hirata,et al.  Higher-order equation-of-motion coupled-cluster methods for electron attachment. , 2007, The Journal of chemical physics.

[11]  N. S. Sariciftci,et al.  Conjugated polymer-based organic solar cells. , 2007, Chemical reviews.

[12]  Rodney J. Bartlett,et al.  Equation of motion coupled cluster method for electron attachment , 1995 .

[13]  S. Hirata,et al.  Higher-order equation-of-motion coupled-cluster methods for ionization processes. , 2006, The Journal of chemical physics.

[14]  M. Nooijen,et al.  Application of double ionization state-specific equation of motion coupled cluster method to organic diradicals. , 2008, The journal of physical chemistry. A.

[15]  Lai‐Sheng Wang,et al.  High resolution photoelectron spectroscopy of C60 , 1999 .

[16]  Anna I Krylov,et al.  Using the charge-stabilization technique in the double ionization potential equation-of-motion calculations with dianion references. , 2011, The Journal of chemical physics.

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

[18]  Christoph J. Brabec,et al.  Design Rules for Donors in Bulk‐Heterojunction Solar Cells—Towards 10 % Energy‐Conversion Efficiency , 2006 .

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

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

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

[22]  John F. Stanton,et al.  Many‐body methods for excited state potential energy surfaces. I. General theory of energy gradients for the equation‐of‐motion coupled‐cluster method , 1993 .

[23]  G. Scuseria,et al.  Calculation of spectra and spin–spin coupling constants using a coupled–cluster polarization propagator method , 1987 .

[24]  K. Hirao,et al.  A generalization of the Davidson's method to large nonsymmetric eigenvalue problems , 1982 .

[25]  J. Chelikowsky,et al.  Structure selection based on high vertical electron affinity for TiO2 clusters. , 2012, Physical review letters.

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

[27]  Uttam Sinha Mahapatra,et al.  A size-consistent state-specific multireference coupled cluster theory: Formal developments and molecular applications , 1999 .

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

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

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

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

[32]  Anna I Krylov,et al.  Perturbative triples correction for the equation-of-motion coupled-cluster wave functions with single and double substitutions for ionized states: Theory, implementation, and examples. , 2009, The Journal of chemical physics.

[33]  Edoardo Aprà,et al.  Implementation of the multireference Brillouin-Wigner and Mukherjee's coupled cluster methods with non-iterative triple excitations utilizing reference-level parallelism. , 2012, The Journal of chemical physics.

[34]  K. Ruud,et al.  The dispersion of the polarizability of C60: A confirmation of recent experimental results through theoretical calculations , 2001 .

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

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

[37]  R J Bartlett,et al.  Parallel implementation of electronic structure energy, gradient, and Hessian calculations. , 2008, The Journal of chemical physics.

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

[39]  D. Lichtenberger,et al.  The He I valence photoelectron spectrum ofC70 in the gas phase , 1992 .

[40]  Marcel Nooijen,et al.  Extended similarity transformed equation-of-motion coupled cluster theory (extended-STEOM-CC): Applications to doubly excited states and transition metal compounds , 2000 .

[41]  Rodney J. Bartlett,et al.  Similarity transformed equation-of-motion coupled-cluster study of ionized, electron attached, and excited states of free base porphin , 1997 .

[42]  Wei Zhao,et al.  Electron affinities of 1,1-diaryl-2,3,4,5-tetraphenylsiloles: direct measurements and comparison with experimental and theoretical estimates. , 2005, Journal of the American Chemical Society.

[43]  H. Monkhorst,et al.  Calculation of properties with the coupled-cluster method , 2009 .

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

[45]  Jarek Nieplocha,et al.  Advances, Applications and Performance of the Global Arrays Shared Memory Programming Toolkit , 2006, Int. J. High Perform. Comput. Appl..

[46]  P. Jørgensen,et al.  Frequency dependent polarizabilities in a multiconfigurational time dependent Hartree-Fock approximation , 1980 .

[47]  P. Piecuch,et al.  Efficient formulation and computer implementation of the active-space electron-attached and ionized equation-of-motion coupled-cluster methods. , 2006, The Journal of chemical physics.

[48]  J. V. Ortiz.  The Electron Propagator Picture of Molecular Electronic Structure , 1997 .

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

[50]  S. Baroni,et al.  Can We Tune the Band Offset at Semiconductor Heterojunctions , 1989 .

[51]  A. Heeger.  Nobel Lecture: Semiconducting and metallic polymers: The fourth generation of polymeric materials* , 2001 .

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

[53]  S. Forrest,et al.  Molecular and morphological influences on the open circuit voltages of organic photovoltaic devices. , 2009, Journal of the American Chemical Society.

[54]  Karol Kowalski,et al.  Excited-state potential energy curves of CH+: a comparison of the EOMCCSDt and full EOMCCSDT results , 2001 .

[55]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[56]  K. Emrich,et al.  An extension of the coupled cluster formalism to excited states (I) , 1981 .

[57]  S. Sze,et al.  Physics of Semiconductor Devices: Sze/Physics , 2006 .

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

[59]  R. Smalley,et al.  UPS of 2–30-atom carbon clusters: Chains and rings , 1988 .

[60]  Vladan Stevanović,et al.  Assessing capability of semiconductors to split water using ionization potentials and electron affinities only. , 2014, Physical chemistry chemical physics : PCCP.

[61]  O. Inganäs,et al.  Mixed C60/C70 based fullerene acceptors in polymer bulk-heterojunction solar cells , 2012 .

[62]  Approximate Inclusion of the T3 and R3 Operators in the Equation-of-motion Coupled-cluster Method , 2004 .

[63]  M. Musiał.  Multireference Fock space coupled cluster method in the effective and intermediate Hamiltonian formulation for the (2,0) sector. , 2012, The Journal of chemical physics.

[64]  R. Needs,et al.  Quantum Monte Carlo simulations of solids , 2001 .

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

[66]  G. Scuseria.  Ab initio theoretical predictions of the equilibrium geometries of C60, C60H60 and C60F60 , 1991 .

[67]  John F. Stanton,et al.  Perturbative treatment of the similarity transformed Hamiltonian in equation‐of‐motion coupled‐cluster approximations , 1995 .

[68]  Steven E. J. Bell,et al.  Reduced‐size polarized basis sets for calculations of molecular electric properties. I. The basis set generation , 2005, J. Comput. Chem..

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

[70]  Lai‐Sheng Wang,et al.  Direct experimental probe of the on-site Coulomb repulsion in the doubly charged fullerene anion C70 2-. , 2006, Physical review letters.

[71]  D. Mukherjee.  Applied many-body methods in spectroscopy and electronic structure , 1992 .

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

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

[74]  J. Fréchet,et al.  Polymer-fullerene composite solar cells. , 2008, Angewandte Chemie.

[75]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

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

[77]  Kenneth W. Nebesny,et al.  Valence and core photoelectron spectroscopy of C60, buckminsterfullerene , 1991 .

[78]  W. Krätschmer,et al.  Solid C60: a new form of carbon , 1990, Nature.

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

[80]  Jenny Nelson,et al.  Polymer:fullerene bulk heterojunction solar cells , 2011 .

[81]  S. Barlow,et al.  Synthesis, ionisation potentials and electron affinities of hexaazatrinaphthylene derivatives. , 2007, Chemistry.

[82]  Active-Space Equation-of-Motion Coupled-Cluster Methods through Quadruples for Excited, Ionized, and Electron-Attached States. , 2007, Journal of chemical theory and computation.

[83]  James B. Anderson,et al.  A random‐walk simulation of the Schrödinger equation: H+3 , 1975 .

[84]  P R C Kent,et al.  Neutral and charged excitations in carbon fullerenes from first-principles many-body theories. , 2008, The Journal of chemical physics.

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

[86]  Josef Paldus,et al.  Correlation Problems in Atomic and Molecular Systems. IV. Extended Coupled-Pair Many-Electron Theory and Its Application to the B H 3 Molecule , 1972 .

[87]  Oreste Villa,et al.  Noniterative Multireference Coupled Cluster Methods on Heterogeneous CPU-GPU Systems. , 2013, Journal of chemical theory and computation.

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

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

[90]  P. Szalay,et al.  Analytic energy derivatives for coupled‐cluster methods describing excited states: General formulas and comparison of computational costs , 1995 .

[91]  Jun Shen,et al.  Doubly electron-attached and doubly ionized equation-of-motion coupled-cluster methods with 4-particle-2-hole and 4-hole-2-particle excitations and their active-space extensions. , 2013, The Journal of chemical physics.

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

[93]  William Weltner,et al.  Carbon molecules, ions, and clusters , 1989 .