Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

First-order properties

Electric multipole moments

The permanent multipole moments characterize the charge distribution in a molecule. The zeroth--order (monopole) moment corresponds to the total charge, the first moment to the dipole, the second to the quadrupole, and so on. The multipole moments can be computed as expectation values of the multipole moment operators discussed elsewhere. As for the operators, we can use the primitive or traceless forms for the quadrupole and higher-order moments. The traceless multipole moments describe the deviation from an isotropic charge distribution. For instance, the traceless quadrupole moment vanishes for a spherical charge distribution when the center of the sphere coincides with the origin, since the transformation from the primitive (traced) quadrupole moment qαβq_{\alpha\beta} to the traceless quadrupole moment Θαβ\Theta_{\alpha\beta} eliminates the isotropic contribution. The traceless multipole moment of order nn is defined by 2n+12n+1 components, whereas the traced tensor involves (n+2)(n+1)/2(n+2)(n+1)/2 independent quantities. A tabulation of non-zero symmetry-independent elements of multipole moments (μ\mu = dipole, Θ\Theta = quadrupole, Ω\Omega = octupole and Φ\Phi = hexadecapole) for different point groups is given elsewhere Kielich (1972) .

In general, when the multipole expansion of the potential is used to describe the influence of a static electric field, only the first non-zero moment is unambiguously defined. The next,higher-order multipole moments depend on the choice of origin, and when computed with respect to the origin of the coordinate system, they depend on the position of the molecule in space. This can be seen from the form of the the electric dipole operator

μ^α=−e∑ir^iα+e∑KZKRKα ,\hat{\mu}_\alpha = -e \sum_{i} {\hat{r}}_{i\alpha} +e \sum_{K}Z_K R_{K \alpha}\ ,

where a shift of the chosen common origin for the expansion of the potential by a vector D{\bf D} for a molecule with the total charge QQ leads to a change of QDQ{\bf D} in the dipole moment; similarly, for a neutral molecule with a non-zero dipole moment, the quadrupole moment will change.

Dipole and quadrupole moments

From the multipole expansion of the energy in the presence of a static electric field given previously we find that the α−\alpha-component of the static permanent dipole moment is defined as the electric field derivative of E(F)E(\mathbf{F})

μα=−∂E(F)∂Fα∣F=0.    \mu_\alpha = -\left. \frac{\partial E\left({\mathbf{F}} \right)}{\partial F_\alpha} \right|_{{\mathbf{F}} = {\bf 0}}. \;\;
(1)

For wave functions that fulfill the Hellmann-Feynman theorem we can compute the dipole moment as the expectation value

μα=⟨0∣μ^α∣0⟩,    \mu_\alpha = \langle 0 \mid \hat{\mu}_\alpha \mid 0 \rangle, \;\;

where μ^α\hat\mu_\alpha is the α−\alpha-component of the dipole moment operator defined just above. Obviously, the expectation value expression is correct for the exact wave function. Furthermore, for most approximate wave functions it is much easier to compute the expectation value than to calculate the energy derivative. However, the results differ when the method used does not fulfill the Hellmann-Feynman theorem.
In this case, one can either derive the analytic expressions which reproduce the correct definition of the electric field derivative of E(F)E(\mathbf{F}) or use the finite field method and numerical differentiation.

The expectation value of Θ^{\hat{\Theta}} - the electric quadrupole operator in its traceless form

Θ^αβ=−e2∑i(3r^iαr^iβ−r^i2δαβ)+e2∑KZK(3RKαRKβ−RK2δαβ).{\hat{\Theta}}_{\alpha\beta} = -\frac{e}{2} \sum_{i} (3{\hat{r}}_{i\alpha}{\hat{r}}_{i\beta}-{\hat{r}}_{i}^{2}\delta_{\alpha \beta}) + \frac{e}{2} \sum_{K} Z_{K} { ( 3 R_{K \alpha} R_{K \beta} - {R}_{K}^{2} \delta_{\alpha \beta} ) }.

as defined following Buckingham Buckingham (1967) ,

Θαβ=⟨0∣Θ^αβ∣0⟩=−3∂E(∇F)∂Fαβ \Theta_{\alpha\beta} = \langle 0 \mid {\hat{\Theta}}_{\alpha \beta} \mid 0 \rangle = -3 \frac{\partial E \left(\nabla {\mathbf{F}} \right)} {\partial F_{\alpha\beta}} \,

can be related to the expectation values of second moments, qαβq_{\alpha\beta} = ⟨0∣q^αβ∣0⟩\langle 0| \hat{q}_{\alpha\beta} |0 \rangle , where the latter are defined as

q^αβ=−e∑ir^iαr^iβ+e∑KZKRKαRKβ.\hat{q}_{\alpha\beta} = -e \sum_i {\hat{r}}_{i\alpha}{\hat{r}}_{i\beta} + {e} \sum_{K} Z_{K} R_{K \alpha} R_{K \beta}.

by a simple linear transformation.

Both Θ\mathbf{\Theta} and q\mathbf{q} are rank two, symmetric tensors.

For molecules with an nn-fold rotation axis with n≥3n \ge 3 :

Θ≡Θzz=−2Θyy=−2Θxx,\Theta \equiv \Theta_{zz} = -2 \Theta_{yy} = -2 \Theta_{xx},
(2)

where zz coincides with the rotation axis. In experiment, two diagonal components of Θ\mathbf{\Theta} are measured.

Units conversion factors for the electric dipole and quadrupole moments are, respectively

1 au of μ\mu = ea0e a_0 = 8.47835×10−30 C m8.47835 \times 10^{-30}\ \mathrm{C\,m} (SI) = 2.5417×10−18 statvolt cm22.5417 \times 10^{-18}\ \mathrm{statvolt}\,\mathrm{cm}^2 (esu)

1 au of Θ\Theta = ea02e a_0^2 = 4.48655×10−40 C m24.48655 \times 10^{-40}\ \mathrm{C\,m}^2 (SI) = 1.3450×10−26 statvolt cm31.3450 \times 10^{-26}\ \mathrm{statvolt}\,\mathrm{cm}^3 (esu)

The unit of Debye is 10−18 statvolt cm210^{-18}\ \mathrm{statvolt}\,\mathrm{cm}^2 and is often used for the electric dipole moment.

References
  1. Kielich, S. (1972). General molecular theory and electric field effects in isotopic dielectrics. In M. Davies (Ed.), Specialist Periodical Report, Dielectric and Related Molecular Processes (Vol. 1, p. 192). Chem. Soc.
  2. Buckingham, A. D. (1967). Permanent and Induced Molecular Moments and Long-Range Intermolecular Forces. Adv. Chem. Phys., 12, 107. 10.1002/9780470143582.ch2
Spectra and properties
Molecular properties
Spectra and properties
Linear Response