Computational methods molecular dynamics

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

During the last few years the progress of computational techniques has made it possible to simulate the dynamic behavior of whole ensembles consisting of several hundred molecules. In this way the limitations of the statistical approach can be at least partly overcome. Two kinds of methods — molecular dynamics and Monte Carlo calculations — were applied to liquids and liquid mixtures and brought new insight into their structure and properties. Even some important characteristics of systems as complicated as associated liquids like water could be... 

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

An overview of the approaches that have been taken to linking different theoretical and computational modeling descriptions is also provided in Fig. 2 The first principles (QC) descriptions are based on the Schrodinger equation and the Bom Oppenheimer approximation as realized in most chemical applications by density functional [14] or Hartree-Fock [15] methods. Molecular dynamics (MD) methods [16], based on classical Newtonian me-... 

Minoru Tanaka and Ryoichi Yamamoto, Computer Physics and Computer Chemistry Molecular Dynamics Method and Monte Carlo Method, Kaibundo Publishing, Tokyo, 1988. 

The basic idea behind an atomistic-level simulation is quite simple. Given an accurate description of the energetic interactions between a collection of atoms and a set of initial atomic coordinates (and in some cases, velocities), the positions (velocities) of these atoms are advanced subject to a set of thermodynamic constraints. If the positions are advanced stochastically, we call the simulation method Monte Carlo or MG [10]. No velocities are required for this technique. If the positions and velocities are advanced deterministically, we call the method molecular dynamics or MD [10]. Other methods exist which are part stochastic and part deterministic, but we need not concern ourselves with these details here. The important point is that statistical mechanics teUs us that the collection of atomic positions that are obtained from such a simulation, subject to certain conditions, is enough to enable aU of the thermophysical properties of the system to be determined. If the velocities are also available (as in an MD simulation), then time-dependent properties may also be computed. If done properly, the numerical method that generates the trajectories... 

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

The generation of El mass spectra is well understood with a rich documentation of the fragmentation processes. But in contrast to NMR, mass spectra still cannot be calculated from a theoretical standpoint. A general method to compute El mass spectra based on a combination of fast quantum chemical methods, molecular dynamics and stochastic preparation of hot ionized species was published by Grimme. It provides mass spectra that compare well with their experimental counterparts, even in subtle details (Grimme, 2013). 

The methodology for microscale heat transfer study can be classified into three categories. The first method attempts to modify the continuum model in such a way that microscale consideration is taken into account. The second method is application of the Boltzmann transport equation (BTE). The third approach is computationally exhaustive molecular dynamics approach, which is typically used when the first two methods fail. The present chapter discusses the introductory concepts of microscale conduction. 

Particle methods (Molecular Dynamics, Dissipative Particle Dynamics, Multi-Particle Collision Dynamics) simulate a system of interacting mass points, and therefore thermal fluctuations are always present. The particles may have size and structure or they may be just point particles. In the former case, the finite solvent size results in an additional potential of mean force between the beads. The solvent structure extends over unphysically large length scales, because the proper separation of scale between solute and solvent is not computationally realizable. In dynamic simulations of systems in thermal equilibrium [43], solvent structure requires that the system be equilibrated with the solvent in place, whereas for a structureless solvent the solute system can be equilibrated by itself, with substantial computational savings [43]. Finally, lattice models have a (rigorously) known solvent viscosity, whereas for particle methods the existing analytical expressions are only approximations (which however usually work quite well). 

---