​
ElectroDynamics, by James Rohlf and Kevin Reiss, Version 1.0, Copyright 2021, All Rights Reserved.
__________________________________________________________________________________

Excerpts

The electric dipole plays a fundamental role in electromagnetism, from statics to radiation, and the electric potential is an invaluable concept for its visualization and interpretation.
Out[193]=
Electric Dipole Potential
By:
Stephen Wolfram
​
type of plot
3D plot
contour plot
range of plot
show field direction
position of charge 1
strength of charge 1
position of charge 2
strength of charge 2
​
__________________________________________________________________________________
Calculate the ratio of electric to gravitational force for an electron-proton interaction.
In[28]:=
ClearAll["Global`*"];UnitConvert[(
2
e
4π
ε
0
)/(
G
proton
(particle)
[
mass
]
electron
(particle)
[
mass
])]
Out[28]=
2.269×
39
10
Draw the line charge integration geometry off the symmetry axis.
Out[495]=
Calculate the electric field at (x, 0, z) from the line of charge above.
In[8]:=
r={x,0,z};r′={x′,0,0};ℛ=r-r′;$Assumptions={L>0,z∈,z≠0,x∈};​​Ε=
λ
4π
ε
0
L/2
∫
-L/2
ℛ
3/2
(ℛ.ℛ)
x′//Together//FullSimplify
Out[9]=

-
2
(L-2x)
+4
2
z
+
2
(L+2x)
+4
2
z
λ
2π
4
L
+8
2
L
(-
2
x
+
2
z
)+16
2
(
2
x
+
2
z
)
ε
0
,0,
2x
2
(L-2x)
+4
2
z
-
2
(L+2x)
+4
2
z
+L
2
(L-2x)
+4
2
z
+
2
(L+2x)
+4
2
z
λ
4πz
4
L
+8
2
L
(-
2
x
+
2
z
)+16
2
(
2
x
+
2
z
)
ε
0

Calculate the electric field for
L=0.1
m
,x=0.4
m
,z=.3
m
,andλ=1
μC

m
.
In[11]:=
UnitConvert[Ε/.{L->0.1
m
,x->0.4
m
,z->.3
m
,λ->1
μC
/
m
,
ε
0

ε
0
,0->0
V
/
m
},
V
/
m
]
Out[11]=
{
2878.75
V/m
,
0
V/m
,
2180.79
V/m
}
Plot the electric field in the x-z plane.
In[26]:=
VectorPlot[Ε[[{1,3}]]/.{L2,λ1,
ε
0
1},{x,-2,2},{z,-2,2},EpilogLine[{{-1,0},{1,0}}],
]
Out[26]=
__________________________________________________________________________________
Draw the geometry for dividing a solid ball into disks.
Out[148]=
Calculate the electric field outside a uniform ball of charge, by adding together disks.
In[194]:=
ClearAll["Global`*"];$Assumptions={z′∈Reals,z′>0,r∈Reals,r>R,R>0};ρ=
Q
(4/3)π
3
R
;
ρ
2
ε
0
R
∫
-R
(1-
r-z′
2
R
-
2
z′
+
2
(r-z′)
)z′
Out[194]=
Q
4π
2
r
ε
0
__________________________________________________________________________________
Calculate the electric field of an ideal dipole.
Transform the electric field to Cartesian coordinates.
Verify that the coordinate transformation produces the same result as taking the gradient in Cartesian.
Verify that the curl is zero.
Draw the electric field lines for the ideal dipole.
Draw the electric field lines for an ideal dipole in 3D.
■
Wolfram Research (2021), StreamPlot3D, Wolfram Language function, https://reference.wolfram.com/language/ref/StreamPlot3D.html.
__________________________________________________________________________________
Draw the boundary conditions for the solution of Laplace’s Equation.
Calculate the full form of the potential inside the above region.
Plot the solution for a = 2 and b = 1, keeping 500 terms.
__________________________________________________________________________________
The following demonstration solves Laplace’s equation for arbitrary boundary conditions on various geometries. A certain sinusoidal boundary condition produces a solution in the shape of a “Pringle” chip.
__________________________________________________________________________________
Calculate the magnetic field due to a circular loop of radius R and current ℐ at a distance z from the center along the axis of symmetry.
__________________________________________________________________________________

2.9 Summary of B, A, and J

Here is a summary of how to calculate one quantity from another in magneto-statics. The hardest conversion to do is to get A from B because there is no easy recipe for undoing a curl. See, for example, Section 2.5.
__________________________________________________________________________________
A dielectric ball placed in a uniform external electric field is an important example that will appear again in the discussion of magnetic fields in matter. The result is a uniform electric field inside the ball, which contributes to the external field as a pure dipole. It is described in Section 3.4.
__________________________________________________________________________________
__________________________________________________________________________________
When squaring the electric field, one may note that the cross term is twice as large as the second term squared, and has the opposite sign. The Poynting vector is
Make a polar plot of the Poynting vector (energy flux).
Plot the 3D radiation pattern.
__________________________________________________________________________________
Make a polar plot of the radiation pattern for β=.9.
Calculate 1-β.
Calculate the radiated power per proton.
Calculate the radiated power per beam.
Calculate the radiated energy per proton per revolution.
__________________________________________________________________________________
Make a polar plot of the electric quadrupole energy flux.
Make a 3D plot of the electric quadrupole energy flux.
__________________________________________________________________________________

8.5 Magnetic Quadrupole Radiation

A magnetic quadruple can be constructed out of 2 out-of-phase magnetic dipoles. Consider 2 current loops of radius d separated by distance b.
This time the first 2 terms of the expansion of the exponential both give zero and we must use the 3rd term. We have separated the current into 2 loops with opposite sign.
Calculate the magnetic quadrupole angular term.
__________________________________________________________________________________
Display the contravariant form of the field tensor.
Display the covariant form of the field tensor.
Transform the field tensor.
Transform back.