Molecular dynamics solvents

Specific solute-solvent interactions involving the first solvation shell only can be treated in detail by discrete solvent models. The various approaches like point charge models, siipennoleciilar calculations, quantum theories of reactions in solution, and their implementations in Monte Carlo methods and molecular dynamics simulations like the Car-Parrinello method are discussed elsewhere in this encyclopedia. Here only some points will be briefly mentioned that seem of relevance for later sections. 

Predicting the solvent or density dependence of rate constants by equation (A3.6.29) or equation (A3.6.31) requires the same ingredients as the calculation of TST rate constants plus an estimate of and a suitable model for the friction coefficient y and its density dependence. While in the framework of molecular dynamics simulations it may be worthwhile to numerically calculate friction coefficients from the average of the relevant time correlation fiinctions, for practical purposes in the analysis of kinetic data it is much more convenient and instructive to use experimentally detemiined macroscopic solvent parameters. 

Wang W, Nelson K A, Xiao L and Coker D F 1994 Molecular dynamics simulation studies of solvent cage effects on photodissociation in condensed phases J. Chem. Phys. 101 9663-71... 

Many of the fiindamental physical and chemical processes at surfaces and interfaces occur on extremely fast time scales. For example, atomic and molecular motions take place on time scales as short as 100 fs, while surface electronic states may have lifetimes as short as 10 fs. With the dramatic recent advances in laser tecluiology, however, such time scales have become increasingly accessible. Surface nonlinear optics provides an attractive approach to capture such events directly in the time domain. Some examples of application of the method include probing the dynamics of melting on the time scale of phonon vibrations [82], photoisomerization of molecules [88], molecular dynamics of adsorbates [89, 90], interfacial solvent dynamics [91], transient band-flattening in semiconductors [92] and laser-induced desorption [93]. A review article discussing such time-resolved studies in metals can be found in... 

Mavri, J., Berendsen, H.J.C., Van Gunsteren, W.F. Influence of solvent on intramolecular proton transfer in hydrogen malonate. Molecular dynamics study of tunneling by density matrix evolution and nonequilibrium solvation. J. Phys. Chem. 97 (1993) 13469-13476. 

To enable an atomic interpretation of the AFM experiments, we have developed a molecular dynamics technique to simulate these experiments [49], Prom such force simulations rupture models at atomic resolution were derived and checked by comparisons of the computed rupture forces with the experimental ones. In order to facilitate such checks, the simulations have been set up to resemble the AFM experiment in as many details as possible (Fig. 4, bottom) the protein-ligand complex was simulated in atomic detail starting from the crystal structure, water solvent was included within the simulation system to account for solvation effects, the protein was held in place by keeping its center of mass fixed (so that internal motions were not hindered), the cantilever was simulated by use of a harmonic spring potential and, finally, the simulated cantilever was connected to the particular atom of the ligand, to which in the AFM experiment the linker molecule was connected. 

D. D. Humphreys, R. A. Friesner, and B. J. Berne. Simulated annealing of a protein in a continuum solvent by multiple-time-step molecular dynamics. J. Phys. Chem., 99 10674-10685, 1995. 

The explicit definition of water molecules seems to be the best way to represent the bulk properties of the solvent correctly. If only a thin layer of explicitly defined solvent molecules is used (due to hmited computational resources), difficulties may rise to reproduce the bulk behavior of water, especially near the border with the vacuum. Even with the definition of a full solvent environment the results depend on the model used for this purpose. In the relative simple case of TIP3P and SPC, which are widely and successfully used, the atoms of the water molecule have fixed charges and fixed relative orientation. Even without internal motions and the charge polarization ability, TIP3P reproduces the bulk properties of water quite well. For a further discussion of other available solvent models, readers are referred to Chapter VII, Section 1.3.2 of the Handbook. Unfortunately, the more sophisticated the water models are (to reproduce the physical properties and thermodynamics of this outstanding solvent correctly), the more impractical they are for being used within molecular dynamics simulations. 

HyperChem uses th e ril 31 water m odel for solvation. You can place th e solute in a box of T1P3P water m oleeules an d impose periodic boun dary eon dition s. You may then turn off the boundary conditions for specific geometry optimi/.aiion or molecular dynamics calculations. However, th is produces undesirable edge effects at the solvent-vacuum interface. 

Often yon need to add solvent molecules to a solute before running a molecular dynamics simiilatmn (see also Solvation and Periodic Boundary Conditions" on page 62). In HyperChem, choose Periodic Box on the Setup m en ii to enclose a soln te in a periodic box filled appropriately with TIP3P models of water inole-cii les. 

The input to a minimisation program consists of a set of initial coordinates for the system. The initial coordinates may come from a variety of sources. They may be obtained from an experimental technique, such as X-ray crystallography or NMR. In other cases a theoretical method is employed, such as a conformational search algorithm. A combination of experimenfal and theoretical approaches may also be used. For example, to study the behaviour of a protein in water one may take an X-ray structure of the protein and immerse it in a solvent bath, where the coordinates of the solvent molecules have been obtained from a Monte Carlo or molecular dynamics simulation. 

Humphreys D D, R A Friesner and B ] Berne 1995. Simulated Annealing of a Protein in a Continuu Solvent by Multiple Time-step Molecular Dynamics. Journal of Physical Chemistry 99 10674-1068... 

In this chapter we shall consider four important problems in molecular n iudelling. First, v discuss the problem of calculating free energies. We then consider continuum solve models, which enable the effects of the solvent to be incorporated into a calculation witho requiring the solvent molecules to be represented explicitly. Third, we shall consider the simi lation of chemical reactions, including the important technique of ab initio molecular dynamic Finally, we consider how to study the nature of defects in solid-state materials. 

---