Molecular dynamics physical simulation

To examine the soUd as it approaches equUibrium (atom energies of 0.025 eV) requires molecular dynamic simulations. Molecular dynamic (MD) simulations foUow the spatial and temporal evolution of atoms in a cascade as the atoms regain thermal equiUbrium in about 10 ps. By use of MD, one can foUow the physical and chemical effects that induence the final cascade state. Molecular dynamics have been used to study a variety of cascade phenomena. These include defect evolution, recombination dynamics, Hquid-like core effects, and final defect states. MD programs have also been used to model sputtering processes. 

Evans and Baranyai [51, 52] have explored what they describe as a nonlinear generalization of Prigogine s principle of minimum entropy production. In their theory the rate of (first) entropy production is equated to the rate of phase space compression. Since phase space is incompressible under Hamilton s equations of motion, which all real systems obey, the compression of phase space that occurs in nonequilibrium molecular dynamics (NEMD) simulations is purely an artifact of the non-Hamiltonian equations of motion that arise in implementing the Evans-Hoover thermostat [53, 54]. (See Section VIIIC for a critical discussion of the NEMD method.) While the NEMD method is a valid simulation approach in the linear regime, the phase space compression induced by the thermostat awaits physical interpretation even if it does turn out to be related to the rate of first entropy production, then the hurdle posed by Question (3) remains to be surmounted. 

Molecular predictions of the properties of interfacial systems are now becoming possible as a result of rapid advances in liquid state chemical physics and computer technology. The objectives of this paper are 1) to review the general approaches and models used in Monte Carlo (MC) and molecular dynamics (MD) simulations of interfacial systems, 2) to describe and discuss results from selected simulation studies of interfacial water, and 3) to discuss the major limitations of these techniques and to offer suggestions for overcoming them. 

One of the major objectives of the physical chemistry studies in water and biomolecules is to fully reproduce the experimentally observed folding/ unfolding behavior of a typical model protein in water by means of molecular simulation. However, the all-atom molecular dynamics (MD) simulation of the folding of a protein from the fully unfolded state to the native structure remains computationally intractable when the size of the target protein is larger than 100 residues and when simulation is carried out with explicit water molecules (i.e., when complete, contextualized simulation is attempted) [1-3]. 

Dynamical simulated annealing (DSA)177 is a variant of restrained molecular dynamics (RMD).178 There are numerous programs available for performing molecular dynamics (MD) simulations, including GROMOS,178 AMBER,179 CHARMM,180 X-PLOR/CNS,181 and OPLS.182 In MD simulations, Newton s equations of motion are solved for all atoms under the influence of a physical force field ( physical), which for a protein has the form183... 

The topic of this article is the use of molecular dynamics (MD) simulations of positive ion-surface interactions for insights into the chemical and physical processes that occur at surfaces immersed in glow discharge plasmas. To understand the signihcance of this topic, it is necessary to have some background in the technology and its major industrial application (Lieberman and Lichtenberg, 1994). The term plasma in this context refers... 

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. 

An alternative to molecular dynamics based simulated annealing is provided by Metropolis importance sampling Monte Carlo (Metropolis et al., 1953) which has been widely exploited in the evaluation of configurational integrals (Ciccotti et al., 1987) and in simulations of the physical properties of liquids and solids (Allen and Tildesley, 1987). Here, as outlined in Chapters 1 and 2, a particle or variable is selected at random and displaced both the direction and magnitude of the applied displacement within standard bounds are randomly selected. The energy of this new state, new, is evaluated and the state accepted if it satisfies either of the following criteria ... 

Gain G and Pasquarello A 1993 First-principles molecular dynamics Computer Simulation in Chemical Physics vol 397 NATO ASI Series C ed M P Allen and D J Tildesley (Dordrecht Kluwer) pp 261-313... 

Molecular dynamics (MD) simulations predict the movement of atoms and molecules with time under some basic laws of physics. Normally, real molecular systems have a larger number of freedoms, because a huge amount of particles are included. It is impractical to solve the movement equations and hnd the properties of such complex systems analytically. Assisted by the development of computational technique, MD simulations overcome the problem using numerical methods, and build an interface between laboratory experiments and theory. 

While molecular dynamics (MD) simulations have proven to be very powerful for studying numerous aspects of protein dynamics and structure [11-13], this technique cannot yet access the millisecond-to-second time-scales required for folding even a small protein. To address this timescale gap, one has to simplify the protein model by reducing the number of degrees of freedom. Such approaches assume that the basic physics could be reproduced in model systems that employ united atoms and effective solvent models. On the basis of recent work, it has become apparent that the crux of the solution to the protein folding problem does not lie in whether a reduced protein model is used, but rather in the development of... 

The advent of high-speed computers, availability of sophisticated algorithms and state-of-the-art computer graphics have made plausible the use of computationally intensive methods such as QM, MM and molecular dynamics (MolD) simulations to determine those physical and structural properties most commonly involved in molecular processes. The power of molecular modeling rests solidly on a variety of well-established scientific disciplines including computer science, theoretical chanistry, biochemistry and biophysics. Molecular modeling has become an indispensable complementary tool for most experimental scientific research. 

A distinctly different approach, which has witnessed much progress recently, is large scale Monte Carlo and molecular dynamics computer simulations [4]. These studies provide many insights regarding the physics of model polymer fluids, and also valuable benchmarks against which approximate theory can be tested. However, an atomistic, off-lattice treatment of high polymer fluids and alloys remains immensely expensive, if not impossible, from a computational point of view. 

Abstract We present in this contribution results from Molecular Dynamics (MD) simulations of a chemically realistic model of 1,4-polybutadiene (PB). The work we will discuss exemplifies the physical questions one can address with these types of simulations. We will specifically compare the results of the computer simulations with nuclear magnetic resonance (NMR) experiments, neutron scattering experiments and dielectric data. These comparisons will show how important it is to understand the torsional dynamics of polymers in the melt to be able to explain the experimental findings. We will then introduce a freely rotating chadn (FRC) model where all torsion potentials have been switched off and show the influence of this procedure on the qualitative properties of local dynamics through comparison with the chemically realistic (CRC) model. 

Because of the very large surface-area-to-volume ratios of micro-devices, adhe-sion/stiction has been considered the most important failure mode and the major obstacle for the commercialization of micro-electromechanical systems (MEMS). In this chapter, most important surface forces are introduced. The physical origin and mathematical models of these surface forces are presented. Then, adhesion effects such as wetting and surface energy, which are related to these surface forces, are extensively discussed. Self-assembled monolayers (SAMs) have recently received considerable attention as molecular-level lubricants in MEMS. The structure and the surface characteristics of SAMs are introduced. Experiments, molecular dynamics (MD) simulations, and theoretical models on the adhesion force between the atomic force microscope (AFM) tip and sample are discussed in detail. Finally, the adhesion problems related to super-hydrophobic films are discussed. 

FIGURE 9.5 Comparison of the equation of state (reduced axial pressure versus reduced numerical density) of Square-Well molecules of = 1.5 confined in cylindrical hard pore with diameter, D/a = 2.2, obtained by isobaric-isothermal Monte Carlo (NPT MC) and molecular dynamic (MD) simulations. Here, squares indicate NPT MC result and circles the MD result. The solid line indicates an analytical fit of the result at the fluid branch, and the dash line is the second order polynomial fit to the solid branch. Error bars are the standard deviation of five independent runs. (From Huang, H. C., J. Chem. Phys., 132, 224504, 2010. With permission. Copyright 2010, American Institute of Physics.)... 

In this section we describe quasielastic neutron scattering studies [129,130] focusing attentions on conformational transitions of polymer chains. For this purpose we first summarize the results of the recent molecular dynamics (MD) simulations on conformational transitions and then discuss the conformational transition mechanism on the basis of neutron data analyzed in terms of a jump diffusion model with damped vibration which has a similar physical picture to that predicted by the MD simulations. 

In spite of their many assets, DPD and SCMF methods are not able to accurately predict physical properties that rely upon time correlation functions, for example, diffusion coefficients. The more feasible mesoscale approach for hydrated ionomer membranes is coarse-grained molecular dynamics (CGMD) simulations. 

---