Decomposition of Intermolecular Interaction Energies within the Local Pair Natural Orbital Coupled Cluster Framework.

The local coupled cluster method DLPNO-CCSD(T) allows calculations on systems containing hundreds of atoms to be performed while typically reproducing canonical CCSD(T) energies with chemical accuracy. In this work, we present a scheme for decomposing the DLPNO-CCSD(T) interaction energy between two molecules into physical meaningful contributions, providing a quantification of the most important components of the chemical interaction. The method, called Local Energy Decomposition (LED), is straightforward and requires negligible additional computing time. Both the Hartree-Fock and the correlation energy are decomposed into contributions from localized or pairs of localized occupied orbitals. Assigning these localized orbitals to fragments allows one to differentiate between intra- and intermolecular contributions to the interaction energy. Accordingly, the interaction energy can be decomposed into electronic promotion, electrostatic, exchange, dynamic charge polarization, and dispersion contributions. The LED scheme is applied to a number of test cases ranging from weakly, dispersively bound complexes to systems with strong ionic interactions. The dependence of the results on the one-particle basis set and various technical aspects, such as the localization scheme, are carefully studied in order to ensure that the results do not suffer from technical artifacts. A numerical comparison between the DLPNO-CCSD(T)/LED and the popular symmetry adapted perturbation theory (DFT-SAPT) is made, and the limitations of the proposed scheme are discussed.

[1]  Y. Mo,et al.  Energy decomposition analysis of intermolecular interactions using a block-localized wave function approach , 2000 .

[2]  Keiji Morokuma,et al.  Why do molecules interact? The origin of electron donor-acceptor complexes, hydrogen bonding and proton affinity , 1977 .

[3]  Wilfried Meyer,et al.  PNO-CI and CEPA studies of electron correlation effects , 1974 .

[4]  Donald G Truhlar,et al.  Density functional theory for transition metals and transition metal chemistry. , 2009, Physical chemistry chemical physics : PCCP.

[5]  W. Kutzelnigg,et al.  Computed potential hypersurface (including electron correlation) of the system Li+/H2 , 1973 .

[6]  Hans-Joachim Werner,et al.  Local treatment of electron correlation in coupled cluster theory , 1996 .

[7]  Donald G Truhlar,et al.  Comparative DFT study of van der Waals complexes: rare-gas dimers, alkaline-earth dimers, zinc dimer, and zinc-rare-gas dimers. , 2006, The journal of physical chemistry. A.

[8]  Georg Jansen,et al.  Intermolecular dispersion energies from time-dependent density functional theory , 2003 .

[9]  Jirí Cerný,et al.  Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs. , 2006, Physical chemistry chemical physics : PCCP.

[10]  Frederick R Manby,et al.  The orbital-specific virtual local triples correction: OSV-L(T). , 2013, The Journal of chemical physics.

[11]  Robert M Parrish,et al.  Tractability gains in symmetry-adapted perturbation theory including coupled double excitations: CCD+ST(CCD) dispersion with natural orbital truncations. , 2013, The Journal of chemical physics.

[12]  J. Herbert,et al.  Energy Decomposition Analysis with a Stable Charge-Transfer Term for Interpreting Intermolecular Interactions. , 2016, Journal of chemical theory and computation.

[13]  P Pulay,et al.  Local Treatment of Electron Correlation , 1993 .

[14]  Pavel Hobza,et al.  Highly accurate CCSD(T) and DFT-SAPT stabilization energies of H-bonded and stacked structures of the uracil dimer. , 2008, Chemphyschem : a European journal of chemical physics and physical chemistry.

[15]  Ł M Mentel,et al.  Can the Counterpoise Correction for Basis Set Superposition Effect Be Justified? , 2014, Journal of chemical theory and computation.

[16]  Donald G Truhlar,et al.  Density functionals with broad applicability in chemistry. , 2008, Accounts of chemical research.

[17]  Frank Neese,et al.  Natural triple excitations in local coupled cluster calculations with pair natural orbitals. , 2013, The Journal of chemical physics.

[18]  S. Grimme,et al.  Density functional theory including dispersion corrections for intermolecular interactions in a large benchmark set of biologically relevant molecules. , 2006, Physical chemistry chemical physics : PCCP.

[19]  H. Lischka,et al.  PNO–CI (pair natural orbital configuration interaction) and CEPA–PNO (coupled electron pair approximation with pair natural orbitals) calculations of molecular systems. II. The molecules BeH2, BH, BH3, CH4, CH−3, NH3 (planar and pyramidal), H2O, OH+3, HF and the Ne atom , 1975 .

[20]  Lori A Burns,et al.  Levels of symmetry adapted perturbation theory (SAPT). I. Efficiency and performance for interaction energies. , 2014, The Journal of chemical physics.

[21]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[22]  Peter Pulay,et al.  Fourth‐order Mo/ller–Plessett perturbation theory in the local correlation treatment. I. Method , 1987 .

[23]  Peter Pulay,et al.  Accurate correlated calculation of the intermolecular potential surface in the coronene dimer , 2010 .

[24]  Martin Schütz,et al.  Molpro: a general‐purpose quantum chemistry program package , 2012 .

[25]  Pavel Hobza,et al.  S66: A Well-balanced Database of Benchmark Interaction Energies Relevant to Biomolecular Structures , 2011, Journal of chemical theory and computation.

[26]  Robert M Parrish,et al.  Spatial assignment of symmetry adapted perturbation theory interaction energy components: The atomic SAPT partition. , 2014, The Journal of chemical physics.

[27]  J. Makarewicz,et al.  Does DFT-SAPT method provide spectroscopic accuracy? , 2015, The Journal of chemical physics.

[28]  V. Staemmler.  Ab initio calculations on small hydrides including electron correlation , 1973 .

[29]  Frank Weinhold,et al.  Natural bond orbital methods , 2012 .

[30]  Donald G Truhlar,et al.  Benchmark Databases for Nonbonded Interactions and Their Use To Test Density Functional Theory. , 2005, Journal of chemical theory and computation.

[31]  Robert Moszynski,et al.  Perturbation Theory Approach to Intermolecular Potential Energy Surfaces of van der Waals Complexes , 1994 .

[32]  Artur Michalak,et al.  A Combined Charge and Energy Decomposition Scheme for Bond Analysis. , 2009, Journal of chemical theory and computation.

[33]  Peter Pulay,et al.  CAN (SEMI) LOCAL DENSITY FUNCTIONAL THEORY ACCOUNT FOR THE LONDON DISPERSION FORCES , 1994 .

[34]  M. Head‐Gordon,et al.  An energy decomposition analysis for intermolecular interactions from an absolutely localized molecular orbital reference at the coupled-cluster singles and doubles level. , 2012, The Journal of chemical physics.

[35]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[36]  Pavel Hobza,et al.  Noncovalent interactions in biochemistry , 2011 .

[37]  Frank Neese,et al.  The ORCA program system , 2012 .

[38]  Georg Hetzer,et al.  Low-order scaling local electron correlation methods. I. Linear scaling local MP2 , 1999 .

[39]  Guntram Rauhut,et al.  Local Treatment of Electron Correlation in Molecular Clusters: Structures and Stabilities of (H2O)n, n = 2−4 , 1998 .

[40]  Peter Pulay,et al.  Local configuration interaction: An efficient approach for larger molecules , 1985 .

[41]  S. Grimme,et al.  Analysis of non-covalent interactions in (bio)organic molecules using orbital-partitioned localized MP2. , 2008, Physical chemistry chemical physics : PCCP.

[42]  F. Neese,et al.  Efficient and accurate local single reference correlation methods for high-spin open-shell molecules using pair natural orbitals. , 2011, The Journal of chemical physics.

[43]  M Pitoňák,et al.  Convergence of the CCSD(T) Correction Term for the Stacked Complex Methyl Adenine-Methyl Thymine: Comparison with Lower-Cost Alternatives. , 2009, Journal of chemical theory and computation.

[44]  L. Belpassi,et al.  Advances in Charge Displacement Analysis. , 2016, Journal of chemical theory and computation.

[45]  D. Truhlar,et al.  Benchmark calculations of reaction energies, barrier heights, and transition-state geometries for hydrogen abstraction from methanol by a hydrogen atom. , 2005, The journal of physical chemistry. A.

[46]  Hans-Joachim Werner,et al.  Low-order scaling local electron correlation methods. IV. Linear scaling local coupled-cluster (LCCSD) , 2001 .

[47]  C. David Sherrill,et al.  Wavefunction methods for noncovalent interactions , 2012 .

[48]  L. Slipchenko,et al.  Accurate Prediction of Noncovalent Interaction Energies with the Effective Fragment Potential Method: Comparison of Energy Components to Symmetry-Adapted Perturbation Theory for the S22 Test Set. , 2012, Journal of chemical theory and computation.

[49]  Stefan Grimme,et al.  Effect of the damping function in dispersion corrected density functional theory , 2011, J. Comput. Chem..

[50]  Georg Jansen,et al.  Symmetry‐adapted perturbation theory based on density functional theory for noncovalent interactions , 2014 .

[51]  Eric D. Glendening,et al.  Natural energy decomposition analysis: An energy partitioning procedure for molecular interactions with application to weak hydrogen bonding, strong ionic, and moderate donor–acceptor interactions , 1994 .

[52]  Donald G. Truhlar,et al.  Hybrid Meta Density Functional Theory Methods for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions: The MPW1B95 and MPWB1K Models and Comparative Assessments for Hydrogen Bonding and van der Waals Interactions , 2004 .

[53]  Rustam Z. Khaliullin,et al.  Unravelling the origin of intermolecular interactions using absolutely localized molecular orbitals. , 2007, The journal of physical chemistry. A.

[54]  S. Grimme,et al.  A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. , 2010, The Journal of chemical physics.

[55]  W. Kutzelnigg,et al.  Ab initio calculations on small hydrides including electron correlation , 1968 .

[56]  Gerald Knizia,et al.  Intrinsic Atomic Orbitals: An Unbiased Bridge between Quantum Theory and Chemical Concepts. , 2013, Journal of chemical theory and computation.

[57]  Dimitrios G Liakos,et al.  Is It Possible To Obtain Coupled Cluster Quality Energies at near Density Functional Theory Cost? Domain-Based Local Pair Natural Orbital Coupled Cluster vs Modern Density Functional Theory. , 2015, Journal of chemical theory and computation.

[58]  Pavel Hobza,et al.  MP2.5 and MP2.X: approaching CCSD(T) quality description of noncovalent interaction at the cost of a single CCSD iteration. , 2013, Chemphyschem : a European journal of chemical physics and physical chemistry.

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

[60]  H. Schaefer,et al.  Structures of protonated benzene dimer and intermolecular interaction decomposition via symmetry-adapted perturbation theory , 2011 .

[61]  A. Hesselmann.  Local correlation potentials from Brueckner coupled-cluster theory. , 2005, The Journal of chemical physics.

[62]  Peter Pulay,et al.  The local correlation treatment. II. Implementation and tests , 1988 .

[63]  C. Ochsenfeld,et al.  Linear-scaling symmetry-adapted perturbation theory with scaled dispersion. , 2013, The Journal of chemical physics.

[64]  José M. Pérez-Jordá,et al.  A density-functional study of van der Waals forces: rare gas diatomics. , 1995 .

[65]  Klaus Ruedenberg,et al.  The Physical Nature of the Chemical Bond , 1962 .

[66]  Frank Neese,et al.  Accurate theoretical chemistry with coupled pair models. , 2009, Accounts of chemical research.

[67]  Krzysztof Szalewicz,et al.  Symmetry‐adapted perturbation theory of intermolecular forces , 2012 .

[68]  Manoj K. Kesharwani,et al.  Exploring the Accuracy Limits of Local Pair Natural Orbital Coupled-Cluster Theory. , 2015, Journal of chemical theory and computation.

[69]  Dimitrios G Liakos,et al.  Efficient and accurate approximations to the local coupled cluster singles doubles method using a truncated pair natural orbital basis. , 2009, The Journal of chemical physics.

[70]  F. Neese,et al.  Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method. , 2009, The Journal of chemical physics.

[71]  F. B. van Duijneveldt,et al.  Convergence to the basis‐set limit in ab initio calculations at the correlated level on the water dimer , 1992 .

[72]  M. Head‐Gordon,et al.  An energy decomposition analysis for second-order Møller-Plesset perturbation theory based on absolutely localized molecular orbitals. , 2015, The Journal of chemical physics.

[73]  F. Neese,et al.  Accurate thermochemistry from a parameterized coupled-cluster singles and doubles model and a local pair natural orbital based implementation for applications to larger systems. , 2012, The Journal of chemical physics.

[74]  C. Sherrill,et al.  Assessment of Empirical Models versus High-Accuracy Ab Initio Methods for Nucleobase Stacking: Evaluating the Importance of Charge Penetration. , 2015, Journal of chemical theory and computation.

[75]  Hans-Joachim Werner,et al.  Local perturbative triples correction (T) with linear cost scaling , 2000 .

[76]  S. Grimme.  Density functional theory with London dispersion corrections , 2011 .

[77]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials , 1995 .

[78]  J. Piquemal,et al.  New intermolecular benchmark calculations on the water dimer: SAPT and supermolecular post‐Hartree–Fock approaches , 2009 .

[79]  S. F. Boys.  Construction of Some Molecular Orbitals to Be Approximately Invariant for Changes from One Molecule to Another , 1960 .

[80]  Dimitrios G Liakos,et al.  Weak Molecular Interactions Studied with Parallel Implementations of the Local Pair Natural Orbital Coupled Pair and Coupled Cluster Methods. , 2011, Journal of chemical theory and computation.

[81]  Frank Neese,et al.  An efficient and near linear scaling pair natural orbital based local coupled cluster method. , 2013, The Journal of chemical physics.

[82]  S. Grimme,et al.  Dispersion-corrected density functional theory for aromatic interactions in complex systems. , 2013, Accounts of chemical research.

[83]  Frank Neese,et al.  Sparse maps--A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory. , 2016, The Journal of chemical physics.

[84]  Paul G. Mezey,et al.  A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions , 1989 .