Equilibrium, molecular dynamics system

Equilibrium Systems. Magda et al (12.) have carried out an equilibrium molecular dynamics (MD) simulation on a 6-12 Lennard-Jones fluid In a silt pore described by Equation 41 with 6 = 1 with fluid particle Interactions given by Equation 42. They used the Monte Carlo results of Snook and van Me gen to set the mean pore density so that the chemical potential was the same In all the simulations. The parameters and conditions set In this work were = 27T , = a, r = 3.5a, kT/e = 1.2, and... 

The static and dynamic properties of polymer-layered silicate nanocomposites are discussed, in the context of polymers in confined spaces and polymer brushes. A wide range of experimental techniques as applied to these systems are reviewed, and the salient results from these are compared with a mean field thermodynamic model and non-equilibrium molecular dynamics simulations. 

Transport coefficients of molecular model systems can be calculated by two methods [8] Equilibrium Green-Kubo (GK) methods where one evaluates the GK-relation for the transport coefficient in question by performing an equilibrium molecular dynamics (EMD) simulation and Nonequilibrium molecular dynamics (NEMD) methods. In the latter case one couples the system to a fictitious mechanical field. The algebraical expression for the field is chosen in such a way that the currents driven by the field are the same as the currents driven by real Navier-Stokes forces such as temperature gradients, chemical potential gradients or velocity gradients. By applying linear response theory one can prove that the zero field limit of the ratio of the current and the field is equal to the transport coefficient in question. 

In a somewhat similar paper, diffusion through a 2D porous solid modeled by a regular array of hard disks was evaluated [65] using non-equilibrium molecular dynamics. It was found that Pick s law is not obeyed in this system unless one takes different diffusion constants for different regions in the flow system. Other non-equilibrium molecular dynamics simulations of diffusion for gases within a membrane have been presented [66]. The membrane was modeled as a randomly... 

One easy extension of equilibrium molecular dynamics is the computation of the properties of condensed matter systems in the presence of external fields. Many experiments are done by applying an external field, whether it be an electric field, heat gradient, or some type of flow. Usually, the experimentalist waits some time for the system under investigation to reach a steady state in the presence of the applied field. Measurements are then performed to deduce structural or dynamical information. Performing this task on a computer is what is commonly known as nonequilibrium molecular dynamics (NEMD). It is nonequilibrium in the sense that, in the presence of an external field, the system at the steady state will be in a state of lower entropy. Upon removing the external field, the system will return to the state of maximum entropy or the equilibrium state. 

Equilibrium molecular dynamics was put on a firm theoretical ground with Andersen s seminal paper on extended system dynamics (to be discussed in great detail later). NEMD found its first success in this area with the advent of the so-called DOLL s (not an acronym) algorithm by Hoover and coworkers. ... 

An immediate use for this conserved quantity is obvious it can (and should) be used to check the NEMD code for algorithmic and programming errors. It is also possible to use the conserved energy in obtaining a knowledge of the phase space. The approach proceeds in the same fashion as presented in the section on equilibrium molecular dynamics. Let F denote the full phase space of the variables, p , q,, ri,, I. We now make the assumption of equal a priori probability for each of the microstates F with energy H. This assumption has traditionally been applied to equilibrium systems only. In the isolated system we consider, this assumption is the most obvious one to make. Thus, one can write the phase space distribution function /(F) as... 

The use of non-Hamiltonian dynamical systems has a long history in mechanics [8] and they have recently been used to study a wide variety of problems in molecular dynamics (MD). In equilibrium molecular dynamics we can exploit non-Hamiltonian systems in order to generate statistical ensembles other than the standard microcanonical ensemble NVE) that is generated by traditional Hamiltonian dynamics. These ensembles, such as the canonical (NVT) and isothermal-isobaric (NPT) ensembles, are much better than the microcanonical ensemble for representing the actual conditions under which experiments are carried out. 

The earliest efforts to develop non-equilibrium molecular dynamics (NEMD) methods used special boundary conditions and/or external fields to induce non-equilibrium behavior in the system. Important contributions to this development include those of Lees and Edwards [72], Gosling et al. [73], Hoover and Ashurst [74] and Ciccotti and Jacucci... 

The same effect in Zeolite-guest systems was demonstrated by Auerbach [26] by equilibrium molecular dynamics and non-equilibrium molecular dynamics after experimental work by Cormer [27]. The energy distributions obtained in Zeolite and Zeolite-Na are shown in Fig. 5.6A. At equilibrium, all the atoms in the system are at the same temperature. When Na-Y Zeolite is exposed to microwave energy, however, the effective steady-state temperature of Na atoms is substantially higher than that of the rest of the framework this is indicative of athermal energy distribution. The steady-state temperatures for binary methanol/benzene mixtures in... 

Equilibrium molecular dynamics is typically applied to an isolated system containing a fixed number of molecules in a fixed volume V. Because the system is isolated, the total energy E is also constant. E is the sum of the molecular kinetic and potential energies. Hence, the variable V, N and E determine the thermodynamic state. 

In non-equilibrium molecular dynamics (NEMD), a driving force is introduced which maintains the system out of equilibrium at steady state, or else a perturbation is introduced and the system studied as it relaxes towards equilibrium. Our simulations are of the latter type. Generation of the initial cells... 

II. DESCRIPTION OF EQUILIBRIUM MOLECULAR DYNAMICS The Ordinary MD System... 

Salanne, M., Marrocchelli, D., Merlet, C. et al. (2011) Thermal conductivity of ionic systems from equilibrium molecular dynamics. J. Phys. Condens. Matter, 23, 102101. 

As in equilibrium molecular dynamics, the equations of motion have to be solved for a system with periodic boundaries. For shear, the boundaries are modified to become the Lees-Edwards sliding brick conditions (Lees Edwards 1972), in which periodic images of the simulation cell above and below the unit cell are moved in opposite directions at a velocity determined by the imposed shear rate (see Fig. 9.9). The properties of the system follow firom the appropriate time averages, <. . >, usually (but not necessarily) after the system has reached the steady state. Given, for example, a system at a number density, n = N/V, under an applied shear rate, the kinetic temperature is constrained with an appropriate thermostat Different properties can then be evaluated, for example, the internal energy. 

It is quite remarquable that a relatively small system exhibits a behaviour which can be understood in terms of macroscopic hydrodynamics. However, it should be stressed that, compared to "usual" equilibrium molecular dynamics, the system described here is an order of magnitude larger and has... 

MARESCHAL - In the Benard problem, the thermal boundaries are simulated along the ways developped in non-equilibrium molecular dynamics, using stochastic boundary conditions (see G. Ciccotti). The boundary layer does not extend over more than a mean free path in the system and can hardly be seen in our measurements. 

Non-equilibrium Molecular Dynamics Simulations of Coarse-Grained Polymer Systems... 

A few theoretical and computational studies have already addressed in some detail the problem of viscosity in ILs.[136] However, a complete microscopic theory of viscosity is currently not available. It is a challenging task to accurately compute the viscosity of a complex system by means of simulation methods. For a system with high viscosity, it is extremely difficult to reach the hydrodynamic limit (zero wave number) where the experimental data is observed. This is because, in order to reach this limit, a very large simulation box is required. Traditional simulation methods normally used for shear viscosity of fluids fall into two categories (a) the evaluation of the transverse-current autocorrelation function (TCAC) through equilibrium molecular dynamics (HMD) trajectories and (b) non-equilibrium molecular dynamics (NEMD) simulations that impose a periodic perturbation. [137] In recent work, Hess[138] compared most of the above methods by performing simulations of Lermard-Jones and water system. They concluded that the NEMD method using a periodic shear perturbation can be the best option. 

A non-equilibrium molecular dynamics procedure is used along with an established fixed charge force field. It is foimd that the simulations quantitatively capture the temperature dependence of the viscosity as well as the drop in viscosity that occurs with increasing water content. Using mixture viscosity models, the authors showed that the relative drop in viscosity with water content is actually less than the predicted values for an ideal system. This means that dissolved water is actually less effective at lowering the viscosity of these mixtures when compared to a solute obeying ideal mixing behaviour. 

Two simulation methods—Monte Carlo and molecular dynamics—allow calculation of the density profile and pressure difference of Eq. III-44 across the vapor-liquid interface [64, 65]. In the former method, the initial system consists of N molecules in assumed positions. An intermolecule potential function is chosen, such as the Lennard-Jones potential, and the positions are randomly varied until the energy of the system is at a minimum. The resulting configuration is taken to be the equilibrium one. In the molecular dynamics approach, the N molecules are given initial positions and velocities and the equations of motion are solved to follow the ensuing collisions until the set shows constant time-average thermodynamic properties. Both methods are computer intensive yet widely used. 

In equilibrium statistical mechanics, one is concerned with the thennodynamic and other macroscopic properties of matter. The aim is to derive these properties from the laws of molecular dynamics and thus create a link between microscopic molecular motion and thennodynamic behaviour. A typical macroscopic system is composed of a large number A of molecules occupying a volume V which is large compared to that occupied by a molecule ... 

The alternative simulation approaches are based on molecular dynamics calculations. This is conceptually simpler that the Monte Carlo method the equations of motion are solved for a system of A molecules, and periodic boundary conditions are again imposed. This method pennits both the equilibrium and transport properties of the system to be evaluated, essentially by numerically solvmg the equations of motion... 

Radiation probes such as neutrons, x-rays and visible light are used to see the structure of physical systems tlirough elastic scattering experunents. Inelastic scattering experiments measure both the structural and dynamical correlations that exist in a physical system. For a system which is in thennodynamic equilibrium, the molecular dynamics create spatio-temporal correlations which are the manifestation of themial fluctuations around the equilibrium state. For a condensed phase system, dynamical correlations are intimately linked to its structure. For systems in equilibrium, linear response tiieory is an appropriate framework to use to inquire on the spatio-temporal correlations resulting from thennodynamic fluctuations. Appropriate response and correlation functions emerge naturally in this framework, and the role of theory is to understand these correlation fiinctions from first principles. This is the subject of section A3.3.2. 

---