We start off with Maxwell's Equation in the Lorentz gauge:
◻
2
A
μ
(
r
,
t
)
=
◻
2
A
μ
(
r
)
=
(
−
μ
1
j
μ
(
r
)
)
{\displaystyle \square ^{2}A^{\mu }(\mathbf {r} ,t)=\square ^{2}A^{\mu }(r)=(-\mu _{1}j^{\mu }(r))}
Where:
A
μ
=
(
A
,
Φ
c
)
,
◻
2
=
∂
μ
∂
μ
=
∇
2
−
1
c
2
∂
∂
t
2
{\displaystyle A^{\mu }=(\mathbf {A} ,{\frac {\Phi }{c}}),\square ^{2}=\partial _{\mu }\partial ^{\mu }=\nabla ^{2}-{\frac {1}{c^{2}}}{\frac {\partial }{\partial t^{2}}}}
j
μ
=
(
j
,
c
ρ
)
,
∂
μ
=
(
∇
,
1
c
∂
∂
t
)
{\displaystyle j^{\mu }=(\mathbf {j} ,c\rho ),\partial _{\mu }=(\mathbf {\nabla } ,{\frac {1}{c}}{\frac {\partial }{\partial t}})}
Lorentz Gauge:
A
μ
=
0
⇒
∇
⋅
A
+
1
c
2
∂
Φ
∂
t
=
0
{\displaystyle A^{\mu }=0\Rightarrow \mathbf {\nabla } \cdot \mathbf {A} +{\frac {1}{c^{2}}}{\frac {\partial \Phi }{\partial t}}=0}
Introduce Green's function at
(
r
,
t
)
=
r
{\displaystyle (\mathbf {r} ,t)=r\quad }
from some impulse source at
r
′
=
(
r
′
,
t
′
)
{\displaystyle r'=(\mathbf {r} ',t')\quad }
◻
r
2
G
(
r
,
r
′
)
=
δ
4
(
r
−
r
′
)
{\displaystyle \square _{r}^{2}G(r,r')=\delta ^{4}(r-r')}
Let
G
~
(
q
)
=
1
(
2
π
)
2
∫
d
4
r
e
i
q
r
G
(
r
,
0
)
{\displaystyle {\tilde {G}}(q)={\frac {1}{(2\pi )^{2}}}\int d^{4}re^{iqr}G(r,0)}
Then
G
(
q
)
=
1
(
2
π
)
2
∫
d
4
q
e
i
q
r
G
~
(
r
,
0
)
{\displaystyle G(q)={\frac {1}{(2\pi )^{2}}}\int d^{4}qe^{iqr}{\tilde {G}}(r,0)}
In free space, translational symmetry implies:
G
(
r
−
r
′
,
0
)
=
G
(
r
,
r
′
)
{\displaystyle G(r-r',0)=G(r,r')\quad }
∴
G
(
r
,
r
′
)
=
1
(
2
π
)
2
∫
d
4
q
e
−
i
q
(
r
−
r
′
)
G
~
(
q
)
{\displaystyle G(r,r')={\frac {1}{(2\pi )^{2}}}\int d^{4}qe^{-iq(r-r')}{\tilde {G}}(q)}
◻
r
2
G
(
r
,
r
′
)
=
1
(
2
π
)
2
|
i
n
t
d
4
q
e
−
i
q
(
r
−
r
′
)
G
~
(
q
)
{\displaystyle \square _{r}^{2}G(r,r')={\frac {1}{(2\pi )^{2}}}|intd^{4}qe^{-iq(r-r')}{\tilde {G}}(q)}
◻
r
2
G
(
r
,
r
′
)
=
1
(
2
π
)
2
∫
d
4
q
e
−
i
q
(
r
−
r
′
)
(
−
k
2
+
ω
2
c
2
)
{\displaystyle \square _{r}^{2}G(r,r')={\frac {1}{(2\pi )^{2}}}\int d^{4}qe^{-iq(r-r')}(-k^{2}+{\frac {\omega ^{2}}{c^{2}}})}
, where
q
=
(
k
,
ω
c
)
{\displaystyle q=(\mathbf {k} ,{\frac {\omega }{c}})\quad }
But,
◻
r
2
G
(
r
,
r
′
)
=
δ
4
(
r
−
r
′
)
=
1
(
2
π
)
4
∫
d
4
q
e
−
i
q
(
r
−
r
′
)
{\displaystyle \square _{r}^{2}G(r,r')=\delta ^{4}(r-r')={\frac {1}{(2\pi )^{4}}}\int d^{4}qe^{-iq(r-r')}}
∴Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \tilde{G}(q)=\frac{(2\pi)^2}{(2\pi)^4}\frac{1}{-q^2}= \frac{-1}{(2\pi)^2q^2}}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G(r,r')=\frac{-1}{(2\pi)^4} \int d^4qe^{-iq(r-r')} \frac{1}{(k+\frac{\omega}{c})(k-\frac{\omega}{c})}}
Chose the "retarded" solution, such that the function is zero unless t>t'
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G(r,r')=\frac{1}{(2\pi)^4}\int d^3ke^{-i\mathbf{k}(r-r')}\int d(\frac{\omega}{c}) \frac{e^{i\omega(t-t')}}{(\frac{\omega}{c}-k)(\frac{\omega}{c}+k)}\Theta}
=
1
(
2
π
)
4
∫
d
3
k
e
−
i
k
(
r
−
r
′
)
(
2
π
i
e
i
c
k
(
t
−
t
′
)
−
e
−
i
c
k
(
t
−
t
′
)
2
k
)
Θ
{\displaystyle ={\frac {1}{(2\pi )^{4}}}\int d^{3}ke^{-i\mathbf {k} (r-r')}(2\pi i{\frac {e^{ick(t-t')}-e^{-ick(t-t')}}{2k}})\Theta }
=
−
2
π
(
2
π
)
4
∫
k
2
d
k
k
sin
(
c
k
(
t
−
t
′
)
)
2
π
∫
−
i
i
d
z
e
−
i
k
|
r
−
r
′
|
z
Θ
{\displaystyle ={\frac {-2\pi }{(2\pi )^{4}}}\int {\frac {k^{2}dk}{k}}\sin \left({ck(t-t')}\right)2\pi \int _{-i}^{i}dze^{-ik|\mathbf {r} -\mathbf {r'} |z}\Theta }
=
−
1
(
2
π
)
2
(
1
2
i
|
r
−
r
′
|
)
2
∫
d
k
s
i
n
(
c
k
(
t
−
t
′
)
)
s
i
n
(
k
|
r
−
r
′
|
)
Θ
{\displaystyle ={\frac {-1}{(2\pi )^{2}}}\left({\frac {1}{2i|\mathbf {r} -\mathbf {r'|} }}\right)2\int dksin(ck(t-t'))sin(k|\mathbf {r} -\mathbf {r'} |)\Theta }
=
1
(
2
π
)
2
2
|
r
−
r
′
|
2
π
4
[
2
δ
(
|
r
−
r
′
|
+
c
(
t
−
t
′
)
)
−
2
δ
(
|
r
−
r
′
|
−
c
(
t
−
t
′
)
)
]
Θ
{\displaystyle ={\frac {1}{(2\pi )^{2}}}{\frac {2}{|\mathbf {r} -\mathbf {r} '|}}{\frac {2\pi }{4}}\left[2\delta (|\mathbf {r} -\mathbf {r} '|+c(t-t'))-2\delta (|\mathbf {r} -\mathbf {r} '|-c(t-t'))\right]\Theta }
But the term
2
δ
(
|
r
−
r
′
|
+
c
(
t
−
t
′
)
)
→
0
∀
t
>
t
′
{\displaystyle 2\delta (|\mathbf {r} -\mathbf {r} '|+c(t-t'))\rightarrow 0\quad \forall \quad t>t'}
∴
G
(
r
,
r
′
)
=
−
1
4
π
δ
(
|
r
−
r
′
|
−
c
(
t
−
t
′
)
)
|
r
−
r
′
|
{\displaystyle G(r,r')={\frac {-1}{4\pi }}\quad {\frac {\delta (|\mathbf {r} -\mathbf {r} '|-c(t-t'))}{|\mathbf {r} -\mathbf {r} '|}}}
Now to get the
G
1
(
r
,
r
′
)
{\displaystyle G_{1}(r,r')\quad }
in the half-space with z>0 with the boundary condition
G
1
{\displaystyle G_{1}\quad }
at
r
3
=
z
=
0
{\displaystyle r_{3}=z=0\quad }
we take the difference:
G
1
(
r
,
r
′
)
=
−
1
4
π
(
δ
(
|
r
−
r
′
|
−
c
(
t
−
t
′
)
)
|
r
−
r
′
|
−
δ
(
|
r
−
r
′
−
2
z
e
3
^
|
−
c
(
t
−
t
′
)
)
|
r
−
r
′
−
2
z
e
3
^
|
)
{\displaystyle G_{1}(r,r')={\frac {-1}{4\pi }}\left({\frac {\delta (|\mathbf {r} -\mathbf {r} '|-c(t-t'))}{|\mathbf {r} -\mathbf {r} '|}}-{\frac {\delta (|\mathbf {r} -\mathbf {r} '-2z{\hat {e_{3}}}|-c(t-t'))}{|\mathbf {r} -\mathbf {r} '-2z{\hat {e_{3}}}|}}\right)}
Now use Green's theorem:
Let
F
=
A
(
r
)
∇
G
1
(
r
,
r
′
)
−
G
1
(
r
,
r
′
)
∇
A
(
r
)
{\displaystyle \mathbf {F} =A(r)\mathbf {\nabla } G_{1}(r,r')-G_{1}(r,r')\mathbf {\nabla } A(r)}
∫
∇
⋅
F
d
4
r
=
∫
c
d
t
∫
d
3
r
[
∇
A
⋅
∇
G
+
A
∇
2
G
1
−
∇
G
⋅
∇
A
−
G
1
∇
2
A
]
{\displaystyle \int \mathbf {\nabla } \cdot \mathbf {F} d^{4}r=\int cdt\int d^{3}r[\mathbf {\nabla } A\cdot \mathbf {\nabla } G+A\nabla ^{2}G_{1}-\mathbf {\nabla } G\cdot \mathbf {\nabla } A-G_{1}\nabla ^{2}A]}
But
∇
2
G
1
(
r
,
r
′
)
=
δ
4
(
r
−
r
′
)
+
1
c
2
∂
2
∂
t
2
G
1
(
r
,
r
′
)
{\displaystyle \nabla ^{2}G_{1}(r,r')=\delta ^{4}(r-r')+{\frac {1}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}}G_{1}(r,r')}
∇
2
A
(
r
)
=
μ
j
(
r
)
+
1
c
2
∂
2
∂
t
2
A
(
r
)
{\displaystyle \nabla ^{2}A(r)=\mu j(r)+{\frac {1}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}}A(r)}
, let
j
(
r
)
=
0
{\displaystyle j(r)=0\quad }
∫
∇
⋅
F
d
4
r
=
A
(
r
′
)
+
1
c
2
∫
d
4
r
[
A
∂
2
∂
t
2
G
1
−
G
1
∂
2
∂
t
2
A
]
{\displaystyle \int \nabla \cdot \mathbf {F} d^{4}r=A(r')+{\frac {1}{c^{2}}}\int d^{4}r\left[A{\frac {\partial ^{2}}{\partial t^{2}}}G_{1}-G_{1}{\frac {\partial ^{2}}{\partial t^{2}}}A\right]}
The last term vanishes if
G
1
(
r
,
r
′
)
a
n
d
A
(
r
)
{\displaystyle G_{1}(r,r')andA(r)\quad }
fall off sufficiently fast at
t
→
∞
{\displaystyle t\rightarrow \infty }
. They do. So:
∫
∇
⋅
F
d
4
r
=
A
(
r
′
)
{\displaystyle \int \nabla \cdot \mathbf {F} d^{4}r=A(r')}
Now invoke the divergence theorem on the half space
z
>
0
{\displaystyle z>0\quad }
:
A
(
r
′
)
=
−
∫
d
2
r
∫
c
d
t
[
A
(
r
)
∂
∂
t
G
1
(
r
,
r
′
)
−
G
1
(
r
,
r
′
)
∂
∂
z
A
(
r
)
]
{\displaystyle A(r')=-\int d^{2}r\int cdt\left[A(r){\frac {\partial }{\partial t}}G_{1}(r,r')-G_{1}(r,r'){\frac {\partial }{\partial z}}A(r)\right]}
, where the last term is zero by the constriction of
G
1
(
z
=
0
)
{\displaystyle G_{1}(z=0)\quad }
A
(
r
′
)
=
−
c
∫
d
t
∫
d
2
r
A
(
r
)
∂
∂
z
G
1
(
r
,
r
′
)
{\displaystyle A(r')=-c\int dt\int d^{2}rA(r){\frac {\partial }{\partial z}}G_{1}(r,r')}
To do the t integral, I need to bring out the z derivative. To do this, I first turn it into a z' derivative, using the relation:
G
1
(
r
,
r
′
)
=
−
1
4
π
(
δ
(
|
r
−
r
′
|
−
c
(
t
−
t
′
)
)
|
r
−
r
′
|
−
δ
(
|
r
−
r
″
|
−
c
(
t
−
t
′
)
)
|
r
−
r
″
|
)
{\displaystyle G_{1}(r,r')={\frac {-1}{4\pi }}\left({\frac {\delta (|\mathbf {r} -\mathbf {r} '|-c(t-t'))}{|\mathbf {r} -\mathbf {r} '|}}-{\frac {\delta (|\mathbf {r} -\mathbf {r} ''|-c(t-t'))}{|\mathbf {r} -\mathbf {r} ''|}}\right)}
, where
r
″
=
r
′
−
2
z
′
e
3
^
{\displaystyle \mathbf {r} ''=\mathbf {r} '-2z'{\hat {e_{3}}}}
∂
∂
z
G
1
(
r
,
r
′
)
=
1
4
π
(
∂
∂
z
(
δ
(
|
r
−
r
′
|
−
c
(
t
−
t
′
)
)
|
r
−
r
′
|
−
δ
(
|
r
−
r
″
|
−
c
(
t
−
t
′
)
)
|
r
−
r
″
|
)
)
{\displaystyle {\frac {\partial }{\partial z}}G_{1}(r,r')={\frac {1}{4\pi }}\left({\frac {\partial }{\partial z}}\left({\frac {\delta (|\mathbf {r} -\mathbf {r} '|-c(t-t'))}{|\mathbf {r} -\mathbf {r} '|}}-{\frac {\delta (|\mathbf {r} -\mathbf {r} ''|-c(t-t'))}{|\mathbf {r} -\mathbf {r} ''|}}\right)\right)}
∴
A
(
r
′
)
=
−
1
4
π
∂
∂
z
′
∫
z
=
0
d
2
r
(
2
A
(
r
,
t
′
−
r
−
r
′
c
r
−
r
′
)
{\displaystyle A(r')={\frac {-1}{4\pi }}{\frac {\partial }{\partial z'}}\int _{z=0}d^{2}r\left(2{\frac {A(\mathbf {r} ,t'-{\frac {\mathbf {r} -\mathbf {r} '}{c}}}{\mathbf {r} -\mathbf {r} '}}\right)}
At
z
=
0
{\displaystyle z=0\quad }
,
|
r
−
r
′
|
=
r
2
+
z
′
2
=
S
,
d
S
=
r
d
r
r
2
+
z
′
2
{\displaystyle |\mathbf {r} -\mathbf {r} '|={\sqrt {r^{2}+z'^{2}}}=S,dS={\frac {rdr}{\sqrt {r^{2}+z'^{2}}}}}
If
A
(
r
,
t
)
{\displaystyle A(\mathbf {r} ,t)\quad }
is independent of
r
{\displaystyle \mathbf {r} \quad }
, then:
A
(
r
′
)
=
−
∂
∂
z
′
∫
z
′
=
0
∞
d
S
A
(
0
,
t
−
S
c
)
=
A
(
∅
,
t
′
−
z
′
c
)
{\displaystyle A(r')={\frac {-\partial }{\partial z'}}\int _{z'=0}^{\infty }dSA\left(\mathbf {0} ,t-{\frac {S}{c}}\right)=A\left(\mathbf {\emptyset } ,t'-{\frac {z'}{c}}\right)}
This gives us uniform translation of waves at velocity c. More generally:
A
(
r
′
)
=
−
1
2
π
∫
z
=
0
d
2
r
∂
∂
z
′
(
A
(
r
,
t
′
−
|
r
−
r
′
|
c
)
|
r
−
r
′
|
)
{\displaystyle A(r')={\frac {-1}{2\pi }}\int _{z=0}d^{2}r{\frac {\partial }{\partial z'}}\left({\frac {A\left(\mathbf {r} ,t'-{\frac {|\mathbf {r} -\mathbf {r} '|}{c}}\right)}{|\mathbf {r} -\mathbf {r} '|}}\right)}
=
−
1
2
π
∫
z
=
0
d
2
r
(
A
(
r
,
t
′
−
|
r
−
r
′
|
c
)
|
r
−
r
′
|
3
(
−
z
′
)
+
A
˙
(
r
,
t
′
−
|
r
−
r
′
|
c
)
|
r
−
r
′
|
c
−
z
′
|
r
−
r
′
|
)
{\displaystyle ={\frac {-1}{2\pi }}\int _{z=0}d^{2}r\left({\frac {A\left(\mathbf {r} ,t'-{\frac {|\mathbf {r} -\mathbf {r} '|}{c}}\right)}{|\mathbf {r} -\mathbf {r} '|^{3}}}(-z')+{\frac {{\dot {A}}\left(\mathbf {r} ,t'-{\frac {|\mathbf {r} -\mathbf {r} '|}{c}}\right)}{|\mathbf {r} -\mathbf {r} '|c}}{\frac {-z'}{|\mathbf {r} -\mathbf {r} '|}}\right)}
A
(
r
′
)
=
−
1
2
π
∫
z
=
0
d
2
r
(
A
(
r
,
t
′
−
|
r
−
r
′
|
c
)
|
r
−
r
′
|
3
(
−
z
′
)
+
1
c
A
˙
(
r
,
t
′
−
|
r
−
r
′
|
c
)
|
r
−
r
′
|
2
(
z
′
)
)
{\displaystyle A(r')={\frac {-1}{2\pi }}\int _{z=0}d^{2}r\left({\frac {A\left(\mathbf {r} ,t'-{\frac {|\mathbf {r} -\mathbf {r} '|}{c}}\right)}{|\mathbf {r} -\mathbf {r} '|^{3}}}(-z')+{\frac {1}{c}}{\frac {{\dot {A}}\left(\mathbf {r} ,t'-{\frac {|\mathbf {r} -\mathbf {r} '|}{c}}\right)}{|\mathbf {r} -\mathbf {r} '|^{2}}}(z')\right)}
In our case, we consider only those waves which degrade as
1
r
{\displaystyle {\frac {1}{r}}\quad }
, so:
A
(
r
′
)
=
−
1
2
π
∫
z
=
0
d
2
r
(
1
c
A
˙
(
r
,
t
′
−
|
r
−
r
′
|
c
)
|
r
−
r
′
|
2
(
z
′
)
)
{\displaystyle A(r')={\frac {-1}{2\pi }}\int _{z=0}d^{2}r\left({\frac {1}{c}}{\frac {{\dot {A}}\left(\mathbf {r} ,t'-{\frac {|\mathbf {r} -\mathbf {r} '|}{c}}\right)}{|\mathbf {r} -\mathbf {r} '|^{2}}}(z')\right)}
A
(
r
′
)
=
−
z
′
2
π
c
∫
z
=
0
d
2
r
(
A
˙
(
r
,
t
′
−
|
r
−
r
′
|
c
)
|
r
−
r
′
|
2
)
{\displaystyle A(r')={\frac {-z'}{2\pi c}}\int _{z=0}d^{2}r\left({\frac {{\dot {A}}\left(\mathbf {r} ,t'-{\frac {|\mathbf {r} -\mathbf {r} '|}{c}}\right)}{|\mathbf {r} -\mathbf {r} '|^{2}}}\right)}
In cylindrical coordinates,
d
2
r
=
r
d
r
d
ϕ
{\displaystyle d^{2}r=rdrd\phi \quad }
. Also,
A
˙
(
r
,
t
′
−
|
r
−
r
′
|
c
)
=
A
0
˙
e
−
i
k
(
t
′
c
−
|
r
−
r
′
|
)
{\displaystyle {\dot {A}}\left(\mathbf {r} ,t'-{\frac {|\mathbf {r} -\mathbf {r} '|}{c}}\right)={\dot {A_{0}}}e^{-ik(t'c-|\mathbf {r} -\mathbf {r} '|)}}
. So:
A
(
r
′
)
=
−
z
′
A
0
˙
2
π
c
∫
z
=
0
r
d
r
d
ϕ
e
−
i
k
(
t
′
c
−
|
r
−
r
′
|
)
|
r
−
r
′
|
2
{\displaystyle A(r')={\frac {-z'{\dot {A_{0}}}}{2\pi c}}\int _{z=0}rdrd\phi {\frac {e^{-ik(t'c-|\mathbf {r} -\mathbf {r} '|)}}{|\mathbf {r} -\mathbf {r} '|^{2}}}}
Let
J
(
r
′
)
=
∫
0
a
r
d
r
∫
0
2
π
d
ϕ
e
i
k
|
r
−
r
′
|
|
r
−
r
′
|
{\displaystyle {\mathcal {J}}(r')=\int _{0}^{a}rdr\int _{0}^{2}\pi d\phi {\frac {e^{ik|\mathbf {r} -\mathbf {r} '|}}{|\mathbf {r} -\mathbf {r} '|}}}
Then
A
(
r
′
)
=
−
z
′
A
0
˙
2
π
c
e
−
i
k
c
t
′
J
(
r
′
)
{\displaystyle A(r')={\frac {-z'{\dot {A_{0}}}}{2\pi c}}e^{-ikct'}{\mathcal {J}}(r')}