Time modeling molecular dynamics

In our implementation of SMD, modified versions of VMD and Sigma communicate with each other using a customized, lightweight protocol. Sigma sends atomic positions resulting from each molecular dynamics time step to VMD for display. When the user specifies restraints on parts of the displayed model, VMD sends them to Sigma, where they are converted into potential-well restraints added to the force field [21]. 

In the classical (CD) or molecular dynamics (MD) model, the movement of all atoms inside a solid is studied as a function of time. The model thus takes the interaction with all neighboring atoms into account. The starting point is Newton s equation for the motion of a single atom i, which is governed by the interaction forces tfom all target atoms in the neighborhood. 

A third approach is to inject particles based on a grand canonical ensemble distribution. At each predetermined molecular dynamics time step, the probability to create or destroy a particle is calculated and a random number is used to determine whether the update is accepted (the probability for both the creation and the destruction of a particle must be equal to ensure reversibility). The probability function depends on the excess chemical potential and must be calculated in a way that is consistent with the microscopic model used to describe the system. In the work of Im et al., a primitive water model is used, and the chemical potential is determined through an analytic solution to the Ornstein-Zernike equation using the hypemetted chain as a closure relation. This method is very accurate from the physical viewpoint, but it has a poorer CPU performance compared with simpler schemes based on... 

The method of molecular dynamics (numerical modeling of the dynamic type) combines both theoretical and experimental approaches. The method specifies the interaction potentials of molecules (atoms), the temperature as a chaotic distribution of velocities, and the laws of motion (in our case, the integration of Newton s laws for each molecule). At the same time, molecular dynamic simulation is an experiment, which enables one to observe the behavior of molecules (atoms) with a resolution of up to 10" s in the time domain and up to 10 m in the space domain, which can t be realized in any other experiment. 

In 2001, Hanke, Price and Lynden-BelP were the first to conduct an atomistic simulation of compovmds that can be called ionic liquids under our definition. They used molecular dynamics to model the crystalline state of 1,3-dimethylimidazolium chloride ([Cimim][Cl]), 1,3-dimethylimidazolium hexafluorophosphate ([CimimllPFg]), l-ethyl-3-methylimidazolium chloride ([C2mim][Cl]), and l-ethyl-3-methylimidazolium hexafluorophosphate ([C2 mimJlPFg]). They also modeled the liquid state of [Cimim][Cl] and [Cimim] [PFg], both of which are relatively high melting substances. Because of this (and the need to speed dynamics and thus limit computation times), the liquid simulations were carried out at temperatures between 400 and 500 K. The form of the potential function they used was... 

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. 

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. 

In molecular mechanics and molecular dynamics studies of proteins, assig-ment of standard, non-dynamical ionization states of protein titratable groups is a common practice. This assumption seems to be well justified because proton exchange times between protein and solution usually far exceed the time range of the MD simulations. We investigated to what extent the assumed protonation state of a protein influences its molecular dynamics trajectory, and how often our titration algorithm predicted ionization states identical to those imposed on the groups, when applied to a set of structures derived from a molecular dynamics trajectory [34]. As a model we took the bovine... 

In molecular dynamics applications there is a growing interest in mixed quantum-classical models various kinds of which have been proposed in the current literature. We will concentrate on two of these models the adiabatic or time-dependent Born-Oppenheimer (BO) model, [8, 13], and the so-called QCMD model. Both models describe most atoms of the molecular system by the means of classical mechanics but an important, small portion of the system by the means of a wavefunction. In the BO model this wavefunction is adiabatically coupled to the classical motion while the QCMD model consists of a singularly perturbed Schrddinger equation nonlinearly coupled to classical Newtonian equations, 2.2. 

Abstract. The overall Hamiltonian structure of the Quantum-Classical Molecular Dynamics model makes - analogously to classical molecular dynamics - symplectic integration schemes the methods of choice for long-term simulations. This has already been demonstrated by the symplectic PICKABACK method [19]. However, this method requires a relatively small step-size due to the high-frequency quantum modes. Therefore, following related ideas from classical molecular dynamics, we investigate symplectic multiple-time-stepping methods and indicate various possibilities to overcome the step-size limitation of PICKABACK. 

In this paper we present a number of time integrators for various problems ranging from classical to quantum molecular dynamics. These integrators share some common features they are new, they are second-order accurate and time-reversible, they improve substantially over standard schemes in well-defined model situations — and none of them has been tested on real applications at the time of this writing. This last feature will hopefully change in the near future [20]. 

Finite difference techniques are used to generate molecular dynamics trajectories with continuous potential models, which we will assume to be pairwise additive. The essential idea is that the integration is broken down into many small stages, each separated in time by a fixed time 6t. The total force on each particle in the configuration at a time t is calculated as the vector sum of its interactions with other particles. From the force we can determine the accelerations of the particles, which are then combined with the positions and velocities at a time t to calculate the positions and velocities at a time t + 6t. The force is assumed to be constant during the time step. The forces on the particles in their new positions are then determined, leading to new positions and velocities at time t - - 2St, and so on. 

Once a PES has been computed, it is often fitted to an analytic function. This is done because there are many ways to analyze analytic functions that require much less computation time than working directly with ah initio calculations. For example, the reaction can be modeled as a molecular dynamics simulation showing the vibrational motion and reaction trajectories as described in Chapter 19. Another technique is to fit ah initio results to a semiempirical model designed for the purpose of describing PES s. 

---