Detailed numerical analysis for dispersion properties and losses has been carried out for a new type of Photonic crystal fiber where the air-holes are arranged in a circular pattern with a silica matrix called as Circular Photonic Crystal Fiber (C-PCF). The dependence of different PCF geometrical parameters namely different circular spacings, air-hole diameter and numbers of air-hole rings are carried out in detail towards practical applications. Our numerical analysis establishes that total dispersion is strongly affected by the interplay between material dispersion and waveguide dispersion. For smaller air-filing fraction, adding extra air-hole rings does not change dispersion much whereas for higher air-filling fraction, the dispersion nature changes significantly. With proper adjustment of the parameters ultra-flattened dispersion could be achieved; though the application can be limited by higher losses. However, the ultra-flat dispersion fibers can be used for practical high power applications like supercontinuum generation (SCG) by reducing the loss at the pumping wavelength by increasing the no of air-hole rings. Broadband smooth SCG can also be achieved with low loss oscillating near-zero dispersion fiber with higher no of air-hole rings. The detail study shows that for realistic dispersion engineering we need to be careful for both loss and dispersion.
1. Introduction
Photonic crystal fibers (PCFs) or microstructured optical fibers (MOFs) [1, 2] are special types of optical fiber where air holes are arranged in a periodic nature in the cladding. These types of fibers possess some novel guiding properties, related to the geometric characteristics of the air holes in their cross-section and have been successfully exploited in different applications [1, 2]. Most of the air holes in the PCFs cladding have been arranged either in a periodic triangular or periodic square orientation. The modal properties, in particular, and the dispersion properties of the above types of PCFs can be altered by varying the hole-to-hole spacing (Λ) and the air hole diameter (d) with air-filling fraction being d/Λ [3, 4]. Both types of PCFs with a silica background can be successfully implemented to compensate the positive dispersion parameter and dispersion slope of the existing inline fibers [3, 4]. These fibers can be engineered for designing ultraflattened near-zero dispersion [5–7] or can be engineered to have ultranegative dispersion values near the communication wavelength [8–10].
Recently, there have been new types of cross-sectional geometry where air holes are arranged in a circular pattern [11–16]. Photonic Crystal with air holes arranged in circular layer arrangements has been investigated for different applications [11–15] such as studying thermal properties for electrical driving [11], high transmission waveguides [12], high quality factor microcavity lasers with isotropic photonic band gap effect [13], localization of electromagnetic waves [14] and achieving ultraflat dispersion [15], and so forth. Argyros et al. [16] realized PCF with circular air holes arrangement with polymer background. Though the above works based on C-PCF demonstrate few researches that have been accomplished, still a detailed guiding and modal properties based on the geometrical parameters of C-PCF with silica background have not yet been explored for its potential applications in the communication wavelength window or infrared (IR) or far-IR applications.
In this work, we have made a detailed numerical analysis of the dispersion and loss property of the C-PCF for the first time. We have considered different circular spacings, air hole diameter, and different numbers of air hole rings (Nr). The detailed study will be useful for practical engineering applications where both loss and dispersion play an integral part.
2. Geometry of the Studied Structure and Modal Analysis
The schematic diagram of the C-PCF is shown in Figure 1 with the air hole diameter as d and the radius of the circles, as well as the spacing between the rings, is R. The C-PCF, where air holes are arranged in a circle of certain radius is constructed by repeating the circular unit around the core center. The circles have 6n (with n=1, 2, 3, and so forth for first, second, and third rings of air holes) number of air holes in a particular air hole ring. The angular spacing between any two consecutive air holes in a circle can be given by θ=360/6n, where n=1, 2, 3, and so forth for first, second, and third layers of the circles forming the cladding cross-section. The structure under investigation is having twofold symmetry (C2v symmetry). We compare the propagation characteristics of circular-lattice PCF with triangular-lattice PCF for the structural parameters; the pitch “Λ” (hole-to-hole distance of regular triangular-lattice PCF) is found to have the direct equivalence with “R” (the radius of the circle and the distance between two consecutive circles). Thus, employing the above analogy under the same air hole diameter d, d/Λ (the air-filling fraction) for triangular PCF can be compared with d/R of C-PCF. Another important geometrical parameter is air-filling fraction f. The corresponding value for triangular-lattice PCF is fΔ=(πd2)/(23Λ2)=0.9069(d2/Λ2) [4], whereas for the C-PCF the value comes out to be f°≈(6πd2)/(8R2)≈0.75(d2/R2). So C-PCF is having an air-filling factor which is 82% compared to that of triangular lattice PCF. With the lower index of the air-filled cladding and the higher index of the silica fiber core (the center air hole being missing), the propagation modes can exist due to the modified total internal reflection (TIR) mechanism in this type of PCF structure. The modal fields are calculated using CUDOS MOF utilities [17] that simulates PCFs using the multipole method [18, 19]. The validity and efficiency of the method has already been established [18, 19]. The numerical calculations, namely, dispersion parameter (D) and losses (L), are calculated with MATLAB. The guiding properties of the C-PCF have been studied in the following section by taking into consideration of the various combinations of radius of the circle and air hole diameter. We discuss the modal properties in the subsequent sections as follows.
The schematic diagram of the studied fiber. Air holes are arranged equally spaced in a circular pattern of radius R in silica background.
3. Propagation Characteristics: Results and Analysis
PCFs, with the ability to change their different geometrical parameters, possess the unique and attractive applications like dispersion tailoring [3], endlessly single-mode operation [20], and higher nonlinearity [21] which are not achievable in standard silica fibers. Controllability of waveguide dispersion in PCFs is an important issue for applications of optical communications, dispersion compensation, nonlinear applications, and so forth. The dispersion property for PCFs with regular triangular lattice and square lattice can be controlled by varying the d and Λ [3, 4]. Similarly, for C-PCF, the dispersion can be controlled by varying the radius of the circle R and air hole diameter d and number of air hole rings (Nr). We have calculated the dispersion parameter (D) of the structures through
(1)D=-λcd2Re[neff]dλ2.
Here Re[neff] is the real part of the effective indices obtained from the simulation, and c is the velocity of light in vacuum.
The confinement loss for the structures has been calculated through
(2)L=2πλ20ln(10)106Im(neff)dBm,
where Im(neff) is the imaginary part of the effective indices (obtained from the simulations) and λ in micrometer. The chromatic dispersion of the background silica material has been taken into account through Sellmeier’s equation.
3.1. Dispersion Properties
We studied the dispersion and loss property of the C-PCF for a broad wavelength range from 0.8 μm to 2.8 μm. Figure 2 shows the dispersion property taking into consideration equal air hole diameter (we have taken it to be 0.80 μm), for different values of circular spacing (R) with Nr=3. The figure shows that for smallest values of the radius of the circle (R=1.50μm) the dispersion property is oscillating in the whole wavelength range and few numbers of zero-dispersion wavelengths can be found. As we increase the value of R, the oscillating nature gradually disappears. Instead, a monotonically increasing nature of the dispersion could be found out. This typical nature of the dispersion can be explained on the basis of the interplay between the material dispersion and waveguide effect upon the dispersion. For smaller values of R, the waveguide nature dominates, which results in an oscillating nature of the dispersion. The figure clearly indicates that the smaller the values of the radius of the circle, the larger the magnitude of oscillation. For all cases of R, for smaller wavelength the dispersion values are always negative. This is because of the fact that, for smaller values of wavelength, the material dispersion is so negative that the waveguide dispersion just cannot compensate the material one. This characteristic justifies that for smaller wavelengths all the dispersion tends towards the material one.
(a) Dispersion and (b) losses for different spacings of the circle (R) for three-ring C-PCF with air hole diameter d=0.8μm.
3.1.1. Ultraflat Dispersion and Its Dependence Upon Nr
An interesting observation that follows by observing from Figure 2 is that, for R=2.40 μm, the dispersion value is almost flat (average D value of 13.5 ps/nm/km around 1.65 μm of wavelength) over a long wavelength range (from 1.38 μm to 1.91 μm, i.e., wavelength band of almost 530 nm) with the dispersion oscillation of 0.35 ps/nm/km. For Nr=4,without changing the other design parameters, we could observe that the oscillation increases to 1.75 ps/nm/km (with wavelength band of 880 nm) with the average D=8.25 ps/nm/km around 1.61 μm as shown in Figure 3(a). With the same parameter, the flatness reduces (i.e., the oscillation increases) with the dispersion oscillation of 3.3 ps/nm/km and 4.5 ps/nm/km for Nr=5 and Nr=6 as shown in Figures 4(a) and 5(a), respectively. However, we could achieve flat dispersion with D=9.75 ps/nm/km (with dispersion fluctuation of 0.45 ps/nm/km around 1.63 μm) with Nr=6 and R=2.6μm. Another interesting observation can be found out from Figure 2(a) regarding the dispersion variation with R=2.10μm. The dispersion oscillates around zero dispersion with D=0±11.2 ps/nm/km for a wavelength bandwidth of 1460 nm. Similar nature can be observed with Nr=4 (D=0±9.1 ps/nm/km, wavelength bandwidth of 1380 nm with R=2.20μm), Nr=5 (D=0±8.3 ps/nm/km, wavelength bandwidth of 1300 nm with R=2.30μm), and Nr=6 (D=0±8.2 ps/nm/km, wavelength bandwidth of 1340 nm with R=2.30μm). This type of dispersion properties (oscillating near zero ultraflat dispersion PCF) is becoming an excellent candidate for broadband smooth supercontinuum generation (SCG) [22, 23] where flatness of the spectrum strongly depends on the dispersion nature [24]. Another interesting fact that can be observed from the dispersion figures is the presence of an almost “R independent dispersion”. For Nr=3, the corresponding wavelength comes around 1.32 μm (at least with the R values from 1.70 μm to 2.5 μm) with the crossing region of approximately 0.1 μm. There are similar regions for other cases with Nr=4, Nr=5, and Nr=6 with the approximate value of the center wavelength of 1.306 μm, 1.305 μm, and 1.311 μm, respectively. So with the increment of Nr the “R independent dispersion” wavelength moves towards smaller wavelengths and then shifts to higher wavelength. However, the particular wavelength remains little perturbed with the change of Nr.
(a) Dispersion and (b) losses for different spacings of the circle (R) for four-ring C-PCF with air hole diameter d=0.8μm.
(a) Dispersion and (b) losses for different spacings of the circle (R) for five-ring C-PCF with air hole diameter d=0.8μm.
(a) Dispersion and (b) losses for different spacings of the circle (R) for six-ring C-PCF with air hole diameter d=0.8μm.
3.2. Loss of the Structure and Its Nr Dependence
For the above discussed almost flattened dispersion cases, the corresponding losses around C-band of wavelength are found out to be approximately 1.95×103 dB/m, 5.13×102 dB/m, 1.55×102 dB/m and 4.82×101 dB/m, for Nr=3, Nr=4, Nr=5, and Nr=6 as shown in Figures 2(b), 3(b), 4(b), and 5(b), respectively. So, though we might be getting a flattened dispersion PCF, the design will not be suitable for practical applications as the losses are very high. However, for high power nonlinear applications like SCG and so forth, the pump wavelength is the wavelength of zero dispersion wavelength (ZDW) [24, 25]. For all the above cases, the values of ZDW come out to be 1.120 μm, 1.132 μm, and 1.134 μm and 1.136 μm, respectively. The corresponding losses around these wavelengths are 3.60×102 dB/m, 6.25×101 dB/m, 1.06×101 dB/m, and 1.81 dB/m, respectively. As high power nonlinear applications like SCG require even less than a meter long of the fiber [25], practical applications of the designed fiber are feasible especially with Nr=6 and with R=2.6μm as shown in Figure 5. The losses corresponding to first ZDW of the oscillating near-zero flat dispersion comes out to be 1.53×102 dB/m, 2.77×101 dB/m, 0.64×101 dB/m, and 0.94 dB/m for Nr=3, Nr=4, Nr=5 and Nr=6 respectively. The above data shows that applications like broadband smooth SCG (where only a meter long of the fiber is required) are possible with Nr=6 and with R=2.3μm. The above results imply that the losses can be reduced by increasing Nr, but the dispersion nature will not remain the same. So to design a proper practical design we need to be careful for both the loss and dispersion. Similar to the case of R independent dispersion, another interesting fact that can be observed is the presence of “R independent loss” which is more prominent. The wavelength for “R independent loss” with 3 air hole rings comes out to be 1.79 μm, and the range for this wavelength is even less than 0.02 μm. Similar phenomenon can be observed for other values of air hole rings as well as with the wavelengths of 1.93 μm, 2.01 μm, and 2.07 μm for Nr=4, Nr=5 and Nr=6 as can be observed from Figures 3(b), 4(b), and 5(b) respectively. Here the “R independent loss” wavelength gets red-shifted with the increase of Nr. Throughout the above designs we have considered the confinement losses of the structure; however, the sources of losses in MOFs can be limited by some other factors like structural imperfections, Rayleigh scattering, and absorption which are not affected by the geometrical losses.
3.3. Effect of Air Hole Diameter (d)
In this section we will be discussing the effect of change of air hole diameter upon dispersion for a fixed value of circular spacing. Figure 5 shows the dispersion and loss characteristics for different values of air hole diameter keeping R (=1.50 μm), and Nr (=3) constant. Figure 6(a) dictates that the oscillation amplitudes in dispersion curves increase with air hole diameter. For d=0.90μm, the oscillation amplitude becomes 250 ps/nm/km, whereas, for the smaller values of d=0.60μm the amplitude of oscillation reduces to 30 ps/nm/km. These typical characteristics can be explained by considering the interplay between material dispersion and waveguide dispersion. For higher values of air hole diameter the waveguide dispersion dominates, whereas for smaller values of air hole diameter, material dispersion takes over the total dispersion. It is interesting to note that the value of Dmax for a particular circular radius increases with the increase of air hole diameter. This is an important finding as the dispersion can be tailored according to the engineer’s requirement for specific applications. The losses for the above-studied structures are presented in Figure 6(b). The loss gets reduced with the increase of air hole diameter as the mode becomes tightly confined.
(a) Dispersion and (b) losses for a three-ring MOF for different air hole diameters d.
3.4. Effect of Number of Air Hole Rings (Nr)
In the following section, we discuss the effect of change of Nr. For this purpose we have subdivided the analysis into two sections. In the first case, we discuss the effect for smaller values of air-filling fraction such that material dispersion dominates the total dispersion and for this obvious reason we will not be having any local maxima in our dispersion curve. The effect has been shown in Figure 7. The dispersion decreases as Nr is increased as clearly shown in Figure 7(a). The difference between successive dispersion curves of two MOFs decreases as Nr increases, as shown in Figure 7(c). This figure clearly shows that the dispersion almost converges, when Nr is increased. This convergence depends upon the operating wavelength. The larger the wavelength, the higher the convergence, as shown in Figure 7(c). The graph clearly indicates that for λ=1.55μm the limit is almost achieved, whereas for smaller wavelength (λ=1.20μm) the convergence is yet to reach its limit even with Nr=8. This factor is different to that of regular triangular PCF [3]. The losses with higher Nr, say for 8 rings, are still higher as shown in Figure 7(b).
(a) Dispersion and (b) losses for different Nr for R=2.0μm and d=0.5μm. (c) Convergence for three different wavelengths for different numbers of air hole rings (Nr) for the above parameters.
We now switch to the second case, where we have considered higher air-filing fraction so that we are having oscillating nature in the dispersion characteristics. For this purpose, we have considered R=1.55μm with d=0.6μm. With the increase of Nr, the oscillation amplitude increases and the negative dispersion value also increases as shown in Figure 8(a). This is due to the increasing effect of waveguide dispersion in the total dispersion. The losses for the above cases are presented in Figure 8(b), which clearly shows that the performance respect losses improve with the increase of air hole rings. The above influence can be explained on the basis of the mode confinement. When the mode is well confined to the core, an additional air hole ring is not going to change the modal properties drastically, whereas for the reverse case for a loosely bound mode an additional air hole ring is going to affect its modal properties significantly and we can expect the field pattern associated with the mode to converges with the increase of air hole ring.
(a) Dispersion and (b) losses for different numbers of air hole rings. R=1.55μm and d=0.60μm.
4. Conclusion and Discussion
We carried out a detail numerical analysis of the dispersion and loss characteristics of a new type of photonic crystal fiber where air holes are arranged in a circular pattern. Influence of different geometrical parameters, namely, the circular spacing, air hole diameter, and number of air hole rings has been carried out in detail. Total dispersion is strongly affected by the interplay between material dispersion and waveguide dispersion. For smaller air-filing fraction, adding extra air hole rings does not change dispersion much whereas, for higher air-filling fraction, the dispersion nature changes significantly. With proper adjustment of the available parameters, ultraflat dispersion (dispersion fluctuation of 0.35 ps/nm/km around 1.65 μm) for a broad wavelength band could be achieved, though the practical applications around C-band are limited due to the higher loss. However, for high power applications like SCG and so forth where less than a meter of the fiber is sufficient, the ultraflat dispersion fiber will be suitable with the low loss at the first ZDW. We also achieved small oscillating (D=0±8.2 ps/nm/km) near-zero dispersion, and the design will be helpful for broadband smooth SCG (where only a meter long of the fiber is required) as the loss at the first ZDW comes out to be 0.95 dB/m. It has been shown that, for increased number of air holes, the dispersion tends to saturate. The saturation is sensitive for different wavelengths, and the saturation improves at higher wavelengths. The detailed study shows that for realistic dispersion engineering, we need to be careful for both the loss and dispersion.
Conflict of Interests
The authors would like to declare that there is no direct financial relation with any commercial identity mentioned in their work that might lead to a conflict of interests.
Acknowledgments
The authors would like to thank Dr. Boris Kuhlmey, University of Sydney, Australia, for providing valuable suggestions in understanding the software for designing and studying the properties of different structures. The authors acknowledge sincerely the Defence Research and Development Organization, Government of India, and CRF of IIT Kharagpur for the financial support to carry out this research.
KnightJ. C.Photonic crystal fibers2003424847851RussellP. S. J.Photonic-crystal fibers20062412472947492-s2.0-3384717063210.1109/JLT.2006.885258KuhlmeyB. T.RenversezG.MaystreD.Chromatic dispersion and losses of microstructured optical fibers20034246346392-s2.0-0037330121BoukA. H.CucinottaA.PoliF.SelleriS.Dispersion properties of square-lattice photonic crystal fibers20041259419462-s2.0-294270987810.1364/OPEX.12.000941ReevesW. H.KnightJ. C.RussellP. S. J.RobertsP. J.Demonstration of ultra-flattened dispersion in photonic crystal fibers200210146096132-s2.0-0001086925SaitohK.KoshibaM.HasegawaT.SasaokaE.Chromatic dispersion control in photonic crystal fibers: application to ultra-flattened dispersion20031188438522-s2.0-0141945705MajiP. S.ChaudhuriP. R.A new design of ultra-flattened near-zero dispersion PCF using selectively liquid infiltration201322532ZhaoX.ZhouG.ShuguangL.LiuZ.WeiD.HouZ.HouL.Photonic crystal fiber for dispersion compensation20084728519051962-s2.0-6074911606910.1364/AO.47.005190YangS.ZhangY.PengX.LuY.XieS.LiJ.ChenW.JiangZ.PengJ.LiH.Theoretical study and experimental fabrication of high negative dispersion photonic crystal fiber with large area mode field2006147301530232-s2.0-3364553191510.1364/OE.14.003015GérômeF.AugusteJ. L.BlondyJ. M.Design of dispersion-compensating fibers based on a dual-concentric-core photonic crystal fiber20042923272527272-s2.0-1064428003110.1364/OL.29.002725HoW. D.LuT. W.HsiaoY. H.LeeP. T.Thermal properties of 12-fold quasi-photonic crystal microcavity laser with size-controlled nano-post for electrical driving20092723530253072-s2.0-7035051354310.1109/JLT.2009.2031501HoriuchiN.SegawaY.NozokidoT.MizunoK.MiyazakiH.High-transmission waveguide with a small radius of curvature at a bend fabricated by use of a circular photonic crystal20053099739752-s2.0-1994437912510.1364/OL.30.000973LeeP. T.LuT. W.FanJ. H.TsaiF. M.High quality factor microcavity lasers realized by circular photonic crystal with isotropic photonic band gap effect200790152-s2.0-3424723327410.1063/1.2724899151125ZhongW.ZhangX.Localized modes in defect-free two-dimensional circular photonic crystals201081162-s2.0-7364908900710.1103/PhysRevA.81.013805013805XuH.WuJ.XuK.DaiY.XuC.LinJ.Ultra-flattened chromatic dispersion control for circular photonic crystal fibers20111352-s2.0-7995582640310.1088/2040-8978/13/5/055405055405ArgyrosA.BassettI. M.van EijkelenborgM. A.LargeM. C. J.ZagariJ.NicoroviciN. A. P.McPhedranR. C.de SterkeC. M.Ring structures in microstructured polymer optical fibres20019138138202-s2.0-0011053599CUDOS MOF Utilitieshttp://sydney.edu.au/science/physics/cudos/research/mofsoftware.shtmlWhiteT. P.KuhlmeyB. T.McPhedranR. C.MaystreD.RenversezG.de SterkeC. M.BottenL. C.Multipole method for microstructured optical fibers. I. Formulation20021910232223302-s2.0-0000992754KuhlmeyB. T.WhiteT. P.RenversezG.MaystreD.BottenL. C.de SterkeC. M.McPhedranR. C.Multipole method for microstructured optical fibers. II. Implementation and results20021910233123402-s2.0-0000992755BirksT. A.KnightJ. C.RussellP. S. J.Endlessly single-mode photonic crystal fiber199722139619632-s2.0-0031185826DudleyJ. M.TaylorJ. R.Ten years of nonlinear optics in photonic crystal fibre20093285902-s2.0-5934909036510.1038/nphoton.2008.285SaitohK.KoshibaM.Highly nonlinear dispersion-flattened photonic crystal fibers for supercontinuum generation in a telecommunication window20041210202720322-s2.0-294272462010.1364/OPEX.12.002027BegumF.NamihiraY.KinjoT.KaijageS.Supercontinuum generation in square photonic crystal fiber with nearly zero ultra-flattened chromatic dispersion and fabrication tolerance analysis201128449659702-s2.0-7864973040010.1016/j.optcom.2010.10.029DudleyJ. M.GentyG.CoenS.Supercontinuum generation in photonic crystal fiber2006784113511842-s2.0-3374946499510.1103/RevModPhys.78.1135RoyS.ChaudhuriP. R.Supercontinuum generation in visible to mid-infrared region in square-lattice photonic crystal fiber made from highly nonlinear glasses200928217344834552-s2.0-6765041047410.1016/j.optcom.2009.05.062