## Liquid-State Physical Chemistry: Fundamentals, Modeling, and Applications (2013)

### 3. Basic Energetics: Intermolecular Interactions

### 3.2. Electrostatic Interaction

The *Coulomb energy* *w* of a charge *q*_{1} in the presence of another charge *q*_{2} is

(3.4)

Here, *ϕ*_{2} represents the *potential* associated with the *force* between the charges *q*_{1} and *q*_{2} with *r* = |** r**| their distance and

*ε*

_{0}the permittivity of vacuum. Molecule 1 can be represented by a charge distribution

__, and the potential__

^{4)}*ϕ*

_{2}due to the charge distribution of molecule 2 by where

*s*= |

_{j}

*s**| = |*

_{j}**−**

*r*

*r**| is the distance of charge*

_{j}*q*to position

_{j}**(**

*r*__Figure 3.2__). The total interaction energy with

*s*= |

_{jk}

*s**| = |*

_{jk}

*r**−*

_{j}

*r**|. After expanding the function*

_{k}*ϕ*

_{2}, evaluating all the derivatives and collecting terms, a procedure with considerable understatement often called “a lengthy but straightforward calculation” (see

__Justification 3.1__), one obtains the desired result. Note that, since this process takes place at constant

*T*and

*V*, the interaction

*W*

_{12}represents the Helmholtz energy.

** Figure 3.2** The interaction

*ρ*

_{1}

*ϕ*

_{2}of two charge distributions

*ρ*

_{1}and

*ρ*

_{2}with common origin O.

Justification 3.1: The multipole expansion*

To find the proper expression for the electrostatic interaction, let us consider an arbitrary charge distribution *ρ* of point charges *q _{j}* at position

*r**from the origin O located at the center of mass (*

_{j}__Figure 3.2__). The potential energy

*ϕ*at a certain point P located at

**outside the sphere containing all charges is**

*r*where *s _{j}* = |

*s**| = |*

_{j}**−**

*r*

*r**| is the distance of charge*

_{j}*q*to P. Developing 1/

_{j}*s*in a Taylor series with respect to

_{j}

*r**, we may write 1/*

_{j}*s*= 1/

_{j}*r*+

*r**(∇1/*

_{j}*r*)

_{O}+ ··· and therefore (using direct notation,

__Appendix B__)

__(3.5)__

where *r* = |** r**| is the distance of the point P to the origin O of

*ρ*. In the second line we defined

*q*= Σ

*as the*

_{j}q_{j}*total charge*,

**= Σ**

*μ*

_{j}q_{j}

*r**as the*

_{j}*dipole moment*, and

**= ½Σ**

*Q*

_{j}q_{j}

*r*

_{j}

*r**as the*

_{j}*quadrupole moment*. Moreover, we abbreviate 1/

*r*by

*ϕ*(0), [∇(1/

*r*)]

_{O}by

*φ′*(0), etc. The minus sign in the last step for the second, fourth, and so on terms arises because the differentiation is now made at P, and not at the origin of the potential O. This follows since for an arbitrary vector

**it holds that ∇(1/**

*x**x*) = −

**/**

*x**x*

^{3}and therefore [∇(1/

*r*)]

_{P}

**=**−[∇(1/

*r*)]

_{O}.

__Equation (3.5)__is usually called the

*multipole expansion*of the potential. The

*interaction W*of a point charge

*q*at

**r**in the potential

*ϕ*of a charge distribution is

*W*=

*qϕ*(

**). Hence, to calculate the interaction energy between two charge distributions**

*r**ρ*

_{1}and

*ρ*

_{2}, separated by

**(**

*r*__Figure 3.2__), we express this energy as the energy of in the field of . The interaction energy

*W*

_{12}=

*ρ*

_{1}

*ϕ*

_{2}is then

If we take the center of mass P of *ρ*_{1} as common origin, the potential *ϕ*_{2} due to *ρ*_{2} is given by the last line of __Eq. (3.5)__, and substitution in the previous equation results in

__(3.6)__

This is the general expression for the electrostatic interaction energy expressed in terms of multipole moments of the charge distributions with respect to their own center of mass [1]. For further evaluation, we need again ∇(1/*r*) = −** r**/

*r*

^{3}and also ∇

^{2}(1/

*r*) = 3

**/**

*rr**r*

^{5}–

**/**

*I**r*

^{3}with

**the unit tensor, easily derived using**

*I**r*

^{2}=

*x*

^{2}+

*y*

^{2}+

*z*

^{2}, and leading to

__Eq. (3.7)__.

The final result reads

__(3.7)__

or explicitly (see __Figure 3.3__ for notation, introducing the abbreviation **ω** = *θ*,*ϕ*)

__(3.8)__

** Figure 3.3** The interaction at a distance

*r*of a point charge

*q*

_{1}with a dipole

*μ*

_{2}in which the charges

*q*are separated by a distance

*δ*(a) and the interaction of a dipole with a dipole at distance

*r*(b).

In this expression is the *charge* of molecule 1, and is the *dipole moment* with *μ*_{1} = |*μ*_{1}|.

The first term (proportional to *q*_{1}, *q*_{2} and 1/*r*) is just the *Coulomb interaction* (more precisely the charge–charge Coulomb interaction) between the molecules at distance *r*. At long distance *r* this interaction term will dominate, but clearly it is zero if one of the molecules is neutral.

The second term (proportional to *q*_{1}, *μ*_{2}, *q*_{2}, *μ*_{1} and 1/*r*^{2}; __Figure 3.3__) is the interaction between the charge of molecule 2 with the dipole moment of molecule 1 and the interaction between the charge of molecule 1 with the dipole moment of molecule 2. This is the *charge–dipole* (*Coulomb*)* interaction*. The charge–dipole interaction decreases as 1/*r*^{2}, as compared to 1/*r* for the charge–charge Coulomb interaction. A simple way to obtain this result is dealt with in __Problem 3.2__. Because the charge–dipole interaction is orientation-dependent and molecules move generally rather rapidly, an orientational average (overbar) must be made (we will learn how to do that in __Chapter 5__). For the moment we just quote the result, valid for sufficiently high temperature,

(3.9)

The third term (proportional to *μ*_{1}, *μ*_{2} and 1/*r*^{3}; __Figure 3.3__) represents the *dipole–dipole* (*Coulomb*)* interaction* between the dipole moments of molecule 1 and molecule 2, and is often addressed as the *Keesom interaction* (it also includes the charge–quadrupole interaction, but we neglect that interaction). We deal with this interaction in a simplified way in __Problem 3.3__. This interaction is also orientation-dependent and for the orientational average one obtains, again for sufficiently high temperature and using **ω** = (*θ*,*ϕ*) with *θ* and *ϕ* the orientation angles

(3.10)

For neutral molecules the dipole–dipole interaction is the leading term. Note (again) that the expressions given represent the Helmholtz energies and that the internal energy expressions *U* are given by *U* = 2*W* (see Problem 2.6).

Problem 3.1

Show that for a charged molecule the dipole moment depends on the choice of the origin but is independent of this choice for a neutral molecule.

Problem 3.2: Charge–dipole interaction

A simple model to estimate the charge–dipole interaction is shown in __Figure 3.3__. The distances between the charge *q*_{1} of molecule 1 and the charges of the dipole moment of molecule 2, +*q*_{2} and −*q*_{2}, respectively, are given by

Show, if it is assumed that the distance *r* is much larger than the distance *δ* by using Coulomb's law and the binomial expansion, that the interaction for fixed direction is where . How large should *r* be as compared to *δ* in order to have an error less than 5%?

Problem 3.3: Dipole–dipole interaction

The dipole–dipole interaction can be estimated in a similar way as for the charge–dipole problem (__Figure 3.3__). Show that the expression for the in-line orientation energy (orientations parallel to the connection line) is

Here, ± indicates whether the orientation is in the same way (attraction) or in the opposite way (repulsion). Also show that the expression for the parallel orientation energy (orientations perpendicular to the connection line) is

again, dependent on whether the orientation is antiparallel (attraction) or parallel (repulsion). How large should *r* be as compared to *δ* in order to have an error less than 5%?

Problem 3.4*

Verify __Equation (3.6)__. Alternatively, derive the expression for the electrostatic interaction up to second order using a double Taylor expansion for with *s _{ij}* = |

**+**

*r*

*r**−*

_{j}

*r**| using the complete expression .*

_{i}