Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic...

51
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, ”Classical Electrodynamics”, 2nd Edition, Section 1-5 Classical Field Theory: Electrostatics-Magnetostatics

Transcript of Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic...

Page 1: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Classical Field Theory:Electrostatics-Magnetostatics

April 27, 20101

1J.D.Jackson, ”Classical Electrodynamics”, 2nd Edition, Section 1-5Classical Field Theory: Electrostatics-Magnetostatics

Page 2: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Electrostatics

The behavior of an electrostatic field can be described by two differentialequations:

~∇ · ~E = 4πρ (1)

(Gauss law) and~∇× ~E = 0 (2)

the latter equation being equivalent to the statement that ~E is thegradient of a scalar function, the scalar potential Φ:

~E = −~∇Φ (3)

Eqns (1) and (3) can be combined into one differential equation for asingle scalar function Φ(~x):

∇2Φ = −4πρ (4)

This equation is called Poisson equation.In the regions of space where there is no charge density, the scalarpotential satisfies the Laplace equation:

∇2Φ = 0 (5)

Classical Field Theory: Electrostatics-Magnetostatics

Page 3: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

For a general distribution ρ(~x ′), the potential is expected to be the sumover all increments of charge d3x ′ρ(~x ′), i.e.,

Φ(~x) =

∫ρ(~x ′)

|~x − ~x ′|d3x ′ (6)

This potential should satisfy Poisson’s equation. But does it? If weoperate with ∇2 on both sides of (6) we get (on ~x not on ~x ′)

∇2Φ(~x) =

∫ρ(~x ′)d3x ′∇2

(1

|~x − ~x ′|

)(7)

But ∇2(1/|~x − ~x ′|) = 0 as long as ~x 6= ~x ′ ! (Why ?)The singular nature of ∇2(1/|~x − ~x ′|) = ∇2(1/r) can be best expressedin terms of the Dirac δ-function.Since ∇2(1/r) = 0 for r 6= 0 and its volume integral is −4π (Why?) wecan write

∇2

(1

|~x − ~x ′|

)= −4πδ(|~x − ~x ′|) (8)

Classical Field Theory: Electrostatics-Magnetostatics

Page 4: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

By definition, if the integration volume contains the point ~x = ~x ′∫δ3(~x − ~x ′)d3x = 1

otherwise is zero. This way we recover Poisson’ s equation

∇2Φ(~x) = −4πρ(~x ′)~x′=~x (9)

Thus, we have not only shown that the potential from Coulomb’s lawsatisfies Poisson’s eqn, but we have established (through the solution ofPoisson’s eqn) the important result that :the potential from a distributed source is the superposition of theindividual potentials from the constituent parcels of charge.We may consider situations in which ρ is comprised of N discrete chargesqi , positioned at ~x ′i so that

ρ(~x ′) =N∑

i=1

qiδ3 (~x ′ − ~x ′i ) (10)

In this case the solution for the potential is a combination of terms

proportional to 1/|~x − ~x ′|.Classical Field Theory: Electrostatics-Magnetostatics

Page 5: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Green Theorem

If in the electrostatic problem involved localized discrete or continuousdistributions of charge with no boundary surfaces, the general solution(6) would be the most convenient and straight forward solution to anyproblem.To handle the boundary conditions it is necessary to develop some newmathematical tools, namely, the identities or theorems due to GeorgeGreen (1824). These follow as simple applications of the divergencetheorem ∫

V

~∇ · ~A d3x =

∮S

~A · ~n da (11)

which applies to any well-behaved vector field ~A defined in the volume Vbounded by the closed surface S .Let ~A = Φ~∇ψ, (Φ and ψ arbitrary scalar fields). Then

~∇ · (Φ~∇ψ) = Φ∇2ψ + ~∇Φ · ~∇ψ (12)

Φ~∇ψ · ~n = Φ∂ψ

∂n(13)

where ∂/∂n is the normal derivative at the surface S .

Classical Field Theory: Electrostatics-Magnetostatics

Page 6: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

When (12) and (13) substituted into the divergence theorem (8)produces the so-called Green’s 1st identity∫

V

(Φ∇2ψ + ~∇Φ · ~∇ψ

)d3x =

∮S

Φ∂ψ

∂nda . (14)

If we rewrite (14) with Φ and ψ interchanged, and then subtract it from(14) we obtain Green’s 2nd identity or Green’s Theorem:∫

V

(Φ∇2ψ − ψ∇2Φ

)d3x =

∮S

[Φ∂ψ

∂n− ψ∂Φ

∂n

]da (15)

Now we can apply Poisson’s equation (8) for discrete charge, substitutingfor ψ = 1/|~x − ~x ′|∫

V

[−4πδ3(~x − ~x ′)Φ(~x ′) +

4πρ(~x ′)

|~x − ~x ′|

]d3x

=

∮S

[Φ∂

∂n′

(1

|~x − ~x ′|

)− 1

|~x − ~x ′|∂Φ

∂n′

]da′ (16)

Classical Field Theory: Electrostatics-Magnetostatics

Page 7: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Integrating the Dirac delta function over all values of ~x ′ within V and for~x within the volume V yields a nonzero result

Φ(~x) =

∫V

ρ(~x ′)

|~x − ~x ′|d3x ′ +

1

∮S

[1

|~x − ~x ′|∂Φ

∂n′− Φ

∂n′

(1

|~x − ~x ′|

)]da′

(17)• The (blue) correction term goes to zero as the surface S goes toinfinity (because S falls of faster than 1/|~x − ~x ′|)• If the integration volume is free of charges, then the first term ofequation (17) becomes zero, and the potential is determined only by thevalues of the potential and the values of its derivatives at the boundaryof the integration region (the surface S).

Classical Field Theory: Electrostatics-Magnetostatics

Page 8: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

• Physical experience leads us to believe that specification of thepotential on a closed surface defines a unique potential problem. This iscalled Dirichlet problem or Dirichlet boundary conditions.• Similarly it is plausible that specification of the electric field (normalderivative of the potential) everywhere on the surface (corresponding to agiven surface-charge density) also defines a unique problem. Thespecification of the normal derivative is known as the Newmannboundary condition.• As it turns out either one of the two conditions results in uniquesolution. It should be clear that a solution to the Poisson eqn with bothΦ and ∂Φ/∂n specified arbitrarily on a closed boundary (Cauchyboundary conditions) does not exist since there are unique solutions forDirichlet and Newmann condition separately.

Classical Field Theory: Electrostatics-Magnetostatics

Page 9: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Green Functions

The solution of the Poisson or Laplace eqn in a finite volume V witheither Dirichlet or Neumann boundary conditions on the boundary surfaceS can be obtained by means of Green’s theorem (15) and the so-calledGreen functions.In obtaining the result (17) we have chosen ψ = 1/|~x − ~x ′| satisfying

∇2

(1

|~x − ~x ′|

)= −4πδ(|~x − ~x ′|) (18)

The function 1/|~x − ~x ′| is only one of a class of functions depending onthe variables ~x and ~x ′ and called Green functions.In general

∇′2G (~x ,~x ′) = −4πδ(|~x − ~x ′|) (19)

where

G (~x ,~x ′) =1

|~x − ~x ′|+ F (~x ,~x ′) (20)

and F satisfying the Laplace equation inside the volume V

∇′2F (~x ,~x ′) = 0 (21)

Classical Field Theory: Electrostatics-Magnetostatics

Page 10: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

If we substitute G (~x ,~x ′) in eqn (17) we get

Φ(~x) =

∫V

ρ(~x ′)G (~x ,~x ′)d3x ′+1

∮S

[G (~x ,~x ′)

∂Φ

∂n′− Φ(~x ′)

∂n′G (~x ,~x ′)

]da′

(22)The freedom in the definition of G means that we can make the surfaceintegral depend only on the chosen type of BC.• For the Dirichlet BC we demand:

GD(~x ,~x ′) = 0 for ~x ′ ∈ S (23)

Then the 1st term on the surface integral of (22) vanishes

Φ(~x) =

∫V

ρ(~x ′)GD(~x ,~x ′)d3x ′ − 1

∮S

Φ(~x ′)∂

∂n′GD(~x ,~x ′)da′ (24)

• For Neumann BC the simplest choice of BC on G is

∂GN

∂n′(~x ,~x ′) = 0 for ~x ′ ∈ S (25)

but application o the Gauss theorem on (19) shows that (how?)∮S

∂GN

∂n′da′ = −4π 6= 0

which is incosistent with ∇′2G (~x ,~x ′) = −4πδ3(~x − ~x ′).

Classical Field Theory: Electrostatics-Magnetostatics

Page 11: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

This will mean that the outflux of G cannot be zero when there is asource enclosed by S . Then the simplest boundary condition that we canuse is

∂GN

∂n′(~x ,~x ′) = −4π

Sfor ~x ′ ∈ S (26)

S is the total area of the boundary surface. Then the solution will be:

Φ(~x) =

∫V

ρ(~x ′)GN(~x ,~x ′)d3x ′ + 〈Φ〉S +1

∮S

∂Φ

∂n′GN(~x ,~x ′)da′ (27)

where 〈Φ〉S is the average value of the potential over the whole surface

〈Φ〉S ≡1

S

∮S

Φ(~x ′)da′ (28)

In most cases S is extremely large (or even infinite), in which case〈Φ〉S → 0.

The physical meaning of F (~x ,~x ′) : it is a solution of the Laplace eqn

inside V and so represents the potential of charges external to the

volume V chosen as to satisfy the homogeneous BC of zero potential on

the surface S .

Classical Field Theory: Electrostatics-Magnetostatics

Page 12: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Green Function

• It is important to understand that no matter how the source isdistributed, finding the Green function is completely independent or ρ(~x ′).• G (~x ,~x ′) depends exclusively on the geometry of the problem, is a“template”, potential and not the actual potential for a given physicalproblem.• In other words G (~x ,~x ′) is the potential due to a unit charge, positionedarbitrarily within the surface S consistent either with GD = 0 or∂GN/∂n = −4π/S on the surface.

• The “true potential” is a convolution of this template with the given

ρ(~x ′).

Classical Field Theory: Electrostatics-Magnetostatics

Page 13: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Laplace Equation in Spherical Coordinates

In spherical coordinates (r , θ, φ) the Laplace equation can be written inthe form

1

r2

∂2

∂r2(rΦ) +

1

r2 sin θ

∂θ

(sin θ

∂Φ

∂θ

)+

1

r2 sin2 θ

∂2Φ

∂φ2= 0 (29)

If we assume

Φ =1

rU(r)P(θ)Q(φ) (30)

Then by substituting in (29) and multiplying with r2 sin2 θ/UPQ weobtain

r2 sin2 θ

[1

U

d2U

dr2+

1

r2 sin θP

d

(sin θ

dP

)]+

1

Q

d2Q

dφ2= 0 (31)

We see that the therms depending on φ have been isolated and we canset

1

Q

d2Q

dφ2= −m2 (32)

with solutionQ = e±imφ (33)

Classical Field Theory: Electrostatics-Magnetostatics

Page 14: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Similarly the remaining terms can be separated as:

1

r2 sin θ

d

(sin θ

dP

)+

[l(l + 1)− m2

sin2 θ

]P = 0 (34)

d2U

dr2− l(l + 1)

r2U = 0 (35)

The radial equation will have a solution

U = Ar l+1 + Br−l (36)

while l is still undetermined.

Classical Field Theory: Electrostatics-Magnetostatics

Page 15: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Legendre Equation and Legendre Polynomials

The equation (34) for P(θ) can be expressed in terms of x = cos θ

d

dx

[(1− x2)

dP

dx

]+

[l(l + 1)− m2

1− x2

]P = 0 (37)

This is the generalized Legendre equation and its solutions are theassociated Legendre functions.We will consider the solution of (54) for m2 = 0

d

dx

[(1− x2)

dP

dx

]+ l(l + 1)P = 0 (38)

The solution should be single valued, finite, and continuous on theinterval −1 ≤ x ≤ 1 in order that it represents a physical potential.The solution can be found in the form of a power series

P(x) = xk∞∑j=0

ajxj (39)

where k is a parameter to be determined.

Classical Field Theory: Electrostatics-Magnetostatics

Page 16: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

By substitution in (54) we get the recurrence relation (how?)

aj+2 =(k + j)(k + j + 1)− l(l + 1)

(k + j + 1)(k + j + 2)aj (40)

while for j = 0, 1 we find that

if a0 6= 0 then k(k − 1) = 0 (41)

if a1 6= 0 then k(k + 1) = 0 (42)

These two relations are equivalent and it is sufficient to choose either a0

or a1 different from zero but not both (why?). We also see that the seriesexpansion is either only odd or only on even powers of x . By choosingeither k = 0 or k = 1 it is possible to prove the following properties:

I The series converges for x2 < 1, independent of the value of l

I The series diverges for x = ±1, unless it terminates.

Since we want solution that is finite at x = ±1, as well as forx2 < 1, we demand that the series terminates.

Since k and j are positive integers or zero, the recurrence relation (40)

will terminate only if l is zero or a positive integer.

Classical Field Theory: Electrostatics-Magnetostatics

Page 17: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

If l is even (odd), then only the k = 0 (k = 1) series terminates.The polynomials in each case have x l as their highest power. Byconvention these polynomials are normalized to have the value unity forx = +1 and are called the Legendre polynomials of order l .The first few are:

P0(x) = 1

P1(x) = x

P2(x) =1

2(3x2 − 1) (43)

P3(x) =1

2(5x3 − 3x)

P4(x) =1

8(35x4 − 30x2 + 3)

The Legendre polynomials can be taken from Rodrigue’s formula:

Pl(x) =1

2l l!

d l

dx l(x2 − 1)l (44)

Classical Field Theory: Electrostatics-Magnetostatics

Page 18: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

The Legendre polynomials form a complete orthonormal set of functionson the interval −1 ≤ x ≤ 1 (prove it)∫ 1

−1

Pl′(x)Pl(x) =2

2l + 1δl′l (45)

Since the Legendre polynomials form a complete set of orthonormalfunctions, any function f (x) on the interval −1 ≤ x ≤ 1 can be expandedin terms of them i.e.

f (x) =∞∑l=0

AlPl(x) (46)

where (how?)

Al =2l + 1

2

∫ 1

−1

f (x)Pl(x)dx (47)

Thus for problems with azimuthal symmetry i.e. m = 0 the generalsolution is:

Φ(r , θ) =∞∑l=0

[Al r

l + Bl r−(l+1)

]Pl(cos θ) (48)

where the coefficients Al [it is not the same as in eqn (47)] and Bl can

be determined from the boundary conditions.

Classical Field Theory: Electrostatics-Magnetostatics

Page 19: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Example

Boundary Value Problems with Azimuthal Symmetry

Let’s specify as V (θ) the potential on the surface of a sphere of radius R,and try to find the potential inside the sphere.If there are no charges at the origin (r = 0) the potential must be finitethere. Consequently Bl = 0 for all l . Then on the surface of the sphere

V (r = R, θ) =∞∑l=0

AlRlPl(cos θ) (49)

and the coefficients Al will be taken via eqn(47)

Al =2l + 1

2R l

∫ π

0

V (θ)Pl(cos θ) sin θdθ (50)

Classical Field Theory: Electrostatics-Magnetostatics

Page 20: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

If, for example V (θ) = ±V on the two hemispheres then the coefficientscan be derived easily and the potential inside the sphere is (how?):

Φ(r , θ) = V

[3

2

( r

R

)P1(cos θ)− 7

8

( r

R

)3

P3(cos θ) +11

16

( r

R

)5

P5(cos θ)− . . .]

(51)To find the potential outside the sphere we merely replace (r/R)l by(R/r)l+1 and the resulting potential will be (how?):

Φ(r , θ) =3

2

(R

r

)2

V

[P1(cos θ)− 7

12

(R

r

)2

P3(cos θ) + . . .

](52)

Classical Field Theory: Electrostatics-Magnetostatics

Page 21: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Legendre Polynomials

An important expansion is that of thepotential at ~x due to a unit point chargeat ~x ′

1

|~x − ~x ′|=∞∑l=0

r l<

r l+1>

Pl(cos γ) (53)

where r<(r>) is the smaller (larger) of |~x |and |~x ′| and γ is the angle between |~x |and |~x ′|.

Show that the potential is :

1

|~x − ~x ′|=∞∑l=0

(Al r

l + Bl r−(l+1)

)Pl(cos γ) on the z-axis

1

|~x − ~x ′|=

1

r>

∞∑l=0

(r<r>

)l

on the x-axis

Classical Field Theory: Electrostatics-Magnetostatics

Page 22: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Associated Legendre Functions and Spherical Harmonics

For problems without azimuthal (axial) symmetry, we need thegeneralization of Pl(cos θ), namely, the solutions of

d

dx

[(1− x2)

dP

dx

]+

[l(l + 1)− m2

1− x2

]P = 0 (54)

for arbitrary l and m.It can be shown that in order to have finite solutions on the interval−1 ≤ x ≤ 1 the parameter l must be zero or a positive integer and thatthe integer m can take on only the values −l , −(l − 1), ..., 0 , ...,(l − 1), l (why?).The solution having these properties is called an associated Legendrefunction Pm

l (x) . For positive m it is defined as

Pml (x) = (−1)m(1− x2)m/2 dm

dxmPl(x) (55)

If Rodrigues’ formula is used an expression valid for both positive andnegative m is obtained:

Pml (x) =

(−1)m

2l l!(1− x2)m/2 dm

dxmPl(x) (56)

Classical Field Theory: Electrostatics-Magnetostatics

Page 23: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

There is a simple relation between P−ml (x) and Pm

l (x) :

P−ml (x) = (−1)m (l −m)!

(l + m)!Pm

l (x) (57)

For fixed m the functions Pml (x) form an orthonormal set in the index l

on the interval −1 ≤ x ≤ 1. The orthogonality relation is∫ 1

−1

Pml′ (x)Pm

l (x)dx =2

2l + 1

(l + m)!

(l −m)!δl′l (58)

We have found that Qm(φ) = e imφ, this function forms a complete set oforthogonal functions in the index m on the interval 0 ≤ φ ≤ 2π.The product Pm

l Qm forms also a complete orthonormal set on the surfaceof the unit sphere in the two indices l ,m.From the normalization condition (58) we can conclude that the suitablynormalized functions, denoted by Ylm(θ, φ), are :

Ylm(θ, φ) =

√2l + 1

(l −m)!

(l + m)!Pm

l (cos θ)e imφ (59)

and alsoYl,−m(θ, φ) = (−1)mY ∗lm(θ, φ) (60)

Classical Field Theory: Electrostatics-Magnetostatics

Page 24: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

The normalization and orthogonality conditions are:∫ 2π

0

∫ π

0

sin θdθY ∗l′m′(θ, φ)Ylm(θ, φ) = δl′lδm′m (61)

An arbitrary function g(θ, φ) can be expanded in spherical harmonics

g(θ, φ) =∞∑l=0

l∑m=−l

AlmYlm(θ, φ) (62)

where the coefficients are

Alm =

∫dΩY ∗lm(θ, φ)g(θ, φ) . (63)

The general solution for a boundary-value problem in sphericalcoordinates can be written in terms of spherical harmonics and powers ofr in a generalization of (48) :

Φ(r , θ, φ) =∞∑l=0

l∑m=−l

[Almr l + Blmr−(l+1)

]Ylm(θ, φ) (64)

If the potential is specified on a spherical surface, the coefficients can bedetermined by evaluating (64) on the surface and using (63).

Classical Field Theory: Electrostatics-Magnetostatics

Page 25: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Spherical Harmonics Ylm(θ, φ)

l = 0 Y00 =

√1

l = 1 Y11 = −√

3

8πsin θe iφ

Y10 =

√3

4πcos θ

l = 2 Y22 =1

4

√15

2πsin2 θe2iφ

Y21 = −√

15

8πsin θ cos θ e iφ

Y20 =

√5

(3

2cos2 θ − 1

2

)

Classical Field Theory: Electrostatics-Magnetostatics

Page 26: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Spherical Harmonics Ylm(θ, φ)

l = 3 Y33 = −1

4

√35

4πsin3 θ e3iφ

Y32 =1

4

√105

2πsin2 θ cos θe2iφ

Y31 = −1

4

√21

4πsin θ(5 cos2 θ − 1) e iφ

Y30 =1

2

√7

4πsin θ

(5 cos3 θ − 3 cos θ

)(65)

Classical Field Theory: Electrostatics-Magnetostatics

Page 27: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Figure: Schematic representation ofYlm on the unit sphere. Ylm is equalto 0 along m great circles passingthrough the poles, and along l −mcircles of equal latitude. Thefunction changes sign each time itcrosses one of these lines.

Classical Field Theory: Electrostatics-Magnetostatics

Page 28: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Addition Theorem for Spherical Harmonics

The spherical harmonics are related toLegendre polynomials Pl by a relation known asthe addition theorem.

The addition theorem expresses a Legendrepolynomial of order l in the angle γ in terms ofproducts of the spherical harmonics of theangles θ, φ and θ′, φ′:

Pl(cos γ) =4π

l(l + 1)

l∑m=−l

Y ∗lm(θ′, φ′)Ylm(θ, φ) (66)

where γ is the angle between the vectors ~x and ~x ′, ~x · ~x ′ = x · x ′ cos γand cos γ = cos θ cos θ′ + sin θ sin θ′ cos(φ− φ′).Pl(cos γ) is a function of the angles θ, φ with the angles θ′, φ′ asparameters and it maybe expanded in a series (63) :

Pl(cos γ) =∞∑

l′=0

l′∑m=−l′

Al′m′(θ′, φ′)Yl′m′(θ, φ) (67)

Classical Field Theory: Electrostatics-Magnetostatics

Page 29: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

The addition theorem offer the possibility to extend the expansion validfor a point charge (axially symmetric distribution) to an arbitrary chargedistribution.By substituting (66) into (53) we obtain

1

|~x − ~x ′|= 4π

∞∑l=0

l∑m=−l

1

2l + 1

r l<

r l+1>

Y ∗lm(θ′, φ′)Ylm(θ, φ) (68)

This equation gives the potential in a completely factorized form in thecoordinates ~x and ~x ′. This is useful in any integration over chargedensities, when one is the variable of integration and the other theobservation point.

Classical Field Theory: Electrostatics-Magnetostatics

Page 30: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Multipole Expansion

A localized distribution of charge is described bythe charge density ρ(~x ′), which is nonvanishingonly inside a sphere a around some origin.The potential outside the sphere can be writtenas an expansion in spherical harmonics

Φ(~x) =∞∑l=0

l∑m=−l

2l + 1qlm

1

r l+1Ylm(θ, φ) (69)

This type of expansion is called multipole expansion;The l = 0 term is called monopole term,The l = 1 are called dipole terms etc.The problem to be solved is the determination of the constants qlm interms of the properties of the charge distribution ρ(~x ′).The solution is very easily obtained from the integral for the potential

Φ(~x) =

∫ρ(~x ′)

|~x − ~x ′|d3x ′ (70)

Classical Field Theory: Electrostatics-Magnetostatics

Page 31: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Using the expansion (68) for 1/|~x − ~x ′| i.e.

Φ(~x) =∞∑l=0

l∑m=−l

2l + 1

1

r l+1Ylm(θ, φ)

∫Y ∗lm(θ′, φ′)ρ(~x ′)r ′ld3x ′ (71)

Consequently the coefficients in (69) are :

qlm =

∫Y ∗lm(θ′, φ′)ρ(~x ′)r ′ld3x ′ (72)

and called the multipole moments of the charge distribution ρ(~x ′)

Classical Field Theory: Electrostatics-Magnetostatics

Page 32: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Monopole moment l = 0

Here, the only component is

q00 =

∫ρ(~x ′)r ′0Y ∗00(θ′, φ′)d3x ′ =

1√4π

∫ρ(~x ′)d3x ′ =

q√4π

(73)

Observed from a large distance r , any charge distribution actsapproximately as if the total charge q (monopole moment) would beconcentrated at one point since the dominant term in (68)

Φ(~x) = 4πq∗00

1

rY00 + · · · =

q

r+ . . . (74)

NOTE:The moments with m ≥ 0 are connected (for real charge density) too themoments with m < 0 through

ql,−m = (−1)mq∗lm (75)

Classical Field Theory: Electrostatics-Magnetostatics

Page 33: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Dipole moment l = 1

In Cartesian coordinates the dipole moment is given by

~p =

∫~x ′ρ(~x ′)d3x ′ (76)

In Spherical representation one obtains

q∗11 =

∫ρ(~x ′)r ′Y ∗11(θ′, φ′)d3x ′

= −√

3

∫ρ(~x ′)r ′ (sin θ′ cosφ′ − i sin θ′ sinφ′) d3x ′

in Cartesian represantation

= −√

3

∫ρ(~x ′)(x ′ − iy ′)d3x ′ = −

√3

8π(px − ipy ) (77)

also

q∗10 =

∫ρ(~x ′)r ′Y ∗10(θ′, φ′)d3x ′ =

√3

∫ρ(~x ′)r ′ cos θ′d3x ′

=

√3

∫z ′ρ(~x ′)d3x =

√3

4πpz (78)

Classical Field Theory: Electrostatics-Magnetostatics

Page 34: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Quadrupole moment l = 2

q∗22 =

∫ρ(~x)r ′2Y ∗22(θ′, φ′)d3x ′

=1

4

√15

∫ρ(~x) [r ′ sin θ′(cosφ′ − i sinφ′)]

2d3x ′

=1

4

√15

∫ρ(~x)(x ′ − iy ′)2d3x ′

because (x ′ − iy ′)2 =1

3

[(3x ′2 − r ′2)− 6ix ′y ′ − (3y ′2 − r ′2)

]=

1

12

√15

2π(Q11 − 2iQ12 − Q22) (79)

where Qij is the traceless (why?) quadrupole moment tensor:

Qij =

∫ (3x ′i x

′j − r2δij

)ρ(~x ′)d3x ′ (80)

Classical Field Theory: Electrostatics-Magnetostatics

Page 35: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Quadrupole moment l = 2

Analogously

q∗21 =

∫ρ(~x ′)r ′2Y ∗21(θ′, φ′)d3x ′ = −

√15

∫ρ(~x ′)z ′(x ′ − iy ′)d3x ′

= −1

3

√15

8π(Q13 − iQ23) (81)

and

q∗20 =

∫ρ(~x ′)r ′2Y ∗20(θ′, φ′)d3x ′ =

1

2

√5

∫ρ(~x ′)(3z ′ − r ′2)d3x ′

=1

2

√5

4πQ33 (82)

From eqn (60) we can get the moments with m < 0 through

ql,−m = (−1)mq∗lm (83)

Classical Field Theory: Electrostatics-Magnetostatics

Page 36: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Multipole Expansion

By direct Taylor expansion of 1/|~x − ~x ′| the expansion of Φ(~x) inrectangular coordinates is:

Φ(~x) =q

r+~p · ~xr3

+1

2

∑i,j

Qijxixj

r5+ . . . (84)

The electric field components for a given multipole can be expressedmost easily in terms of spherical coordinates. From ~E = −~∇Φ and (69)for fixed (l ,m) we get(how?):

Er =4π(l + 1)

2l + 1qlm

1

r l+2Ylm

Eθ = − 4π

2l + 1qlm

1

r l+2

∂θYlm (85)

Eφ = − 4π

2l + 1qlm

1

r l+2

1

sin θ

∂φYlm

For a dipole ~p along the z-axis they reduce to:

Er =2p cos θ

r3, Eθ =

p sin θ

r3, Eφ = 0 (86)

Classical Field Theory: Electrostatics-Magnetostatics

Page 37: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Multipole Expansion

The field intensity at a point ~x due to a dipole ~p at a point ~x0

(r = |~x − ~x0|) is

~E (~x) =3~n (~p · ~n)− ~p|~x − ~x0|3

with ~n =~r

r(87)

because Φ(~x) = ~p · ~x/r3 and

~E (~x) = −~∇Φ = −∇(~p · ~x)

r3− (~p · ~x)∇

(1

r3

)=

3~n (~p · ~n)− ~p|~x − ~x0|3

(88)

Classical Field Theory: Electrostatics-Magnetostatics

Page 38: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Energy of a Charge Distribution in an External Field

The multipole expansion of the potential of acharge distribution can also be used todescribe the interaction of the chargedistribution with an external field.The electrostatic energy of the chargedistribution ρ(~x) placed in an external fieldΦ(~x) is given by

W =

∫V

ρ(~x)Φ(~x)d3x (89)

The external field (if its is slowly varying over the region where ρ(~x) isnon-negligible) may be expanded in a Taylor series:

Φ(~x) = Φ(0) + ~x · ~∇Φ(0) +1

2

3∑i=1

3∑j=1

xixj∂2Φ

∂xi∂xj(0) + . . . (90)

Since ~E = −~∇Φ for the external field

Φ(~x) = Φ(0)− ~x · ~E (0)− 1

2

3∑i=1

3∑j=1

xixj∂Ej

∂xi(0) + . . . (91)

Classical Field Theory: Electrostatics-Magnetostatics

Page 39: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Since ~∇ · ~E = 0 for the external field we can substract

1

6r2~∇ · ~E (0)

from the last term to obtain finally the expansion:

Φ(~x) = Φ(0)− ~x · ~E (0)− 1

6

3∑i=1

3∑j=1

(3xixj − r2δij

) ∂Ej

∂xi(0) + . . . (92)

When this inserted into (89) the energy takes the form:

W = qΦ(0)− ~p · ~E (0)− 1

6

3∑i=1

3∑j=1

Qij∂Ej

∂xi(0) + . . . (93)

Notice, that:

I the total charge interacts with the potential,

I the dipole moment with the electric field,

I the quadrupole with the electric field gradient, etc

Classical Field Theory: Electrostatics-Magnetostatics

Page 40: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Examples

EXAMPLE 1 : Show that for a uniform charged sphere all multipolemoments vanish except q00.

If the sphere has a radius R0 and constant charge density ρ(~x) = ρ then

q∗lm = ρ

∫Ω′

∫ R0

0

r ′lY ∗lm(θ′, φ′)r ′2dr ′dΩ′ = ρR l+3

0

l + 3

∫Ω′

Y ∗lm(θ′, φ′)dΩ′

but since Y00 = 1/√

4π we have from the orthogonality relation

q∗lm = ρR l+3

0

l + 3

√4π

∫Ω′

Y ∗lmY00dΩ′ = ρR l+3

0

l + 3

√4πδl0δm0

Classical Field Theory: Electrostatics-Magnetostatics

Page 41: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

EXAMPLE 2: Perform multipole decomposition of a uniform chargedistribution whose surface is a weakly deformed sphere:

R = R0

(1 +

2∑m=−2

a∗2mY2m(θ, φ)

), |a2m| << 1

The multipole moments are

q∗lm =

∫Ω′

∫ R(θ′,φ′)

0

ρr ′lY ∗lmr ′2dr ′dΩ′

= ρR l+3

0

l + 3

∫Ω′

Y ∗lm(θ′, φ′)

(1 +

2∑m=−2

a∗2mY ∗2m(θ′, φ′)

)l+3

dΩ′

Classical Field Theory: Electrostatics-Magnetostatics

Page 42: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

q∗lm ≈ ρR l+3

0

l + 3

∫Ω′

Y ∗lm(θ′, φ′)

(1 + (l + 3)

2∑m=−2

a∗2mY ∗2m(θ′, φ′)

)dΩ′

Apart from the monopole moment q00 of the previous example, we find

q∗lm = ρR l+3

0

l + 3(l+3)

2∑µ=−2

a∗2µ

∫Ω′

Y ∗lmY2µdΩ′ = ρR l+30

∑µ

a∗2µδmµδl2 (l > 0)

Thus the nonvanishing multipole moments are : q∗2m = ρR50 a∗2m.

Classical Field Theory: Electrostatics-Magnetostatics

Page 43: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Magnetostatics

• There is a radical difference between magnetostatics andelectrostatics : there are no free magnetic charges• The basic entity in magnetic studies is the magnetic dipole.• In the presence of magnetic materials the dipole tends to align itselfin a certain direction. That is by definition the direction of the magneticflux density , denoted by ~B (also called magnetic induction).• The magnitude of the flux density can be defined by the magnetictorque ~N exerted on the magnetic dipole:

~N = ~µ× ~B (94)

where ~µ is the magnetic moment of the dipole.• In electrostatics the conservation of charge demands that the chargedensity at any point in space be related to the current density in thatneighborhood by the continuity equation

∂ρ/∂t + ~∇ · ~J = 0 (95)

• Steady-state magnetic phenomena are characterized by no change inthe net charge density anywhere in space, consequently in magnetostatics

~∇ · ~J = 0 (96)Classical Field Theory: Electrostatics-Magnetostatics

Page 44: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Biot & Savart Law

If d ~ is an element of length (pointing in thedirection of current flow) of a wire which carriesa current I and ~x is the coordinate vector fromthe element of length to an observation pointP, then the elementary flux density d~B at thepoint P is given in magnitude and direction by

d~B = kId ~× ~x|~x |3

(97)

NOTE: (97) is an inverse square law, just as is Coulomb’s law ofelectrostatics. However, the vector character is very different.

• If we consider that current is charge in motion and replace Id ~ by q~vwhere q is the charge and ~v the velocity. The flux density for such acharge in motion would be

~B = kq~v × ~x|~x |3

(98)

This expression is time-dependent and valid only for charges with small

velocities compared to the speed of light.Classical Field Theory: Electrostatics-Magnetostatics

Page 45: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Diff. Equations for Magnetostatics & Ampere’s Law

The basic law (97) for the magnetic induction can be written down in

general form for a current density ~J(~x):

~B(~x) =1

c

∫~J(~x ′)× (~x − ~x ′)

|~x − ~x ′|3d3x ′ (99)

This expression for ~B(~x) is the magnetic analog of electric field in termsof the charge density:

~E (~x) =1

c

∫ρ(~x ′)

(~x − ~x ′)|~x − ~x ′|3

d3x ′ (100)

In order to obtain DE equivalent to (99) we use the relation

(~x − ~x ′)|~x − ~x ′|3

= −~∇(

1

|~x − ~x ′|

)as ~∇

(1

r

)= −

~r

r3(101)

and (99) transforms into

~B(~x) =1

c~∇×

∫ ~J(~x ′)

|~x − ~x ′|d3x ′ (102)

Classical Field Theory: Electrostatics-Magnetostatics

Page 46: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Classical Field Theory: Electrostatics-Magnetostatics

Page 47: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

From (102) follows immediately:

~∇ · ~B = 0 (103)

• This is the 1st equation of magnetostatics and corresponds to~∇× ~E = 0 in electrostatics.By analogy with electrostatics we now calculate the curl of ~B

~∇× ~B =1

c~∇× ~∇×

∫ ~J(~x ′)

|~x − ~x ′|d3x ′ (104)

which for steady-state phenomena (~∇ · ~J = 0) reduces to (how?) 2

~∇× ~B =4π

c~J (105)

• This is the 2nd equation of magnetostatics and corresponds to~∇ · ~E = 4πρ in electrostatics.

2Remember : ~∇ד

~∇× ~A”

= ~∇“

~∇ · ~A”−∇2~A

Classical Field Theory: Electrostatics-Magnetostatics

Page 48: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

The integral equivalent of (105) is called Ampere’s law

It can be obtained by applying Stokes’stheorem to the integral of the normalcomponent of (105) over the open surface Sbounded by a closed curve C . Thus∫

S

~∇× ~B · ~n da =4π

c

∫S

~J · ~nda (106)

is transformed into∮C

~B · d ~=4π

c

∫S

~J · ~nda (107)

Since the surface integral of the current density is the total current Ipassing through the closed curve C , Ampere’s law can be written in theform: ∮

C

~B · d ~=4π

cI (108)

Classical Field Theory: Electrostatics-Magnetostatics

Page 49: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Vector Potential

The basic differential laws in magnetostatics are

~∇× ~B =4π

c~J and ~∇ · ~B = 0 (109)

The problems is how to solve them.• If the current density is zero in the region of interest, ~∇× ~B = 0permits the expression of the vector magnetic induction ~B as the gradientof a magnetic scalar potential ~B = −~∇ΦM , then (109) reduces to theLaplace equation for ΦM .• If ~∇ · ~B = 0 everywhere, ~B must be the curl of some vector field~A(~x), called the vector potential

~B(~x) = ~∇× ~A(~x) (110)

and from (102) the general form of ~A is:

~A(~x) =1

c

∫ ~J(~x ′)

|~x − ~x ′|d3x ′ + ~∇Ψ(~x) (111)

Classical Field Theory: Electrostatics-Magnetostatics

Page 50: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

The added gradient of an arbitrary scalar function Ψ shows that, for agiven magnetic induction ~B, the vector potential can be freelytransformed according to

~A→ ~A + ~∇Ψ (112)

This transformation is called a gauge transformation. Suchtransformations are possible because (110) specifies only the curl of ~A.

• If (110) is substituted into the first equation in (109), we find

~∇×(~∇× ~A

)=

c~J or ~∇

(~∇ · ~A

)−∇2~A =

c~J (113)

If we exploit the freedom implied by (112) we can make the convenient

choice of gauge (Coulomb gauge) ~∇ · ~A = 0. Then each component ofthe vector potential satisfies the Poisson equation

∇2~A = −4π

c~J (114)

The solution for ~A in unbounded space is (111) with Ψ =constant:

~A(~x) =1

c

∫ ~J(~x ′)

|~x − ~x ′|d3x ′ (115)

Classical Field Theory: Electrostatics-Magnetostatics

Page 51: Classical Field Theory: Electrostatics-MagnetostaticsElectrostatics The behavior of an electrostatic eld can be described by two di erential equations: r~ E~= 4ˇˆ (1) (Gauss law)

Further Reading

I Green Function

I Greiner W. ”Classical Electrodynamics”, Chapter 1, 2, 3, 4I Jackson J.K ”Classical Electrodynamics”, 2nd Edition,

Sections 1.7-1.10, 3.1-3.6, 4.1-4.2, 5.1-5.4

Classical Field Theory: Electrostatics-Magnetostatics