Electrostatic potential of a uniformly charged disk through Green’s theorem

Abstract

All existing derivations of the electrostatic potential of a uniformly charged disk are technically rather involved. In an old and now almost forgotten publication, Duffin and McWhirter proposed a method for calculating the electrostatic potentials of planar bodies by a skillful application of Green’s theorem. It is shown that this method significantly simplifies the problem of calculating the electrostatic potential of a uniformly charged disk and makes it an almost trivial task.

Keywords: Electrostatic potential; Uniformly charged disk; Green’s theorem.

1 Introduction

Calculating the electrostatic potential of a uniformly charged disk (or its gravitational analogue) is a classic problem, and there are a variety of ways to solve this problem [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. However, all existing solution methods are quite complex from a technical point of view, which perhaps explains why this interesting electrostatic problem is not reflected in introductory textbooks on classical electrodynamics.

When the observation point is in the plane of the disk, in [4, 5] the electrostatic potential was calculated by skillful use of Green’s theorem [11, 12]

$$\iint\limits_{D}\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dxdy=\oint\limits_{\partial D}\left(Pdx+Qdy\right).$$ (1)

This neat trick replaces the double integral over the surface of the disk with a line integral over the boundary of the disk and really makes calculating the electrostatic potential much easier. However, it cannot be extended in the form in which it is used in [5] to observation points located outside the plane of the disk.

It is well known that Green’s theorem (1) is equivalent to the two-dimensional divergence theorem [11, 12]

$$\iint\limits_{D}\left(\frac{\partial A_{x}}{\partial x}+\frac{\partial A_{y}}{\partial y}\right)dxdy=\oint\limits_{\partial D}{\bf{A}}\cdot{\bf{n}}\,dl.$$ (2)

Indeed, it is clear from Fig.1 that ${\bf n}\,dl=dy{\bf i}-dx{\bf j}$, and if we choose ${\bf A}=(Q,-P)$ it is easy to see that (2) transforms into (1).

Figure 1: A planar region $D$ bounded by a closed curve $\partial D$, the tangent vector $d{\bf{l}}$ and the outer normal vector ${\bf{n}}\,dl$ at a given point.

Although (1) and (2) are equivalent, Duffin and McWhirter [13] showed that the divergence theorem (2) is more useful and suggestive when used in electrostatic problems for planar charged bodies (they attribute the origin of the idea to Collie [14]). Crucial observation of [13] is that, if we take a two-variable scalar function ($z$ considered as a parameter)

$$\psi(x,y)=r-|z|\ln{(r+|z|)},\;\;\;r=\sqrt{(x-X)^{2}+(y-Y)^{2}+z^{2}},$$ (3)

then

$$\frac{\partial^{2}\psi}{\partial x^{2}}+\frac{\partial^{2}\psi}{\partial y^{2}}=\frac{r(r+|z|)-(x-X)^{2}}{r(r+|z|)^{2}}+\frac{r(r+|z|)-(y-Y)^{2}}{r(r+|z|)^{2}}=\frac{1}{r}.$$ (4)

Therefore, taking

$${\bf A}={\bf\nabla}\psi=\left(\frac{\partial\psi}{\partial x},\frac{\partial\psi}{\partial y}\right),$$ (5)

in (2), we find that the electrostatic potential created by a uniformly charged planar region $D$ with a surface charge density $\sigma$ at the observation point $(X,Y,z)$ is given by

$$\phi(X,Y,z)=\frac{\sigma}{4\pi\epsilon_{0}}\iint\limits_{D}\frac{dxdy}{r}=\frac{\sigma}{4\pi\epsilon_{0}}\oint\limits_{\partial D}\nabla\psi\cdot{\bf n}\,dl.$$ (6)

Note that the absolute value for $z$ in (3) is needed to avoid logarithmic singularity when $z<0,x=X,y=Y$.

2 Calculation of the electrostatic potential

Let point $P$ in Fig.2 be the projection onto the disk of the point at which we want to find the electrostatic potential (observation point). The coordinate axes can be selected so that $P=(X,Y)=(\eta R,0)$, $R$ being the disk radius and $\eta\geq 0$. In view of the axial symmetry of the problem, to characterize the integration point on the disk, we choose polar coordinates $x=\rho\cos{\alpha},\,y=\rho\sin{\alpha}$. Then $r=\sqrt{(\eta R)^{2}-2\eta R\rho\cos{\alpha}+\rho^{2}+z^{2}}$ and

$$\left.{\bf n}\cdot{\bf\nabla}\psi\right|_{\rho=R}=\left.\frac{\partial\psi}{\partial\rho}\right|_{\rho=R}=\frac{1-\eta\cos{\alpha}}{\sqrt{1+\eta^{2}-2\eta\cos{\alpha}+z^{2}/R^{2}}+|z|/R}.$$ (7)

Figure 2: If $P$ is projection of the observation point to the disk plane, $D_{1}$ and $D_{2}$ subdomains of the disk contribute equally to the electrostatic potential.

From the symmetry it is clear that the $D_{1}$ and $D_{2}$ halves of the disk make equal contributions to the potential (see Fig.2). In a line integral, the linear segments of $D_{1}$ and $D_{2}$ go in opposite directions, so their contributions cancel each other (alternatively, the symmetry argument can be applied to the line integral itself, stating that the two halves of the circle contribute equally). Therefore, since $dl=Rd\alpha$,

$$\phi(\eta R,z)=\frac{\sigma R}{2\pi\epsilon_{0}}\int\limits_{0}^{\pi}\frac{1-\eta\cos{\alpha}}{\sqrt{1+\eta^{2}-2\eta\cos{\alpha}+z^{2}/R^{2}}+|z|/R}d\alpha=\frac{\sigma R}{2\pi\epsilon_{0}}\left[I_{1}-\frac{|z|}{R}\,I_{2}\right],$$ (8)

where

$$I_{1}=\int\limits_{0}^{\pi}\frac{(1-\eta\cos{\alpha})\sqrt{1+\eta^{2}-2\eta\cos{\alpha}+z^{2}/R^{2}}}{1+\eta^{2}-2\eta\cos{\alpha}}d\alpha,$$ (9)

and

$$I_{2}=\int\limits_{0}^{\pi}\frac{1-\eta\cos{\alpha}}{1+\eta^{2}-2\eta\cos{\alpha}}d\alpha.$$ (10)

The integral $I_{2}$ can be calculated by a standard substitution $t=\tan{\frac{\alpha}{2}}$ which gives

$$I_{2}=2\int\limits_{0}^{\infty}\frac{1-\eta+(1+\eta)t^{2}}{(1+t^{2})[(1-\eta)^{2}+(1+\eta)^{2}t^{2}]}dt=\int\limits_{0}^{\infty}\left[\frac{1}{1+t^{2}}+\frac{1+\eta}{1-\eta}\;\frac{1}{1+\left(\frac{1+\eta}{1-\eta}t\right)^{2}}\right]dt.$$ (11)

In the second integral, we make a substitution $s=\frac{1+\eta}{|1-\eta|}t$, where the absolute value is needed to have limits $(0,\infty)$ in $s$. Finally,

$$I_{2}=\frac{\pi}{2}\left[1+\frac{|1-\eta|}{1-\eta}\right]=\pi\theta(1-\eta),$$ (12)

where $\theta(x)=\left\{\begin{array}[]{c}1,x>0,\\ 0,x<0.\end{array}\right.$ is the Heaviside step function.

Let $\alpha=\pi-2\varphi$. Then (9) can be rewritten in the following way

$$I_{1}=2\int\limits_{0}^{\pi/2}\frac{(1+\eta-2\eta\sin^{2}{\varphi})\sqrt{(1+\eta)^{2}+z^{2}/R^{2}-4\eta\sin^{2}{\varphi}}}{(1+\eta)^{2}-4\eta\sin^{2}{\varphi}}\,d\varphi$$ (13)

Denoting temporarily $A=(1+\eta)^{2}+z^{2}/R^{2}-4\eta\sin^{2}{\varphi}$, we have

$$1+\eta-2\eta\sin^{2}{\varphi}=\frac{1}{2}\left[1-\eta^{2}+A-z^{2}/R^{2}\right],$$ (14)

and

$$I_{1}=\int\limits_{0}^{\pi/2}\frac{\left[1-\eta^{2}+A-z^{2}/R^{2}\right]A}{\left[A-z^{2}/R^{2}\right]\sqrt{A}}\;d\varphi.$$ (15)

But

$$\frac{\left[1-\eta^{2}+A-z^{2}/R^{2}\right]A}{A-z^{2}/R^{2}}=A+(1-\eta^{2})\frac{A}{A-z^{2}/R^{2}}=A+1-\eta^{2}+z^{2}/R^{2}\frac{1-\eta^{2}}{A-z^{2}/R^{2}}.$$ (16)

Therefore, $I_{1}=J_{1}+J_{2}+J_{3}$, where

$$J_{1}=\int\limits_{0}^{\pi/2}\sqrt{A}\,d\varphi=\sqrt{(1+\eta)^{2}+z^{2}/R^{2}}\int\limits_{0}^{\pi/2}\sqrt{1-k^{2}\sin^{2}{\varphi}}\;d\varphi,$$

$$J_{2}=(1-\eta^{2})\int\limits_{0}^{\pi/2}\frac{d\varphi}{\sqrt{A}}=\frac{1-\eta^{2}}{\sqrt{(1+\eta)^{2}+z^{2}/R^{2}}}\int\limits_{0}^{\pi/2}\frac{d\varphi}{\sqrt{1-k^{2}\sin^{2}{\varphi}}},$$ (17)

with

$$k^{2}=\frac{4\eta}{(1+\eta)^{2}+z^{2}/R^{2}},$$ (18)

and

$$J_{3}=\frac{z^{2}}{R^{2}}(1-\eta^{2})\int\limits_{0}^{\pi/2}\frac{d\varphi}{\sqrt{A}\,(A-z^{2}/R^{2})}=$$

$$\frac{1-\eta}{1+\eta}\,\frac{z^{2}/R^{2}}{\sqrt{(1+\eta)^{2}+z^{2}/R^{2}}}\int\limits_{0}^{\pi/2}\frac{d\varphi}{(1-n^{2}\sin^{2}{\varphi})\sqrt{1-k^{2}\sin^{2}{\varphi}}},$$ (19)

with

$$n^{2}=\frac{4\eta}{(1+\eta)^{2}}.$$ (20)

From the definitions of complete elliptic integrals of the first, second and third types [15] it is clear that the final expression for the electrostatic potential of a uniformly charged disk, valid for all observation points, takes the form

$$\eqalign{\;\;\;\;\;\phi(\eta R,z)=\frac{\sigma R}{2\pi\epsilon_{0}}\left[\frac{1-\eta^{2}}{\sqrt{(1+\eta)^{2}+z^{2}/R^{2}}}\,K(k)+\sqrt{(1+\eta)^{2}+z^{2}/R^{2}}\,E(k)+\right.\cr\left.\;\;\;\;\;\frac{1-\eta}{1+\eta}\frac{z^{2}/R^{2}}{\sqrt{(1+\eta)^{2}+z^{2}/R^{2}}}\,\Pi(n^{2},k)-\pi\,\theta(1-\eta)\,\frac{|z|}{R}\right].}$$ (21)

Because of the presence of the $\theta$-function, when calculating the electric field from (21), care must be taken not to get a spurious $\delta$-function singularity [16].

3 Concluding remarks

If $z=0$, then

$$\frac{\partial}{\partial x}\frac{x-X}{r}+\frac{\partial}{\partial y}\frac{y-Y}{r}=\frac{1}{r},$$ (22)

and we can use Green’s theorem (1) to get

$$\phi(X,Y,z)=\frac{\sigma}{4\pi\epsilon_{0}}\iint\limits_{D}\frac{dxdy}{r}=\frac{\sigma}{4\pi\epsilon_{0}}\oint\limits_{\partial D}\frac{(x-X)dy-(y-Y)dx}{r}.$$ (23)

Using this formula, Kondratyev in [4, 5] considered a number of electrostatic problems, including uniformly charged circular and elliptical disks. Many other interesting results on potential theory with and without the constraint $z=0$ can also be found in [4, 5].

The use of Green’s theorem for planar electrostatic problems in the form proposed in [13] makes calculating the electrostatic potential of a uniformly charged disk in terms of complete elliptic integrals an almost trivial task.

The reason for this simplification is that the method fully exploits the cylindrical symmetry of the problem. If this symmetry is broken, such as in the case of an elliptical disk, the calculations will not be as simple.

The electrostatic potential of a uniformly charged elliptical disk was found in [4, 5]. Of course, it will be interesting to reproduce this result using Green’s theorem in the form (2).

A drawback of the method under consideration is that its application is limited only to the case of uniform charge density.

Nevertheless, we hope that students will benefit from this approach and that this interesting problem will finally find its place in introductory classical electrodynamics textbooks.

References