* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download The Field and Energy Distributions of the Fundamental Mode in the
Survey
Document related concepts
Transcript
Hindawi Publishing Corporation ISRN Optics Volume 2013, Article ID 567501, 8 pages http://dx.doi.org/10.1155/2013/567501 Research Article The Field and Energy Distributions of the Fundamental Mode in the Solid-Core Photonic Crystal Fibers for Different Geometric Parameters Halime Demir Inci and Sedat Ozsoy Department of Physics, Faculty of Science and Arts, Erciyes University, 38039 Kayseri, Turkey Correspondence should be addressed to Sedat Ozsoy; [email protected] Received 16 June 2013; Accepted 12 July 2013 Academic Editors: Y. S. Kivshar and S. Ponomarenko Copyright © 2013 H. Demir Inci and S. Ozsoy. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. For a solid core photonic crystal fiber with the triangular lattice, the field and energy distributions of fundamental mode are considered at 1.55 πm wavelength. The silica core is constituted by removing the 7 air holes. The cladding consists of the twodimensional silica-air photonic crystal with the 4 rings of air holes. The field and energy distributions were investigated for three different values of d/Ξ. Here, d and Ξ represent the diameter of air holes and the pitch length, respectively. The simulations show that, for the fixed Ξ, the increase in d/Ξ ratio does not cause the considerable changes in the field and energy intensity distributions, but for the fixed d, the increase in this ratio affects the intensity distributions reasonably. 1. Introduction In recent years, the photonic crystal fibers (PCFs) have been of significant interest due to their unique structures and new properties [1β6]. Generally, photonic crystal fibers consist of an arrangement of air holes in the cladding extending the whole length of the fiber. Photonic crystal fibers are categorized into two groups according to light guiding mechanism. One is the index guiding photonic crystal fiber, and the other is the photonic band gap PCF. In the index guiding PCFs, the core region is solid and the light is confined in the central core as in the conventional fibers. The photonic crystal fiber consists of the pure silica fiber with an array of the air holes along the length of the fiber. The core is constituted by removing the central hole from the structure. The higher effective refractive index of the surrounding holes forms cladding in which leading the index guiding mechanism analogous to total internal reflection. Consequently, the light guiding can be explained by the total internal reflection which is also the way light is guided in step index fibers. If around the central air hole there is the two-dimensional photonic crystal cladding which consists of a periodic array of the silica and air, such fibers are called the photonic band gap photonic crystal fibers. The central hole in these fibers acts as the core, and in this place the light is guided by the photonic band gap effect. The frequencies within the band gap of the structure will expose to the multiple Bragg reflection that leads to the destructive interference of the light trying to propagate away from the air core. The function of the air core is to provide a defect at the periodic structure that the propagation of the frequencies within the photonic band gap is really allowed. Therefore, it is not needed for the core indexes of these fibers to be bigger than the effective index of cladding. PCFs have been shown to possess many important properties as the single mode operation over wide range of wavelength, the highly tunable dispersion, the propagation of high power densities without exciting unwanted nonlinear effects, and the high birefringence. These properties have the practical importance in the design of sophisticated broadband optical telecommunication networks [7] and active sensor systems [8]. In this work, the software based on the finite element method (FEM) is used for the field and energy distribution simulations of the fundamental mode at the solid core photonic crystal fiber structures with the different geometrical parameters. 2 ISRN Optics 2. Theoretical Model The large index contrast and complex structure in the photonic crystal fibers make them difficult to be treated mathematically. The standard optical fiber analysis does not help, and moreover in the majority of photonic crystal fibers, it is practically impossible to perform modal analysis analytically; therefore the Maxwell equations must be solved numerically. The main idea consists of transforming this complicated problem, which can be described by the curl-curl equation [9]: β × (ππβ1 β × πΈ)β β π02 ππ πΈβ = 0. d (1) This is an eigenvalue equation, the eigenvalues are π½/π0 the effective indexes, and the eigenvectors are the electric field components (πΈπ₯ , πΈπ¦ , πΈπ§ ). In (1), πΈβ is the electric field, π0 is the wavenumber in vacuum, and ππ and ππ are the dielectric permittivity and magnetic permeability tensors, respectively. Here in the analysis of the properties of PCFs, the FEM is used. The FEM allows the photonic crystal fiber crosssection in the π₯-π¦ plane to be divided into a patchwork of triangular elements in which they can possess different sizes, shapes, and refractive indexes. In this way, any geometry including the PCF air hole and the medium characteristics can be accurately described. In particular, the FEM is suited for studying fibers with the nonperiodic arrangements of the air holes. However, it provides a full vectorial analysis that is necessary to model the PCFs with the large air holes and high index variations and to accurately predict the properties of these fibers [10]. The FEM is basically divided in to four steps [11]. The domain discretization is the first and perhaps the most important step in any finite element method. The first step consists of dividing the solution domain Ξ© into the subdomains of finite number in which form a patchwork of fundamental elements that can possess different sizes, shapes, and physical properties. The second step is choosing the interpolation functions that provide an approximation of the unknown solution within each element. In this specific case, the solution domain is the transverse cross-section of optical waveguide that is divided in triangular finite elements, and the nodal approach is followed using secondorder polynomials as interpolating functions in each finite element (each triangular is characterized by 6 nodes, and the unknowns are the magnetic field components). In the next step, the curl-curl equation is transformed into a generalized eigenvalue problem by applying the variational formulation as follows: [π΄] {Ξ¦} β π [π΅] {Ξ¦} = {0} , Ξ (2) where the eigenvalue π is the mode effective refractive index (πeff ) and the eigenvector Ξ¦ is the full vectorial electric distribution (πΈπ₯ , πΈπ¦ , πΈπ§ ) in the nodal points. Finally, the fourth and final step in the finite element analysis is to solve the eigenvalue system. The matrices [π΄] and [π΅] are sparse and symmetric; therefore, the computational time can be effectively minimized by using the sparse matrix solver. Figure 1: The cross-section of the solid core PCF considered. Ξ is the pitch length, and π is the diameter of air holes. Table 1: The effective refractive indexes of the fundamental mode for different π/Ξ ratios. π/Ξ Fixed pitch (Ξ = 4.2 πm) 0.08 0.09 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 Fixed diameter (π = 0.84 πm) 0.08 0.09 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 π/Ξ πeff (π = 1.55 πm, πsilica = 1.44409) 0.37 1.442911 β 2.534351π β 8π 1.442822 β 1.378536π β 8π 1.442747 β 5.387208π β 9π 1.442334 β 4.984438π β 12π 1.442108 β 2.353948π β 14π 1.441914 β 8.43704π β 17π 1.441721 β 2.22923π β 19π 1.44152 β 3.398141π β 19π 1.441303 β 5.284044π β 19π 1.441068 β 4.492512π β 20π 1.440809 β 3.595418π β 19π 0.15 0.17 0.19 0.37 0.56 0.71 0.91 1.11 1.25 1.48 1.67 1.443843 β 7.982422π β 11π 1.443773 β 9.149398π β 11π 1.443693 β 4.61459π β 11π 1.442334 β 4.984438π β 12π 1.439848 β 1.198296π β 12π 1.436078 β 1.976832π β 13π 1.43085 β 8.180793π β 14π 1.423967 β 3.413292π β 14π 1.41519 β 2.279162π β 14π 1.404232 β 2.485407π β 14π 1.390737 β 3.793315π β 14π 3. Simulation Results For the simulations, the photonic crystal fiber with crosssection shown in Figure 1 is investigated at 1.55 πm wavelength. Here Ξ is the pitch length, and π is the diameter of air holes. The fiber core is silica, and it is formed by removing ISRN Optics 3 Max: 2.103e β 4 ×10β4 2 Max: 0.915 Max: 233.945 200 0.8 150 0.6 100 0.4 50 0.2 0 0 β50 β0.2 1 0.8 β100 β0.4 0.6 β150 β0.6 0.4 β200 β0.8 0.2 Min: β234.822 (a) 1.8 1.6 1.4 1.2 Min: β0.917 (b) Min: 1.715e β 13 (c) Figure 2: For Ξ = 4.2 πm and π/Ξ = 0.2, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. Max: 4.995e β 5 ×10β5 Max: 0.578 Max: 149.414 4.5 100 0.4 50 0.2 0 0 β50 β0.2 1.5 β100 β0.4 1 0.5 4 3.5 3 2.5 2 Min: β149.38 (a) Min: 3.07e β 8 Min: β0.579 (b) (c) Figure 3: For Ξ = 4.2 πm and π/Ξ = 0.08, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. Max: 202.371 200 150 0.6 8 100 0.4 7 50 0.2 6 0 0 β50 β0.2 β100 β0.4 2 β0.6 1 β150 β200 Min: β201.783 (a) Max: 9.003 e β 5 ×10β5 9 Max: 0.770 5 4 3 Min: 4.689e β 11 Min: β0.772 (b) (c) Figure 4: For Ξ = 4.2 πm and π/Ξ = 0.1, (a) the distribution of z-component of the electric field, (b) the distribution of z-component of the magnetic field, and (c) the total energy distribution. 4 ISRN Optics 200 150 0.6 100 0.4 50 0.2 0 0 β50 β0.2 β100 β0.4 β150 β0.6 β200 2 1.5 1 0.5 β0.8 Min: β0.879 Min: β224.829 (a) Max: 2.444e β 4 ×10β4 Max: 0.879 0.8 Max: 225.835 (b) Min: 5.523 e β 16 (c) Figure 5: For Ξ = 4.2 πm and π/Ξ = 0.3, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. Max: 0.814 0.8 Max: 221.453 200 0.6 150 100 0.4 50 0.2 0 0 β50 β0.2 β100 β0.4 β150 β0.6 β200 2.5 2 1.5 1 0.5 β0.8 Min: β0.814 Min: β222.651 (a) Max: 2.788e β 4 ×10β4 (b) Min: 1.143e β 26 (c) Figure 6: For Ξ = 4.2 πm and π/Ξ = 0.6, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. Max: 217.578 200 Max: 2.607 e β 4 ×10β4 2.5 Max: 0.792 150 0.6 100 0.4 50 0.2 0 0 β50 β0.2 β100 2 1.5 1 β0.4 β150 0.5 β0.6 β200 Min: β0.792 Min: β217.289 (a) (b) Min: 2.713e β 39 (c) Figure 7: For Ξ = 4.2 πm and π/Ξ = 0.9, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. ISRN Optics 5 Max: 1.823 e β 4 ×10β4 1.8 Max: 0.609 0.6 Max: 155.133 150 1.6 100 0.4 50 0.2 1.2 0 0 1 β50 β0.2 β100 β0.4 1.4 0.8 0.6 0.4 0.2 β150 Min: β154.856 (a) β0.6 Min: β0.607 (b) Min: 1.555e β 12 (c) Figure 8: For π = 0.84 and π/Ξ = 0.08 πm, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. 150 0.6 100 0.4 50 0.2 0 0 β50 β0.2 β100 β0.4 β150 β0.6 Min: β169.009 (a) Max: 1.888e β 4 ×10β4 1.8 Max: 0.668 Max: 168.202 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 Min: 6.405e β 13 Min: β0.665 (b) (c) Figure 9: For π = 0.84 and π/Ξ = 0.1 πm, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. Max: 1.141 Max: 284.171 Max: 2.166e β 4 ×10β4 2 1.8 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 1 200 0.5 100 0 0 β100 β0.5 β200 β1 Min: β285.155 (a) Min: 2.611 e β 13 Min: β1.141 (b) (c) Figure 10: For π = 0.84 πm and π/Ξ = 0.3 πm, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. 6 ISRN Optics Max: 365.285 Max: 1.921e β 4 ×10β4 1.8 Max: 1.484 300 1 200 100 0.5 0 0 1.6 1.4 1.2 1 0.8 β100 β0.5 β200 0.6 0.4 β1 0.2 β300 Min: β364.21 Min: β1.481 (a) Min: 5.172e β 16 (b) (c) Figure 11: For π = 0.84 πm and π/Ξ = 0.6 πm, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. Max: 1.578e β 4 ×10β4 Max: 1.967 Max: 610.695 600 1.5 1.4 1 1.2 200 0.5 1 0 0 0.8 400 β200 β400 β600 Min: β610.153 β0.5 0.6 β1 0.4 β1.5 0.2 Min: β1.969 (a) Min: 5.851e β 17 (b) (c) 1.2 1.05 1 1 0.95 0.8 Normalized w Normalized w Figure 12: For π = 0.84 πm and π/Ξ = 0.9 πm, (a) the distribution of π§-component of the electric field, (b) the distribution of π§-component of the magnetic field, and (c) the total energy distribution. 0.6 0.4 0.9 0.85 0.8 0.2 0.75 0 0.7 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 d/Ξ (a) 0.6 0.8 1 d/Ξ (b) Figure 13: The normalized total energy of the fundamental mode at the central of the core for different π/Ξ ratios for (a) fixed pitch Ξ = 4.2 πm and (b) fixed diameter π = 0.84 πm. ISRN Optics 7 Normalized z-component of the Poynting ο¬ux 1.2 Normalized z-component of the Poynting flux 1 0.8 0.6 0.4 0.2 0 0 0.2 0.4 0.6 0.8 1.05 1 0.95 0.9 0.85 0.8 0.75 0.7 1 0 0.2 0.4 0.6 0.8 1 d/Ξ d/Ξ (a) (b) Figure 14: The normalized z-component of the Poynting flux of the fundamental mode at the central of the core for different π/Ξ ratios for (a) fixed pitch Ξ = 4.2 πm and (b) fixed diameter π = 0.84 πm. 1.4435 1.45 1.443 1.44 1.43 1.4425 neff neff 1.42 1.442 1.41 1.4415 1.4 1.441 1.39 1.38 1.4405 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 d/Ξ (a) 0.6 0.8 1 d/Ξ (b) Figure 15: The effective refractive indexes of the fundamental mode for different π/Ξ ratios for (a) fixed pitch Ξ = 4.2 πm and (b) fixed diameter π = 0.84 πm. 7 air holes from the structure. The cladding is the twodimensional photonic crystal with 4 rings of the triangular lattice air holes in the silica matrix. Firstly, for the fixed pitch length Ξ = 4.2 πm, the field and energy distributions are investigated by varying the diameters of air holes for the π/Ξ values of 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9. Later, a similar investigation is executed for the fixed diameter of air-hole π = 0.84 πm, by varying the pitch length. For Ξ = 4.2 πm and the cross-section with the radius of 23.1 πm, the electric and magnetic fields and the energy distributions of the structure for π/Ξ = 0.2, 0.08, 0.1, 0.3, 0.6, and 0.9 are shown in Figures 2, 3, 4, 5, 6, and 7, respectively. For a given pitch length, Ξ = 4.2 πm, the electric and magnetic field distributions are more extended into the cladding region (i.e., guiding is weakened) when the π/Ξ ratio reduces (Figures 2β9). In other words, the light is well confined in the core when the π/Ξ ratio increases. For π/Ξ = 0.08, the effective index of the fundamental mode becomes complex, and thus the mode becomes leaky. As a result of this, the energy intensity transported in the core region is also reduced. The reason of this behavior is that the effective cladding index approaches the core index, and hence the index difference decreases. For a fixed hole diameter, π = 0.84 πm, Figures 9 and 10 show the fields and energy distributions of the structure for π/Ξ = 0.1 and the cross-section with the radius of 46.2 πm and π/Ξ = 0.3 and the cross-section with the radius of 15.4 πm, respectively. In addition, Figures 11 and 12 show the fields and energy distributions of the structure for π/Ξ = 0.6 and the cross-section with the radius of 7.7 πm and π/Ξ = 0.9 and the cross-section with the radius of 5.13 πm, respectively. Figures 2, 8, and 9 show that, for a given diameter of the air holes, the confinement and the total energy confined in 8 the core region increase when the π/Ξ ratio increases up to a certain π/Ξ value, for which the effective index becomes complex (see Table 1). However, when the π/Ξ ratio increases above this value, it can be seen that the confinement is lower and the energy involved in the core region decreases (see Figures 2 and 10). This behaviour of the energy can be seen more clearly from Figure 13 and attributed to the mode becoming leaky due to approaching its cutoff. In addition, π§-component of the Poynting flux of the fundamental mode versus π/Ξ values is illustrated in Figure 14. For the different π/Ξ ratios, the effective index variation of the fundamental mode is given in Table 1. The effective indexes are obtained by using cylindrical PML condition. Figure 15 shows the real part of the effective refractive index of the fundamental mode versus π/Ξ ratios for both fixed pitch and fixed diameter. 4. Conclusion In this study, the effect of the diameters of air holes and the pitch length on the field and energy distributions of the fundamental mode in a solid core PCF with the triangular lattice of air holes in the silica matrix is investigated. For a given pitch length, the confinement and the energy density for the fundamental mode decrease when π/Ξ decreases. For a given diameter of the air holes, the confinement and the total energy confined in the core region increase when the π/Ξ ratio increases up to a certain π/Ξ value, for which the effective index becomes complex. However, when the π/Ξ ratio increases above this value, the confinement is lower and the energy involved in the core region decreases. This behavior of the energy can be attributed to the mode becoming leaky due to approaching its cutoff. References [1] J. C. Knight, T. A. Birks, P. S. J. Russell, and D. M. Atkin, βAllsilica single-mode optical fiber with photonic crystal cladding,β Optics Letters, vol. 21, no. 19, pp. 1547β1549, 1996. [2] J. C. Knight, βPhotonic crystal fibres,β Nature, vol. 424, no. 6950, pp. 847β851, 2003. [3] S. Arismar Cerqueira Jr., F. Luan, C. M. B. Cordeiro, A. K. George, and J. C. Knight, βHybrid photonic crystal fiber,β Optics Express, vol. 14, no. 2, pp. 926β931, 2006. [4] T. A. Birks, J. C. Knight, and P. S. J. Russell, βEndlessly singlemode photonic crystal fiber,β Optics Letters, vol. 22, no. 13, pp. 961β963, 1997. [5] A. Ortigosa-Blanch, J. C. Knight, W. J. Wadsworth et al., βHighly birefringent photonic crystal fibers,β Optics Letters, vol. 25, no. 18, pp. 1325β1327, 2000. [6] W. J. Wadsworth, A. Ortigosa-Blanch, J. C. Knight, T. A. Birks, T.-P. M. Man, and P. S. J. Russell, βSupercontinuum generation in photonic crystal fibers and optical fiber tapers: a novel light source,β Journal of the Optical Society of America B, vol. 19, no. 9, pp. 2148β2155, 2002. [7] M. D. Nielsen, J. R. Folkenberg, N. A. Mortensen, and A. Bjarklev, βBandwidth comparison of photonic crystal fibers and conventional single-mode fibers,β Optics Express, vol. 12, no. 3, pp. 430β435, 2004. ISRN Optics [8] S. O. Konorov, A. M. Zheltikov, and M. Scalora, βPhotoniccrystal fiber as a multifunctional optical sensor and sample collector,β Optics Express, vol. 13, no. 9, pp. 3454β3459, 2005. [9] K. Saitoh and M. Koshiba, βSingle-polarization single-mode photonic crystal fibers,β IEEE Photonics Technology Letters, vol. 15, no. 10, pp. 1384β1386, 2003. [10] T. M. Monro, D. J. Richardson, N. G. R. Broderick, and P. J. Bennett, βModeling large air fraction holey optical fibers,β Journal of Lightwave Technology, vol. 18, no. 1, pp. 50β56, 2000. [11] J. Jin, The Finite Element Method in Electromagnetics, Wiley, New York, NY, USA, 2nd edition, 2002. Hindawi Publishing Corporation http://www.hindawi.com The Scientific World Journal Journal of Journal of Gravity Photonics Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Journal of Advances in Condensed Matter Physics Soft Matter Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Journal of Aerodynamics Journal of Fluids Hindawi Publishing Corporation http://www.hindawi.com Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Volume 2014 International Journal of Optics Submit your manuscripts at http://www.hindawi.com International Journal of Statistical Mechanics Hindawi Publishing Corporation http://www.hindawi.com Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Volume 2014 Journal of Thermodynamics Journal of ComputationalβMethodsβinβPhysics Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Journal of Solid State Physics Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Physics Research International Advances in High Energy Physics Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Volume 2014 Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Journal of Astrophysics Biophysics Hindawi Publishing Corporation http://www.hindawi.com Atomic and Molecular Physics Journal of Journal of Hindawi Publishing Corporation http://www.hindawi.com Hindawi Publishing Corporation http://www.hindawi.com Advances in Astronomy Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com International Journal of Superconductivity Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014