Molecular dynamics motion

What is the mechanism of spins dropping down from the state to the a state and fanning out around the two cones, and what determines the rates (R = HT and/ 2 = 1/72) of NMR relaxation These processes are intimately tied to the motion of molecules as they tumble ( reorient ) in solution in their rapid Brownian motion, and measurement of the NMR relaxation parameters T and T2 can even give us detailed information about molecular dynamics (motion) from the point of view of each spin in the molecule. A simplified model... 

We adopted the method by Car and Parrinello [22] to ensure a self-consistent evolution of the electronic structure during molecular dynamics motion. The electronic structure was described in the framework of density functional theory (DFT) with the generalized gradient approximation (GGA) due to Becke (B) for the exchange energy and Lee, Yang and Parr (LYP) for the correlation energy... 

Two simulation methods—Monte Carlo and molecular dynamics—allow calculation of the density profile and pressure difference of Eq. III-44 across the vapor-liquid interface [64, 65]. In the former method, the initial system consists of N molecules in assumed positions. An intermolecule potential function is chosen, such as the Lennard-Jones potential, and the positions are randomly varied until the energy of the system is at a minimum. The resulting configuration is taken to be the equilibrium one. In the molecular dynamics approach, the N molecules are given initial positions and velocities and the equations of motion are solved to follow the ensuing collisions until the set shows constant time-average thermodynamic properties. Both methods are computer intensive yet widely used. 

Because of limitations of space, this section concentrates very little on rotational motion and its interaction with the vibrations of a molecule. However, this is an extremely important aspect of molecular dynamics of long-standing interest, and with development of new methods it is the focus of mtense investigation [18, 19, 20. 21. 22 and 23]. One very interesting aspect of rotation-vibration dynamics involving geometric phases is addressed in section A1.2.20. 

In equilibrium statistical mechanics, one is concerned with the thennodynamic and other macroscopic properties of matter. The aim is to derive these properties from the laws of molecular dynamics and thus create a link between microscopic molecular motion and thennodynamic behaviour. A typical macroscopic system is composed of a large number A of molecules occupying a volume V which is large compared to that occupied by a molecule ... 

The alternative simulation approaches are based on molecular dynamics calculations. This is conceptually simpler that the Monte Carlo method the equations of motion are solved for a system of A molecules, and periodic boundary conditions are again imposed. This method pennits both the equilibrium and transport properties of the system to be evaluated, essentially by numerically solvmg the equations of motion... 

If we wish to know the number of (VpV)-collisions that actually take place in this small time interval, we need to know exactly where each particle is located and then follow the motion of all the particles from time tto time t+ bt. In fact, this is what is done in computer simulated molecular dynamics. We wish to avoid this exact specification of the particle trajectories, and instead carry out a plausible argument for the computation of r To do this, Boltzmann made the following assumption, called the Stosszahlansatz, which we encountered already in the calculation of the mean free path ... 

If a fluid is placed between two concentric cylinders, and the inner cylinder rotated, a complex fluid dynamical motion known as Taylor-Couette flow is established. Mass transport is then by exchange between eddy vortices which can, under some conditions, be imagmed as a substantially enlranced diflfiisivity (typically with effective diflfiision coefficients several orders of magnitude above molecular difhision coefficients) that can be altered by varying the rotation rate, and with all species having the same diffusivity. Studies of the BZ and CIMA/CDIMA systems in such a Couette reactor [45] have revealed bifiircation tlirough a complex sequence of front patterns, see figure A3.14.16. 

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... 

Molecular dynamics consists of the brute-force solution of Newton s equations of motion. It is necessary to encode in the program the potential energy and force law of interaction between molecules the equations of motion are solved numerically, by finite difference techniques. The system evolution corresponds closely to what happens in real life and allows us to calculate dynamical properties, as well as thennodynamic and structural fiinctions. For a range of molecular models, packaged routines are available, either connnercially or tlirough the academic conmuinity. 

Ryckaert J-P, Ciccotti G and Berendsen H J C 1977 Numerical integration of the Cartesian equations of motion of a system with constraints molecular dynamics of n-alkanes J. Comput. Phys. 23 327-41... 

The study of quautum effects associated with nuclear motion is a distinct field of chemistry, known as quantum molecular dynamics. This section gives an overview of the methodology of the field for fiirtlier reading, consult [1, 2, 3, 4 and 5,]. 

In its most fiindamental fonn, quantum molecular dynamics is associated with solving the Sclirodinger equation for molecular motion, whether using a single electronic surface (as in the Bom-Oppenlieimer approximation— section B3.4.2 or with the inclusion of multiple electronic states, which is important when discussing non-adiabatic effects, in which tire electronic state is changed [15,16, YL, 18 and 19]. 

As in any field, it is usefiil to clarify tenninology. Tliroughout this section an atom more specifically refers to its nuclear centre. Also, for most of the section the /)= 1 convention is used. Finally, it should be noted that in the literature the label quantum molecular dynamics is also sometimes used for a purely classical description of atomic motion under the potential created by tlie electronic distribution. 

The complexity of polymeric systems make tire development of an analytical model to predict tlieir stmctural and dynamical properties difficult. Therefore, numerical computer simulations of polymers are widely used to bridge tire gap between tire tlieoretical concepts and the experimental results. Computer simulations can also help tire prediction of material properties and provide detailed insights into tire behaviour of polymer systems. A simulation is based on two elements a more or less detailed model of tire polymer and a related force field which allows tire calculation of tire energy and tire motion of tire system using molecular mechanisms, molecular dynamics, or Monte Carlo teclmiques 1631. 

Molecular dynamics tracks tire temporal evolution of a microscopic model system tlirough numerical integration of tire equations of motion for tire degrees of freedom considered. The main asset of molecular dynamics is tliat it provides directly a wealtli of detailed infonnation on dynamical processes. 

The method of molecular dynamics (MD), described earlier in this book, is a powerful approach for simulating the dynamics and predicting the rates of chemical reactions. In the MD approach most commonly used, the potential of interaction is specified between atoms participating in the reaction, and the time evolution of their positions is obtained by solving Hamilton s equations for the classical motions of the nuclei. Because MD simulations of etching reactions must include a significant number of atoms from the substrate as well as the gaseous etchant species, the calculations become computationally intensive, and the time scale of the simulation is limited to the... 

The starting point for the theory of molecular dynamics, and indeed the basis for most of theoretical chemistry, is the separation of the nuclear and electionic motion. In the standard, adiabatic, picture this leads to the concept of nuclei moving over PES corresponding to the electronic states of a system. 

Creveld, L., Amadei, A., Van Schaik, C., Pepermans, R., De Vlieg, J., Berendsen, H.J.C. Identification of functional and unfolding motions of cutinase as obtained from molecular dynamics computer simulations. Submitted (1998). 

Hayward, S., Kitao, A., Berendsen, H.J.C. Model-free methods to analyze domain motions in proteins from simulation A comparison of normal mode analysis and molecular dynamics simulation of lysozyme. Proteins 27 (1997) 425-437. 

The simulations also revealed that flapping motions of one of the loops of the avidin monomer play a crucial role in the mechanism of the unbinding of biotin. The fluctuation time for this loop as well as the relaxation time for many of the processes in proteins can be on the order of microseconds and longer (Eaton et al., 1997). The loop has enough time to fluctuate into an open state on experimental time scales (1 ms), but the fluctuation time is too long for this event to take place on the nanosecond time scale of simulations. To facilitate the exit of biotin from its binding pocket, the conformation of this loop was altered (Izrailev et al., 1997) using the interactive molecular dynamics features of MDScope (Nelson et al., 1995 Nelson et al., 1996 Humphrey et al., 1996). 

Molecular dynamics simulations ([McCammon and Harvey 1987]) propagate an atomistic system by iteratively solving Newton s equation of motion for each atomic particle. Due to computational constraints, simulations can only be extended to a typical time scale of 1 ns currently, and conformational transitions such as protein domains movements are unlikely to be observed. 

The important underlying components of protein motion during a simulation can be extracted by a Principal Component Analysis (PGA). It stands for a diagonalization of the variance-covariance matrix R of the mass-weighted internal displacements during a molecular dynamics simulation. 

An interesting approach has recently been chosen in the MBO(N)D program ([Moldyn 1997]). Structural elements of different size varying from individual peptide planes up to protein domains can be defined to be rigid. During an atomistic molecular dynamics simulation, all fast motion orthogonal to the lowest normal modes is removed. This allows use of ca. 20 times longer time steps than in standard simulations. 

Related to the previous method, a simulation scheme was recently derived from the Onsager-Machlup action that combines atomistic simulations with a reaction path approach ([Oleander and Elber 1996]). Here, time steps up to 100 times larger than in standard molecular dynamics simulations were used to produce approximate trajectories by the following equations of motion ... 

Elamrani et al. 1996] Elamrani, S., Berry, M.B., Phillips Jr., G.N., McCammon, J.A. Study of Global Motions in Proteins by Weighted Masses Molecular Dynamics Adenylate Kinase as a Test Case. Proteins 25 (1996) 79-88 [Elcock et al. 1997] Elcock, A.H., Potter, M.J., McCammon, J.A. Application of Poisson-Boltzmann Solvation Forces to Macromolecular Simulations. In Computer Simulation of Biomoleeular Systems, Vol. 3, A.J. Wilkinson et al. eds., ESCOM Science Publishers B.V., Leiden... 

Abstract. Molecular dynamics (MD) simulations of proteins provide descriptions of atomic motions, which allow to relate observable properties of proteins to microscopic processes. Unfortunately, such MD simulations require an enormous amount of computer time and, therefore, are limited to time scales of nanoseconds. We describe first a fast multiple time step structure adapted multipole method (FA-MUSAMM) to speed up the evaluation of the computationally most demanding Coulomb interactions in solvated protein models, secondly an application of this method aiming at a microscopic understanding of single molecule atomic force microscopy experiments, and, thirdly, a new method to predict slow conformational motions at microsecond time scales. 

---