Molecular dynamics chain molecules

In the next section we describe the basic models that have been used in simulations so far and summarize the Monte Carlo and molecular dynamics techniques that are used. Some principal results from the scaling analysis of EP are given in Sec. 3, and in Sec. 4 we focus on simulational results concerning various aspects of static properties the MWD of EP, the conformational properties of the chain molecules, and their behavior in constrained geometries. The fifth section concentrates on the specific properties of relaxation towards equilibrium in GM and LP as well as on the first numerical simulations of transport properties in such systems. The final section then concludes with summary and outlook on open problems. 

The particular models used to demonstrate the theory obviously have many drawbacks as true representations of polymer crystals. These could include the lack of a fold energy, no distinction between new molecules and those already attached, neglect of chain ends, a somewhat arbitrary choice of pinning rules etc. However, they all serve their purpose in that they show that an energetic free energy barrier is not necessary to obtain the experimental curves. A truly representative growth picture can probably only be achieved via molecular dynamics. 

While thin polymer films may be very smooth and homogeneous, the chain conformation may be largely distorted due to the influence of the interfaces. Since the size of the polymer molecules is comparable to the film thickness those effects may play a significant role with ultra-thin polymer films. Several recent theoretical treatments are available [136-144,127,128] based on Monte Carlo [137-141,127, 128], molecular dynamics [142], variable density [143], cooperative motion [144], and bond fluctuation [136] model calculations. The distortion of the chain conformation near the interface, the segment orientation distribution, end distribution etc. are calculated as a function of film thickness and distance from the surface. In the limit of two-dimensional systems chains segregate and specific power laws are predicted [136, 137]. In 2D-blends of polymers a particular microdomain morphology may be expected [139]. Experiments on polymers in this area are presently, however, not available on a molecular level. Indications of order on an... 

Using the first-principles molecular-dynamics simulation, Munejiri, Shimojo and Hoshino studied the structure of liquid sulfur at 400 K, below the polymerization temperature [79]. They found that some of the Ss ring molecules homolytically open up on excitation of one electron from the HOMO to the LUMO. The chain-like diradicals S " thus generated partly recombine intramolecularly with formation of a branched Sy=S species rather than cyclo-Ss- Furthermore, the authors showed that photo-induced polymerization occurs in liquid sulfur when the Ss chains or Sy=S species are close to each other at their end. The mechanism of polymerization of sulfur remains a challenging problem for further theoretical work. 

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

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. 

Studies of the effect of permeant s size on the translational diffusion in membranes suggest that a free-volume model is appropriate for the description of diffusion processes in the bilayers [93]. The dynamic motion of the chains of the membrane lipids and proteins may result in the formation of transient pockets of free volume or cavities into which a permeant molecule can enter. Diffusion occurs when a permeant jumps from a donor to an acceptor cavity. Results from recent molecular dynamics simulations suggest that the free volume transport mechanism is more likely to be operative in the core of the bilayer [84]. In the more ordered region of the bilayer, a kink shift diffusion mechanism is more likely to occur [84,94]. Kinks may be pictured as dynamic structural defects representing small, mobile free volumes in the hydrocarbon phase of the membrane, i.e., conformational kink g tg ) isomers of the hydrocarbon chains resulting from thermal motion [52] (Fig. 8). Small molecules can enter the small free volumes of the kinks and migrate across the membrane together with the kinks. 

---