Cylindrical multipole moments explained

Cylindrical multipole moments are the coefficients in a series expansion of a potential that varies logarithmically with the distance to a source, i.e., as

lnR

. Such potentials arise in the electric potential of long line charges, and the analogous sources for the magnetic potential and gravitational potential.

For clarity, we illustrate the expansion for a single line charge, then generalize to an arbitrary distribution of line charges. Through this article, the primed coordinates such as

(\rho\prime,\theta\prime)

refer to the position of the line charge(s), whereas the unprimed coordinates such as

(\rho,\theta)

refer to the point at which the potential is being observed. We use cylindrical coordinates throughout, e.g., an arbitrary vector

r

has coordinates

(\rho,\theta,z)

where

\rho

is the radius from the

z

axis,

\theta

is the azimuthal angle and

z

is the normal Cartesian coordinate. By assumption, the line charges are infinitely long and aligned with the

z

axis.

Cylindrical multipole moments of a line charge

The electric potential of a line charge

λ

located at

(\rho',\theta')

is given by\Phi(\rho, \theta) = \frac \ln R = \frac \ln \left| \rho^ + \left(\rho' \right)^ - 2\rho\rho'\cos (\theta - \theta') \right|where

R

is the shortest distance between the line charge and the observation point.

By symmetry, the electric potential of an infinite line charge has no

z

-dependence. The line charge

λ

is the charge per unit length in the

z

-direction, and has units of (charge/length). If the radius

\rho

of the observation point is greater than the radius

\rho'

of the line charge, we may factor out

\rho2

\Phi(\rho, \theta) =\frac \left\and expand the logarithms in powers of

(\rho'/\rho)<1

\Phi(\rho, \theta) =\frac \left\which may be written as\Phi(\rho, \theta) =\frac \ln \rho +\frac \sum_^ \fracwhere the multipole moments are defined as\beginQ &= \lambda,\\C_k &= \frac \left(\rho' \right)^k \cos k\theta', \\S_ &= \frac \left(\rho' \right)^k \sin k\theta'.\end

Conversely, if the radius

\rho

of the observation point is less than the radius

\rho'

of the line charge, we may factor out

\left(\rho'\right)2

and expand the logarithms in powers of

(\rho/\rho')<1

\Phi(\rho, \theta) =\frac \left\which may be written as\Phi(\rho, \theta) =\frac \ln \rho' +\frac \sum_^ \rho^ \left[I_{k} \cos k\theta + J_{k} \sin k\theta \right] where the interior multipole moments are defined as\beginQ &= \lambda, \\I_k &= \frac \frac, \\J_k &= \frac \frac.\end

General cylindrical multipole moments

The generalization to an arbitrary distribution of line charges

λ(\rho',\theta')

is straightforward. The functional form is the same \Phi(\mathbf) = \frac \ln \rho + \frac \sum_^ \fracand the moments can be written\beginQ &= \int d\theta' \, d\rho' \, \rho' \lambda(\rho', \theta') \\C_k &= \frac \int d\theta' \, d\rho' \left(\rho'\right)^ \lambda(\rho', \theta') \cos k\theta' \\S_k &= \frac \int d\theta' \, d\rho' \left(\rho'\right)^ \lambda(\rho', \theta') \sin k\theta'\endNote that the

λ(\rho',\theta')

represents the line charge per unit area in the

(\rho-\theta)

plane.

Interior cylindrical multipole moments

Similarly, the interior cylindrical multipole expansion has the functional form\Phi(\rho, \theta) =\frac \ln \rho' +\frac \sum_^ \rho^ \left[I_{k} \cos k\theta + J_{k} \sin k\theta \right] where the moments are defined\beginQ &= \int d\theta' \, d\rho' \, \rho' \lambda(\rho', \theta') \\I_ &= \frac\int d\theta' \, d\rho' \frac\lambda(\rho', \theta') \\J_ &= \frac\int d\theta' \, d\rho' \frac\lambda(\rho', \theta') \end

Interaction energies of cylindrical multipoles

A simple formula for the interaction energy of cylindrical multipoles (charge density 1) with a second charge density can be derived. Let

f(r\prime)

be the second charge density, and define

λ(\rho,\theta)

as its integral over z\lambda(\rho, \theta) = \int dz \, f(\rho, \theta, z)

The electrostatic energy is given by the integral of the charge multiplied by the potential due to the cylindrical multipolesU = \int d\theta \, d\rho \, \rho \, \lambda(\rho, \theta) \Phi(\rho, \theta)

If the cylindrical multipoles are exterior, this equation becomesU = \frac \int d\rho \, \rho \, \lambda(\rho, \theta) \ln \rho + \frac \sum_^ \int d\theta \, d\rho \left[C_{1k} \frac{\cos k\theta}{\rho^{k-1}} + S_{1k} \frac{\sin k\theta}{\rho^{k-1}}\right] \lambda(\rho, \theta)where

Q1

,

C1k

and

S1k

are the cylindrical multipole moments of charge distribution 1. This energy formula can be reduced to a remarkably simple formU = \frac \int d\rho \, \rho \, \lambda(\rho, \theta) \ln \rho + \frac \sum_^ k \left(C_ I_ + S_ J_ \right)where

I2k

and

J2k

are the interior cylindrical multipoles of the second charge density.

The analogous formula holds if charge density 1 is composed of interior cylindrical multipolesU = \frac \int d\rho \, \rho \, \lambda(\rho, \theta) + \frac \sum_^ k \left(C_ I_ + S_ J_ \right)where

I1k

and

J1k

are the interior cylindrical multipole moments of charge distribution 1, and

C2k

and

S2k

are the exterior cylindrical multipoles of the second charge density.

As an example, these formulae could be used to determine the interaction energy of a small protein in the electrostatic field of a double-stranded DNA molecule; the latter is relatively straight and bears a constant linear charge density due to the phosphate groups of its backbone.

See also