An example of such surface of discontinuity is the initial wave front emanating from a source that starts radiating at a certain instant of time.
The surfaces of field discontinuity thus become geometrical optics wave fronts with the corresponding geometrical optics fields defined as:
Those fields obey transport equations consistent with the transport equations of the Sommerfeld-Runge approach. Light rays in Luneburg's theory are defined as trajectories orthogonal to the discontinuity surfaces and with the right parametrisation they can be shown to obey Fermat's principle of least time thus establishing the identity of those rays with light rays of standard optics.
The above developments can be generalised to anisotropic media.[11]
The proof of Luneburg's theorem is based on investigating how Maxwell's equations govern the propagation of discontinuities of solutions. The basic technical lemma is as follows:
A technical lemma. Let
be a hypersurface (a 3-dimensional manifold) in spacetime
on which one or more of:
,
-
∮
Γ
0
(
∇
φ
⋅
[
ε
E
→
]
)
d
S
‖
∇
4
D
φ
‖
=
0
{\displaystyle \oint _{\Gamma _{0}}(\nabla \varphi \cdot [\varepsilon \mathbf {\vec {E}} ])\,{dS \over \|\nabla ^{4D}\varphi \|}=0}
![{\displaystyle \oint _{\Gamma _{0}}(\nabla \varphi \cdot [\varepsilon \mathbf {\vec {E}} ])\,{dS \over \|\nabla ^{4D}\varphi \|}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/765fd0c151f60c43f357280e247d260dc0581ee6)
-
∮
Γ
0
(
∇
φ
⋅
[
μ
H
→
]
)
d
S
‖
∇
4
D
φ
‖
=
0
{\displaystyle \oint _{\Gamma _{0}}(\nabla \varphi \cdot [\mu \mathbf {\vec {H}} ])\,{dS \over \|\nabla ^{4D}\varphi \|}=0}
![{\displaystyle \oint _{\Gamma _{0}}(\nabla \varphi \cdot [\mu \mathbf {\vec {H}} ])\,{dS \over \|\nabla ^{4D}\varphi \|}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/13ce0be2f64e43744b3df2b6db3883789c3a4846)
-
∮
Γ
0
(
∇
φ
×
[
H
→
]
−
1
c
φ
t
[
ε
∇
4
D
{\displaystyle \nabla ^{4D}}
denotes the gradient in the 4D
x
y
z
t
{\displaystyle xyzt}
-space. And since
Γ
0
{\displaystyle \Gamma _{0}}
is arbitrary, the integrands must be equal to 0 which proves the lemma.
It's now easy to show that as they propagate through a continuous medium, the discontinuity surfaces obey the eikonal equation. Specifically, if
ε
{\displaystyle \varepsilon }
and
μ
{\displaystyle \mu }
are continuous, then the discontinuities of
E
→
{\displaystyle \mathbf {\vec {E}} }
and
H
→
{\displaystyle \mathbf {\vec {H}} }
satisfy:
[
ε
E
→
]
=
ε
[
E
→
]
{\displaystyle [\varepsilon \mathbf {\vec {E}} ]=\varepsilon [\mathbf {\vec {E}} ]}
and
[
μ
H
→
]
=
μ
[
H
→
]
{\displaystyle [\mu \mathbf {\vec {H}} ]=\mu [\mathbf {\vec {H}} ]}
. In this case the first two equations of the lemma can be written as:
-
∇
φ
×
[
H
→
]
−
ε
c
φ
t
[
E
→
]
=
0
{\displaystyle \nabla \varphi \times [\mathbf {\vec {H}} ]-{\varepsilon \over c}\,\varphi _{t}\,[\mathbf {\vec {E}} ]=0}
![{\displaystyle \nabla \varphi \times [\mathbf {\vec {H}} ]-{\varepsilon \over c}\,\varphi _{t}\,[\mathbf {\vec {E}} ]=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/73dddb90e4cf5597724bc161094ddbe008d97d02)
-
∇
φ
×
[
E
→
]
+
μ
c
φ
t
[
H
→
]
=
0
{\displaystyle \nabla \varphi \times [\mathbf {\vec {E}} ]+{\mu \over c}\,\varphi _{t}\,[\mathbf {\vec {H}} ]=0}
![{\displaystyle \nabla \varphi \times [\mathbf {\vec {E}} ]+{\mu \over c}\,\varphi _{t}\,[\mathbf {\vec {H}} ]=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f3789a6ee64532429305477a7a4ee4af8ce77655)
Taking cross product of the first equation with
∇
φ
{\displaystyle \nabla \varphi }
and substituting the second yields:
-
ε
{\displaystyle \varepsilon }
and
μ
{\displaystyle \mu }
are continuous, then the discontinuities of
E
→
{\displaystyle \mathbf {\vec {E}} }
and
H
→
{\displaystyle \mathbf {\vec {H}} }
satisfy:
[
ε
E
→
]
=
ε
[
E
→
]
{\displaystyle [\varepsilon \mathbf {\vec {E}} ]=\varepsilon [\mathbf {\vec {E}} ]}
and
[
μ
H
→
]
=
μ
[
H
→
]
{\displaystyle [\mu \mathbf {\vec {H}} ]=\mu [\mathbf {\vec {H}} ]}
. In this case the first two equations of the lemma can be written as:
Taking cross product of the first equation with
∇
φ
{\displaystyle \nabla \varphi }
and substituting the second yields:
-
∇
φ
×
(
∇
φ
×
[
H
→
]
)
−
ε
c
φ
t
(
∇
φ
×
[
E
→
]
)
=
(
∇
φ
⋅
[
H
→
]
)
∇
φ
−
‖
∇
φ
‖
2
[
H
→
]
+
ε
μ
c
2
φ
t
2
[
By the second of Maxwell's equations,
∇
φ
⋅
[
H
→
]
=
0
{\displaystyle \nabla \varphi \cdot [\mathbf {\vec {H}} ]=0}
, hence, for points lying on the surface
φ
=
0
{\displaystyle \varphi =0}
only:
-
‖
∇
φ
‖
2
=
ε
μ
c
2
φ
t
2
{\displaystyle \|\nabla \varphi \|^{2}={\varepsilon \mu \over c^{2}}\varphi _{t}^{2}}

(Notice the presence of the discontinuity is essential in this step as we'd be dividing by zero otherwise.)
Because of the physical considerations one can assume without loss of generality that
φ
{\displaystyle \varphi }
is of the following form:
φ
(
x
,
y
,
z
,
t
)
=
ψ
(
x
,
y
,
z
)
−
c
t
{\displaystyle \varphi (x,y,z,t)=\psi (x,y,z)-ct}
, i.e. a 2D surface moving through space, modelled as level surfaces of
ψ
{\displaystyle \psi }
. (Mathematically
ψ
{\displaystyle \psi }
exists if (Notice the presence of the discontinuity is essential in this step as we'd be dividing by zero otherwise.)
Because of the physical considerations one can assume without loss of generality that
φ
{\displaystyle \varphi }
is of the following form:
φ
(
x
,
y
,
z
,
t
)
=
ψ
(
x
,
y
,
z
)
−
c
t
{\displaystyle \varphi (x,y,z,t)=\psi (x,y,z)-ct}
φ
{\displaystyle \varphi }
is of the following form:
φ
(
x
,
y
,
z
,
t
)
=
ψ
(
x
,
y
,
z
)
−
c
t
{\displaystyle \varphi (x,y,z,t)=\psi (x,y,z)-ct}
, i.e. a 2D surface moving through space, modelled as level surfaces of
ψ
{\displaystyle \psi }
. (Mathematically
ψ
{\displaystyle \psi }
exists if
φ
t
≠
0
{\displaystyle \varphi _{t}\neq 0}
by the implicit function theorem.)
The above equation written in terms of
ψ
{\displaystyle \psi }
becomes:
i.e.,
-
ψ
x
2
+
ψ
y
2
+
ψ
z
2
=
n
2
{\displaystyle \psi _{x}^{2}+\psi _{y}^{2}+\psi _{z}^{2}=n^{2}}

which is the eikonal equation and it holds for all
x
{\displaystyle x}
which is the eikonal equation and it holds for all
x
{\displaystyle x}
,
y
{\displaystyle y}
,
z
{\displaystyle z}
, since the variable
t
{\displaystyle t}
is absent. Other laws of optics like Snell's law and Fresnel formulae can be similarly obtained by considering discontinuities in
ε
{\displaystyle \varepsilon }
and
μ
{\displaystyle \mu }
.
General equation using four-vector notation
In four-vector notation used in special relativity, the wave equation can be written as
-
∂
2
ψ
∂
x
i
∂
x
In four-vector notation used in special relativity, the wave equation can be written as
-
and the substitution
ψ
=
A
e
i
S
/
ϵ
{\displaystyle \psi =Ae^{iS/\epsilon }}
leads to[12]
-
−
A
ϵ
2
∂
S
∂
x
i
∂
S
∂
x
∂
S
∂
x
i
∂
S
∂
x
i
=
0.
{\displaystyle {\frac {\partial S}{\partial x_{i}}}{\frac {\partial S}{\partial x^{i}}}=0.}

Once eikonal is found by solving the above equation, the wave four-vector can be found from
-
k
i
=
−
∂
S
∂
x
i
.
{\displaystyle k_{i}=-{\frac {\partial S}{\partial x^{i}}}.}

See also
References
- ^ Arthur Schuster, An Introduction to the Theory of Optics, London: Edward Arnold, 1904 online.
- ^ Greivenkamp, John E. (2004). Field Guide to Geometrical Optics. SPIE Field Guides. 1. SPIE. pp. 19–20. ISBN 0-8194-5294-7.