Molecular dynamics simulation calculation techniques

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. 

The rest of this paper will be organized as follows In the next section we describe a set of interatomic potentials which we have developed to model iron-nickel alloys and report the calculational techniques of molecular dynamics simulations which we... 

In our calculations we make use of several standard techniques of molecular dynamics simulations. The integration of the equations of motions is done by the velocity form of the Verlet-algorithm with a time step of 1.5 The temperature is controlled... 

Molecular dynamics simulations, with quantum-mechanically derived energy and forces, can provide valuable insights into the dynamics and structure of systems in which electronic excitations or bond breaking processes are important. In these cases, conventional techniques with classical analytical potentials, are not appropriate. Since the quantum mechanical calculation has to be performed many times, one at each time step, the choice of a computationally fast method is crucial. Moreover, the method should be able to simulate electronic excitations and breaking or forming of bonds, in order to provide a proper treatment of those properties for which classical potentials fail. 

An interesting combined use of discrete molecular and continuum techniques was demonstrated by Floris et al.181,182 They used the PCM to develop effective pair potentials and then applied these to molecular dynamics simulations of metal ion hydration. Another approach to such systems is to do an ab initio cluster calculation for the first hydration shell, which would typically involve four to eight water molecules, and then to depict the remainder of the solvent as a continuum. This was done by Sanchez Marcos et al. for a group of five cations 183 the continuum model was that developed by Rivail, Rinaldi et al.14,108-112 (Section III.2.ii). Their results are compared in Table 14 with those of Floris et al.,139 who used a similar procedure but PCM-based. In... 

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. 

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. 

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. 

What is next Several examples were given of modem experimental electrochemical techniques used to characterize electrode-electrolyte interactions. However, we did not mention theoretical methods used for the same purpose. Computer simulations of the dynamic processes occurring in the double layer are found abundantly in the literature of electrochemistry. Examples of topics explored in this area are investigation of lateral adsorbate-adsorbate interactions by the formulation of lattice-gas models and their solution by analytical and numerical techniques (Monte Carlo simulations) [Fig. 6.107(a)] determination of potential-energy curves for metal-ion and lateral-lateral interaction by quantum-chemical studies [Fig. 6.107(b)] and calculation of the electrostatic field and potential drop across an electric double layer by molecular dynamic simulations [Fig. 6.107(c)]. 

The development of multiscale simulation techniques that involve the atomistic modeling of various structures and processes still remains at its early stage. There are many problems to be solved associated with more accurate and detailed description of these structures and processes. These problems include the development of efficient and fast methods for quantum calculations at the atomistic level, the development of transferable interatomic potentials (especially, reactive potentials) for molecular dynamic simulations, and the development of strategies for the application of multiscale simulation methods to other important processes and materials (optical, magnetic, sensing, etc.). 

---