Optimization Search Space of Optical and Geometrical Parameters.
1. Introduction
The demand for data communication is growing rapidly due to the increasing popularity of the Internet and other factors [1]. The unprecedented success of information technology in recent years ushers an explosive growth in demand for the Internet access in the 21st century. Ultra-wideband transmission media are needed in order to provide high-speed communications for a much larger number of users. A promising solution to the capacity crunch can come from wavelength-division-multiplexed optical fiber communication systems that are shown to provide enormous capacities on the order of terabit per second over long distances. These systems utilize single-mode fibers, in conjunction with erbium-doped fiber amplifiers, as the transmission medium [2]. The optical silica fiber might be the only proper choice to realize this task. Optical fiber based communication is the excellent alternative for these purposes which needs low dispersion as well as dispersion slope and large bandwidth supported by optical physical medium[3]. Optical transmission began to be used in trunk cables about 1990; the capacity of those systems was several hundred Mbit/s per fiber. The capacity jumped to 2.5 (5.0) Gbit/s per fiber with the introduction of optical repeaters using erbium-doped fiber amplifiers in 1995. It jumped again in 1998, to 10 (20) Gbit/s per fiber, with the introduction of wavelength division multiplexing (WDM) [5]. Overall transmission capacity now exceeds 100 Gbit/s per fiber due to improvements in WDM techniques [6]. Usually, those techniques are called dense WDM (DWDM). In recent years, the increasing demands for transmission capacity have led to intense research activities on high capacity DWDM communication system [7].
Nowadays, applications such as optical time division multiplexing (OTDM) and dense wavelength division multiplexing (DWDM) are usual tasks in industry [1]. Therefore by considering these applications, providing a large bandwidth and high-speed communication possibility using optical fibers is highly interesting [3].
In the following, we review requirements for DWDM, Dispersion properties, optical nonlinearity, loss properties, and design of optical fiber for DWDM.
2. Requirements for DWDM
The number of wavelengths (channels) in the fibers of DWDM systems is increasing. As discussed earlier, a WDM signal typically occupies a bandwidth of 30 nm or more, although it is bunched in spectral packets of bandwidth
3. Dispersion properties
Silica fibers suffer from some disadvantages especially dispersion and dispersion slope. Meanwhile, these two factors cause severe restrictions for high-speed pulse propagation [1]. The dispersion value becomes larger by the wavelength increasing in the conventional optical fibers. So owing to the dissimilar broadening for different channels, the multi-channel application realization would be hard. A suitable optical fiber should meet the small dispersion as well as the small dispersion slope in the predefined wavelength interval. The dispersion properties are the dispersion itself and the dispersion slope of the optical fiber. The dispersion value cannot be in the zero-value region because FWM causes interaction between signals (optical channels) in DWDM systems when there is phase matching between the optical channels due to zero dispersion. Therefore, the dispersion value in the signal wavelength region must have the proper non-zero value. The sign of the dispersion value should be positive for short-distance transmission and negative for ultra-long-distance transmission [12] because of modulation instability in the positive dispersion in a long link. When the signal wavelength band becomes wider, the difference in the dispersion values at the edges of the wavelength band becomes larger. Dispersion compensation thus becomes difficult for long distance DWDM transmission. To achieve both long distance and high speed transmission with easy dispersion compensation for a wide wavelength band, the dispersion slope should be reduced. For single wavelength communication, dispersion shifted fiber is enough. But for applications such as DWDM this method cannot provide high speed possibility. In these applications, the physical media should provide the flat, minimum, and uniform dispersion as well as dispersion slope ideally. An important limitation induced by chromatic dispersion and its slope is broadening factor which restricts the bit rate parameter.
To minimize pulse broadening in an optical fiber, the chromatic dispersion should be low over the wavelength range used. A fiber in which the chromatic dispersion is low over a broad wavelength range is called a dispersion-flattened fiber.
4. Optical nonlinearity
The response of any dielectric to light becomes nonlinear for intense electromagnetic fields, and optical fibers are no exception. Even though silica is intrinsically not a highly nonlinear material, the waveguide geometry that confines light to a small cross section over long fiber lengths makes nonlinear effects quite important in the design of modern light wave systems [1]. The much higher power level due to simultaneous transmission of multiplexed channels and propagation over much longer distances made possible with the utilization of fiber amplifiers, cause the otherwise weak and negligible fiber nonlinearities to affect the signal transmission significantly [2]. When the number of signal wavelengths carried in an optical fiber increases, the average transmission power density becomes larger than that in conventional systems. Consequently, optical-fiber nonlinearities have emerged as a main issue. This nonlinearity seriously limits transmission capacity with various nonlinear interactions, which are generally categorized as scattering effects and optical signal interactive modulation. Because the signal power density is stronger due to the greater number of channels in DWDM systems, optical fiber nonlinearity limits the number/spacing of the channels and the length/speed of the transmission. In general, the refractive index of optical fiber has a weak dependence on optical intensity (equal to signal power (p) per effective area (Aeff) in the fiber. Optical fiber nonlinearity arises from modulation of the refractive index caused by changes in the optical intensity of the signal. This cause four wave mixing (FWM), self-phase modulation (SPM) and cross phase modulation (XPM) can be observed in the fiber. The XPM and SPM distort the signals. Therefore, optical fiber nonlinearity must be reduced. The most practical way to do this is to enlarge Aeff [13]. The relationship between Aeff and mode field diameter (MFD) is direct and proportional. As a result, enlarging MFD is a practical solution for low nonlinearity. The choice of dispersion shifted fibers (DSFs) along with erbium-doped fiber amplifiers (EDFAs), for operation at 1550nm window, would be an ideal one to achieve greater transmission distance and utilize full capacity of transmission system [5,7,15]. However, when the system is operated at the zero dispersion wavelengths, the nonlinear interaction between the channels and noise components is increased. The system working slightly away from the zero dispersion wavelengths can reduce these unwanted interactions. The WDM system reduces the nonlinear effects and enables multi-wavelength transmission through non-zero dispersion shifted fibers having very small dispersion in duration 1530-1610 nm. In order to increase the information carrying capacity, latest high speed communication system is based on the dense wavelength division multiplexing/demultiplexing (DWDM) [16, 17]. In such systems, nonlinear effects like four wave mixing (FWM), which arise due to simultaneous transmission at many closely spaced wavelengths and high optical gain from EDFA, imposes serious limitations on the use of a DSF with zero dispersion wavelength at 1550 nm [18,19]. To overcome this difficulty, the nonzero dispersion shifted fibers having small dispersion in the range ~ 2−4 ps
5. Loss properties
Progress in optical fiber fabrication technologies has resulted in a routine production of low loss single mode fibers. This enables us to apply the single mode fibers promisingly in high bit rate and long haul optical transmission systems. Structural optimization must be established so as to provide desirable transmission characteristics for given operating conditions. A basic design consideration has been made by taking into account transmission characteristics such as fiber intrinsic loss, bending loss, splice loss, and launching efficiency [22, 23]. Use of commercially available erbium doped fiber amplifiers (EDFA), which forces optical communication systems to be operated in the 1550 nm window, has significantly reduced the link length limitation imposed by attenuation in the optical fiber [15]. The fiber loss is one of the significant restrictions in the optical fiber communication links. It is one of some reasons limit the maximum distance that information can be sent without presence of the repeaters. Meanwhile, due to the loss, the pulse amplitude reduces so that the initial information cannot be restored in the noisy conditions. Seeing that, in the fiber design one likes to shift the zero dispersion wavelength to the region that the fiber has the lowest level attenuation. The combination of natural attenuations has a global minimum around 1.55 μm and that is why most optical communication systems are operated at this wavelength [4, 18]. A kind of loss which must be taken into account in fiber design is the bending loss. Every time an optical fiber is bent, radiation occurs. When a bent occurs, a portion of the power propagating in the cladding is lost through radiation.
6. Design of DWDM fiber
There are three methods to increase the capacity of a DWDM transmission system, using a broad wavelength range, narrowing channel spacing and increasing a bit rate per channel. However, one of disadvantages for the last two methods is the degradation of the transmission performance due to optical nonlinear effects. In this area, there are three categories which cover all designs. There are based on using zero dispersion shifted fibers (ZDSFs), non-zero dispersion shifted fiber (NZDSFs) and dispersion flattened fibers (DFFs).
7. Zero Dispersion Shifted Fibers (ZDSFs)
Use of commercially available erbium doped fiber amplifiers (EDFA), which forces optical communication systems to be operated in the 1550 nm window, has significantly reduced the link length limitation imposed by attenuation in the optical fiber. However, high bit rate (~10 Gb
For calculation of dispersion and dispersion slope the following parameters are used.
where P and Q are geometrical parameters. Also, the optical parameters for the structure are defined as follows.
For evaluating of the index of refraction difference between core and cladding the following definition is done.
Here, we propose a novel methodology to make design procedure systematic. It is done by the aim of optimization technique and based on the Genetic Algorithm. A GA belongs to a class of evolutionary computation techniques [28] based on models of biological evolution. This method has been proved useful in the domains that are not understood well; search spaces that are too large to be searched efficiently through standard methods. Here, we concentrate on dispersion and dispersion slope simultaneously to achieve to the small dispersion and its slope in the predefined wavelength duration. Our goal is to propose a special fitness function that optimizes the pulse broadening factor. To achieve this, we have defined a weighted fitness function. In fact, the weighting function is necessary to describe the relative importance of each subset in the fitness function [24]; in other words, we let the pulse broadening factor have different coefficient in each wavelength. To weight the mentioned factor in the predefined wavelength interval, we have used the Gaussian weighting function. The central wavelength (λ0) and the Gaussian parameter (σ) are used for the manipulation of the proposed fitness function and their effects on system dispersion and dispersion slope. To express the fiber optic structure, we considered three optical and geometrical parameters. According to the GA technique, the problem will have six genes, which explain those parameters. It should be mentioned that the initial range of parameters are chosen after some conceptual examinations. The initial population has 50 chromosomes, which cover the search space approximately. By using the initial population, the dispersion (β2) and dispersion slope (β3), which are the important parameters in the proposed fitness function, can be calculated. Consequently elites are selected to survive in the next generation. Gradually the fitness function leads to the minimum point of the search zone with an appropriate dispersion and slope. Equation (6) shows our proposal for the weighted fitness function of the pulse broadening factor.
where
The flowchart given in Figure 2 explains the foregoing design strategy clearly.
To illustrate capability of the suggested technique and weighted fitness function, the MII triple-clad optical fiber is studied, and the simulated results are demonstrated below. In the presented figures, we consider four simulation categories including dispersion related quantities, nonlinear behavior of the proposed fibers, electrical field distribution in the structures, and fiber losses.
For all the simulations, we consider λ0=1500, 1550 nm and σ = 0, 0.027869 and 0.036935 µm as design constants. To apply the GA for optimization, we consider the search space illustrated in Table 1 for each parameter as a gene. The choice of these intervals is done according to two items. The designed structure must be practical in terms of manufacturing and have high probability of supporting only one propagating mode [24].
Parameter | a (µm) | p | Q | R 1 | R 2 | Δ |
duration | [2-2.6] | [0.4-0.9] | [0.1-0.7] | [0.05-0.99] | [(-0.99)- (-0.05)] | [2×~10 -3 - 1×~10 -2] |
The wavelength and distance durations for optimization are selected as follows. For λ0=1550nm: 1500 nm<λ< 1600 nm, for λ0=1500 nm: 1450 nm <λ< 1550 nm, and 0 <
λ -0(µm) | a (µm) | Δ | R 1 | R 2 | p | Q | |
σ=0 | 1.55 | 2.0883 | 8.042e-3 | 0.5761 | -0.4212 | 0.7116 | 0.3070 |
1.5 | 2.1109 | 7.036e-3 | 0.6758 | -0.2785 | 0.8356 | 0.2389 | |
|
1.55 | 2.0592 | 9.899e-3 | 0.7320 | -0.2670 | 0.7552 | 0.2599 |
1.5 | 2.5822 | 9.111e-3 | 0.5457 | -0.4237 | 0.7425 | 0.2880 | |
|
1.55 | 2.2753 | 9.933e-3 | 0.5779 | -0.4218 | 0.6666 | 0.3428 |
1.5 | 2.5203 | 9.965e-3 | 0.4867 | -0.3841 | 0.6819 | 0.3324 |
It is found that optimization method for precise tuning of the zero dispersion wavelengths as well as the small dispersion slope requires large value for the index of refraction difference (Δ). That is to say that providing large index of refraction is excellent for the simultaneous optimization of zero dispersion wavelength and dispersion slope. First, we consider the dispersion behavior of the structures. To demonstrate the capability of the proposed algorithm for the assumed data, the obtained dispersion characteristics of the structures are illustrated in Figure 3. It shows that the zero dispersion wavelengths can be controlled precisely by controlling the central wavelength. Meanwhile, the Gaussian parameters are used to manipulate the dispersion slope of the profile. Considering Figure 3 and Table 3, it is found that the zero value for the Gaussian parameter can tune the zero dispersion wavelengths accurately (~100 times better than other cases).
Second, the dispersion slope is examined. The presented curves say that by increasing the Gaussian parameter the dispersion slope becomes smaller, and it is going to be smooth in large wavelengths. Furthermore it is clear that there is a trade-off between tuning the zero dispersion wavelengths and decreasing the dispersion slope as shown in Figure 3, 4, and Table 3.
type |
|
Dispersion |
Dispersion Slope |
Effective Area |
Mode Field Diameter |
Quality Factor |
|
1.55 | -2.57e-4 | 0.0695 | 191.92 | 7.95 | 3.04 |
1.5 | 2.55e-5 | 0.0828 | 344.15 | 9.76 | 3.61 | |
|
1.55 | -0.013 | 0.0647 | 194.79 | 7.12 | 3.85 |
1.5 | 0.008 | 0.0597 | 209.95 | 6.70 | 4.68 | |
|
1.55 | -0.085 | 0.0592 | 150.05 | 6.82 | 3.22 |
1.5 | -0.089 | 0.0564 | 164.21 | 6.55 | 3.82 |
The normalized field distribution of the MII based designed structures is illustrated in Figure 5 and 6. Because of the special structure, the field distribution peak has fallen in region III. As such most of the field distribution displaces to the cladding region. In addition it is observed that the field distribution peak is shifted toward the core, and its tail is depressed in the cladding region by increasing the Gaussian parameter (except σ=0). On the other hand the field distribution slope increases inside the cladding region by increasing of the Gaussian parameter.
The effective area or nonlinear behavior of the suggested structures is illustrated in Figure 7. It is observed that the effective area becomes smaller by increasing the Gaussian parameter. Figure 5–7, and Table 3 indicate a trade-off between the large effective area and the small dispersion slope. The results illustrated in Figure 4 show that the dispersion slope reduces by increasing the Gaussian parameter. However the field distribution shifts toward the core, which concludes the small effective area in this case. Foregoing points show that there is an inherent trade-off between these two important quantities.
The mode field diameter that corresponds to the bend loss is illustrated in Figure 8 and 9 for both central wavelengths. It is clearly observed that the mode field diameter decreases by increasing the Gaussian parameter. In other words, the Gaussian parameter is suitable for the bend loss manipulation in these structures. Furthermore, Table 3 shows that the mode field diameter is ~7µm in the designed structure.
As another concept to consider, Table 3 says that the mode field diameter is not affected noticeably by increasing the effective area. This is the origin of raising the quality factor in these structures. This is a key point why the average amount of the quality factor in the proposed structures is increased in Figure 9. The quality factor of the designed fibers is illustrated in Figure 10. The calculations show that the quality factor is generally larger than 3. It is mentionable that the quality factor is smaller than unity in the inner depressed clad fibers (
As another result the dispersion length is illustrated in Figure 11 for the given Gaussian parameter and two central wavelengths. The narrow peaks at λ=1500nm and 1550nm imply the precise tuning of the zero dispersion wavelengths. The higher-order dispersion length of the designed fibers is demonstrated in Figure 12. It is clear that the higher-order dispersion length increases by raising the Gaussian parameter.
In the following, the nonlinear effect length for 1 mW input power is illustrated in Figure 13. First, it can be extracted that the suggested structures have the high nonlinear effect length. For the general distances, these simulations show that the fiber input power can become some hundred times greater to have the nonlinear effect length comparable with the fiber dispersion length. Second, the nonlinear effect length decreases and increases, respectively, by raising the Gaussian parameter and wavelength.
The amount of the fiber bending loss strongly depends on the bend radius and the mode field diameter. Figure 14 and 15, respectively, illustrate the bending loss (dB/m) versus the bending radius (mm) at λ0= 1550 nm and 1500nm with variance of the weighting function (σ) as a parameter. According to Figure 8, 9, 14, and 15, it is clear that smaller mode field diameter yields to the greater tolerance to the bending loss.
All of the presented outcomes show that the suggested idea has capability to introduce a fiber including higher performance. We have presented a novel method that includes the small dispersion, its slope, high effective area, and small mode field diameter simultaneously [24]. So all options required for the zero dispersion shifted communication system are achieved successfully. This advantage is obtained owing to the selection of the basic fiber structure as well as the method of optimization. Our selected fiber structure is the MII, and we use the weighted fitness function applied in the GA for optimization. By combining the suitable structure and the novel optimization method, all of the stated advantages can be gathered simultaneously. The features of the proposed method are capable of being extended to all of fiber structures, introduce two parameters instead of multi-designing parameters, and tune the zero dispersion wavelengths precisely.
The ring index profiles fibers have been closely paid attentions because it has the larger effective-areas that can minimize the harmful effects of fiber nonlinearity [29]. For the proposed MII fiber structures, the small dispersion and its slope have been obtained thanks to a design method based on genetic algorithm. But there is not any concentration on the bending loss characteristic at the design process. Here we want to enter bending loss effect on the fitness function directly and attempt to present an optimized RII triple-clad optical fiber to obtain the wondering performance from dispersion, its slope, and bending loss points of view. The index refraction profile of the RII fiber structure is shown in Figure 16.
To calculate the dispersion, its slope and bending loss characteristics of the structure, the geometrical and optical parameters are defined as follows.
The design method is based on the combination of the Genetic Algorithm (GA) and Coordinate Descent (CD) approaches. It is well known that the GA is the scatter-shot and the CD is the single-shot searching technique. The single-shot search is very quick compared to the scatter-shot type, but depends critically on the guessed initial parameter values. This description indicates that for the CD search, there is a considerable emphasis on the initial search position. In this method, it is possible to define a fitness function and evaluate every individuals of the population with it. So we have combined the CD and GA methods to improve the initial point selection with the help of generation elite and inherit the quick convergence of coordinate descent [30]. In other words, we cover and evaluate the answer zone by initial population and deriving few generations and use the elite of the latest generation as an initial search position in the CD (Figure 17).
To derive the suggested design methodology, the following weighted cost function is introduced. We have normalized the pulse broadening factor in the manner to be comparable with bending loss. This normalization is essential to optimize the pulse broadening factor and bending loss simultaneously. If not, the bending loss impact will be imperceptible and be lost in the broadening factor term.
The bending radius is set on
To show the capability of the proposed algorithm, Table 4 is presented to clarify the different characteristics of these three structures. By considering on Figure 18 and Table 4, it is clear that there is a trade-off between the zero dispersion wavelength tuning and the dispersion slope decreasing. In other words, it is found out that the zero value for the
The effective area or nonlinear behavior of the suggested structures is listed in Table 4. These values are high enough for the optical transmission applications. Owing to the special structure of the RII type fiber, the field distribution peak has fallen in the first cladding layer. As such most of the field distribution displaces to the cladding region. This is the origin of large effective area in the designed structures. The normalized field distribution of the RII based designed structures is illustrated in Figure 19.
type | D ( λ=1.55 μm ) (ps/km/nm) | S ( λ=1.55 μm ) (ps/km/nm 2 ) | BL ( λ=1.55 μm ) (dB/m) | A eff ( λ=1.55 μm ) (μm 2 ) |
σ = 0.0 | 1.38e-4 | 0.048 | 1.90e-2 | 86.84 |
σ =1.12e-8 | -6.15e-4 | 0.041 | 1.67e-1 | 82.53 |
σ =3.69e-8 | 4.50e-2 | 0.035 | 4.66e-2 | 86.01 |
Due to the refractive index thermo-optic coefficient and the thermal expansion coefficient, the optical and geometrical parameters are altered. Consequently, the optical transmission characteristics of the optical fiber such as dispersion, its slope and bending loss are confronted to change. In order to evaluate the thermal stability of the designed structures, the following results are extracted and presented in Table 5. The
In the least design we have focused on RII depressed core triple clad single mode optical fiber and presented a combined optimization approach to obtain desirable design goals. Furthermore, we have used the special fitness function including dispersion, its slope and bending loss impacts simultaneously. With application of this fitness function in the case of higher
Type | dD/dT (ps/km/nm/°C) | dS/dT (ps/km/nm 2 /°C) | dλ 0 /dT (nm/°C) | dBL/dT (dBL/m/°C) |
σ = 0.0 | -1.22 ×~ 10 -3 | +2.83 ×~ 10 -6 | +2.5 ×~ 10 -2 | +3.97 ×~ 10 -6 |
σ =1.12e-8 | -1.21×~ 10 -3 | +2.93 ×~ 10 -6 | +3.33 ×~ 10 -2 | +2.70 ×~ 10 -5 |
σ =3.69e-8 | -1.21 ×~ 10 -3 | +2.93 ×~ 10 -6 | +2.5 ×~ 10 -2 | +8.79 ×~ 10 -6 |
[ ( -1.77 ) - ( +1.77 ) ] ps/km/nm and [ ( 0.037 ) – ( 0.033 ) ] ps/km/nm2. Also the amount of bending loss at 1.55 μm with 1cm radius of curvature and effective area are 4.66e-2 dB/m and 86.01 μm2 respectively. In the meantime, the thermal stabilities of the designed structures are evaluated. It is possible to design zero dispersion shifted by using graded index structure. The main options are dispersion value and the effective area at 1550nm to minimize pulse broadening and nonlinearity effects. Excess investigation of large mode area fibers show that there is not serious focusing on design of zero dispersion shifted fiber based on the graded index refractive structures. The index refraction profile of the triangular-core graded index optical fiber structure, which is suggested by us for the first time, is shown in Figure 20. It is clear that the proposed graded index fiber has a linear variation in core region. According to the TMM approach, it is assumed that the refractive index of the fiber with an arbitrary but axially symmetric profile is approximately expressed by a staircase function. So the field distribution and guided modes are calculated [26].
Also for easy handling of the problem and calculating the dispersion and its slope of the proposed fiber, following optical and geometrical parameters are defined.
In order to explain layers’ refractive index,
In the defined fitness function, the summation is proposed to include optimum broadening factor for each length up to 200 km. The short glance on eq.(10) shows that the above fitness function is a limited version of the past one [24], which the wavelength duration optimizing is abbreviated. Also, it is not weighted because the fitness function is evaluated at single wavelength (
Core radius | Δ | p | η | μ |
*1.8μm | 8×~10 -3 | 0.5 | 0.25 | 0.44 |
The value of the dispersion is –0.04549ps/km/nm at 1.55μm, which is properly small. The nonlinear effects in a single mode fiber are the ultimate restricting factors for bit rate and distance in long haul optical fiber communication system. Therefore, the large effective area single mode fibers have been the subject of considerable studies recently. The effective area of the suggested structure is 65.3 μm2 at 1.55μm, which is acceptable for this application. The associated mode field diameter at aforesaid wavelength is 9.03μm.
8. Non-Zero Dispersion Shifted Fibers (NZDSFs)
Use of commercially available erbium doped fiber amplifiers (EDFA), which forces optical communication systems to be operated in the 1550 nm window, has significantly reduced the link length limitation imposed by attenuation in the optical fiber. However, high bit rate (~ 10 Gb
where
From Eq. (11) it is found that the exact refractive index profile are controlled by seven parameters (i.e.,
The Gaussian approximation method is used for calculating electrical field distribution for the designed refractive index profile. The designed fiber has a dispersion about 4 ps/km/nm and dispersion slope about 0.06 ps/km/nm2 at 1.55 µm operating wavelength, which can be used to avoid four-wave mixing (FWM). Also the zero dispersion wavelength is adjusted near 1480nm. In addition, calculation also shown that the designed fiber not only has a large effective area over 100 µm2 but also has low bending loss (<1.3×10-3 dB with 30 mm bending radius and 100 turns) and low splice loss (<6.38×10-3 dB) with conventional fiber. In order to broaden the wavelength range – not only the conventional C-band (1530–1565 nm) – but also the L-band (1565–1620 nm) and S-band (1460–1530 nm), the dispersion slope should be as low as possible, which can decrease the cost of dispersion compensation and suppress the self-phase modulation, especially for the 40 Gb/s communication system. Now we present a new depressed core index dual ring profile, which can provide a large Aeff and a low dispersion slope simultaneously, and the zero dispersion wavelength is less than 1430 nm [32]. The fibers with this refractive index profile have the low bending loss and low intrinsic loss. The splice loss also can reach the accepted value as it is spliced with conventional single mode fiber. The principal design requirements for the fibers are large Aeff, small dispersion slope, low bending and intrinsic loss, low polarization-mode dispersion, and zero dispersion wavelength that should be lower than 1430 nm.
Relative refractive index Δni is defined by the equation: Δni=(ni-nc)/nc, where nc is the outer cladding layer’s refractive index. The dopant in the glass can decrease the glass viscosity, i.e, more dopant concentration means less glass viscosity. High difference value between Δn1 and Δn2 may cause high difference value in viscosity property of the depressed core layer and the first raised ring. Therefore, more mechanical stress will be built in the optical fibers during the drawing process [33]. On the other hand, high difference value between layers’ refractive index can increase the compositional variation in the optical fiber. Therefore, more thermal stress can be also caused by the radial variation of thermal expansion coefficient due to the big compositional variation in optical fiber [34]. The residual stress in optical fibers can not only weaken the strength of optical fibers, but also increase the fiber’s attenuation and polarization mode dispersion values. According to the above description, the refractive index profile shown in Figure 22 is proposed to overcome problems. Every parameters of the fiber profile are set out in Table 6, where Δni are the relative refractive index of different layer from the depressed core layer to the cladding, respectively.
Δn 1 (%) | Δn 2 (%) | Δn 3 (%) | Δn 4 (%) | Δn 5 (%) | r 1 (µm) | r 2 (µm) | r 3 (µm) | r 4 (µm) | r 5 (µm) |
0.14 | 0.57 | -0.27 | 0.30 | -0.18 | 2.50 | 4.10 | 6.88 | 9.98 | 12.41 |
Table 7 shows the optical characteristic of fabricated fiber designed according to the refractive index profile parameters as Table 6, where MFD, RDS are the mode field diameter and relative dispersion slope, respectively. It is noted that the fiber has a large Aeff of 105 µm2 and a small dispersion slope of about 0.065 ps/ km /nm2 simultaneously. Macro bending loss at 1550 nm is less than 0.05 dB/km (100 turns on the 60 mm diameter mandrel).
Parameters item | Wavelength (nm) | Typical value |
Dispersion (ps/km/nm) | 1460 | 3.465 |
1550 | 9.569 | |
1625 | 14.324 | |
Dispersion Slope(ps/km/nm 2 ) | 1550 | 0.0648 |
Zero dispersion wavelength (nm) | - | 1413 |
Attenuation (dB/km) | 1550 | 0.210 |
MFD (μm) | 1550 | 10.2 |
A eff (μm 2 ) | 1550 | 105.6 |
A eff ×~dispersion | 1550 | 1010.5 |
RDS (nm -1 ) | 1550 | 0.0068 |
PMD(ps/km 0.5 ) | 1550 | 0.04 |
Macro bending loss for 100 turns on the 60mm diameter (dB/km) | 1550 | 0.006 |
1625 | 0.015 |
From the Table 7 we can see that the zero dispersion wavelength is below 1430 nm, the dispersion at 1460, 1550 and 1625 nm are 3.465, 9.569 and 12.324 ps/km/nm, respectively. Therefore, this fiber not only can be used at the conventional C band for transmission link, but also can be suited for S-band and L-band. With the progress being made in the practical application of the Raman amplifier, this fiber also can be applicable in transmission links using distributed Raman amplifier in the future. Furthermore, the value of Aeff × dispersion is also large enough to suppress the dispersion-related non-linear effects in the transmission system [35]. It is obvious that a large effective area fiber with non-zero dispersion about 4 ps/km/nm at 1.55 µm wavelength band is a good approach to avoid four-wave mixing (FWM) effect, which in turn enhances the performance of wavelength division multiplexing (WDM) system [31]. However, such large effective area fibers have relatively large dispersion slope which also restricts the numbers of different wavelength signals [21, 36-38]. Therefore, considerable efforts are being made to reduce the dispersion slope of such fibers to deal with a rapid progress of DWDM system [39,40]. The refractive index profile of the fiber is shown in Figure 23. This profile has been named RI type in the literatures.
In order to obtain the flat modal field over the entire central dip region, the effective index (
The values of the MFD and
a(μm) | b(μm) | d(μm) | ∆ 1 (%) | ∆ 2 (%) | ∆ 3 (%) |
1.0 | 3.1 | 6.5 | 0.03 | 0.48 | -0.20 |
The dispersion and dispersion values of the designed fiber at
9. Dispersion Flattened Fibers (DFFs)
The dispersion value becomes larger by the wavelength increasing in the conventional optical fibers. So owing to the dissimilar broadening for different channels, the multi-channel application realization would be hard. A suitable optical fiber should meet the small dispersion as well as the small dispersion slope in the predefined wavelength interval [42]. The concept of providing the attractive option of low dispersion over a range of wavelengths was first suggested by Kawakami and Nishida in 1974 [43, 44]. They proposed the original
where C is the chromatic dispersion. A normalized W profile is given by
where
This procedure is repeated for different combinations of N1 and N2, and the result is given in Table 9. The first column of this table, i.e.,
According to Table 9, the global minimum is 0.9 ps/km /nm, and the corresponding optimal W fiber is
N 2 | |||||
N 1 | 1.000 | 0.995 | 0.990 | 0.985 | 0.980 |
1.002 | 12 | 12 | 12 | 12 | 12 |
1.005 | 8.6 | 8.6 | 8.5 | 8.5 | 8.5 |
1.010 | 4.8 | 4.3 | 4.2 | 4.1 | 4.1 |
1.015 | 7.9 | 2.3 | 1.8 | 1.6 | 1.5 |
1.020 | 14 | 1.3 | 0.9 | 1.0 | 1.2 |
1.025 | 21 | 1.1 | 1.7 | 2.4 | 2.9 |
1.030 | 27 | 4.0 | 2.4 | 3.7 | 4.5 |
An exhaustive method for calculating the minimum rms chromatic dispersion in W fibers has been presented [45]. The procedure is to generate all W fibers with a certain cutoff wavelength and then find the minimum rms dispersion by one-dimensional minimization followed by direct inspection. It was found, in the case investigated, that the W fiber is capable of dispersion flattening only if high doping levels are used. Dispersion and its slope are responsible on the pulse broadening [4]. So we believe that mixing and gathering these parameters on the design procedure would lead us to a fiber with exciting performances. In other words, we concentrate on dispersion and slope simultaneously to achieve the small dispersion and its slope in the predefined wavelength interval, the band which we want to have flat dispersion behavior [42]. We use WII-type optical fiber structure which its refractive index profile is shown in Figure 27.
Also, for easy handling of the problem the following optical parameters are defined as follows.
For this structure the geometrical parameters are introduced in the following.
To achieve the flattening purpose, we have defined a weighted fitness function. Equation (17) shows our proposal for the weighted fitness function of the un-chirped pulse broadening factor.
It is useful to say that at first we applied this fitness function to design dispersion shifted fiber. But outcomes presentation will show that this function is appropriate in Flattening application too [42]. By using this fitness function, the zero dispersion wavelength can be shifted to the central wavelength (λ0). Since, the weight of the pulse broadening factor at λ0 is greater than others in the weighted fitness function; it is more likely to find the zero dispersion wavelength at λ0 compared to the other wavelengths. In the meantime, the flattening of the dispersion curve is controlled by Gaussian parameter (σ). To put it in another way, the weighting Gaussian function becomes broader in the predefined wavelength interval by increasing the Gaussian parameter (σ). As a result, the effect of the pulse broadening factor with greater value is regarded in different wavelengths, which causes a considerable decrease in the dispersion slope in the interval. Consequently, the zero dispersion wavelength and dispersion slope can be tuned by λ0 and σ, respectively. The wavelength and the distance durations for the design are defined as follows: 1.50 µm < λ < 1.60 µm; 0 < z < 200 km. In the simulations un-chirped initial pulse with 5 ps as full-width at half-maximum is used. From Figure 28, it is clear that the zero dispersion wavelength is successfully set at λ0 which is equal to 1.55 µm. Furthermore, the dispersion curve becomes so flat by adding the Gaussian weighting term to the fitness function. In other words, in the absence of weighting function, the optimized dispersion has higher slope compared to its presence.
The impact of sigma parameter on the dispersion and its slope is illustrated in Figure 29 and 30. It is obvious that the dispersion value reduces by the sigma parameter increase in the predefined wavelength interval and the curve is the smoothest in the highest sigma case. This event can be described based on the fact that the weighting function (Gaussian function) has large values around the central wavelength by the increase in the Gaussian parameter. So, a large band of the wavelength around the central wavelength has almost the same chance for optimization and thus the dispersion will be small and uniform for this band. In other words, the duration of this band can be controlled by the sigma parameter.
The dispersion slope is strongly affected by the presence of σ in such a manner that its increase has the considerable influence on the dispersion slope and decreases it obviously. This result is easily visible in Figure 30 which shows the dispersion slope versus wavelength with variance of the σ as a parameter. According to the presented weighted function based GA optimization, the following optical and geometrical parameters are obtained. We find out that the optimal value of
type |
|
P | Q | R 1 | R 2 |
|
|
2.5462 | 0.6942 | 0.3046 | 7.0897 | -0.9517 | 4.886e-3 |
|
2.4450 | 0.7189 | 0.4049 | 3.0497 | -0.9854 | 8.159e-3 |
|
2.5763 | 0.8478 | 0.4774 | 2.2437 | -0.9877 | 7.178e-3 |
|
2.4374 | 0.8461 | 0.4098 | 1.7966 | -0.7076 | 8.064e-3 |
|
2.5914 | 0.8942 | 0.4728 | 2.2583 | -0.9780 | 6.812e-3 |
In Table 11 and 12 the simulated numerical values of dispersion and dispersion slope based on the presented algorithm for wavelength duration in
type |
|
|
|
|
|
-2.039 | -0.032e-3 | 2.002 | 4.041 |
|
-0.899 | 0.0238 | 0.765 | 1.664 |
|
-0.618 | 0.0032 | 0.363 | 0.981 |
|
-0.606 | 0.0051 | 0.346 | 0.952 |
|
-0.510 | 0.0034 | 0.191 | 0.701 |
type |
|
|
|
|
0.0416 | 0.0402 | 0.0401 |
|
0.0206 | 0.0164 | 0.0133 |
|
0.0153 | 0.0096 | 0.0488 |
|
0.0152 | 0.0093 | 0.0044 |
|
0.0137 | 0.0068 | 0.0007 |
As a final result of these simulations, we should point out that for zero value of the Gaussian parameter zero-dispersion wavelength has high accuracy compared to nonzero-values of the Gaussian parameters cases. By comparing presented results and the ones demonstrated earlier as a dispersion flattened optical fiber, it is clear that the least design has considerable band width, the band between to zero dispersion wavelength. Moreover, the dispersion value tolerance in this interval is so small which is a direct result of its small slope.
10. Conclusion
In this chapter some special fiber structures for covering broadband optical fiber communications were reviewed. For these three cases R, W and M two different types for each of them were considered and discussed in detail. We have been shown that using the proposed design method in this chapter systematic approach for broadband applications can be found and considering the fibers in this chapter ultra broadband communications are available.
References
- 1.
Kazumasa Ohsono Tomoyuki Nishio Takahiro Yamazaki Tomomi Onose Kotaro Tan 2000 “Low Non-linear Dispersion-shifted Fiber for DWDM Transmission”, HITACHI CABLE REVIW,19 . - 2.
Baghdadi J. A. Safaai-Jazi A. Hattori H. T. 2001 Optical“.fiberswith.lownonlinearity.lowpolarization-mode.dispersionfor.terabitcommunications,”.OpticsLaserTechnology.Vol33 285 291 ,. - 3.
Makouei S. savadi-Oskouei M. Rostami A. Kanani Z. D. K. 2002 “DISPERSION FLATTENED OPTICAL FIBER DESIGN FOR LARGE BANDWIDTH AND HIGH-SPEED OPTICAL COMMUNICATIONS USING OPTIMIZATION TECHNIQUE”, Progress In Electromagnetics Research B,13 21 40 , 2009. - 4.
Govind P. - 5.
Kato T. Hirano M. Tada A. Fukuada K. Fujii T. Ooishi T. Yokoyama Y. Yoshida M. Onishi M. 2002 “Dispersion flattened transmission line consisting of wide-band non-zero dispersion shifted fiber and dispersion compensating fiber module ”, ,8 231 239 ,. - 6.
Tewari R. Pal B. P. Das U. K. 1992 “Dispersion shifted dual shape core fibers: Optimization based on spot size definitions”, ,10 1 5 , - 7.
Dipankar Ghosh ,Debashri Ghosh Mousumi Basu 2008 “Designing a graded index depressed clad non-zero dispersion shifted optical fiber for wide band transmission system”, , vol.119 63 68 ,. - 8.
Mikkelsen B. Raybon G. Zhu B. Essiambre R. J. Bernasconi P. G. Dreyer K. Stulz L. W. Knudsen S. N. 2001 “High spectral efficiency (0.53 bit/s/Hz) WDM transmission of 160 Gb/s per wavelength over 400 km of Fiber ”, Technical Digest of OFC 2001,, Paper ThF2. - 9.
Ito T. Fukuchi K. Sekiya K. Ogasawara D. Ohhira R. Ono T. 2000 “6.4 Tb/s (160 × 40 Gb/s) WDM transmission experiment with 0.8 bit/s/Hz spectral efficiency ”, ECOC, 2000, Postdeadline paper PD1.1. - 10.
Kani J. Hattori K. Jinno M. Kanamori T. Oguchi K. 1999 “Tripple-wavelength-band WDM transmission over cascaded dispersion-shifted fibers ”, Technical Digest of OAA’99,, Paper WC2. - 11.
Fukuchi K. Kasamatsu T. Morie M. Ohhira R. Ito T. Sekiya K. Ogasawara D. Ono T. 2001 “10.92 -Tb/s (273 × 40-Gb/s) triple-band/ultra-dense WDM optical-repeatered transmission experiment ”, OFC 2001, Postdeadline paper PD24,. - 12.
Naka A. et al. Lightwave Technol. 12 12 2 February (1994)280 287 - 13.
Liu Y. et al. 1996 OFC’96, WK15, - 14.
Varshney R. K. Ghatak A. K. Goyal I. C. siny C. Antony 2003 “Design of a flat field fiber with very small dispersion slope ,” ,9 189 198 , - 15.
Sakamoto J. Kani J. Jinno M. Aisawa S. Fukui M. Yamada M. Oguchi K. 1998 “Wide wavelength band (1535-1560 nm and 1574-1600 nm), 28 × 10 Gbit/s WDM transmission over 320 km dispersion shifted fiber ”, Electron. Letter,34 392 394 , - 16.
Yashida S. Kuwano S. Iwashita K. 1996 “10 Gbit/s ×10 channel WDM transmission experiment over 1200 km with repeater spacing of 100 km without gain equalization or pre-emphasis ”, Optical Fiber Communication (OFC’) 96, San Jose, CA, TuD6, - 17.
Chraplyvy A. R. 1990 “Limitation on lightwave communications imposed by optical fiber nonlinearities ”, Lightwave Technology,8 1548 1557 , - 18.
Tkach R. W. Chraplyvy A. R. Fabrizio F. Gnauck A. H. Derosier R. M. 1995 “Four photon mixing and high speed WDM systems ”, Lightwave Technology,13 841 849 , - 19.
Akasaka Y. 1999 “New optical fibers for high bit rate and high capacity transmission ”, SPIE Proc.,3666 23 29 , - 20.
Liu Y. Antos A. J. Bhagavatula V. A. Newhouse M. A. 1995 Single mode dispersion shifted fiber with effective area larger than 80 μm2 and good bending performance ”, Proc. of ECOC’95, TuL2.4, - 21.
Tateda M. Kato Y. Seikai S. Uchida N. 1982 Design“.considerationon.singlemode.fiberparameters,”.PansI. E. C.IECE Japan, (in Japanese)J56-B 56 324 331 ,. - 22.
Kitayama K. I. Kato Y. Ohashi M. Ishida Y. Uchida N. 1983 “Design considerations for the structural optimization of a single mode fiber ,” Ligthwave Technology,1 363 369 , - 23.
Savadi-Oskouei M. Makouei S. Rostami A. Koozeh Z. D. Kanani 2007 “Proposal for optical fiber designs with ultrahigh effective area and small bending loss applicable to long haul communications” , ,46 6330 6339 , - 24.
Namihira Y. 1994 “Relationship between nonlinear effective area and mode-diameter for dispersion shifted fiber,” ,30 262 264 , - 25.
Hautakorpi M. Kaivola M. 2005 “Modal analysis of M-type dielectric-profile optical fibers in the weakly guiding approximation, ” Optical Society of America A.,22 1163 1169 ,. - 26.
Nunes F. D. de Souza C. A. 1996 Melo“.Theoreticalstudy.ofcoaxial.fibers,”Applied.Optics vol. 35388 398 ,. - 27.
Goldberg D. E. 1989 Genetic Algorithms in Search, Optimization and Machine Learning (Addison-Wesley, 1989). - 28.
Tian X. Zhang X. 2004 “ Dispersion flattened design of large effective area single mode fibers with ring index profiles”, ,230 105 113 , - 29.
Tseng Z. O. 1992 On“.theConvergence.ofthe.CoordinateDescent.Methodfor.ConvexDifferentiable.Minimization” vol72 7 35 ,. - 30.
Yin S. K.w Chung. H. Liu H. Kurtz P. Reichard K. 2000 “ A new design for non-zero dispersion shifted fiber (NZ-DSF) with a large effective area over 100µm2 and low bending and splice loss”, Optics Communications,177 225 232 , - 31.
Jiang X. Wang R. 2004 “Non-zero dispersion shifted optical fiber with ultra large effective area and low dispersion slope for terabit communication systems”, optics Communications,236 69 74 , - 32.
Kim B. H. Park Y. Kim D. Y. Paek U. C. Han W. T. 2002 in: OFC_2002, Technical Digest,,173 174 . - 33.
Pack U. C. Kurkjian C. R. 1975 58 330. - 34.
Matsuo S. Tanigawa S. Himeno K. Harada K. 2002 in: OFC_2002, Technical Digest,,329 330 . - 35.
Safaai-Jazi A. 1996 “Large effective area fibers: propagation properties and optimum designs”, SPIE Proc.,3666 30 39 , 1999. - 36.
Nouchi P. Sansonetti P. Von Wirth J. Le Sergent, C. - 37.
Silva V. 1998 new“. A.designfor.dispersionshifted.fiberwith.aneffective.arealarger.than 1002 μm2 and good bending characteristics,” Proc. of OFC’98, ThK1,301 . - 38.
Kumano N. Mukasa K. Matsushita S. Yagi T. 2002 “Zero dispersion slope NZ-DSF with ultra wide bandwidth over 300 nm”, Proc. of ECOC’2002, - 39.
Zhu B. Nelson L. E. Leng L. Stulz S. Knudsen S. Peckham D. 2002 “1.6 Tbits/s (40 × 42.7 Gbit/s) WDM transmission over 2400 Km of fiber with 100 Km dispersion managed spans”, ,38 647 648 ,. - 40.
Nouchi P. Sansonetti P. Von Wirth J. Le Sergent,” C. 1996 New dispersion shifted fiber with effective area larger than 90 μm2”, Proc. of ECOC’96, MoB3.2,. - 41.
Savadi-Oskouei M. Rostami A. Makouei,” S. 2009 A novel fiber design strategy for simultaneously introducing ultra small dispersion and dispersion slope using genetic algorithm”, European Transaction on Telecommunications,20 37 47 ,. - 42.
Kawakami S. Nishida S. 1974 “Anomalous dispersion of new doubly clad optical fibers,” ,10 38 40 ,. - 43.
Kawakami S. Nishida S. Sumi M. 1976 Transmission“.characteristicsof.W-typeoptical.fibers”Proc. Inst. Elec. Eng.,123 586 590 ,. - 44.
James B. Day C. R. 1986 “A review of single mode fibers with modified dispersion characteristics”, Lightwave Technology,4 967 979 ,