Lect4 Optical Waveguides
Lect4 Optical Waveguides
Lect4 Optical Waveguides
Waveguide structures
Waveguide modes
Field equations
Wave equations
Guided modes in symmetric slab waveguides
General formalisms for step-index planar
waveguides
1
Waveguide structures
2
Waveguide structures
• The basic structure of a dielectric waveguide consists of a
longitudinally extended high-index optical medium, called the core,
which is transversely surrounded by low-index media, called the
cladding. A guided optical wave propagates in the waveguide along its
longitudinal direction.
x x
n1 n(x)
y z y z
n2 n2
n(x) n(x)
6
Channel waveguides
Most waveguides used in device applications are nonplanar
waveguides.
For a nonplanar waveguide, the index profile n(x, y) is a
function of both transverse coordinates x and y.
There are many different types of nonplanar waveguides that
are differentiated by the distinctive features of their index
profiles.
One very unique group is the circular optical fibers (to be
discussed in Lecture 5).
Another important group of nonplanar waveguides is the
channel waveguides, which include
The buried channel waveguides
The strip-loaded waveguides
The ridge waveguides
The rib waveguides
The diffused waveguides.
7
Representative channel waveguides
w
n3 w
n2 d n1 d n1
d n1 n2 n2
w
Buried channel Strip-loaded Ridge
waveguide waveguide waveguide
w x
h
d n1
n1 y
d
n2 n2
w
z
rib waveguide Diffused waveguide 8
Representative channel waveguides
A buried channel waveguide is formed with a high-index
waveguiding core buried in a low-index surrounding medium.
The waveguiding core can have any cross-sectional geometry
though it is often a rectangular shape.
A strip-loaded waveguide is formed by loading a planar
waveguide, which already provides optical confinement in
the x direction, with a dielectric strip of index n3<n1 or a
metal strip to facilitate optical confinement in the y direction.
The waveguiding core of a strip waveguide is the n1 region
under the loading strip, with its thickness d determined by the
thickness of the n1 layer and its width w defined by the width
of the loading strip.
A ridge waveguide has a structure that looks like a strip
waveguide, but the strip, or the ridge, on top of its planar
structure has a high index and is actually the waveguiding
core. A ridge waveguide has strong optical confinement
because it is surrounded on three sides by low-index air (or 9
cladding material).
Representative channel waveguides
A rib waveguide has a structure similar to that of a strip or
ridge waveguide, but the strip has the same index as the high-
index planar layer beneath it and is part of the waveguiding
core.
These four types of waveguides are usually termed
rectangular waveguides with a thickness d in the x direction
and a width w in the y direction, though their shapes are
normally not exactly rectangular.
A diffused waveguide is formed by creating a high-index
region in a substrate through diffusion of dopants, such as
LiNbO3 waveguide with a core formed by Ti diffusion.
Because of the diffusion process, the core boundaries in the
substrate are not sharply defined.
A diffused waveguide also has a thickness d defined by the
diffusion depth of the dopant in the x direction and a width w
defined by the distribution of the dopant in the y direction. 10
Silicon optical waveguides (nanophotonic wires)
air w
Si h d
SiO2
Si substrate
silica silica
13
Waveguide modes
14
Waveguide modes
Waveguide modes exist that are characteristic of a particular
waveguide structure.
A waveguide mode is a transverse field pattern whose
amplitude and polarization profiles remain constant along
the longitudinal z coordinate.
Therefore, the electric and magnetic fields of a mode can be
written as follows
17
Intuitive ray and wavefront picture
wavefronts
x
n3
λ/n1
θ θ
n1 d z
n2
There are two critical angles associated with the internal
reflections at the lower and upper interfaces:
If θ > θc2 > θc3, the wave inside the core is totally reflected at
both interfaces and is trapped by the core, resulting in guided 18
modes.
Guided modes
As the wave is reflected back and forth between the two
interfaces, it interferes with itself.
19
k-vector triangle
• The orthogonal components of the propagation constant, β and
kx, are related by the “k-vector triangle.”
x
k1 θ
kx
core n1
β
z
d k1
z
β
kx
d + d
kx
where m is an integer = 0, 1, 2, …
Because m can assume only integral values, only certain
discrete values of θ can satisfy the transverse resonance 22
condition.
Transverse electric polarization
• For slab waveguides, we define the x-z plane as the plane of incidence.
• Waves with this polarization are labeled transverse electric (TE) fields
because the electric field vector lies entirely in the x y plane (i.e. Ez = 0)
that is transverse to the direction of net travel (the z direction).
n1k x
E
H z
core n1 TE
cladding
23
Transverse magnetic polarization
• For the parallel, or p, polarization, the electric field is no longer
purely transverse. It has a component along the z direction.
n1k x
E
H z
TM
core n1
cladding
24
Discrete guided modes
The transverse resonance condition results in discrete values
of the propagation constant βm for guided modes identified by
the mode number m.
Although the critical angles, θc2 and θc3, do not depend on the
polarization of the wave, the phase shifts, ϕ2(θ) and ϕ3(θ),
caused by the internal reflection at a given angle θ depend on
the polarization.
Therefore, TE and TM waves have different solutions for
the transverse resonance condition, resulting in different βm
and different mode characteristics for a given mode number
m.
For a given polarization, solution of the transverse resonance
condition yields a smaller value of θ and a correspondingly
smaller value of β for a larger value of m. Therefore, β0 > β1
> β2 > …
The guided mode with m = 0 is called the fundamental mode
and those with m ≠ 0 are higher-order modes. 25
Qualitative picture of a waveguide mode
Standing wave oscillation (in the x direction)
cladding x
core n1 + - kx + - z
θ1 k1
k1 h x
d n1 β
z
n2
d z
28
Discrete waveguide modes
• Because m can assume only integral values, only certain discrete
values of θ = θm can satisfy the resonance condition.
β1
• As the spatial rate of change of phase at the boundary (or the
projection of the wavefront propagation) on the n1 side must match with
that on the n2 side, we have β1 = β2 = β. This condition is known as
“phase-matching.” Phase-matching allows coupling of oscillating field
between two media.
β1 = β2 => n1k sin θ1 = n2k sin θ2 (Snell’s Law) 31
Evanescent field in total internal reflection
• Phase-matching at TIR θ1 > θc (i.e. sin θ1 > n2/n1)
β1 = n1 k sin θ1 = β2 > n2k
• k-vector triangle in n2
k2x = [(n2k)2 – (n1 k sin θ1)2]1/2
where E2 is the peak value of the electric field at the lower (x = -d/2) and
upper (x = d/2) boundaries. 33
Waveguide effective index
• We can define the waveguide phase velocity vp as
vp = ω / β
neff = c / vp
Or neff = cβ / ω = β / k
=> neff = n1 sin θ
• The effective refractive index is a key parameter in guided
propagation, just as the refractive index is in unguided wave travel.
34
• For waveguiding at n1-n2 interface, we see that n2 ≤ neff ≤ n1
λz = λ/neff 35
Field equations
36
Field equations
For a linear, isotropic dielectric waveguide characterized by a
spatial permittivity distribution ε(x, y), Maxwell’s equations
can be written as
37
Field equations
These two Maxwell’s equations can be written as
38
Field equations
k2 = ω2µ0ε(x, y) is a function of x and y to account for the
transverse spatial inhomogeneity of the waveguide structure.
41
Wave equations
The common approach to finding Ez and Hz is to solve the
wave equations together with boundary conditions.
For the case of a linear, isotropic waveguide with a spatially
dependent ε(x, y), the two Gauss’ laws for E and H can be
written as
42
Wave equations
The three components Ex, Ey and Ez for the electric field are
generally coupled together because ∇ε ≠ 0 in a waveguide.
For the same reason Hx, Hy and Hz are also coupled.
This fact indicates that the vectorial characteristics of a mode
field in a waveguide are strongly dependent on the geometry
and index profile of the waveguide.
43
Wave equations for step-index waveguides
The index profile of a step-index waveguide is piecewise
constant.
We can write a homogeneous wave equation separately for
each region of constant ε because ∇ε = 0 within each region.
Assuming E and H of the harmonic guided wave form, and
with ∇ε = 0 for each region of constant ε, we obtain for the
longitudinal components
However, note that this approach does not work for graded-
index waveguides because ∇ε ≠ 0 for such waveguides.
45
Wave equations for planar waveguides
Homogeneous wave equations exist for planar waveguides of
any index profile n(x).
For a planar waveguide, the modes are either TE or TM.
Furthermore, we consider ∂/∂y = 0 because the index profile
is independent of the y coordinate. The wave equations are
substantially simplified.
x
n3
y n1
z (x, y, z)
n2 (n1 > n2, n3)
Planar (slab) waveguides for
integrated photonics (e.g. laser chips) 46
Wave equations for planar waveguides: TE modes
For any TE mode of a planar waveguide, Ez = 0.
Such that from relations of the transverse components to the
longitudinal components, we see that Ex = Hy = 0 because ∂Hz/
∂y = 0.
The only nonvanishing field components are Hx, Ey and Hz.
Because there is only one nonvanishing electric field
component Ey, the wave equation for Ey is naturally decoupled
from the other field components. Therefore, we have
47
Wave equations for planar waveguides: TM modes
For any TM mode of a planar waveguide, Hz = 0.
Then, Hx = Ey = 0 because ∂Ez/∂y = 0.
The only nonvanishing field components are Ex, Hy and Ez.
The wave equation for Hy
48
Guided modes in symmetric
slab waveguides
49
The solutions
The solutions to the wave equation for the three regions have
the form (assuming TE polarization)
50
Symmetric slab waveguide solutions
The solutions are sinusoidal or exponential according to
Sinusoidal ki2 > β2
Exponential ki2 < β2
Guided modes: cladding region has exponential solutions β
> k2 > k3 while core region has sinusoidal solutions β < k1
For x → ±∞, we require the solutions to remain finite.
The solutions have the general form (assuming n2 = n3):
Cladding (x ≥ 0)
Core (-d ≤ x ≤ 0)
Cladding (x ≤ -d)
51
Symmetric slab waveguide solutions
The transverse propagation constants k2
iκ
κ = (β2 – k22)1/2 β
x=0
k1
θ
h
h = (k12 – β2)1/2 = k1cosθ
β
x=-d
In order to determine the allowed β and the unspecified
constants A, B, C and D, we need to match the solution in
cladding with the solution in core.
Therefore, boundary conditions must be specified at the
core-cladding interfaces.
We expect at least one arbitrary constant in the final solution
given by the overall field strength.
52
Boundary conditions
Ey continuous at x = 0 → A = B
The first three results can be used to solve the electric field
distributions within the waveguide core and cladding.
x ≥ 0
-d ≤ x ≤ 0
x ≤ -d
55
Eigenvalue equations for symmetric slab waveguides
Recall the reflection coefficient rTE = exp –iϕTE=> ϕ(θ) = -ϕTE(θ)
rTM = -exp –iϕTM=> ϕ(θ) = -ϕTM(θ)
TE:
tan [(k1d cos θ)/2 - mπ/2] = (n12 sin2θ - n22)1/2 / (n1 cos θ)
TM:
tan [(k1d cos θ)/2 - mπ/2] = (n12/n22) (n12 sin2θ - n22)1/2 / (n1 cos θ)
m = 0, 1, 2, …, 56
Normalized waveguide parameters
The mode properties of a waveguide are commonly characterized
in terms of dimensionless normalized waveguide parameters.
The normalized frequency, also known as the V number, of a
step-index planar waveguide is defined as
58
Example #1: Symmetric weakly guiding slab waveguides
59
Graphic solutions for the eigenvalues of guided TE and TM
modes of a weakly guiding symmetric slab waveguide
n1 = 3.6, n2 = 3.55, V = 5π
tan (hd/2 - mπ/2)
V=5π
TM
m=0 1 2 3 4 5
TE
hd
60
Mode chart for the first six TE and TM modes (m = 0 – 5) of
symmetric slab waveguides in AlGaAs (n1 = 3.6, n2 = 3.55)
(neff = n1)
m=0
Normalized guide index b
1
2
TE
3
TM
4 5
(neff = n2)
(V = π)
V V = 5π
Cut-off for m = 1 61
Example #2: Symmetric strongly guiding slab waveguides
62
Graphic solutions for the eigenvalues of guided TE and TM
modes of a strongly guiding symmetric slab waveguide
n1 = 3.5, n2 = 1.45, V = 5π
V=5π
tan (hd/2 - mπ/2)
TM
m=0 1 2 3 4 5
TE
hd
63
Mode chart for the first six TE and TM modes (m = 0 – 5) of
symmetric slab waveguides in SOI (n1 = 3.5, n2 = 1.45)
(neff = n1)
m=0
Normalized guide index b
TE 2
3
TM
4 5
(neff = n2)
(V = π)
V V = 5π
Cut-off for m = 1 64
Normalized guide index vs. V number
• When the V number is very small (e.g. d/λ << 1), the guided ray
travels close to the critical angle (b << 1) . The effective index is close
to that of the cladding layer n2.
=>The wave penetrates deeply into the cladding layers, because the rays
are near the critical angle. The evanescent decay is slow.
• As the V number increases, the ray travels more nearly parallel to the
waveguide axis, and the effective refractive index lies between n1 and n2.
• For a very large V number (e.g. d/λ >> 1) the effective index is near
that of the core index n1. The wave in the cladding layer decays very
rapidly for evanescent waves traveling at angles far above the critical
angle. 65
Cutoff conditions
• For example, consider V = 15 on the mode chart, the TE5/TM5 modes
could not propagate because V was not large enough to intersect with the
b vs. V curves.
=> The TE5/TM5 modes, and all higher-ordered modes, are cut off.
• Cutoff occurs when the propagation angle for a given mode just equals
the critical angle θc --- a guided mode transits to an unguided radiation
mode.
• The fields extend to infinity for κ = 0 (i.e. the fields become unguided!).
This defines the cutoff condition for guided modes.
66
Cutoff conditions
The cutoff value V = Vc for a particular guided mode is the
value of V at the point where b = 0 (i.e. the b vs. V relation
intersects with the axis b = 0).
For κ = 0,
⇒ Vc = hd
=> Vc = mπ, m = 0, 1, 2, …
67
Number of modes
TE and TM modes of a symmetric waveguide have the same
cutoff condition:
Vc = mπ
for the mth TE and TM modes
0
1
TE 2
TM
3
4
5
69
V
Singlemode vs. multimode waveguides
• A multimode waveguide is one that supports more than one
propagating mode.
• For a multimode waveguide, at a fixed thickness the higher-ordered
modes propagate with smaller β values than the lower-ordered modes.
Higher-order
d lower-order
*If we wish to propagate only the TE0/TM0 mode, then we must have
V < π
single-mode condition
=> This cuts off the m = 1 mode and all higher-order modes. 70
Further remarks on TE/TM mode charts
• For the weakly guiding AlGaAs waveguide, TE and TM modes have
about the same effective index and propagation angle, but their electric-
field vectors point in orthogonal directions.
73
Step-index planar waveguides
A step-index planar (slab) waveguide
x
n3 k3
d/2
β
n1 k1
-d/2
n2 k2
z
(assume n1 > n2 > n3)
74
Normalized waveguide parameters
The measure of the asymmetry of the waveguide is represented
by an asymmetry factor a, which depends on the polarization of
the mode. For TE modes, we have
77
Longitudinal propagation constant β
Because k1, k2 and k3 are well-defined parameters of a given
waveguide, the only parameter that has to be determined for a
particular waveguide mode is the longitudinal propagation
constant β.
x < -d/2 79
Guided TE modes
And the following eigenvalue equations:
x < -d/2
and
82
Modal dispersion
Guided modes have discrete allowed values of β.
They are determined by the allowed values of h1 because β
and h1 are directly related to each other through k12 – β2 = h12.
Because κ2 and κ3 are uniquely determined by β, they are also
uniquely determined by h1.
In terms of the normalized waveguide parameters, we have
83
Modal dispersion
The solutions of the eigenvalue equation for TE modes yield
the allowed parameters for guided TE modes, while those of
the eigenvalue equation for TM modes yield the allowed
parameters for guided TM modes.
A transcendental equation is usually solved graphically by
plotting its left- and right-hand sides as a function of h1d.
The solutions yield the allowed values of β, or the normalized
guide index b, as a function of the parameters a and V.
For a given waveguide, a guided TE mode has a larger
propagation constant than the corresponding TM mode of the
same order:
87
These equations yield only discrete eigenvalues for given
values of waveguide parameters n1, n2, n3, V and aE.
They are transcendental equations that have to be solved
graphically or numerically.
88
Cutoff conditions
κ2 and κ3 are real and positive for a guided mode, s.t. the
fields of the mode decay exponentially in the transverse
direction outside the core region and remain bound to the
core.
This is equivalent to the condition that θ > θc2 > θc3 in the ray
optics picture s.t. the ray in the core is totally reflected by
both interfaces.
Because θc2 > θc3, the transition from a guided mode to an
unguided radiation mode occurs when θ = θc2.
This corresponds to the condition that β = k2 and κ2 = 0.
The fields extend to infinity on the substrate side for κ2 = 0.
This defines the cutoff condition for guided modes.
89
Cutoff conditions
At cutoff, V = Vc. The cutoff value Vc for a particular guided
mode is the value of V at the point where the curve of its b vs.
V dispersion relation intersects with the horizontal axis b = 0.
By setting κ2 = 0, we find from the normalized waveguide
parameters relations that
and
m = 0, 1, 2, …
90
Cutoff conditions
A similar mathematical procedure yields the cutoff condition
for the mth guided TM mode:
m = 0, 1, 2, …
93