VirialVLE:Background:Model Potentials

From Wiketomica
Basic Layout
Run Simulation
Etomica Modules

The interaction potential energy surface (PES) between two molecules depends on their distance and relative orientation. These days the PES can be calculated ab–initio by using high–quality electronic structure methods. However, the PES may have a complicated form. Simple models have been devised that capture the general features of intermolecular interaction. With suitable adjustments and refinements, these potentials are also capable of yielding good agreement with experiment when used in molecular simulations. An example for a simple interaction model potential U between two neutral molecules is the Lennard–Jones (LJ, or 6/12) potential:

..... (Equation 3)

Here, quantifies the strength of the interaction, R is the distance between the molecules, and is the location of the minimum, as shown in the plot below:

VLE MinLoc.jpg

Another way of writing this potential is :

..... (Equation 4)

where is the separation where the potential is zero. In the plot above, the potential U exhibits an attraction between two molecules at larger distances and a strong repulsion at short distances, which is typical for neutral molecules. In the VLE simulation module we attempt to obtain reasonable agreement with experimental data. In this case, for a molecule such as CO2 the simple LJ potential is not accurate enough. The reason is that the atoms in CO2 carry partial charges, i.e. the charge distribution in the molecule is non–spherical. As a consequence, the interaction potential depends not only on the distance R but also on the relative orientation of the two molecules and has a somewhat different distance dependence. For molecules without permanent dipole moments (example: CO2) the lowest–order electrostatic term relevant here is the quadrupole–quadrupole interaction. Given a charge density of the molecule, the elements of a traceless electric quadrupole moment tensor can be defined as [1]

..... (Equation 5)

where . "Traceless" means that . For a molecule, we have two contributions: one from the electron density (counts negative), and one from the nuclei A with charges , positions and point–charge density . The quadrupole generates an electric potential at a position R (with distance R) of :

..... (Equation 6)

along with a non–uniform electric field . The field is seen to decay with the 4th inverse power of R at large distances R. Another quadrupolar molecule can interact with the gradient components of this electric field [2]; the interaction energy U is proportional to at large distances and therefore provides an important contribution to the PES of neutral unpolar molecules such as CO2. Note that without a quadrupole–quadrupole term we would be missing some longer–range interactions in the PES; the term of the LJ potential does not reach as far out.