Notebooks and calculation files for: Modeling of the 3-Coupled-Core Fiber: Comparison Between Scalar and Vector Random Coupling Models

The files with simulation results for JLT submission "Modeling of the 3-Coupled-Core Fiber: Comparison Between Scalar and Vector Random Coupling Modelsr".

"3CCF_supermodes" file is the Mathematica code which enables to calculate supermodes (eigenvectors of M(w)) and their propagation constants of 3-coupled-core fiber (4CCF). These results are uploaded to the python notebook "3CCF_modelingJLTPaper" in order to plot them to get Fig. 3 in the paper. "TransferMatrix" is the python file with functions used for modeling, simulation and plotting. It is also uploaded in the python notebook "3CCF_modelingJLTPaper", where all the calculations for figures in the paper are presented.

! UPD 25.09.2023: There is an error in the formula of birefringence calculation. It is in the function "CouplingCoefficients" in  "TransferMatrix" file. There the variable "birefringence" has to be calculated according to the formula (19) [A. Ankiewicz, A. Snyder, and X.-H. Zheng, “Coupling between parallel optical fiber cores–critical examination”, Journal of Lightwave Technology, vol. 4, no. 9,pp. 1317–1323, 1986]:

(4*U**2*W*spec.k0(W)*spec.kn(2, W_)/(spec.k1(W)*V**4))*((spec.iv(1, W)/spec.k1(W))-(spec.iv(2, W)/spec.k0(W)))

The correct formula gives almost the same result (the difference is 10^-5), but one has to use a correct formula anyway.

P.s. In case of any questions or suggestions or if you need more explanations, you are welcome to write me an email ekader@chalmers.se. If it seems like the code does not work or mistakes in simulations are found, I also appreciate letting me know.