Application of the Electrostatically Embedded Many-Body Expansion to Microsolvation of Ammonia in Water Clusters.

The electrostatically embedded many-body expansion (EE-MB), at both the second and third order, that is, the electrostatically embedded pairwise additive (EE-PA) approximation and the electrostatically embedded three-body (EE-3B) approximation, are tested for mixed ammonia-water clusters. We examine tetramers, pentamers, and hexamers for three different density functionals and two levels of wave function theory, We compare the many-body results to the results of full calculations performed without many-body expansions. Because of the differing charge distributions in the two kinds of monomers, this provides a different kind of test of the usefulness of the EE-MB method than was provided by previous tests on pure water clusters. We find only small errors due to the truncation of the many-body expansion for the mixed clusters. In particular, for tests on tetramers and pentamers, the mean absolute deviations for truncation at second order are 0.36-0.98 kcal/mol (average: 0.66 kcal/mol), and the mean absolute deviations for truncation at third order are 0.04-0.28 (average: 0.16 kcal/mol). These may be compared to a spread of energies as large as 4.24 kcal/mol in the relative energies of various structures of pentamers and to deviations of up to 8.57 kcal/mol of the full calculations of relative energies from the best estimates of the relative energies. When the methods are tested on hexamers, the mean unsigned deviation per monomer remains below 0.10 kcal/mol for EE-PA and below 0.03 kcal/mol for EE-3B. Thus the additional error due to the truncation of the expansion is small compared to the accuracy needed or the other approximations involved in practical calculations. This means that the EE-MB expansion in combination with density functional theory or wave function theory for the oligomers provides a useful practical model chemistry for making electronic structure calculations and simulations more affordable by improving the scaling with respect to system size.

[1]  Martin Head-Gordon,et al.  Quantum chemistry and molecular processes , 1996 .

[2]  T. Dyke,et al.  Rotational spectra and structure of the ammonia–water complex , 1985 .

[3]  D. R. Hanson,et al.  First measurement of prenucleation molecular clusters , 2000 .

[4]  Yuichi Inadomi,et al.  Fragment molecular orbital method: application to molecular dynamics simulation, ‘ab initio FMO-MD’ , 2003 .

[5]  P. Mcmurry,et al.  Measurements of new particle formation and ultrafine particle growth rates at a clean continental site , 1997 .

[6]  M. Canagaratna,et al.  The Nitric Acid−Water Complex: Microwave Spectrum, Structure, and Tunneling , 1998 .

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

[8]  Roger E Bumgarner,et al.  Microwave and tunable far‐infrared laser spectroscopy of the ammonia–water dimer , 1992 .

[9]  P. Åstrand,et al.  Novel model for calculating the intermolecular part of the infrared spectrum for molecular complexes , 1995 .

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

[11]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[12]  Donald G Truhlar,et al.  Assessment of Model Chemistries for Noncovalent Interactions. , 2006, Journal of chemical theory and computation.

[13]  Donald G Truhlar,et al.  SM6:  A Density Functional Theory Continuum Solvation Model for Calculating Aqueous Solvation Free Energies of Neutrals, Ions, and Solute-Water Clusters. , 2005, Journal of chemical theory and computation.

[14]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[15]  J. Novoa,et al.  Ionization of Bases in Water: Structure and Stability of the NH4+···OH- Ionic Forms in Ammonia−Water Clusters , 1996 .

[16]  Kazuo Kitaura,et al.  Extending the power of quantum chemistry to large systems with the fragment molecular orbital method. , 2007, The journal of physical chemistry. A.

[17]  D. Donaldson.  Adsorption of Atmospheric Gases at the Air−Water Interface. I. NH3 , 1999 .

[18]  D. Truhlar,et al.  Electrostatically Embedded Many-Body Expansion for Simulations. , 2008, Journal of chemical theory and computation.

[19]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[20]  Donald G Truhlar,et al.  Electrostatically Embedded Many-Body Expansion for Large Systems, with Applications to Water Clusters. , 2007, Journal of chemical theory and computation.

[21]  John A. Pople,et al.  Nobel Lecture: Quantum chemical models , 1999 .

[22]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

[23]  P. Mcmurry,et al.  Measurement of Expected Nucleation Precursor Species and 3–500-nm Diameter Particles at Mauna Loa Observatory, Hawaii , 1995 .

[24]  Marat Valiev,et al.  Fast electron correlation methods for molecular clusters in the ground and excited states , 2005 .

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