Numerical Estimation of Multimode Fiber Modes and Propagation Constants: 

Part 1: Straight Fibers

Under the weakly guided approximation, analytical solutions for the mode profiles of step-index (SI) and graded-index (GRIN) multimode fibers (MMF) can be found [1]. It also gives a semi-analytical solution for the dispersion relation in SI MMFs, and, by adding stronger approximations, an analytical solution for the parabolic profile GRIN MMFs [2] (note that those approximations do fail for lower order modes). An arbitrary index profile requires numerical simulations to estimate the mode profiles and the corresponding propagation constants of the modes. I present in this tutorial how to numerically estimate the scalar solution for the profiles and propagation constants of guides modes in multimode circular waveguide with arbitrary index profile and in the presence of bending. I released a beta version of the Python module pyMMF based on such an approach [3]. It relies on expressing the transverse Helmholtz equation as an eigenvalue problem. Solutions are found by finding the eigenvectors of a large but sparse matrix representing the equation on the discretized space.

Modes of step-index multimode fibers

Scattering media were the first type of "complex media" for which wavefront shaping techniques were applied. Quickly, applications were developed for multimode fibers as well. One can consider multimode fiber as a complex media; because of its inherent modal dispersion (different modes travel at different speeds) and also because of the possible coupling between modes, the output field of the fiber does not resemble its input one. Wavefront shaping in multimode fibers has had a fast development because of its applications in biomedical endoscopic imaging and for telecommunications, where the exploitation of the spatial modes in multimode fibers offers a promising way to increase data rates compared to single-mode fibers. 

I present here quickly the expression of the modes of a step-index multimode fiber and the so-called linearly polarized modes, that are convenient for manipulation using shaping techniques.