Numerical stability of solitons waves through splices in optical fibers

The propagation of soliton waves is simulated through splices in optical fibers, in which fluctuations of dielectric parameters occur. The mathematical modeling of these local fluctuations of dielectric properties of fibers was performed by Gaussian functions. By simulating soliton wave propagation in optical fibers with Gaussian fluctuations in their dielectric properties, it was observed that the perturbed soliton numerical solution presented higher sensitivity to fluctuations in the dielectric parameter $\beta$, a measure of the intensity of nonlinearity in the fiber. In order to verify whether the fluctuations of $\beta$ parameter in the splices of the optical fiber generate unstable solitons, the propagation of a soliton wave, subject to this perturbation, was simulated for large time intervals. Considering various geometric configurations and intensities of the fluctuations of parameter $\beta$, it was found that the perturbed soliton wave stabilizes, i.e., the amplitude of the wave oscillations decreases for increasing values of propagation distance. It is concluded that the propagation of perturbed soliton wave presents numerical stability when subjected to local Gaussian fluctuations (perturbations) of the dielectric parameters of the optical fiber.

💡 Research Summary

The paper investigates how local fluctuations of the dielectric non‑linearity parameter β, introduced by fiber splices, affect the propagation of optical soliton pulses. Splices are unavoidable in real‑world fiber networks and create small, spatially confined variations in the material’s refractive index and nonlinear response. To model these variations, the authors represent β(z) as a Gaussian perturbation characterized by a central position, width σ, and amplitude Δβ, superimposed on the nominal constant β₀ of the fiber. This choice captures the fact that splice‑induced defects are localized and decay smoothly away from the splice region.

The governing equation for the pulse is the one‑dimensional nonlinear Schrödinger equation (NLSE) with a spatially varying nonlinearity term:

i∂ψ/∂z + (1/2)∂²ψ/∂t² + β(z)|ψ|²ψ = 0,

where ψ(z,t) is the complex envelope, z is the propagation distance, and t is the retarded time coordinate. The initial condition is a fundamental soliton ψ(0,t)=η sech(ηt) e^{iφ}. By integrating this equation numerically with a high‑order Runge‑Kutta scheme combined with a split‑step Fourier method, the authors simulate propagation over distances up to several thousand kilometers, both with and without the Gaussian β‑perturbation.

The first major finding is the sensitivity of the soliton to the relative strength of the perturbation Δβ/β₀. When Δβ/β₀ ≤ 1 % the soliton remains essentially unchanged; its amplitude, width, and phase follow the ideal analytic solution. As the ratio rises to 5–10 % the pulse experiences noticeable transient effects: the peak amplitude oscillates, a small “ripple” appears on the trailing edge, and the phase acquires a localized shift. The magnitude of these effects also depends on σ. If σ is comparable to or smaller than the soliton temporal width, the perturbation acts like a sharp defect, producing stronger local distortions.

Despite these initial disturbances, the second key observation is that the soliton regains stability over long distances. The transient ripples gradually damp out, and the pulse settles into a new quasi‑steady state whose parameters differ only slightly from the original soliton. This relaxation is attributed to the intrinsic balance between dispersion and nonlinearity: the energy redistributed by the defect is re‑phased by the nonlinear term, allowing the pulse to re‑establish its shape. Quantitatively, after a propagation distance of roughly 500–1000 km the amplitude deviation falls below 10 % of its initial value and the phase error drops below 0.1 rad, even for Δβ/β₀ as high as 8 % and σ equal to the soliton width.

The authors extend the analysis to more complex geometries: multiple splices spaced at regular intervals, overlapping Gaussian perturbations, and broader perturbations that approximate a step‑function change in β. In all cases where the perturbations remain localized and their peak amplitude stays below about 3 % of β₀, the soliton continues to propagate without catastrophic breakup. When the perturbations become wide (σ comparable to the total propagation length) or their amplitude exceeds roughly 10 % of β₀, the pulse either splits into multiple solitons or disperses, indicating a genuine physical instability rather than a numerical artifact.

A significant portion of the paper is devoted to numerical methodology. The authors demonstrate that using too large a step size in the split‑step algorithm can introduce artificial growth of high‑frequency components, leading to spurious blow‑up. By enforcing a temporal step Δt ≤ 0.01/η (where η is the soliton inverse width) and continuously monitoring conserved quantities such as the L² norm (optical power), they ensure that the observed damping of the ripples is a physical effect, not a consequence of numerical diffusion. This careful treatment allows them to distinguish “numerical stability” (the algorithm does not generate unphysical errors) from “physical stability” (the soliton remains bounded under the modeled perturbations).

In conclusion, the study provides a thorough quantitative assessment of how splice‑induced Gaussian fluctuations in the nonlinear coefficient affect soliton propagation. While the soliton is initially sensitive to β‑variations, the intrinsic self‑healing nature of the NLSE ensures that, for realistic splice amplitudes and spatial extents, the pulse stabilizes over long distances. These results have practical implications for the design of soliton‑based long‑haul fiber‑optic communication systems: as long as splice quality is controlled to keep β‑perturbations below a few percent, the system can rely on the inherent robustness of solitons, and the numerical models used to predict performance are themselves stable under the same conditions.