Molecular dynamics pressure

Fundamental frequency of the resonator Correlation function for surface roughness Root mean square height of a roughness Wave vector of shear waves in quartz, (Uy pq//rq Correlation length of surface roughness Thickness of the liquid film Thickness of interfacial layer Molecular dynamics Pressure in a liquid Quartz crystal microbalance Hydrodynamic roughness factor Electrochemical roughness factor Coordinates (normal and lateral)... 

The three formulas (29), (35), and (42) can all be used to determine the molecular dynamics pressure. The resulting values from a molecular dynamics calculation for 12 hard disks at a reduced volume r = 2 are shown in Table 1 along with the value reported by Hoover and Alder from (29) and the value... 

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 Fig. III-7 we show a molecular dynamics computation for the density profile and pressure difference P - p across the interface of an argonlike system [66] (see also Refs. 67, 68 and citations therein). Similar calculations have been made of 5 in Eq. III-20 [69, 70]. Monte Carlo calculations of the density profile of the vapor-liquid interface of magnesium how stratification penetrating about three atomic diameters into the liquid [71]. Experimental measurement of the transverse structure of the vapor-liquid interface of mercury and gallium showed structures that were indistinguishable from that of the bulk fluids [72, 73]. 

The Langmuir-Hinshelwood picture is essentially that of Fig. XVIII-14. If the process is unimolecular, the species meanders around on the surface until it receives the activation energy to go over to product(s), which then desorb. If the process is bimolecular, two species diffuse around until a reactive encounter occurs. The reaction will be diffusion controlled if it occurs on every encounter (see Ref. 211) the theory of surface diffusional encounters has been treated (see Ref. 212) the subject may also be approached by means of Monte Carlo/molecular dynamics techniques [213]. In the case of activated bimolecular reactions, however, there will in general be many encounters before the reactive one, and the rate law for the surface reaction is generally written by analogy to the mass action law for solutions. That is, for a bimolecular process, the rate is taken to be proportional to the product of the two surface concentrations. It is interesting, however, that essentially the same rate law is obtained if the adsorption is strictly localized and species react only if they happen to adsorb on adjacent sites (note Ref. 214). (The apparent rate law, that is, the rate law in terms of gas pressures, depends on the form of the adsorption isotherm, as discussed in the next section.)... 

AokI M I and Tsumuraya K 1997 Ab initio molecular-dynamics study of pressure-induced glass-to-crystal transitions In the sodium system Pbys. Rev. B 56 2962-8... 

Andersen H C 1980 Molecular dynamics simulations at constant pressure and/or temperature J. Chem. 

A typical molecular dynamics simulation comprises an equflibration and a production phase. The former is necessary, as the name imphes, to ensure that the system is in equilibrium before data acquisition starts. It is useful to check the time evolution of several simulation parameters such as temperature (which is directly connected to the kinetic energy), potential energy, total energy, density (when periodic boundary conditions with constant pressure are apphed), and their root-mean-square deviations. Having these and other variables constant at the end of the equilibration phase is the prerequisite for the statistically meaningful sampling of data in the following production phase. 

The forces are calculated as part of a molecular dynamics simulation, cind so little additional effort is required to calculate the virial and thus the pressure. The forces are not routinely calculated during a Monte Carlo simulation, and so some additional effort is required to determine the pressure by this route. When calculating the pressure it is also important to check that the components of the pressure in all three directions are equal. 

The thermodynamic properties that we have considered so far, such as the internal energy, the pressure and the heat capacity are collectively known as the mechanical properties and can be routinely obtained from a Monte Carlo or molecular dynamics simulation. Other thermodynamic properties are difficult to determine accurately without resorting to special techniques. These are the so-called entropic or thermal properties the free energy, the chemical potential and the entropy itself. The difference between the mechanical emd thermal properties is that the mechanical properties are related to the derivative of the partition function whereas the thermal properties are directly related to the partition function itself. To illustrate the difference between these two classes of properties, let us consider the internal energy, U, and the Fielmholtz free energy, A. These are related to the partition function by ... 

Calculating the statistical efficiency, a. A plot of tj,a A)i,la A) against tj, shows a steep rise before j off. Here the property A corresponds to the pressure calculated from the molecular dynamics simulation of... 

Just as one may wish to specify the temperature in a molecular dynamics simulation, so may be desired to maintain the system at a constant pressure. This enables the behavior of the system to be explored as a function of the pressure, enabling one to study phenomer such as the onset of pressure-induced phase transitions. Many experimental measuremen are made under conditions of constant temperature and pressure, and so simulations in tl isothermal-isobaric ensemble are most directly relevant to experimental data. Certai structural rearrangements may be achieved more easily in an isobaric simulation than i a simulation at constant volume. Constant pressure conditions may also be importai when the number of particles in the system changes (as in some of the test particle methoc for calculating free energies and chemical potentials see Section 8.9). 

The pressure often fluctuates much more than quantities such as the total energy in constant NVE molecular dynamics simulation. This is as expected because the pressure related to the virial, which is obtained as the product of the positions and the derivativ of the potential energy function. This product, rijdf rij)/drij, changes more quickly with than does the internal energy, hence the greater fluctuation in the pressure. 

Statistical mechanics computations are often tacked onto the end of ah initio vibrational frequency calculations for gas-phase properties at low pressure. For condensed-phase properties, often molecular dynamics or Monte Carlo calculations are necessary in order to obtain statistical data. The following are the principles that make this possible. 

In a normal molecular dynamics simulation with repeating boundary conditions (i.e., periodic boundary condition), the volume is held fixed, whereas at constant pressure the volume of the system must fluemate. In some simulation cases, such as simulations dealing with membranes, it is more advantageous to use the constant-pressure MD than the regular MD. Various schemes for prescribing the pressure of a molecular dynamics simulation have also been proposed and applied [23,24,28,29]. In all of these approaches it is inevitable that the system box must change its volume. 

An algorithm for performing a constant-pressure molecular dynamics simulation that resolves some unphysical observations in the extended system (Andersen s) method and Berendsen s methods was developed by Feller et al. [29]. This approach replaces the deterministic equations of motion with the piston degree of freedom added to the Langevin equations of motion. This eliminates the unphysical fluctuation of the volume associated with the piston mass. In addition, Klein and coworkers [30] present an advanced constant-pressure method to overcome an unphysical dependence of the choice of lattice in generated trajectories. 

Whenever the polymer crystal assumes a loosely packed hexagonal structure at high pressure, the ECC structure is found to be realized. Hikosaka [165] then proposed the sliding diffusion of a polymer chain as dominant transport process. Molecular dynamics simulations will be helpful for the understanding of this shding diffusion. Folding phenomena of chains are also studied intensively by Monte Carlo methods and generalizations [166,167]. 

Numerical simulation and molecular dynamics simulation on the removal action of CMP have been widely studied in recent years. In 1927, Preston [125] presented the mechanical model which relates the removal rate to the down pressure and relative velocity as follows ... 

FIG. 23 Surface pressure vs. area/molecule isotherms at 300 K from molecular dynamics simulations of Karaborni et al. (Refs. 362-365). All are for hydrocarbon chains with carboxylate-like head groups, (a) (filled squares) A 20-carbon chain, (b) (filled circles) A 16-carbon chain with a square simulation box the curve is shifted 5 A to the right, (c) (open squares) A 16-carbon chain with a nonsquare box with dimensions in the ratio xly = (3/4) to fit a hexagonal lattice the curve is shifted 5 A to the right. (Reproduced with permission from Ref. 365. Copyright 1993 American Chemical Society.)... 

In recent years, a class of methods has been developed for molecular dynamics simulations to be performed with an external pH parameter, like temperature or pressure [18, 43, 44, 70], These methods treat the solution as an infinite proton bath, and are thus referred to as constant pH molecular dynamics (PHMD). In PHMD, conformational dynamics of a protein is sampled simultaneously with the protonation states as a function of pH. As a result, protein dielectric response to the... 

An interesting extension of the original methodology was proposed by Lopes and Tildesley to allow the study of more than two phases at equilibrium [21], The extension is based on setting up a simulation with as many boxes as the maximum number of phases expected to be present. Kristof and Liszi [22, 23] have proposed an implementation of the Gibbs ensemble in which the total enthalpy, pressure and number of particles in the total system are kept constant. Molecular dynamics versions of the Gibbs ensemble algorithm are also available [24-26]. 

---