Table Of ContentCharged line segments and ellipsoidal equipotentials
T L Curtright§, N M Aden, X Chen, M J Haddad, S Karayev, D B Khadka, and J Li
Department of Physics, University of Miami
Coral Gables, FL 33124-8046, USA
§
6 [email protected]
1
0 January 18, 2016
2
n
a Abstract
J
Thisisasurveyoftheelectrostaticpotentialsproducedbychargedstraight-linesegments,invarious
5
numbers of spatial dimensions, with comparisons between uniformly charged segments and those hav-
1
ing non-uniform linear charge distributions that give rise to ellipsoidal equipotentials surrounding the
] segments. A uniform linear distribution of charge is compatible with ellipsoidal equipotentials only for
h
threedimensions. Inhigherdimensions,thelinearchargedensitygivingrisetoellipsoidalequipotentials
p
iscounter-intuitive—thechargedistributionhasamaximumatthecenterofthesegmentandvanishes
-
s attheendsofthesegment. Onlyintwodimensionsisthecontinuouschargedistributionintuitive—for
s thatonecaseofellipsoidalequipotentials,thechargeispeakedattheendsofthesegmentandminimized
a
at the center.
l
c
.
s
c Contents
i
s
y 1 Introduction 1
h
p 2 Electrostatics in two dimensions 2
[
1 3 Point charge in D dimensions 13
v
7 4 Uniformly charged line segments for all dimensions D >2 14
4
0 5 Ellipsoidal equipotentials for all dimensions D >2 19
4
0
6 Generalizations 25
.
1
0 7 Summary 26
6
1
:
v
i
X
r
a
1 Introduction
Ithasbecomeawidespreadpracticetostudythephysicsofsystemsinvariousnumbersofspatialdimensions,
not necessarily D =3. For instance, graphene with D =2, and string or membrane theory with D as high
as25, aretwoexamplesthatimmediatelycometomindforboththeirexperimentalandtheoreticalinterest.
Moreover, geometric ideas provide a common framework used to pursue such studies. A pedagogical goal
of this paper is to encourage students to think along these lines in the context of a familiar subject —
electrostatics.
1
Forexample,inthreespatialdimensionsauniformlychargedstraight-linesegmentgivesrisetoanelectric
potentialΦwhoseequipotentialsurfacesareprolateellipsoidsofrevolutionaboutthesegment,withtheends
of the segment providing the foci of the ellipsoid. As an immediate consequence of the geometry for these
→− →−
prolate ellipsoidal equipotentials, the associated electric field E = −∇Φ — always normal to surfaces of
constant Φ — has at any observation point a direction that bisects the angle formed by the pair of lines
from the observation point to each of the two foci of the ellipsoid. This beautiful electrostatic example was
presentedbyGeorgeGreenin1828[1], andithasbeendiscussedinmanybookssincethen[2]-[10]including
at least two texts from this century [11]. While the straight-line segment is an idealization, nevertheless
it provides insight into the behavior of real thin-wire conductors, especially upon approximating those real
wires as very narrow, “needle-like” ellipsoids.
Butasitturnsout,thelinesegmentprobleminthreedimensionsisaveryspecialcase,insomesensethe
most ideal of all possible electrostatic worlds. In any other dimension of space, uniformly charged segments
do not produce ellipsoidal equipotentials. Conversely, in any other dimension, if the equipotentials are
ellipsoidalaboutalinearlydistributedstraight-linesegmentofcharge,thenthatchargedistributioncannot
be uniform.
One need look no farther than two dimensional systems to see clearly that there are differences between
uniformly charged segments and those with charge distributed so as to produce ellipsoidal equipotentials.
Indeed, D = 2 is the only ellipsoidal equipotential case which is intuitive in the sense that the associated
linear charge distribution has maxima at the ends of the segment, as one might naively expect from the
repulsive force between like charges placed on a segment of a real, thin conductor at a finite potential. In
contrast,thedistributionofchargeneededtoproduceellipsoidalequipotentialsiscounter-intuitiveinhigher
dimensions. ToproducesuchequipotentialsforD >3thechargedistributionmusthaveamaximumatthe
center and vanish at the ends of the segment.
We begin our discussion in Section 2 in two dimensions, where the two types of charged segments are
readily analyzed. Then we compare and contrast uniformly charged line segments with those admitting
ellipsoidal equipotentials for any number of spatial dimensions. As a preliminary, we first discuss briefly in
Section 3 the potential of a point charge in D dimensions. We then use this information in Section 4 to
compute the potentials and electric fields for uniformly charged line segments. We continue in Section 5 by
considering systems with ellipsoidal equipotentials in D dimensions. We then determine the linear charge
distributions that produce such potential configurations, and we find the remarkable result that the linear
charge density giving rise to ellipsoidal equipotentials is peaked at the center of the segment, for D >3.
The discussion of general D affords the opportunity to illustrate how continuous D can be used as a
mathematical device to regulate singular behavior. This too is a widespread practice in theoretical physics.
We use D in this way in Section 5 of the paper to interpolate continuously between intuitive and counter-
intuitive charge distributions for segments with ellipsoidal equipotentials.
Finally, inSection6, weinvokewell-knownmethodstoindicatehowthevariouslinesegmentresultsalso
give solutions to a class of electrostatic boundary value problems where the charge is moved outward and
distributed on one of the equipotential surfaces that surrounds the original segment.
2 Electrostatics in two dimensions
The point-particle electric potential in 2D is well-known to be logarithmic. For a point charge Q located at
the origin, up to a constant R that sets the distance scale,
Φ (→−r)=kQln(R/r) , ∇2Φ (→−r)=−2πkQ δ2(→−r) , (1)
point point
where k is the two dimensional analogue of Coulomb’s constant.
By linear superposition, a finite length, uniformly charged line segment, with constant charge/length λ
distributed on the x axis for −L/2≤x≤L/2, produces a potential
(cid:18) (cid:112) u(cid:19)(cid:12)(cid:12)u=−x+L/2
Φline(x,y)=−kλ uln u2+y2−u+yarctan y (cid:12)(cid:12) +kλLln(R) . (2)
u=−x−L/2
2
This result may be established by integrating the contributions of infinitesimal point-like bits of charge that
make up the segment, using (1) and the indefinite integral
(cid:90) (cid:16)(cid:112) (cid:17) (cid:16)(cid:112) (cid:17) u
ln u2+y2 du=uln u2+y2 −u+yarctan . (3)
y
Written out in full,
(cid:18) (cid:18)x− 1L(cid:19) (cid:18)x+ 1L(cid:19)(cid:19)
Φ (x,y)=kλ L+yarctan 2 −yarctan 2 (4)
line y y
(cid:115) (cid:115)
(cid:18) 1 (cid:19) 1 (cid:18) 1 (cid:19)2 (cid:18) 1 (cid:19) 1 (cid:18) 1 (cid:19)2
+ x− L kλln x− L +y2− x+ L kλln x+ L +y2 .
2 R 2 2 R 2
A plot of the potential surface shows the essential features.
x y
Figure 1: Potential surface for a uniformly charged line segment in 2D.
The top of the potential surface is curved and not a straight line, indicating that the charged segment
itself is not an equipotential. This follows analytically from (4). Although points on the segment are not
at the same potential, they are all at finite values of the potential, for D =2. Explicitly, for y =0 and all
x,
Φline(x,0)=kλL−kλ(cid:0)21L+x(cid:1)ln(cid:0)(cid:12)(cid:12)12L+x(cid:12)(cid:12)/R(cid:1)−kλ(cid:0)12L−x(cid:1)ln(cid:0)(cid:12)(cid:12)12L−x(cid:12)(cid:12)/R(cid:1) . (5)
We plot f(x)=Φ (x,0)/kλ versus x to show the shape of the potential along the x-axis, for L=2 and
line
R=1.
3
f
2
4 2 2 4
x
2
Figure 2: Uniformly charged line segment potential for D =2, along the x-axis.
Ontheotherhand,transversetothex-axisthepotentialhasadiscontinuousslopeacrossthelinesegment.
That is to say, the electric field normal to the line segment, E , is discontinuous due to the presence of the
y
charge density on the segment. For example, for x=0 this transverse profile is given by
(cid:20) (cid:18)(cid:113) (cid:19) (cid:21)
Φ(0,y)=kλ L−Lln 1L2+y2/R −2yarctan(cid:0)1L/y(cid:1) . (6)
4 2
We also plot g(y) = Φ (0,y)/kλ versus y to show the shape of the potential along the y-axis, for L = 2
line
and R=1.
g
2
4 2 2 4
y
2
Figure 3: Uniformly charged line segment potential for D =2, along the y-axis.
Moreover, the equipotentials are not ellipses surrounding the segment on the xy-plane, as is especially
clear for points close to the segment. The potential contours actually intersect the segment. A view of the
potentialcontoursfrombelow thepotentialsurfaceshowsthesefeaturesgraphically. (Seetheplottofollow.
But note the view in that plot is actually an orthogonal projection of the contours onto the xy-plane, and
not the true perspective of an observer on the potential axis a finite distance below that plane.)
4
x
2 1 0 1 2
y
Figure 4: Potential contours for a uniformly charged line segment in 2D.
The components of the electric field produced by the segment are given by
(cid:118)
∂ (cid:117)(cid:117)(cid:0)x+ 1L(cid:1)2+y2
Ex(x,y)=−∂xΦline(x,y)=kλln(cid:116)(cid:0)x− 21L(cid:1)2+y2 , (7)
2
(cid:18) (cid:19) (cid:18) (cid:19)
∂ L+2x L−2x
E (x,y)=− Φ (x,y)=kλarctan +kλarctan , (8)
y ∂y line 2y 2y
→−
where we have used d arctanu = 1 . Note that E (x,y) is independent of the scale R, since changing
du u2+1
→−
R just amounts to adding a constant to the potential. But of course E (x,y) does depend on L as this sets
the physical length scale for the system.
(cid:82) →− →−
As a check, in the units we have chosen the charge density along the line follows from ∇ · E dxdy =
(cid:82)
2πk ρ dxdy, integrated over a horizontal rectangle containing an infinitesimal portion of the x-axis, as the
height of the rectangle is taken to zero. Thus, using lim arctanz =±π/2,
z→±∞
0 if x>L/2
2πk λ(x)=2limE (x,y)=2kλ× π if −L/2<x<L/2 . (9)
y
y→0 0 if x<−L/2
That is to say, the charge/length between −L/2 and L/2 is the constant λ, as expected. Elsewhere, the
→−
charge density vanishes, as follows from ∇2Φ (r)=0 for all points not coincident with the segment.
line
It is instructive to make a vector field plot of (E (x,y),E (x,y)), especially near the charged segment.
x y
Again it is evident graphically that the segment itself is not an equipotential, since the electric field lines
are not perpendicular to the segment as the x-axis is approached for −L/2<x<L/2, except at the single
point x=0. Here are vector field plots for L=2.
5
2
y
1
2 1 1 2
x
1
2
y 0.1
1.2 1.0 0.8 0.6 0.4 0.2 0.2 0.4 0.6 0.8 1.0 1.2
x
0.1
→−
Figures 5 & 6: Vector field plots of E for a uniformly charged line segment situated between x=−1 and
x=+1. The field is evaluated at the center of each arrow.
6
In contrast to the uniformly charged segment in 2D, consider a distribution of charge along the segment
such that equipotentials are ellipsoidal. The relevant charge distribution turns out to be
0 if x>L/2
2λL 1
λ(x)= √ 1 if −L/2<x<L/2 , (10)
π L2−4x2 0 if x<−L/2
as we shall confirm in the following. Since (cid:82)1 √ 1 ds = π, the total charge on the segment is still
−1 1−s2
(cid:82)L/2
Q= λ(x) dx=λL, the same as for the uniformly charged case.
−L/2
The corresponding potential is now
(cid:18) (cid:19)
4R
Φ (x,y)=kλLln √ , (11)
line
s+ s2−L2
where s is a sum of two distances, from the observation point (x,y) to each of the two ends of the segment.
That is,
(cid:115)
(cid:18) 1 (cid:19)2
s=r +r , r = x± L +y2 . (12)
− + ± 2
Because the positional dependence of the potential is given entirely by s, the equipotentials are ellipses on
the xy-plane, with the ends of the segment serving as the foci of each equipotential ellipse. Again, we plot
the potential surface.
x y
Figure 7: Potential surface for a non-uniformly charged line segment in 2D, with λ(x)= √2λL/π .
L2−4x2
The top of the potential surface is now a straight line, indicating that the charged segment itself is an
equipotential,andthepointsonthesegmentareatafinite potentialforD =2,namely,Φ=kλLln(4R/L).
7
This follows analytically from (11). Explicitly, for y =0 and all x,
4R
Φline(x,0)=kλLln(cid:12)(cid:12)x− 1 L(cid:12)(cid:12)+(cid:12)(cid:12)x+ 1 L(cid:12)(cid:12)+(cid:113)2(cid:0)x2− 1L2+(cid:12)(cid:12)x2− 1 L2(cid:12)(cid:12)(cid:1) . (13)
2 2 4 4
We plot f(x) = 1 Φ (x,0) versus x to show the shape of the potential along the x-axis, for L = 2 and
kλ line
R=1.
f
1
4 2 2 4
x
1
2
3
Figure 8: Non-uniformly charged line segment potential for D =2, along the x-axis.
Transverse to the x-axis, the potential again has a discontinuous slope across the line segment due to the
presence of the charge density on the segment. For example, for x=0 this transverse profile is given by
(cid:32) (cid:33)
4R
Φ (0,y)=Lkλln . (14)
line (cid:112)
L2+4y2+2|y|
We plot g(y) = 1 Φ (0,y) versus y to show the shape of the potential along the y-axis, for L = 2 and
kλ line
R=1.
g
1
4 2 2 4
y
1
2
3
Figure 9: Non-uniformly charged line segment potential for D =2, along the y-axis.
8
A view of the potential contours from below the potential surface shows the equipotential ellipses. (But
again note the view in that plot is actually an orthogonal projection of the contours onto the xy-plane, and
not the true perspective of an observer on the potential axis a finite distance below that plane.)
x
2 1 0 1 2
y
Figure 10: Ellipsoidal potential contours for a non-uniformly charged line segment in 2D.
The components of the electric field are now
∂ kλL ∂s
E (x,y)=− Φ (x,y)= √ , (15)
x ∂x line s2−L2∂x
∂ kλL ∂s
E (x,y)=− Φ (x,y)= √ , (16)
y ∂y line s2−L2∂y
whose final form in rectangular coordinates follows from
∂s x− 1 L x+ 1 L 1 (cid:18) 1 (cid:19) ∂s y y 1
= 2 + 2 = sx− (r −r )L , = + = sy . (17)
∂x r r r r 2 + − ∂y r r r r
− + + − − + + −
In terms of x, y, r , and r ,
+ −
(cid:18)(cid:18) (cid:19) (cid:18) (cid:19) (cid:19)
kλL 1 1
E (x,y)= x− L r + x+ L r , (18)
x r r (cid:113)(r +r )2−L2 2 + 2 −
+ − − +
kλL
E (x,y)= (r +r )y . (19)
y (cid:113) − +
r r (r +r )2−L2
+ − − +
9
Writing everything out in terms of x and y gives
(cid:115) (cid:115)
(cid:18)1 (cid:19)2 1 (cid:18)1 (cid:19)2
r+r− = 4 L2+x2+y2 −L2x2, (r−+r+)2−L2 =2x2+y2− 4L2+ 4 L2+x2+y2 −L2x2 .
(20)
→−
The components of E are not elegant in these coordinates [13] but their properties are fully encoded and
amenable to machine computation.
(cid:113) (cid:113)
kλL (cid:0)x+ 1 L(cid:1) (cid:0)x− 1 L(cid:1)2+y2+(cid:0)x− 1 L(cid:1) (cid:0)x+ 1 L(cid:1)2+y2
E (x,y)= √ 2 2 2 2 , (21)
x (cid:114)
2 (cid:113) (cid:113)
(cid:0)1 L2+x2+y2(cid:1)2−L2x2 x2+y2− 1L2+ (cid:0)1 L2+x2+y2(cid:1)2−L2x2
4 4 4
(cid:18)(cid:113) (cid:113) (cid:19)
(cid:0)x− 1 L(cid:1)2+y2+ (cid:0)x+ 1 L(cid:1)2+y2 y
kλL 2 2
E (x,y)= √ . (22)
y (cid:114)
2 (cid:113) (cid:113)
(cid:0)1 L2+x2+y2(cid:1)2−L2x2 x2+y2− 1L2+ (cid:0)1 L2+x2+y2(cid:1)2−L2x2
4 4 4
Although the geometric features of the electric field may not be transparent from these expressions, never-
→−
theless the fact that the equipotentials are ellipsoidal allows one to immediately visualize the direction of E
as normal to those surfaces of constant Φ.
(cid:82) →− →− (cid:82)
To confirm the charge density along the line, we again use ∇ · E dxdy =2πk ρ dxdy, and integrate
over a horizontal rectangle containing an infinitesimal portion of the x-axis, as the height of the rectangle is
taken to zero. The crucial features here are s−L → |2x−L| > 0 as y → 0 for |x| > L/2 (away from the
segment), but s−L→0 as y →0 for −L/2<x<L/2 (on the segment). More precisely, as points on the
segment are approached transversely,
(cid:18)1 (cid:19)(cid:115) y2 (cid:18)1 (cid:19)(cid:115) y2
s= L−x 1+ + L+x 1+
2 (cid:0)1 L−x(cid:1)2 2 (cid:0)1 L+x(cid:1)2
2 2
(cid:18) (cid:19) (cid:18) (cid:19)
∼ L+ 1 1 + 1 y2+O(cid:0)y4(cid:1)=L+ 1 L 1 y2+O(cid:0)y4(cid:1) . (23)
y→0 2 1 L−x 1 L+x 2 1 L2−x2
2 2 4
Thus for −L/2<x<L/2,
(cid:18) (cid:19) (cid:114)
y 1 y 1 y 1 1
lim √ = √ lim √ = √ lim = L2−x2 . (24)
(cid:114)
y→0 s2−L2 2Ly→0 s−L 2Ly→0 (cid:16) (cid:17) L 4
1 L 1 y2
2 1 L2−x2
4
Combining this with the more obvious lim (r r )= 1 L2−x2 for −L/2<x<L/2 gives
y→0 + − 4
(cid:18)2kλLs y (cid:19) 2kλL 0 if x>L/2
2πk λ(x)=2limE (x,y)= lim √ = × 1 if −L/2<x<L/2 .
y→0 y y→0 r+r− s2−L2 (cid:113)1 L2−x2 0 if x<−L/2
4
(25)
That is to say, the charge/length between −L/2 and L/2 is
2λL
λ(x)= √ , (26)
π L2−4x2
→−
as anticipated above in (10). Elsewhere, the charge density vanishes, again as follows from ∇2Φ (r)=0
line
for all points not coincident with the segment.
Oncemoreitisinstructivetomakeavectorfieldplotof(E (x,y),E (x,y)),especiallynearthecharged
x y
segment. The fact that the segment itself is an equipotential is evident graphically since the electric field
lines are perpendicular to the segment as the x-axis is approached for −L/2 < x < L/2. Here are vector
field plots for L=2.
10