Molecular dynamics simulation 54,1-order

Since the development of grazing incidence x-ray diffraction, much of the convincing evidence for long-range positional order in layers has come from this technique. Structural relaxations from distorted hexagonal structure toward a relaxed array have been seen in heneicosanol [215]. Rice and co-workers combine grazing incidence x-ray diffraction with molecular dynamics simulations to understand several ordering transitions [178,215-219]. 

Prenkel, D. Pree energy computation and first order phase transitions. In Molecular Dynamic Simulation of Statistical Mechanical Systems, Enrico Fermi Summer School, Varenna 1985, G. Ciccotti and W. Hoover, eds. North Holland, Amsterdam (1986) 43-65. 

The principal idea behind the CSP approach is to use input from Classical Molecular Dynamics simulations, carried out for the process of interest as a first preliminary step, in order to simplify a quantum mechanical calculation, implemented in a subsequent, second step. This takes advantage of the fact that classical dynamics offers a reasonable description of many properties of molecular systems, in particular of average quantities. More specifically, the method uses classical MD simulations in order to determine effective... 

One drawback to a molecular dynamics simulation is that the trajectory length calculated in a reasonable time is several orders of magnitude shorter than any chemical process and most physical processes, which occur in nanoseconds or longer. This allows yon to study properties that change w ithin shorter time periods (such as energy finctnations and atomic positions), but not long-term processes like protein folding. 

Variation in Verlet order parameter during the equilibration phase of a molecular dynamics simulation of... 

Molecular dynamics calculations are more time-consuming than Monte Carlo calculations. This is because energy derivatives must be computed and used to solve the equations of motion. Molecular dynamics simulations are capable of yielding all the same properties as are obtained from Monte Carlo calculations. The advantage of molecular dynamics is that it is capable of modeling time-dependent properties, which can not be computed with Monte Carlo simulations. This is how diffusion coefficients must be computed. It is also possible to use shearing boundaries in order to obtain a viscosity. Molec-... 

Molecular Dynamics and Monte Carlo Simulations. At the heart of the method of molecular dynamics is a simulation model consisting of potential energy functions, or force fields. Molecular dynamics calculations represent a deterministic method, ie, one based on the assumption that atoms move according to laws of Newtonian mechanics. Molecular dynamics simulations can be performed for short time-periods, eg, 50—100 picoseconds, to examine localized very high frequency motions, such as bond length distortions, or, over much longer periods of time, eg, 500—2000 ps, in order to derive equiUbrium properties. It is worthwhile to summarize what properties researchers can expect to evaluate by performing molecular simulations ... 

In a different context, conformational analysis is essential for the analysis of molecular dynamics simulations. As discussed in Chapter 3, the direct output of a molecular dynamics simulation is a set of confonnations ( snapshots ) that were saved along the trajectory. These conformations are subsequently analyzed in order to extract information about the system. However, if, during a long simulation, the molecule moves from one... 

The concentration of salt in physiological systems is on the order of 150 mM, which corresponds to approximately 350 water molecules for each cation-anion pair. Eor this reason, investigations of salt effects in biological systems using detailed atomic models and molecular dynamic simulations become rapidly prohibitive, and mean-field treatments based on continuum electrostatics are advantageous. Such approximations, which were pioneered by Debye and Huckel [11], are valid at moderately low ionic concentration when core-core interactions between the mobile ions can be neglected. Briefly, the spatial density throughout the solvent is assumed to depend only on the local electrostatic poten-... 

Fig. 7 gives an example of such a comparison between a number of different polymer simulations and an experiment. The data contain a variety of Monte Carlo simulations employing different models, molecular dynamics simulations, as well as experimental results for polyethylene. Within the error bars this universal analysis of the diffusion constant is independent of the chemical species, be they simple computer models or real chemical materials. Thus, on this level, the simplified models are the most suitable models for investigating polymer materials. (For polymers with side branches or more complicated monomers, the situation is not that clear cut.) It also shows that the so-called entanglement length or entanglement molecular mass Mg is the universal scaling variable which allows one to compare different polymeric melts in order to interpret their viscoelastic behavior. 

The rapid rise in computer speed over recent years has led to atom-based simulations of liquid crystals becoming an important new area of research. Molecular mechanics and Monte Carlo studies of isolated liquid crystal molecules are now routine. However, care must be taken to model properly the influence of a nematic mean field if information about molecular structure in a mesophase is required. The current state-of-the-art consists of studies of (in the order of) 100 molecules in the bulk, in contact with a surface, or in a bilayer in contact with a solvent. Current simulation times can extend to around 10 ns and are sufficient to observe the growth of mesophases from an isotropic liquid. The results from a number of studies look very promising, and a wealth of structural and dynamic data now exists for bulk phases, monolayers and bilayers. Continued development of force fields for liquid crystals will be particularly important in the next few years, and particular emphasis must be placed on the development of all-atom force fields that are able to reproduce liquid phase densities for small molecules. Without these it will be difficult to obtain accurate phase transition temperatures. It will also be necessary to extend atomistic models to several thousand molecules to remove major system size effects which are present in all current work. This will be greatly facilitated by modern parallel simulation methods that allow molecular dynamics simulations to be carried out in parallel on multi-processor systems [115]. 

The simulation of a first-order phase transition, especially one where the two phases have a significant difference in molecular area, can be difficult in the context of a molecular dynamics simulation some of the works already described are examples of this problem. In a molecular dynamics simulation it can be hard to see coexistence of phases, especially when the molecules are fairly complicated so that a relatively small system size is necessary. One approach to this problem, described by Siepmann et al. [369] to model the LE-G transition, is to perform Monte Carlo simulations in the Gibbs ensemble. In this approach, the two phases are simulated in two separate but coupled boxes. One of the possible MC moves is to move a molecule from one box to the other in this manner two coexisting phases may be simulated without an interface. Siepmann et al. used the chain and interface potentials described in the Karaborni et al. works [362-365] for a 15-carbon carboxylic acid (i.e. pen-tadecanoic acid) on water. They found reasonable coexistence conditions from their simulations, implying, among other things, the existence of a stable LE state in the Karaborni model, though the LE phase is substantially denser than that seen experimentally. The re-... 

Tarek et al. [388] studied a system with some similarities to the work of Bocker et al. described earlier—a monolayer of n-tetradecyltrimethylammonium bromide. They also used explicit representations of the water molecules in a slab orientation, with the mono-layer on either side, in a molecular dynamics simulation. Their goal was to model more disordered, liquid states, so they chose two larger molecular areas, 0.45 and 0.67 nm molecule Density profiles normal to the interface were calculated and compared to neutron reflectivity data, with good agreement reported. The hydrocarbon chains were seen as highly disordered, and the diffusion was seen at both areas, with a factor of about 2.5 increase from the smaller molecular area to the larger area. They report no evidence of a tendency for the chains to aggregate into ordered islands, so perhaps this work can be seen as a realistic computer simulation depiction of a monolayer in an LE state. 

FIG. 26 Optimized structure of a water monolayer on mica obtained from molecular dynamic simulations by Odelius et al. The water molecules and the first layer of sihca tetrahedra of the mica substrate are shown in a side view in the top. K ions are the large dark balls. The bottom drawing shows a top view of the water. Oxygen atoms are dark, hydrogen atoms light. Notice the ordered icelike structure and the absence of free OH groups. All the H atoms in the water are involved in a hydrogen bond to another water molecule or to the mica substrate. (From Ref. 73.)... 

---