Lennard-Jones 6-12 function/term

June et al. (12) used TST as an alternative method to investigate Xe diffusion in silicalite. Interactions between the zeolite oxygen atoms and the Xe atoms were modeled with a 6-12 Lennard-Jones function, with potential parameters similar to those used in previous MD simulations (11). Simulations were performed with both a rigid and a flexible zeolite lattice, and those that included flexibility of the zeolite framework employed a harmonic term to describe the motion of the zeolite atoms, with a force constant and bond length data taken from previous simulations (26). 

The interaction parameters for the water molecules were taken from nonempirical configuration interaction calculations for water dimers (41) that have been shown to give good agreement between experimental radial distribution functions and simulations at low sorbate densities. The potential terms for the water-ferrierite interaction consisted of repulsion, dispersion, and electrostatic terms. The first two of these terms are the components of the 6-12 Lennard-Jones function, and the electrostatic term accounts for long-range contributions and is evaluated by an Ewald summation. The... 

Auerbach et al. (101) used a variant of the TST model of diffusion to characterize the motion of benzene in NaY zeolite. The computational efficiency of this method, as already discussed for the diffusion of Xe in NaY zeolite (72), means that long-time-scale motions such as intercage jumps can be investigated. Auerbach et al. used a zeolite-hydrocarbon potential energy surface that they recently developed themselves. A Si/Al ratio of 3.0 was assumed and the potential parameters were fitted to reproduce crystallographic and thermodynamic data for the benzene-NaY zeolite system. The functional form of the potential was similar to all others, including a Lennard-Jones function to describe the short-range interactions and a Coulombic repulsion term calculated by Ewald summation. 

The intermolecular potential term is represented by a simple Lennard-Jones function that is attenuated at short interatomic distances by a cubic spline so that at small (covalent) intemuclear distances, the description of the interaction is that of the intramolecular term only. The original form of... 

The interaction between the solvent and solute in these solutions can, in many cases, be described simply in terms of the Van der Waals forces. The corresponding total potential, W, can then be given by the Lennard-Jones function ... 

A later work dealt with two questions (1) what effect does a change from a three-parameter Buckingham to a two-parameter Lennard-Jones function have, nothing else being touched (2) when Coulomb terms are included, what is the effect of using the same charge for all carbon atoms, also the anomeric ... 

As a third test case, we wanted an Ionic compound, where the Coulomb Interaction Is dominating. We chose KCl because we wanted to use Ar parameters In the van der Waals terms (the Lennard-Jones functions). NaCl was also used In the convergence tests. 

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. 

The first two terms on the right-hand side of Eq. (83) are usually assumed to be harmonic, as given for example by Eq. (6-74). The third term is often developed in a Fourier series, as given by Eq. (82). The potential function appropriate to the interaction between nonbonded atoms is taken to be of the Lennard-Jones type (Section 6.7.3). In all of these cases the necessary force constants are estimated by comparing the results obtained from a large number of similar molecules. If electrostatic interactions are to be considered, effective atomic charges must be suggested and Coulomb s law applied directly [see Eq. (6-81)]. 

Here, avaw is a positive constant, and and ctJ are the usual Lennard-Jones parameters found in macromolecular force fields. The role played by the term avdw (1 — A)2 in the denominator is to eliminate the singularity of the van der Waals interaction. Introduction of this soft-core potential results in bounded derivatives of the potential energy function when A tends towards 0. 

A better approach than non-linear scaling is to attempt to reduce and/or eliminate the singularity in the function that occurs on the step when a noninteracting group starts to interact. A clever approach has been described that reduces the problem by modifying the Lennard-Jones van der Waals term in the potential function.30,31 For a pair of atoms where one group vanishes at the X =1 endpoint, the modified Lennard-Jones 6-12 function takes the form ... 

In a statistical Monte Carlo simulation the pair potentials are introduced by means of analytical functions. In the election of that analytical form for the pair potential, it must be considered that when a Monte Carlo calculation is performed, the more time consuming step is the evaluation of the energy for the different configurations. Given that this calculation must be done millions of times, the chosen analytic functions must be of enough accuracy and flexibility but also they must be fastly computed. In this way it is wise to avoid exponential terms and to minimize the number of interatomic distances to be calculated at each configuration which depends on the quantity of interaction centers chosen for each molecule. A very commonly used function consists of a sum of rn terms, r being the distance between the different interaction centers, usually, situated at the nuclei. In particular, non-bonded interactions are usually represented by an atom-atom centered monopole expression (Coulomb term) plus a Lennard-Jones 6-12 term, as indicated in equation (51). 

The OPLS model is an example of pair potential where non-bonded interactions are represented through Coulomb and Lennard-Jones terms interacting between sites centred on nuclei (equation (51). Within this model, each atomic nucleus has an interaction site, except CH groups that are treated as united atoms centered on the carbon. It is important to note that no special functions were found to be needed to describe hydrogen bonding and there are no additional interaction sites for lone pairs. Another important point is that standard combining rules are used for the Lennard-Jones interactions such that An = (Ai As )1/2 and Cu = (C Cy)1/2. The A and C parameters may also be expressed in terms of Lennard-Jones o s and e s as A = 4ei Oi and C ... 

In the study of reactivity, Jorgensen and col. have normally used both, the OPLS model and potential functions derived from ab initio calculations. As we have already indicated, when intermolecular pair potentials are applied to the study of a chemical process, the evolution of charges, as well as the Lennard-Jones terms, along the reaction coordinate, has to be considered. For the SN2 reaction in water between chloride anion... 

Here the atoms in the system are numbered by i, j, k, l = 1,..., N. The distance between two atoms i, j is ry, q is the (partial) charge on an atom, 6 is the angle defined by the coordinates (i, j, k) of three consecutive atoms, and 4> is the dihedral angle defined by the positions of four consecutive atoms, e0 is the dielectric permittivity of vacuum, n is the dihedral multiplicity. The potential function, as given in equation (6), has many parameters that depend on the atoms involved. The first term accounts for Coulombic interactions. The second term is the Lennard-Jones interaction energy. It is composed of a strongly repulsive term and a van der Waals-like attractive term. The form of the repulsive term is chosen ad hoc and has the function of defining the size of the atom. The Ay coefficients are a function of the van der Waals radii of the... 

The first three terms, stretch, bend and torsion, are common to most force fields although their explicit form may vary. The nonbonded terms may be further divided into contributions from Van der Waals (VdW), electrostatic and hydrogen-bond interactions. Most force fields include potential functions for the first two interaction types (Lennard-Jones type or Buckingham type functions for VdW interactions and charge-charge or dipole-dipole terms for the electrostatic interactions). Explicit hydrogen-bond functions are less common and such interactions are often modeled by the VdW expression with special parameters for the atoms which participate in the hydrogen bond (see below). 

---