Water model

Water model

effects.

An alternative to the explicit water models is to use an implicit solvation model, also known as a continuum model.

imple 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 potential for models such as TIP3P and TIP4P is represented by

:E_{ab} = sum_{i} ^{on a} sum_{j} ^{on b} frac {k_Cq_iq_j}{r_{ij + frac {A}{r_{OO}^{12 - frac {B}{r_{OO}^6}

where "kC", the electrostatic constant, has a value of 332.1 Å·kcal/mol in the units commonly used in molecular modeling; "qi" are the partial charges relative to the charge of the electron; "rij" is the distance between two atoms or charged sites; and "A" and "B" are the Lennard-Jones parameters. The charged sites may be on the atoms or on dummy sites (such as lone pairs). In most water models, the Lennard-Jones term applies only to the interaction between the oxygen atoms.

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.

:

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.

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)":

:S(r_{ij}) = egin{cases} 0, & mbox{if }r_{ij} le R_L \ frac{(r_{12} - R_L)^2(3R_U - R_L - 2r_{ij})}{(R_U - R_L)^2}, & mbox{if }R_L le r_{ij} le R_U \ 1, & mbox{if }R_U le r_{ij}end{cases}

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 [H. Nada, J.P.J.M. van der Eerden, "J. Chem. Phys." 2003, "118", 7401. doi|10.1063/1.1562610] . 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. [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]
*Coarse-grained models. One- and two-site models of water have also been developed. [ S. Izvekov, G. A. Voth. Multiscale coarse graining of liquid-state systems "J. Chem. Phys. 2005, "123", 134105. doi|10.1063/1.2038787] In coarse grain models, each site can represent several water molecules.

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).

ee also

*Force field (chemistry)

References


Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать реферат

Look at other dictionaries:

  • Water — This article is about general aspects of water. For a detailed discussion of its properties, see Properties of water. For other uses, see Water (disambiguation) …   Wikipedia

  • Water cluster — In chemistry a water cluster is a discrete hydrogen bonded assembly or cluster of molecules of water [cite journal | title = Water: From Clusters to the Bulk | author = Ralf Ludwig | journal = Angew. Chem. Int. Ed. | year = 2001 | volume = 40 |… …   Wikipedia

  • Water (molecule) — Chembox new Name = Water (H2O) ImageFileL1 = Water 2D labelled.png ImageSizeL1 = 120px ImageNameL1 = The water molecule has this basic geometric structure ImageFileR1 = Water molecule 3D.svg ImageSizeR1 = 100px ImageNameR1 = Space filling model… …   Wikipedia

  • Water privatization — is a short hand for private sector participation in the provision of water services and sanitation, although more rarely it refers to privatization of water resources themselves. Because water services are seen as such a key public service,… …   Wikipedia

  • Water trading — is the process of buying and selling of water access entitlements, also often called water rights . The terms of the trade can be either permanent or temporary, depending on the legal status of the water rights. Some of the western states of the… …   Wikipedia

  • Water supply and sanitation in Latin America — is characterized by insufficient access and in many cases by poor service quality, with detrimental impacts on public health. Water and sanitation services are provided by a vast array of mostly local service providers under an often fragmented… …   Wikipedia

  • Water-line model — Water line Wa ter line 1. (Shipbuilding) Any one of certain lines of a vessel, model, or plan, parallel with the surface of the water at various heights from the keel. [1913 Webster] Note: In a half breadth plan, the water lines are outward… …   The Collaborative International Dictionary of English

  • Water heating — is a thermodynamic process using an energy source to heat water above its initial temperature. Typical domestic uses of hot water are for cooking, cleaning, bathing, and space heating. In industry, both hot water and water heated to steam have… …   Wikipedia

  • Water line — Wa ter line 1. (Shipbuilding) Any one of certain lines of a vessel, model, or plan, parallel with the surface of the water at various heights from the keel. [1913 Webster] Note: In a half breadth plan, the water lines are outward curves showing… …   The Collaborative International Dictionary of English

  • Model Town, Lahore — Model Town Lahore is one of residential suburbs of Lahore, Punjab, Pakistan. It is located next to Faisal Town, Gulberg, Garden Town and Township suburbs of Lahore. Contents 1 History 1.1 Model Town Society 1.2 Evacuee property …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”