Molecular dynamics type simulations

When addressing problems in computational chemistry, the choice of computational scheme depends on the applicability of the method (i.e. the types of atoms and/or molecules, and the type of property, that can be treated satisfactorily) and the size of the system to be investigated. In biochemical applications the method of choice - if we are interested in the dynamics and effects of temperature on an entire protein with, say, 10,000 atoms - will be to run a classical molecular dynamics (MD) simulation. The key problem then becomes that of choosing a relevant force field in which the different atomic interactions are described. If, on the other hand, we are interested in electronic and/or spectroscopic properties or explicit bond breaking and bond formation in an enzymatic active site, we must resort to a quantum chemical methodology in which electrons are treated explicitly. These phenomena are usually highly localized, and thus only involve a small number of chemical groups compared with the complete macromolecule. 

We have reviewed above the GH approach to reaction rate constants in solution, together with simple models that give a deeper perspective on the reaction dynamics and various aspects of the generalized frictional influence on the rates. The fact that the theory has always been found to agree with Molecular Dynamics computer simulation results for realistic models of many and varied reaction types gives confidence that it may be used to analyze real experimental results. 

The microscopic structure of water at the solution/metal interface has been the focus of a large body of literature, and excellent reviews have been published summarizing the extensive knowledge gained from experiments, statistical mechanical theories of varied sophistication, and Monte Carlo and molecular dynamics computer simulations. To keep this chapter to a reasonable size, we limit ourselves to a brief summary of the main results and to a sample of the type of information that can be gained from computer simulations. 

Following the early studies on the pure interface, chemical and electrochemical processes at the interface between two immiscible liquids have been studied using the molecular dynamics method. The most important processes for electrochemical research involve charge transfer reactions. Molecular dynamics computer simulations have been used to study the rate and the mechanism of ion transfer across the water/1,2-dichloroethane interface and of ion transfer across a simple model of a liquid/liquid interface, where a direct comparison of the rate with the prediction of simple diffusion models has been made. ° ° Charge transfer of several types has also been studied, including the calculations of free energy curves for electron transfer reactions at a model liquid/liquid... 

Molecular dynamics (MD) simulations are a class of molecular mechanics calculation which directly model the motions of molecular systems, often providing considerable information which cannot be obtained by any other technique, theoretical or experimental. MD simulations have only recently been applied to problems of carbohydrate conformation and motions, but it is likely that this technique will be widely used for modeling carbohydrates in the future. This paper introduces the basic techniques of MD simulations and illustrates the types of information which can be gained from such simulations by discussing the results of several simulations of sugars. The importance of solvation in carbohydrate systems will also be discussed, and procedures for including solvation in molecular dynamics simulations will be introduced and again illustrated from carbohydrate studies. 

It should be mentionned that the result of q-dependent chain orientation in shear was also found by nonequilibrium molecular dynamics (NEMD) simulations by Kroger et al [32], The increasing power of computational techniques will surely result in increased accuracy and usefulness of this type of numerical simulation. 

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. 

Complementing the information conveyed by Ty, a fugacity f is introduced and defined as f = Nj/Nfib), where Nj is the expected number of water molecules in Dj at equilibrium and Nj(b), the number associated with the same volume in bulk solvent. Thus, the chemical potential p,j of water in Dj becomes prj = k Tln [A//A/(b)] ( b= Boltzmann constant). The 0-dependence of T, /-values (Fig. 14.1) is obtained at equilibrium determined from classical trajectories generated by molecular dynamics (MD) simulations within an NPT ensemble of the type described in Chap. 4. The computations start with the protein... 

We note that in molecular-dynamics (MD) simulations we make no approximations other than the ones implied in the interatomic potentials and the fact that the d5mamics of the atoms is purely classical (no quantum effects on the atomie motion). For example, no approximation is made sis to what type of ehemical reaction can or can not occur complex phenomena such as pressure effects, multi-molecular reactions, and relaxation are explicitly described in NEMD. In this sense, the simulations presented here provide a full-physics, full-chemistry description of energetic materials. 

The microkinetic interpretation of the origin of these different patterns may be confirmed by molecular dynamics (MD) simulations. As an example. Fig. 11 presents the result of an MD study with methane in a cation-free zeolite of type LTA [117,118]. By increasing the Lennard-Jones distance a between the methane molecules and the oxygen of the zeolite lattice one is able to simulate the influence of a reduction of the window diameter on the dif-... 

Figure 1 Schematic drawing of the various types of surfactants used in the molecular dynamics (MD) simulations discussed in this chapter. The lighter spheres represent solvophilic beads (solvent bead or head bead), and the darker spheres are the solvophobic beads (tail bead or oil bead).

---