Quantum-mechanical explanation of intermolecular interactions


In the natural sciences, an intermolecular force is an attraction between two molecules or atoms. They occur from either momentary interactions between molecules or permanent electrostatic attractions between dipoles. They can be explained using a simple phenomenological approach, or using a quantum mechanical approach.

Perturbation theory

Hydrogen bonding, dipole-dipole interactions, and London forces are most naturally accounted for by Rayleigh-Schrödinger perturbation theory. In this theory—applied to two monomers A and B—one uses as unperturbed Hamiltonian the sum of two monomer Hamiltonians,
In the present case the unperturbed states are products
with and

Supermolecular approach

The early theoretical work on intermolecular forces was invariably based on RS-PT and its antisymmetrized variants. However, since the beginning of the 1990s it has become possible to apply standard quantum chemical methods to pairs of molecules. This approach is referred to as the supermolecule method. In order to obtain reliable results one must include electronic correlation in the supermolecule method, and take care of the basis set superposition error. This is the effect that the atomic orbital basis of one molecule improves the basis of the other. Since this improvement is distance dependent, it easily gives rise to artifacts.

Exchange

The monomer functions ΦnA and ΦmB are antisymmetric under permutation of electron coordinates, but the product states are not antisymmetric under intermolecular exchange of the electrons. An obvious way to proceed would be to introduce the intermolecular antisymmetrizer. But, as already noticed in 1930 by Eisenschitz and London, this causes two major problems. In the first place the antisymmetrized unperturbed states are no longer eigenfunctions of H, which follows from the non-commutation
In the second place the projected excited states
become linearly dependent and the choice of a linearly independent subset is not apparent. In the late 1960s the Eisenschitz-London approach was revived and different rigorous variants of symmetry adapted perturbation theory were developed.. The different approaches shared a major drawback: they were very difficult to apply in practice. Hence a somewhat less rigorous approach was introduced: apply ordinary RS-PT and introduce the intermolecular antisymmetrizer at appropriate places in the RS-PT equations. This approach leads to feasible equations, and, when electronically correlated monomer functions are used, weak symmetry forcing is known to give reliable results.
The first-order energy including exchange is in almost all symmetry-adapted perturbation theories given by the following expression
The main difference between covalent and non-covalent forces is the sign of this expression. In the case of chemical bonding this interaction is attractive and responsible for large bonding energies—on the order of a hundred kcal/mol. In the case of intermolecular forces between closed shell systems, the interaction is strongly repulsive and responsible for the "volume" of the molecule.
Roughly speaking, the exchange interaction is proportional to the differential overlap between Φ0A and Φ0B. Since the wavefunctions decay exponentially as a function of distance, the exchange interaction does too. Hence the range of action is relatively short, which is why exchange interactions are referred to as short range interactions.

Electrostatic interactions

By definition the electrostatic interaction is given by the first-order Rayleigh-Schrödinger perturbation energy :
Let the clamped nucleus α on A have position vector Rα,
then its charge times the Dirac delta function, Zα δ, is the charge density of this nucleus.
The total charge density of monomer A is given by
with the electronic charge density given by an integral over nA − 1 primed electron coordinates:
An analogous definition holds for the charge density of monomer B. It can be shown that the first-order quantum mechanical expression can be written as
which is nothing but the classical expression for the electrostatic interaction between two charge distributions. This shows that the first-order RS-PT energy is indeed equal to the electrostatic interaction between A and B.

Multipole expansion

At present it is feasible to compute the electrostatic energy without any further approximations other than those applied in the computation of the monomer wavefunctions. In the past this was different and a further approximation was commonly introduced: VAB was expanded in a series in inverse powers of the intermolecular distance R. This yields the multipole EXPANSION of the electrostatic energy. Since its concepts still pervade the theory of intermolecular forces, it will be presented here. In this article the following expansion is proved
with the Clebsch-Gordan series defined by
and the irregular solid harmonic is defined by
The function YL,M is a normalized spherical harmonic, while
and are spherical multipole moment operators. This expansion is manifestly in powers of 1/RAB.
Insertion of this expansion into the first-order expression gives a very similar expansion for the electrostatic energy, because the matrix element factorizes,
with the permanent multipole moments defined by
We see that the series is of infinite length, and, indeed, most molecules have an infinite number of non-vanishing multipoles. In the past, when computer calculations for the permanent moments were not yet feasible, it was common to truncate this series after the first non-vanishing term.
Which term is non-vanishing, depends very much on the symmetry of the molecules constituting the dimer. For instance, molecules with an inversion center such as a homonuclear diatomic, or an organic molecule like ethene do not possess a permanent dipole moment, but do carry a quadrupole moment.
Molecules such a hydrogen chloride and water lack an inversion center and hence do have a permanent dipole. So, the first non-vanishing electrostatic term in, e.g., the N2—H2O dimer, is the lA=2, lB=1 term. From the formula above follows that this term contains the irregular solid harmonic of order L = lA + lB = 3, which has an R−4 dependence. But in this dimer the quadrupole-quadrupole interaction is not unimportant either, because the water molecule carries a non-vanishing quadrupole as well.
When computer calculations of permanent multipole moments of any order became possible, the matter of the convergence of the multipole series became urgent. It can be shown that, if the charge distributions of the two monomers overlap, the multipole expansion is formally divergent.

Ionic interactions

It is debatable whether ionic interactions are to be seen as intermolecular forces, some workers consider them rather as special kind of chemical bonding.
The forces occur between charged atoms or molecules. Ionic bonds are formed when the difference between the electron affinity of one monomer and the ionization potential of the other is so large that electron transfer from the one monomer to the other is energetically favorable. Since a transfer of an electron is never complete there is always a degree of covalent bonding.
Once the ions are formed, the interaction between them can be seen as a special case of multipolar attraction, with a 1/RAB distance dependence. Indeed, the ionic interaction is the electrostatic term with lA = 0 and lB = 0. Using that the irregular harmonics for L = 0 is simply
and that the monopole moments and their Clebsch-Gordan coupling are
we recover—as to be expected—Coulomb's law
For shorter distances, where the charge distributions of the monomers overlap, the ions will repel each other because of inter-monomer exchange of the electrons.
Ionic compounds have high melting and boiling points due to the large amount of energy required to break the forces between the charged ions. When molten they are also good conductors of heat and electricity, due to the free or delocalized ions.
Writing
and similarly for B, we get the well-known expression
As a numerical example we consider the HCl dimer depicted above. We assume that the left molecule is A and the right B, so that the z-axis is along the molecules and points to the right. Our convention of the dipole moment is such that it points from negative to positive charge.
Note parenthetically that in organic chemistry the opposite convention is used. Since organic chemists hardly ever perform vector computations with dipoles, confusion hardly ever arises. In organic chemistry dipoles are mainly used as a measure of charge separation in a molecule.
So,
The value of μHCl is 0.43, so that at a distance of 10 bohr the dipole-dipole attraction is −3.698 10−4 hartree.
If one of the molecules is neutral and freely rotating, the total electrostatic interaction energy becomes zero.. In gases and liquids molecules are not rotating completely freely—the rotation is weighted by the Boltzmann factor exp, where k is the Boltzmann constant and T the absolute temperature.
It was first shown by John Lennard-Jones that the temperature-averaged dipole-dipole interaction is
Since the averaged energy has an R−6 dependence, it is evidently much weaker than the unaveraged one, but it is not completely zero. It is attractive, because the Boltzmann weighting favors somewhat the attractive regions of space. In HCl-HCl we find for T = 300 K and RAB = 10 bohr the averaged attraction −62 J/mol, which shows a weakening of the interaction by a factor of about 16 due to thermal rotational motion.

Anisotropy and non-additivity of intermolecular forces

Consider the interaction between two electric point charges at position
and. By
Coulomb's law the interaction potential
depends only on the distance between
the particles. For molecules this is different. If we see a molecule as
a rigid 3-D body, it has 6 degrees of freedom. The
interaction energy of two molecules in isotropic and
homogeneous space is in general a function of 2×6−6=6 degrees of freedom
.
The analytic description of the interaction of two arbitrarily shaped rigid
molecules requires therefore 6 parameters. The
fact that the intermolecular interaction depends on the orientation of the molecules
is expressed by stating that the potential is anisotropic. Since point charges are
by definition spherical symmetric, their interaction is isotropic. Especially in the
older literature, intermolecular interactions are regularly assumed to be
isotropic.
Consider three arbitrary point charges at distances
r12, r13, and r23 apart.
The total interaction U is additive; i.e., it is the sum
Again for molecules this can be different. Pretending that the
interaction depends on distances only—but see above—the
interaction of three molecules takes in general the form
where is a non-additive three-body interaction.
Such an interaction can be caused by exchange interactions, by induction, and by
dispersion.