Midinfrared optical rogue waves in soft glass photonic crystal fiber

: We investigate numerically the formation of extreme events or rogue waves in soft glass tellurite ﬁbers and demonstrate that optical loss drastically diminishes shot-to-shot ﬂuctuations characteristic of picosecond pumped supercontinuum (SC). When loss is neglected these ﬂuctuations include extreme events such as formation of highly energetic pulses located at the red end of the spectrum and we obtain right-skewed heavy-tailed distributions characteristic of extreme events statistics. On the other hand, when loss is included bandwidth ﬂuctuations follow Gaussian-like statistical distributions. Our results thus implicitly show that rogue waves will not occur in any SC spectrum that is limited by loss, such as commercial silica ﬁber based SC sources.


Introduction
Ever since it was first reported [1], SC generation has been widely investigated and demonstrated in several pump regimes, from continuous wave to ultrashort femtosecond [2]. Invariably, novel physics are uncovered and in particular an intriguing phenomenon in the form of statistically rare highly energetic pulses located at the red end of the spectra has been predicted numerically [3,4] and observed experimentally [5] in optical fibers.
Such extreme events occur since pulses or solitons generated in a fiber can interact and transfer energy between each other during collisions. On average there will be a net energy transfer to the larger soliton and the mechanism can generate very high amplitude solitons [3,6,7]. This energy transfer and generation of strong nonlinear excitations, is in fact a generic nonlinear phenomenon well-known from biology [7][8][9], optical waveguide arrays [10], and discrete lattices [11].
The term optical "rogue wave" (RW) was coined [5] due to the high resemblance to the infamous so-called oceanic RWs. Moreover, these statistically rare events exhibited a shifting towards longer wavelengths and they distinctly appeared in the red end of the spectrum as large bandwidth fluctuations following right-skewed heavy-tailed statistics [5].
It should be noted, that in recent studies [12,13] the authors pointed to the possibility of a RW occuring inside the SC spectrum, however in the following we limit the investigations to RWs as giant solitons located on the long wavelength side of the SC.
The above results were obtained in silica fibers at low powers, with spectra far from the loss edge, however here we investigate midinfrared (MIR) RW formation in soft glass fibers. Soft glasses, which among others include chalcogenide (n 2 ∼0.8×10 −17 m 2 /W) [14] and tellurite (n 2 ∼0.6×10 −18 m 2 /W) [15], exhibit much higher nonlinearities than silica (n 2 ∼0.3×10 −19 m 2 /W) [16], and thus offer a great reduction in the required fiber length or pump power. In addition soft glass fibers allow low loss guidance at much longer wavelengths, well into the MIR regime.
Hence we consider a 60 cm tellurite fiber with pitch (Λ=3 μm), relative hole size (d/Λ=0.4) and zero-dispersion wavelength (ZDW) ∼1780 nm [17] see Fig. 1(a), in which we show the measured loss and calculated dispersion. Now, since thulium mode-locked lasers exhibit a broad and high gain spectrum in the 1.8 to 2.1 μm wavelength regime [18], we consider pumping at the thulium wavelength λ =1930 nm, i.e., anomalous dispersion region to generate a large number of solitons as in the aforementioned studies in silica.
As seen in Fig. 1(a), this fiber entails substantial loss and we demonstrate that generation of RWs can be greatly reduced when the band edge of the generated SC approaches the loss edge at ∼3 μm, similar to SC fluctuations limited by a photonic bandgap in a solid-core fiber [19]. We stress here however, that in a bandgap fiber the dispersion deviates to infinity as it approaches the edge of the gap. The rapid increase in dispersion alone will halt the solitons and this is exactly what stops them red-shifting in this case, before entering the actual regime of loss.
Therefore this has no resemblance whatsoever with the material loss edge of a conventional fiber, at which the dispersion is still finite and the solitons are still allowed to red-shift into the loss edge. Hence the dynamics is different in many ways, e.g., solitons/rogue waves approaching the loss edge will decay in power, whereas the rogue waves in a bandgap fiber will maintain their power, but spread out due to the increasing dispersion.
The broad absorption in the region 2.5-5 μm is due to water in the glass, note however that this glass was melted in dry atmosphere and thus has already a considerably reduced water content compared with glasses melted in ambient atmosphere [17]. In addition, improvement of the glass melting procedure is expected to further reduce the water content, where complete removal of water would shift the loss edge to ∼4 μm, at which confinement loss becomes significant. In comparison, the material IR edge loss is <1 dB/m for λ <4.9 μm [20] which is negligible compared to confinement loss.
We consider pulse propagation both with and without including loss and demonstrate, for the first time to our knowledge, MIR RW formation, combining a tellurite fiber with a thulium pump. We base our simulations on the generalized nonlinear Schrödinger equation (GNLSE) Fig. 1. (a) Optical loss for tellurite fiber with Λ=3 μm and d/Λ=0.4. Shown are the combined contribution from both material and confinement loss, as well as confinement loss alone. (b) Raman gain spectrumh R (Ω) of tellurite with maximum peak at Ω R ∼22.5 THz. and silica with peak at Ω R ∼13.2 THz. The dashed curve shows the approximate Raman gain spectrum of tellurite.
[16] which we solve in the so-called interaction picture [21,22], and we consider both higher order dispersion, loss, and stimulated Raman scattering, see [17] for details. In optical fibers the generation of RWs is assisted by the Raman effect [23], leading to an on-average energy transfer to the larger soliton, which is also the one red-shifting the fastest [24]. The fact that the Raman red shift increases with soliton power is exactly what enables spectral measurements of RWs [5].
We implement the Raman response function in the GNLSE, where f R =0.064 and h R (t) is derived from the measured Raman gain spectrumh R (Ω) [15], see Fig. 1(b). The tellurite Raman profile is markedly different from silica and the width (bandwidth) is substantially larger. Moreover, it consists of several peaks and thus the Raman profile h R (t) should be approximated as a sum of decaying harmonic oscillators and we use two oscillators. The approximated gain profile is shown in Fig. 1(b) (dashed curve) which is more physical since it is non-zero for all Ω>0, whereas the measured profile, due to filtering of the pump, drops off abruptly at Ω∼5 THz. The approximated Raman gain thus allows us to model small frequency shifts, which is important for narrowband pulses, e.g., ps pulses.
We pump at λ =1930 nm, with a 5 ps pulse duration, and various peak powers. The nonlinear coefficient is ∼0.14 W −1 m −1 and noise is included by adding one random phase photon per mode in the frequency domain. To capture the statistics of the RWs, we perform 200 simulations for each peak power, i.e., 1000, 1500 and 2000 W and the resulting spectra are shown in Fig. 2. In all cases a SC is developed. However, the detrimental effect of loss is evident, giving rise to much narrower SC than in the corresponding lossless case. In particular, when pumping with a peak power of 1000 W, the SC is not formed at all when loss is included (therefore not shown in Fig. 2). Without loss the SC are very broad, and the ensemble of each peak power contains several MIR RWs at ∼3.4 μm, ∼3.8 μm, and ∼4.2 μm for peaks power of 1000, 1500 and 2000 W, respectively. From the inserts on the right it follows that shot-to-shot fluctuations can be substantial, which is manifested in noticeable bandwidth deviations from the most red-shifted RW, to the narrowest spectra.
When loss is included, the bandwidth fluctuations are much smaller. Moreover, the fluctuations decrease with peak power, since increasing the peak power, forces a larger fraction of the pulses to fully develop into the maximum allowed SC bandwidth, effectively limited by optical loss.
To investigate the statistics of each ensemble, we use a long-pass filter and select components above a particular wavelength, depending on the spectral extent of the ensemble median. In particular, cut-off wavelengths from top to bottom in Fig. 2 are 3.2 μm, 3.6 μm, 2.65 μm, 4 μm, and 2.7 μm. Fourier transforming the remaining spectral components reveals a series of short pulses of varying peak powers, see Fig. 3.