Mode solver semivectorial

import numpy as np
import matplotlib.pyplot as plt
import modes as ms

ms.mode_solver_semi(n_modes=2, thickness=0.22, width=0.5, wavelength=1.55)

<modes._mode_solver_semi_vectorial.ModeSolverSemiVectorial at 0x7f090c463610>

ms.mode_solver_semi(n_modes=2, thickness=0.22, width=0.5, wavelength=1.55, logscale=True)

<modes._mode_solver_semi_vectorial.ModeSolverSemiVectorial at 0x7f08d900fd90>

ms.mode_solver_semi?

Convergence tests

Lets change the simulations size and see how the neff changes

sub_widths = np.linspace(1.5, 3, 3)
n_effs = [ms.mode_solver_semi(plot=False, n_modes=1, sub_width=sub_width).results['n_effs'][0].real for sub_width in sub_widths]

plt.plot(sub_widths, n_effs)
plt.xlabel('simulation width (um)')
plt.ylabel('neff')

Text(0, 0.5, 'neff')

ds = np.linspace(0.01, 0.02, 3)
n_effs = [ms.mode_solver_semi(plot=False, n_modes=1, y_step=d, x_step=d).results['n_effs'][0].real for d in ds]

plt.plot(ds*1e3, n_effs)
plt.xlabel('mesh size (nm)')
plt.ylabel('neff')

Text(0, 0.5, 'neff')