Lecture Slab Waveguide Analysis
Lecture Slab Waveguide Analysis
Lecture Slab Waveguide Analysis
Computational Science:
Computational Methods in Engineering
Outline
• Slab Waveguides
• Formulation
• Solution
• Implementation in MATLAB
• More About Resolution and Spacer Regions
1
2/18/2020
Slab Waveguides
Slide 3
Refractive Index n
Light travels at different speeds when it is inside different materials.
Frequency is constant.
Speed changes.
Wavelength changes.
2
2/18/2020
Snell’s Law
z
Snell’s law quantifies the angles
of light rays at an interface.
1 n1
x
n1 sin 1 n2 sin 2
n2
2
Slide 5
Critical Angle c
z There exists a special angle, the
critical angle, where the ray in
the low-index medium is at 90°.
1 90 n1 n1 sin 1 n2 sin 2
n1 sin 90 n2 sin c
x n1 n2 sin c
sin c n1 n2
c sin 1 n1 n2
n2 2 c
c sin 1 n1 n2
where n2 n1
Slide 6
3
2/18/2020
Slide 7
n1 Conditions
n2 n1
TIR
and
n2 n2 n3
TIR
n3
Slide 8
4
2/18/2020
k0 neff k0 n sin
Slide 9
Rigorous Analysis
k0 neff k0 n sin
Slide 10
5
2/18/2020
Confinement
Confinement
Slide 11
Wave oscillations
y z phase constant
Slide 12
6
2/18/2020
Formulation
Slide 13
What is Formulation?
Formulation is the initial analytical work we do before implementing a computer code.
Usually we start with the governing equation(s) and end with the matrix equation to be
solved.
Slide 14
7
2/18/2020
Governing Equations
Since this is an electrodynamics problem, we start with Maxwell’s curl equations.
E j H
H j E
Vector Curl
The curl of a vector is a measure of the vector field’s tendency to circulate about an axis. The curl quantity is
directly along this axis and the magnitude measures the strength of the circulation.
Slide 15
8
2/18/2020
H z H y
H j E y
z
j Ex
H x H z
j E y
z x
H y H x
j Ez
x y
Slide 17
0
y
y z
z
Device is uniform. Wave
propagates in this direction so j
wave phase is increasing. z
Slide 18
9
2/18/2020
y 0 Apply
Since nothing is changing in the y direction, any derivative with respect to y must be zero.
Ez E y E y
j H x j H x
y z z
Ex Ez Ex Ez
j H y j H y
z x z x
E y Ex E y
j H z j H z
x y x
H z H y H y
j E x j E x
y z z
H x H z H x H z
j E y j E y
z x z x
H y H x H y
j Ez j E z
x y x
Slide 19
10
2/18/2020
What About ?
The guided mode has the following mathematical form
E x, y , z A x e j z
Now calculate the partial derivative with respect to z and see what happens.
E x, y, z A x e j z A x e j z e j z A x
z z z z
j A x e j z j E x, y , z
E x, y , z
1D Governing Equations
The equations for the Ey mode were
H z
Hx j E y
z x
E y j H x
z
E y
j H z
x
Slide 22
11
2/18/2020
x k0 x 2
k0
0 n
Slide 23
12
2/18/2020
Normalized Equations
Applying the normalizations to all equations to get
dH z
jneff H x j 0 r Ey
dx 0
0
jneff E y j H
0 r x
dE y
j 0 r H z
dx 0
These are substituted into the first equation to get a single equation containing only Ey.
This is why it was called the Ey mode.
dH z
jneff H x j 0 r Ey
dx 0
0 d 0 dE y 0
jneff neff E y j j r Ey
0
dx 0 dx 0
d 2Ey d 2 Ey
2
neff Ey r Ey r E y neff
2
Ey
dx 2 dx 2
Slide 26
13
2/18/2020
Eigen-Value Problem
For optical problems, people like to put everything in terms of refractive index. This is
related to the relative permittivity through r = n2.
d 2 Ey
n 2 E y neff
2
Ey
dx 2
Slide 27
Matrix Form
We go term-by-term to write the equation in matrix form.
d2
dx 2 n x E y x neff E y x
2 2
D 2
x n 2 e y neff
2
ey
or
D 2
x ε e y neff
2
ey
Slide 28
14
2/18/2020
Solution
Slide 29
e1 2
2 M
y ey 2
2
ey 2
M
e1 3
ey 3 e 3 M = # modes
V y y Usually M = Nx
1
ey N x 1 ey N x 1 eyM N x 1
2
e1 N 2 M
y ey N x ey N x
x
n
2
1
eff
n Eigen-vectors and eigen-values come in pairs.
2
2
D x eff
Do not mix up their pairing!
n
2
M
eff
Slide 30
15
2/18/2020
Slide 31
Implementation
in MATLAB
Slide 32
16
2/18/2020
Implementation Outline
1. Initialize MATLAB
2. Dashboard (materials, dimensions, etc.)
3. Calculate Grid
4. Build Device on Grid
5. Perform Finite-Difference Analysis
6. Visualize the Results
Slide 33
Dashboard
% slabdemo.m
% INITIALIZE MATLAB
close all;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
make b? %% DASHBOARD
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Enough to allow
b n2 n2 = 1.0;
a = 3*lam0;
% GRID
b = 5*lam0;
What grid resolution NRES = 10;
should we use? dx = lam0/NRES;
17
2/18/2020
% COMPUTE GRID
b n2 Sx = a + 2*b;
Nx = ceil(Sx/dx);
Sx = Nx*dx;
xa = [0.5:Nx-0.5]*dx;
nx1
a n1 Sx nx
xa = xa - mean(xa);
b n2 % BUILD N
N = zeros(Nx,1);
N(1:nx1-1) = n2;
N(nx1:nx2) = n1;
N(nx2+1:Nx) = n2;
Slide 35
D
N e
% SOLVE EIGEN-VALUE PROBLEM
2
x
2
neff
2
ey A = DX2 + N^2;
y [V,D] = eig(full(A));
NEFF = sqrt(diag(D));
A
Slide 36
18
2/18/2020
% SORT MODES
[~,ind] = sort(real(NEFF),'descend');
V = V(:,ind);
NEFF = NEFF(ind);
Slide 38
19
2/18/2020
Notes
• Higher-order modes converge slower.
• Higher-order modes have a smaller neff.
0
x
NRES
Slide 39
Spacer Region b
Metal
Remember the Dirichlet boundary conditions?
Values outside of the grid are forced to zero.
b n2
This means we really are simulating a slab
waveguide inside of a large metal waveguide.
a n1
It is only possible to get an accurate simulation
of the slab waveguide when the metal
b n2 waveguide is large enough.
Slide 40
20
2/18/2020
Slide 41
Slide 42
21
2/18/2020
Slide 43
Notes
• Under normal circumstances, the
spacer region size can be ~0.250.
• Modes near cutoff require larger
spacer regions to resolve.
• Thin waveguides may require larger
spacer regions.
• Always check for convergence of
spacer region size.
Slide 44
22