Spatial assignment of symmetry adapted perturbation theory interaction energy components: The atomic SAPT partition.

We develop a physically-motivated assignment of symmetry adapted perturbation theory for intermolecular interactions (SAPT) into atom-pairwise contributions (the A-SAPT partition). The basic precept of A-SAPT is that the many-body interaction energy components are computed normally under the formalism of SAPT, following which a spatially-localized two-body quasiparticle interaction is extracted from the many-body interaction terms. For electrostatics and induction source terms, the relevant quasiparticles are atoms, which are obtained in this work through the iterative stockholder analysis (ISA) procedure. For the exchange, induction response, and dispersion terms, the relevant quasiparticles are local occupied orbitals, which are obtained in this work through the Pipek-Mezey procedure. The local orbital atomic charges obtained from ISA additionally allow the terms involving local orbitals to be assigned in an atom-pairwise manner. Further summation over the atoms of one or the other monomer allows for a chemically intuitive visualization of the contribution of each atom and interaction component to the overall noncovalent interaction strength. Herein, we present the intuitive development and mathematical form for A-SAPT applied in the SAPT0 approximation (the A-SAPT0 partition). We also provide an efficient series of algorithms for the computation of the A-SAPT0 partition with essentially the same computational cost as the corresponding SAPT0 decomposition. We probe the sensitivity of the A-SAPT0 partition to the ISA grid and convergence parameter, orbital localization metric, and induction coupling treatment, and recommend a set of practical choices which closes the definition of the A-SAPT0 partition. We demonstrate the utility and computational tractability of the A-SAPT0 partition in the context of side-on cation-π interactions and the intercalation of DNA by proflavine. A-SAPT0 clearly shows the key processes in these complicated noncovalent interactions, in systems with up to 220 atoms and 2845 basis functions.

[1]  R. Ahlrichs,et al.  Efficient molecular numerical integration schemes , 1995 .

[2]  Tatiana Korona,et al.  Exchange-Dispersion Energy: A Formulation in Terms of Monomer Properties and Coupled Cluster Treatment of Intramonomer Correlation. , 2009, Journal of chemical theory and computation.

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

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

[5]  R. Wheatley,et al.  Atomic charge densities generated using an iterative stockholder procedure. , 2009, The Journal of chemical physics.

[6]  Rick A. Kendall,et al.  The impact of the resolution of the identity approximate integral method on modern ab initio algorithm development , 1997 .

[7]  Florian Weigend,et al.  A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency , 2002 .

[8]  Alistair P. Rendell,et al.  COUPLED-CLUSTER THEORY EMPLOYING APPROXIMATE INTEGRALS : AN APPROACH TO AVOID THE INPUT/OUTPUT AND STORAGE BOTTLENECKS , 1994 .

[9]  A. van der Avoird,et al.  Symmetry‐adapted perturbation theory of nonadditive three‐body interactions in van der Waals molecules. I. General theory , 1995 .

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

[11]  Patrick Bultinck,et al.  Critical analysis and extension of the Hirshfeld atoms in molecules. , 2007, The Journal of chemical physics.

[12]  Jan Almlöf,et al.  Laplace transform techniques in Mo/ller–Plesset perturbation theory , 1992 .

[13]  R. Wheatley,et al.  Redefining the atom: atomic charge densities produced by an iterative stockholder approach. , 2008, Chemical communications.

[14]  J. Šponer,et al.  Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations? , 2010, The journal of physical chemistry. B.

[15]  D. Langreth,et al.  Stacking interactions and DNA intercalation. , 2009, The journal of physical chemistry. B.

[16]  Ashley L. Ringer,et al.  First principles computation of lattice energies of organic solids: the benzene crystal. , 2008, Chemistry.

[17]  Jean-Philip Piquemal,et al.  NCIPLOT: a program for plotting non-covalent interaction regions. , 2011, Journal of chemical theory and computation.

[18]  Andrew S. Dixon,et al.  Charge-Transfer Interactions in Macromolecular Systems: A New View of the Protein/Water Interface , 1998 .

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

[20]  J. Recio,et al.  Electron delocalization and bond formation under the ELF framework , 2011 .

[21]  A. Becke.  A multicenter numerical integration scheme for polyatomic molecules , 1988 .

[22]  G. Beran,et al.  Predicting Organic Crystal Lattice Energies with Chemical Accuracy , 2010 .

[23]  Crystal structure prediction for cyclotrimethylene trinitramine (RDX) from first principles. , 2009, Physical chemistry chemical physics : PCCP.

[24]  Martin W. Feyereisen,et al.  Use of approximate integrals in ab initio theory. An application in MP2 energy calculations , 1993 .

[25]  J. Almlöf,et al.  Integral approximations for LCAO-SCF calculations , 1993 .

[26]  Tatiana Korona,et al.  First-order exchange energy of intermolecular interactions from coupled cluster density matrices and their cumulants. , 2008, The Journal of chemical physics.

[27]  Kevin E. Riley,et al.  Pairwise decomposition of residue interaction energies using semiempirical quantum mechanical methods in studies of protein-ligand interaction. , 2005, Journal of the American Chemical Society.

[28]  M. Alderton,et al.  Distributed multipole analysis , 2006 .

[29]  C. Chabalowski,et al.  Using Kohn−Sham Orbitals in Symmetry-Adapted Perturbation Theory to Investigate Intermolecular Interactions , 2001 .

[30]  B. Rice,et al.  Predicting structure of molecular crystals from first principles. , 2008, Physical review letters.

[31]  H. Jónsson,et al.  Pipek-Mezey Orbital Localization Using Various Partial Charge Estimates. , 2014, Journal of chemical theory and computation.

[32]  J. Leszczynski,et al.  Physical nature of ethidium and proflavine interactions with nucleic acid bases in the intercalation plane. , 2006, The journal of physical chemistry. B.

[33]  J. Platts,et al.  Local electron correlation descriptions of the intermolecular stacking interactions between aromatic intercalators and nucleic acids , 2009 .

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

[35]  Tatiana Korona,et al.  Second-order exchange-induction energy of intermolecular interactions from coupled cluster density matrices and their cumulants. , 2008, Physical chemistry chemical physics : PCCP.

[36]  Kenneth M. Merz,et al.  The Role of Polarization and Charge Transfer in the Solvation of Biomolecules , 1999 .

[37]  Richard A. Friesner,et al.  Solution of the Hartree–Fock equations by a pseudospectral method: Application to diatomic molecules , 1986 .

[38]  Gilbert N. Lewis,et al.  The Atom and the Molecule , 1916, Resonance.

[39]  Richard A. Friesner,et al.  Solution of the Hartree–Fock equations for polyatomic molecules by a pseudospectral method , 1987 .

[40]  Pavel Hobza,et al.  Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations. , 2002, Journal of the American Chemical Society.

[41]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[42]  Frank Weinhold,et al.  Natural hybrid orbitals , 1980 .

[43]  D. Bondarev,et al.  Nature of Intercalator Amiloride−Nucelobase Stacking. An Empirical Potential and ab Initio Electron Correlation Study , 2000 .

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

[45]  Robert M Parrish,et al.  Tensor hypercontraction. II. Least-squares renormalization. , 2012, The Journal of chemical physics.

[46]  Tatiana Korona,et al.  A coupled cluster treatment of intramonomer electron correlation within symmetry-adapted perturbation theory: benchmark calculations and a comparison with a density-functional theory description , 2013 .

[47]  A. Stone,et al.  Distributed polarizabilities obtained using a constrained density-fitting algorithm. , 2006, The Journal of chemical physics.

[48]  Anthony J Stone,et al.  Distributed Multipole Analysis:  Stability for Large Basis Sets. , 2005, Journal of chemical theory and computation.

[49]  Krzysztof Szalewicz,et al.  Intermolecular potentials based on symmetry-adapted perturbation theory with dispersion energies from time-dependent density-functional calculations. , 2005, The Journal of chemical physics.

[50]  Jan Almlöf,et al.  Elimination of energy denominators in Møller—Plesset perturbation theory by a Laplace transform approach , 1991 .

[51]  Kazuo Kitaura,et al.  A new energy decomposition scheme for molecular interactions within the Hartree‐Fock approximation , 1976 .

[52]  R. Podeszwa,et al.  Three-body symmetry-adapted perturbation theory based on Kohn-Sham description of the monomers. , 2007, The Journal of chemical physics.

[53]  Robert M Parrish,et al.  Large-scale symmetry-adapted perturbation theory computations via density fitting and Laplace transformation techniques: investigating the fundamental forces of DNA-intercalator interactions. , 2011, The Journal of chemical physics.

[54]  J. L. Whitten,et al.  Coulombic potential energy integrals and approximations , 1973 .

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

[56]  K. Dill,et al.  Predicting absolute ligand binding free energies to a simple model site. , 2007, Journal of molecular biology.

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

[58]  C. David Sherrill,et al.  Potential energy curves for cation-pi interactions: off-axis configurations are also attractive. , 2009, The journal of physical chemistry. A.

[59]  A. Savin,et al.  Classification of chemical bonds based on topological analysis of electron localization functions , 1994, Nature.

[60]  Stanisl,et al.  Many‐body symmetry‐adapted perturbation theory of intermolecular interactions. H2O and HF dimers , 1991 .

[61]  C. Sherrill,et al.  Quantum-mechanical analysis of the energetic contributions to π stacking in nucleic acids versus rise, twist, and slide. , 2013, Journal of the American Chemical Society.

[62]  Edward G Hohenstein,et al.  Density fitting of intramonomer correlation effects in symmetry-adapted perturbation theory. , 2010, The Journal of chemical physics.

[63]  Charles L. Brooks,et al.  Are Many-Body Effects Important in Protein Folding? , 2000 .

[64]  Kenneth M Merz,et al.  Interaction energy decomposition in protein-protein association: a quantum mechanical study of barnase-barstar complex. , 2007, Biophysical chemistry.

[65]  T. Nakano,et al.  Fragment interaction analysis based on local MP2 , 2007 .

[66]  Edward G Hohenstein,et al.  Accurate Interaction Energies for Problematic Dispersion-Bound Complexes: Homogeneous Dimers of NCCN, P2, and PCCP. , 2011, Journal of chemical theory and computation.

[67]  Robert M Parrish,et al.  Exact tensor hypercontraction: a universal technique for the resolution of matrix elements of local finite-range N-body potentials in many-body quantum problems. , 2013, Physical review letters.

[68]  Richard F. W. Bader.  A quantum theory of molecular structure and its applications , 1991 .

[69]  Sequence-dependent configurational entropy change of DNA upon intercalation. , 2010, The journal of physical chemistry. B.

[70]  Matthew L. Leininger,et al.  Psi4: an open‐source ab initio electronic structure program , 2012 .

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

[72]  A. Stone,et al.  Localization methods for distributed polarizabilities , 1994 .

[73]  Huub J. J. Van Dam,et al.  Starting SCF calculations by superposition of atomic densities , 2006, J. Comput. Chem..

[74]  Paul S. Bagus,et al.  A new analysis of charge transfer and polarization for ligand–metal bonding: Model studies of Al4CO and Al4NH3 , 1984 .

[75]  K. Szalewicz,et al.  Dispersion Energy in the Coupled Pair Approximation with Noniterative Inclusion of Single and Triple Excitations , 1995 .

[76]  Y. Shao,et al.  Accelerating MP2C dispersion corrections for dimers and molecular crystals. , 2013, The Journal of chemical physics.

[77]  Robert M Parrish,et al.  Tensor hypercontraction density fitting. I. Quartic scaling second- and third-order Møller-Plesset perturbation theory. , 2012, The Journal of chemical physics.

[78]  Julia Contreras-García,et al.  Revealing noncovalent interactions. , 2010, Journal of the American Chemical Society.

[79]  Kenneth M. Merz,et al.  Divide and Conquer Interaction Energy Decomposition , 1999 .

[80]  J. Herbert,et al.  An efficient, fragment-based electronic structure method for molecular systems: self-consistent polarization with perturbative two-body exchange and dispersion. , 2011, The Journal of chemical physics.

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

[82]  Richard A. Friesner,et al.  Solution of self-consistent field electronic structure equations by a pseudospectral method , 1985 .

[83]  Dietrich Braess,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig Approximation of 1/x by Exponential Sums in [1, ∞) , 2022 .

[84]  John R. Sabin,et al.  On some approximations in applications of Xα theory , 1979 .

[85]  Donald G Truhlar,et al.  Perspectives on Basis Sets Beautiful: Seasonal Plantings of Diffuse Basis Functions. , 2011, Journal of chemical theory and computation.

[86]  C David Sherrill,et al.  Formal Estimation of Errors in Computed Absolute Interaction Energies of Protein-ligand Complexes. , 2011, Journal of chemical theory and computation.

[87]  Georg Jansen,et al.  The helium dimer potential from a combined density functional theory and symmetry-adapted perturbation theory approach using an exact exchange–correlation potential , 2003 .

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

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

[90]  Wolfgang Hackbusch,et al.  Minimax approximation for the decomposition of energy denominators in Laplace-transformed Møller-Plesset perturbation theories. , 2008, The Journal of chemical physics.

[91]  Anthony J. Stone,et al.  Distributed multipole analysis, or how to describe a molecular charge distribution , 1981 .

[92]  Sándor Suhai,et al.  SCC-DFTB-D study of intercalating carcinogens: Benzo(a)pyrene and its metabolites complexed with the G-C base pair , 2003 .

[93]  P. Jørgensen,et al.  Orbital localization using fourth central moment minimization. , 2012, The Journal of chemical physics.

[94]  K. Hongo,et al.  The Importance of Electron Correlation on Stacking Interaction of Adenine-Thymine Base-Pair Step in B-DNA: A Quantum Monte Carlo Study. , 2013, Journal of chemical theory and computation.

[95]  Wei Ji,et al.  Real-Space Identification of Intermolecular Bonding with Atomic Force Microscopy , 2013, Science.

[96]  R. Podeszwa,et al.  Efficient Calculations of Dispersion Energies for Nanoscale Systems from Coupled Density Response Functions. , 2012, Journal of chemical theory and computation.

[97]  F. J. Luque,et al.  Nature of base stacking: reference quantum-chemical stacking energies in ten unique B-DNA base-pair steps. , 2006, Chemistry.

[98]  Tatiana Korona,et al.  Dispersion energy from density-fitted density susceptibilities of singles and doubles coupled cluster theory. , 2008, The Journal of chemical physics.

[99]  P. Ayers,et al.  The conformational sensitivity of iterative stockholder partitioning schemes , 2012 .

[100]  Victor F. Lotrich,et al.  Symmetry-adapted perturbation theory of three-body nonadditivity of intermolecular interaction energy , 1997 .

[101]  P. Żuchowski,et al.  Symmetry-adapted perturbation theory based on unrestricted Kohn-Sham orbitals for high-spin open-shell van der Waals complexes. , 2012, The Journal of chemical physics.

[102]  K. Szalewicz,et al.  Third-order interactions in symmetry-adapted perturbation theory. , 2006, The Journal of chemical physics.

[103]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[104]  K. Szalewicz,et al.  Many‐body theory of exchange effects in intermolecular interactions. Second‐quantization approach and comparison with full configuration interaction results , 1994 .

[105]  John M Herbert,et al.  An improved treatment of empirical dispersion and a many-body energy decomposition scheme for the explicit polarization plus symmetry-adapted perturbation theory (XSAPT) method. , 2013, The Journal of chemical physics.

[106]  Axel D. Becke,et al.  A Simple Measure of Electron Localization in Atomic and Molecular-Systems , 1990 .

[107]  Covalent bond orders and atomic anisotropies from iterated stockholder atoms. , 2012, Physical chemistry chemical physics : PCCP.

[108]  Francesco Aquilante,et al.  Fast noniterative orbital localization for large molecules. , 2006, The Journal of chemical physics.

[109]  F. Weinhold,et al.  Natural population analysis , 1985 .

[110]  M. Schütz,et al.  Density-functional theory-symmetry-adapted intermolecular perturbation theory with density fitting: a new efficient method to study intermolecular interaction energies. , 2005, The Journal of chemical physics.

[111]  B. Silvi,et al.  What Can Tell the Topological Analysis of ELF on Hydrogen Bonding? , 2005 .

[112]  C. David Sherrill,et al.  Density fitting and Cholesky decomposition approximations in symmetry-adapted perturbation theory: Implementation and application to probe the nature of π-π interactions in linear acenes , 2010 .

[113]  Klaus Ruedenberg,et al.  Localized Atomic and Molecular Orbitals , 1963 .

[114]  E. Schröder,et al.  Stacking interactions and the twist of DNA. , 2008, Journal of the American Chemical Society.

[115]  J. Connolly,et al.  On the applicability of LCAO‐Xα methods to molecules containing transition metal atoms: The nickel atom and nickel hydride , 2009 .

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