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 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(q)=\frac{1}{(2\pi)^2} \int d^4qe^{iqr} \tilde{G}(r,0)}
In free space, translational symmetry implies:
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',0)=G(r,r') \quad }
∴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)^2}\int d^4q e^{-iq(r-r')} \tilde{G} (q)}
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 \square^2_rG(r,r')=\frac{1}{(2\pi)^2}|int d^4qe^{-iq(r-r')}\tilde{G}(q)}
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 \square^2_rG(r,r')=\frac{1}{(2\pi)^2}\int d^4qe^{-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')}}
∴
G
~
(
q
)
=
(
2
π
)
2
(
2
π
)
4
1
−
q
2
=
−
1
(
2
π
)
2
q
2
{\displaystyle {\tilde {G}}(q)={\frac {(2\pi )^{2}}{(2\pi )^{4}}}{\frac {1}{-q^{2}}}={\frac {-1}{(2\pi )^{2}q^{2}}}}
G
(
r
,
r
′
)
=
−
1
(
2
π
)
4
∫
d
4
q
e
−
i
q
(
r
−
r
′
)
1
(
k
+
ω
c
)
(
k
−
ω
c
)
{\displaystyle G(r,r')={\frac {-1}{(2\pi )^{4}}}\int d^{4}qe^{-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'
G
(
r
,
r
′
)
=
1
(
2
π
)
4
∫
d
3
k
e
−
i
k
(
r
−
r
′
)
∫
d
(
ω
c
)
e
i
ω
(
t
−
t
′
)
(
ω
c
−
k
)
(
ω
c
+
k
)
Θ
{\displaystyle G(r,r')={\frac {1}{(2\pi )^{4}}}\int d^{3}ke^{-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, 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 d^2r=rdrd\phi \quad}
. Also, 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 \dot{A}\left(\mathbf{r},t'-\frac{|\mathbf{r}-\mathbf{r}'|}{c}\right)=\dot{A_0}\frac{e^{-ik(t'c-|\mathbf{r}-\mathbf{r}'|)}}{|\mathbf{r}-\mathbf{r}'}}