Wideband Gain Flattened Hybrid Erbiumdoped Fiber Amplifier/Fiber Raman Amplifier
 Author: Afkhami Hossein, Mowla Alireza, Granpayeh Nosrat, Hormozi Azadeh Rastegari
 Organization: Afkhami Hossein; Mowla Alireza; Granpayeh Nosrat; Hormozi Azadeh Rastegari
 Publish: Journal of the Optical Society of Korea Volume 14, Issue4, p342~350, 25 Dec 2010

ABSTRACT
An optimal wideband gain flattened hybrid erbiumdoped fiber amplifier/fiber Raman amplifier (EDFA/FRA) has been introduced. A new and effective optimization method called particle swarm optimization (PSO) is employed to find the optimized parameters of the EDFA/FRA. Numerous parameters which are the parameters of the erbiumdoped fiber amplifier (EDFA) and the fiber Raman amplifier (FRA)define the gain spectrum of a hybrid EDFA/FRA. Here, we optimize the length, Er^{3+} concentration, and pump power and wavelength of the EDFA and also pump powers and wavelengths of the FRA to obtain the flattest operating gain spectrum. Hybrid EDFA/FRA with 6pumped and 10pumpedFRAs have been studied. Gain spectrum variations are 1.392 and 1.043 dB for the 6pumped and 10pumpedFRAs,respectively, in the 108.5 km hybrid EDFA/FRAs, with 1 mW of input signal powers. Dense wavelength division multiplexing (DWDM) system with 60 signal channels in the wavelength range of 1529.21627.1 nm, i.e. the wide bandwidth of 98 nm, is studied. In this work, we have added FRA’s pump wavelengths to the optimization parameters to obtain better results in comparison with the results presented in our previous works.

KEYWORD
Hybrid fiber amplifier , Erbiumdoped fiber amplifier , Fiber Raman amplifier , Optimization , Particle swarm optimization , (060.2320) Fiber optics amplifiers and oscillators , (140.4480) Optical amplifiers , (060.0060) Fiber optics and optical communications , (060.2410) Fibers , erbium

I. INTRODUCTION
After the invention of EDFAs, data transfer capacities of optical communication systems exploded based on wavelength division multiplexing (WDM) systems [1]. Then more demands for higher capacity systems, in addition to the progress in optical device technology, encouraged the researchers to focus on the multiwavelength pump FRAs which have broader bandwidths [2]. Hybrid EDFA/FRAs could have the benefits of both EDFAs and FRAs, hence some efforts have been carried on to offer efficient constructions for hybrid EDFA/FRAs [3]. Hybrid EDFA/FRAs have long lengths and consequently lower numbers of link repeaters are needed. To build these hybrid amplifiers, low price EDFAs are used. Their costs are relatively lower than FRAs and their noise figures (NFs) are lower than those of the EDFAs [2]. A WDM system with hybrid EDFA/FRAs has lots of elements and the combination of these elements and the value of their parameters are of great importance. The module combination is the art of engineering work. But choosing the module parameters could be difficult,because these parameters affect the characteristics of the hybrid EDFA/FRAs. Here, a global optimization method (GOM)is quite necessary. GOMs help us to find the optimum values for parameters of the modules. PSO is a GOM that is effective and fast and can deal with a high number of variables [4, 5]. Using PSO, the parameters of the EDFA and FRA used in hybrid EDFA/FRA have been optimized to obtain a wideband gain flattened amplifier. In this paper,hybrid EDFA/FRA configurations are introduced, the method of simulation of hybrid EDFA/FRA is discussed and then PSO is presented and the results of optimization are given for different pumping schemes and configurations. In addition to optimization of the length, Er^{3+} concentration, and pump power and wavelength of the EDFA and also pump powers of the FRA all of which have been carried on in our previous work [6], the pump wavelengths of the FRA are optimized in order to have a flattergain hybrid EDFA/FRA. Adding the pump wavelengths of the FRA to the other optimization variables increases the number of optimization variables to 24 in the 10 pumped FRA hybrid EDFA/FRA.Dealing with this high number of variables is quite a challenge. Here PSO shows its strength to manipulate a large number of variables in a fast and straightforward manner in comparison with the other optimization methods such as genetic algorithm (GA) [714]. The results are discussed in section V.
II. HYBRID EDFA/FRA CONFIGURATION
FIG. 1 depicts the configuration which we have chosen for the hybrid EDFA/FRA. There are 60 signal channels with frequencies between 184.25 and 196.05 THz, placed equally in 200 GHz of frequency distances launched to the single mode fiber (SMF). SMF allows the propagation of only one mode, retains the fidelity of each light ray and has higher bandwidths over longer links. TABLE 1 presents the characteristics of the designed WDM system. If we utilize communication systems with the rate of 40 Gbit/s at each of these channels, the capacity of this WDM link could reach 2.4 Tbit/s. Choosing the values of the input signal powers is of great importance and will affect the gain spectrum variations of the amplifier. We select the input signal powers to be 1 or 10 mW to have a comparison
between the results. Here, SMF acts as the transmission medium. A dispersion compensating fiber (DCF) is located after the SMF to compensate the dispersion that took place in the SMF. DCF has a high negative value of dispersion which is supposed to neutralize the SMF lower positive value of dispersion. So, shorter length of DCF compared to the length of SMF is needed and the length of DCF must be chosen appropriately. Afterward, a dispersion shifted fiber (DSF) is placed which has zero dispersion at 1550 nm wavelength and is pumped by 6 or 10 backward pump LD sources which acts as an FRA. FRA pump powers are propagated counter to the direction of the signal powers. In this paper, we have chosen 2 different sets of fiber lengths. Accordingly, SMF, DCF, and DSF could be 50, 8.5, and 50 km in the first configuration and 80, 14,and 80 km in the second configuration, respectively. An isolator (ISO) is laid between the DCF and DSF to suppress the residue of the FRA backward pumps and noise, raised in the DSF. A WDM coupler is set up to couple the power of pumps to the DSF in the counterpropagating direction. Pump powers and wavelengths are chosen through the optimization process discussed in section IV. Next module is a Cband EDFA. Principal range of amplification in the Cband EDFA is 15301570 nm and its gain falls at wavelengths over 1570 nm. On the other hand, the amplification at the higher range of 15701630 takes place mostly through the DSF as an FRA. TABLE 2 lists the characteristics of different modules used in our hybrid EDFA/FRA configurations.Lots of variables exist in the design of EDFA. In this paper we have chosen 4 of the most important variables to be optimized, the length, the Er^{3+} concentration, the pump power and wavelength of the EDFA. Finally, we have a detector at the end of the link span to receive the output powers. Gain spectrum variations will be computed at this point. This span can be repeated if longer links are needed.As indicated, four different configurations have been proposed.These configurations are 108.5 or 174 km in length and each could adopt 6 or 10pumpedFRA. Additionally, the amplifier characteristics are calculated in the cases of different input signal powers of 1 or 10 mW for each channel. As discussed in section V, gain spectrum variations of 1.794 dB can be achieved for the 174 km link span,10pumpedFRA hybrid EDFA/FRA with 1 mW of input signal power.
III. HYBRID EDFA/FRA SIMULATION
Simulation of EDFA and FRA should be done separately.Also, attenuation, dispersion and noise of the other components such as SMF, DCF, DSF, erbiumdoped fiber (EDF), and ISO must be taken into account. Therefore, the simulation of hybrid EDFA/FRA has to be a combination of simulations of different modules and components of the hybrid amplifier.In the following sections we discuss the simulation of EDFAs and FRAs separately.
1. Simulation of the EDFA : Giles model [15] is the most common model used for the optical amplification that occurs in a silica fiber doped by Er^{3+}. We assume a twolevel,homogeneously broadened model for Er^{3+} that is sufficient to have a simple and accurate EDFA model [15, 16]. In this model ^{4}I_{15/2} is the ground level and ^{4}I_{13/2} is the excited level. We suppose that all the radiative transitions occur between these two levels. Computer analysis and simulation is discrete, so we should break the continuous spectrum of the optical waves into discrete optical beams with frequency bandwidth of Δv_{k} and optical power ofP_{k} centered at frequencyν_{k} or wavelengthλ_{k} . We should solve an ordinary differential equation for each one of the optical beams to find the propagation of optical powers along the fiber. Therefore,N ordinary differential equations should be solved, whereN is the number of the pump and signal channels. Rate equations of the excited level of Er^{3+} ions are computed as a summation of the emission and absorption contributions of all theN optical beams [17]. We have used Al/Psilica fiber as the EDFA, the emission and absorption cross sections of which are depicted in FIG.2(a) [18]. All fibers are assumed to be SMFs, hence theV number should be less than 2.405. Also we have assumed the steady state condition of∂ N_{2}(r, z, t)/∂t = 0 M whereN_{2} is the upper level population, and also the weakly guiding approximation of？ = (n_{1} n_{2})/n_{2} << 1 wheren_{1} andn_{2} are the refractive indices of the fiber core and cladding,respectively. The flow of signal and pump power through EDFA is formulated by ordinary differential equations[17]. Rate equations in the steadystate condition have to be solved to acquire the metastable level population of Er^{3+} ions [17].Since the EDFA pump is in the copropagating direction,we need to solve an initial value problem. We ignore the counter propagating pumps and eliminate the backward amplified spontaneous emission (ASE) from our simulation to avoid solving the boundary value problem which is a bit more challenging. In this paper, a onestep predictorcorrector method [19] is used to simulate the EDFA.Using this method, the propagated powers at channels are calculated step by step using their values at the previous steps and their predicted values at the next steps [20].
The pumping wavelength is around 1480 nm. Step index fiber with uniform distribution of Er^{3+} ions are used in this simulation. EDFA characteristics which have been used are listed in TABLE 3.
2. Simulation of the FRA : Coupled ordinary differential equations govern the propagation of the average optical power of pumps, signals, ASEs, and Rayleigh backscattering beams in a FRA [19, 21]. For copropagating pumps, an initial value problem must be solved and for the case of counterpropagating or bidirectional pumps, twopoint boundary value problems must be solved. Here, we have counterpropagating pumps, hence a twopoint boundary value problem has to be solved. Fourstep predictorcorrector method based on the AdamsBashforthMoulton formula [19] and shooting algorithm [22] are used to solve this problem numerically[23]. FIG. 2(b) shows the Raman gain efficiency of the DSF [2] and FIG. 2(c) shows attenuation spectra of the DSF [2], SMF, and DCF [24] used in our simulation.We have used C^{++} programming to model the hybrid EDFA/FRA. The time needed to accomplish the simulation, depends
on the configuration of the setup. It could be 0.2 to 1 second when run by an AMD Athlon 2.71 GHz CPU with 1.00 Gb of RAM.
IV. PARTICLE SWARM OPTIMIZATION IMPLEMENTATION
PSO is a populationbased stochastic method for solving optimization problems which is classified as a swarm intelligence technique. In this method, a population of particles which are bunches of variables that present candidate solutions,will be selected randomly from the possible solution space and these particles will move according to their social learning from each other toward the optimum solution.The rules and constraints which lead the flying particles in the search space toward the global optimized solution are originated from the behavioral models of bird flocks [4,5]. This method was introduced and developed by Kennedy and Eberhart [25, 26].
Let us assume that F is the fitness function. There are n particles P_{i} where
i = 1, …, n, x_{i} is the current position ofP_{i}, Gb is the global best which is the best fitness value that has ever happened,Pb_{i} is the personal best ofP_{i} which is the best previous fitness value ofP_{i} ,x_{Gb} andx_{Pbi} are the positions of theGb andPb_{i} , respectively,v_{i} is the velocity ofP_{i} andk is the time interval. We can present PSO method in the following steps:Step 1. Pis are randomly chosen from the search space to form the swarm P(k)
Step 2. Calculate the fitness values F(xi(k)) for each Pi.
Step 3. If F(xi(k))>Pbi, then Pbi=F(xi(k)) and xPbi=xi(k).
Step 4. If F(xi(k))max>Gb, then Gb=F(xi(k))max and xGb=xi(k).
Step 5. Set the new velocities
where
φ is the inertia weight, k is the iteration number and kmax=300 is the total number of iterations,C_{f} =2 is the accelerator constant.ρ_{1} andρ_{2} are uniform random values between 0 and 1.Step 6. Set the positions
Step 7. Check the convergence criterion, if fulfilled, the optimization is over, if not go to step 2.
To perform the optimization, PSO acts as the main program and the simulation of the hybrid EDFA/FRA will be the subprogram. Fitness function is the maximum signal network gain deviation from zero which is supposed to be minimum,formulated in the following:
where
P_{s} ^{k}(0) andP_{s} ^{k}(L) are thek th input and output signal powers, respectively. Taking into account the designcriteria, variables of optimization and their ranges in this work have been listed in TABLE 4. As it was mentioned,we used C^{++} programming to simulate the hybrid EDFA/FRA and execute the optimization process. The necessary time to accomplish one round of optimization depends on the configuration which is to be optimized and could be between 2 and 3 hours, by an AMD Athlon 2.71 GHz CPU with 1.00 Gb of RAM. The number of particles and iterations in our PSO model are 40 and 300, respectively.
V. RESULTS AND DISCUSSIONS
In this paper we have simulated 4 different configurations for hybrid EDFA/FRA. Two configurations for 108.5 km and two for 174 km link with 6 and 10pumpedFRA.Also, we survey the results for both input signal powers of 1 or 10 mW.
A group of optimization variables or the group of simulation inputs which are introduced in TABLE 4 are known as a particle
P_{i} in the PSO. We haven=40 particles in our optimization process and it means that we have 40 groups of optimization variables which are randomly spread in our optimization space. Each one of the optimization variables makes a dimension of the optimization space. Particle positionsx_{i} are arrays which are composed of the optimization variable values. We also haven=40 particle positions xi in our optimization process that need to be optimized.At each iteration, we insert the inputs which are the optimization variable values or the particle positions to the main hybrid EDFA/FRA simulation computer program or fitness function. Then the simulation will result in the output or the fitness valueF(x_{i}) which is the maximum signal network gain deviation from zero. We consider the deviation from zero to take into account both the negative and positive gain values. According to the PSO steps which are introduced in the previous section, new values for optimization variables are selected at the end of each iteration. We have 24 and 16 optimization variables for 10 and 6 pumped FRA hybrid amplifiers. So, the particle positions are composed of 24 or 16 elements.To reach the flattest gain spectrum for the hybrid EDFA/FRA, it is better to set the EDFA pump wavelength to have a Cband EDFA. Because of the shape of the EDFA gain spectrum, it is preferred to have amplification at lower wavelengths by EDFA and at higher wavelengths by FRA. In the optimization process, the EDFA pump power is inclined to attain high values. So, it is better to choose 1480 nm pumping wavelength for the EDFA,although it has higher noise compared to 980 nm ones.Also, in high power schemes it is more efficient to use 1480 nm pumping wavelength. The radius of the Er^{3+}doped area at EDFA is equal to the radius of the fiber used as EDFA that is 6.5 ㎛. An EDFA with a small core radius will result in better gain spectrum characteristics,but its fabrication is hard and it is more susceptible to damage in facing high EDFA pump power. Also a small core radius EDFA is rare in the market. Er^{3+} ion concentration has to be chosen in a confined area, since low Er^{3+}concentration would elongate the EDFA and consequently create higher variation of the gain spectrum. Also, high Er^{3+}concentration will increase unfavorable intra ion interactions[1]. The main source of noise of EDFAs is ASE. The level of noise in FRAs is lower than that in EDFAs. In this paper, we have assumed an acceptably low level for noise, and to speed up the optimization process, we have not simulated noises. Lower input signal power in EDFAs will result in higher and flatter gain spectra, but input signal power should not be so low that it could not be retrieved from the noise. The designed amplifier is a distributed amplifier, and it is better to define net gain for it. We prefer the signal powers at the end of the link to be the same as their values at the beginning of the link.So, the ideal net gain of our hybrid EDFA/FRA is zero.We have assumed the EDFA to be polarization independent as it almost is, but we consider the polarization dependency of FRA by including a factor in our simulations.
Amplification in the FRA is due to the nonlinear phenomenon of stimulated Raman scattering (SRS). By using a high number of pumps in FRAs with different wavelengths and powers, the gain spectrum is tailored and the EDFA gain variation is compensated. The amplification band can be chosen simply by adjusting the pump wavelengths at appropriate distances. The noise transition from pumps to signals occurs less in the counterpropagating pump scheme, hence it is more advantageous to place the pumps at the end of the FRA and increase the length of link span. We have to appoint the input signal powers to somehow avoid harmful interference with noise due to the low level of the input signal power, and unfavorable nonlinear effects due to the high level of the input signal powers.
All the optimization results, obtained by PSO, for different schemes and configurations are given in TABLES 58. It is obvious that gain spectrum variations will be lower if we choose 10pumped configurations instead of 6pumped ones, or if we choose 108.5 km links instead of 174 km ones, and also if we choose 1mW of input signal power instead of 10 mW. As listed in TABLE 5, gain spectrum variations of 1.392 and 1.043 dB are attained for 108.5
km, 6 and 10pumpedFRA hybrid EDFA/FRA, respectively with 1mW input signal powers. Also, as listed in TABLE 6, gain spectrum variations of 1.864 and 1.794 dB are attained for 174 km, 6 and 10pumpedFRA hybrid EDFA/FRA, respectively with 1mW input signal powers.
Using the results listed in the TABLES 5 and 6, optimized signal gain spectra of 108.5 and 174 km 10pumpedFRA hybrid EDFA/FRA are shown in FIG. 3(a), 108.5 and 174 km 6pumpedFRA hybrid EDFA/FRA are depicted in FIG. 3(b), 6 and 10pumpedFRA 108.5 km hybrid EDFA/FRA
is pictured in FIG. 4(a), and 6 and 10pumpedFRA 174 km hybrid EDFA/FRA is pictured in FIG. 4(b), all for 1 mW input signal powers. These figures present good visual comparison of the optimized gain spectra of different schemes.
It is common to select the FRA pump wavelengths in an equidistant form, but here we search the FRA pump
wavelengths through an optimization process which will result in better outcomes for gain spectrum variations. For example, gain spectrum variations of 2.91 dB for a 130.5 km, 6pumpedFRA hybrid EDFA/FRA, equidistant FRA pump wavelengths [6] has become 1.864 dB for 174 km,6pumpedFRA hybrid EDFA/FRA, by optimization of FRA pump wavelengths and also gain spectrum variations of 2.03 dB for a 130.5 km, 10pumpedFRA hybrid EDFA/FRA, by equidistant FRA pump wavelengths [6] has become 1.794 dB for 174 km, 10pumpedFRA hybrid EDFA/FRA, by optimization of FRA pump wavelengths,which are notable improvements.
For most cases, the wavelength spectral distance between the FRA pump wavelengths, presented inTABLES 58,are longer than 8 nm. At the closest point it will reach to 2.3 nm. By using narrow spectral width laser diodes at the places where the pump wavelengths are near the adjacent pump wavelengths, the pump powers will not interfere with each other and sound operation is guaranteed. Typical spectral width of a laser diode is about 0.110 nm.
FIG. 5 shows a comparison between the convergence speeds of PSO and GA optimization processes which are carried on the optimization of one of the configurations of our WDM systems. It is apparent that PSO acts faster and more efficiently in searching for the optimum values for a system under test.
VI. CONCLUSION
In this paper, we have used an effective optimization method, called particle swarm optimization (PSO) to offer optimum constructions for hybrid erbiumdoped fiber amplifier/fiber Raman amplifiers (EDFA/FRAs). PSO allowed us to include a high number of parameters in the optimization process, so that we could optimize as much as 24 parameters simultaneously. Using the efficient PSO algorithm, pump wavelengths of FRA are optimized in addition to the FRA pump powers, EDFA pump power and wavelength, Er^{3+}concentration, and EDFA length. This results in obtaining hybrid amplifiers with gain spectrum variations of as low as 1.043 dB for a 108.5 km and 1.794 dB for a 174 km,10pumpedFRA hybrid EDFA/FRA, respectively. A comparison between the convergence speeds of PSO and genetic algorithm is done that shows that PSO is faster than genetic algorithm in converging toward the best values of amplifier parameters.

[FIG. 1.] Configurations of the hybrid EDFA/FRA with 6 or10pumpedFRA in backward direction.

[TABLE 1.] Characteristics of the WDM transmission system

[TABLE 2.] Characteristics of hybrid EDFA/FRA fiber modules

[TABLE 3.] EDFA and FRA parameters used in the hybrid amplifier modeling

[FIG. 2.] Characteristics of the fibers (a) emission andabsorption cross sections of Er+3 doped in Al/Psilica SMF (b)Raman gain efficiency of DSF (c) DCF SMF and DSFattenuation spectrum. [2 23 24]

[TABLE 4.] Variables of the optimization of the EDFA/FRAsand their ranges

[TABLE 5.] Values of the optimized parameters of 6 & 10pumpedFRA 108.5 km hybrid EDFA/FRA with 1mW of input signalpower obtained by PSO

[TABLE 6.] Values of the optimized parameters of 6 & 10pumpedFRA 174 km hybrid EDFA/FRA with 1mW of input signal powerobtained by PSO

[FIG. 3.] Optimum signal gain spectra of (a) 108.5 & 174 km10pumpedFRA hybrid EDFA/FRA and (b) 108.5 & 174 km6pumpedFRA hybrid EDFA/FRA (both for input signalpower of 1mW).

[FIG. 4.] Optimum signal gain spectra of (a) 6 &10pumpedFRA 108.5 km hybrid EDFA/FRA and (b) 6 &10pumpedFRA174 km hybrid EDFA/FRA (both for input signal power of1mW).

[TABLE 7.] Values of the optimized parameters of 6 & 10pumpedFRA 108.5 km hybrid EDFA/FRA with 10mW of input signalpower obtained by PSO

[TABLE 8.] Values of the optimized parameters of 6 & 10pumpedFRA 174 km hybrid EDFA/FRA with 10mW of input signalpower obtained by PSO

[FIG. 5.] Comparison of the convergence speeds of the PSOand GA toward the optimum values versus the normalizedexecuting time.