Relaxation time molecular dynamics simulation

Molecular dynamics simulations of proteins often begin with a known structure (such as an X-ray diffraction structure) that you want to maintain during equilibration. Since the solvent may contain high energy hot spots, equilibration of the protein and solvent at the same time can change the protein conformation. To avoid this, select only the water molecules and run a molecular dynamics equilibration. This relaxes the water while fixing the protein structure. Then deselect the water and equilibrate the whole system. 

One would prefer to be able to calculate aU of them by molecular dynamics simulations, exclusively. This is unfortunately not possible at present. In fact, some indices p, v of Eq. (6) refer to electronically excited molecules, which decay through population relaxation on the pico- and nanosecond time scales. The other indices p, v denote molecules that remain in their electronic ground state, and hydrodynamic time scales beyond microseconds intervene. The presence of these long times precludes the exclusive use of molecular dynamics, and a recourse to hydrodynamics of continuous media is inevitable. This concession has a high price. Macroscopic hydrodynamics assume a local thermodynamic equilibrium, which does not exist at times prior to 100 ps. These times are thus excluded from these studies. 

Even if we consider a single solvent, e g., water, at a single temperature, say 298K, depends on the solute and in fact on the coordinate of the solute which is under consideration, and we cannot take xF as a constant. Nevertheless, in the absence of a molecular dynamics simulation for the solute motion of interest, XF for polar solvents like water is often approximated by the Debye model. In this model, the dielectric polarization of the solvent relaxes as a single exponential with a relaxation time equal to the rotational (i.e., reorientational) relaxation time of a single molecule, which is called Tp) or the Debye time [32, 347], The Debye time may be associated with the relaxation of the transverse component of the polarization field. However the solvent fluctuations and frictional relaxation occur on a faster scale given by [348,349]... 

Monte Carlo and Molecular Dynamics simulations of water near hydrophobic surfaces have yielded a wealth of information about the structure, thermodynamics and transport properties of interfacial water. In particular, they have demonstrated the presence of molecular layering and density oscillations which extend many Angstroms away from the surfaces. These oscillations have recently been verified experimentally. Ordered dipolar orientations and reduced dipole relaxation times are observed in most of the simulations, indicating that interfacial water is not a uniform dielectric continuum. Reduced dipole relaxation times near the surfaces indicate that interfacial water experiences hindered rotation. The majority of simulation results indicate that water near hydrophobic surfaces exhibits fewer hydrogen bonds than water near the midplane. 

Figure 1 The shear stress relaxation function, C(t), obtained from a molecular dynamics simulation of500 SRP spheres at a reduced temperature of 1.0 and effective volume fraction of 0.45. Note that n = 144 and 1152 (from Equation (1)) cases are superimposable with the analytic function of Equation (4) ( Algebraic on the figure) for short times, t (or nt here)...

Engstrdm et al. [112] used molecular dynamics simulations to study quadrupole relaxation mechanism for Li+, Na+, and Cl ions in dilute aqueous solutions. They found that NMR relaxation rate for these ions was determined by the relaxation of water molecules in the first solvation shell. The simulations show nonexponential solvation dynamics which can be modeled by two relaxation time constants < 0.1 ps and x2 lps (see Table 4). 

The hydration dynamics were also studied by Karim et al. [58], These authors used the TIP4P model of water in their molecular dynamics simulations. The observed hydration dynamics was nonexponential with average time constant in the range of 0.4-0.7ps. In this case simulated relaxation of the first solvation shell was also faster than that of the other shells. 

The behavior of VACF and of D in one-dimensional systems are, therefore, of special interest. The transverse current mode of course does not exist here, and the decay of the longitudinal current mode (related to the dynamic structure factor by a trivial time differentiation) is sufficiently fast to preclude the existence of any "dangerous" long-time tail. Actually, Jepsen [181] was the first to derive die closed-form expression for the VACF and the diffusion coeffident for hard rods. His study showed that in the long time VACF decays as 1/f3, in contrast to the t d 2 dependence reported for the two and three dimensions. Lebowitz and Percus [182] studied the short-time behavior of VACF and made an exponential approximation for VACF [i.e, Cv(f) = e 2 ], for the short times. Haus and Raveche [183] carried out the extensive molecular dynamic simulations to study relaxation of an initially ordered array in one dimension. This study also investigated the 1/f3 behavior of VACF. However, none of the above studies provides a physical explanation of the 1/f3 dependence of VACF at long times, of the type that exists for two and three dimensions. 

At all events, the role of Coulombic forces for VET in solution was first examined in a molecular dynamics simulation for the 680 cm-1 C-Cl vibration of the CH3CI molecule (modeled as a diatomic) in water solvent by Whitnell et al. (19). The (classical) relaxation time Ti was determined both by nonequilibrium simulations and by use of the classical Landau-Teller (LT) formula (1,3,19,20). 

The nature of the excitation has a profound influence on the subsequent relaxation of molecular Uquid systems, as the molecular dynamics simulations show. This influence can be exerted at field-on equiUbrium and in decay transients (the deexdtation effect). GrigoUni has shown that the effect of high-intensity excitation is to slow the time decay of the envelope of such oscillatory functions as the angular velocity autocorrelation function. The effect of high-intensity pulses is the same as that of ultrafast (subpicosecond laser) pulses. The computer simulation by Abbot and Oxtoby shows that... 

In our previous papers , we have shown that collective jump motions of atoms take place in highly supercooled fluid states, mainly contributing to the a relaxation, and therefore represents the molecular-level mechanisms. The main purpose of this paper is to study both a and / relaxations from S q,u>) and x (9,w) in a supercooled fluid by a super-long-time molecular dynamics (MD) simulation for a model fluid of binary soft-sphere mixtures. In particular, we focus on studying the type of each relaxation (Debye or non-Debye ) and the molecular-level processes for the / relaxation. 

The calculated results agreed qualitatively with that by molecular dynamics simulations (Fig. 1). In the long-time region, solvent relaxation for a change in a solute charge from 0 to e (z=0- l) was slower than that obtained by the linearized equation. In addition, relaxation for z=0- 1 was slower than that for z= 1 - 0. On the other hand, in the short-time region, solvent relaxation for z=0—1 was similar to that by the linearized equation. 

---