Analytical comparison of photonic crystal fibers for dispersion compensation with different structures using FDTD method

doi:10.15406/paij.2018.02.00078

eISSN: 2576-4543

Physics & Astronomy International Journal

Review Article Volume 2 Issue 2

Analytical comparison of photonic crystal fibers for dispersion compensation with different structures using FDTD method

Faramarz E Seraji,

Verify Captcha

Regret for the inconvenience: we are taking measures to prevent fraudulent form submissions by extractors and page crawlers. Please type the correct Captcha word to see email ID.

Vajieh Arsang

Department of Communication Technology, Iran Telecom Research Center, Iran

Correspondence: Faramarz E Seraji, Faculty member Optical Communication Group, Department of Communication Technology, Iran Telecom Research Center, Tehran, Iran, Tel 98?21?84977723

Received: March 26, 2017 | Published: April 25, 2018

Citation: Seraji FE Arsang V. Analytical comparison of photonic crystal fibers for dispersion compensation with different structures using FDTD method. Phys Astron Int J. 2018;2(2):155-158. DOI: 10.15406/paij.2018.02.00078

Download PDF

Abstract

In new generation of optical communication networks, simultaneous transmission of several information channels has become possible over a single optical fiber. But dispersion and dispersion slope of optical fibers are important factors, which limit high capacity optical fiber communication systems to perform high speed transmission.

In this paper, silica core photonic crystal fibers (PCFs) with three and six air–hole rings in the cladding are investigated for effective refractive index, propagation constant and dispersion characteristics, using FDTD method. It is shown that by increasing the number of the air–hole rings in the cladding around the core, the negative dispersion of the PCFs will increase, too. The obtained negative dispersion of –300 ps/nm.km is comparatively high with respect to the published results in the related literature.

Keywords: dispersion, modeling, FDTD method, photonic crystal fibers

Introduction

Photonic crystal fibers (PCFs) have attained popularity in designing of compact optical devices utilized in optical networks^1,2 and optical sensing systems.^3,4 The widely used applications are in designing of dispersion compensators commonly used in WDM optical networks reducing the number of required repeaters.^5,6 One of the comparative characteristic features of the PCFs is their flexibilities in response to dispersion effects on propagating pulses present in a long–haul transmission optical fiber line.^7,8

Till date, the photonic crystal fibers (PCFs) are known to be new types of optical fibers with different structures and profiles as compared to conventional optical fibers that are used in optical communication systems.^9,10 Having different optical characteristics, the PCFs as optical devices have been used in various applications in optical systems, such as super continuum spectrum generation,¹¹ high optical power transfer,¹² dispersion compensations,¹³ tunable optical filters,¹⁴ Erbium–doped optical amplifiers,¹⁵ Bragg grating–based devices,¹⁶ optical nonlinear applications.^17,18

Different methods are used to evaluate dispersion properties of the PCFs with various structures.¹⁹ In this paper, by using FDTD analysis, two different PCF structures are investigated for utilization in compensations of positive dispersion in optical fiber as a transmission link.

For determination of optimal parameters values of the PCF, analytical and numerical methods, such as the scalar effective index method,²⁰ the vectorial effective index method (VEIM),²¹ the improved vectorial effective index method (IVEIM),²² and the finite difference frequency domain (FDFD) method,^23,24 are usually utilized. Couple of reports on designs of dispersion compensating fibers (DCFs) are published on optimization of compensation of the dispersions using PCFs in optical transmission bands ranging from 1460 nm to 1675 nm.^25,26

A survey shows a design of a pentagonal PCF with large flattened negative dispersion by using the full vector finite element method (FEM) where the average negative dispersions for two optimized designs were −611.9 ps/nm/km over 1,460–1,625 nm and −474 ps/nm.km over 1425–1675 nm wavelength bands, respectively.²⁷ In another report a porous–core circular PCF with circular arrangement of air holes, both in the periodic cladding and the porous core is simulated by using an efficient FEM. It has shown a flattened dispersion of $\pm$ $\pm$ 0.09 ps/THz/cm of 0.9–1.3 THz.²⁸

An ultra–flattened dispersion over the range of –210 – 15 ps/nm.km is reported in²⁹ by doping the core of a PCF with Germanium. In another report, a design of a hybrid PCF with elliptical and circular air holes of hexagonal layout with triangular lattice of five rings around the solid core is presented in,³⁰ where they obtained large flattened dispersion is of the order of 4.88 ps/nm.km over a wavelength range of 1200–1800 nm.

In our previous report, by using two–dimensional finite difference time domain (2D–FDTD), three PCF structures were proposed to optimize the dispersions by considering the effects of the geometrical parameters, such as air–hole diameters and the center–to–center spacing between the holes (pitch). The air–holes were arranged in the cladding in the form of triangular lattice by constituting 11 rings around the core.^31,32

FDTD numerical method

For a linear, isotropic, non–dispersive material with no source, the time–variant Maxwell’s equations are expressed as follows:^9,33,34

$\frac{\partial \vec{H}}{\partial t} = - \frac{1}{μ (r)} \nabla \times \vec{E}$ $\frac{\partial \vec{H}}{\partial t} = - \frac{1}{μ (r)} \nabla \times \vec{E}$ (1)

$\frac{\partial \vec{E}}{\partial t} = \frac{1}{ε (r)} \nabla \times \vec{H} - \frac{σ (r)}{ε (r)} \vec{E}$ $\frac{\partial \vec{E}}{\partial t} = \frac{1}{ε (r)} \nabla \times \vec{H} - \frac{σ (r)}{ε (r)} \vec{E}$ (2)

Where $ε (r), μ (r), σ (r)$ $ε (r), μ (r), σ (r)$ are susceptibility, permeability, and conductivity of the dielectric material, respectively, and $\vec{E}$ $\vec{E}$ and $\vec{H}$ $\vec{H}$ are respective electric and magnetic fields vectors. Using Yee cell approach, the Maxwell’s equations can be written in the form of partial differential equations with respect to time and space.³⁵

If we indicate $β$ $β$ as the propagation constant along $z$ $z$ direction in a PCF, then the expression $φ (x, y, z) = φ (x, y) e^{j β z}$ $φ (x, y, z) = φ (x, y) e^{j β z}$ will denote component of each field where $j^{2} = - 1$ $j^{2} = - 1$ . Therefore, the differentiation of Maxwell’s equations with respect to $z$ $z$ may be replaced by $j β$ $j β$ and the fields would be expressed only in terms of transverse components.

A 2D Yee unit cell across the fiber cross section is shown in Figure 1.³⁵ The $x -$ $x -$ component of first curl of the Maxwell’s equation for magnetic field is derived as:³³

$\frac{\partial H_{x}}{\partial t} = - \frac{1}{μ} (\frac{\partial E_{z}}{\partial y} - \frac{\partial E_{y}}{\partial z})$ $\frac{\partial H_{x}}{\partial t} = - \frac{1}{μ} (\frac{\partial E_{z}}{\partial y} - \frac{\partial E_{y}}{\partial z})$ (3)

Figure 1 2D yee unit cell across the fiber cross section.³⁵

By discretization of Equation (3) in terms of time and space and by using Yee cell technique, we obtain:

$\vec{H_{x}} |_{i, k}^{n + 1 / 2} = \vec{H_{x}} |_{i, k}^{n - 1 / 2} - \frac{Δ t}{μ_{i, k}} (\frac{\vec{E_{z}} |_{i, k + 1}^{n} - \vec{E_{z}} |_{i, k}^{n}}{Δ y} - j β \vec{E_{y}} |_{i, k}^{n})$ $\vec{H_{x}} |_{i, k}^{n + 1 / 2} = \vec{H_{x}} |_{i, k}^{n - 1 / 2} - \frac{Δ t}{μ_{i, k}} (\frac{\vec{E_{z}} |_{i, k + 1}^{n} - \vec{E_{z}} |_{i, k}^{n}}{Δ y} - j β \vec{E_{y}} |_{i, k}^{n})$ (4)

Where $n$ $n$ is the discrete time step, $i$ $i$ and $k$ $k$ denote grid points of discretization in $x y -$ $x y -$ plane, and $Δ t$ $Δ t$ , $Δ x$ $Δ x$ , and $Δ y$ $Δ y$ denote the time increment, distances between two adjacent points along $x$ $x$ and $y$ $y$ directions, respectively. Similarly, the components of the other fields can also be derived.

To reduce the computational time, one prefers to consider the real numbers in above equation. For this particular reason, let us consider the real parts of $E_{z}, H_{x}, and H_{y}$ $E_{z}, H_{x}, and H_{y}$ to be $\cos (β z + φ)$ $\cos (β z + φ)$ and that of $H_{z}, E_{x}, and E_{y}$ $H_{z}, E_{x}, and E_{y}$ to be $\sin (β z + φ)$ $\sin (β z + φ)$ . Therefore, Equation (4) reduces to:³⁵

$H_{x} |_{i, k}^{n + 1 / 2} = H_{x} |_{i, k}^{n - 1 / 2} - \frac{Δ t}{μ_{i, k}} (\frac{E_{z} |_{i, k + 1}^{n} - E_{z} |_{i, k}^{n}}{Δ y} - β E_{y} |_{i, k}^{n})$ $H_{x} |_{i, k}^{n + 1 / 2} = H_{x} |_{i, k}^{n - 1 / 2} - \frac{Δ t}{μ_{i, k}} (\frac{E_{z} |_{i, k + 1}^{n} - E_{z} |_{i, k}^{n}}{Δ y} - β E_{y} |_{i, k}^{n})$ (5)

The $y$ $y$ and $z$ $z$ components of magnetic fields are obtained as follows:

$H_{y} |_{i, k}^{n + 1 / 2} = H_{y} |_{i, k}^{n - 1 / 2} + \frac{Δ t}{μ_{i, k}} (\frac{E_{z} |_{i + 1, k}^{n} - E_{z} |_{i, k}^{n}}{Δ x} - β E_{x} |_{i, k}^{n})$ $H_{y} |_{i, k}^{n + 1 / 2} = H_{y} |_{i, k}^{n - 1 / 2} + \frac{Δ t}{μ_{i, k}} (\frac{E_{z} |_{i + 1, k}^{n} - E_{z} |_{i, k}^{n}}{Δ x} - β E_{x} |_{i, k}^{n})$ (6)

$H_{z} |_{i, k}^{n + 1 / 2} = H_{z} |_{i, k}^{n - 1 / 2} + \frac{Δ t}{μ_{i, k}} (\frac{E_{x} |_{i, k + 1}^{n} - E_{x} |_{i, k}^{n}}{Δ y} - \frac{E_{y} |_{i + 1, k}^{n} - E_{y} |_{i, k}^{n}}{Δ x})$ $H_{z} |_{i, k}^{n + 1 / 2} = H_{z} |_{i, k}^{n - 1 / 2} + \frac{Δ t}{μ_{i, k}} (\frac{E_{x} |_{i, k + 1}^{n} - E_{x} |_{i, k}^{n}}{Δ y} - \frac{E_{y} |_{i + 1, k}^{n} - E_{y} |_{i, k}^{n}}{Δ x})$ (7)

Similarly, one can obtain the components of electrical fields. The stability of the FDTD numerical method is determined by the following expression:³⁵

$Δ t \leq \frac{1}{c \sqrt{Δ x^{- 2} + Δ y^{- 2} + {(β / 2)}^{- 2}}}$ $Δ t \leq \frac{1}{c \sqrt{Δ x^{- 2} + Δ y^{- 2} + {(β / 2)}^{- 2}}}$ (8)

For our boundary computational treatment, we use perfectly matched layers method. One of the key characteristic parameters of PCF is chromatic dispersion, which is summation of waveguide and material dispersions. The expression for total dispersion, consisting of material and waveguide dispersion is given as:

$D_{t o t} = - \frac{λ}{2 π c} (2 \frac{d β}{d λ} + λ \frac{d^{2} β}{d λ^{2}})$ $D_{t o t} = - \frac{λ}{2 π c} (2 \frac{d β}{d λ} + λ \frac{d^{2} β}{d λ^{2}})$ (9)

To simplify the above equation, we define normalized propagation constant $β_{N} (= β / k_{0} =2 π β / λ)$ $β_{N} (= β / k_{0} =2 π β / λ)$ and substitute it in (9) to obtain:^9,18,10

$D_{t o t} = - \frac{λ}{c} (\frac{d^{2} β_{N}}{d λ^{2}}) π$ $D_{t o t} = - \frac{λ}{c} (\frac{d^{2} β_{N}}{d λ^{2}}) π$ (10)

Where $λ$ $λ$ and $c$ $c$ are wavelength and velocity of light in a vacuum, respectively. The effective refractive index of fundamental mode is defined as $n_{e f f} = β (λ / 2 π)$ $n_{e f f} = β (λ / 2 π)$ .³⁶

To evaluate Equation (10) in terms of $β_{N}$ $β_{N}$ , $d β_{N} / d λ$ $d β_{N} / d λ$ , $d^{2} β_{N} / d λ^{2}$ $d^{2} β_{N} / d λ^{2}$ , we use FDTD method. The refractive index of the core material is determined by using Sellmeier formula as:^18,37

$n_{s i l i c a}^{2} (λ) = 1 + \sum_{k = 1}^{3} \frac{a_{k} λ^{2}}{λ^{2} - b_{k}^{2}}$ $n_{s i l i c a}^{2} (λ) = 1 + \sum_{k = 1}^{3} \frac{a_{k} λ^{2}}{λ^{2} - b_{k}^{2}}$ (11)

Where $λ$ $λ$ is the wavelength of light, and $a_{k}$ $a_{k}$ and $b_{k}$ $b_{k}$ are the Sellmeier coefficients. The practical values of $a_{k}$ $a_{k}$ and $b_{k}$ $b_{k}$ for $k = 1, 2, 3$ $k = 1, 2, 3$ are given as:³⁸

$a_{1} = 0.6961663, a_{2} = 0.4079426, a_{3} = 0.8974794,$ $a_{1} = 0.6961663, a_{2} = 0.4079426, a_{3} = 0.8974794,$

$b_{1} = 0.0684043, b_{2} = 0.1162414, b_{3} = 9.896161.$ $b_{1} = 0.0684043, b_{2} = 0.1162414, b_{3} = 9.896161.$

In this paper, by using the above equations and other related expressions in FDTD method, two PCF samples with different numbers of air–hole rings are selected for design in order to study their dispersion behaviors.

PCF with three rings of air–holes

The selected PCF structure with air–hole of equal diameters $d = 0.75 μ m$ $d = 0.75 μ m$ and air–hole spacing of $Λ = 2 μ m$ $Λ = 2 μ m$ arranged in three rings in the form of triangular lattice is illustrated in Figure 2. The core diameter is taken as $1.45 μ m$ $1.45 μ m$ .

The spectrum of chromatic dispersion of fundamental mode in the wavelength range of 1100–2000 nm is illustrated in Figure 3 (Red curve). The dispersion at 1510 nm is obtained as –140 ps/nm.km. The dispersion of the proposed PCF in this case is nearly twice the dispersion of conventional dispersion compensating fiber (i.e.,–80 ps/nm.km).

Figure 2 The designed PCF with three rings of air–holes.

Figure 3 Dispersion spectra of the PCFs with three and six circular air–hole rings in the cladding region.

The PCF with six rings of air–holes

As of previous analysis, a PCF structure is considered with the same air–hole diameter of $d = 0.75 μ m$ $d = 0.75 μ m$ , core diameter of $1.45 μ m$ $1.45 μ m$ and a different air–hole spacing of $Λ = 1.1 μ m$ $Λ = 1.1 μ m$ . In this case, six air–hole rings are considered in the cladding region (in this case, figure not shown).³⁶

By using the home–made programming code, in this case as well, the chromatic dispersion spectrum has been depicted in Figure 3 (Blue curve). In this case, the dispersion value of –300 ps/nm.km is obtained, which twice the value is obtained in the previous design of the PCF with three circular air–hole rings. The results obtained here is comparable to the ones reported,¹⁸ in which the proposed structure with seven air–hole rings were used. The reported dispersion value was –150 ps/nm.km.

In the simulations performed by FDTD method, by varying $β$ $β$ values in each case, the resonance fundamental frequencies, wavelengths, and effective refractive indices are determined and tabulated in Table 1.^17,18

$\times 1 0^{6} β$ $\times 1 0^{6} β$	$n_{e f f}$ $n_{e f f}$	$λ (μ m)$ $λ (μ m)$	Fundamental Frequency (THz)
5	1.3010	1.6349	183.50
6	1.3343	1.3973	214.70
7	1.3575	1.2185	246.20
8	1.3745	1.0795	277.90
9	1.3884	0.9631	309.50
12	1.4116	0.7391	405.90
15	1.4239	0.59641	503.00

Table 1 Calculated effective refractive index used for determination of dispersion of PCF with six air-hole rings

Conclusion

In this paper, by using FDTD method, two proposed photonic crystal fibers have been analyzed with three and six air–hole rings in the cladding arranged in triangular lattices. In the case of three air–hole rings in the cladding of the PCF, for air–hole with equal diameters of $d = 0.75 μ m$ $d = 0.75 μ m$ and air–hole spacing of $Λ = 2 μ m$ $Λ = 2 μ m$ , the dispersion value is determined to be –140ps/nm.km at wavelength 1510nm.

In case of the PCF with six air–hole rings, air–hole spacing of $Λ = 1.1 μ m$ $Λ = 1.1 μ m$ , and with the same air–hole of diameter of $d = 0.75 μ m$ $d = 0.75 μ m$ , the dispersion value is found to be –300ps/nm.km. In both the cases, the diameters of the PCF core are taken as $1.45 μ m$ $1.45 μ m$ .

It is further shown that when the air–hole spacing $Λ$ $Λ$ is decreased by a ratio of 55% and by doubling the number of air–hole rings, the dispersion will increase more than twice the value. The PCF with high negative dispersion value is suitable for design of high dispersion compensating fiber used in long–haul optical transmission links.

Acknowledgements

The authors greatly acknowledge the allotment of the post–graduate student by GJK Institute of Higher Education, Qazvin, Iran, for collaboration in the running academic project of optical group at ITRC.