Potential: “effective 156 Lennard-Jones

Fig. 10. Effective Lennard-Jones potential of neon in the solid state. The solid curve represents the potential before averaging over the vibrational motion. The dashed curve results from the averaging.

More recently, Khoury and Robinson calculated virial coefficients for ethane + hydrogen sulphide mixtures from an effective Lennard-Jones 6—12 potential that includes a contribution from the dipole-induced-dipole interaction. The parameters were evaluated by fitting their measurements on CgHg and HgS to a 6—12 potential and to a Stockmayer potential respectively. Bradley and King have used several potential functions to analyze interactions of phen-anthrene with a number of small polyatomic molecules. 

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. 

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. 

The simplest atomistic model for the formation of a crystal in continuous space requires the definition of some effective attractive potential between any two atoms, which is defined independently of the other atoms in the cluster or crystal. The most frequently studied potential is the Lennard-Jones potential... 

Kihara20 used a core model in which the Lennard-Jones potential is assumed to hold for the shortest distance between the molecular cores instead of molecular centers. By use of linear, tetrahedral, and other shapes of cores, various molecules can be approximated. Thomaes,41 Rowlinson,35 Hamann, McManamey, and Pearse,14 Atoji and Lipscomb,1 Pitzer,30 and Balescu,4 have used other models of attracting centers and other mathemtical methods, but obtain similar conclusions. The primary effect is to steepen the potential curve so that in terms of inverse powers of the inter-... 

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

Multiparticle collision dynamics provides an ideal way to simulate the motion of small self-propelled objects since the interaction between the solvent and the motor can be specified and hydrodynamic effects are taken into account automatically. It has been used to investigate the self-propelled motion of swimmers composed of linked beads that undergo non-time-reversible cyclic motion [116] and chemically powered nanodimers [117]. The chemically powered nanodimers can serve as models for the motions of the bimetallic nanodimers discussed earlier. The nanodimers are made from two spheres separated by a fixed distance R dissolved in a solvent of A and B molecules. One dimer sphere (C) catalyzes the irreversible reaction A + C B I C, while nonreactive interactions occur with the noncatalytic sphere (N). The nanodimer and reactive events are shown in Fig. 22. The A and B species interact with the nanodimer spheres through repulsive Lennard-Jones (LJ) potentials in Eq. (76). The MPC simulations assume that the potentials satisfy Vca = Vcb = Vna, with c.,t and Vnb with 3- The A molecules react to form B molecules when they approach the catalytic sphere within the interaction distance r < rc. The B molecules produced in the reaction interact differently with the catalytic and noncatalytic spheres. 

A typical Lennard-Jones (6-12) function is plotted in Fig. 9. Often, the second team on the right-hand side of Eq. (100) is added to represent an effective potential function, viz. 

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.

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

The QM/MM interactions (Eqm/mm) are taken to include bonded and non-bonded interactions. For the non-bonded interactions, the subsystems interact with each other through Lennard-Jones and point charge interaction potentials. When the electronic structure is determined for the QM subsystem, the charges in the MM subsystem are included as a collection of fixed point charges in an effective Hamiltonian, which describes the QM subsystem. That is, in the calculation of the QM subsystem we determine the contributions from the QM subsystem (Eqm) and the electrostatic contributions from the interaction between the QM and MM subsystems as explained by Zhang et al. [13],... 

Fig. 3.8. Transition matrix of move proposal probabilities for the Lennard-Jones fluid at p = 0.88 and N = 110. The energy range of —700 to —500 in Lennard-Jones units has been discretized into 100 bins. Due to the adjustment of the random displacement moves to achieve 50% acceptance, the transition probabilities are highly banded. The tuned moves change the potential energy by only a small amount, and as a result, each energy level is effectively only connected to a few neighbors...

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

---