Abstract
Mid-infrared supercontinuum generation in dispersion-engineered ZnSe rib waveguides was investigated for the first time. Numerical results showed that the zero-dispersion wavelength can be shifted to a shorter wavelength by adjusting structural parameters and refractive index contrast between the core and cladding layers in the waveguide. The optical field can be well confined in the 4- and 8-μm wide waveguides with a 2-μm thick cladding layer of Ge5As10S85 glass. The effect of waveguide parameters on the bandwidth of the supercontinuum spectrum at a 5-cm-long waveguide was also simulated to understand the effect of the pump wavelength and structure parameters on the supercontinuum generation. Our results showed that supercontinuum output could vary over a wide range depending on structural parameters of the waveguide, the pump power and wavelength. An ultrabroad supercontinuum spectrum from 3.0 up to 12.2 μm (> 2 octaves) was confirmed in a 4 μm-width waveguide with a peak pump power of 20 kW and a pump wavelength of 4.5 μm, which is promising as one of the on-chip supercontinuum light sources for many applications such as biomedical imaging, and environmental and industrial sensing in the mid-infrared.
Supercontinuum (SC) generation in nonlinear materials excited by ultrashort pulses of peak power has raised great research interest, due to its wide applications in various fields
Except for the adjustable ZDW in optical fiber, planar waveguide structures are also compatible with the well-developed semiconductor processing, which have better scalability and low fabrication cost. It has been demonstrated that supercontinuum spectra can be excited in ChG planar waveguides using a high peak power of few kW
In this paper, we used ZnSe as core materials and simulated the dispersion, optical field distribution, the nonlinear coefficient, and effective mode area in the waveguide consisting of bottom and cladding layers with different thicknesses and refractive index contrast. Different cladding materials and waveguide structures affect the zero-dispersion wavelength. Furthermore, we numerically simulated SC generation in a 5-cm-long dispersion-engineered ZnSe rib waveguide pumped at different wavelength from 3 to 4.5 μm and different peak power from 100 W to 20 kW, and demonstrated a broadband SC spectrum from 3.0 to 12.2 μm in a 4-μm-width waveguide.
The schematic diagram of the designed ZnSe rib waveguide is shown in

Fig. 1 (a) ZnSe waveguide structure, (b) refractive index of the core layer ZnSe, the cladding layer Ga2As30S68 and Ge5As10S85.
图1 (a) ZnSe波导结构,(b)ZnSe芯层、Ga2As30S68和Ge5As10S85覆盖层折射率
The Sellmeier equation used in the simulation of the wavelength-dependent linear refractive index n over the entire wavelength range of ZnSe, Ga2As30S68, and Ge5As10S85 glass is given by
. | (1) |
The waveguide structures were optimized by using commercial software (COMSOL). With numerical analysis, the effective refractive index can be calculated by the Finite Element Modeling solver. Subsequently, the effective index was used for calculating the dispersion parameter as well as all other higher-order dispersion parameters.
The dispersion parameter curves of the modes can be calculated using:
, | (2) |
here is the wavelength in micrometers, D is the dispersion parameter of the transmission mode, is the effective refractive index of the fundamental mode and is the light speed in vacuum. We investigated the dispersion of the waveguide, in which the top and bottom cladding layers were fixed at 2 µm, the total thickness of ZnSe film was 3 µm, and the ratio of to was variable.
Dispersion of the waveguide plays a significant role in SC generation. Ideally, dispersion near the pump wavelength should be anomalous as well as relatively fla

Fig. 2 Calculated dispersion curves of the fundamental quasi-TE (a) the dispersion parameter curves for the fundamental quasi-TE mode calculated from for eight waveguide geometries employing Ga2As30S68 glass for both the upper and lower claddings, (b) and (c) Map of the dispersion parameter of w=4 and 8 μm ZnSe rib waveguides as a function of core thickness and wavelength, respectively, employing Ge5As10S85 glass for both the upper and lower claddings. The dash lines show the change of the ZDWs
图2 TE基模的色散曲线 (a)采用Ga2As30S68做包层的8种波导结构的有效折射率计算TE基模的色散,(b-c)上下包层为Ge5As10S85 4和8 μm宽不同结构的ZnSe脊波导色散图,虚线表示0色散点位置
Therefore, all simulations in the rest part of the paper were performed in the waveguide with a refractive index contrast of 0.4, since the larger refractive index contrast can tune the ZDW to a shorter wavelength as shown in

Fig. 3 The optical filed distribution for quasi-TE polarization in the waveguide (a-c) for w=4 μm, and (d-f) for w=8 μm waveguide with =2 μm and=1 μm at a wavelength of 2, 6, and 10 μm, respectively
图3 (a-c)宽4 μm和(d-f)宽8 μm, H1=2 μm和H2=1 μm的波导在2、6和10 μm波长处TE光场分布
The nonlinear coefficient (Kerr effect), γ, which is determined by the effective mode area of the waveguide and the nonlinear refractive index of the material, can be calculated using the following formula,
, | (3) |
where is the nonlinear refractive index of ZnSe
, | (4) |
where E is the electric field’s transverse component propagating inside the waveguide.
The effective mode area and Kerr nonlinearity coefficient were calculated and shown in

Fig. 4 Effective area and nonlinear coefficient of the fundamental mode calculated in the waveguides (a) w = 4 μm, =2 μm, and=1 μm (b) w = 8μm, =2 μm, and=1 μm. (c) Dispersion distribution of the waveguides with w = 4 and 8μm. (d) The second-order dispersion of the waveguides with w = 4 and 8 μm.
图4 计算波导基模有效面积和非线性系数(a) 宽4 μm (b) 宽8 μm, H1=2 μm, H2=1 μm。宽4和8 μm波导的(c)色散分布 (d) 二阶色散
We performed the simulations of SC generation by using a generalized nonlinear Schrdinger equation (GNLSE) with a chirp-free Gaussian-shaped pump pulse as the initial condition
, | (5) |
where A is the electrical field amplitude, A(z,T) is the electric field wave amplitude as a function of propagation distance and time, T=t- is the retarded time frame moving at the group velocity vg = , α is the linear propagation loss of the waveguide including a wavelength-independent propagation loss of 0.6 dB/cm
, | (6) |
where the delayed Raman contribution is given by:
. | (7) |
We calculated the Raman gain from the data available, and τ1 and τ2 from the linewidth of Raman spectrum of ZnSe, the response function coefficient is 0.08, τ1 is 21.06 fs, and τ2 is 4.4 ps for ZnSe waveguide
The GNLSE for the fundamental quasi-TE mode of the waveguides was calculated by using commercial software (MATLAB) to simulate SC generation. For numerical analysis of SC generation in a novel 5-cm-long dispersion engineered ZnSe rib waveguide, sub-femtosecond pulses with 150 fs duration and a repetition rate of 1 kHz were used as an exciting source.

Fig. 5 Simulated SC spectra at a pump wavelength of (a) 3.0 µm, (b) 4.5 µm, (c) 3.0 µm and (d) 4.5 µm for the two waveguides at different peak power up to 20 kW, respectively.
图5 在泵浦波长为(a)3.0 μm (b)4.5 μm (c)3.0 μm和(d)4.5 μm的两个波导在不同峰值功率时的SC仿真光谱
For SC generation, self-phase modulation (SPM) alters the broadening rate imposed on the pulse by the group-velocity dispersion (GVD), and this has a correlation with the optical solitons in the abnormal dispersion region of the waveguide. The group velocity dispersion (β2) in

Fig. 6 Simulated SC spectra at different pump wavelengths of (a) 3.0 µm, (b) 3.5 µm, (c) 4.0 µm, and (d) 4.5 µm for the two waveguides with a peak power of 20 kW, respectively.
图6 两个波导在20 kW峰值功率不同泵浦波长(a)3.0 μm (b)3.5 μm (c)4.0 μm (d)4.5 μm的SC光谱
The spectral evolution corresponding to two curves in

Fig. 7 The spectral evolution plots and temporal density plots corresponding to two curves in Fig. 5.2 (d) at a peak power of 20 kW for (a) the spectral evolution plot of 4 μm waveguide, (b) the spectral evolution plot of 8 μm waveguide, (c) the temporal density plot of 4 μm waveguide, (d) the temporal density plot of 8 μm waveguide.
图7 对应图5.2 (d) 20千瓦峰值功率的SC光谱和时间谱,(a)4 μm波导SC光谱,(b)8 μm波导SC光谱,(c)4 μm波导时间谱,(d)8 μm波导时间谱
When the dispersion length is comparable to the length of the waveguide, both GVD and SPM contribute to the formation of the solitons in the abnormal dispersion region
We have numerically analyzed dispersion parameters, optical field distribution, nonlinear coefficient, and SC generation in ZnSe rib waveguide. It was found that ZDW can be shifted to shorter wavelength with increasing refractive index contrast between the core and cladding layer in the waveguide, and the optical field distribution can be well confined in 4- and 8-μm-width waveguide employing Ge5As10S85 glass as both the upper and lower claddings. With increasing pump wavelength from 3.0 μm to 4.5 μm, SC spectrum broadens, and an ultrabroad SC spectrum can be obtained up to 9.2 μm (> 2 octaves) in a waveguide pumped by a peak power of 20 kW and a wavelength of 4.5 μm. Furthermore, the simulations confirm that the SC generation is initiated by self-phase modulation, followed by soliton dynamics and soliton self-frequency shift, both of which increase with increasing pump wavelength.
References
Dudley J M, Taylor J R. Supercontinuum generation in optical fibers [M]. New York: Cambridge University Press, 2010. 10.1017/cbo9780511750465.005 [Baidu Scholar]
Eggleton B J, Lutherdavies B, Richardson K. Chalcogenide photonics [J]. Nature Photonics, 2011, 5:141-148. 10.1038/nphoton.2011.309 [Baidu Scholar]
Petersen C R, Møller U, Kubat I, et al. Mid-infrared supercontinuum covering the 1.4–13.3 μm molecular fingerprint region using ultra-high NA chalcogenide step-index fibre [J]. Nature Photonics, 2014, 8:830. 10.1038/nphoton.2014.213 [Baidu Scholar]
Gao W Q, Amraoui E M, Liao M S,et al. Mid-infrared supercontinuum generation in a suspended-core As2S3 chalcogenide microstructured optical fiber [J]. Optics Express, 2013, 21:9573. 10.1364/oe.21.009573 [Baidu Scholar]
Karim M R, Rahman B M , Agrawal G P. Mid-infrared supercontinuum generation using dispersion -engineered Ge(11.5)As(24)Se(64.5) chalcogenide channel waveguide [J]. Opt Express, 2015, 23:6903-6914. 10.1364/oe.23.006903 [Baidu Scholar]
Liao M S, Gao W Q, Cheng T L, et al.Five-octave-spanning supercontinuum generation in fluride glass[J]. App.Phy.Exp, 2013,6 (032503):1-3 . 10.7567/apex.6.032503 [Baidu Scholar]
Pigeon J J, Tochitsky S Y, Gong Cet al. Supercontinuum generation from 2 to 20 mum in GaAs pumped by picosecond CO(2) laser pulses[J]. Opt Lett, 2014, 39:3246-3249. 10.1364/ol.39.003246 [Baidu Scholar]
Kulkarni O P, Alexander V V, Kumar M, et al. 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]. J. Opt. Soc. Am. B,2011 28(10):2486–2498. 10.1364/josab.28.002486 [Baidu Scholar]
Porcel M A G, Schepers F, Epping J P, et al. Two-octave spanning supercontinuum generation in stoichiometric silicon nitride waveguides pumped at telecom wavelengths[J]. Optics Express,2017, 25:1542-1554. 10.1364/oe.25.001542 [Baidu Scholar]
Saini T S, Kumar A, Sinha R K. Design and modeling of dispersion engineered rib waveguide for ultra-broadband M-IR supercontinuum generation[J]. Journal of Modern Optics, 2017, 64(2):143-149. 10.1080/09500340.2016.1216190 [Baidu Scholar]
Karim M R, Ahmad H, Ghosh Set al. Design of dispersion-engineered As2Se3 channel waveguide for mid-infrared region supercontinuum generation editors-pick[J]. J. Appl. Phys, 2018, 123: 213101. 10.1063/1.5033494 [Baidu Scholar]
Zhao Z M, Wang X S, Dai S X, et al. 1.5–14 μm midinfrared supercontinuum generation in a low-loss Te-based chalcogenide step-index fiber[J]. Opt. Lett, 2016, 41:5222-5225. 10.1364/ol.41.005222 [Baidu Scholar]
Zhao Z M, Wu B, Wang X S, et al. Mid-infrared supercontinuum covering 2.0-16 μm in a low-loss telluride single-mode fiber[J]. Laser & Photonics Reviews, 2017,11:1700005. 10.1002/lpor.201700005 [Baidu Scholar]
Yu Y, Gai X, Wang T, et al. Mid-infrared supercontinuum generation in chalcogenides[J]. Optical Materials Express, 2013, 3:1075. 10.1364/ome.3.001075 [Baidu Scholar]
Dantanarayana H G, Abdel-Moneim N, Tang Z, et al. Refractive index dispersion of chalcogenide glasses for ultra-high numerical-aperture fiber for mid-infrared supercontinuum generation[J]. Optical Materials Express, 2014, 4:1444. 10.1364/OME.4.001444 [Baidu Scholar]
Kubat I, Agger C S, Møller U, et al. Mid-infrared supercontinuum generation to 12.5μm in large NA chalcogenide step-index fibres pumped at 4.5μm[J]. Opt. Express,2014, 22:19169-19182. 10.1364/OE.22.019169 [Baidu Scholar]
Al-kadry A, Baker C, Amraoui M E, et al. Broadband supercontinuum generation in As2Se3 chalcogenide wires by avoiding the two-photon absorption effects[J]. Opt. Lett, 2013, 38:1185-1187. 10.1364/ol.38.001185 [Baidu Scholar]
Yu Y, Gai X, Ma P, et al. Experimental demonstration of linearly polarized 2-10 mum supercontinuum generation in a chalcogenide rib waveguide[J]. Opt. Let , 2016 41:958-961. 10.1364/ol.41.000958 [Baidu Scholar]
Gai X, Choi D-Y, Madden S, Y, et al. Supercontinuum generation in the mid-infrared from a dispersion-engineered As2S3 glass rib waveguide[J]. Opt. Lett, 2012, 37: 3870-3872. 10.1364/ol.37.003870 [Baidu Scholar]
Al-kadry A, Amraoui M E, Messaddeq Y, et al. Two octaves mid-infrared supercontinuum generation in As2Se 3 microwires[J]. Opt. Express, 2014, 22: 31131-31137. 10.1364/oe.22.031131 [Baidu Scholar]
Wang R P, Yan K L, Yang Z Y, et al. Structural and physical properties of Ge11.5As24S64.5·xSe64.5·(1 - X) glasses[J]. Journal of Non Crystalline Solids, 2015, 427: 16-19. 10.1016/j.jnoncrysol.2015.07.034 [Baidu Scholar]
Li Q L, Qi D F, Wang X S, et al. Femto- and nano-second laser-induced damages in chalcogenide glasses[J]. Japanese Journal of Applied, 2019, 58: 080911. 10.7567/1347-4065/ab3368 [Baidu Scholar]
Krola H, Grezes-Besset C, Gallais L, et al. Study of laser-induced damage at 2 microns on coated and uncoated ZnSe substrates[J]. SPIE,2006, 6403:640316. 10.1117/12.695742 [Baidu Scholar]
Ma W Q, Wang L L, Zhang P Q, et al. Surface damage and threshold determination of Ge–As–Se glasses in femtosecond pulsed laser micromachining[J]. Journal of the American Ceramic Society , 2020, 103(1) 94-102. 10.1111/jace.16699 [Baidu Scholar]
The refractive index of ZnSe.Online at https://refractiveindex.info/?shelf=main&book=ZnSe&page=Connolly. 10.1016/0020-0891(83)90011-8 [Baidu Scholar]
Durand M, Houard A, Lim K, et al. Study of filamentation threshold in zinc selenide[J]. Opt Express, 2014, 22:5852-5858. 10.1364/oe.22.005852 [Baidu Scholar]
Wang T, Gai X, Wei W H, et al. Systematic z-scan measurements of the third order nonlinearity of chalcogenide glasses[J]. Opt. Mater. Express, 2014, 4: 1011. 10.1364/ome.4.001011 [Baidu Scholar]
Suminas R , Tamosauskas G, Valiulis G, et al. Multi-octave spanning nonlinear interactions induced by femtosecond filamentation in polycrystalline ZnSe[J]. Appl Phys Lett, 2017, 110(241106):1-4. 10.1063/1.4986440 [Baidu Scholar]
Mouawad O, Bejot P, Billard F, et al. Filament-induced visible-to-mid-IR supercontinuum in a ZnSe crystal: Towards multi-octave supercontinuum absorption spectroscopy[J]. Optical Materials, 2016, 60: 355. 10.1016/j.optmat.2016.08.009 [Baidu Scholar]
Werner K, Hastings M G, Schweinsberg A, et al. Ultrafast mid-infrared high harmonic and supercontinuum generation with n2 characterization in zinc selenide[J]. Optics Express, 2019, 27(3):2867. 10.1364/oe.27.002867 [Baidu Scholar]
Suminas R, Marcinkeviciute A, Tamosauskas Get al. Even and odd harmonics-enhanced supercontinuum generation in zinc-blende semiconductors[J]. Journal of the Optical Society of America B, 2019,36(2):A22. 10.1364/josab.36.000a22 [Baidu Scholar]
Mittal V, Sessions N P, Wilkinson J S, et al. Optical quality ZnSe films and low loss waveguides on Si substrates for mid-infrared applications[J]. Optical Materials Express, 2017, 7: 712-745. 10.1364/ome.7.000712 [Baidu Scholar]
Mittal V, Nedeljkovic M, Rowe D J, et al. Chalcogenide glass waveguides with paper-based fluidics for mid-infrared absorption spectroscopy[J]. Opt. Lett. 2018, 43: 2913. 10.1364/ol.43.002913 [Baidu Scholar]
Borisova Z U. Glassy Semiconductors [M]. New York: Plenum Press,1981. 10.1007/978-1-4757-0851-6 [Baidu Scholar]
Zhu M, Liu H, Li X, et al. Ultrabroadband flat dispersion tailoring of dual-slot silicon waveguides[J]. Optics Express, 2012, 20:15899-15907. 10.1364/OE.20.015899 [Baidu Scholar]
Karim M R, Rahman B M, Agrawal G P. Dispersion engineered Ge(1)(1).(5)As(2)(4) Se(6)(4).(5) nanowire for supercontinuum generation: a parametric study[J]. Opt. Express, 2014, 22: 3102931040. 10.1364/oe.22.031029 [Baidu Scholar]
Gai X, Madden S, Choi D-Y, et al. Dispersion engineered Ge11.5As24Se64.5 nanowires with a nonlinear parameter of 136W-1m-1 at 1 550 nm[J]. Optics Express, 2010,18: 18866-18874. 10.1364/oe.18.018866 [Baidu Scholar]
Hutchings D C, Van Stryland E W. Nondegenerate two-photon absorption in zinc blende semiconductors[J]. Journal of the Optical Society of America B, 1992, 9:2065-2074. 10.1364/josab.9.002065 [Baidu Scholar]
Anand S, Verma P, Jain K P, et al. Temperature dependence of optical phonon lifetimes in ZnSe[J]. Physica B, 1996, 226:331-337. 10.1016/0921-4526(96)84974-x [Baidu Scholar]
Agrawal G. Nonlinear fiber optics[M]. New York: Elsevier Press, 2013. 10.1016/b978-0-12-397023-7.00012-7 [Baidu Scholar]