Molecular dynamics mixtures

It was noted in connection with Eq. III-56 that molecular dynamics calculations can be made for a liquid mixture of rare gas-like atoms to obtain surface tension versus composition. The same calculation also gives the variation of density for each species across the interface [88], as illustrated in Fig. Ill-13b. The density profiles allow a calculation, of course, of the surface excess quantities. 

Surface tension is usually predicted using group additivity methods for neat liquids. It is much more difficult to predict the surface tension of a mixture, especially when surfactants are involved. Very large molecular dynamics or Monte Carlo simulations can also be used. Often, it is easier to measure surface tension in the laboratory than to compute it. 

FIQ. 3 Diffusion coefficient of benzene molecules in benzene-polystyrene mixtures normalized by the diffusion coefficient of neat benzene molecular dynamics results, NMR measurements and prediction by the Mackie-Meares model [26]. 

The only feasible procedure at the moment is molecular dynamics computer simulation, which can be used since most systems are currently essentially controlled by classical dynamics even though the intermolecular potentials are often quantum mechanical in origin. There are indeed many intermolecular potentials available which are remarkably reliable for most liquids, and even for liquid mixtures, of scientific and technical importance. However potentials for the design of membranes and of the interaction of fluid molecules with membranes on the atomic scale are less well developed. 

J. D. Weeks, D. Chandler, H. C. Andersen. Role of repulsive forces in determining the equilibrium structure of simple liquids. J Chem Phys 54 5237, 1971. R. L. Rowley, M. W. Schuck, J. Perry. A direct method for determination of chemical potential with molecular dynamics simulations. 2. Mixtures. Mol Phys 55 125, 1995. 

Dielectric relaxation measurements of polyethylene grafted with acrylic acid(AA), 2-hydroxyethyl methacrylate (HEMA) and their binary mixture were carried out in a trial to explore the molecular dynamics of the grafted samples [125]. Such measurements provide information about their molecular packing and interaction. It was possible to predict that the binary mixture used yields a random copolymer PE—g—P(AA/HEMA), which is greatly enriched with HEMA. This method of characterization is very interesting and is going to be developed in different polymer/monomer systems. 

The lattice gas has been used as a model for a variety of physical and chemical systems. Its application to simple mixtures is routinely treated in textbooks on statistical mechanics, so it is natural to use it as a starting point for the modeling of liquid-liquid interfaces. In the simplest case the system contains two kinds of solvent particles that occupy positions on a lattice, and with an appropriate choice of the interaction parameters it separates into two phases. This simple version is mainly of didactical value [1], since molecular dynamics allows the study of much more realistic models of the interface between two pure liquids [2,3]. However, even with the fastest computers available today, molecular dynamics is limited to comparatively small ensembles, too small to contain more than a few ions, so that the space-charge regions cannot be included. In contrast, Monte Carlo simulations for the lattice gas can be performed with 10 to 10 particles, so that modeling of the space charge poses no problem. In addition, analytical methods such as the quasichemical approximation allow the treatment of infinite ensembles. 

Yu HB, Geerke DP, Liu HY, van Gunsteren WE (2006) Molecular dynamics simulations of liquid methanol and methanol-water mixtures with polarizable models. J Comput Chem 27(13) 1494-1504... 

Noskov S Y, Lamoureux G, Roux B (2005) Molecular dynamics study of hydration in ethanol-water mixtures using a polarizable force field. J Phys Chem B 109(14) 6705-6713... 

Reaction Ensemble Molecular Dynamics Direct Simulation of the Dynamic Equilibrium Properties of Chemically Reacting Mixtures. 

D. R Geerke and W. F. Van Gunsteren, The performance of non polarizable and polarizable force field parameter sets for ethylene glycol in molecular dynamics simulations of the pure liquid and its aqueous mixtures. Mol. Phys. 105, 1861 1881 (2007). 

During the last few years the progress of computational techniques has made it possible to simulate the dynamic behavior of whole ensembles consisting of several hundred molecules. In this way the limitations of the statistical approach can be at least partly overcome. Two kinds of methods — molecular dynamics and Monte Carlo calculations — were applied to liquids and liquid mixtures and brought new insight into their structure and properties. Even some important characteristics of systems as complicated as associated liquids like water could be... 

G. B. Dutt, S. Doraiswamy and N. Periasamy, Molecular dynamics of polar dye probes in tertiary-butyl alcohol-water mixtures,/. Chem. Phys. 94, 5360-5368 (1991). 

Zhang, Z., Bhide, S.Y., Berkowitz, M.L. Molecular dynamics simulations of bilayers containing mixtures of sphingomyelin with cholesterol and phosphatidylcholine with cholesterol. J. Phys. Chem. B 2007, 111, 12888-97. 

Bhargava, B.L., and Balasubramanian, S., Insights into the structure and dynamics of a room-temperature ionic liquid Ab initio molecular dynamics simulation studies of l-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PFj]) and the [bmim][PFj]-C02 mixture, /. Phys. Chem. B, 111, 4477-4487, 2007. 

Sprik et al, 1993 Signorini et al, 1990), a typical example being the orientational disorder associated with NH in NH Br. Detailed simulations have been reported on (NaCN),, t(KCN),t and other mixed alkali halides and alkali cyanides. Other systems studied include potassium and calcium nitrate crystals and their mixtures. The transition from the crystalline to the superionic conductor phase in solid electrolytes has also been successfully investigated. Molecular dynamics studies of Agl were carried out by Parrinello, Rahman Vashishta (1983). LijSO has been investigated by molecular dynamics by Impey et al. (1985). Here, the Li ions become mobile at high temperatures. The ions exhibit orientational disorder and the orientational... 

---