Symmetry-adapted Perturbation Theory Calculations Research Articles

Transferable, Polarizable Force Field for Ionic Liquids.

A general, transferable, polarizable force field for molecular simulation of ionic liquids (ILs) and their mixtures with molecular compounds is developed. This polarizable model is derived from the widely used CL&P fixed-charge force field that describes most families of ILs, in a form compatible with OPLS-AA, one of the major force fields for organic compounds. Models for ILs with fixed, integer-ionic charges lead to pathologically slow dynamics, a problem that is corrected when polarization effects are included explicitly. In the model proposed here, Drude-induced dipoles are used with parameters determined from atomic polarizabilities. The CL&P force field is modified upon inclusion of the Drude dipoles to avoid double-counting of polarization effects. This modification is based on first-principles calculations of the dispersion and induction contributions to the van der Waals interactions using symmetry-adapted perturbation theory (SAPT) for a set of dimers composed of positive, negative, and neutral fragments representative of a wide variety of ILs. The fragment approach provides transferability, allowing the representation of a multitude of cation and anion families, including different functional groups, without the need to reparametrize. Because SAPT calculations are expensive, an alternative predictive scheme was devised, requiring only molecular properties with a clear physical meaning, namely, dipole moments and atomic polarizabilities. The new polarizable force field, CL&Pol, describes a broad set of ILs and their mixtures with molecular compounds and is validated by comparisons with experimental data on density, ion diffusion coefficients, and viscosity. The approaches proposed here can also be applied to the conversion of other fixed-charge force fields into polarizable versions.

Modeling Carbon Dioxide Vibrational Frequencies in Ionic Liquids: II. Spectroscopic Map.

The primary challenge for connecting molecular dynamics (MD) simulations to linear and two-dimensional infrared measurements is the calculation of the vibrational frequency for the chromophore of interest. Computing the vibrational frequency at each time step of the simulation with a quantum mechanical method like density functional theory (DFT) is generally prohibitively expensive. One approach to circumnavigate this problem is the use of spectroscopic maps. Spectroscopic maps are empirical relationships that correlate the frequency of interest to properties of the surrounding solvent that are readily accessible in the MD simulation. Here, we develop a spectroscopic map for the asymmetric stretch of CO2 in the 1-butyl-3-methylimidazolium hexafluorophosphate ([C4C1im][PF6]) ionic liquid (IL). DFT is used to compute the vibrational frequency of 500 statistically independent CO2-[C4C1im][PF6] clusters extracted from an MD simulation. When the map was tested on 500 different CO2-[C4C1im][PF6] clusters, the correlation coefficient between the benchmark frequencies and the predicted frequencies was R = 0.94, and the root-mean-square error was 2.7 cm-1. The calculated distribution of frequencies also agrees well with experiment. The spectroscopic map required information about the CO2 angle, the electrostatics of the surrounding solvent, and the Lennard-Jones interaction between the CO2 and the IL. The contribution of each term in the map was investigated using symmetry-adapted perturbation theory calculations.

Symmetry-Adapted Perturbation Theory Energy Analysis of Alkyl Fluorine-Aromatic Interactions in Torsion Balance Systems.

Symmetry-adapted perturbation theory (SAPT) calculations are carried out to elucidate the intermolecular interactions present between fluorinated and nonfluorinated alkyl chain groups and aromatic π systems in the folded and unfolded conformers of Wilcox torsion balance systems. The calculations predict the folded conformers to be 2.0-2.3 kcal/mol more stable than the unfolded conformers, with the preference for the folded conformer being greater in the fluorinated alkyl chain case. We also establish that a simple electrostatic analysis, based on atomic charges, is inadequate for understanding the conformational preferences of these systems. In the folded conformers, there are sizable charge penetration contributions that are not recovered by point charge models. Additionally, the SAPT analysis reveals that exchange-repulsion interactions make a significant contribution to the relative stability of the folded and unfolded conformer.

Rotational Spectrum, Structure, and Interaction Energy of the Trifluoroethylene···Carbon Dioxide Complex.

Rotational spectra for four isotopologues of the 1:1 weakly bound complex between trifluoroethylene (HFC═CF2) and carbon dioxide (CO2) were recorded using 480 MHz bandwidth chirped-pulse and resonant cavity Fourier transform microwave spectroscopy between 5.0 and 18.5 GHz. Two planar forms are possible: experimental rotational constants, planar moments, and dipole moment components are consistent with the form in which CO2 is positioned at the CHF end of the TFE subunit and is approximately perpendicular to the C═C bond; the other form, with CO2 aligned roughly parallel to the C═C bond, is not observed, consistent with ab initio relative energy predictions. Symmetry-adapted perturbation theory (SAPT) calculations provided interaction energies for possible structural forms of this complex, and comparisons are made with this and other members of the series of carbon dioxide complexes with fluorinated ethylenes (vinyl fluoride, 1,1-difluoroethylene, cis- and trans-1,2-difluoroethylene, and trifluoroethylene).

Theoretical investigation of the nature and strength of simultaneous interactions of π–π stacking and halogen bond including NMR, SAPT, AIM and NBO analysis

Ab initio MP2/aug-cc-pVDZ calculations were performed to investigate mutual effect between π–π stacking and halogen bond interactions in X-ben||pyr···Cl–F complexes (X = CN, F, Cl, Br, CH3, OH and H where || and ··· denote π–π stacking and halogen bonds). The results indicate the cooperativity of π–π stacking and halogen bonds in these complexes. This effect was discussed in terms of the energetic, geometrical parameters and charge-transfer properties of the complexes. To explore on the two-bonded spin–spin coupling constant 2X J(N–F) across 15N···35Cl–19F halogen bond in X-ben||pyr···Cl–F complexes, NMR calculations were performed at PBE0/aug-cc-pVDZ levels of theory. To get more insight into the physical nature of the binding energies, Symmetry Adapted Perturbation Theory calculations were carried out. Energy decomposition indicates that the percentage of the electrostatic term in the halogen bonding system constitutes approximately half of the total attractive binding energies, while the percentage of the dispersion term in the π–π stacking complexes constitutes approximately half of the attractive binding energies. In addition, atoms in molecules, natural bond orbital and molecular electrostatic potential were also used to probe the π–π stacking interactions and halogen bonding strengths.

First-Principles United Atom Force Field for the Ionic Liquid BMIM(+)BF4(-): An Alternative to Charge Scaling.

Molecular dynamics study of ionic liquids (ILs) is a challenging task. While accurate fully polarizable atomistic models exist, they are computationally too demanding for routine use. Most nonpolarizable atomistic models predict diffusion constants that are much lower than experiment. Scaled charge atomistic models are cost-effective and give good results for single component ILs but are in qualitative error for the phase behavior of mixtures, due to inaccurate prediction of the IL cohesive energy. In this work, we present an alternative approach for developing computationally efficient models that importantly preserves both the correct dynamics and cohesive energy of the IL. Employing a "top-down" approach, a hierarchy of coarse-grained models for BMIM(+)BF4(-) are developed by systematically varying the polarization/atomic resolution of the distinct functional groups. Parametrization is based on symmetry-adapted perturbation theory (SAPT) calculations involving the homomolecular species; all cross interactions are obtained from mixing rules, and there are no adjustable parameters. We find that enhanced dynamics from a united-atom description counteracts the effect of reduced polarization, enabling computationally efficient models that exhibit quantitative agreement with experiment for both static and dynamic properties. We give explicit suggestions for reduced-description models that are computationally more efficient, more accurate, and more fundamentally sound than existing nonpolarizable atomistic models.

Next-Generation Force Fields from Symmetry-Adapted Perturbation Theory.

Symmetry-adapted perturbation theory (SAPT) provides a unique set of advantages for parameterizing next-generation force fields from first principles. SAPT provides a direct, basis-set superposition error free estimate of molecular interaction energies, a physically intuitive energy decomposition, and a seamless transition to an asymptotic picture of intermolecular interactions. These properties have been exploited throughout the literature to develop next-generation force fields for a variety of applications, including classical molecular dynamics simulations, crystal structure prediction, and quantum dynamics/spectroscopy. This review provides a brief overview of the formalism and theory of SAPT, along with a practical discussion of the various methodologies utilized to parameterize force fields from SAPT calculations. It also highlights a number of applications of SAPT-based force fields for chemical systems of particular interest. Finally, the review ends with a brief outlook on the future opportunities and challenges that remain for next-generation force fields based on SAPT.

Theoretical Characterization of the Minimum-Energy Structure of (SF6)2.

MP2 and symmetry-adapted perturbation theory calculations are used in conjunction with the aug-cc-pVQZ basis set to characterize the SF6 dimer. Both theoretical methods predict the global minimum structure to be of C2 symmetry, lying 0.07-0.16 kJ/mol below a C2h saddle point structure, which, in turn, is predicted to lie energetically 0.4-0.5 kJ/mol below the lowest-energy D2d structure. This is in contrast with IR spectroscopic studies that infer an equilibrium D2d structure. It is proposed that the inclusion of vibrational zero-point motion gives an averaged structure of D2d symmetry.

Anion Recognition by Pyrylium Cations and Thio-, Seleno- and Telluro- Analogues: A Combined Theoretical and Cambridge Structural Database Study

Pyrylium salts are a very important class of organic molecules containing a trivalent oxygen atom in a six-membered aromatic ring. In this manuscript, we report a theoretical study of pyrylium salts and their thio-, seleno- and telluro- analogues by means of DFT calculations. For this purpose, unsubstituted 2,4,6-trimethyl and 2,4,6-triphenyl cations and anions with different morphologies were chosen (Cl–, NO3– and BF4–). The complexes were characterized by means of natural bond orbital and “atoms-in-molecules” theories, and the physical nature of the interactions has been analyzed by means of symmetry-adapted perturbation theory calculations. Our results indicate the presence of anion-π interactions and chalcogen bonds based on both σ- and π-hole interactions and the existence of very favorable σ-complexes, especially for unsubstituted cations. The electrostatic component is dominant in the interactions, although the induction contributions are important, particularly for chloride complexes. The geometrical features of the complexes have been compared with experimental data retrieved from the Cambridge Structural Database.

Accurate description of intermolecular interactions involving ions using symmetry-adapted perturbation theory.

Three new data sets for intermolecular interactions, AHB21 for anion-neutral dimers, CHB6 for cation-neutral dimers, and IL16 for ion pairs, are assembled here, with complete-basis CCSD(T) results for each. These benchmarks are then used to evaluate the accuracy of the single-exchange approximation that is used for exchange energies in symmetry-adapted perturbation theory (SAPT), and the accuracy of SAPT based on wave function and density-functional descriptions of the monomers is evaluated. High-level SAPT calculations afford poor results for these data sets, and this includes the recently proposed "gold", "silver", and "bronze standards" of SAPT, namely, SAPT2+(3)-δMP2/aug-cc-pVTZ, SAPT2+/aug-cc-pVDZ, and sSAPT0/jun-cc-pVDZ, respectively [ Parker , T. M. , et al. , J. Chem. Phys. 2014 , 140 , 094106 ]. Especially poor results are obtained for symmetric shared-proton systems of the form X(-)···H(+)···X(-), for X = F, Cl, or OH. For the anionic data set, the SAPT2+(CCD)-δMP2/aug-cc-pVTZ method exhibits the best performance, with a mean absolute error (MAE) of 0.3 kcal/mol and a maximum error of 0.7 kcal/mol. For the cationic data set, the highest-level SAPT method, SAPT2+3-δMP2/aug-cc-pVQZ, outperforms the rest of the SAPT methods, with a MAE of 0.2 kcal/mol and a maximum error of 0.4 kcal/mol. For the ion-pair data set, the SAPT2+3-δMP2/aug-cc-pVTZ performs the best among all SAPT methods with a MAE of 0.3 kcal/mol and a maximum error of 0.9 kcal/mol. Overall, SAPT2+3-δMP2/aug-cc-pVTZ affords a small and balanced MAE (<0.5 kcal/mol) for all three data sets, with an overall MAE of 0.4 kcal/mol. Despite the breakdown of perturbation theory for ionic systems at short-range, SAPT can still be saved given two corrections: a "δHF" correction, which requires a supermolecular Hartree-Fock calculation to incorporate polarization effects beyond second order, and a "δMP2" correction, which requires a supermolecular MP2 calculation to account for higher-order induction-dispersion coupling. These corrections serve to remove artifacts introduced by the single exchange approximation in the exchange-induction and exchange-dispersion interactions, and obviate the need for ad hoc scaling of the first- and second-order exchange energies. Finally, some density-functional and MP2-based electronic structure methods are assessed as well, and we find that the best density-functional method for computing binding energies in these data sets is B97M-V/aug-cc-pVTZ, which affords a MAE of 0.4 kcal/mol, whereas complete-basis MP2 affords an MAE of 0.3 kcal/mol.

Assessment of hydrophobic interactions and their contributions through the analysis of the methane dimer.

Hydrophobic Interactions (HIs) are important in many phenomena of molecular recognition in chemistry and biology. Still, the relevance of HIs is sometimes difficult to evaluate particularly in large systems and intramolecular interactions. We put forward a method to estimate the magnitude and the different contributions of a given HI of the C···C, H-C···H, and H···H type through (i) the analysis of the electron density in the intermolecular region for eleven relative orientations of the methane dimer and (ii) the subsequent decomposition of the corresponding interaction energy in physically significant contributions using Symmetry Adapted Perturbation Theory (SAPT). Strong correlations were found between the topological properties of ρ(r) calculated at intermolecular bond critical points and Eint(SAPT) plus its different contributions with the C···C distance of the considered orientations of (CH4 )2 . These correlations were used to construct Mollier-like diagrams of Eint(SAPT) and its components as a function of the separation between two carbons and the orientation of the groups bonded to these atoms. The ethane dimer and tert-butylcyclohexane are used as representative examples of this new approach. Overall, we anticipate that this new method might prove useful in the study of both intramolecular and intermolecular HIs particularly of those within large systems wherein SAPT or electronic structure calculations are computationally expensive or even prohibitive.

Determination of structure and properties of molecular crystals from first principles.

CONSPECTUS: Until recently, it had been impossible to predict structures of molecular crystals just from the knowledge of the chemical formula for the constituent molecule(s). A solution of this problem has been achieved using intermolecular force fields computed from first principles. These fields were developed by calculating interaction energies of molecular dimers and trimers using an ab initio method called symmetry-adapted perturbation theory (SAPT) based on density-functional theory (DFT) description of monomers [SAPT(DFT)]. For clusters containing up to a dozen or so atoms, interaction energies computed using SAPT(DFT) are comparable in accuracy to the results of the best wave function-based methods, whereas the former approach can be applied to systems an order of magnitude larger than the latter. In fact, for monomers with a couple dozen atoms, SAPT(DFT) is about equally time-consuming as the supermolecular DFT approach. To develop a force field, SAPT(DFT) calculations are performed for a large number of dimer and possibly also trimer configurations (grid points in intermolecular coordinates), and the interaction energies are then fitted by analytic functions. The resulting force fields can be used to determine crystal structures and properties by applying them in molecular packing, lattice energy minimization, and molecular dynamics calculations. In this way, some of the first successful determinations of crystal structures were achieved from first principles, with crystal densities and lattice parameters agreeing with experimental values to within about 1%. Crystal properties obtained using similar procedures but empirical force fields fitted to crystal data have typical errors of several percent due to low sensitivity of empirical fits to interactions beyond those of the nearest neighbors. The first-principles approach has additional advantages over the empirical approach for notional crystals and cocrystals since empirical force fields can only be extrapolated to such cases. As an alternative to applying SAPT(DFT) in crystal structure calculations, one can use supermolecular DFT interaction energies combined with scaled dispersion energies computed from simple atom-atom functions, that is, use the so-called DFT+D approach. Whereas the standard DFT methods fail for intermolecular interactions, DFT+D performs reasonably well since the dispersion correction is used not only to provide the missing dispersion contribution but also to fix other deficiencies of DFT. The latter cancellation of errors is unphysical and can be avoided by applying the so-called dispersionless density functional, dlDF. In this case, the dispersion energies are added without any scaling. The dlDF+D method is also one of the best performing DFT+D methods. The SAPT(DFT)-based approach has been applied so far only to crystals with rigid monomers. It can be extended to partly flexible monomers, that is, to monomers with only a few internal coordinates allowed to vary. However, the costs will increase relative to rigid monomer cases since the number of grid points increases exponentially with the number of dimensions. One way around this problem is to construct force fields with approximate couplings between inter- and intramonomer degrees of freedom. Another way is to calculate interaction energies (and possibly forces) "on the fly", i.e., in each step of lattice energy minimization procedure. Such an approach would be prohibitively expensive if it replaced analytic force fields at all stages of the crystal predictions procedure, but it can be used to optimize a few dozen candidate structures determined by other methods.

Theoretical study of the binding energy of a methane molecule in a (H2O)20 dodecahedral cage.

The interaction energy of a methane molecule encapsulated in a dodecahedral water cage is calculated using the MP2, MP2C, various dispersion-corrected DFT, and diffusion Monte Carlo (DMC) methods. The MP2, MP2C, and DMC methods give binding energies of -5.04, -4.60, and -5.3 ± 0.5 kcal/mol, respectively. In addition, the two- and three-body contributions are evaluated using the DFT, MP2, and CCSD(T) methods. All of the DFT methods considered appreciably overestimate the magnitude of the three-body contribution to the interaction energy. The two- and three-body energies are further analyzed by use of symmetry-adapted perturbation theory (SAPT) which allows decomposition into electrostatics, exchange, induction, and dispersion contributions. The SAPT calculations reveal that the induction, dispersion, and exchange three-body contributions to the methane-cage binding energy are all sizable, with the net three-body contribution to the binding energy being about 1 kcal/mol.

Symmetry-adapted perturbation theory with Kohn-Sham orbitals using non-empirically tuned, long-range-corrected density functionals.

The performance of second-order symmetry-adapted perturbation theory (SAPT) calculations using Kohn-Sham (KS) orbitals is evaluated against benchmark results for intermolecular interactions. Unlike previous studies of this "SAPT(KS)" methodology, the present study uses non-empirically tuned long-range corrected (LRC) functionals for the monomers. The proper v xc (r)→0 asymptotic limit is achieved by tuning the range separation parameter in order to satisfy the condition that the highest occupied KS energy level equals minus the molecule's ionization energy, for each monomer unit. Tests for He2, Ne2, and the S22 and S66 data sets reveal that this condition is important for accurate prediction of the non-dispersion components of the energy, although errors in SAPT(KS) dispersion energies remain unacceptably large. In conjunction with an empirical dispersion potential, however, the SAPT(KS) method affords good results for S22 and S66, and also accurately predicts the whole potential energy curve for the sandwich isomer of the benzene dimer. Tuned LRC functionals represent an attractive alternative to other asymptotic corrections that have been employed in density-functional-based SAPT calculations, and we recommend the use of tuned LRC functionals in both coupled-perturbed SAPT(DFT) calculations and dispersion-corrected SAPT(KS) calculations.

Energetic factors determining the binding of type I inhibitors to c-Met kinase: experimental studies and quantum mechanical calculations

To decipher the molecular interactions between c-Met and its type I inhibitors and to facilitate the design of novel c-Met inhibitors. Based on the prototype model inhibitor 1, four ligands with subtle differences in the fused aromatic rings were synthesized. Quantum chemistry was employed to calculate the binding free energy for each ligand. Symmetry-adapted perturbation theory (SAPT) was used to decompose the binding energy into several fundamental forces to elucidate the determinant factors. Binding free energies calculated from quantum chemistry were correlated well with experimental data. SAPT calculations showed that the predominant driving force for binding was derived from a sandwich π-π interaction with Tyr-1230. Arg-1208 was the differentiating factor, interacting with the 6-position of the fused aromatic ring system through the backbone carbonyl with a force pattern similar to hydrogen bonding. Therefore, a hydrogen atom must be attached at the 6-position, and changing the carbon atom to nitrogen caused unfavorable electrostatic interactions. The theoretical studies have elucidated the determinant factors involved in the binding of type I inhibitors to c-Met.

Predicting DNA–intercalator binding: the development of an arene–arene stacking parameter from SAPT analysis of benzene‐substituted benzene complexes

A series of substituted naphthalimides were synthesized and intercalated into the DNA sequence d(GCGCGCGC)2, and an experimental ΔTm value was obtained. Two‐parameter QSAR analyses were performed to generate a theoretical ΔTm value. Although by no means exhaustive in terms of parameter selection, the correlations did not yield statistics that indicated the models met the threshold for significance at the 95% confidence level. Rather than continue with an exhaustive search of all possible QSAR parameters, a one‐parameter QSAR analysis was performed utilizing a novel arene–arene stacking parameter, designated Ππ, developed from Symmetry‐Adapted Perturbation Theory (SAPT) energy decomposition studies of calculated benzene‐substituted benzene dimer binding energies. The QSAR analysis using the Ππ stacking parameter yielded statistics suggesting the model was significant at the 95% confidence level. The approach of developing a novel QSAR parameter via SAPT calculations, rather than exhaustively searching all traditional QSAR parameters, is presented both as a new approach for QSAR studies and as a unique application of SAPT. Copyright © 2013 John Wiley &amp; Sons, Ltd.

Adsorption of Small Organic Molecules on Graphene

We present a combined experimental and theoretical quantification of the adsorption enthalpies of seven organic molecules (acetone, acetonitrile, dichloromethane, ethanol, ethyl acetate, hexane, and toluene) on graphene. Adsorption enthalpies were measured by inverse gas chromatography and ranged from -5.9 kcal/mol for dichloromethane to -13.5 kcal/mol for toluene. The strength of interaction between graphene and the organic molecules was estimated by density functional theory (PBE, B97D, M06-2X, and optB88-vdW), wave function theory (MP2, SCS(MI)-MP2, MP2.5, MP2.X, and CCSD(T)), and empirical calculations (OPLS-AA) using two graphene models: coronene and infinite graphene. Symmetry-adapted perturbation theory calculations indicated that the interactions were governed by London dispersive forces (amounting to ∼60% of attractive interactions), even for the polar molecules. The results also showed that the adsorption enthalpies were largely controlled by the interaction energy. Adsorption enthalpies obtained from ab initio molecular dynamics employing non-local optB88-vdW functional were in excellent agreement with the experimental data, indicating that the functional can cover physical phenomena behind adsorption of organic molecules on graphene sufficiently well.

Ab initio intermolecular potential energy surfaces of He–CS2, Ne–CS2 and Ar–CS2 complexes

The potential energy surfaces of the He–CS2, Ne–CS2 and Ar–CS2 van der Waals complexes were calculated for the first time at the CCSD(T) level of theory using the aug-cc-pVDZ basis set augmented with a set of midbond functions (3s3p2d1f1g). It was found that the calculated interaction potential, using the applied basis set, readily converges to the complete basis set limit. For a broad range of intermolecular separations and configurations, the interaction energies were obtained by the supermolecular approach with the full counterpoise correction for the basis set superposition error (BSSE). In addition, symmetry-adapted perturbation theory (SAPT) calculations were performed with the same basis set in order to determine the character of the interaction energy of the most stable configuration of each complex at different intermolecular separations in order to make a comparison with the CCSD(T) results. The CCSD(T) calculated potential energy surface of each complex was fitted to an analytic expression to obtain the values of the isotropic dipole–dipole () and dipole–quadruple () dispersion coefficients of each complex. Finally, the interaction second virial coefficients (B12) were obtained using the calculated potential energy surface and used together with the experimental second virial coefficients of pure gases (CS2, Ar, Ne and He) to obtain the second virial coefficient of mixtures of CS2 with rare gas at different temperatures and mole fractions.

Breakdown of the Single-Exchange Approximation in Third-Order Symmetry-Adapted Perturbation Theory

We report third-order symmetry-adapted perturbation theory (SAPT) calculations for several dimers whose intermolecular interactions are dominated by induction. We demonstrate that the single-exchange approximation (SEA) employed to derive the third-order exchange-induction correction (E(exch-ind)((30))) fails to quench the attractive nature of the third-order induction (E(ind)((30))), leading to one-dimensional potential curves that become attractive rather than repulsive at short intermolecular separations. A scaling equation for (E(exch-ind)((30))), based on an exact formula for the first-order exchange correction, is introduced to approximate exchange effects beyond the SEA, and qualitatively correct potential energy curves that include third-order induction are thereby obtained. For induction-dominated systems, our results indicate that a "hybrid" SAPT approach, in which a dimer Hartree-Fock calculation is performed in order to obtain a correction for higher-order induction, is necessary not only to obtain quantitative binding energies but also to obtain qualitatively correct potential energy surfaces. These results underscore the need to develop higher-order exchange-induction formulas that go beyond the SEA.

Adsorption of a water molecule on the MgO(100) surface as described by cluster and slab models

The interaction of a water molecule with the (100) surface of MgO as described by cluster models is studied using MP2, coupled MP2 (MP2C) and symmetry-adapted perturbation theory (SAPT) methods. In addition, diffusion Monte Carlo (DMC) results are presented for several slab models as well as for the smallest, 2X2 cluster model. For the 2X2 model it is found that the MP2C, DMC, and CCSD(T) methods give nearly the same potential energy curve for the water-cluster interaction, whereas the potential energy curve from the SAPT calculations differs slightly from those of the other methods. The interaction of the water molecule with the cluster models of the MgO(100) surface is weakened upon expanding the number of layers from one to two and also upon expanding the description of the layers from 2X2 to 4X4 to 6X6. The SAPT calculations reveal that both these expansions of the cluster model are accompanied by reductions in the magnitudes of the induction and dispersion constributions. The best estimate of the energy for binding an isolated water molecule to the surface obtained from the cluster model calculations is in good agreement with that obtained from the DMC calculations using a 2-layer slab model with periodic boundary conditions.