Robust split-step Fourier methods for simulating the propagation of ultra-short pulses in single- and two-mode optical communication fibers

Ralf Deiterding,Stephen W. Poole

DOI: https://doi.org/10.48550/arXiv.1504.01331

2015-04-07

Abstract:Extensions of the split-step Fourier method (SSFM) for Schrödinger-type pulse propagation equations for simulating femto-second pulses in single- and two-mode optical communication fibers are developed and tested for Gaussian pulses. The core idea of the proposed numerical methods is to adopt an operator splitting approach, in which the nonlinear sub-operator, consisting of Kerr nonlinearity, the self-steepening and stimulated Raman scattering terms, is reformulated using Madelung transformation into a quasilinear first-order system of signal intensity and phase. A second-order accurate upwind numerical method is derived rigorously for the resulting system in the single-mode case; a straightforward extension of this method is used to approximate the four-dimensional system resulting from the nonlinearities of the chosen two-mode model. Benchmark SSFM computations of prototypical ultra-fast communication pulses in idealized single- and two-mode fibers with homogeneous and alternating dispersion parameters and also high nonlinearity demonstrate the reliable convergence behavior and robustness of the proposed approach.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: to develop and test a robust split - step Fourier method (SSFM) for simulating the propagation of ultrashort pulses in single - mode and two - mode fibers. Specifically, this research aims to deal with the propagation problems of femtosecond - level ultrashort pulses in single - mode and two - mode fibers. Due to their extremely short time widths (T0 < 100 fs), such ultrashort pulses will be significantly affected by nonlinear effects during propagation, such as Kerr nonlinearity, self - steepening, and stimulated Raman scattering, etc. Therefore, traditional numerical methods may not be able to accurately simulate these complex phenomena.

To solve these problems, the author proposes a new class of SSFM based on the operator splitting method and reformulates the nonlinear sub - operator as a first - order quasilinear system of signal intensity and phase through the Madelung transformation. For the single - mode case, an upwind numerical method with second - order accuracy is derived; while for the two - mode case, this method is extended to approximately handle the four - dimensional system caused by the selected two - mode model. Finally, the reliable convergence behavior and robustness of this method are verified through benchmark calculations in idealized single - mode and two - mode fibers.

In summary, the main objective of this paper is to improve the existing numerical methods so that they can more accurately simulate the propagation characteristics of ultrashort pulses in fiber - optic communication systems, especially in the presence of high nonlinearity and dispersion parameters. This helps to improve the design and performance optimization of future high - speed optical communication networks.

### Key formulas involved

1. **Basic equation for the propagation of ultrashort pulses in single - mode fibers**:
\[
\frac{\partial A}{\partial z}+\frac{\alpha}{2}A+\left(\sum_{k \geq 1}\frac{\beta_{k}i^{k - 1}}{k!}\frac{\partial^{k}}{\partial t^{k}}\right)A = i\gamma\left(1 + i\frac{1}{\omega_{0}}\frac{\partial}{\partial t}\right)\left[A\int_{-\infty}^{\infty}R(t')|A(t - t')|^{2}dt'\right]
\]

2. **Simplified equation considering dispersion and nonlinear effects**:
\[
\frac{\partial A}{\partial z}+\frac{\alpha}{2}A+\beta_{1}\frac{\partial A}{\partial t}+i\frac{\beta_{2}}{2}\frac{\partial^{2}A}{\partial t^{2}}-\frac{\beta_{3}}{6}\frac{\partial^{3}A}{\partial t^{3}}=i\gamma\left(A|A|^{2}+i\frac{1}{\omega_{0}}\frac{\partial}{\partial t}(A|A|^{2})-T_{R}A\frac{\partial|A|^{2}}{\partial t}\right)
\]

3. **Equation after introducing group velocity**:
\[
\frac{\partial A}{\partial z}+\frac{\alpha}{2}A+i\frac{\beta_{2}}{2}\frac{\partial^{2}A}{\partial T^{2}}-\frac{\beta_{3}}{6}\frac{\partial^{3}A}{\partial T^{3}}=i\gamma\left(A|A|^{2}+iS\frac{\partial}{\partial T}(A|A|^{2})-T_{R}A\frac{\partial|A|^{2}}{\partial T}\right)
\]

4. **Propagation of two interacting ultrashort pulses in two - mode fibers**