Many-body interaction analysis: algorithm development and application to large molecular clusters.

A completely automated algorithm for performing many-body interaction energy analysis of clusters (MBAC) [M. J. Elrodt and R. J. Saykally, Chem. Rev. 94, 1975 (1994); S. S. Xantheas, J. Chem. Phys. 104, 8821 (1996)] at restricted Hartree-Fock (RHF)/MA Plesset 2nd order perturbation theory (MP2)/density functional theory (DFT) level of theory is reported. Use of superior guess density matrices (DM's) for smaller fragments generated from DM of the parent system and elimination of energetically insignificant higher-body combinations, leads to a more efficient performance (speed-up up to 2) compared to the conventional procedure. MBAC approach has been tested out on several large-sized weakly bound molecular clusters such as (H(2)O)(n), n=8, 12, 16, 20 and hydrated clusters of amides and aldehydes. The MBAC results indicate that the amides interact more strongly with water than aldehydes in these clusters. It also reconfirms minimization of the basis set superposition error for large cluster on using superior quality basis set. In case of larger weakly bound clusters, the contributions higher than four body are found to be repulsive in nature and smaller in magnitude. The reason for this may be attributed to the increased random orientations of the interacting molecules separated from each other by large distances.

[1]  Sotiris S. Xantheas,et al.  Ab initio studies of cyclic water clusters (H2O)n, n=1–6. II. Analysis of many‐body interactions , 1994 .

[2]  Sudhir A. Kulkarni,et al.  Ab initio investigations on neutral hydrogen peroxide clusters: (H2O2)n (n=2–4) , 2003 .

[3]  R. Mcweeny.  Perturbation Theory for the Fock-Dirac Density Matrix , 1962 .

[4]  Jules W. Moskowitz,et al.  Water Molecule Interactions , 1970 .

[5]  P. Rosen.  The Nonadditivity of the Repulsive Potential of Helium , 1953 .

[6]  L. Jansen,et al.  Stability of Crystals of Rare-Gas Atoms and Alkali Halides in Terms of Three-Body Interactions. II. Alkali-Halide Crystals , 1964 .

[7]  István Mayer,et al.  Hierarchy of counterpoise corrections for N-body clusters: generalization of the Boys-Bernardi scheme , 1997 .

[8]  A. Shostak.  Interaction Energy among Three Helium Atoms , 1955 .

[9]  M J Elrod,et al.  Many-body effects in intermolecular forces. , 1994, Chemical reviews.

[10]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[11]  Edward Teller,et al.  Interaction of the van der Waals Type Between Three Atoms , 1943 .

[12]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[13]  S. Xantheas,et al.  First Principles Examination of the Acetylene-Water Clusters, HCCH-(H2O)x, x ) 2, 3, and 4 , 2002 .

[14]  Sotiris S. Xantheas,et al.  On the importance of the fragment relaxation energy terms in the estimation of the basis set superposition error correction to the intermolecular interaction energy , 1996 .

[15]  K. Mierzwicki,et al.  Basis set superposition error in N-body clusters , 2003 .

[16]  E. Davidson,et al.  An analysis of the hydrogen bond in ice , 1990 .

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

[18]  Granino A. Korn,et al.  Mathematical handbook for scientists and engineers , 1961 .

[19]  Shridhar R. Gadre,et al.  Tailoring approach for exploring electron densities and electrostatic potentials of molecular crystals , 2004 .

[20]  L. Jansen,et al.  Perturbation Analysis of Three- and Four-Atom Exchange Interactions in a Gaussian Effective-Electron Model , 1968 .

[21]  K. Babu,et al.  Exploring hydration patterns of aldehydes and amides: ab initio investigations , 2004 .

[22]  Hideaki Umeyama,et al.  The origin of hydrogen bonding. An energy decomposition study , 1977 .

[23]  Shridhar R. Gadre,et al.  Structure and Stability of Water Clusters (H2O)n, n ) 8-20: An Ab Initio Investigation , 2001 .