Abstract
Background: The quantumchemical description of the interactions in water clusters is an essential basis for deriving accurate and physically sound models of the interaction potential for water to be used in molecular simulations. In particular, the role of manybody interactions beyond the twobody interactions, which are often not explicitly taken into account by empirical force fields, can be accurately described by quantum chemistry methods on an adequate level, e.g., randomphase approximation electron correlation methods. The relative magnitudes of the different interaction energy contributions obtained by accurate ab initio calculations can therefore provide useful insights that can be exploited to develop enhanced force field methods.
Results: In line with earlier theoretical studies of the interactions in water clusters, it has been found that the main contribution to the manybody interactions in clusters with a size of up to N = 13 molecules are higherorder polarisation interaction terms. Compared to this, manybody dispersion interactions are practically negligible for all studied sytems. The twobody dispersion interaction, however, plays a significant role in the formation of the structures of the water clusters and their stability, since it leads to a distinct compression of the cluster sizes compared to the structures optimized on an uncorrelated level. Overall, the manybody interactions amount to about 13% of the total interaction energy, irrespective of the cluster size. The electron correlation contribution to these, however, amounts to only about 30% to the total manybody interactions for the largest clusters studied and is repulsive for all structures considered in this work.
Conclusion: While this shows that three and higherbody interactions can not be neglected in the description of water complexes, the electron correlation contributions to these are much smaller in comparison to the twobody electron correlation effects. Efficient quantum chemistry approaches for describing intermolecular interactions between water molecules may therefore describe higherbody interactions on an uncorrelated Hartree–Fock level without a serious loss in accuracy.
Introduction
The description of the intermolecular interactions between water molecules is essential for an understanding of the structures and properties of water through the different stages of assemblies, from the dimer over the liquid phase to the bulk phase. Moreover, many chemical processes are explicitly or implicitly influenced by a water environment. An example for this is the hydrogenbond cooperativity effect that can have a significant impact on the properties of the bare solute molecules [1]. In order to describe such phenomena, computer simulations have become an indispensable tool, since they enable a description of water on a molecular level that often can provide further insights than are accessible from spectroscopic measurements.
The basis for such simulations are the socalled force fields that describe both the covalent as well as the noncovalent interactions within the system. These commonly depend on a number of empirical parameters that are determined either by a fit to experimentally known liquid or bulk properties, or by fitting to energies from ab initio quantum chemistry methods. The most popular potentials for water are the TIP3P [2], TIP4P [2,3] and TIP5P [4,5] force fields, which are based on a modeling of the water pair potential using an electrostatic contribution described by interacting point charges and a van der Waals interaction contribution using LennardJones potentials. In more advanced ab initio water pair potentials the force field is fitted to highlevel quantum chemistry results for the water dimer. Force fields belonging to this category are, e.g., the TTM3F and TTM4F models [6,7], the AMOEBA force field [8,9], the DPP2 model [10] and various force fields derived by Szalewicz and coworkers [1113]. There also exist a number of pair potentials that go beyond the point charge approximation [14].
A comparison with highlevel coupledcluster energies for a large number of water dimers and tetramers has revealed, however, that polarisation effects, which are not accounted for in the classical point charge potentials, are essential to describe the structures of (H_{2}O)_{2} and (H_{2}O)_{3} in many different conformations [15]. A rather good correlation between the coupledcluster energies and the forcefield energies is found for the polarisable AMOEBA2003 [9,16] and TTM4F [17] potentials both for the dimer and for the trimer. Both these methods are based on the induced dipole scheme in which polarisable point dipoles, which are assigned to the molecules, interact with the surrounding electric field and are computed in a selfconsistent manner. While these force fields, too, rely on a certain degree of empiricism, a number of other force field exist that aim at a more physically sound decomposition of the interaction energy into distinct contributions. Examples for such force fields are the sum of interaction between fragments (SIBFA) [1820] and the effective fragment potential (EFP) [21] method. The most recent version of the latter, EFP2, can describe both (longrange) polarisation as well as chargetransfer interactions. The latter was found to yield a significant contribution to the interaction energy of the water dimer [22].
However, one of the most significant results of [22] was that the dipole–quadrupole polarisability term of the multipole expansion of the dispersion interaction, which is usually neglected in force fields or dispersion corrected DFT methods, is a quite large positive (for the clusters considered) and anisotropic contribution to the interaction energy of small water clusters [22]. It was found to be almost half of the magnitude of the leadingorder dipole–dipole term. It was therefore concluded by Guidez et al. that this term should not be neglected in the description of the interactions in water.
Another challenge for water models is the description of nonadditive manybody terms to the interaction energy [2325]. It has been found that these contribute 15% to the total interactions in the condensed phase [26] and this amount even increases to 17–30% for small water clusters [2628]. Explicit evaluations of threebody interaction energy terms using symmetryadapted perturbation theory (SAPT) for the water trimer have revealed that the strongest contribution to the threebody energy originates from the polarisation (induction) energy while the threebody dispersion interaction is rather small [28,29]. Moreover, manybody exchange effects, including exchange–induction and exchange–dispersion, are relatively large yet cancel each other due to opposite signs at the uud minimum configuration [29] (“uud” indicates that two “free” hydrogen atoms point above the plane formed by the three oxygen atoms (u→up) and one below it (d→down) [30]). In the study by Hodges et al. it has been shown for several structures of the water tetramer that total fourbody interactions are in most cases much smaller or even negligible compared to the threebody interactions [27]. The only exception to this was observed for the squared geometry in which the hydrogen bonds act cooperatively to enhance the induction energy. In a recent work by Hapka et al. it was shown that also standard density functional theory methods are able to describe nonadditive effects to the interaction energy quite well for hydrogenbonded clusters, yet, fail to do so for dispersion bound complexes [31]. The common conclusion from the quantum chemistry studies of small water clusters was that damped classical polarization models should be able to accurately capture the nonadditive effects for larger clusters of water, because of the fact that shortrange contributions, including manybody exchange effects, grow less strongly with the size of the system than induction, dispersion or electrostatic interaction energies [32].
In this work we will study the impact of electron correlation effects and manybody interactions on the structures and energies of water clusters (H_{2}O)_{n} with cluster sizes ranging between n = 2 and n = 13. In doing so, we focus on the importance of dispersion energy contributions to the interaction energy, including twobody and manybody dispersion effects. It will be shown that overall electron correlation effects have only a minor impact on the orientation of the water molecules in the various minimum structures of the clusters. Dispersion interactions, however, make up a significant contribution of about 60% to the total interaction energies in larger water clusters and lead to a compression of the cluster sizes on average. The total magnitude of the electron correlation contribution to the interaction energy is, however, only about half the size of the sum of the two, three and fourbody dispersion interactions, i.e., the large dispersion interaction contribution is strongly quenched by further repulsive correlation contributions. Noting that the correlation effect to the molecular dipole moment of the water molecule reduces the dipole moment from the Hartree–Fock method by about 0.13 Debye [33], it can be assessed that a fraction of this repulsive correlation contribution originates from the reduction of the electrostatic interaction energy [34]. Empirical models for water that are based on a fitting to ab initio results therefore have to take the influences of the different correlation effects carefully into account.
ManyBody Expansion of the Interaction Energy
Consider a cluster system containing N molecules. The total interaction energy of this system is given by
with E(123…N) denoting the total energy of the system and E(A) (with A = 1,2,…) denoting the onebody clusters (monomers). The idea of the manybody expansion is to decompose the total interaction energy into terms arising from the twobody, threebody, fourbody, etc. interactions, so that the total energy is given by
with Δ^{M} denoting the Mth body interaction contribution. These describe the change of the total energy of the system due to the interactions of the Mth body clusters. The first three terms are defined as
As shown in Figure 1, the number of three and fourbody clusters increases much more rapidly with the size of the system than the number of twobody clusters. In spite of this, however, a truncation of the expansion of Equation 2 after the twobody term will usually capture 90% and more of the total interaction energy. This originates in part from the fact that for larger manybody systems many threebody and fourbody clusters possess structures with far distant molecules. Moreover, the individual interaction energy contributions to the three and higherbody interaction energy tend to be much smaller than the sum of the twobody interaction energies, as will be shown in this work.
Method
The calculations in this work have been performed using the EXXRPA (exactexchange randomphase approximation) electron correlation method [3537], which is based on an exactexchange Kohn–Sham reference determinant [3843]. The total energy of the EXXRPA method is given by
where E^{EXX} denotes the energy of the EXX reference determinant and is the EXXRPA correlation energy. Like the Hartree–Fock method, the EXX method lacks the description of electron correlation effects, i.e., electron–electron interactions beyond zeroth and first order in the Coulomb interaction. Thus, any inter and intramolecular correlation effects of the EXXRPA method are described by the correlation term of Equation 4.
In order to analyze individual interaction energy contributions of the twobody clusters, the DFTSAPT method [4452] method has been used. With this method, the interaction energy of a dimer system AB is given by
The interaction energy terms in Equation 5 are: electrostatic interaction energy, firstorder exchange interaction energy, induction energy, exchangeinduction energy, (twobody) dispersion energy and exchangedispersion energy. The exchange interaction energy terms in Equation 5 are shortrange contributions to the interaction energy and stem from a tunneling of the electrons between the two monomers. They quickly decay exponentially with the distance of the interacting systems. The superscripts (1) and (2) denote the order of the individual terms with respect to the interaction energy operator.
In this work all terms in Equation 5 have been computed using EXX monomer wave functions. Moreover, a timedependent EXX (TDEXX) response approach was used to compute the secondorder interaction energy contributions [5359]. Due to this choice, the subtotal
approximates the EXX interaction energy between the monomers A and B. The difference between Δ^{2}E^{EXX}(AB) and these interaction energy terms can be interpreted as higherorder exchange–induction interaction terms not accounted for by the DFTSAPT method (when truncated at second order).
Moreover, the sum of the dispersion and exchange–dispersion energy
is a fraction of the EXXRPA correlation energy contribution to the intermolecular interaction energy, since both the DFTSAPT terms as well as the EXXRPA correlation energy are computed with the exactexchange response kernel [53,54]. Here, the main difference of the quantities on the lefthand and righthand side of Equation 7 stems from the correlation energy contributions to the electrostatic, induction and their exchange interaction energy counterparts. Note that an analogous decomposition of the supermolecular interaction energy into distinct terms is also possible for the secondorder MøllerPlesset perturbation theory method [60].
Finally, we have also computed threebody and fourbody dispersion interaction energies that contribute to the threebody and fourbody interaction energy terms:
Thus, it will be possible to evaluate the importance of dispersion interactions both for the total manybody interaction energy as well as for its correlation interaction contribution.
Computational Details
The structures of the water clusters have been taken from the work of Maheshwary et al., see Figure 2. This set of structures contains both the prism form as well as the cage form of the water hexamer. These are almost isoenergetic and form the first noncyclic global minimum structures of (H_{2}O)_{x}. In addition, the cyclicchair conformation of the hexamer is considered, too, in this work, since it is known that this structure is characterised by strong manybody interactions [6163]. For consistency, the cyclicchair geometry has been optimized on the same level that was used in [23], namely Hartree–Fock employing the 631G(d,p) basis set [64]. The resulting structure is shown in Figure 3.
All calculations of this work were performed using the augccpVTZ basis set by Dunning and coworkers [65]. In order to correct the basis set error of the correlation energy terms, these were extrapolated using a doubleζ to tripleζ twopoint extrapolation using the formula from Bak and coworkers [66]. Core electrons have been kept frozen in these calculations.
The counterpoise correction of Boys and Bernadi [67] has been employed in all subcluster calculations to reduce the basis set superposition error. All calculations have been performed using a developers version of the Molpro quantum chemistry package [68,69].
Results and Discussion
Influence of correlation effects on the structures
To investigate the influence of electron correlation effects (including intermolecular dispersion interactions) on the structure of the water clusters, we have reoptimized the structures from [23] on the Hartree–Fock (HF) level using the def2TZVP basis set [70]. These modified HF geometries were then further optimized with the MP2 (secondorder MøllerPlesset perturbation theory) method, which takes electron correlation effects into account at the second order. The resulting structures for the dimer up to the heptamer are shown in Figure 4. Here, for better visibility, the HF structures are colored in blue and the MP2 structures in red. They are superpositioned in such a way that the average distances of the atoms of the respective geometries are minimised.
As can be seen in Figure 4, the electron correlation effects described by the MP2 method hardly change the global structures of the various clusters. This indicates that they mainly originate from electrostatic and induction interactions that can already be described reasonably well with the HF method. However, except for the trimer, where HF and MP2 geometries hardly differ from each other, one can observe that the correlation effects lead to a compression of the structures relative to the ones obtained with HF. This can be attributed to dispersion interactions between the molecules which is an additional attractive interaction energy contribution not accounted for by the HF method, see also below.
The dependence of the electrostatic and dispersion energy on the structure is highlighted in Figure 5 for three different conformations of the water dimer. These three structures have in common the distance of the oxygen atoms (2.98 Å in this example), but they possess different orientations of the hydrogen atoms. The first structure in Figure 5 corresponds to the equilibrium. As can be seen in the figure, the electrostatic interaction energy is strongly influenced by the orientation of the molecules and changes by almost +10 kcal/mol from the hydrogenbonded structure to the second one that is characterised by parallel dipoles of the two water molecules. Compared to this, the dispersion interaction hardly changes upon a disordering of the hydrogen atoms. While it has a minimum, too, at the equilibrium structure, for the other two structures it lies only +0.7 kcal/mol higher in energy. In line with the structure changes displayed in Figure 4, one can thus conclude that dispersion interactions almost act isotropically and therefore will generally try to maximise the contacts of the interacting sites.
This is also illustrated in Figure 6, which shows the twobody and threebody dispersion energies for three different structures of the water trimer. Here, again structure 1 corresponds to the equilibrium, which is characterised by an equilateral triangle formed by the three oxygen atoms. Compared to this, in structure 2 one of the sides of the triangle is extended, fixing however the other two at a length of 2.87 Å as in the first one. Finally, structure 3 is a linear shaped structure, see Figure 6. As can be seen in the diagram, when the trimer transforms from the most compact equilibrium structure to the linear form, the twobody dispersion interaction strongly reduces by +6 kcal/mol, yet remains attractive having a magnitude of about −6 kcal/mol for the structures 2 and 3. The blue horizontal bars in Figure 6 show the energy levels of the corresponding threebody dispersion energy for the three conformations. As can be seen, it possesses just the opposite dependency on the structure as the twobody dispersion energy. That is, it is repulsive at the triangularly shaped equilibrium structure but turns into an attractive contribution when the structure changes to the linearly shaped form (see the scale on the righthand side of the diagram in the figure). One can readily describe this anisotropic behavior of the threebody dispersion energy by the simplified Axilrod–Teller form of the interaction energy between three atoms [71]:
where R_{ij} and θ_{i} denote the sides and the angles of the ABC triangle and is a constant coefficient. Equation 10 shows that when ABC is in a linear configuration, the threebody dispersion energy is negative (an attractive contribution), while the equilateral triangular configuration leads to a repulsive interaction. Figure 7 shows the number of threebody subclusters that possess a stabilising threebody dispersion energy contribution (determined from the results of the calculations performed in this work). It can be seen that this number grows much less strongly with the cluster size than the number of subclusters with a destabilising threebody dispersion interaction. This qualitatively explains that the total threebody dispersion energy for the respective water clusters considered in this work is always repulsive, too (see below). Furthermore, a comparison of the total magnitudes of the twobody and threebody dispersion energies shows that the threebody dispersion interaction is much weaker than the twobody dispersion interaction. Therefore, its effect on the shape and energies of the larger water clusters can almost be neglected, see also below.
Manybody interactions in the water trimer and tetramer structures
In this section we analyse the three and fourbody interactions in the cyclic water trimer and tetramer structures, see Figure 2. More precisely, we here want to identify the main interaction energy contributions to the manybody interaction in these two cases. Note, though, that the interaction energy terms of the DFTSAPT method, see Equation (Equation 5), can only describe the twobody interactions between two subsystems A and B. While threebody contributions to the DFTSAPT method have been developed by Podeszwa and Szalewicz [29] manybody effects can, however, be also described by the twobody DFTSAPT terms with the aid of the pseudodimer technique. For this, recall that the interaction energy of a trimer ABC can be approximated by the sum of all mutual twobody interactions:
The twobody interaction terms in Equation 11 describe the interactions between two isolated monomers neglecting, however, the perturbation by the third one. Alternatively, twobody interactions can be determined by combining two monomers to one single (pseudo)monomer and calculating the interaction with the remaining one. In case of the trimer three different possibilities exist:
where the term Δ^{2}(AB − C) now denotes the interaction of the combined system AB with monomer C. As can be easily understood, the term Δ^{2}(AB − C) contains the two twobody interactions Δ^{2}(AC) and Δ^{2}(BC) and a remainder that describes the change of these two interactions due to the perturbation of A by B and vice versa. This precisely is the contribution that is not described by Equation 11 and, thus, is a contribution to the threebody interaction energy Δ^{3}(ABC). According to the three different possible pseudodimers, see Equation 12, one can extract the following three terms:
Recall that using the supermolecule method the term Δ^{3}(AB − C) is given by Δ^{3}(AB − C) = E(ABC) − E(AB) − E(C). An insertion of the corresponding total energy expressions into the other terms in Equation 13–Equation 15 shows that
and therefore the sum of the terms Δ^{3}(AB − C), Δ^{3}(AC − B) and Δ^{3}(BC − A) scaled by a factor of one third can be identified as threebody interaction energy.
We have calculated these terms using the DFTSAPT method, which allows us to analyse the contribution of the different interaction energy terms of Equation 5 to the threebody interaction energy. The results for the water trimer are shown in Figure 8. In addition, the diagram also contains the threebody interaction energy of the EXXRPA method as well as its correlation contribution. We have also compared the sum of the terms (the overlines are used here to distinguish the terms from the twobody SAPT interaction terms) to the total EXX threebody interaction energy and found a good agreement (DFTSAPT: −1.63 kcal/mol, EXX: −1.70 kcal/mol). This shows that higherorder interaction contributions to the threebody interaction energy of the EXX method are small.
The diagram in Figure 8 shows that the main contributions to the threebody interaction energy stem from the electrostatic and induction energies. These should not be confused with the standard twobody interaction energy counterparts but should rather be interpreted as the changes of the electrostatic and induction interactions due to the perturbations by another monomer in the trimer. Accordingly, one may interpret the electrostatic contribution in Figure 8 as an induction effect and the induction contribution as a higherorder polarisation effect contained in the threebody interaction energy. The exchange energy counterparts to and are quite small and strongly cancel each other, making up only a marginal net contribution to the threebody energy. This holds true also for the correlation contribution to the threebody interaction, including the dispersion energy, which is found to be repulsive. Therefore, in line with the findings of the previous section, the manybody correlation effects are negligible in the sum of the different threebody terms, see also Figure 6.
The pseudodimer scheme can also be applied to the tetramer. While this can be done in various ways, in this work we used the following terms:
One finds
and therefore, analogously to the procedure described above for the trimer, the sum has to be scaled by a factor of one third to reproduce the sum of the threebody and fourbody interactions. The DFTSAPT interaction energy decomposition of the manybody interactions of the water tetramer is shown in Figure 9. While compared to the trimer case the magnitudes of the individual components are distinctly larger, qualitatively the situation is similar to the trimer case. Namely, the main contributions to the manybody interactions in the tetramer stem from the (changes in the) electrostatic and induction energy while again correlation effects are comparably small. The difference between the total EXX manybody interaction energy and the sum of the first and second order energies (excluding the (exchange–)dispersion energy) amounts to −0.15 kcal/mol, which is slightly larger than in case of the trimer. This indicates an increasing importance of higher order interaction energy terms to the manybody interaction energy for larger cluster sizes.
Dependence of energy contributions on the cluster size
Various contributions to the interaction energy of the water clusters are presented in Table 1 up to the water tridecamer. The second to fourth column show the twobody, threebody and fourbody dispersion energies, the fifth column contains the sum of all twobody interactions () and the seventh column the total (allbody) interaction energy, i.e., E_{N} = E(123…N) − . In addition, the sixth and the eighth column display the correlation contributions to the twobody and allbody interaction energies, respectively.
Table 1: Two, three and fourbody dispersion energies and twobody and total (allbody) interaction energies of the water clusters. All energies are in kcal/mol.
(H_{2}O)_{x}  dispersion  twobody  allbody  
2  −2.28  −4.47  −1.24  −4.47  −1.24  
3  −8.02  0.06  −12.22  −4.47  −13.76  −4.31  
4  −13.21  0.07  −0.01  −20.27  −6.88  −24.27  −6.61 
5  −16.89  0.04  −0.03  −25.62  −8.33  −32.05  −8.17 
6 (cage)  −24.09  0.31  −0.04  −35.12  −13.79  −40.35  −12.78 
6 (prism)  −24.69  0.33  −0.05  −35.53  −14.25  −40.72  −13.17 
6 (ring)  −20.33  0.01  −0.01  −31.63  −10.06  −39.72  −9.60 
7  −29.86  0.36  −0.06  −43.64  −17.00  −50.71  −15.74 
8  −37.72  0.48  −0.09  −55.15  −21.64  −64.06  −19.99 
9  −41.59  0.44  −0.08  −61.60  −23.53  −72.36  −21.72 
10  −47.41  0.49  −0.12  −70.27  −27.06  −82.32  −24.82 
11  −53.43  0.79  −0.12  −76.72  −31.97  −87.45  −28.86 
12  −61.96  0.91  −0.18  −89.79  −37.14  −102.94  −33.33 
13  −64.94  0.92  −0.19  −94.26  −38.88  −107.82  −34.93 
As can be seen in Table 1, even for large cluster sizes the three and fourbody dispersion energies are fairly small compared to the twobody dispersion interaction. Moreover, since they possess alternate signs, they also partially cancel such that their sum amounts to only 1% of the twobody dispersion at a cluster size of N = 13. One can therefore conclude that manybody dispersion effects are negligible for the description of water clusters. This is also illustrated in the diagram in Figure 10 in which the three and fourbody dispersion interaction is plotted along with the total manybody correlation interaction , defined by the difference of the and terms from Table 1. Here one can see that the total manybody correlation interaction, like the manybody dispersion energy, is repulsive and is a significant contribution to the total manybody interaction energy at larger cluster sizes. For instance, for N = 13 the correlation contribution reduces the overall attractive manybody interaction energy by 22%.
Yet, how significant are manybody interactions compared to the sum of the twobody interactions in the water clusters? From the results shown in Table 1 one finds that the manybody interactions lead to a lowering of the interaction energy of about 13% for all clusters on average. The only exception to this is found for the cyclic structures of the water tetramer, pentamer and hexamer for which the manybody interaction contributes even 20–25% to the total interaction energy. This strong increase of the manybody interactions in the cyclic structures of water clusters is well known [61]. More recently, Bates et al. [72], Hincapie et al. [73] and Chen et al. [74] have performed highlevel coupledcluster calculations for various isomers of (H_{2}O)_{6}. These more recent investigations support the findings of earlier studies of the water hexamer [62,63] that the cyclic structures are less favorable than the prism and cage forms in spite of the strong manybody interaction contribution. The results for the three hexamer structures studied in this work indicate why the ring form is less stable than the other two structures. We find a considerably lower stabilisation of the ringhexamer due to twobody dispersion and (thereby) total twobody interactions by about 4 kcal/mol compared to the cage and prism isomers, see Table 1. This result agrees well with the local molecular orbital energy decomposition analysis of the MP2 interaction energies for the corresponding hexamer structures by Chen and coworkers [74]. Yet, the total interaction energies for the three structures are within a range of 1 kcal/mol (note that compared to this Bates et al. found that the cyclicchair structure is more unstable by +1.83 kcal/mol than the prism structure using the CCSD(T) method [72]). Thus, the decrease of the twobody interaction energy and the increase of the manybody interaction when switching from the cage/prism form to the ring form of (H_{2}O)_{6} almost cancel. Apparently, the water hexamer is the first cluster where the twobody interactions start to dominate the global shape of the cluster geometry, favoring more compact structures than noncompact ring structures.
The above analyis shows that while the manybody interactions in the water clusters are a significant contribution to the total interaction their contribution is scalable and could be modeled, e.g., by scaling the sum of the pair interactions to describe the interactions in the system effectively. This is an approach that is used also in many force field parametrisations for water.
In Figure 11 the ratio of the twobody dispersion energy over the twobody and total interaction energies (including the twobody energies) is plotted. One can see that for all cluster sizes the twobody dispersion amounts to about 60% to the total interaction energy and even almost 70% to the total twobody interaction. This clearly demonstrates the significance of dispersion interactions for the stabilisation of the water clusters. The almost constant dependence of the ratio on the cluster size again demonstrates that the increase of the dispersion interaction energy is very similar to the increase of the total interaction energy. Again, this shows that the magnitude of the total interaction energy of the water clusters (in their equilibrium) could be approximated well by a pairinteraction model.
Conclusion
Different interaction energy contributions to the total interaction energy of water clusters have been analyzed in this work. It has been shown that the main orientation of the water clusters, characterised by the formation of hydrogen bonds, can be well reproduced already on an uncorrelated level using the Hartree–Fock method. However, electron correlation effects to the interaction energy, including the twobody dispersion interaction, lead to a compression of the cluster sizes relative to the structures optimized with the Hartree–Fock method. This global effect originates from the almost isotropic character of the twobody dispersion energy.
Compared to this, the threebody dispersion interaction energy is more strongly dependent on the orientation of the water molecules. However, its was found that dispersion interactions beyond the twobody level are negligible for the description of the stability of water clusters.
The main contributions to the manybody interactions (beyond the twobody level) are described by higher order polarisation interactions, in line with previous studies of the interactions in water clusters [28,29]. This was found through a decomposition of the interaction energies using the DFTSAPT method with the aid of the pseudodimer technique. Overall, manybody interactions are quite significant and contribute about 13% to the total interaction energy of the water clusters. This amount was found to be almost independent on the cluster size. Because of this, manybody interactions in water should be accurately reproduceable by properly scaled twobody interaction energy terms. Since, however, the manybody polarisation interactions, like their twobody counterparts, may be very anisotropic, this approach needs to be carefully tested also for nonequilibrium structures not considered in this work.
It can generally be concluded that electron correlation effects, including dispersion interactions, are crucial for the description of the twobody interactions in water clusters, yet they yield comparably smaller contributions to the manybody interactions. Efficient computational approaches that are based on the manybody expansion to describe the interactions in water may therefore restrict the description of electron correlation to the twobody level without a severe loss in accuracy.
References

Ludwig, R. Phys. Chem. Chem. Phys. 2002, 4, 5481. doi:10.1039/B207000F
Return to citation in text: [1] 
Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. J. Chem. Phys. 1983, 79, 926. doi:10.1063/1.445869
Return to citation in text: [1] [2] 
Horn, H. W.; Swope, W. C.; Pitera, J. W.; Madura, J. D.; Dick, T. J.; Hura, G. L.; HeadGordon, T. J. Chem. Phys. 2004, 120, 9665. doi:10.1063/1.1683075
Return to citation in text: [1] 
Mahoney, M. W.; Jorgensen, W. L. J. Chem. Phys. 2000, 112, 8910. doi:10.1063/1.481505
Return to citation in text: [1] 
Rick, S. W. J. Chem. Phys. 2004, 120, 6085. doi:10.1063/1.1652434
Return to citation in text: [1] 
Fanourgakis, G. S.; Xantheas, S. S. J. Chem. Phys. 2008, 128, 074506. doi:10.1063/1.2837299
Return to citation in text: [1] 
Burnham, C. J.; Anick, D. J.; Mankoo, P. K.; Reiter, G. F. J. Chem. Phys. 2008, 128, 154519. doi:10.1063/1.2895750
Return to citation in text: [1] 
Ren, P. Y.; Ponder, J. W. J. Phys. Chem. B 2004, 108, 13427. doi:10.1021/jp0484332
Return to citation in text: [1] 
Ponder, J. W.; Wu, C.; Ren, P.; Pande, V. S.; Chodera, J. D.; Schnieders, M. J.; Haque, I.; Mobley, D. L.; Lambrecht, D. S.; DiStasio, R. A., Jr.; HeadGordon, M.; Clark, G. N. I.; Johnson, M. E.; HeadGordon, T. J. Phys. Chem. B 2010, 114, 2549. doi:10.1021/jp910674d
Return to citation in text: [1] [2] 
Kumar, R.; Wang, F.F.; Jenness, G. R.; Jordan, K. D. J. Chem. Phys. 2010, 132, 014309. doi:10.1063/1.3276460
Return to citation in text: [1] 
Mas, E. M.; Szalewicz, K.; Bukowski, R.; Jeziorski, B. J. Chem. Phys. 1997, 107, 4207. doi:10.1063/1.474795
Return to citation in text: [1] 
Mas, E. M.; Bukowski, R.; Szalewicz, K.; Groenenboom, G. C.; Wormer, P. E. S.; van der Avoird, A. J. Chem. Phys. 2000, 113, 6687. doi:10.1063/1.1311289
Return to citation in text: [1] 
Cencek, W.; Szalewicz, K.; Leforestier, C.; van Harrevelt, R.; van der Avoird, A. Phys. Chem. Chem. Phys. 2008, 10, 4716. doi:10.1039/b809435g
Return to citation in text: [1] 
Torheyden, M.; Jansen, G. Mol. Phys. 2006, 104, 2101. doi:10.1080/00268970600679188
Return to citation in text: [1] 
Cisneros, G. A.; Wikfeldt, K. T.; Ojamäe, L.; Lu, J.; Xu, Y.; Torabifard, H.; Bartok, A. P.; Csányi, G.; Molinero, V.; Paesani, F. Chem. Rev. 2016, 116, 7501. doi:10.1021/acs.chemrev.5b00644
Return to citation in text: [1] 
Ren, P.; Ponder, J. W. J. Comput. Chem. 2002, 23, 1497. doi:10.1002/jcc.10127
Return to citation in text: [1] 
Burnham, C. J.; Li, J.; Xantheas, S. S.; Leslie, M. J. Chem. Phys. 1999, 110, 4566. doi:10.1063/1.478797
Return to citation in text: [1] 
Gresh, N.; Claverie, P.; Pullman, A. Theor. Chim. Acta 1984, 66, 1. doi:10.1007/BF00577135
Return to citation in text: [1] 
Piquemal, J.P.; Cisneros, G. A.; Reinhardt, P.; Gresh, N.; Darden, T. A. J. Chem. Phys. 2006, 124, 104101. doi:10.1063/1.2173256
Return to citation in text: [1] 
Piquemal, J.P.; Chevreau, H.; Gresh, N. J. Chem. Theory Comput. 2007, 3, 824. doi:10.1021/ct7000182
Return to citation in text: [1] 
Day, P. N.; Jensen, J. H.; Gordon, M. S.; Webb, S. P.; Stevens, W. J.; Krauss, M.; Garmer, D.; Basch, H.; Cohen, D. J. Chem. Phys. 1996, 105, 1968. doi:10.1063/1.472045
Return to citation in text: [1] 
Guidez, E. B.; Gordon, M. S. J. Phys. Chem. A 2017, 121, 3736. doi:10.1021/acs.jpca.6b11403
Return to citation in text: [1] [2] [3] 
Maheshwary, S.; Patel, N.; Sathyamurthy, N.; Kulkarni, A. D.; Gadre, S. R. J. Phys. Chem. A 2001, 105, 10525. doi:10.1021/jp013141b
Return to citation in text: [1] [2] [3] [4] 
Kulkarni, A. D.; Ganesh, V.; Gadre, S. R. J. Chem. Phys. 2004, 121, 5043. doi:10.1063/1.1780156
Return to citation in text: [1] 
Stone, A. The Theory of Intermolecular Forces, 2nd ed.; Oxford University Press, 2013.
Return to citation in text: [1] 
Mas, E. M.; Bukowski, R.; Szalewicz, K. J. Chem. Phys. 2003, 118, 4404. doi:10.1063/1.1542872
Return to citation in text: [1] [2] 
Hodges, M. P.; Stone, A. J.; Xantheas, S. S. J. Phys. Chem. A 1997, 101, 9163. doi:10.1021/jp9716851
Return to citation in text: [1] [2] 
Milet, A.; Moszynski, R.; Wormer, P. E. S.; van der Avoird, A. J. Phys. Chem. A 1999, 103, 6811. doi:10.1021/jp990773d
Return to citation in text: [1] [2] [3] 
Podeszwa, R.; Szalewicz, K. J. Chem. Phys. 2007, 126, 194101. doi:10.1063/1.2733648
Return to citation in text: [1] [2] [3] [4] 
Schütz, M.; Bürgi, T.; Leutwyler, S.; Bürgi, H. B. J. Chem. Phys. 1993, 99, 5228. doi:10.1063/1.465991
Return to citation in text: [1] 
Hapka, M.; Rajchel, Ł.; Modrzejewski, M.; Schäffer, R.; Chałasiński, G.; Szczęśniak, M. M. J. Chem. Phys. 2017, 147, 084106. doi:10.1063/1.4986291
Return to citation in text: [1] 
Misquitta, A. J.; Stone, A. J. J. Chem. Theory Comput. 2008, 4, 7. doi:10.1021/ct700104t
Return to citation in text: [1] 
Heßelmann, A.; Jansen, G. Chem. Phys. Lett. 1999, 315, 248. doi:10.1016/S00092614(99)012518
Return to citation in text: [1] 
Heßelmann, A.; Jansen, G. J. Chem. Phys. 2000, 112, 6949. doi:10.1063/1.481309
Return to citation in text: [1] 
Heßelmann, A.; Görling, A. Mol. Phys. 2010, 108, 359. doi:10.1080/00268970903476662
Return to citation in text: [1] 
Heßelmann, A.; Görling, A. Mol. Phys. 2011, 109, 2473. doi:10.1080/00268976.2011.614282
Return to citation in text: [1] 
Bleiziffer, P.; Heßelmann, A.; Görling, A. J. Chem. Phys. 2012, 136, 134102. doi:10.1063/1.3697845
Return to citation in text: [1] 
Talman, J. D.; Shadwick, W. F. Phys. Rev. A 1976, 14, 36. doi:10.1103/PhysRevA.14.36
Return to citation in text: [1] 
Görling, A. Phys. Rev. Lett. 1999, 83, 5459. doi:10.1103/PhysRevLett.83.5459
Return to citation in text: [1] 
Hirata, S.; Ivanov, S.; Grabowski, I.; Bartlett, R. J.; Burke, K.; Talman, J. D. J. Chem. Phys. 2001, 115, 1635. doi:10.1063/1.1381013
Return to citation in text: [1] 
Holas, A.; Cinal, M. Phys. Rev. A 2005, 72, 032504. doi:10.1103/PhysRevA.72.032504
Return to citation in text: [1] 
HeatonBurgess, T.; Bulat, F. A.; Yang, W. Phys. Rev. Lett. 2007, 98, 256401. doi:10.1103/PhysRevLett.98.256401
Return to citation in text: [1] 
Heßelmann, A.; Götz, A. W.; Della Sala, F.; Görling, A. J. Chem. Phys. 2007, 127, 054102. doi:10.1063/1.2751159
Return to citation in text: [1] 
Heßelmann, A.; Jansen, G. Chem. Phys. Lett. 2002, 357, 464. doi:10.1016/S00092614(02)005389
Return to citation in text: [1] 
Heßelmann, A.; Jansen, G. Chem. Phys. Lett. 2002, 362, 319. doi:10.1016/S00092614(02)010977
Return to citation in text: [1] 
Misquitta, A.; Szalewicz, K. Chem. Phys. Lett. 2002, 357, 301. doi:10.1016/S00092614(02)00533X
Return to citation in text: [1] 
Heßelmann, A.; Jansen, G. Chem. Phys. Lett. 2003, 367, 778. doi:10.1016/S00092614(02)017967
Return to citation in text: [1] 
Misquitta, A. J.; Jeziorski, B.; Szalewicz, K. Phys. Rev. Lett. 2003, 91, 033201. doi:10.1103/PhysRevLett.91.033201
Return to citation in text: [1] 
Heßelmann, A.; Jansen, G. Phys. Chem. Chem. Phys. 2003, 5, 5010. doi:10.1039/b310529f
Return to citation in text: [1] 
Heßelmann, A.; Jansen, G.; Schütz, M. J. Chem. Phys. 2005, 122, 014103. doi:10.1063/1.1824898
Return to citation in text: [1] 
Misquitta, A. J.; Podeszwa, R.; Jeziorski, B.; Szalewicz, K. J. Chem. Phys. 2005, 123, 214103. doi:10.1063/1.2135288
Return to citation in text: [1] 
Jansen, G. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2013, 4, 127. doi:10.1002/wcms.1164
Return to citation in text: [1] 
Görling, A. Int. J. Quantum Chem. 1998, 69, 265. doi:10.1002/(SICI)1097461X(1998)69:3<265::AIDQUA6>3.0.CO;2T
Return to citation in text: [1] [2] 
Görling, A. Phys. Rev. A 1998, 57, 3433. doi:10.1103/PhysRevA.57.3433
Return to citation in text: [1] [2] 
Hirata, S.; Ivanov, S.; Grabowski, I.; Bartlett, R. J. J. Chem. Phys. 2002, 116, 6468. doi:10.1063/1.1460869
Return to citation in text: [1] 
Shigeta, Y.; Hirao, K.; Hirata, S. Phys. Rev. A 2006, 73, 010502. doi:10.1103/PhysRevA.73.010502
Return to citation in text: [1] 
Heßelmann, A.; Ipatov, A.; Görling, A. Phys. Rev. A 2009, 80, 012507. doi:10.1103/PhysRevA.80.012507
Return to citation in text: [1] 
Ipatov, A.; Heßelmann, A.; Görling, A. Int. J. Quantum Chem. 2010, 110, 2202. doi:10.1002/qua.22561
Return to citation in text: [1] 
Heßelmann, A.; Görling, A. J. Chem. Phys. 2011, 134, 034120. doi:10.1063/1.3517312
Return to citation in text: [1] 
Chałasiński, G.; Szcześniak, M. M. Mol. Phys. 1988, 63, 205. doi:10.1080/00268978800100171
Return to citation in text: [1] 
Xantheas, S. S. J. Chem. Phys. 1994, 100, 7523. doi:10.1063/1.466846
Return to citation in text: [1] [2] 
Mhin, B. J.; Kim, J.; Lee, S.; Lee, J. Y.; Kim, K. S. J. Chem. Phys. 1994, 100, 4484. doi:10.1063/1.466279
Return to citation in text: [1] [2] 
Gregory, J. K.; Clary, D. C. J. Phys. Chem. C 1996, 100, 18014. doi:10.1021/jp9616019
Return to citation in text: [1] [2] 
Curtiss, L. A.; Raghavachari, K.; Redfern, P. C.; Rassolov, V.; Pople, J. A. J. Chem. Phys. 1998, 109, 7764. doi:10.1063/1.477422
Return to citation in text: [1] 
Dunning, T. H., Jr. J. Chem. Phys. 1989, 90, 1007. doi:10.1063/1.456153
Return to citation in text: [1] 
Bak, K. L.; Jørgensen, P.; Olsen, J.; Helgaker, T.; Klopper, W. J. Chem. Phys. 2000, 112, 9229. doi:10.1063/1.481544
Return to citation in text: [1] 
Boys, S. F.; Bernardi, F. Mol. Phys. 1970, 19, 553. doi:10.1080/00268977000101561
Return to citation in text: [1] 
MOLPRO, a package of ab initio programs, 2015.1;TTI  TechnologieTransferInitiative GmbH an der Universität Stuttgart: Stuttgart, Germany, 2015. see http://www.molpro.net
Return to citation in text: [1] 
Werner, H.J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 242. doi:10.1002/wcms.82
Return to citation in text: [1] 
Weigend, F.; Ahlrichs, R. Phys. Chem. Chem. Phys. 2005, 7, 3297. doi:10.1039/b508541a
Return to citation in text: [1] 
Axilrod, B. M.; Teller, E. J. Chem. Phys. 1943, 11, 299. doi:10.1063/1.1723844
Return to citation in text: [1] 
Bates, D. M.; Tschumper, G. S. J. Phys. Chem. A 2009, 113, 3555. doi:10.1021/jp8105919
Return to citation in text: [1] [2] 
Hincapie, G.; Acelas, N.; Castaño, M.; David, J.; Restrepo, A. J. Phys. Chem. A 2010, 114, 7809. doi:10.1021/jp103683m
Return to citation in text: [1] 
Chen, Y.; Li, H. J. Phys. Chem. A 2010, 114, 11719. doi:10.1021/jp104822e
Return to citation in text: [1] [2]
34.  Heßelmann, A.; Jansen, G. J. Chem. Phys. 2000, 112, 6949. doi:10.1063/1.481309 
35.  Heßelmann, A.; Görling, A. Mol. Phys. 2010, 108, 359. doi:10.1080/00268970903476662 
36.  Heßelmann, A.; Görling, A. Mol. Phys. 2011, 109, 2473. doi:10.1080/00268976.2011.614282 
37.  Bleiziffer, P.; Heßelmann, A.; Görling, A. J. Chem. Phys. 2012, 136, 134102. doi:10.1063/1.3697845 
38.  Talman, J. D.; Shadwick, W. F. Phys. Rev. A 1976, 14, 36. doi:10.1103/PhysRevA.14.36 
39.  Görling, A. Phys. Rev. Lett. 1999, 83, 5459. doi:10.1103/PhysRevLett.83.5459 
40.  Hirata, S.; Ivanov, S.; Grabowski, I.; Bartlett, R. J.; Burke, K.; Talman, J. D. J. Chem. Phys. 2001, 115, 1635. doi:10.1063/1.1381013 
41.  Holas, A.; Cinal, M. Phys. Rev. A 2005, 72, 032504. doi:10.1103/PhysRevA.72.032504 
42.  HeatonBurgess, T.; Bulat, F. A.; Yang, W. Phys. Rev. Lett. 2007, 98, 256401. doi:10.1103/PhysRevLett.98.256401 
43.  Heßelmann, A.; Götz, A. W.; Della Sala, F.; Görling, A. J. Chem. Phys. 2007, 127, 054102. doi:10.1063/1.2751159 
64.  Curtiss, L. A.; Raghavachari, K.; Redfern, P. C.; Rassolov, V.; Pople, J. A. J. Chem. Phys. 1998, 109, 7764. doi:10.1063/1.477422 
23.  Maheshwary, S.; Patel, N.; Sathyamurthy, N.; Kulkarni, A. D.; Gadre, S. R. J. Phys. Chem. A 2001, 105, 10525. doi:10.1021/jp013141b 
61.  Xantheas, S. S. J. Chem. Phys. 1994, 100, 7523. doi:10.1063/1.466846 
62.  Mhin, B. J.; Kim, J.; Lee, S.; Lee, J. Y.; Kim, K. S. J. Chem. Phys. 1994, 100, 4484. doi:10.1063/1.466279 
63.  Gregory, J. K.; Clary, D. C. J. Phys. Chem. C 1996, 100, 18014. doi:10.1021/jp9616019 
23.  Maheshwary, S.; Patel, N.; Sathyamurthy, N.; Kulkarni, A. D.; Gadre, S. R. J. Phys. Chem. A 2001, 105, 10525. doi:10.1021/jp013141b 
53.  Görling, A. Int. J. Quantum Chem. 1998, 69, 265. doi:10.1002/(SICI)1097461X(1998)69:3<265::AIDQUA6>3.0.CO;2T 
54.  Görling, A. Phys. Rev. A 1998, 57, 3433. doi:10.1103/PhysRevA.57.3433 
60.  Chałasiński, G.; Szcześniak, M. M. Mol. Phys. 1988, 63, 205. doi:10.1080/00268978800100171 
44.  Heßelmann, A.; Jansen, G. Chem. Phys. Lett. 2002, 357, 464. doi:10.1016/S00092614(02)005389 
45.  Heßelmann, A.; Jansen, G. Chem. Phys. Lett. 2002, 362, 319. doi:10.1016/S00092614(02)010977 
46.  Misquitta, A.; Szalewicz, K. Chem. Phys. Lett. 2002, 357, 301. doi:10.1016/S00092614(02)00533X 
47.  Heßelmann, A.; Jansen, G. Chem. Phys. Lett. 2003, 367, 778. doi:10.1016/S00092614(02)017967 
48.  Misquitta, A. J.; Jeziorski, B.; Szalewicz, K. Phys. Rev. Lett. 2003, 91, 033201. doi:10.1103/PhysRevLett.91.033201 
49.  Heßelmann, A.; Jansen, G. Phys. Chem. Chem. Phys. 2003, 5, 5010. doi:10.1039/b310529f 
50.  Heßelmann, A.; Jansen, G.; Schütz, M. J. Chem. Phys. 2005, 122, 014103. doi:10.1063/1.1824898 
51.  Misquitta, A. J.; Podeszwa, R.; Jeziorski, B.; Szalewicz, K. J. Chem. Phys. 2005, 123, 214103. doi:10.1063/1.2135288 
52.  Jansen, G. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2013, 4, 127. doi:10.1002/wcms.1164 
53.  Görling, A. Int. J. Quantum Chem. 1998, 69, 265. doi:10.1002/(SICI)1097461X(1998)69:3<265::AIDQUA6>3.0.CO;2T 
54.  Görling, A. Phys. Rev. A 1998, 57, 3433. doi:10.1103/PhysRevA.57.3433 
55.  Hirata, S.; Ivanov, S.; Grabowski, I.; Bartlett, R. J. J. Chem. Phys. 2002, 116, 6468. doi:10.1063/1.1460869 
56.  Shigeta, Y.; Hirao, K.; Hirata, S. Phys. Rev. A 2006, 73, 010502. doi:10.1103/PhysRevA.73.010502 
57.  Heßelmann, A.; Ipatov, A.; Görling, A. Phys. Rev. A 2009, 80, 012507. doi:10.1103/PhysRevA.80.012507 
58.  Ipatov, A.; Heßelmann, A.; Görling, A. Int. J. Quantum Chem. 2010, 110, 2202. doi:10.1002/qua.22561 
59.  Heßelmann, A.; Görling, A. J. Chem. Phys. 2011, 134, 034120. doi:10.1063/1.3517312 
66.  Bak, K. L.; Jørgensen, P.; Olsen, J.; Helgaker, T.; Klopper, W. J. Chem. Phys. 2000, 112, 9229. doi:10.1063/1.481544 
67.  Boys, S. F.; Bernardi, F. Mol. Phys. 1970, 19, 553. doi:10.1080/00268977000101561 
72.  Bates, D. M.; Tschumper, G. S. J. Phys. Chem. A 2009, 113, 3555. doi:10.1021/jp8105919 
73.  Hincapie, G.; Acelas, N.; Castaño, M.; David, J.; Restrepo, A. J. Phys. Chem. A 2010, 114, 7809. doi:10.1021/jp103683m 
29.  Podeszwa, R.; Szalewicz, K. J. Chem. Phys. 2007, 126, 194101. doi:10.1063/1.2733648 
70.  Weigend, F.; Ahlrichs, R. Phys. Chem. Chem. Phys. 2005, 7, 3297. doi:10.1039/b508541a 
71.  Axilrod, B. M.; Teller, E. J. Chem. Phys. 1943, 11, 299. doi:10.1063/1.1723844 
68.  MOLPRO, a package of ab initio programs, 2015.1;TTI  TechnologieTransferInitiative GmbH an der Universität Stuttgart: Stuttgart, Germany, 2015. see http://www.molpro.net 
69.  Werner, H.J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 242. doi:10.1002/wcms.82 
23.  Maheshwary, S.; Patel, N.; Sathyamurthy, N.; Kulkarni, A. D.; Gadre, S. R. J. Phys. Chem. A 2001, 105, 10525. doi:10.1021/jp013141b 
62.  Mhin, B. J.; Kim, J.; Lee, S.; Lee, J. Y.; Kim, K. S. J. Chem. Phys. 1994, 100, 4484. doi:10.1063/1.466279 
63.  Gregory, J. K.; Clary, D. C. J. Phys. Chem. C 1996, 100, 18014. doi:10.1021/jp9616019 
6.  Fanourgakis, G. S.; Xantheas, S. S. J. Chem. Phys. 2008, 128, 074506. doi:10.1063/1.2837299 
7.  Burnham, C. J.; Anick, D. J.; Mankoo, P. K.; Reiter, G. F. J. Chem. Phys. 2008, 128, 154519. doi:10.1063/1.2895750 
22.  Guidez, E. B.; Gordon, M. S. J. Phys. Chem. A 2017, 121, 3736. doi:10.1021/acs.jpca.6b11403 
4.  Mahoney, M. W.; Jorgensen, W. L. J. Chem. Phys. 2000, 112, 8910. doi:10.1063/1.481505 
5.  Rick, S. W. J. Chem. Phys. 2004, 120, 6085. doi:10.1063/1.1652434 
22.  Guidez, E. B.; Gordon, M. S. J. Phys. Chem. A 2017, 121, 3736. doi:10.1021/acs.jpca.6b11403 
2.  Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. J. Chem. Phys. 1983, 79, 926. doi:10.1063/1.445869 
3.  Horn, H. W.; Swope, W. C.; Pitera, J. W.; Madura, J. D.; Dick, T. J.; Hura, G. L.; HeadGordon, T. J. Chem. Phys. 2004, 120, 9665. doi:10.1063/1.1683075 
18.  Gresh, N.; Claverie, P.; Pullman, A. Theor. Chim. Acta 1984, 66, 1. doi:10.1007/BF00577135 
19.  Piquemal, J.P.; Cisneros, G. A.; Reinhardt, P.; Gresh, N.; Darden, T. A. J. Chem. Phys. 2006, 124, 104101. doi:10.1063/1.2173256 
20.  Piquemal, J.P.; Chevreau, H.; Gresh, N. J. Chem. Theory Comput. 2007, 3, 824. doi:10.1021/ct7000182 
2.  Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. J. Chem. Phys. 1983, 79, 926. doi:10.1063/1.445869 
21.  Day, P. N.; Jensen, J. H.; Gordon, M. S.; Webb, S. P.; Stevens, W. J.; Krauss, M.; Garmer, D.; Basch, H.; Cohen, D. J. Chem. Phys. 1996, 105, 1968. doi:10.1063/1.472045 
14.  Torheyden, M.; Jansen, G. Mol. Phys. 2006, 104, 2101. doi:10.1080/00268970600679188 
9.  Ponder, J. W.; Wu, C.; Ren, P.; Pande, V. S.; Chodera, J. D.; Schnieders, M. J.; Haque, I.; Mobley, D. L.; Lambrecht, D. S.; DiStasio, R. A., Jr.; HeadGordon, M.; Clark, G. N. I.; Johnson, M. E.; HeadGordon, T. J. Phys. Chem. B 2010, 114, 2549. doi:10.1021/jp910674d 
16.  Ren, P.; Ponder, J. W. J. Comput. Chem. 2002, 23, 1497. doi:10.1002/jcc.10127 
11.  Mas, E. M.; Szalewicz, K.; Bukowski, R.; Jeziorski, B. J. Chem. Phys. 1997, 107, 4207. doi:10.1063/1.474795 
12.  Mas, E. M.; Bukowski, R.; Szalewicz, K.; Groenenboom, G. C.; Wormer, P. E. S.; van der Avoird, A. J. Chem. Phys. 2000, 113, 6687. doi:10.1063/1.1311289 
13.  Cencek, W.; Szalewicz, K.; Leforestier, C.; van Harrevelt, R.; van der Avoird, A. Phys. Chem. Chem. Phys. 2008, 10, 4716. doi:10.1039/b809435g 
17.  Burnham, C. J.; Li, J.; Xantheas, S. S.; Leslie, M. J. Chem. Phys. 1999, 110, 4566. doi:10.1063/1.478797 
10.  Kumar, R.; Wang, F.F.; Jenness, G. R.; Jordan, K. D. J. Chem. Phys. 2010, 132, 014309. doi:10.1063/1.3276460 
72.  Bates, D. M.; Tschumper, G. S. J. Phys. Chem. A 2009, 113, 3555. doi:10.1021/jp8105919 
8.  Ren, P. Y.; Ponder, J. W. J. Phys. Chem. B 2004, 108, 13427. doi:10.1021/jp0484332 
9.  Ponder, J. W.; Wu, C.; Ren, P.; Pande, V. S.; Chodera, J. D.; Schnieders, M. J.; Haque, I.; Mobley, D. L.; Lambrecht, D. S.; DiStasio, R. A., Jr.; HeadGordon, M.; Clark, G. N. I.; Johnson, M. E.; HeadGordon, T. J. Phys. Chem. B 2010, 114, 2549. doi:10.1021/jp910674d 
15.  Cisneros, G. A.; Wikfeldt, K. T.; Ojamäe, L.; Lu, J.; Xu, Y.; Torabifard, H.; Bartok, A. P.; Csányi, G.; Molinero, V.; Paesani, F. Chem. Rev. 2016, 116, 7501. doi:10.1021/acs.chemrev.5b00644 
28.  Milet, A.; Moszynski, R.; Wormer, P. E. S.; van der Avoird, A. J. Phys. Chem. A 1999, 103, 6811. doi:10.1021/jp990773d 
29.  Podeszwa, R.; Szalewicz, K. J. Chem. Phys. 2007, 126, 194101. doi:10.1063/1.2733648 
26.  Mas, E. M.; Bukowski, R.; Szalewicz, K. J. Chem. Phys. 2003, 118, 4404. doi:10.1063/1.1542872 
22.  Guidez, E. B.; Gordon, M. S. J. Phys. Chem. A 2017, 121, 3736. doi:10.1021/acs.jpca.6b11403 
23.  Maheshwary, S.; Patel, N.; Sathyamurthy, N.; Kulkarni, A. D.; Gadre, S. R. J. Phys. Chem. A 2001, 105, 10525. doi:10.1021/jp013141b 
24.  Kulkarni, A. D.; Ganesh, V.; Gadre, S. R. J. Chem. Phys. 2004, 121, 5043. doi:10.1063/1.1780156 
25.  Stone, A. The Theory of Intermolecular Forces, 2nd ed.; Oxford University Press, 2013. 
32.  Misquitta, A. J.; Stone, A. J. J. Chem. Theory Comput. 2008, 4, 7. doi:10.1021/ct700104t 
33.  Heßelmann, A.; Jansen, G. Chem. Phys. Lett. 1999, 315, 248. doi:10.1016/S00092614(99)012518 
27.  Hodges, M. P.; Stone, A. J.; Xantheas, S. S. J. Phys. Chem. A 1997, 101, 9163. doi:10.1021/jp9716851 
31.  Hapka, M.; Rajchel, Ł.; Modrzejewski, M.; Schäffer, R.; Chałasiński, G.; Szczęśniak, M. M. J. Chem. Phys. 2017, 147, 084106. doi:10.1063/1.4986291 
29.  Podeszwa, R.; Szalewicz, K. J. Chem. Phys. 2007, 126, 194101. doi:10.1063/1.2733648 
30.  Schütz, M.; Bürgi, T.; Leutwyler, S.; Bürgi, H. B. J. Chem. Phys. 1993, 99, 5228. doi:10.1063/1.465991 
26.  Mas, E. M.; Bukowski, R.; Szalewicz, K. J. Chem. Phys. 2003, 118, 4404. doi:10.1063/1.1542872 
27.  Hodges, M. P.; Stone, A. J.; Xantheas, S. S. J. Phys. Chem. A 1997, 101, 9163. doi:10.1021/jp9716851 
28.  Milet, A.; Moszynski, R.; Wormer, P. E. S.; van der Avoird, A. J. Phys. Chem. A 1999, 103, 6811. doi:10.1021/jp990773d 
28.  Milet, A.; Moszynski, R.; Wormer, P. E. S.; van der Avoird, A. J. Phys. Chem. A 1999, 103, 6811. doi:10.1021/jp990773d 
29.  Podeszwa, R.; Szalewicz, K. J. Chem. Phys. 2007, 126, 194101. doi:10.1063/1.2733648 
© 2018 Heßelmann; licensee BeilsteinInstitut.
This is an Open Access article under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Organic Chemistry terms and conditions: (https://www.beilsteinjournals.org/bjoc)