Molecular dynamics velocity correlation function

Figure 3 Normalised single particle velocity correlation functions and their normal mode approximations ip (<)/i/) (0) and i/ ° (<)/0 (O). Left panel centre-of-rnass velocity correlation function for CO, 80 K, 1 bar right panel rotational velocity correlation function for CSj, 293 K, 10 kbar. Thick lines Molecular dynamics result dotted lines normal mode result deluding unstable modes thin lines normal mode result including unstable modes.

For translational motion, the velocity correlation function is particularly useful and, as we will show, can be utilized to provide a relationship between the echo amplitude and the molecular dynamics in the case of general modulation wave forms. Its Fourier spectrum is simply the selfdiffusion tensor (Lenk, 1977 Stepisnik, 1981) D (w), where a and /3 may take each of the Cartesian directions, x, y, z, that is. 

The dynamical analysis of the isothermal molecular dynamics simulations [359] was presented in Ref. 360. The motion perpendicular to the surface for the first-layer molecules is for all coverages of the damped oscillatory type around the minimum of the total potential normal to the plane the frequency increases slightly as the number of second-layer molecules increases. The translational velocity correlation function parallel to the surface stems from an apparently liquid-like phase of the first-layer molecules. Inspection of trajectory plots suggests that much of the observed first-layer motion is in directions parallel to the rows of molecules and can be inter-... 

In the molecular dynamics comparison of the full and the purely repulsive Lennard-Jones potentials, the velocity correlation function was calculated using both potentials in a variety of systems. However, the extremely close... 

MD simulation is advantageous for obtaining dynamic properties directly, since the MD technique provides not only particle positions but also particle velocities that enable us to utilize the response theory (e.g., the Kubo formula [175,176]) to calculate the transport coefficients from time-dependent correlation functions. For example, we will examine the self-diffusion process of a tagged PFPE molecular center of mass (Fig. 1.49) from the simulation to gain insight into the excitation of translational motion, specifically, spreading and replenishment. The squared displacement of the center mass of a molecule or a bead is used as a measure of translational movement. The self-diffusion coefficient D can be represented as a velocity autocorrelation function... 

Figure 6. (a) 2-Gilorobutane at 50 K, 6x6 site-site potential, angular velocity autocorrelation functions. Crosshatdiing indicates computer noise difference between R and S enantiomers. (—) Racemic mixture, (b) As for (a), under the influence of a strong field E, producing a torque — 63 XE in each molecule of the molecular dynamics sample. (1) (—) Racemic mixture (2) (—) R enantiomer. Ordinate Normalized correlation function abscissa time, ps. 

Computer simulations of the molecular dynamics of the liquid state (see also Chapter VI) clearly show that the correlation function of the velocity variable is not exponential rather it usually exhibits a sort of damped oscillatory behavior. This means that the Markovian assumption is often invalid. This makes it n sary, when studying a chemical reaction in a liquid phase, to replace the standard Kramers condition [see Eq. (4b)] with a more realistic correlation function having a finite lifetime. Recall the rate expression obtained by Kramers for moderate to high frictions, Eq. (6). This may be cast into the form k = tst/(" >y) where given by Eq. (7), is essentially an equilibrium property depraiding on the thermodynamic equilibrium inside the well. As a canonical equilibrium property, it is not afifected by whether or not the system is Markovian. The calculation of the factor fiui, y) depends, however, on the dynamics of the system and will thus be modified when non-Markovian behavior is allowed for. 

Correlations by Computation of Molecular Dynamics. The power of modem computing systems has made it possible to solve the dynamical equations of motion of a model system of several hundred molecules, with fairly realistic interaction potentials, and hence by direct calculation obtain correlation functions for linear velocity, angular velodty, dipole orientation, etc. Rahman s classic paper on the motion of 864 atoms of model argon has stimulated a great amount of further work, of which we cite particularly that of Beme and Harpon nitrogen and carbon monoxide, and that of Rahman himself and Stillinger on water. ... 

Computer simulation of molecular dynamics is concerned with solving numerically the simultaneous equations of motion for a few hundred atoms or molecules that interact via specified potentials. One thus obtains the coordinates and velocities of the ensemble as a function of time that describe the structure and correlations of the sample. If a model of the induced polarizabilities is adopted, the spectral lineshapes can be obtained, often with certain quantum corrections [425,426]. One primary concern is, of course, to account as accurately as possible for the pairwise interactions so that by carefully comparing the calculated with the measured band shapes, new information concerning the effects of irreducible contributions of inter-molecular potential and cluster polarizabilities can be identified eventually. Pioneering work has pointed out significant effects of irreducible long-range forces of the Axilrod-Teller triple-dipole type [10]. Very recently, on the basis of combined computer simulation and experimental CILS studies, claims have been made that irreducible three-body contributions are observable, for example, in dense krypton [221]. 

Our preliminary MD simulations provided similar results as the MC method for the calculation of the static properties of confined PFPE nanofilms, especially the radius of gyration and end-bead density profiles. The anisotropic molecular conformation and experimental layering structures in the film were also verified. MD simulations provide a powerful tool for examining the dynamics of nanofilms through correlation functions by tracking the trajectories of molecules, including the space and velocity coordinates. MD simulations, therefore, are suitable for... 

The correlation function, <-P2[am(0) ( )]>. provides a measure of the internal motions of particular residues in the protein.324 333 Figure 46 shows the results obtained for Trp-62 and Trp-63 from the stochastic boundary molecular dynamics simulations of lysozyme used to analyze the displacement and velocity autocorrelation functions. The net influence of the solvent for both Trp-62 and Trp-63 is to cause a slower decay in the anisotropy than occurs in vacuum. In vacuum, the anisotropy decays to a plateau value of 0.36 to 0.37 (relative to the initial value of 0.4) for both residues within a picosecond. In solution there is an initial rapid decay, corresponding to that found in vacuum, followed by a slower decay (without reaching a plateau value) that continues beyond the period (10 ps) over which the correlation function is ex-... 

G. J. Martyna, J. Chem. Phys. (in press, 1996). In this paper, an effective set of molecular dynamics equations are specified that provide an alternative path-integral approach to the calculation of position and velocity time correlation functions. This approach is essentially based on the Wigner phase-space function. For general nonlinear systems, the appropriate MD mass in this approach is not the physical mass, but it must instead be a position-dependent effective mass. 

Space does not permit us to give here a detailed discussion of the effects of finite-system size on molecular dynamic calculations of time correlations functions. We have given elsewhere a discussion of such effects on the velocity autocorrelation function from a hydrodynamic point of view. This reference can also be consulted for a more extensive discussion of results for both the velocity autocorrelation function and the super-Burnett self-diffusion coefficient, including comparisons with theoretical predictions. 

Consider a dense gas of hard spheres, all with mass m and diameter a. Since the collisions of hard-sphere molecules are instantaneous, the probability is zero that any particle will collide with more than one particle at a time. Hence we still suppose that the dynamical events taking place in the gas are made up of binary collisions, and that to derive an equation for the single-particle distribution function /(r, v, /) we need only take binary collisions into account. However, the Stosszahlansatz used in deriving the Boltzmann equation for a dilute gas should be modified to take into account any spatial and velocity correlations that may exist between the colliding spheres. The Enskog theory continues to ignore the possibility of correlations in the velocities before collision, but attempts to take into account the spatial correlations. In addition, the Enskog theory takes into account the variation of the distribution function over distances of the order of the molecular diameter, which also leads to corrections to the Boltzmann equation. 

Finally, we mention the fact that the molecular dynamics calculation of Po t) for a gas of hard disks exhibited a vortex type of velocity correlation between the tagged molecule and the surrounding molecules that is very similar to the hydrodynamic flow field surrounding a moving volume element in a fluid initially at rest. This vortex pattern, illustrated in Fig. 26, suggests that a fraction of the momentum transferred by the tagged particle to the particles in front of it is eventually returned to it from behind. This process causes the velocity autocorrelation function to be larger than it would be if these vortices did not occur, and it is connected with the slow decay of the velocity autocorrelation function. 

---