Water force field models

Force Field Models for the Simulation of Liquid Water... 

Flexible force field models require 30 or more beads in Cartesian path integral descriptions to obtain converged intramolecular energies. Rigid water models require only about five beads due to the weaker intermolecular quantum effects. 

The literature force field models were used with the TIP3P water model to calculate the interaction of the CO2 with a single water molecule, as was done in the ab initio calculations presented above [22-24], When the literature CO2 complexes were compared to the ab initio complexes, the Steele model provided the best match with the ab initio data, with an interaction energy and distance of... 

We conclude this section by giving a topical example of the utility of conditional averages in considering molecularly complex systems (Ashbaugh et al, 2004). We considered the RPLC system discussed above (p. 5), but without methanol n-Ci8 alkyl chains, tethered to a planar support, with water as the mobile phase. The backside of the liquid water phase contacts a dilute water vapor truncated by a repulsive wall see Fig. 1.2, p. 7. Thus, it is appropriate to characterize the system as consistent with aqueous liquid-vapor coexistence at low pressure. A standard CHARMM force-field model (MacKerell Jr. et al, 1998) is used, as are standard molecular dynamics procedures - including periodic boimdary conditions - to acquire the data considered here. Our interest is in the interface between the stationary alkyl and the mobile liquid water phases at 300 K. 

In order to obtain more information on the solvation process Yang and Cui performed a so called natural energy decomposition analysis (NEDA) on monomethyl phosphate ester (MMP) solvated in water. They used a supermolecular approach where the solute plus a number of water molecules (up to 34) were treated quantum-mechanically. A further set of water molecules was treated with a force-field model. Their results indicate that there is a substantial charge transfer between the solute and the nearest solvent molecules. The interaction energy due to this transfer was found to amount to some 70-80% of that of the electric interactions. Since MMP forms hydrogen bonds with the water molecules, all results together suggests that for such a system it is important to include the nearest solvent molecules in the quantum-mechanical treatment, whereas a continuum approximation or a force field may not be sufficiently accurate. 

TD-DFT) they calculated the transition energies and dipole moments for NMA both in vacuum and in an aqueous solution. Moreover, in the treatment of the solvent they compared two different approaches, i.e., a polarizable-continuum method (COSMO) and a supermolecule approach. For the latter, the authors performed molecular-dynamics calculations using a force-field model and, subsequently, extracted a cluster containing the solute and 3 water molecules that form hydrogen bonds to the solute. Averages over 90 such configurations were ultimately determined. 

Chevrot et al. reported MD simulations of the interface between water and the partly hygroscopic IL [C4mim][PF6] comparing different force field models for water and two IL models where the ions were charged 1 and 0.9 [97], For the IL-humidity better agreement was found with experiment when the reduced charges were applied. The authors also studied demixing and found the phases to... 

In the present microscopic investigation, it is true that some results seem to be contradictory. However, they may also be viewed as evidence that there is still room for improvement in the molecular force field in the estimation of preferential exclusion parameters with the use of MD simulations. For example, it may be con-jectnred that improvanent of the partial atomic charge distribution of ectoine and alteration of the water force field, e.g., from TIP3P to some other water model, would bring abont an increase in the water accumulation around the ectoine molecules. This would improve the preferential exclusion parameters Vc so that they are in mnch better agreanent with the experimental ones. 

Recently, we studied the water dependence in mixtures of water and the protonated form of Nafion [53] using both standard force field models and an empirical valence bond model to account for the Grotthuss structural diffusion mechanism of aqueous proton transport. Results showed a transition of an irregularly shaped filamentous (cylindrical) structure in the case of low water (A = 5, where A is defined as the ratio of water molecules to sulfonate groups) to a structure more in accord with the above-discussed models of nano-separation, where larger clusters form which are connected by narrow bridges. A comparison of aqueous cluster sizes indicated, for a simulation time of 30 ns, no percolating clusters for A = 5, whereas at A = 10 most water molecules were located in a connected cluster (see [53]). Other structural... 

Gerig has evaluated force field models for trifluoroethanol-water mixture in order to perform a reliable MD simulation of the system. In this study it was very important to obtain good description of solvent fluorine/solute hydrogen NMR cross-relaxation obtained from solute-solvent nuclear Overhauser effects (NOE). An additional requirement was good agreement between experimental and simulated diffusion coefficients as NMR cross relaxation depends on the diffusion of the components in the mfacture. It was found that OPLS-AA (TFE) and TIP5P-E (H2O) did do a reasonably good work. The mixture will be used as a solvent for small peptides in the future. 

In this model of electrostatic in teraction s, two atoms (i and j) have poin t charges tq and qj. The magnitude of the electrostatic energy (V[. , [ ) varies inversely with the distance between the atoms, Rjj. fh e effective dielectric constant is . For in vacuo simulations or simulation s with explicit water rn olecules, the den om in a tor equals uRjj, In some force fields, a distance-dependent dielectric, where the denominator is uRjj Rjj, represen is solvent implicitly. 

---