Application of the Molecular Dynamics Method

The most successful and complete study of the behavior of a sample of waterlike particles in three dimensions has been carried out by Rahman 

There is a wealth of other information, on both equilibrium and kinetic properties, of this sample of waterlike particles obtained by the molecular dynamics computations which are not reproduced here. The interested reader is referred to the original articles for further details. 

We have described two main lines of development in the theory of liquid water. The first, and the older, was founded on the mixture-model approach (Chapter 5) to liquids, which offers certain approximate or ad hoc models for the fluids as a whole. The second approach may be referred to as the ab initio method, based on first principles of statistical mechanics. In the past, these two lines of development were thought to be conflicting and a vigorous debate has taken place on this issue. As we have stressed throughout this and the previous chapter, both approaches can be developed from first principles, and, in fact, provide complementary information on this liquid. Once we attempt to pursue this theory along either route, we must introduce serious approximations. Therefore, it is very difficult to establish a clear-cut preference for one approach or the other. 

Another general comment regarding the two approaches is their ability to explain the anomalous behavior of liquid water and its solutions. We 

---

Both molecular dynamics studies and femtosecond laser spectroscopy results show that molecules with a sufficient amount of energy to react often vibrate until the nuclei follow a path that leads to the reaction coordinate. Dynamical calculations, called trajectory calculations, are an application of the molecular dynamics method that can be performed at semiempirical or ah initio levels of theory. See Chapter 19 for further details. 

Lantelme, R, Turq, R, Quentrec, B., and Lewis, J.W.E., Application of the molecular dynamics method to a liquid system with long range forces (Molten NaCl), Mol. Phys., 28,1537-1549, 1974. 

In this and the following three sections we discuss the principal types of application of the molecular dynamics method. We begin with the problem of the approach to equilibrium. 

An application of the molecular dynamics method to simulate the liquid-vapor surface of molecular fluids is described. A predictor-corrector algorithm is used to solve the equations of translational and rotational motion, where the orientations of molecules are expressed in quaternions. The method is illustrated with simulations of 216 homonuclear (N2 and Clz) diatomic molecules. Properties calculated include surface tensions and density-orientation profiles. 

The essence of the molecular dynamics method is a straightforward application of Newton s second law. To obtain the time evolution of the system, we integrate the set of coupled second order differential equations... 

In tfiis chapter we address first the electrochemical application of the more familiar method of molecular (or atom) dynamics, and later turn to consider Monte Carlo methods, in each case giving a short introduction that should motivate the reader to pursue reading more specific works. Although the present research field is relatively new, the investigations are already too extensive to review in detail in a single chapter. For this reason, we discuss here the more extended research branches in the field and present a few representative examples. The application of simulations applied to nanostructuring problems is discussed in Chapter 36 liquid-liquid interfaces have been addressed by I. Benjamin (1997). 

Some of the major areas of activity in this field have been the application of the method to more complex materials, molecular dynamics, [28] and the treatment of excited states. [29] We will deal with some of the new materials in the next section. Two major goals of the molecular dynamics calculations are to determine crystal structures from first principles and to include finite temperature effects. By combining molecular dynamics techniques and ah initio pseudopotentials within the local density approximation, it becomes possible to consider complex, large, and disordered solids. 

Recent rapid developments in ultrashort pulse laser [1-5] make it possible to probe not only the dynamics of population of the system but also the coherence (or phase) of the system. To treat these problems, the density matrix method is an ideal approach. The main purpose of this paper is to briefly describe the application of the density matrix method in molecular terms and show how to apply it to study the photochemistry and photophysics [6-9]. Ultrafast radiationless transactions taking place in bacterial photosynthetic reaction centers (RCs) are very important examples to which the proposed theoretical approach can be applied. 

In this article, we shall discuss studies of phase changes in solids carried out by the application of the generalized Monte Carlo and the molecular dynamics methods. This topic is of particular significance because of the recent modifications of the method to include both variation in size and shape of the simulation cell. What is especially gratifying is that we are now able to make meaningful predictions of phase transitions in real solids by employing reliable pair potentials. Besides phase transitions of molecular solids, we shall examine phase transi-... 

In these past 10 years, it has been demonstrated that the TR-QELS method is a versatile technique that can provide much information on interfacial molecular dynamics [1-11]. In this chapter, we intend to show interfacial behaviour of molecules elucidated by the TR-QELS method. In Section 3.2, we present the principle, the historical background and the experimental apparatus for TR-QELS. The dynamic collective behaviour of molecules at liquid/liquid interfaces was first obtained by improving the time resolution of the TR-QELS method. In Section 3.3, we present an application of the TR-QELS method to a phase transfer catalyst system and describe results on the scheme of the catalytic reactions. This is the first application of the TR-QELS method to a practical liquid/liquid interface system. In Section 3.4, we show chemical oscillations of interfacial tension and interfacial electric potential. In this way, the TR-QELS method allows us to analyze non-linear adsorption/desorption behaviour of surfactant molecules in the system. 

The first structure of human renin was obtained from prorenin produced by expression of its cDNA in transfected mammalian cells. Prorenin was cleaved in the laboratory to renin using the protease trypsin. Because the carbohydrates in renin are not required for bioactivity, oligosaccharides were removed enzymatically. This process facilitates crystallization in some cases and also removes the contribution of the heterogeneous sugar chains to the diffraction pattern. The structure was determined without the use of heavy-atom derivatives, by application of molecular replacement techniques based on the atomic coordinates of porcine pepsinogen as the model. The molecular dynamic method of refinement was used extensively to arrive at a 2.5 A resolution structure. However, some of the loop regions were not well resolved in this structure (Sielecki et al, 1989 Sail et al, 1990). 

Pak Y, Wang S. Application of a molecular dynamics simulation method with a generalized effective potential to the flexible molecular docking problems. J Phys Chem B 2000 104 354-359. 

The molecular dynamics methods that we have discussed in this chapter, and the examples that have been used to illustrate them, fall into the category of atomistic simulations, in that all of the actual atoms (or at least the non-hydrogen atoms) in the core system are represented explicitly. Atomistic simulations can provide very detailed information about the behaviour of the system, but as we have discussed this typically limits a simulation to the nanosecond timescale. Many processes of interest occur over a longer timescale. In the case of processes which occur on a macroscopic timescale (i.e. of the order of seconds) then rather simple models may often be applicable. Between these two extremes are phenomena that occur on an intermediate scale (of the order of microseconds). This is the realm of the mesoscale Dissipative particle dynamics (DPD) is particularly useful in this region, examples include complex fluids such as surfactants and polymer melts. 

Increases in computer power and improvements in algorithms have greatly extended the range of applicability of classical molecular simulation methods. In addition, the recent development of Internal Coordinate Quantum Monte Carlo (ICQMC) has allowed the direct comparison of classical simulations and quantum mechanical results for some systems. In particular, it has provided new insights into the zero point energy problem in many body systems. Classical studies of non-linear dynamics and chaos will be compared to ICQMC results for several systems of interest to nanotechnology applications. The ramifications of these studies for nanotechnology applications will be discussed. 

The second line of inquiry alluded to was the use of modem computational power, both hardware and software, for the evaluation of pair distribution functions. When the Kirkwood-Buff paper was published, the use of computers for this kind of scientific computation was in its infancy. Indeed, one can say that it was in its prenatal stage. It is difficult to put a date on the time when computers became powerful enough to compute pair correlation functions and, consequently, KB integrals with sufficient accuracy for application to real systems. They have certainly reached that stage at the time of the writing of these words. The computational method of choice in carrying out these calculations is the molecular dynamics method. Since this kind of calculation is discussed in detail in several of the later chapters of this work, we eschew discussion here. 

In this chapter, we concentrate on luminescence approaches that are most suitable for stilbenes to be molecular probes and are based on specific and nonspecific labeling, competition, solvatochromism, experimental molecular dynamics, and singlet-singlet and triplet-triplet energy transfer. A general survey is made of the physical principles and application of the fluorescence probe methods stressing on latest developments in this area. 

Analyses of the data generated by the molecular dynamics method is an important part of the overall simulation. We have developed a number of very useful tools for the examination of the results such as to extract the maximum information. These include high resolution spectral estimators and the application of neural/fuzzy systems. A more detailed description of the methods and examples of application can be found in a recent review article [11]. 

Over the past 10-15 years computational methods have been developed which permit the study of protein and nucleic acid motions and structure, as well as some aspects of their reactivity. These techniques, known as biopolymer dynamics and mechanics [1,2], evolved from pioneering work by Alder and Wainwright [3] and Rahman [4] on the classical simulation of condensed phase systems. They were solidified by the first application of classical molecular dynamics to proteins by McCammon, Gelin and Karplus in 1977 [5]. Today a broad range of biophysical processes are explored using molecular simulation methods [1, 2]. 

---