Supercontinuum generation in photonic crystal fibers infiltrated with tetrachloroethylene

This study proposes a photonic crystal fiber made of fused silica glass, with the core infiltrated with tetrachloroethylene (C2Cl4) as a new source of supercontinuum (SC) spectrum. We studied numerically the guiding properties of the several different fiber structures in terms of characteristic dispersion, mode area, and attenuation of the fundamental mode. Based on the results, the structural geometries of three C2Cl4-core photonic crystal fibers were optimized in order to support the broadband SC generations. The first fiber structure with lattice constant 1.5 μm and filling factor 0.4 operates in all-normal dispersion. The SC with a broadened spectral bandwidth of 0.8–2 μm is generated by a pump pulse with a central wavelength of 1.56 μm, 90 fs duration and energy of 1.5 nJ. The second proposed structure, with lattice constant 4.0 μm and filling factor 0.45, performs an anomalous dispersion for wavelengths longer than 1.55 μm. With the same pump pulse as the first fiber, we obtained the coherence SC spectrum in an anomalous dispersion range with wavelength range from 1 to 2 μm. Meanwhile, the third selected fiber (lattice constant 1.5 μm, filling factor 0.55) has two zero dispersion wavelengths at 1.04 μm and 1.82 μm. The octave-spanning of the SC spectrum formed in this fiber was achieved in the wavelength range of 0.7–2.4 μm with an input pulse whose optical properties are 1.03 μm wavelength, 120 fs duration and energy of 2 nJ. Those fibers would be good candidates for all-fiber SC sources as cost-effective alternatives to glass core fibers.


Introduction
Supercontinuum (SC) lasers, i.e. light sources exhibiting an ultra-broad band of frequencies and good quality beam (Liu et al. 2019), have been a subject of intense research over the past decade. Their applications cover many areas such as spectroscopy, frequency metrology, biomedicine, and optical coherence tomography (Holzwarth et al. 2000;Udem et al. 2002;Tu and Boppart 2013;Sharma et al. 2008).
A number of techniques have been proposed for the generation of SC in which the high nonlinearities and low dispersion profiles play important roles. A typical approach for this purpose is utilizing photonic crystal fibers (PCFs) made of silica or highly nonlinear soft glasses (Dudley et al. 2006;Dai et al. 2018;Nguyen et al. 2019;Klimczak et al. 2019;Buczynski et al. 2010). Silica fibers work efficiently for SC generations from the visible to the near-infrared (IR) range, but not for the mid-IR range due to their limited transmission bandwidth (Dudley et al. 2006). Moreover, silica shows relatively low nonlinearity, so silica fibers require input pulses with high peak power for broad SC spectra. Although the nonlinear coefficient of the silica-core PCFs can be enhanced by reducing its effective core area, this approach will lead to difficulties in light coupling. A second common approach is to use highly nonlinear soft glass materials such as chalcogenide (Dai et al. 2018), tellurite (Nguyen et al. 2019;Klimczak et al. 2019) or lead-bismuth-gallate glass (Buczynski et al. 2010). However, these fibers are incompatible with silica for fusion splicing, limited in power handling, and usually require complex pump systems. Other disadvantages of these glasses include the following. First, they might be mechanically weak. Secondly, more complex fiber structures may be difficult to fabricate, and may require more thermal processes which can result in crystallization.
An alternative solution is liquid-infiltrated optical fibers, as demonstrated in previous research (Zhang et al. 2006;Fanjoux et al. 2017; Van et al. 2018a;Maji and Chaudhuri 2014;Pniewski et al. 2016;Ding et al. 2018;Lanh et al. 2019; Van et al. 2020). Unlike silica, nonlinear liquids show high nonlinearity that is comparable with the nonlinear soft glasses (Kedenburg et al. 2014). Some of them, e.g. carbon disulfide (CS 2 ), carbon tetrachloride (CCl 4 ) are highly transparent in the near-IR range (Churin et al. 2013;Kedenburg et al. 2012). In addition, liquids are able to be mixed to optimize their optical properties, e.g. linear refractive index for practical applications (Wohlfarth 2008). Another advantage of using liquid core PCFs is that the generated SC spectra can be controlled by the modification of temperature or pressure since it is possible to alter the PCFs' dispersion properties via these factors (Van et al. 2018b;Du et al. 2004). Therefore, the infiltration of PCFs with liquids adds an extra degree of freedom for broadening applications in potential fields (Ebnali-Heidari et al. 2012), especially, with use of mid-IR SC laser sources in medical diagnostics (Chemnit et al. 2018), biomedical technologies (Guo et al. 2004) and remote sensing (Ambrico et al. 2000).
Up to date, most of the experimental results of liquid-core fibers SC generations are focused on the use of step-index liquid-core optical fibers since they are straight forward and convenient to fill micro-capillaries with liquids, as an equivalent of a step-index optical fiber. However, this approach lacks flexibility in designing the dispersion and modal properties of the fibers (Churin et al. 2013). Recently, use of liquid-filled PCFs has been received more attention because it overcomes the disadvantages of liquids capillary approach. The obtained results indicated that it is possible to shift the ZDW and match it with a pump wavelength of the high-power commercial laser, as well as obtaining allnormal and flat dispersion regime in the expected spectra (Van et al. 2018a;Maji and Chaudhuri 2014;Ding et al. 2018). Due to the higher nonlinear refractive index of selected liquids (Ding et al. 2018;Kedenburg et al. 2014), it is possible to achieve the broadband, coherent SC spectrum either in a few centimeters long liquid-core fiber sample, or with a lower peak power (Van et al. 2020).
Most of the works on SC generation considered with carbon disulfide (CS 2 ), nitrobenzene (C 6 H 5 NO 2 ), toluene (C 7 H 8 ), carbon tetrachloride (CCl 4 ) or chloroform (CHCl 3 ) (Ding et al. 2018;Lanh et al. 2019; Van et al. 2020;Churin et al. 2013;Hoang et al. 2019Hoang et al. , 2018. Although, CS 2 has a high nonlinear refractive index n 2 = 220 × 10 -20 m 2 /W at the wavelength of 1030 nm (Kedenburg et al. 2014), 100 times higher than for fused silica n 2 = 2.19 × 10 -20 m 2 /W (Kabaciński et al. 2019) and are transparent in the visible to near-IR (0.5-2.2 µm) range (Churin et al. 2013), it is also a highly toxic, carcinogenic and explosive (Challenor 2002). C 6 H 5 NO 2 and C 7 H 8 are relatively low transparent in the near-IR range (Zhang et al. 2006;Hoang et al. 2018), resulting in a reduction of spectral broadening in the PCFs (Raei et al. 2018). Meanwhile, the advantage of CHCl 3 is its moderate toxicity when compared to CS 2 and C 7 H 8 . The upper limit of the acceptable concentration of CHCl 3 in air is 249 mg m −3 (IPCS 1994), while those limits are 188 mg m −3 and 31 mg m −3 for C 7 H 8 , CS 2 , respectively (Van et al. 2017). Its linear refractive index is lower than silica (Lanh et al. 2019). Therefore, total reflection mechanism cannot be achieved in PCFs. As a consequence, CHCl 3 was not taken in consideration. Among those liquids, CCl 4 is a better choice because of the moderate toxic property and highly transparent in the near-IR range (Kedenburg et al. 2012). However, for CCl 4 , the linear refractive index is similar to that of silica which leads to the low index contrast between the core and cladding and then results in the decrease of nonlinear coefficient (Ding et al. 2018).
In this paper, we numerically study on PCFs made of fused silica, having core infiltrated with tetrachloroethylene (C 2 Cl 4 ) as a propagating medium for SC generations. We have selected C 2 Cl 4 because it has high nonlinear refractive index n 2 = 16.75 × 10 -20 m 2 /W at λ = 1.55 µm (Chemnit et al. 2018) (seven times higher than for fused silica n 2 = 2.19 × 10 -20 m 2 /W (Kedenburg et al. 2014). In particular, C 2 Cl 4 provides a super-wide transparent from the visible to mid-IR (0.5-3.39 µm) (Takahashi et al. 1985). In addition, the upper limit of the acceptable concentration of C 2 Cl 4 in air is 345 mg m −3 (DHHS 1993). C 2 Cl 4 is then less toxic when compared to CHCl 3 , C 7 H 8 , and CS 2 which is extremely important for practical applications. In our work, dispersion properties were optimized by modifying geometrical parameters of PCF, i.e. lattice constant and filling factor in order to obtain the dispersion flatness and location of ZDW in the analyzed wavelength range. For those optimum fiber structures, we utilized the generalized nonlinear Schrödinger equation (GNLSE) to demonstrate SC generations. The advantages of the selected PCFs would also be discussed.

Numerical modeling of the PCF
The cross-section of the geometrical structure of the modeled PCF is shown in Fig. 1. We assumed that the cladding of the fiber is made of fused silica glass and consists of eight rings of air holes arranged in regular hexagonal lattice defined by the lattice pitch Λ and air holes diameter d. The linear filling factor of the cladding is defined as f = d/Λ and a constant filling factor is used for all rings to simplify future fiber development. The diameter of the central air hole has to be large enough and fulfill the relation: Dc = 2Λ-1.2d, to ensure the technical feasibility of the considered structures during the fabrication process. During numerical simulations process, the diameters of cladding holes increase with the growth of the filling factors and lattice pitches, while the core diameters decrease.
We assumed that the central hole is filled with C 2 Cl 4 . The characteristics of the real and imaginary parts of the refractive index of C 2 Cl 4 showed in Fig. 2. The refractive index characteristics modeled obeying Sellmeier equation and given by the formula: where B i and C i are Sellmeier coefficients and λ is the wavelength in μm. The Sellmeier's coefficients for fused silica and C 2 Cl 4 are presented in Table 1.
The numerical analysis was carried out with the use of Lumerical Mode Solution software (Mode Solution 2020). The perfectly matched layer is the boundary condition used for the analysis of the PCF which makes no reflection at the boundary and reduces the loss.
(1) In addition, conformal meshing was applied on all material interfaces to reduce meshing errors and increase the numerical accuracy of the simulations.

Modeling modal and dispersion properties of the PCF
Preliminary numerical simulations were conducted for the PCF having Λ changing from 1.0 to 4.5 μm with a step of 0.5 μm and f changing from 0.3 to 0.8 with a step of 0.05. A number of simulations were carried out to observe the flatness of dispersion and the ZDW shift. Each simulation was done for wavelengths in the range of 0.5 μm to 2.0 μm. Figure 3a shows the maximum dispersion of the fundamental mode for various Λ and f. For the case of the PCF with Λ is greater than or equal to 2.0 μm, the maximum dispersion is greater than zero that means we got both normal and anomalous dispersion. In this case, for a given Λ value, the maximum dispersion will increase with increasing f. Meanwhile, for a given f, the maximum dispersion values will decrease when we increase Λ. The example dispersion characteristic for this region is shown in Fig. 3b. In the rest cases, the dispersion takes values only in the normal dispersion regime or in both normal and anomalous dispersion regimes. Figure 3c shows an example of the dispersion characteristic of the  PCF in the all normal dispersion regime at locations '2' and '3' in Fig. 3a, respectively. In this case, the maximum dispersion increases with increasing f when Λ is constant. Similar behavior can also be observed when increase the lattice pitch with the filling factor is constant. In addition, increasing the lattice pitch makes dispersion flatter ultimately becomes monotonous. Figure 4 shows the ZDW shift with various Λ and f. For a given Λ value, ZDW is usually shifted toward the smaller wavelengths with the increase of filling factor f. Meanwhile, for a given f, the ZDW will be shifted toward longer wavelengths if increasing Λ. In particular, for the case of red area (corresponding with Λ = 1.0 μm and f changing from 0.65 to 0.8) and green area (corresponding with Λ = 1.5 μm and f changing from 0.45 to 0.8), the PCF has two ZDWs. In the rest cases, for the left of the red area, corresponding the fibers with Λ = 1.0 μm, it is fulfilled for f from 0.3 to 0.6 and from 0.3 to 0.4 for fiber with Λ = 1.5 μm left of the green area, the PCF with all normal dispersion have no ZDW.
The optimization criteria aiming at generations of SC at selected wavelengths (1.03 μm and 1.56 μm), include the flatness, indication of the dispersion characteristic, and the amount of difference between ZDW and pump wavelengths. On the basis of preliminary calculations, we proposed three PCFs, namely #F 1 , #F 2 and #F 3 . The geometrical parameters of the fibers are introduced in Table 2. The first fiber (#F 1 ) is designed to work in the normal dispersion range with the pump wavelength of 1.56 μm. PCF with normal dispersions does not have any ZDW which, in Fig. 4, is in the left of the red area and left of the green area as described above. Although, fibers having Λ = 1.0 μm (an example as '2'  in Fig. 3a) have all normal dispersion with f changing from 0.3 to 0.6, but their dispersion curves are less flat than those of fibers with Λ = 1.5 μm (see Fig. 3c). In addition, the broadest SC spectrum would be achieved by pumping near the local maximum of the dispersion curve. Therefore, PCF with parameters Λ = 1.5 μm, f = 0.4 was selected as optimal (point #F 1 in Fig. 4).
For fiber #F 2 with parameters Λ = 4.0 μm, f = 0.45 (point #F 2 in Fig. 4), it is expected to generate the high coherent SC in the anomalous regime at the pump wavelength of 1.56 μm. This fiber was selected because it has ZDW at 1.55 μm which is the closest one to the pumping wavelength. In addition, this fiber also has the flattest dispersion for all dispersion characteristic with Λ = 4.0 μm.
The fiber #F 3 with parameters Λ = 1.5 μm, f = 0.55 (point #F 3 in Fig. 4), is expected to generate high coherent SC in the normal regime at the pump wavelength of 1.03 μm. The reason for this selection is that the fiber has the first ZDW at 1.04 μm, which is the closest wavelength compared to the one pumped. Another reason is that the dispersion at 1.03 μm equals -2.763 ps/nm/km, which is the closest value to the horizontal axis for all dispersion characteristic with Λ = 1.5 μm. This was beneficial to spectrally broaden, effectively suppress, and the noise seeded dynamics.
The properties of these three fiber structures are studied in detail. The material losses shown in Fig. 2 were incorporated in calculations. Figure 5 depicts the dispersion characteristics for these fibers. Fiber #F 1 exhibits all-normal dispersion in the whole analyzed wavelength range and the dispersion at 1.56 μm equals − 15 ps/nm/km, while fiber #F 2 has ZDW at 1.55 μm and the dispersion at the pump wavelength is equal to 3.20 ps/nm/km. In the rest case, fiber #F 3 has first ZDW at 1.04 μm and the dispersion at 1.03 μm is equal to − 4.85 ps/nm/km.
The effective mode area and nonlinear coefficient of these structures are presented in Fig. 6. The nonlinear coefficient γ is calculated by the equation γ = n 2 ω 0 /cA eff , where ω 0 is the central frequency of the pulse, c is the speed of light in the vacuum, and A eff is the effective mode area for the fundamental mode of the fiber. The obtained results show that the mode area of the fundamental mode increases almost linearly with the wavelength. The effective mode area for the fundamental mode of fibers #F 1 and #F 2 are 4332 μm 2 , and 16.670 μm 2 for a pump wavelength of 1.56 μm, while their nonlinear coefficient is 156.9 W −1 km −1 and 40.79 W −1 km −1 , respectively. In the case of the Fig. 5 Characteristics of the fundamental mode dispersion for the fiber #F 1 , #F 2 , and #F 3 pump wavelength at 1.03 μm, the effective mode area of the fiber #F 3 equals 3591 μm 2 , while its nonlinear coefficient is 189.3 W −1 km −1 .
Attenuation characteristics of the fundamental mode as a function of the wavelength for these three fibers are shown in Fig. 7. It is clear that the attenuation characteristics follow the bulk material losses of C 2 Cl 4 , playing the major role. The highest losses are observed for fiber #F 2 because this fiber has the largest core. For the wavelength of 1.88 μm the losses reach the highest values with 0.04 dB cm −1 , 0.042 dB cm −1 , and 0.053 dB·cm −1 for #F 1 , #F 3 , and #F 2 , respectively.
For optimized fibers, we also verify whether the selected fiber was single mode or multimode at the pump wavelengths. The normalized frequency V of fiber #F 1 and fiber #F 2 calculated at 1.56 μm are around 2.24 and 2.49, respectively, so the fiber #F 1 is a single mode and fiber #F 2 is multimode until this wavelength. Meanwhile, V parameter of the fiber #F 3 calculated at 1.03 μm is 2.54, which indicates this multimode behavior.

Supercontinuum generation in optimized structures
The following discussions are dedicated to study on SC formation in those optimized structured fibers only for the fundamental mode. For this purpose, we solved the GNLSE with use of the split-step Fourier method (Agrawal 2000): where Ã ( ) is the frequency domain of an envelope pulse A(t). The influence of dispersion characteristics ( ) on the nonlinear properties of the fiber is described by a product of the complex spectral envelope Ã ( ) and ̃( ) = ( ) − − 0 1 − 0 , where 0 and 1 are the dispersion coefficients associated with the Taylor series expansion of the propagation constant ( ) about 0 . The term ̃ ( ) is the frequency domain of the attenuation coefficient α (t). In simulations, the attenuation term is approximated by the sum of confinement loss and material loss, n 2 is the nonlinear refractive index of C 2 Cl 4 , while A eff ( ) is the effective mode area.
The response function R t ′ depends on the contribution mechanism e.g. molecular reorientation, molecular interactions, collision-induced polarizability, electronic response, etc. It can be expressed as (Chemnit et al. 2018;Agrawal 2000): where the subscripts el, d, l, c indicates the bound-electronic, molecular reorientation, molecular interaction and collision-induced, respectively, while the coefficients C 2l,2d,2c are normalized constant and N = (n 2el + n 2c + n 2d + n 2l ). n 2,R is the Raman contribution and the term g( ) stands for the distribution function of vibration motion.
The all-normal SC generation in #F 1 was investigated with the input pulses having a central wavelength at 1.56 μm, pulse duration of 90 fs and various energy pulses. With input pulse energy of 1.5 nJ, spectral bandwidth covers from 0.8 to 2.0 μm within 10 dB dynamic range. Figure 8b depicts the evolution of the SC spectrum along fiber #F 1 , and temporal profile at various propagation lengths. At the beginning of propagation length, the spectral broadening is induced by a self-phase modulation (SPM), which is characterized by the temporal spectrogram in S-shape. Due to the effect of dispersion nonlinear, i.e. selfsteepening, the spectrum is asymmetry toward shorter wavelengths. The short wavelengths generated by SPM experience the high amount of normal dispersion and travel slower than the tail component at the center of the pulse. The overlap between the slower and tail components would generate the new wavelengths by four-wave mixing (FWM), which are Fig. 8 SC spectrum: a for various input pulse energy in 5 cm length when pumped with pulses of 90 fs duration, 1.56 μm pump wavelength and the first-order coherence is calculated with 20 individual pairs with random single-photon-per-mode noise seed for 1.5 nJ; b evolution of the SC along fiber and temporal profile at various propagation lengths, for #F 1 attributed to optical wave breaking (OWB). As shown in Fig. 8b, OWB firstly occurs in the trailing edge of the pulse at 1.5 cm of propagation and generates the new wavelength band around 0.9 μm. On the leading edge, OWB occurs only after 2.5 cm of propagation and creates a new wavelength band around 2.0 μm. During further propagation, the spectral broadening at the trailing edge of the pulse is restricted by the high slope of the dispersion profile. At the leading edge of the pulse, the further broadening is also limited due to the large effective mode area, and thus low nonlinear coefficient at the long wavelength. As shown in Fig. 8b, the spectral broadening is completed around 2.5 cm. For further propagation, the energy of the pulse is redistributed from the central area to the edges resulting in the further flat on the wings of the output spectrum and stretching out of group delay trace, as well.
The first-order degree of coherence is calculated with 20 individual pairs of the random noise seeds. Its value is approx. 1 over the spectrum (Fig. 8a). It means that the all-normal dispersion SC generation in #F 1 has the potential for high temporal coherence.
In the case of the fiber #F 2 , the SC generation is analyzed with central wavelength of input pulse at 1.56 μm, and pulse duration of 90 fs. The pump wavelength at the anomalous dispersion regime, soliton dynamics, e.g. soliton fission, self-frequency shifting (Dudley Fig. 9 SC spectrum: a for various input pulse energy in 10 cm length when pumped with pulses of 90 fs duration, 1.56 μm pump wavelength and the first-order coherence is calculated with 20 individual pairs with random single-photon-per-mode noise seed for 1.5 nJ; b evolution of the SC along fiber and temporal profile at various propagation length, for #F 2 et al. 2006) would play a significant role in spectral broadening. The input pulse transforms into higher-order solitons which can fall apart into a series of fundamental solitons (Zhang et al. 2006) due to a small perturbation, such as high-order dispersion, or stimulated Raman scattering (SRS). During further propagation, the solitons would shift toward the longer wavelengths due to soliton self-frequency shifts, causing a considerable spectral expansion toward the side of red light. For the spectral broadening across the ZDWs, the dispersive wave can play a significant role to generate the new wavelength band in the normal dispersion regime (Dudley et al. 2006).
The characteristic lengths of #F 2 fiber are calculated according to equation: where L D , L NL , L fiss define the dispersive, nonlinear, and soliton fission characteristic length scales, respectively. P 0 and t 0 are the peak power and the pulse duration of the input pulse. β 2 represents the group velocity dispersion at the central wavelength of the input laser pulse.
The spectral bandwidth of the SC spectrum extends from 1.0 to 2.0 μm within 30 dB dynamics range (Fig. 9a). The evolution of the SC spectrum in #F 2 is shown in Fig. 9b. At (c) (b) Fig. 10 SC spectrum: a for various input pulse energy in 10 cm length when pumped with pulses of 120 fs duration, 1.03 μm pump wavelength; b the first-order coherence is calculated with 20 individual pairs with random single-photon-per-mode noise seed for 2.0 nJ and 2.5 nJ; c evolution of the SC along fiber and temporal profile at various propagation length for 2.0 nJ, for #F 3 the beginning of the propagation, SPM is the main contribution factor to spectral broadening. Since pump wavelength (1.56 μm) is near the ZDW (1.55 μm), a fraction of the spectrum is induced by SPM across the ZDW and experiences in normal-dispersion regime before soliton fission occurs. Moreover, with the small value of dispersion (β 2 = − 0.0041 ps 2 /m), the characteristic length 5L D , over which the ejected solitons begin separating, is extremely larger than the length of fiber (10 cm). The amplification of input noise caused by modulation instability can be reduced in the normal dispersion and near-zero anomalous dispersion regimes resulting in the improvement of the coherent characteristics. The complex degree of first-order coherence is calculated from 20 individual pairs of spectra generated with the random noise seed. Its value is approx. 1 over the spectrum as shown in Fig. 9a. It indicates that the SC spectrum is high coherence.
Fiber #F 3 has two ZDWs at 1.04 μm (λ D1 ) and 1.82 μm (λ D2 ); the pump wavelength is at normal dispersion regime and near the zero-dispersion wavelength λ D1 . For the case of 2.0 nJ input pulse energy, the spectral bandwidth of the SC spectrum extends from 0.7 to 2.4 μm. At the beginning of propagation, the spectral broadening is induced by SPM with S-shape characteristics of temporal profile. After across λ D1 , the spectrum broadens rapidly due to soliton dynamics, e.g. soliton fission. The red-dispersive wave appears at 2 cm of propagation when the spectrum across zero-dispersion wavelength λ D2 and create a spectrum at wavelengths around 2.4 μm (detail in Fig. 10c). C 2 Cl 4 has a majority contribution for the nonlinear response from the non-instantaneous mechanism, e.g. diffusive reorientation libration. This leads the spectral bandwidth induced by soliton dynamics to be reduced when compared to the one in a media with a purely instantaneous Kerr response. The high fraction of noninstantaneous nonlinear response also suppresses the effect of modulation instability on the SC generation of the anomalous dispersion regime (Maji and Chaudhuri 2014). Therefore, the SC spectrum in #F 3 is high coherence with random single-photon-per-made noise (Fig. 10b). However, unfortunately, this still comes at a cost of the complex temporal structure of the pulse, as can be seen in Fig. 10a. Furthermore, the coherent of SC spectral will decrease when input peak power increasing. For the case of energy, when it is higher than 2.5 nJ, the SC spectrum becomes low coherent since the modulation instability frequency components are not coherent with the pump and fission of each soliton occur randomly (as shown in Fig. 10b).
The comparison between the properties of proposed PCFs and some other liquid-filled core PCFs are presented in Table 3. All the cases, the refractive index of liquids has higher than for fused silica. The obtained results indicated that the proposed fibers can be achieved at the same level of bandwidth as previous works. However, unlike most of other liquid-filled PCFs (except for C 6 H 5 NO 2 -filled PCFs), fiber #F 1 having all normal dispersion can support a singlemode at pump wavelength (1.56 μm). It is worth mentioning that a single-mode fiber can offer the SC generation with further broad spectral bandwidth similarly/like in multimode fiber with the same input peak power. The reason for this phenomenon is that the single-mode fibers allow avoiding the interaction between the modes existing in the fiber, which would transfer light energy from fundamental mode into high-order modes resulting in the reduction of the spectral bandwidth of SC spectrum (Poletti and Horak 2009).

Conclusion
In this paper, we reported numerical studies on three optimum structures of PCFs made of fused silica with a tetrachloroethylene-filled core to obtain both all-normal and anomalous dispersion characteristics. Due to the modifications of their micro-structured geometries, the large scope optimization process of the PCF structure has been carried out. This was done to achieve the flat dispersion, zero dispersion wavelength matched with a pump wavelength for broad SC spectrum. According to the conducted simulations, optimized fiber #F 1 with lattice constant Λ = 1.5 μm and filling factor f = 0.4 exhibited an all-normal dispersion, and its peak equals − 15 ps/nm/km at 1.56 μm. Fiber #F 2 with Λ = 4.0 μm and f = 0.45 has anomalous dispersion for wavelengths longer than 1.55 μm, and the dispersion at 1.56 μm is equal to 0.6 ps/nm/km. Fiber #F 3 with Λ = 1.5 μm, f = 0.55 has two zero dispersion wavelengths at 1.04 μm and 1.82 μm. For this fiber, the dispersion at 1.03 μm is equal to − 4.85 ps/nm/km. Coherent, octave-spanning of the SC spectrum in #F 1 was shown in the wavelength range of 0.8 to 2.0 μm when the input pulse had a pump wavelength at 1.56 μm, 90 fs duration and 1.5 nJ pulse energy. Meanwhile, fiber #F 2 could achieve a high coherence of the SC with pumped wavelength at an anomalous dispersion regime. The SC in #F 2 induced by soliton dynamics and octave-spanning of the SC spectrum (from 1.0 to 2.0 μm) was generated with the same input pulse as fiber #F 1 . The last fiber (#F 3 ) demonstrated an over octave-spanning of the SC spectrum (0.7 ÷ 2.4 μm) with the following initial pulse parameters: pump wavelength at 1.03 μm, and 120 fs duration and 2.0 nJ pulse energy, which correspond to the peak power of 16.66 kW.
Our optimized fibers have certain advantages over other types of optical fibers dedicated to SC. Due to the higher nonlinearity of C 2 Cl 4 than that of fused silica, lower power of input pulses is required rather than in the case of silica PCFs ( Van et al. 2018a;Ferhata et al. 2018). In comparison with other liquid-filled PCFs in previous works (Zhang et al. 2006;Ding et al. 2018; Van et al. 2020;Hoang et al. 2019; Van et al. 2017;Guo et al. 2016), fiber #F 1 with an all-normal dispersion can offer a singlemode guidance at the pump wavelength, but this is not obtained in the cases of other liquids mentioned above, except for C 6 H 5 NO 2 -filled PCFs (Van et al. 2020). However, C 6 H 5 NO 2 -filled PCFs has relatively high attenuation in the near-IR range (Zhang et al. 2006). Meanwhile, fibers #F 2 and #F 3 are also multi-mode, which is the same as C 7 H 8 , CCl 4 filled-core PCFs. Nevertheless, the attenuation of the fibers investigated in this work (around 0.02 dB/cm) is smaller than that of C 7 H 8 -filled (around 0.4 dB/ cm) ( Van et al. 2017) and equivalent to CCl 4 -filled PCFs (around 0.0185 dB/cm) (Ding et al. 2018) at pump wavelengths. Another advantage of our fibers is that C 2 Cl 4 is hightransparent from the visible to the mid-IR range (0.5-3.39 μm) (Takahashi et al. 1985). Therefore, PCFs infiltrated with C 2 Cl 4 can be used to generate SC in the mid-IR range, which has a significant role in biological imaging and optical coherent tomography (Holzwarth et al. 2000;Udem et al. 2002;Tu and Boppart 2013;Israelsen et al. 2019).