VirialVLE:Background:Quadrupole Term

From Wiketomica
Basic Layout
Run Simulation
Etomica Modules

For a small molecule such as CO2 molecule it is possible to perform a fairly high–level ab–initio computation of the quadrupole moment tensor elements . As an example, for an isolated CO2 molecule at zero Kelvin a calculation at the ab–initio CCSD/cc-pVTZ level of theory (including geometry optimization) with the Gaussian 03 code yielded the following result:

Quadrupole moment (field-independent basis, debyeÅ):
XX= -14.6739 YY= -14.6739 ZZ= -19.8172
XY= 0.0000 XZ= 0.0000 YZ= 0.0000

for the CO2 molecule aligned with the z–axis and the carbon atom at the coordinate origin. The optimized C=O bond length was 1.159 Å which is very close to the experimental bond length of 1.162 Å obtained by Herzberg (1966) as retreived from the CCCBDB. (accessed in Feb. 2008, it was not explicitly stated whether the experiment refers to or , we assume it is the equilibrium value ).

The Gaussian code outputs the Cartesian quadrupole moments, i.e. instead of Eq. (5) it calculates and prints

..... (Equation 8)

This q–tensor is not traceless. The conversion to Eq. (5) is done via

..... (Equation 9)


For CO2, because of rotational symmetry we have , therefore . From the computed data above, we obtain = -5.1433. For the linear CO2 molecule the magnitude of this value is the quadrupole moment Q, the numerical parameter that is adjusted in the module Virial along with the LJ parameters and to fit the experimental B(T) data. The quadrupole moment is also one of the parameters that can be supplied as input for the VLE module to simulate the phase equilibrium. The units of Q are usually given in debyeÅ, i.e. a dipole moment units multiplied with a length unit. Both are “atomic scale” units appropriate to describe the quadrupole moment of a single molecule.