Lennard-Jones potential energy function

The reaction of trimethylene biradical was successfully treated by means of dynamics simulations by two groups with different PESs as described above.11 15 The success led one of the groups to extend the study to analyze the collisional and frictional effects in the trimethylene decomposition in an argon bath.16 A mixed QM/MM direct dynamics trajectory method was used with argon as buffer medium. Trimethylene intramolecular potential was treated by AM1-SRP fitted to CASSCF as before, and intermolecular forces were determined from Lennard-Jones 12-6 potential energy functions. 

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.

Lennard-Jones gas dimensionless thermal diffusion ratio, 25 307 Lennard-Jones molecules, 25 302 Lennard-Jones potential, 7 620 23 94 Lennard-Jones potential energy function, 25 302... 

In PLPscore, both the steric and hydrogen bonding terms are calculated from a piecewise linear potential function (see Fig. 1), instead of a smooth 6-12 Lennard-Jones potential energy function. The difference between the two terms is simply in the parameter values chosen for each term. 

The Lennard-Jones potential energy function systematically describes the attractive and repulsive forces at all distances in terms of s the well depth and a the location of the potential minimum. 

AB > AB = constants in the Lennard-Jones potential-energy function for the molecular pair AB Oab is in A. 

Since the Lennard-Jones potential-energy function is used, the equation is strictly valid only for nonpolar gases. The Lennard-Jones constants for the unlike molecular pair AB can be estimated from the constants for like pairs A A and BB ... 

Fig. 2.10. Comparison of Dymond-Alder (—) and Lennard-Jones (--) potential-energy functions for argon. The insert compares them in the domain where repulsive forces dominate.

The second simulation technique is molecular dynamics. In this technique, which was pioneered by Alder, initial positions of theparticles of a system of several hundred particles are assigned in some way. Displacements of the particles are determined by numerically simulating the classical equations of motion. Periodic boundary conditions are applied as in the Monte Carlo method. The first molecular dynamics calculations were done on systems of hard spheres, but the method has been applied to monatomic systems having intermolecular forces represented by the square-well and Lennard-Jones potential energy functions, as well as on model systems representing molecular substances. Commercial software is now available to carry out molecular dynamics simulations on desktop computers. ... 

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 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 8-1. A plot of Eq. (8-16), the Lennard-Jones 6-12 potential energy function.

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

We now describe a relatively simple MD model of a low-index crystal surface, which was conceived for the purpose of studying the rate of mass transport (8). The effect of temperature on surface transport involves several competing processes. A rough surface structure complicates the trajectories somewhat, and the diffusion of clusters of atoms must be considered. In order to simplify the model as much as possible, but retain the essential dynamics of the mobile atoms, we will consider a model in which the atoms move on a "substrate" represented by an analytic potential energy function that is adjusted to match that of a surface of a (100) face-centered cubic crystal composed of atoms interacting with a Lennard-Jones... 

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. 

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

A considerable amount of work has been done on the development of water-ion potential energy functions." " Most of these functions are of the standard Lennard-Jones plus Coulomb form, with parameters selected to give the experimental free energy or enthalpy of solvation. ... 

---