Mid-infrared supercontinuum generation using dispersion-engineered Ge 11 . 5 As 24 Se 64 . 5 chalcogenide channel waveguide

We numerically investigate mid-infrared supercontinuum (SC) generation in dispersion-engineered, air-clad, Ge11.5As24Se64.5 chalcogenide-glass channel waveguides employing two different materials, Ge11.5As24S64.5 or MgF2 glass for their lower cladding. We study the effect of waveguide parameters on the bandwidth of the SC at the output of 1-cm-long waveguide. Our results show that output can vary over a wide range depending on its design and the pump wavelength employed. At the pump wavelength of 2 μm the SC never extended beyond 4.5 μm for any of our designs. However, supercontinuum could be extended to beyond 5 μm for a pump wavelength of 3.1 μm. A broadband SC spanning from 2 μm to 6 μm and extending over 1.5 octave could be generated with a moderate peak power of 500 W at a pump wavelength of 3.1 μm using an air-clad, all-chalcogenide, channel waveguide. We show that SC can be extended even further when MgF2 glass is used for the lower cladding of chalcogenide waveguide. Our numerical simulations produced SC spectra covering the wavelength range 1.8-7.7 μm (> two octaves) by using this geometry. Both ranges exceed the broadest SC bandwidths reported so far. Moreover, we realize it using 3.1 μm pump source and relatively low peak power pulses. By employing the same pump source, we show that SC spectra can cover a wavelength range of 1.8-11 μm (> 2.5 octaves) in a channel waveguide employing MgF2 glass for its lower cladding with a moderate peak power of 3000 W. © 2015 Optical Society of America OCIS codes: (000.4430) Numerical approximation and analysis; (130.2755) Glass waveguides; (190.4390) Nonlinear optics, devices; (320.6629) Dispersion; (320.6629) Supercontinuum generation. References and links 1. X. Gai, T. Han, A. Prasad, S. Madden, D. Y. Choi, R. Wang, D. Bulla, and B. Luther-Davies, “Progress in optical waveguides fabricated from chalcogenide glasses,” Opt. Express 18(25), 26635–26646 (2010). 2. J. Hu, C. R. Menyuk, L. B. Shaw, J. S. Sanghera, and I. D. Aggarwal, “Maximizing the bandwidth of supercontinuum generation in As2Se3 chalcogenide fibers,” Opt. Express 18(3), 6722–6739 (2010). 3. I. D. Aggarwal and J. S. Sanghera, “Development and applications of chalcogenide glass optical fibers at NRL,” J. Optoelectron. Adv. Mater. 4(3), 665–678 (2002). #233410 $15.00 USD Received 28 Jan 2015; revised 26 Feb 2015; accepted 26 Feb 2015; published 5 Mar 2015 (C) 2015 OSA 9 Mar 2015 | Vol. 23, No. 5 | DOI:10.1364/OE.23.006903 | OPTICS EXPRESS 6903 4. B. J. Eggleton, B. Luther-Davies, and K. Richardson, “Chalcogenide photonics,” Nat. Photonics 5, 141–148 (2011). 5. F. Luan, M. D. Pelusi, M. R. E. Lamont, D. Y. Choi, S. J. Madden, B. Luther-Davies, and B. J. Eggleton, “Dispersion engineered As2S3 planar waveguides for broadband four-wave mixing based wavelength conversion of 49 Gb/s signals,” Opt. Express 17(5), 3514–3520 (2009). 6. S. J. Madden, D. Y. Choi, D. A. Bulla, A. V. Rode, B. Luther-Davies, V. G. Ta’eed, M. D. Pelusi, and B. J. Eggleton, “Long, low loss etched As2S3 chalcogenide for all-optical signal regeneration,” Opt. Express 15(22), 14414–14421 (2007). 7. M. R. E. Lamont, C. M. Sterke, and B. J. Eggleton, “Dispersion engineering of highly nonlinear As2S3 waveguides for parametric gain and wavelength conversion,” Opt. Express 15(15), 9458–9463 (2007). 8. M. R. Karim, B. M. A. Rahman, and G.P. Agrawal, “Dispersion engineered Ge11.5As24Se64.5 nanowire for supercontinuum generation: A parametric study,” Opt. Express 22(25), 31029–31040 (2014). 9. D. D. Hudson, E. C. Mägi, A. C. Judge, S. A. Dekker, and B. J. Eggleton, “Highly nonlinera chalcogenide glass micro/nanofiber devices: Design, theory, and octave-spanning spectral generation,” Opt. Commun. 285, 4660– 4669 (2012). 10. P. Ma, D. Y. Choi, Y. Yu, X. Gai, Z. Yang, S. Debbarma, S. Madden, and B. Luther-Davies, “Low-loss chalcogenide waveguides for chemical sensing in the mid-infrared,” Opt. Express 21(24), 29927–29937 (2013). 11. M. R. E. Lamont, B. Luther-Davies, D. Y. Choi, S. Madden, and B. J. Eggleton, “Supercontinuum generation in dispersion engineered highly nonlinear (γ = 10 /W/m) As2S3 chalcogenide planar waveguide,” Opt. Express 16(19), 14938–14944 (2008). 12. X. Gai, S. Madden, D. Y. Choi, D. Bulla, and B. Luther-Davies, “Dispersion engineered Ge11.5As24Se64.5 nanowires with a nonlinear parameter of 136 W−1m−1 at 1550 nm,” Opt. Express 18(18), 18866–18874 (2010). 13. X. Gai, D. Choi, S. Madden, Z. Yang, R. Wang, and B. Luther-Devies, “Supercontinuum generation in the midinfrared from a dispersion-engineered As2S3 glass rib waveguide,” Opt. Lett. 37(18), 3870–3872 (2012). 14. Y. Yu, X. Gai, T. Wang, P. Ma, R. Wang, Z. Yang, D. Choi, S. Madden, and B. Luther-Davies, “Mid-infrared supercontinuum generation in chalcogenides,” Opt. Express 3(8), 1075–1086 (2013). 15. Y. Yu, X. Gai, P. Ma, D. Choi, Z. Yang, R. Wang, S. Debbarma, S. J. Madden, and B. Luther-Davies, “A broadband, quasi-continuous, mid-infrared supercontinuum generated in a chalcogenide glass waveguide,” Laser Photonics Rev., 1–7 (2014). 16. R. J. Weiblen, A. Docherty, J. Hu, and C. R. Menyuk, “Calculation of the expected bandwidth for a mid-infrared supercontinuum source based on As2S3 chalcogenide photonic crystal fibers,” Opt. Express 18(25), 26666– 26674 (2010). 17. L. B. Shaw, R. R. Gattass, J. S. Sanghera, and I. D. Aggarwal, “All-fiber mid-IR supercontinuum source from 1.5 to 5 μm,” Proc. SPIE 7914 (79140P), 1–5 (2011). 18. A. Marandi, C. W. Rudy, V. G. Plotnichenko, E. M. Dianov, K. L. Vodopyanov, and R. L. Byer, “Mid-infrared supercontinuum generation in tapered chalcogenide fiber for producing octave spanning frequency comb around 3 μm,” Opt. Express 20(22), 24218–24225 (2012). 19. I. Savelii, O. Mouawad, J. Fatome, B. Kibler, F. Desevedavy, G. Gadret, J. C. Jules, P. Y. Bony, H. Kawashima, W. Gao, T. Kohoutek, T. Suzuki, Y. Ohishi, and F. Smektala, “Mid-infrared 2000-nm bandwidth supercontinuum generation in suspended-core microstructured Sulphide and Tellurite optical fibers,” Opt. Express 20(24), 27083– 27093 (2012). 20. W. Gao, M. E. Amraoui, M. Liao, H. Kawashima, Z. Duan, D. Deng, T. Cheng, T. Suzuki, Y. Messaddeq, and Y. Ohishi, “Mid-infrared supecontinuum generation in a suspended-core As2S3 chalcogenide microstructured optical fiber,” Opt. Express 21(8), 9573–9583 (2013). 21. N. Granzow, M.A. Schmidt, W. Chang, L. Wang, Q. Coulombier, J. Troles, P. Toupin, I. Hartl, K. F. Lee, M. E. Fermann, L. Wondraczek, and P. St. J. Russell, “Mid-infrared supercontinuum generation in As2S3 “nano-spike” step index waveguide,” Opt. Express 21(9), 10969–10977 (2013). 22. C. Wei, X. Zhu, R. A. Norwood, F. Seng, and N. Peyghambarian, “Numerical investigation on high power midinfrared supercontinuum fiber lasers pumped at 3 μm,” Opt. Express 21(24), 29488–29504 (2013). 23. I. Kubat, C. R. Petersen, U. V. Moller, A. B. Seddon, T. M. Benson, L. Brilland, D. Mechin, P. M. Moselund, and O. Bang, “Thulium pumped mid-infrared 0.9-9 μm supercontinuum generation in concatenated fluoride and chalcogenide glass fibers,” Opt. Express 22(4), 3959–3967 (2014). 24. I. Kubat, C. S. Agger, U. Moller, A. B. Seddon, Z. Tang, S. Sujecki, T. M. Benson, D. Furniss, S. Lamarini, K. Scholle, P. Fuhrberg, B. Napier, M. Farries, J. Ward, P. M. Moselund, and O. Bang, “Mid-infrared supercontinuum generation to 12.5 μm in large NA chalcogenide step-index fibers pumped at 4.5μm,” Opt. Express 22(16), 19169–19182 (2014). 25. D. D. Hudson, M. Baudisch, D. Werdehausen, B. J. Eggleton, and J. Biegert, “1.9 octave supercontinuum generation in a As2S3 step-index fiber driven by mid-IR OPCPA,” Opt. Lett. 39(19), 5752–5755 (2014). 26. A. Al-Kadry, M. E. Amraoui, Y. Messaddeq, and M. Rochette, “Two octaves mid-infrared supercontinuum generation in As2Se3 microwires,” Opt. Express 22(25), 31131–31137 (2014). 27. C. R. Petersen, U. Mφ ller, I. Kubat, B. Zhou, S. Dupont, J. Ramsay, T. Benson, S. Sujecki, M. Abdel-Moneim, Z. Tang, D. Furniss, A. Seddon, and O. Bang, “Mid-infrared supercontinuum covering the 1.4-13.3 μm molecular #233410 $15.00 USD Received 28 Jan 2015; revised 26 Feb 2015; accepted 26 Feb 2015; published 5 Mar 2015 (C) 2015 OSA 9 Mar 2015 | Vol. 23, No. 5 | DOI:10.1364/OE.23.006903 | OPTICS EXPRESS 6904 fingerprint region using ultra-high NA chalcogenide step-index fibre,” Nat. Photonics 8, 830–834 (2014). 28. U. Mφ ller, Y. Yu, I. Kubat, C. R. Petersen, X. Gai, L. Brilland, D. Mechin, C. Caillaud, J. Troles, B. LutherDavies, and O. Bang, “Multi-milliwatt mid-infrared supercontinuum generation in a suspended core chalcogenide fiber,” Opt. Express 23(3), 3282–3291 (2015). 29. J. H. Kim, M. Chen, C. Yang, J. Lee, S. Yin, P. Ruffin, E. Edwards, C. Brantley, and C. Luo, “Broadband IR supercontinuum generation using single crystal sapphire fibers,” Opt. Express 16(6), 4085–4093 (2008). 30. B. Kuyken, X, Liu, R. M. Osgood Jr., R. Baets, G. Roelkens, and W. M. J. Green, “Mid-infrared to telecom-band supercontinuum generation in highly nonlinear silicon-on-insulator wire waveguides ,” Opt. Express 19(21), 20172–20181 (2011). 31. O. P. Kulkarni, V. V. Alexander, M. Kumar, M. J. Freeman, M. N. Islam, F. L. Terry Jr., M. Neelakandan, and A. Chan “Supercontinuum generation from 1.9 to 4.5 μm in ZBLAN fiber with high average power generation beyond 3.8 μm using a thulium-doped fiber amplifier,” J. Opt. Soc. Am. B 28(10), 2486–2498 (2011). 32. M. Liao, W. Gao, T. Cheng, X. Xue, Z. Duan, D. Deng, H. Kawashima, T. Suzuki, and Y. Ohishi, “Five-octavespanning supercontinuum generation in fluride glass,” App. Phy. Exp. 6(032503) 1–3 (2013). 33. E. A. Anashkina, A. V. Andrianov, M. Y. Koptev, S. V. Muravyev, and A. V. Kim, “Towards mid-infrared supercontinuum generation with germano-silicate fibers,” IEEE J. of Sel. Top. in Quan. Elect. 20 (5), 7600608 (2014). 34. M. Bass, G. Li, and E. V. Stryland, Hand Book of Optics Vol-IV 3rd ed. (The McGraw-Hill, New York, 2010). 35. D. Yeom, E. C. Mägi, M. R. E. Lamont, M. A. F. Ro


Introduction
Recently chalcogenide glasses (ChGs) have emerged as promising nonlinear materials having a number of unique properties that makes them attractive for fabricating planar optical waveguides and using them for application such as mid-infrared (MIR) supercontinuum (SC) generation and optical sensing [1].SC generation in MIR has increasingly become a focus for research because bright MIR light sources can be used for molecular fingerprint spectroscopy, frequency metrology, optical coherent tomography and microscopy [2,3].Chalcogenide glasses can provide MIR transparency with sulphides transmitting to beyond 8.5 µm, selenides to 14 µm and tellurites to around 20 µm [4].A major advantage of chalcogenide glasses is their large ultra-fast nonlinearity amongst all glasses, making them good materials for planar waveguides to generate SC in the MIR [5][6][7][8][9].A few of these ChG materials such as As 2 S 3 , As 2 Se 3 , Ge 11.5 As 24 S 64.5 and Ge 11.5 As 24 Se 64.5 glasses are highly suitable for making active and passive devices in the MIR region.Amongst them Ge 11.5 As 24 Se 64.5 glass has excellent film-forming properties with high thermal and optical stability under intense illumination [10].Recently interest has grown in designing and optimizing planar waveguides made from Ge 11.5 As 24 Se 64.5 chalcogenide glass for broadband MIR SC generation with suitably tailored group-velocity dispersion (GVD), including a zero-dispersion wavelength (ZDW) close to the central wavelength of the pump [11,12].
In recent years, several experimental and theoretical investigations on mid-infrared SC generation were reported in ChG planar waveguides [13][14][15], ChG fibers [16][17][18][19][20][21][22][23][24][25][26][27][28], single crystal sapphire fibers [29], silicon-on-insulator (SOI) wire waveguides [30], ZBLAN fiber [31], fluoride glass [32] and germano-silicate fibers [33].Gai et al. [13] reported SC generation from 2.9 µm to 4.2 µm in dispersion-engineered As 2 S 3 glass rib waveguide (6.6 cm long) pumped with 7.5 ps duration pulses at a wavelength of 3.26 µm with a pulse peak power of around 2 kW.Yu et al. [14] reported SC generation up to 4.7 µm using a 4.7-cm-long As 2 S 3 glass rib waveguide employing MgF 2 glass as a substrate and pumped with 7.5 ps duration with a peak power of 1 kW pulses at a wavelength of 3.26 µm.They observed a flat SC extending from 2.5 µm to 7.5 µm in 5-mm-thick bulk sample of Ge 11.5 As 24 Se 64.5 glass pumped with 150 fs duration pulses with up to 20 MW peak power at a wavelength of 5.3 µm.They also reported theoretically that SC could be generated beyond 10 µm in a dispersion-engineered all-chalcogenide Ge 11.5 As 24 Se 64.5 glass rib waveguide pumped with 250 fs duration pulses at a wavelength of 4 µm or longer.Yu et al. [15] recently reported the generation of broadband SC spanning from 1.8 µm to 7.5 µm in dispersion-engineered Ge 11.5 As 24 Se 64.5 glass 1-cm-long rib waveguide pumped at 4 µm using 320 fs pulses from a periodically-poled lithium niobate optical parametric amplifier pumped with commercially available mode-locked Yb laser with a peak power of 3260 W. For step-index fiber based on As 2 S 3 , Hudson et al. [25] reported 1.9-octave SC extending from 1.6 µm to 5.9 µm pumped at a wavelength of 3.1 µm with a peak power of 520 kW.Al-Kadry et al. [26] reported MIR SC generation from 1.1 µm to 4.4 µm in 10-cm-long As 2 Se 3 microwires pumped with 800 fs duration pulses at a wavelength of 1.94 µm with energy of 500 pJ.Petersen et al. [27] used 100 fs pulses at 6.3 µm to generate SC covering the range 1.4-13.3µm in 85 mm long As 2 S 3 step-index fiber with a 16 µm core.Mφ ller et al. [28] used a 18 cm long, low-loss, suspended core As 38 S 62 fiber and generated a SC spanning from 1.7 to 7.5 µm by pumping this fiber with 320 fs pulses with a peak power of 5.2 kW at a wavelength of 4.4 µm.The SC generation in As 2 S 3 fibers and waveguides has been limited to 5 µm due to their increased losses in the long wavelength region [13,17,18].Although As 2 Se 3 chalcogenide glass has its loss edge in the long wavelength region around 16-17 µm, the use of such materials for SC generation has so far been limited owing to high peak power pump sources required in the MIR region [4,26,27].
In this paper, we demonstrate numerically that a 1-cm-long dispersion-engineered Ge 11.5 As 24 Se 64.5 glass rectangular channel waveguide can generate broadband SC in the MIR regime.In our previous work [8], a dispersion-engineered Ge 11.5 As 24 Se 64.5 nanowire was designed, with polymer as top and silica as bottom claddings, for low-power wideband SC generation.To extend the SC to the MIR region and to ensure that the claddings as well as the core materials are transparent at such long wavelengths, rectangular channel waveguides were designed with air on top and either Ge 11.5 As 24 S 64.5 or MgF 2 glass as the lower cladding material.To tailor the dispersion to obtain the ZDW close to the pump wavelength, two separate geometries were employed, one with a lower cladding of Ge 11.5 As 24 S 64.5 glass having lower index contrast and another with a lower cladding of MgF 2 glass having larger index contrast.Both were designed and optimized by adjusting the waveguides dimensions.In earlier experiments it was found that the extension of SC depends on the pump wavelength, among many other factors.To achieve sufficiently long wavelength extension in MIR region, the pump wavelength will need to be used around 4-5 µm or longer.The nonlinear parameter, γ = 2πn 2 /λ A eff , where n 2 is the nonlinear refractive index and A eff is the effective area of the mode, of the ChG waveguide decreases with the increasing pump wavelength, λ .As A eff also increases, high pump powers would be required to generate SC in the long-wavelength region, which can cause damage to the ChG waveguides if relatively wider pump pulses are employed [14].For this reason we choose, for our numerical simulation, sub-picosecond pulses of 85 fs (FWHM) duration at a pump wavelength of 3.1 µm with a repetition rate of 160 kHz [45] and with low to moderate peak powers.The propagation loss (α) and nonlinear refractive index (n 2 ) were taken to be 0.5 dB/cm and 4.3 × 10 −18 m 2 /W, respectively [10,15].Since the actual n 2 of ChG glasses has not been measured in the MIR region, its measured value at 1.55 µm [12] was reduced by a factor of two at the pump wavelength of 3.1 µm [15].
We carry out simulations for two different structures of ChG channel waveguides assuming that the input pulse excites the fundamental quasi-TE mode.The peak power varies in the range of 25-500 W and 100-3000 W for two pump sources at wavelength of 2 µm and 3.1 µm, respectively.Using the pump source at 2 µm with a peak power of 500 W, the SC bandwidth extended in the range 1.3-3.3µm and 1.3-3.5 µm for the waveguides designed and optimized with Ge 11.5 As 24 S 64.5 glass or MgF 2 glass acting as lower cladding, respectively.A broadband MIR SC extending over 1.5 octave from 2 µm to 6 µm could be generated with a lower cladding of Ge 11.5 As 24 S 64.5 glass.The SC extended over more than 2 octave (1.8 µm to 7.7 µm) when a lower cladding of MgF 2 glass was employed.We calculate the SC bandwidth at the -30 dB level from the peak.Our calculated bandwidths are the largest reported so far for SC generated using chalcogenide glass waveguides pumped at a wavelength of 3.1 µm with a moderate peak power of 500 W.

Theory
The schematic diagram of the Ge 11.5 As 24 Se 64.5 chalcogenide glass channel waveguide used in our simulations is shown in Fig. 1.The wavelength-dependent linear refractive index of Ge 11.5 As 24 Se 64.5 , Ge 11.5 As 24 S 64.5 and MgF 2 glasses over the entire wavelength range used in the simulation was obtained using the Sellmeier equation, Here λ is the wavelength in micrometers.We use the Sellmeier coefficients calculated by fitting smooth curves to the measured refractive index data for the ChGs [10] and the MgF 2 [34] glasses.The values of the integer m and the fitting coefficients are given in Table 1.Dispersion of the waveguide plays an important role in determining the SC spectrum.Ideally the dispersion near the pump wavelength should be anomalous with a small value that does not change much [35][36][37][38][39][40].Small flattened anomalous dispersion with high nonlinearity is required for large frequency shifts of the Raman soliton.ChG glasses generally have smaller material dispersion than silicon and to obtain the ZDW close to the pump wavelength, relatively larger waveguide dispersion is required to offset the material dispersion [12,41].We use a finite element (FE) based mode-solver [42] to obtain the propagation constant β (ω) of the fundamental mode over a range of frequencies and calculate the effective index, n eff = λ β (ω)/2π from this β (ω), which is subsequently used for numerically calculating the GVD parameter, dλ 2 (ps/nm/km) as well as all other higher-order dispersion parameters.Spectral broadening of a SC mainly depends on these dispersion parameters and the nonlinear coefficient, γ which in turn depends on the nonlinear refractive index of the material, n 2 and effective mode area of the waveguide.Pump wavelength is also an important factor for the extension of SC spectrum in the long wavelength region.Optimized mode area, A eff can be obtained with our FE mode-solver and γ of the waveguide can be calculated using A eff .
The FE approach used here is based on the vector-H-field formulation, since it is one of the most accurate and numerically efficient approaches to obtain the modal field profiles and mode propagation constants β (ω) of various quasi-TE and quasi-TM modes.The full-vectorial formulation is based on the minimization of the full H-field energy functional [42], where H is the full-vectorial magnetic field, * denotes a complex conjugate and transpose, ω 2 is the eigenvalue (ω being the angular frequency), p is a weighting factor for the penalty term to eliminate spurious modes and ε and μ are the permittivity and permeability tensors, respectively.The two-dimensional cross-section of the waveguide is discretized by using a large number of triangular elements.All three components of the magnetic fields can be represented as piece-wise polynomials within the elements.With a proper choice of waveguide discretization we can accurately calculate the energy functional by integrating it over each element.
To study SC generation in the MIR region, simulations were performed using a generalized nonlinear Schrödinger equation (GNLSE) for the slowly varying envelope of the pulse [8,12]: where A is the electrical field amplitude, α is the linear propagation loss of the waveguide, is the k th order dispersion parameter, and T = t − z v g is the retarded time frame moving with the group velocity v g = 1 β 1 (ω 0 ) at the pump frequency ω 0 .The nonlinear coefficient is γ = n 2 ω 0 cA eff (ω 0 ) , where n 2 is the nonlinear refractive index and c is the speed of light in vacuum, A eff (ω 0 ) is the effective area of the mode at the pump frequency ω 0 , and α 2 = 9.3 × 10 −14 m/W is the two-photon absorption coefficient [12].Finally the material response function includes both the instantaneous electronic response (Kerr type) and the delayed Raman response and has the form The f R represents the fractional contribution of the delayed Raman response.The other parameters used for the simulation are listed in Table 2.

Numerical results
As the MIR SC generation requires the cladding as well as core materials of a waveguide to be transparent at long wavelengths, we first choose a channel waveguide fabricated using only chalcogenide materials with a rectangular core made of Ge 11.5 As 24 Se 64.5 glass and using As 36 S 64 glass for both the upper and lower claddings, whose refractive index of 2.37 at 1.55 µm provides an index contrast of 0.3.By varying waveguide dimensions for realizing ZDW in the range of 2-4 µm or longer, we numerically calculated GVD for three different structures as a function of wavelength for the fundamental quasi-TE mode and plot them with material dispersion in Fig. 2. It can be observed from this figure that the material dispersion, shown by a To study SC generation, we solve the generalized nonlinear Schrödinger equation (GNLSE) given in Eq. ( 3).We used the same computational method that was used in Karim at el. [8].That study has indicated that there is a possibility to obtain spurious results if the adequate number of dispersion terms is not included during SC simulations.We solved Eq. ( 3) with the split-step Fourier method (SSFM) including up to 10 th order dispersion.We tested accuracy of the finiteelement modal solution for a channel waveguide by a powerful extrapolation technique called Aitken's extrapolation [44].We observed the convergence between the raw FEM results and extrapolated values as the number of FE elements increases along the transverse dimensions of the waveguide.We also tested the accuracy of numerically calculated GVD parameters for different number of elements used along the transverse dimensions of a channel waveguide.We fitted the dispersion data with a Taylor series expansion including up to 10 th order dispersion and found that this method accurately reproduced the GVD curves obtained by FE mode-solver over the whole wavelength range of calculations.A sufficiently small time step was chosen to accommodate the spectral expansion of the generated SC during numerical simulations.
We consider a waveguide containing a core of dimensions W = 1.5 µm and H = 1.1 µm with air and Ge 11.5 As 24 S 64.5 glass as upper and lower claddings, respectively.Another geometry has a core with W = 1.35 µm and H = 0.4 µm but MgF 2 glass as the lower cladding.Using FE mode-solver, we obtain A eff = 1.09 µm 2 which yields γ = 24.79/W/m for the structure with Ge 11.5 As 24 S 64.5 lower cladding and A eff = 0.51 µm 2 which yields γ = 53.39/W/m for the structure with MgF 2 lower cladding at a pump wavelength of 2 µm.The GVD parameter calculated at the pump wavelength for both waveguides has a value of 13 ps/nm/km.These waveguides supported propagation of the fundamental TE mode up to the cut-off wavelength near 3.1 µm for the air-clad all-ChG structure and around 4 µm for the structure employing MgF 2 glass as lower cladding.A sech pulse of 150 fs duration (FWHM) was launched with peak power between 25 W and 500 W for numerical simulations.We included a wavelength-independent propagation loss of 2.5 dB/cm [10] for our 1-cm-long rectangular channel waveguides.Figure 4 shows the predicted SC spectra for the two waveguides at three different power levels at a pump wavelength of 2 µm.In the case of 100 W input power the SC spectrum extends from 1.5 µm to around 3 µm, producing a -30 dB bandwidth of 1500 nm for the waveguide using Ge 11.5 As 24 S 64.5 glass for its lower cladding.At the same power level the SC spectrum extends from 1.4 µm to around 3.1 µm producing a -30 dB bandwidth of 1700 nm for the waveguide with MgF 2 glass for its lower cladding.After increasing input peak power level at 500 W, the SC spectra broadened from 1.3 µm to 3.3 µm (output bandwidth of 2000 nm) and from 1.3 µm to 3.5 µm (output bandwidth of 2200 nm) for these two waveguides, respectively.The spectral evolution plots corresponding to Fig. 4(c) are shown in Fig. 5.It is apparent that a larger output bandwidth can be realized by using a waveguide employing MgF 2 glass for its lower cladding.The reason behind this bandwidth enhancement solely related to the higher index contrast between the core and cladding materials when MgF 2 glass is used for the lower cladding.However, it was not possible to extend the SC spectrum to beyond 3.5 µm by using a pump wavelength of 2 µm.
To extend the SC in the MIR regime, we need to shift the pump wavelength toward longer wavelengths [14].We employ a pump wavelength of 3.1 µm since such a pump source has been realized at a repetition rate of 160 kHz [45].Using 85 fs duration pulses at a wavelength of 3.1 µm, we focus on a waveguide with W = 4 µm and H = 1.6 µm with Ge 11.5 As 24 S 64.5 glass for its lower lower cladding.The second waveguide structure had W = 5 µm and H = 0.95 µm but employed MgF 2 glass for its lower cladding.Using the FE mode-solver, we obtain A eff = 4.25 µm 2 and γ = 2.05 /W/m for the first structure and A eff = 3.32 µm 2 and γ = 2.63 /W/m for the second structure.The GVD parameter calculated at the pump wavelength for these waveguides has values of 10.22 ps/nm/km and 21 ps/nm/km, respectively.These waveguides supported propagation of the fundamental TE mode up to the cut-off wavelength around 6.5 µm for the air-clad all-ChG structure and beyond 12 µm for the structure employing MgF 2 as lower cladding.
After evaluating higher-order dispersion terms up to tenth-order from GVD curves shown in Fig. 3(b), we performed numerical simulations for SC generation in both the waveguide geometries at a power level between 100 W and 3000 W. The dispersion length, L D = T 2 P / |β 2 | for he 85 fs pump pulse for these two structures is 45 mm and 28 mm and the nonlinear length, L NL = 1/γP at a peak power of 500 W is 0.98 mm and 0.76 mm, respectively.The soliton order, N = L D /L NL was 7 and 6 for the two waveguides, respectively, and the soliton fission length, L fiss ≈ L D /N was found to be 6.6 mm and 4.1 mm, respectively.Figure 6 shows the predicted spectra at various power level up to 3000 W for the two waveguide geometries at a pump wavelength of 3.1 µm.The spectral evolution plots corresponding to Fig. 6(c) and 6(d) are shown in Fig. 7.It is observed from Fig. 7 that the SC generation is dominated by soliton fission, which results in many short pulses generated through the soliton fission process whose spectra shifted towards the long wavelength side of the input spectrum.Raman-induced frequency shift (RIFS) is reducing gradually as solitons moved towards the second ZDWs located near 4.15 µm and 4.4 µm for the two waveguide geometries, respectively.Due to the spectral recoil effect [46,47], RIFS were completely suppressed near the second ZDWs.At the same time, nonsolitonic radiation in the form of dispersive wave is produced at a wavelength that lies beyond the second ZDW.For the wavguide geometry employing Ge 11.5 As 24 S 64.5 glass as a lower cladding, it can be observed from Fig. 6 that the SC extends over 4 µm covering a wavelength range from 2 µm to 6 µm and for the structure employing MgF 2 as a bottom cladding, the SC extends over 6 µm covering a wavelength range from 1.8 µm to around 7.7 µm both at a peak power of 500 W. By increasing power level up to 3000 W, one can generate a SC spectrum that extends from 2 µm to 10 µm producing a bandwidth of 8 µm (> 2 octave) for the first waveguide and bandwidth of 9.2 µm (> 2.5 octave) for the second waveguide employing MgF 2 glass for its lower cladding.
There may be a disadvantage for ChG waveguides fabricated using MgF 2 glass for the lower cladding owing to its fragility.This problem becomes severe for long waveguides but is manageable for short waveguides around 1-cm-long [14].Ma et al. [10] have shown experimentally that, if the top surface of the waveguides made from chalcogenide materials is left uncoated, the bare waveguide surface becomes rapidly contaminated by absorbing water and hydrocarbons from surrounding environment, resulting in increased losses.To prevent the surface contamination, a thin 10-nm protective coating layer of fluoro-polymer may be placed on top of the waveguide geometries proposed in this paper.We have tested numerically placing of such a thin layer and found that the dispersions at a pump wavelength of 3.1 µm for the two waveguides used in our work was 9.85 ps/nm/km (10.23 ps/nm/km without coating) and 21.68 ps/nm/km (21 ps/nm/km without coating).Such relatively small changes in the β 2 value do not produce noticeable changes on the SC generated at the waveguide output.
Recently Yu et al. [15] proposed an air-clad rib waveguide whose core was made with Ge 11.5 As 24 Se 64.5 glass and whose lower cladding was made with Ge 11.5 As 24 S 64.5 glass.They were able to generate SC covering the wavelength range 1.8-7.5 µm when pumped with an optical parametric amplifier at 4 µm with a peak power of 3260 W. By rigorous numerical simulations, it is shown here that MIR SC can be generated with dispersion-engineered airclad rectangular channel waveguide employing the same materials and covering the wavelength range of 1.9-7 µm by employing a pump source at a wavelength of 3.1 µm with a peak power of 3000 W. It is observed from Fig. 6 that MIR SC can be extended in the long wavelength regime up to 11 µm for a waveguide employing MgF 2 glass for its lower cladding when pumped with the same power and pump source although MgF 2 glass would be expected to absorb beyond 9 µm [10].

Conclusions
We have numerically demonstrated MIR SC generation by using dispersion-engineered, airclad, channel waveguides designed and optimized such that they use either Ge 11.5 As 24 S 64.5 ChG glass or MgF 2 glass for its lower cladding material.The SC is generated with the waveguides proposed here by using pump pulses with low to moderate peak power at wavelength near 2 µm or 3.1 µm.Although the nonlinear parameter has larger values at a pump wavelength of 2 µm, the SC spectra were extended only over 1.3-3.3µm and 1.3-3.5 µm, respectively even at the highest but moderate peak power of 500 W for the two proposed structures.To extend SC to beyond the 5 µm, it is necessary to choose a longer pump wavelength around 4-5 µm.
To realize the ZDW of the waveguide around this wavelength the core size of the waveguide must be increased which increases the effective mode area and hence reduces the nonlinear parameter.The only solution is to increase the peak power of the input pulse toward MW levels which can damage the input facet of the waveguide.Considering these factors, we designed and optimized air-clad two rectangular channel waveguides by employing two different lower cladding materials and choose 3.1 µm for the pump wavelength such that broadband SC could be generated at moderate peak power levels.
Using pump source at a wavelength of 3.1 µm with a relatively low peak power of 500 W we obtained a SC spectrum extended over 1.5 octave and covering the wavelength range from 2-6 µm with an air-clad all-ChG structure.However, when we used an optimized waveguide using MgF 2 glass for its lower cladding, the SC spectrum extended over more than two-octaves covered a wavelength range of 1.8-7.7 µm.This, to the best of our knowledge, is the broadest MID SC at such a low peak power and at a lower pump wavelength of 3.1 µm using a chalcogenide waveguide.We have also found that, MID SC can be extended over in the wavelength range 1.9-7 µm and 1.8-11 µm with a moderate peak power of 3000 W when the air-clad channel waveguide is designed with a Ge 11.5 As 24 Se 64.5 glass core and employs Ge 11.5 As 24 S 64.5 or MgF 2 glass for its lower cladding, respectively.

Fig. 2 .
Fig. 2. GVD curves for the fundamental quasi-TE mode calculated from n eff for three waveguides geometries employing As 36 S 64 glass for both the upper and lower claddings.The black solid line curve shows the material dispersion curve for comparison.

Fig. 3 .
Fig. 3. GVD curves for the waveguide geometries employing two different lower claddings (solid black curve for Ge 11.5 As 24 S 64.5 and red dashed curve for MgF 2 ) for the fundamental quasi-TE mode (a) at a pump wavelength of 2 µm and (b) at a pump wavelength of 3.1 µm.Vertical dotted line indicates the position of pump wavelength.

Fig. 4 .
Fig. 4. Simulated SC spectra at a pump wavelength of 2 µm for (a) air-clad all-chalcogenide waveguide at peak power from 25, 100, and 500 W; (b) air-clad chalcogenide core employing MgF 2 for its lower cladding at the same power levels;(c) waveguides with two different lower claddings at a peak power of 500 W only.

Fig. 6 .
Fig. 6.Simulated SC spectra at a pump wavelength of 3.1 µm for (a) air-clad allchalcogenide waveguide at peak power between 100 W and 3000 W; (b) air-clad chalcogenide core employing MgF 2 for its lower cladding for the same power levels; (c) waveguides employing with two different lower claddings at a peak power of 500 W only; (d) waveguides with two different lower claddings at peak power of 3000 W only.