TIP4P

TIk experimentally determined dipole moment of a water molecule in the gas phase is 1.85 D. The dipole moment of an individual water molecule calculated with any of thv se simple models is significantly higher for example, the SPC dipole moment is 2.27 D and that for TIP4P is 2.18 D. These values are much closer to the effective dipole moment of liquid water, which is approximately 2.6 D. These models are thus all effective pairwise models. The simple water models are usually parametrised by calculating various pmperties using molecular dynamics or Monte Carlo simulations and then modifying the... 

The simulated free surface of liquid water is relatively stable for several nanoseconds [68-72] because of the strong hydrogen bonds formed by liquid water. The density decrease near the interface is smooth it is possible to describe it by a hyperbolic tangent function [70]. The width of the interface, measured by the distance between the positions where the density equals 90% and 10% of the bulk density, is about 5 A at room temperature [70,71]. The left side of Fig. 3 shows a typical density profile of the free interface for the TIP4P water model [73]. 

FIG. 3 Left density profile, p z), from a 500 ps simulation of a thin film consisting of 200 TIP4P water molecules at room temperature. Right orientational distribution, p cos d), with 3 the angle between the molecular dipole moment p and the surface normal z. The vertical lines in the left plot indicate the boundary z-ranges,... 

Fig. 5 shows data from a simulation of TIP4P water that is confined on both sides by a rhombohedral mercury crystal with (111) surface structure. Bosio et al. [135] deduce from their X-ray studies that a solid o-mercury lattice with a larger lattice constant in the z direction may be used as a good structural model for liquid mercury. Thus, the mercury phase was modeled as a rigid crystal in order to simplify the simulations. The surface of such a crystal shows rather low corrugation. 

A complete set of intermolecular potential functions has been developed for use in computer simulations of proteins in their native environment. Parameters have been reported for 25 peptide residues as well as the common neutral and charged terminal groups. The potential functions have the simple Coulomb plus Lennard-Jones form and are compatible with the widely used models for water, TIP4P, TIP3P and SPC. The parameters were obtained and tested primarily in conjunction with Monte Carlo statistical mechanics simulations of 36 pure organic liquids and numerous aqueous solutions of organic ions representative of subunits in the side chains and backbones of proteins... 

The polarizable fluctuating charge model in CHARMM results from the work of Patel, Brooks and co-workers [92, 214], The water model is based on the TIP4P-FQ model of Rick, Stuart and Berne [17], In the development of the force field the electronegativities and hardnesses were treated as empirical parameters and do not have any association with experimental or QM values, for example, from ionization energies and electron affinities of single atoms. 

W. L. Jorgensen, J. F. Blake, and J. K. Buckner, Free energy of TIP4P water and the free... 

---