Information about Water Model
)
A water model is defined by its geometry, together with other parameters such as the atomic charges and Lennard-Jones parameters.
Simple water models
The simplest water models treat the water molecule as rigid and rely only on non-bonded interactions. The electrostatic interaction is modeled using Coulomb's law and the dispersion and repulsion forces using the Lennard-Jones potential.
The charged sites may be on the atoms or on dummy sites (such as lone pairs). The Lennard-Jones term typically applies to the oxygen atoms only.
The figure below shows the general shape of the 3- to 6-site water models. The exact geometric parameters (the OH distance and the HOH angle) vary depending on the model.
){kind=spacer}
3-site
The 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 Lennard-Jones parameters. The 3-site 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 3-site models.
| TIPS[1] | SPC[2] | TIP3P[3] | SPC/E[4] | |
|---|---|---|---|---|
| r(OH), Ã… | 0.9572 | 1.0 | 0.9572 | 1.0 |
| HOH, deg | 104.52 | 109.47 | 104.52 | 109.47 |
| A × 10−3, kcal Ã…12/mol | 580.0 | 629.4 | 582.0 | 629.4 |
| B, kcal Ã…6/mol | 525.0 | 625.5 | 595.0 | 625.5 |
| q(O) | −0.80 | −0.82 | −0.834 | −0.8476 |
| q(H) | +0.40 | +0.41 | +0.417 | +0.4238 |
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.
Other models:
- Fergunson (flex. SPC)
- CVFF (flex.)
4-site
The 4-site 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 Bernal-Fowler 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 TIP4P-Ew 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.
| BF[5] | TIPS2[6] | TIP4P<ref name="tip3p" /> | TIP4P-Ew[7] | TIP4P/Ice[8] | TIP4P/2005[9] | |
|---|---|---|---|---|---|---|
| r(OH), Ã… | 0.96 | 0.9572 | 0.9572 | 0.9572 | 0.9572 | 0.9572 |
| HOH, deg | 105.7 | 104.52 | 104.52 | 104.52 | 104.52 | 104.52 |
| r(OM), Ã… | 0.15 | 0.15 | 0.15 | 0.125 | 0.1577 | 0.1546 |
| A × 10−3, kcal Ã…12/mol | 560.4 | 695.0 | 600.0 | 656.1 | 857.9 | 731.3 |
| B, kcal Ã…6/mol | 837.0 | 600.0 | 610.0 | 653.5 | 850.5 | 736.0 |
| q(M) | −0.98 | −1.07 | −1.04 | −1.04844 | −1.1794 | −1.1128 |
| q(H) | +0.49 | +0.535 | +0.52 | +0.52422 | +0.5897 | +0.5564 |
Others:
- TIP4PF (flexible)
5-site
The 5-site models place the negative charge on dummy atoms (labeled L) representing the lone pairs of the oxygen atom, with a tetrahedral-like geometry. An early model of these types was the BNS model of Ben-Naim and Stillinger, proposed in 1971, soon succeeded by the ST2 model of Stillinger and Rahman in 1974. Mainly due to their higher computational cost, five-site 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 TIP5P-E model is a reparameterization of TIP5P for use with Ewald sums.| BNS[10] | ST2<ref name="st2" /> | TIP5P[11] | TIP5P-E[12] | |
|---|---|---|---|---|
| r(OH), Ã… | 1.0 | 1.0 | 0.9572 | 0.9572 |
| HOH, deg | 109.47 | 109.47 | 104.52 | 104.52 |
| r(OL), Ã… | 1.0 | 0.8 | 0.70 | 0.70 |
| LOL, deg | 109.47 | 109.47 | 109.47 | 109.47 |
| A × 10−3, kcal Ã…12/mol | 77.4 | 238.7 | 544.5 | 590.3 |
| B, kcal Ã…6/mol | 153.8 | 268.9 | 554.3 | 628.2 |
| q(L) | −0.19562 | −0.2357 | −0.241 | −0.241 |
| q(H) | +0.19562 | +0.2357 | +0.241 | +0.241 |
| RL, Ã… | 2.0379 | 2.0160 | ||
| RU, Ã… | 3.1877 | 3.1287 |
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 RL and RU parameters only apply to BNS and ST2.
6-site
A 6-site model that combines all the sites of the 4- and 5-site models was developed by Nada and van der Eerden[13]. It was found to reproduce the structure and melting of ice better than other models.Other
- MB model. A more abstract model resembling the Mercedes-Benz logo that reproduces some features of water in two-dimensional systems. It is not used as such for simulations of "real" (i.e., three-dimensional) systems, but it is useful for qualitative studies and for educational purposes.[14]
- Coarse-grained models. One- and two-site models of water have also been developed.[15]
Computational cost
The 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 3-site 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 4-site model, 10 distances are required (every charged site with every charged site, plus the O-O interaction, or 3 × 3 + 1). For the 5-site model, 17 distances are required (4 × 4 + 1). Finally, for the 6-site 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 also
References
1. ^ Jorgensen, W. L. Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water. J. Am. Chem. Soc. 1981, 103, 335-340.
2. ^ H.J.C. Berendsen, J.P.M. Postma, W.F. van Gunsteren, and J. Hermans, In Intermolecular Forces, edited by B. Pullman (Reidel, Dordrecht, 1981), p. 331.
3. ^ Jorgensen, W. L.; Chandrasekhar, J.; Madura, J.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys 1983, 79, 926-935. doi:10.1063/1.445869
4. ^ H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma. The Missing Term in Effective Pair Potentials. J. Phys. Chem 1987, 91, 6269-6271. doi:10.1021/j100308a038
5. ^ Bernal, J. D.; Fowler, R.H. J. Chem. Phys. 1933, 1, 515. doi:10.1063/1.1749327
6. ^ Jorgensen, W. L. Revised TIPS for simulations of liquid water and aqueous solutions. J. Chem. Phys 1982, 77, 4156-4163. doi:10.1063/1.444325
7. ^ H. W. Horn, W. C. Swope, J. W. Pitera, J. D. Madura, T. J. Dick, G. L. Hura, and T. Head-Gordon. Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. J. Chem. Phys. 2004, 120, 9665-9678. doi:10.1063/1.1683075
8. ^ J. L. F. Abascal, E. Sanz, R. GarcÃa Fernández, and C. Vega. A potential model for the study of ices and amorphous water: TIP4P/Ice. J. Chem. Phys. 2005, 122, 234511. doi:10.1063/1.1931662
9. ^ J. L. F. Abascal and C. Vega. A general purpose model for the condensed phases of water: TIP4P/2005. J. Chem. Phys. 2005, 123, 234505. doi:10.1063/1.2121687
10. ^ F.H. Stillinger, A. Rahman, Improved simulation of liquid water by molecular dynamics. J. Chem. Phys. 1974, 60, 1545-1557. doi:10.1063/1.1681229
11. ^ Mahoney, M. W.; Jorgensen, W. L. A five-site model liquid water and the reproduction of the density anomaly by rigid, non-polarizable models. J. Chem. Phys. 2000, 112, 8910-8922. doi:10.1063/1.481505
12. ^ Rick, S. W. A reoptimization of the five-site water potential (TIP5P) for use with Ewald sums. J. Chem. Phys. 2004, 120, 6085-6093. doi:10.1063/1.1652434
13. ^ H. Nada, J.P.J.M. van der Eerden, J. Chem. Phys. 2003, 118, 7401. doi:10.1063/1.1562610
14. ^ K. A. T. Silverstein, A. D. J. Haymet, and K. A. Dill. A Simple Model of Water and the Hydrophobic Effect. J. Am. Chem. Soc. 1998, 120, 3166-3175. doi:10.1021/ja973029k
15. ^ S. Izvekov, G. A. Voth. Multiscale coarse graining of liquid-state systems J. Chem. Phys. 2005, 123'', 134105. doi:10.1063/1.2038787
2. ^ H.J.C. Berendsen, J.P.M. Postma, W.F. van Gunsteren, and J. Hermans, In Intermolecular Forces, edited by B. Pullman (Reidel, Dordrecht, 1981), p. 331.
3. ^ Jorgensen, W. L.; Chandrasekhar, J.; Madura, J.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys 1983, 79, 926-935. doi:10.1063/1.445869
4. ^ H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma. The Missing Term in Effective Pair Potentials. J. Phys. Chem 1987, 91, 6269-6271. doi:10.1021/j100308a038
5. ^ Bernal, J. D.; Fowler, R.H. J. Chem. Phys. 1933, 1, 515. doi:10.1063/1.1749327
6. ^ Jorgensen, W. L. Revised TIPS for simulations of liquid water and aqueous solutions. J. Chem. Phys 1982, 77, 4156-4163. doi:10.1063/1.444325
7. ^ H. W. Horn, W. C. Swope, J. W. Pitera, J. D. Madura, T. J. Dick, G. L. Hura, and T. Head-Gordon. Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. J. Chem. Phys. 2004, 120, 9665-9678. doi:10.1063/1.1683075
8. ^ J. L. F. Abascal, E. Sanz, R. GarcÃa Fernández, and C. Vega. A potential model for the study of ices and amorphous water: TIP4P/Ice. J. Chem. Phys. 2005, 122, 234511. doi:10.1063/1.1931662
9. ^ J. L. F. Abascal and C. Vega. A general purpose model for the condensed phases of water: TIP4P/2005. J. Chem. Phys. 2005, 123, 234505. doi:10.1063/1.2121687
10. ^ F.H. Stillinger, A. Rahman, Improved simulation of liquid water by molecular dynamics. J. Chem. Phys. 1974, 60, 1545-1557. doi:10.1063/1.1681229
11. ^ Mahoney, M. W.; Jorgensen, W. L. A five-site model liquid water and the reproduction of the density anomaly by rigid, non-polarizable models. J. Chem. Phys. 2000, 112, 8910-8922. doi:10.1063/1.481505
12. ^ Rick, S. W. A reoptimization of the five-site water potential (TIP5P) for use with Ewald sums. J. Chem. Phys. 2004, 120, 6085-6093. doi:10.1063/1.1652434
13. ^ H. Nada, J.P.J.M. van der Eerden, J. Chem. Phys. 2003, 118, 7401. doi:10.1063/1.1562610
14. ^ K. A. T. Silverstein, A. D. J. Haymet, and K. A. Dill. A Simple Model of Water and the Hydrophobic Effect. J. Am. Chem. Soc. 1998, 120, 3166-3175. doi:10.1021/ja973029k
15. ^ S. Izvekov, G. A. Voth. Multiscale coarse graining of liquid-state systems J. Chem. Phys. 2005, 123'', 134105. doi:10.1063/1.2038787
Computational chemistry is a branch of chemistry that uses computers to assist in solving chemical problems. It uses the results of theoretical chemistry, incorporated into efficient computer programs, to calculate the structures and properties of molecules and solids.
..... Click the link for more information.
..... Click the link for more information.
Water (H2O, HOH) is the most abundant molecule on Earth's surface, composing of about 70% of the Earth's surface as liquid and solid state in addition to being found in the atmosphere as a vapor.
..... Click the link for more information.
..... Click the link for more information.
Implicit solvation (sometimes known as continuum solvation) is a method of representing solvent as a continuous medium instead of individual “explicit” solvent molecules most often used in molecular dynamics simulations and in other applications of molecular
..... Click the link for more information.
..... Click the link for more information.
molecular mechanics refers to the use of Newtonian mechanics to model molecular systems. The potential energy of all systems in molecular mechanics is calculated using force fields.
..... Click the link for more information.
..... Click the link for more information.
polarization (Brit., polarisation) is the property of electromagnetic waves, such as light, that describes the direction of the transverse electric field. More generally, the polarization of a transverse wave describes the direction of oscillation in the plane
..... Click the link for more information.
..... Click the link for more information.
In physics, chemistry, and biology, intermolecular forces are forces that act between stable molecules or between functional groups of macromolecules. These non-covalent forces, which give rise to bonding energies of less than a few kcal/mol, are generally much weaker than the
..... Click the link for more information.
..... Click the link for more information.
Coulomb's law, developed in the 1780s by French physicist Charles Augustin de Coulomb, may be stated as follows:
..... Click the link for more information.
- The magnitude of the electrostatic force between two points electric charges is directly proportional to the product of the magnitudes of each
..... Click the link for more information.
Lennard-Jones potential (also referred to as the L-J potential, 6-12 potential or, less commonly, 12-6 potential) is a simple mathematical model that represents this behavior. It was proposed in 1931 by John Lennard-Jones of Bristol University.
..... Click the link for more information.
..... Click the link for more information.
A partial charge is a charge with an absolute value of less than one elementary charge unit (that is, smaller than the charge of the electron).
..... Click the link for more information.
Partial atomic charges
Partial charges are created due to the asymmetric distribution of electrons in chemical bonds...... Click the link for more information.
Molecular dynamics (MD) is a form of computer simulation wherein atoms and molecules are allowed to interact for a period of time under known laws of physics, giving a view of the motion of the atoms.
..... Click the link for more information.
..... Click the link for more information.
bisection is the division of something into two equal parts, usually by a line, which is then called a bisector. The most often considered types of bisectors are segment bisectors and angle bisectors.
..... Click the link for more information.
..... Click the link for more information.
In physics, density is mass m per unit volume V—how heavy something is compared to its size. A small, heavy object, such as a rock or a lump of lead, is denser than a lighter object of the same size or a larger object of the same weight, such as pieces of
..... Click the link for more information.
..... Click the link for more information.
The enthalpy of vaporization, (symbol ), also known as the heat of vaporization or heat of evaporation, is the energy required to transform a given quantity of a substance into a gas.
..... Click the link for more information.
..... Click the link for more information.
Metropolis-Hastings algorithm is a rejection sampling algorithm used to generate a sequence of samples from a probability distribution that is difficult to sample from directly. This sequence can be used in Markov chain Monte Carlo simulation to approximate the distribution (i.e.
..... Click the link for more information.
..... Click the link for more information.
In physical chemistry, mineralogy, and materials science, a phase diagram is a type of graph used to show the equilibrium conditions between the thermodynamically-distinct phases.
..... Click the link for more information.
..... Click the link for more information.
A lone pair is a (valence) electron pair without bonding or sharing with other atoms. They are found in the outermost electron shell of an atom, so lone pairs are a subset of a molecule's valence electrons.
..... Click the link for more information.
..... Click the link for more information.
The water dimer consists of two water molecules loosely bound by a hydrogen bond. It is the smallest water cluster. Because it is the simplest model system for studying hydrogen bonding in water, it has been the target of so many theoretical (and later experimental) studies that
..... Click the link for more information.
..... Click the link for more information.
In computational mechanics and statistical mechanics, a radial distribution function (RDF), g(r), describes how the density of surrounding matter varies as a function of the distance from a distinguished point.
..... Click the link for more information.
..... Click the link for more information.
Neutron diffraction is a crystallographic method for the determination of the atomic structure of a material. This is a form of elastic scattering where the neutrons exiting the experiment have more or less the same energy as the incident neutrons.
..... Click the link for more information.
..... Click the link for more information.
This article may be too long.
Please discuss this issue on the talk page and help summarize or split the content into subarticles of an article series.
..... Click the link for more information.
Please discuss this issue on the talk page and help summarize or split the content into subarticles of an article series.
- For other uses of the name Mercedes, see Mercedes.
..... Click the link for more information.
In mechanics, a constraint algorithm is a method for satisfying constraints for bodies that obey Newton's equations of motion. There are three basic approaches to satisfying such constraints: choosing novel unconstrained coordinates ("internal coordinates"), introducing explicit
..... Click the link for more information.
..... Click the link for more information.
force field (also called a forcefield) refers to the functional form and parameter sets used to describe the potential energy of a system of particles (typically but not necessarily atoms).
..... Click the link for more information.
..... Click the link for more information.
digital object identifier (or DOI) is a permanent identifier given to a document, which is not related to its current location. A typical use of a DOI is to give a scientific paper or article a unique identifying number that can be used by anyone to locate details of the paper, and
..... Click the link for more information.
..... Click the link for more information.
digital object identifier (or DOI) is a permanent identifier given to a document, which is not related to its current location. A typical use of a DOI is to give a scientific paper or article a unique identifying number that can be used by anyone to locate details of the paper, and
..... Click the link for more information.
..... Click the link for more information.
digital object identifier (or DOI) is a permanent identifier given to a document, which is not related to its current location. A typical use of a DOI is to give a scientific paper or article a unique identifying number that can be used by anyone to locate details of the paper, and
..... Click the link for more information.
..... Click the link for more information.
digital object identifier (or DOI) is a permanent identifier given to a document, which is not related to its current location. A typical use of a DOI is to give a scientific paper or article a unique identifying number that can be used by anyone to locate details of the paper, and
..... Click the link for more information.
..... Click the link for more information.
digital object identifier (or DOI) is a permanent identifier given to a document, which is not related to its current location. A typical use of a DOI is to give a scientific paper or article a unique identifying number that can be used by anyone to locate details of the paper, and
..... Click the link for more information.
..... Click the link for more information.
digital object identifier (or DOI) is a permanent identifier given to a document, which is not related to its current location. A typical use of a DOI is to give a scientific paper or article a unique identifying number that can be used by anyone to locate details of the paper, and
..... Click the link for more information.
..... Click the link for more information.
digital object identifier (or DOI) is a permanent identifier given to a document, which is not related to its current location. A typical use of a DOI is to give a scientific paper or article a unique identifying number that can be used by anyone to locate details of the paper, and
..... Click the link for more information.
..... Click the link for more information.
digital object identifier (or DOI) is a permanent identifier given to a document, which is not related to its current location. A typical use of a DOI is to give a scientific paper or article a unique identifying number that can be used by anyone to locate details of the paper, and
..... Click the link for more information.
..... Click the link for more information.
This article is copied from an article on Wikipedia.org - the free encyclopedia created and edited by online user community. The text was not checked or edited by anyone on our staff. Although the vast majority of the wikipedia encyclopedia articles provide accurate and timely information please do not assume the accuracy of any particular article. This article is distributed under the terms of GNU Free Documentation License.
Herod_Archelaus
