Molecular dynamics simulation principles

A molecular dynamics simulation samples the phase space of a molecule (defined by the position of the atoms and their velocities) by integrating Newton s equations of motion. Because MD accounts for thermal motion, the molecules simulated may possess enough thermal energy to overcome potential barriers, which makes the technique suitable in principle for conformational analysis of especially large molecules. In the case of small molecules, other techniques such as systematic, random. Genetic Algorithm-based, or Monte Carlo searches may be better suited for effectively sampling conformational space. 

The principle of the molecular dynamics simulation approach is the movement of atoms under the action of a force field. 

Using the first-principles molecular-dynamics simulation, Munejiri, Shimojo and Hoshino studied the structure of liquid sulfur at 400 K, below the polymerization temperature [79]. They found that some of the Ss ring molecules homolytically open up on excitation of one electron from the HOMO to the LUMO. The chain-like diradicals S " thus generated partly recombine intramolecularly with formation of a branched Sy=S species rather than cyclo-Ss- Furthermore, the authors showed that photo-induced polymerization occurs in liquid sulfur when the Ss chains or Sy=S species are close to each other at their end. The mechanism of polymerization of sulfur remains a challenging problem for further theoretical work. 

In addition to enhancing surface reactions, water can also facilitate surface transport processes. First-principles ab initio molecular dynamics simulations of the aqueous/ metal interface for Rh(l 11) [Vassilev et al., 2002] and PtRu(OOOl) alloy [Desai et al., 2003b] surfaces showed that the aqueous interface enhanced the apparent transport or diffusion of OH intermediates across the metal surface. Adsorbed OH and H2O molecules engage in fast proton transfer, such that OH appears to diffuse across the surface. The oxygen atoms, however, remained fixed at the same positions, and it is only the proton that transfers. Transport occurs via the symmetric reaction... 

Beginning with values of ]r) and ]v) at time 0, one calculates the new positions and then the new velocities. This method is second-order in At, too. For additional details, see Allen, M. R, and D. J. Tildesley, Computer Simulation of Liquids, Clarendon Press, Oxford (1989) Frenkel, D., and B. Smit, Understanding Molecular Simulation, Academic Press (2002) Haile, J. M., Molecular Dynamics Simulation, Wiley (1992) Leach, A. R., Molecular Modelling Principles and Applications, Prentice-Hall (2001) Schlick, T., Molecular Modeling and Simulations, Springer, New York (2002). 

Molecular dynamics simulations are capable of addressing the self-assembly process at a rudimentary, but often impressive, level. These calculations can be used to study the secondary structure (and some tertiary structure) of large complex molecules. Present computers and codes can handle massive calculations but cannot eliminate concerns that boundary conditions may affect the result. Eventually, continued improvements in computer hardware will provide this added capacity in serial computers development of parallel computer codes is likely to accomplish the goal more quickly. In addition, the development of realistic, time-efficient potentials will accelerate the useful application of dynamic simulation to the self-assembly process. In addition, principles are needed to guide the selec-... 

Temperature effects are included explicitly in molecular dynamics simulations by including kinetic energy terms - the balls representing the atoms are now on the move The principles are simple. In the microcanonical ensemble (NVE) ... 

The strategy in a molecular dynamics simulation is conceptually fairly simple. The first step is to consider a set of molecules. Then it is necessary to choose initial positions of all atoms, such that they do not physically overlap, and that all bonds between the atoms have a reasonable length. Subsequently, it is necessary to specify the initial velocities of all the atoms. The velocities must preferably be consistent with the temperature in the system. Finally, and most importantly, it is necessary to define the force-field parameters. In effect the force field defines the potential energy of each atom. This value is a complicated sum of many contributions that can be computed when the distances of a given atom to all other atoms in the system are known. In the simulation, the spatial evolution as well as the velocity evolution of all molecules is found by solving the classical Newton equations of mechanics. The basic outcome of the simulation comprises the coordinates and velocities of all atoms as a function of the time. Thus, structural information, such as lipid conformations or membrane thickness, is readily available. Thermodynamic information is more expensive to obtain, but in principle this can be extracted from a long simulation trajectory. 

In principle, MC algorithms can be tuned for particular systems and can thus be more efficient than MD for obtaining equilibrium distributions. An interesting idea is to use MC simulations to obtain accurate initial guesses for subsequent MD simulations. Already as early as 1993, Venable and co-workers [68] used a scheme for efficiently sampling configurations of individual lipids in a mean field. These configurations were then used to develop the initial conditions for the molecular dynamic simulations. 

Enantioselective Hydrosilylation by Chiral Pd Based Homogeneous Catalysts with First-Principles and Combined QM/MM Molecular Dynamics Simulations... 

The main utility of recent molecular dynamics calculations has been first to determine if the assumption of quadratic free energy curves holds and second, to calculate the reorganization free energy in order to test Eqs. (24) and (25). In principle, a molecular dynamic simulation of the free... 

Significant progress has recently been made in several areas which will have a profound effect on the ability of molecular dynamics to handle more complex problems. In this section we speculate on several areas which appear to hold promise for advancing computer modeling studies. In section 4.1, recent progress in both analytic potential energy expressions and first principles calculations are briefly mentioned. Recent advances in computational techniques are discussed in section 4.2. These include the use of constraints within the classical equations of motion to model thermostats in the surface region, and the incorporation of Monte Carlo techniques into molecular dynamics simulations. 

The results shown in this study are limited to the KMC algorithm. In principle, to model a realistic system the set k, can be found using molecular dynamics simulations, or other similar techniques these can then be used as input into KMC. Here, our purposes are more general. A square lattice is examined, and there are two simplified rate laws of interest. The first is often used in KMC simulations of deposition, and is termed i-kinetics. A second set represents a kinetics that is analogous to the hopping probabilities used in the Metropolis simulations. We term the latter A/-kinetics. 

Rovira, C., and Parrinello, M. 2000. First-principles molecular dynamics simulations of models of the myoglobin active center. Int. J. Quantum Chem. 80 1172-80. 

The credit load for die computational chemistry laboratory course requires that the average student should be able to complete almost all of the work required for the course within die time constraint of one four-hour laboratory period per week. This constraint limits the material covered in the course. Four principal computational methods have been identified as being of primary importance in the practice of chemistry and thus in the education of chemistry students (1) Monte Carlo Methods, (2) Molecular Mechanics Methods, (3) Molecular Dynamics Simulations, and (4) Quantum Chemical Calculations. Clearly, other important topics could be added when time permits. These four methods are developed as separate units, in each case beginning with die fundamental principles including simple programming and visualization, and building to the sophisticated application of the technique to a chemical problem. 

---