To use all functions of this page, please activate cookies in your browser.
my.chemeurope.com
With an accout for my.chemeurope.com you can always see everything at a glance – and you can configure your own website and individual newsletter.
 My watch list
 My saved searches
 My saved topics
 My newsletter
Water modelIn computational chemistry, classical water models are used for the simulation of liquid water and aqueous solutions (with socalled explicit solvent, as opposed to the implicit solvation models). These models generally use the approximations of molecular mechanics. Many different models have been proposed; they can be classified by the number of points used to define the model (atoms plus dummy sites), whether the structure is rigid or flexible, and whether the model includes polarization effects. Additional recommended knowledge
Simple water modelsThe simplest water models treat the water molecule as rigid and rely only on nonbonded interactions. The electrostatic interaction is modeled using Coulomb's law and the dispersion and repulsion forces using the LennardJones potential. The charged sites may be on the atoms or on dummy sites (such as lone pairs). The LennardJones term typically applies to the oxygen atoms only. The figure below shows the general shape of the 3 to 6site water models. The exact geometric parameters (the OH distance and the HOH angle) vary depending on the model. 3siteThe simplest models have three interaction sites, corresponding to the three atoms of the water molecule. Each atom gets assigned a point charge, and the oxygen atom also gets the LennardJones parameters. The 3site models are very popular for molecular dynamics simulations because of their simplicity and computational efficiency. Most models use a rigid geometry matching the known geometry of the water molecule. An exception is the SPC model, which assumes an ideal tetrahedral shape (HOH angle of 109.47°) instead of the observed angle of 104.5°. The table below lists the parameters for some 3site models.
The SPC/E model adds an average polarization correction to the potential energy function: where μ is the dipole of the effectively polarized water molecule (2.35 D for the SPC/E model), μ^{0} is the dipole moment of an isolated water molecule (1.85 D from experiment), and α_{i} is an isotropic polarizability constant, with a value of 1.608 × 10^{−40} F m. Since the charges in the model are constant, this correction just results in adding 1.25 kcal/mol (5.22 kJ/mol) to the total energy. The SPC/E model results in a better density and diffusion constant than the SPC model. The TIP3P model implemented in the CHARMM force field is a slightly modified version of the original. The difference lies in the LennardJones parameters: unlike TIP3P, the CHARMM version of the model places LennardJones parameters on the hydrogen atoms. The charges are not modified.^{[5]} Other models:
4siteThe 4site models place the negative charge on a dummy atom (labeled M in the figure) placed near the oxygen along the bisector of the HOH angle. This improves the electrostatic distribution around the water molecule. The first model to use this approach was the BernalFowler model published in 1933, which may also be the earliest water model. However, the BF model doesn't reproduce well the bulk properties of water, such as density and heat of vaporization, and is therefore only of historical interest. This is a consequence of the parameterization method; newer models, developed after modern computers became available, were parameterized by running Metropolis Monte Carlo or molecular dynamics simulations and adjusting the parameters until the bulk properties are reproduced well enough. The TIP4P model, first published in 1983, is widely implemented in computational chemistry software packages and often used for the simulation of biomolecular systems. There have been subsequent reparameterizations of the TIP4P model for specific uses: the TIP4PEw model, for use with Ewald summation methods; the TIP4P/Ice, for simulation of solid water ice; and TIP4P/2005, a general parameterization for simulating the entire phase diagram of water.
Others:
5siteThe 5site models place the negative charge on dummy atoms (labeled L) representing the lone pairs of the oxygen atom, with a tetrahedrallike geometry. An early model of these types was the BNS model of BenNaim and Stillinger, proposed in 1971, soon succeeded by the ST2 model of Stillinger and Rahman in 1974. Mainly due to their higher computational cost, fivesite models were not developed much until 2000, when the TIP5P model of Mahoney and Jorgensen was published. When compared with earlier models, the TIP5P model results in improvements in the geometry for the water dimer, a more "tetrahedral" water structure that better reproduces the experimental radial distribution functions from neutron diffraction, and the temperature of maximum density of water. The TIP5PE model is a reparameterization of TIP5P for use with Ewald sums.
Note, however, that the BNS and ST2 models do not use Coulomb's law directly for the electrostatic terms, but a modified version that is scaled down at short distances by multiplying it by the switching function S(r): Therefore the R_{L} and R_{U} parameters only apply to BNS and ST2. 6siteA 6site model that combines all the sites of the 4 and 5site models was developed by Nada and van der Eerden^{[14]}. It was found to reproduce the structure and melting of ice better than other models. Other
Computational costThe computational cost of a water simulation increases with the number of interaction sites in the water model. The CPU time is approximately proportional to the number of interatomic distances that need to be computed. For the 3site model, 9 distances are required for each pair of water molecules (every atom of one molecule against every atom of the other molecule, or 3 × 3). For the 4site model, 10 distances are required (every charged site with every charged site, plus the OO interaction, or 3 × 3 + 1). For the 5site model, 17 distances are required (4 × 4 + 1). Finally, for the 6site model, 26 distances are required (5 × 5 + 1). When using rigid water models in molecular dynamics, there is an additional cost associated with keeping the structure constrained, using constraint algorithms (although with bond lengths constrained it is often possible to increase the time step). See alsoReferences


This article is licensed under the GNU Free Documentation License. It uses material from the Wikipedia article "Water_model". A list of authors is available in Wikipedia. 