Lennard function

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. 

One fascinating feature of the physical chemistry of surfaces is the direct influence of intermolecular forces on interfacial phenomena. The calculation of surface tension in section III-2B, for example, is based on the Lennard-Jones potential function illustrated in Fig. III-6. The wide use of this model potential is based in physical analysis of intermolecular forces that we summarize in this chapter. In this chapter, we briefly discuss the fundamental electromagnetic forces. The electrostatic forces between charged species are covered in Chapter V. 

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. 

The biasing function is applied to spread the range of configurations sampled such that the trajectory contains configurations appropriate to both the initial and final states. For the creation or deletion of atoms a softcore interaction function may be used. The standard Lennard-Jones (LJ) function used to model van der Waals interactions between atoms is strongly repulsive at short distances and contains a singularity at r = 0. This precludes two atoms from occupying the same position. A so-called softcore potential in contrast approaches a finite value at short distances. This removes the sin-... 

Fig. 3. Curves calculated using (8) for a series of increasing a values. The curves were calculated using tr = 0.6 nm and e = 0.4 kj/mol. Note that for a = 0.0 the normal 6-12 Lennard Jones potential energy function is recovered.

If computing time does not play the major role that it did in the early 1980s, the [12-6] Lennard-Jones potential is substituted by a variety of alternatives meant to represent the real situation much better. MM3 and MM4 use a so-called Buckingham potential (Eq. (28)), where the repulsive part is substituted by an exponential function ... 

A 6-12 function (also known as a Lennard-Jones function) frequently simulates van der Waats in tcraction s in force fields (ec iia-tion t 1). 

The Lennard-Jones 12-6 potential contains just two adjustable parameters the collision diameter a (the separation for which the energy is zero) and the well depth s. These parameters are graphically illustrated in Figure 4.34. The Lennard-Jones equation may also be expressed in terms of the separation at which the energy passes through a minimum, (also written f ). At this separation, the first derivative of the energy with respect to the internuclear distance is zero (i.e. dvjdr = 0), from which it can easily be shown that v = 2 / cr. We can thus also write the Lennard-Jones 12-6 potential function as follows ... 

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. 

Fig. 6.20 A switching function that applies over a narrow range near the cutoff and its effect on the Lennard-Janes potential.

The classical kinetic theoty of gases treats a system of non-interacting particles, but in real gases there is a short-range interaction which has an effect on the physical properties of gases. The most simple description of this interaction uses the Lennard-Jones potential which postulates a central force between molecules, giving an energy of interaction as a function of the inter-nuclear distance, r. 

Figure 3.7 The Lennard-Jones potential of the interaction of gaseous atoms as a function of the internuclear distance...

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. 

Fig. 11(a) displays plots of the in-plane pair correlation function for s = 2. and 3.0 well outside the regime where K exhibits its first maximum (see Fig. 12). The plots indicate that the transverse structures of one- and two-layer fluids (see Fig. 10) are essentially identical and typical of dense Lennard-Jones fluids. However, the transverse structure of a two-layer fluid is significantly affected as the peak of K is approached, as can be seen in Fig. 11(b) where g (zi,pi2) is plotted for s = 2.55 and 2.75, which points...

The integrals are over the full two-dimensional volume F. For the classical contribution to the free energy /3/d([p]) the Ramakrishnan-Yussouff functional has been used in the form recently introduced by Ebner et al. [314] which is known to reproduce accurately the phase diagram of the Lennard-Jones system in three dimensions. In the classical part of the free energy functional, as an input the Ornstein-Zernike direct correlation function for the hard disc fluid is required. For the DFT calculations reported, the accurate and convenient analytic form due to Rosenfeld [315] has been used for this quantity. 

The calculations have been carried out for a one-component, Lennard-Jones associating fluid with one associating site. The nonassociative van der Waals potential is thus given by Eq. (87) with = 2.5a, whereas the associative forces are described by means of Eq. (60), with d = 0.5contact with an attracting wall. The fluid-wall potential is given by the Lennard-Jones (9-3) function... 

The density functional approach has also been used to study capillary condensation in slit-like pores [148,149]. As in the previous section, a simple model of the Lennard-Jones associating fluid with a single associative site is considered. All the parameters of the interparticle potentials are chosen the same as in the previous section. Our attention has been focused on the influence of association on capillary condensation and the evaluation of the phase diagram [42]. 

Note that, for = 0, the potential given above does not reduce to the Lennard-Jones (12-6) function, because the soft Lennard-Jones repulsive branch is replaced by a hard-sphere potential, located at r = cr. The results for the nonassociating Lennard-Jones fluid can be found in Ref. 159. 

We have studied, by MD, pure water [22] and electrolyte solutions [23] in cylindrical model pores with pore diameters ranging from 0.8 to more than 4nm. In the nonpolar model pores the surface is a smooth cylinder, which interacts only weakly with water molecules and ions by a Lennard-Jones potential the polar pore surface contains additional point charges, which model the polar groups in functionalized polymer membranes. 

Figure 8-1. A plot of Eq. (8-16), the Lennard-Jones 6-12 potential energy function.

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

For small systems, where accurate interaction energy profiles are available, it has been shown that the Morse function actually gives a slightly better description than an Exp.-6, which again performs significantly better than a Lennard-Jones 12-6 potential. This is illustrated for the H2-He interaction in Figure 2.9. 

The main difference between the three functions is in the repulsive part at short distances the Lennard-Jones potential is much too hard, and the Exp.-6 also tends to overestimate the repulsion. It furthermore has the problem of inverting at short distances. For chemical purposes these problems are irrelevant, energies in excess of lOOkcal/mol are sufficient to break most bonds, and will never be sampled in actual calculations. The behaviour in the attractive part of the potential, which is essential for intermolecular interactions, is very similar for the three functions, as shown in... 

---