Optical waveguide modes

Optical waveguide modes#

Hide code cell source
from collections import OrderedDict

import matplotlib.pyplot as plt
import numpy as np
import shapely
import shapely.affinity
from scipy.constants import epsilon_0, speed_of_light
from shapely.ops import clip_by_rect
from skfem import Basis, ElementTriP0
from skfem.io.meshio import from_meshio

from femwell.maxwell.waveguide import compute_modes
from femwell.mesh import mesh_from_OrderedDict
from femwell.visualization import plot_domains

We describe the geometry using shapely. In this case it’s simple: we use a shapely.box for the waveguide. For the surrounding we buffer the core and clip it to the part below the waveguide for the box. The remaining buffer is used as the clad. For the core we set the resolution to 30nm and let it fall of over 500nm

wg_width = 2.5
wg_thickness = 0.3
core = shapely.geometry.box(-wg_width / 2, 0, +wg_width / 2, wg_thickness)
env = shapely.affinity.scale(core.buffer(5, resolution=8), xfact=0.5)

polygons = OrderedDict(
    core=core,
    box=clip_by_rect(env, -np.inf, -np.inf, np.inf, 0),
    clad=clip_by_rect(env, -np.inf, 0, np.inf, np.inf),
)

resolutions = dict(core={"resolution": 0.03, "distance": 0.5})

mesh = from_meshio(mesh_from_OrderedDict(polygons, resolutions, default_resolution_max=10))
mesh.draw().show()
../../_images/cd1563b677ee80171a44c39cb63fd86144941a7dbd99296c7962bd4e2aed607a.png

Let’s also plot the domains

plot_domains(mesh)
plt.show()
../../_images/0a59406e2e9ea1896fd6a2c7244e55e16559f59edb894c6eb83c67b8c0996098.png

On this mesh, we define the epsilon. We do this by setting domainwise the epsilon to the squared refractive index.

basis0 = Basis(mesh, ElementTriP0())
epsilon = basis0.zeros()
for subdomain, n in {"core": 1.9963, "box": 1.444, "clad": 1}.items():
    epsilon[basis0.get_dofs(elements=subdomain)] = n**2
basis0.plot(epsilon, colorbar=True).show()
../../_images/fae81c866ccc0239c69f3e4cd211cde3520fb55f3771085ae9138363efb5f84e.png

And now we call compute_modes to calculate the modes of the waveguide we set up. As modes can have complex fields as soon as the epsilon gets complex, so we get a complex field for each mode. Here we show only the real part of the mode.

wavelength = 1.55

modes = compute_modes(basis0, epsilon, wavelength=wavelength, num_modes=2, order=2)
for mode in modes:
    print(f"Effective refractive index: {mode.n_eff:.4f}")
    mode.show("E", part="real", colorbar=True)
    mode.show("E", part="imag", colorbar=True)
Effective refractive index: 1.5996+0.0000j
../../_images/30ca8625177dd76171e0822959b233b244310f33fe2f01a1dd1b30498490afe3.png ../../_images/ea3263367698aaae1f1cd3793a19879a45cf14b729144f5f2f3daf890aef1309.png
Effective refractive index: 1.5202+0.0000j
../../_images/a822a2472129eaf88f7d2bd37bac6506b9a2f2fcba9a104c4a723e791a0f039f.png ../../_images/38380aac667d812d51908b0b269bf7a006c0fe8b573d6d2b5ac5210282a1d7b1.png

The intensity can be plotted directly from the mode object +

modes[0].show("I", colorbar=True)
../../_images/7fa6f2684de1e504f85f0c2e553ff194d51ddbfc2c0bddc357e1ebdcf512ebce.png

Now, let’s calculate with the modes: What percentage of the mode is within the core for the calculated modes?

powers_in_waveguide = []
confinement_factors_waveguide = []

for mode in modes:
    powers_in_waveguide.append(mode.calculate_power(elements="core"))
    confinement_factors_waveguide.append(mode.calculate_confinement_factor(elements="core"))
print(powers_in_waveguide)
print(confinement_factors_waveguide)
[(0.604909264262037+0j), (0.5754481622844256+0j)]
[(0.7520342958493351+0j), (0.7441578331749016+0j)]