Lennard solid

Molecular dynamics calculations have been made on atomic crystals using a Lennard-Jones potential. These have to be done near the melting point in order for the iterations not to be too lengthy and have yielded density functioi). as one passes through the solid-vapor interface (see Ref. 45). The calculations showed considerable mobility in the surface region, amounting to the presence of a... 

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. 

Agrawal R and Kofke D A 1995 Thermodynamio and struotural properties of model systems at solid-fluid ooexistenoe. II. Melting and sublimation of the Lennard-Jones system Mol. Phys. 85 43-59... 

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. 

Extensive computer simulations have been caiTied out on the near-surface and surface behaviour of materials having a simple cubic lattice structure. The interaction potential between pairs of atoms which has frequently been used for inert gas solids, such as solid argon, takes die Lennard-Jones form where d is the inter-nuclear distance, is the potential interaction energy at the minimum conesponding to the point of... 

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. 

Here Tq are coordinates in a reference volume Vq and r = potential energy of Ar crystals has been computed [288] as well as lattice constants, thermal expansion coefficients, and isotope effects in other Lennard-Jones solids. In Fig. 4 we show the kinetic and potential energy of an Ar crystal in the canonical ensemble versus temperature for different values of P we note that in the classical hmit (P = 1) the low temperature specific heat does not decrease to zero however, with increasing P values the quantum limit is approached. In Fig. 5 the isotope effect on the lattice constant (at / = 0) in a Lennard-Jones system with parameters suitable for Ne atoms is presented, and a comparison with experimental data is made. Please note that in a classical system no isotope effect can be observed, x "" and the deviations between simulations and experiments are mainly caused by non-optimized potential parameters. 

Lead, excess entropy of solution of noble metals in, 133 Lead-thalium, solid solution, 126 Lead-tin, system, energy of solution, 143 solution, enthalpy of formation, 143 Lead-zinc, alloy (Pb8Zn2), calculation of thermodynamic quantities, 136 Legendre expansion in total ground state wave function of helium, 294 Lennard-Jones 6-12 potential, in analy-... 

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

The system used in the simulations usually consists of solid walls and lubricant molecules, but the specific arrangement of the system depends on the problem under investigation. In early studies, hard spherical molecules, interacting with each other through the Lennard-Jones (L-J) potential, were adopted to model the lubricant [27], but recently we tend to take more realistic models for describing the lubricant molecules. The alkane molecules with flexible linear chains [28,29] and bead-spring chains [7,30] are the examples for the most commonly used molecular architectures. The inter- and intra-molecular potentials, as well as the interactions between the lubricant molecule and solid wall, have to be properly defined in order to get reliable results. Readers who intend to learn more about the specific techniques of the simulations are referred to Refs. [27-29]. 

Fig. 9 Solid line A typical Lennard-Jones (6-12) potential dotted line the effective potential for given values of the Initial parameters g and b.

Fig. 2.2. Average electrostatic potential mc at the position of the methane-like Lennard-Jones particle Me as a function of its charge q. mc contains corrections for the finite system size. Results are shown from Monte Carlo simulations using Ewald summation with N = 256 (plus) and N = 128 (cross) as well as GRF calculations with N = 256 water molecules (square). Statistical errors are smaller than the size of the symbols. Also included are linear tits to the data with q < 0 and q > 0 (solid lines). The fit to the tanh-weighted model of two Gaussian distributions is shown with a dashed line. Reproduced with permission of the American Chemical Society...

Fig. 10.12. Vapor-liquid phase behavior for the Lennard-Jones fluid. Solid triangles and hollow squares indicate the results of the particle addition/deletion and volume scaling variants of the flat-histogram simulation using the Wang-Landau algorithm. Crosses are from a histogram reweighting study based on grand-canonical measurements at seven state points. The solid line is from Lotfi, et al. [76], Reprinted figure with permission from [75]. 2002 by the American Physical Society...

Transition matrix estimators have received less attention than the multicanonical and Wang-Landau methods, but have been applied to a small collection of informative examples. Smith and Bruce [111, 112] applied the transition probability approach to the determination of solid-solid phase coexistence in a square-well model of colloids. Erring ton and coworkers [113, 114] have also used the method to determine liquid-vapor and solid-liquid [115] equilibria in the Lennard-Jones system. Transition matrices have also been used to generate high-quality data for the evaluation of surface tension [114, 116] and for the estimation of order parameter weights in phase-switch simulations [117]. 

Errington, J. R., Solid-liquid phase coexistence of the Lennard-Jones system through phase-switch Monte Carlo simulation, J. Chem. Phys. 2004,120, 3130-3141... 

Figure 3 Comparison of Lennard-Jones (solid line) and hard-sphere (circles) model ordering maps for the equilibrium fluid and the fee crystal. (Adapted from Ref. 32.)... 

The first MC (16) and MD (17) studies were used to simulate the properties of single particle fluids. Although the basic MC (11,12) and MD (12,13) methods have changed little since the earliest simulations, the systems simulated have continually increased in complexity. The ability to simulate complex interfacial systems has resulted partly from improvements in simulation algorithms (15,18) or in the interaction potentials used to model solid surfaces (19). The major reason, however, for this ability has resulted from the increasing sophistication of the interaction potentials used to model liquid-liquid interactions. These advances have involved the use of the following potentials Lennard-Jones 12-6 (20), Rowlinson (21), BNS... 

One of the first studies of multiple ions at the water/solid interface was by Spohr and Heinzinger, who carried out a simulation of a system of 8 Li" and 81" ions dissolved in 200 water molecules between uncharged flat Lennard-Jones walls.However, the issues discussed in their paper involved water structure and dynamics and the single-ion properties mentioned earlier. No attempt was made to consider the ions distributions and ion-ion correlations. This work has recently been repeated using more realistic water-metal potentials. ... 

Molecular dynamics has been used extensively to explore the solid-liquid interface. In one such study, a modified Lennard-Jones potential has been used to model this interaction in the spreading of a droplet [4], of the form... 

Theoretical calculations of the surface free energy of solids date back to 1928 with the work of Lennard-Jones and Dent 10). Displacements of the positive and negative ions when a given interior layer becomes a surface layer were allowed for by Verwey 11). Calculations by Shuttleworth 12) showed that van der Waals terms make a significant contribution to surface energy. Benson and his co-workers have made an extensive study of alkali-halides 13-16) and of magnesium oxide 17). 

---