# The Field and Energy Distributions of the Fundamental Mode in the Solid-Core Photonic Crystal Fibers for Different Geometric Parameters

Authors: Halime Demir Inci; Sedat Ozsoy
Journal: ISRN Optics
 (2013)

---

## Abstract

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 two-dimensional 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.


---

## Body

## 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. 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)∇×(μr-1∇×E→)-k02εrE→=0.This is an eigenvalue equation, the eigenvalues areβ/k0 the effective indexes, and the eigenvectors are the electric field components (Ex,Ey,Ez). In (1), E→ is the electric field, k0 is the wavenumber in vacuum, and εr and μr 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 cross-section in thex-y 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 second-order 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:
(2)[A]{Φ}-λ[B]{Φ}={0},
where the eigenvalue λ is the mode effective refractive index (neff) and the eigenvector Φ is the full vectorial electric distribution (Ex,Ey,Ez) in the nodal points.Finally, the fourth and final step in the finite element analysis is to solve the eigenvalue system. The matrices[A] and [B] are sparse and symmetric; therefore, the computational time can be effectively minimized by using the sparse matrix solver.


## 3. Simulation Results

For the simulations, the photonic crystal fiber with cross-section shown in Figure1 is investigated at 1.55 μm wavelength. Here Λ is the pitch length, and d is the diameter of air holes. The fiber core is silica, and it is formed by removing 7 air holes from the structure. The cladding is the two-dimensional photonic crystal with 4 rings of the triangular lattice air holes in the silica matrix.Figure 1
The cross-section of the solid core PCF considered.Λ is the pitch length, and d is the diameter of air holes.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 d/Λ 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 d=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 d/Λ = 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Λ=4.2 μm and d/Λ=0.2, (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.
(a)
(b)
(c)ForΛ=4.2 μm and d/Λ=0.08, (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.
(a)
(b)
(c)ForΛ=4.2μm and d/Λ=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.
(a)
(b)
(c)ForΛ=4.2 μm and d/Λ=0.3, (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.
(a)
(b)
(c)ForΛ=4.2 μm and d/Λ=0.6, (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.
(a)
(b)
(c)ForΛ=4.2 μm and d/Λ=0.9, (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.
(a)
(b)
(c)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 d/Λ ratio reduces (Figures 2–9). In other words, the light is well confined in the core when the d/Λ ratio increases. For d/Λ=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.Ford=0.84 and d/Λ=0.08 μm, (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.
(a)
(b)
(c)Ford=0.84 and d/Λ=0.1 μm, (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.
(a)
(b)
(c)For a fixed hole diameter,d=0.84 μm, Figures 9 and 10 show the fields and energy distributions of the structure for d/Λ=0.1 and the cross-section with the radius of 46.2 μm and d/Λ=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 d/Λ=0.6 and the cross-section with the radius of 7.7 μm and d/Λ=0.9 and the cross-section with the radius of 5.13 μm, respectively.Ford=0.84 μm and d/Λ=0.3 μm, (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.
(a)
(b)
(c)Ford=0.84 μm and d/Λ=0.6 μm, (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.
(a)
(b)
(c)Ford=0.84 μm and d/Λ=0.9μm, (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.
(a)
(b)
(c)Figures2, 8, and 9 show that, for a given diameter of the air holes, the confinement and the total energy confined in the core region increase when the d/Λ ratio increases up to a certain d/Λ value, for which the effective index becomes complex (see Table 1). However, when the d/Λ 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, z-component of the Poynting flux of the fundamental mode versus d/Λ values is illustrated in Figure 14.Table 1
The effective refractive indexes of the fundamental mode for differentd/Λ ratios.
d
/
Λ
λ
/
Λ
n
eff (λ = 1.55 µm, nsilica = 1.44409)
Fixed pitch(Λ = 4.2 μm)
0.08
0.37
1.442911
-
2.534351
e
-
8
i
0.09
1.442822
-
1.378536
e
-
8
i
0.10
1.442747
-
5.387208
e
-
9
i
0.20
1.442334
-
4.984438
e
-
12
i
0.30
1.442108
-
2.353948
e
-
14
i
0.40
1.441914
-
8.43704
e
-
17
i
0.50
1.441721
-
2.22923
e
-
19
i
0.60
1.44152
-
3.398141
e
-
19
i
0.70
1.441303
-
5.284044
e
-
19
i
0.80
1.441068
-
4.492512
e
-
20
i
0.90
1.440809
-
3.595418
e
-
19
i
Fixed diameter(d = 0.84 μm)
0.08
0.15
1.443843
-
7.982422
e
-
11
i
0.09
0.17
1.443773
-
9.149398
e
-
11
i
0.10
0.19
1.443693
-
4.61459
e
-
11
i
0.20
0.37
1.442334
-
4.984438
e
-
12
i
0.30
0.56
1.439848
-
1.198296
e
-
12
i
0.40
0.71
1.436078
-
1.976832
e
-
13
i
0.50
0.91
1.43085
-
8.180793
e
-
14
i
0.60
1.11
1.423967
-
3.413292
e
-
14
i
0.70
1.25
1.41519
-
2.279162
e
-
14
i
0.80
1.48
1.404232
-
2.485407
e
-
14
i
0.90
1.67
1.390737
-
3.793315
e
-
14
iThe normalized total energy of the fundamental mode at the central of the core for differentd/Λ ratios for (a) fixed pitch Λ=4.2 μm and (b) fixed diameter d=0.84 μm.
(a)
(b)The normalizedz-component of the Poynting flux of the fundamental mode at the central of the core for different d/Λ ratios for (a) fixed pitch Λ=4.2 μm and (b) fixed diameter d=0.84 μm.
(a)
(b)For the differentd/Λ 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 d/Λ ratios for both fixed pitch and fixed diameter.The effective refractive indexes of the fundamental mode for differentd/Λ ratios for (a) fixed pitch Λ=4.2 μm and (b) fixed diameter d=0.84 μm.
(a)
(b)


## 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 whend/Λ decreases. For a given diameter of the air holes, the confinement and the total energy confined in the core region increase when the d/Λ ratio increases up to a certain d/Λ value, for which the effective index becomes complex. However, when the d/Λ 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.



---

*Source: 567501-2013-08-06.xml*