Lennard-jones liquid

Sese, L. M., A quantum Monte Carlo study of liquid Lennard-Jones methane, path-integral and effective potentials, Mol. Phys. 1992, 76, 1335-1346... 

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. 

The gradient model for interfacial tension described in Eqs. III-42 and III-43 is limited to interaction potentials that decay more rapidly than r. Thus it can be applied to the Lennard-Jones potential but not to a longer range interaction such as dipole-dipole interaction. Where does this limitation come from, and what does it imply for interfacial tensions of various liquids ... 

Molecular dynamics and density functional theory studies (see Section IX-2) of the Lennard-Jones 6-12 system determine the interfacial tension for the solid-liquid and solid-vapor interfaces [47-49]. The dimensionless interfacial tension ya /kT, where a is the Lennard-Jones molecular size, increases from about 0.83 for the solid-liquid interface to 2.38 for the solid-vapor at the triple point [49], reflecting the large energy associated with a solid-vapor interface. 

Wilson M R 1997 Molecular dynamics simulations of flexible liquid crystal molecules using a Gay-Berne/Lennard-Jones model J. Chem. Phys. 107 8654-63... 

McDonald I R and Singer K 1967 Calculation of thermodynamic properties of liquid argon from Lennard-Jones parameters by a Monte Carlo method Discuss. Faraday Soc. 43 40-9... 

Lotfi A, Vrabeo J and Fisoher J 1992 Vapour liquid equilibria of the Lennard-Jones fluid from the NpT plus test partiole method Mol. Phys. 76 1319-33... 

Holcomb C D, Clancy P and Zollweg J A 1993 A critical study of the simulation of the liquid-vapour interface of a Lennard-Jones fluid Mol. Phys. 78 437-59... 

Rare-gas clusters can be produced easily using supersonic expansion. They are attractive to study theoretically because the interaction potentials are relatively simple and dominated by the van der Waals interactions. The Lennard-Jones pair potential describes the stmctures of the rare-gas clusters well and predicts magic clusters with icosahedral stmctures [139, 140]. The first five icosahedral clusters occur at 13, 55, 147, 309 and 561 atoms and are observed in experiments of Ar, Kr and Xe clusters [1411. Small helium clusters are difficult to produce because of the extremely weak interactions between helium atoms. Due to the large zero-point energy, bulk helium is a quantum fluid and does not solidify under standard pressure. Large helium clusters, which are liquid-like, have been produced and studied by Toennies and coworkers [142]. Recent experiments have provided evidence of... 

The range of systems that have been studied by force field methods is extremely varied. Some force fields liave been developed to study just one atomic or molecular sp>ecies under a wider range of conditions. For example, the chlorine model of Rodger, Stone and TUdesley [Rodger et al 1988] can be used to study the solid, liquid and gaseous phases. This is an anisotropic site model, in which the interaction between a pair of sites on two molecules dep>ends not only upon the separation between the sites (as in an isotropic model such as the Lennard-Jones model) but also upon the orientation of the site-site vector with resp>ect to the bond vectors of the two molecules. The model includes an electrostatic component which contciins dipwle-dipole, dipole-quadrupole and quadrupole-quadrupole terms, and the van der Waals contribution is modelled using a Buckingham-like function. 

The Morse function which is given above was obtained from a study of bonding in gaseous systems, and dris part of Swalin s derivation should probably be replaced with a Lennard-Jones potential as a better approximation. The general idea of a variable diffusion step in liquids which is more nearly akin to diffusion in gases than the earlier treatment, which was based on the notion of vacant sites as in solids, remains as a valuable suggestion. 

Recently, many experiments have been performed on the structure and dynamics of liquids in porous glasses [175-190]. These studies are difficult to interpret because of the inhomogeneity of the sample. Simulations of water in a cylindrical cavity inside a block of hydrophilic Vycor glass have recently been performed [24,191,192] to facilitate the analysis of experimental results. Water molecules interact with Vycor atoms, using an empirical potential model which consists of (12-6) Lennard-Jones and Coulomb interactions. All atoms in the Vycor block are immobile. For details see Ref. 191. We have simulated samples at room temperature, which are filled with water to between 19 and 96 percent of the maximum possible amount. Because of the hydrophilicity of the glass, water molecules cover the surface already in nearly empty pores no molecules are found in the pore center in this case, although the density distribution is rather wide. When the amount of water increases, the center of the pore fills. Only in the case of 96 percent filling, a continuous aqueous phase without a cavity in the center of the pore is observed. 

P. Bryk, A. Patrykiejew, O. Pizio, S. Sokolowski. The chemical potential of Lennard-Jones associating fluids from osmotic Monte Carlo simulations. Mol Phys 92 949, 1997 A method for the determination of chemical potential for associating liquids. Mol Phys 90 665, 1997. 

A complete set of intermolecular potential functions has been developed for use in computer simulations of proteins in their native environment. Parameters have been reported for 25 peptide residues as well as the common neutral and charged terminal groups. The potential functions have the simple Coulomb plus Lennard-Jones form and are compatible with the widely used models for water, TIP4P, TIP3P and SPC. The parameters were obtained and tested primarily in conjunction with Monte Carlo statistical mechanics simulations of 36 pure organic liquids and numerous aqueous solutions of organic ions representative of subunits in the side chains and backbones of proteins... 

Thermodynamic information can also be obtained from simulations. Currently we are measuring the differences in chemical potential of various small molecules in dimethylimidazolium chloride. This involves gradually transforming one molecule into another and is a computationally intensive process. One preliminary result is that the difference in chemical potential of propane and dimethyl ether is about 17.5 kj/mol. These molecules are similar in size, but differ in their polarity. Not surprisingly, the polar ether is stabilized relative to the non-polar propane in the presence of the ionic liquid. One can also investigate the local arrangement of the ions around the solute and the contribution of different parts of the interaction to the energy. Thus, while both molecules have a favorable Lennard-Jones interaction with the cation, the main electrostatic interaction is that between the chloride ion and the ether molecule. 

The theory introduced by Lennard-Jones and Devonshire13 17 for the study of liquids provides a powerful method for the quantitative evaluation of the partition function of a solute molecule within its cavity.51 Because the application of this method to the present problem has been described in detail,62 we shall restrict ourselves to its most essential features. 

In a review of the subject, Ubbelohde [3] points out that there is only a relatively small amount of data available concerning the properties of solids and also of the (product) liquids in the immediate vicinity of the melting point. In an early theory of melting, Lindemann [4] considered that when the amplitude of the vibrational displacements of the atoms of a particular solid increased with temperature to the point of attainment of a particular fraction (possibly 10%) of the lattice spacing, their mutual influences resulted in a loss of stability. The Lennard-Jones—Devonshire [5] theory considers the energy requirement for interchange of lattice constituents between occupation of site and interstitial positions. Subsequent developments of both these models, and, indeed, the numerous contributions in the field, are discussed in Ubbelohde s book [3]. 

---