Computer simulation molecular dynamics method

The principal advantage of the time correlation function method is that it provides a new set of microscopic functions for a fluid, the time correlation functions, which can be studied directly by experimental observations of the fluidt or by computer-simulated molecular dynamics. The time correlation functions depend even more sensitively on the microscopic properties of the fluid molecules than the transport coefficients, which are expressed as time integrals of the correlation functions. Thus, a further test of kinetic theory has been found it must not only lead to expressions for the transport coefficients for dilute and dense gases that are in agreement with experiment, but also describe the dependence of the time correlation functions on both time and the density of the gas. One of the principal successes of kinetic theory is that it provides a quantitatively correct description of the short- and long-time... 

As we discussed earlier, the generalized Boltzmann equation leads to a density expansion of the transport coefficients of a dense gas. However, general expressions for transport coefficients of a fluid that are not in the form of an expansion can be derived by another technique, the time correlation function method. This approach has provided a general framework by means of which one can make detailed comparisons between theoretical results, the results of computer-simulated molecular dynamics,and experimental results. ... 

Being a d5mamic property, the diffusion coefficient can only be computed by molecular dynamics methods and not by Monte Carlo techniques. As we shall see, the difficult part for an ionic liquid is that the long enough simulation time can be of the order of 10-100 ns which is hardly attainable for current hardware and sophisticated molecular entities. Diffusive regimes where... 

Sometimes the theoretical or computational approach to description of molecular structure, properties, and reactivity cannot be based on deterministic equations that can be solved by analytical or computational methods. The properties of a molecule or assembly of molecules may be known or describable only in a statistical sense. Molecules and assemblies of molecules exist in distributions of configuration, composition, momentum, and energy. Sometimes, this statistical character is best captured and studied by computer experiments molecular dynamics, Brownian dynamics, Stokesian dynamics, and Monte Carlo methods. Interaction potentials based on quantum mechanics, classical particle mechanics, continuum mechanics, or empiricism are specified and the evolution of the system is then followed in time by simulation of motions resulting from these direct... 

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... 

The computational description of a large biomolecular complex such as the chromatin fiber requires techniques that are difierent from the widely applied molecular dynamics methods used to simulate biopolymers at atomic resolution. 

Interestingly, in the experiments devoted solely to computational chemistry, molecular dynamics calculations had the highest representation (96-98). The method was used in simulations of simple liquids, (96), in simulations of chemical reactions (97), and in studies of molecular clusters (98). One experiment was devoted to the use of Monte Carlo methods to distinguish between first and second-order kinetic rate laws (99). One experiment used DFT theory to study two isomerization reactions (100). 

The method presented in this chapter serves as a link between molecular properties (e.g., cavities and their occupants as measured by diffraction and spectroscopy) and macroscopic properties (e.g., pressure, temperature, and density as measured by pressure guages, thermocouples, etc.) As such Section 5.3 includes a brief overview of molecular simulation [molecular dynamics (MD) and Monte Carlo (MC)] methods which enable calculation of macroscopic properties from microscopic parameters. Chapter 2 indicated some results of such methods for structural properties. In Section 5.3 molecular simulation is shown to predict qualitative trends (and in a few cases quantitative trends) in thermodynamic properties. Quantitative simulation of kinetic phenomena such as nucleation, while tenable in principle, is prevented by the capacity and speed of current computers however, trends may be observed. 

We do not enter into the description of Monte Carlo of Molecular Dynamics methods, as these details are not important for our discussion. There are other more general aspects of computer simulations to consider here. 

Computer simulation is being used increasingly in diverse areas of science in the past few years. It has also emerged to become one of the powerful means for investigating condensed matter (/). The principal tools employed in computer simulation are the Monte Carlo and the molecular dynamics methods. In these methods, properties of a collection of particles, usually between 30 and 1000 in number, interacting via a potential />(r) are obtained numerically. Reliable estimates of equilibrium and transport properties as well as microscopic properties can be obtained from such calculations. 

The problems being addressed in recent work carried out in various laboratories include the fundamental nature of the solute-water intermolecular forces, the aqueous hydration of biological molecules, the effect of solvent on biomolecular conformational equilibria, the effect of biomolecule - water interactions on the dynamics of the waters of hydration, and the effect of desolvation on biomolecular association 17]. The advent of present generation computers have allowed the study of the structure and statistical thermodynamics of the solute in these systems at new levels of rigor. Two methods of computer simulation have been used to achieve this fundamental level of inquiry, the Monte Carlo and the molecular dynamics methods. 

Recent single-molecule experimental studies of proteins provide more detailed views of protein motions, and confirm that a wide variety of timescales is involved in, e.g., catalytic action of enzymes [7,14,15,19,33], Of course, molecular dynamics simulations have been used to probe motions in single proteins for many years, and advances in both theory and computational science have made simulations a powerful approach to building theoretical understanding of protein dynamics [1], The recent introduction of accelerated molecular dynamics methods is helpful in this context [11]. Although detailed dynamical information is sacrificed to the enhanced sampling of conformational space in these methods, which have been shown to access conformational fluctuations that are revealed by nuclear magnetic resonance experiments on the millisecond... 

Computer simulations (Mmite Carlo method and molecular dynamics calculations have been used to analyze the stmcture of water. According to these studies, water is made up of a continuous distribution of associates. The pair interaction between nearest neighbor molecules should not be thought of as either... 

Recent years have seen the extensive application of computer simulation techniques to the study of condensed phases of matter. The two techniques of major importance are the Monte Carlo method and the method of molecular dynamics. Monte Carlo methods are ways of evaluating the partition function of a many-particle system through sampling the multidimensional integral that defines it, and can be used only for the study of equilibrium quantities such as thermodynamic properties and average local structure. Molecular dynamics methods solve Newton s classical equations of motion for a system of particles placed in a box with periodic boundary conditions, and can be used to study both equilibrium and nonequilibrium properties such as time correlation functions. 

An appealing approach to the study of nucleation is to observe it directly in a computer simulation using the method of molecular dynamics. It is evident that one cannot closely mimic experimental conditions, since computer time scales extend only over tens to hundreds of picoseconds. Thus deep quenches carried out at high quench rates are necessary to form supercooled liquids that have some reasonable chance of nucleating. Under such conditions glass formation is also observed, and in fact the first observation of nucleation was a chance event, although since that time more systematic studies have been carried out. 

Finite temperature effects are accessible by ab initio molecular dynamics methods, but the simulation times yet achieved are still at least an order of magnitude too small to give accurate numbers [36]. Again, particularly the development of massive parallel computer hardware and algorithms will change this situation in the near future. 

---