Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Optimization criteria and design of few-mode Erbium-doped fibers for cladding-pumped amplifiers

Open Access Open Access

Abstract

We propose a novel optimization method that combines two design criteria to reduce the differential modal gain (DMG) in few-mode cladding-pumped erbium-doped fiber amplifiers (FM-EDFAs). In addition to the standard criterion that considers the mode intensity and dopant profile overlap, we introduce a second criterion that ensures that all doped regions have the same saturation behavior. With these two criteria, we define a figure-of-merit (FOM) that allows the design of FM-EDFAs with low DMG without high computational cost. We illustrate this method with the design of six-mode erbium-doped fibers (EDFs) for amplification over the C-Band targeting designs that are compatible with standard fabrication processes. The fibers have either a step-index or a staircase refractive index profile (RIP), with two ring-shaped erbium-doped regions in the core. With a staircase RIP, a fiber length of 29 m and 20 W of pump power injected in the cladding, our best design leads to a minimum gain of 22.6 dB while maintaining a DMGmax under 0.18 dB. We further show that the FOM optimization achieves a robust design with low DMG over a wide range of variations in signal power, pump power and fiber length.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

With the unrelenting need to transfer large amounts of data, optical fiber communication networks have experienced tremendous growth over the past decade. As we approach the fundamental Shannon limit of standard single-mode fiber transmission, space-division multiplexing (SDM) technology is increasingly considered for delivering high-capacity [13]. SDM can be accomplished either by increasing the number of cores in a single fiber strand or by multiplying channels through transmission on different modes of a few-mode fiber (FMF), the latter technique being also referred to as mode division multiplexing (MDM) [46]. To support long haul transmission [7], FMFs require amplification to compensate for the inevitable background loss of silica-based fibers. Few-mode SDM communication systems must thus be able to rely on optical amplifiers that deliver high performance for all channels.

Although erbium-doped fiber amplifiers (EDFAs) are a well-established optical amplification technology [8,9], the design of SDM amplifiers is challenging. Such amplifiers must be able to simultaneously amplify all modes with high gain, while providing low differential modal gain (DMG) and low noise figure (NF). The architecture of a few-mode erbium-doped fiber amplifier (FM-EDFA) is akin to that of a single-mode EDFA, where the signal to be amplified is injected into the doped core of the erbium-doped fiber (EDF). However, in FM-EDFAs, the gain must be balanced across both signal wavelengths and modes, this means that both differential spectral gain (DSG) and DMG should be considered when designing a FM-EDFA [10]. To achieve these goals, the design of few-mode erbium-doped fibers (FM-EDFs) entails the optimization of two key parameters, namely the erbium doping profile and the refractive index profile.

With uniform doping in the core, amplifier performance depends on the overlap of the signal mode intensity and pump mode intensity, as determined by the refractive index profile and injection conditions. Some amplifier designs optimize the pump modal content to shape the pump transverse intensity distribution and produce a population inversion in the EDF that renders the small signal gain equal for the various signal modes. This is possible if the amplified mode closely matches the pump intensity profile [1113]. This approach is difficult to implement as the coupling efficiency of the pump to the various modes must be precisely controlled. Selective excitation of pump modes is challenging and complex to implement experimentally. For example, spatial light modulators (SLM) can be used but they have low damage threshold. Consequently, in this paper, we consider cladding pumping that can be implemented using a high-power low-cost multimode pump laser [14].

Tailoring of the refractive index profile (RIP) of the fiber to adjust the signal mode field distribution is the second approach proposed in the literature to minimize DMG. For example, the FM-EDF proposed in [15] has a dip at the center of its RIP that is engineered to modify the intensity distribution of the LP01 mode. By adjusting the depth of the refractive index dip, the LP01 mode profile resembles more closely the LP11 mode and LP21 profiles, which minimizes the DMG between these spatial modes. The first experimental demonstration of a FM-EDFA was accomplished with a fiber having a simple uniformly doped core, yielding a DMG between 5 and 10 dB [15]. Subsequently, Wada et al. demonstrated an EDF with a dip in the centre of its RIP and obtained a gain, in the L-band, of about 20 dB for five modes with a 3 dB DMG using an LP01 pump mode [16]. This approach can lead to substantial splice loss unless mode-field adaptors are used. Pushing this approach further leads to the design of a ring-core fiber in which the central dip reaches a refractive index that is equal, or even lower, than that of the cladding. Since all the signal spatial modes will have a similar overlap factor with the erbium-doped core, the ring-core fiber amplifier provides similar gain for all guided modes [1722].

In another approach, several interesting fibers with tailored erbium doping profiles (EDP) were proposed to reduce the DMG of FM-EDFAs. A fiber with central dips in both the EDP and the RIP was proposed and fabricated by Jung et al. [15]. In this fiber, the erbium ion concentration increased from the center to the edge of the fiber core. Wakayama et al. further extended the doped region in the cladding to provide more gain to the higher order fiber modes when using cladding pumping, which decreased DMG to 3.3 dB. Such a cladding-pumped six-mode EDFA was experimentally demonstrated and a DMG of less than 3.3 dB was obtained over the C-band [23]. Similarly, by intentionally oversizing the fiber core to allow a greater number of modes than required, the desired modes are more confined in the doped core, which lowers the DMG. In [24], a 10-mode EDFA was demonstrated with a fiber having a 24 µm core diameter and a 2 dB DMG was reported. This method also requires mode field adaptors to reduce splice loss.

In most of this previous work, the criteria used in the fiber design optimization targeted equal overlap of the signal modes with the EDP, which amounts to equalizing the small signal gain. When the EDFA is operated in saturation, as in a fully loaded scenario, this condition is no longer valid. To find a more general solution, other optimization techniques, such as genetic algorithms, can be used. If all possible variations of the erbium doping profile and refractive index of the fiber are examined, an excessively large number of possible designs must be simulated. The spectral gain and noise figure (NF) of each mode are usually computed using coupled differential equations accounting for the erbium population rate equations of a two-level model, and the pump, signal and amplified spontaneous emission (ASE) mode intensity propagation. These coupled equations are solved using algorithms such as the 4th order Runge-Kutta method. In addition to the high computational costs when examining all possible designs, we do not gain physical insights [2527]. For example, a one-ring EDP design was optimized to minimize the DMG in a three-mode EDF using a full propagation model [25]. In [26], a genetic algorithm is proposed to examine an enlarged two-layer EDP of a cladding-pumped fiber. The algorithm minimized the difference in the signal mode overlap with the doped regions and, simulations of the resulting 6-mode cladding pumped EDF design predicted an average gain of 24.8 dB and a DMG of less than 0.64 dB over the C-Band.

In this paper, we present a novel design method that is based on a novel figure of merit (FOM) combining two optimization criteria to achieve low DMG in cladding-pumped FM-EDFAs that are fully loaded. This FOM requires only the computation of the mode intensity profile, the total intensity profile, and their overlap with the EDP. This significantly decreases the total computational time required to find a solution, as compared to a Runge-Kutta method-based solution. Our method allows exploration of a large number of fiber designs. Furthermore, the low DMG property of these designs is robust with respect to variations in pump power, fiber length and input signal power.

The paper is organized as follows: in section 2, we present the optimization criteria and method. We start with the standard mode intensity-EDP overlap criteria to which we add a second criteria that takes saturation of the EDP into account. We then detail the simulation method, which is based on generating RIP, calculating mode intensity profiles, computing overlaps with different EDP, and finally evaluating the FOM. In section 3, we apply the method to a step-index refractive index profile and, in section 4, we compare the results to a staircase RIP with a raised refractive index region near the center of the fiber core. Finally, we examine the robustness of both designs in section 5 and a conclusion follows.

2. Optimization criteria and design method

To reduce the computation time, we propose a FOM combining two criteria that do not require solution of the full simulation model of EDF. The optimization of these two criteria is necessary to achieve low DMG in cladding-pumped FM-EDFAs. We assume that cladding pumping results in a uniform pump intensity distribution over the fiber cross-section. The cladding diameter used for all the designs presented in this paper is the standard 125 µm that is mechanically and optically compatible with current fiber technologies. Under the assumption of uniform power distribution, the required pump power directly scales with the cladding-to-core area ratio. Uniform pump power distribution is a usual assumption in the modeling of cladding pumped amplifiers [2830] that was validated by comparing experimental measurements to simulation results in the case of a 100 µm cladding in [30]. Thus, as long as the assumption holds, it is possible to reduce the necessary pump power by reducing the cladding diameter, e.g. the pump power would be reduced by a factor of nearly two for a 90 µm cladding diameter. In the examples we examine, we restrict our demonstration to designs that are simple and compatible with EDF fabrication that uses the modified chemical vapor deposition technique (MCVD) with solution doping. Firstly, we present the standard criterion and show that two fibers meeting this criterion can have very different DMG. We then introduce the second criterion and discuss the physical insight on which it is based. Lastly, we combine these criteria in a new FOM and present the steps in the numerical optimization technique. Throughout our discussion, the DMG is defined as

$$ DMG(\lambda )= \mathop {\max }\limits_{m{\; }\epsilon {\; }M} ({G({m_i},\lambda } )- G({m_j},\lambda ),$$
where the gain G, in dB, is a function of the wavelength λ and of the LP mode order mj $\epsilon {\; }M$; the maximization is taken over modes with ji. The DMG represents the largest modal gain difference at each wavelength. The value DMGmax is defined as the maximum of DMG(λ) over the C-band. With our optimization, we aim to minimize DMGmax.

2.1 First optimization criterion

Equalizing modal gain in FM-EDFAs is commonly based on the calculation of the overlap between each mode and the doped regions [26] and, as such, it corresponds to optimizing the design under a small signal gain condition. Dividing the fiber core into rings, as represented in Fig. 1, the criterion is defined as

$$ {Z_1} = \mathop {\max }\limits_{{\lambda _i}} \left( {\frac{{\mathop {\max }\limits_{m \in M} (\sum _{k = 1}^K {{\Gamma }_{s,(\textrm{m},{\lambda _i},k)}}{\rho _{Er,(k)}})}}{{\mathop {\min }\limits_{m \in M} (\sum _{k = 1}^K {{\Gamma }_{s,(\textrm{m},{\lambda _i},k)}}{\rho _{Er,(k)}})}}} \right),$$
where ${\varGamma _{s,({m,{\lambda_i},k} )}}$ is the distribution of the mth mode intensity over a specific region k of the fiber, λi is the wavelength and ${\rho _{Er,(k )}}$ is the erbium ion concentration over the same k region of the fiber. In the present case, each k region is defined as a 0.1 µm thick ring (annular regions), in the fiber core.

 figure: Fig. 1.

Fig. 1. Fiber core divided into K rings (annular regions), each ring k characterized by its erbium ion concentration ${\rho _{Er,(k )}}$ and mode intensity distribution ${\mathrm{\Gamma }_{s,({\textrm{m},{\mathrm{\lambda }_i},k} )}}$.

Download Full Size | PDF

Minimizing Z1 ensures that the different signal modes have access to a similar quantity of erbium ions for amplification, and therefore that a low DMG can be obtained. However, when used alone, this criterion is not sufficient since it does not consider the fact that the saturation of the erbium ion population in each ring is a function of the total signal intensity in that ring. To illustrate that this criterion alone is not sufficient, we examine the EDP profiles EDF-A and EDF-B shown in Fig. 2. The fibers have the same step-index RIP but two different EDP. The fiber mode profiles, simulated with COMSOL, considering the step-index RIP with a core index, ncore, of 1.449, a numerical aperture, NA, of 0.12 and a core radius of 10 µm, are displayed in Fig. 2. The graph also shows the total signal intensity, Ik, defined as the sum of the intensity of all modes in a given region k.

 figure: Fig. 2.

Fig. 2. The mode profile distributions and the total signal intensity of EDF-A and EDF-B across the fiber core. The green solid and dashed lines are the two-ring erbium doping profile of EDF-A and EDF-B respectively. The inset shows the RIP used for both designs.

Download Full Size | PDF

Now, calculating Z1 for both profiles, we obtain a value of 1.012 for EDF-A and 1.022 for EDF-B, indicating that EDF-A should provide better gain equalization between the modes considering the more equal overlap of each mode with the erbium ion distribution. We now simulate both EDF using the full simulation model with the two-level erbium ion model and propagation equations that are solved with a 4th Runge-Kutta method. The simulations consider the radial dependency of the signal power and population inversion across the fiber core as indicated by the summation over k in the equations presented in [27]. In the simulations, we neglect clustering and use typical absorption and emission cross-sections of aluminosilicate erbium doped fibers. We also consider that the LP11a, LP11b, LP21a and LP21b are excited with equal power and propagate independently, i.e. the coupling between the modes is neglected. These mode pairs thus result in a power distribution with no azimuthal dependence. This model is described in more detail in [27] and the simulation parameters can be found in Table 1. In this configuration, the EDFA is cladding-pumped and we use 36 wavelength channels per mode, which amounts to a total input power of -9 dBm per mode or -1.22 dBm total input signal power. That is, we are examining performance outside the small signal gain region. The resulting gain and DMG, calculated over the C-band, are shown in Fig. 3. and Fig. 4. Results reveal that EDF-B has much lower DMG, with DMGmax = 0.464 dB, than EDF-A, with DMGmax = 2.065 dB. EDF-A and EDF-B respectively display a maximum noise figure of 3.96 dB and 4.12 dB over the C-band. In these simulations, the optimal length is the fiber length for which the minimal gain over the C-band for all the modes is the highest. We used the minimum gain to determine the optimal fiber length since, when optimizing for gain flatness, the process would sometimes lead to fiber designs with low DMG accompanied by very low gain values.

 figure: Fig. 3.

Fig. 3. Gain (left) and DMG (right) over the C-band in EDF-A for the six modes LP01, LP11a, LP11b, LP21a, LP21b and LP02 (P­p = 20 W, Ps = -9 dBm/mode, fiber length = 23 m).

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. Gain (left) and DMG (right) over the C-band in EDF-B for the six modes LP01, LP11a, LP11b, LP21a, LP21b and LP02 (P­p = 20 W, Ps = -9 dBm/mode, fiber length = 30 m).

Download Full Size | PDF

Tables Icon

Table 1. FM-EDFA simulation parameters for evaluation of gain, DMG and noise figure

The difference in DMG observed between the two FM-EDF can be understood when examining the distribution of the total signal intensity in the core as represented by the red curve in Fig. 2. In the case of EDF-A, there is an important difference between the total signal intensity overlapping with the two doped regions. This means that, as the signal propagates along the doped fiber, the two doped rings will have different inversion levels. Because the overlap of each mode with the two rings are different, the modes will experience different gain saturation, resulting in an increased DMG with respect to the small signal criterion. In the case of EDF-B, the total signal intensity at the position of both doped rings is similar and, consequently, the doped rings will experience the same level of gain saturation. This means that all modes will also saturate similarly since how the mode overlap with the doped region is divided between the two rings no longer introduces a difference in the modal gain. It is thus necessary to place erbium doping at locations that experience similar level of total signal intensity in order to achieve low DMG.

2.2 Second optimization criterion

The results presented in the previous section highlight the need for a second optimization criterion, Z2, that will take saturation of the population inversion into account to correctly estimate if a fiber design, with a specific RIP and EDP combination, will provide a low DMG. We are considering cladding-pumped FM-EDFAs and assume a uniform distribution of the pump power over the whole cladding area. We here further assume that the amplifier has uniform input signal conditions, which means that all signal modes are present with the same number of channel wavelengths and the same input power. We now define a new optimization criterion that will direct solutions towards designs for which the erbium doped regions will overlap with similar levels of total signal intensity, for example by placing erbium doping in the region where the total signal intensity, Ik, is flatter, and preferably equal (see Fig. 2). The criterion is defined as

$$ {Z_2} = \frac{{\mathop {\max }\limits_{k{\; }\epsilon {\; }EDP} ({{I_{({k,\rho \neq 0} )}}} )}}{{\mathop {\min }\limits_{k{\; }\epsilon {\; }EDP} ({{I_{({k,\rho \neq 0} )}}} )}},$$
where ${I_{({k,\; \; 0} )}}$ is the total signal intensity over a specific k region of the fiber where erbium doping is present.

To achieve the ideal DMG of 0 dB, a few mode EDF must thus meet the following three conditions: it has a Z1 = 1, i.e. all the modes have the same overlap with erbium ion population; it has Z2 = 1, i.e. the total signal intensity is uniform over the doped region such that saturation will be uniform for all modes; and the same signal power spectrum is injected for each guided mode in the fiber (although the absolute power level can change provided that the change is uniform across channels). Respecting these conditions, the signal power spectrum will be modified for each mode exactly the same way between zi and zj = zi+ Δz, starting with zi = 0, such that the three conditions remain true at any other longitudinal position along the fiber. Now, evaluating the Z2 criterion for EDF-A and EDF-B, we respectively obtain 1.582 and 1.053, which shows that the signal intensity inside the doped regions of EDF-B is more uniform than in EDF-A, thereby explaining the previously demonstrated DMG difference between these fibers.

Thus, considering that a fiber design that respects the conditions mentioned above could theoretically lead to a DMG of 0 dB, we propose a novel FOM and design method that targets these solutions. Optimization of fiber designs with this method involves parameter sweeping of combinations of RIP and EDP, but is much less computationally demanding than using the full simulation model with Runge-Kutta solutions for each possible combination of parameters.

2.3 Figure-of-merit

Our goal is to examine a wide variety of designs by sweeping parameters describing the RIP and EDP of the fiber. The criteria Z1 and Z2 for each of these designs are evaluated from overlaps of the EDP with the mode intensity profiles and total intensity profile. Then, we choose the optimal design between all the possible designs for which Z1 and Z2 were calculated by evaluating the figure of merit Ztot, which is defined as

$$ {Z_{tot}} = ({Z_1} - 1)({Z_2} - 1).$$
We examine below the usefulness of this FOM that minimizes both Z1 and Z2. In the following method, we calculated Z1 and Z2 individually to gain more physical insight on the designs. However, one could also simply use Ztot as the sole optimization criteria to obtain a low DMG FM-EDF, without first calculating the Z1 and Z2 values.

2.4 Design method

As an example, in order to target designs compatible with common fabrication processes based on MCVD with solution doping, we focus on designs that have two-doped ring regions in the core. Furthermore, we consider at most two different values of refractive index in the core, one near the center of the core and one around. Although challenging to fabricate, fiber designs with similar small features were successfully fabricated with a MCVD-based manufacturing process in [31]. Such a generic profile is presented in Fig. 5. and the overall design process is illustrated by the schematic in Fig. 6. Despite its simplicity, the choice of the two-ring EDP offers great flexibility and provides a large number of explorable designs.

 figure: Fig. 5.

Fig. 5. EDP (green) and RIP (black) with the optimization parameters used in the simulations.

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. Schematic of the methodology used to design low-DMG few-mode EDFs.

Download Full Size | PDF

The method consists of five steps. The first step consists of defining the RIP such that the fiber core supports only the desired number of modes. We wish to explore a space of RIPs defined by a set of four parameters, namely the inner core index and radius (nic and ric), and the outer core index and radius (noc and roc) as represented in Fig. 6. We first choose to limit our investigation to a fixed outer core radius, roc = 10 µm, as it is a reasonable value for fabrication purposes, power coupling and bending loss. The second step consists of calculating the mode profiles using COMSOL while sweeping values of nic and ric. For every parameter set (nic, ric), we need to determine the highest possible noc. This value is found by identifying the noc for which the effective index of the lowest order unwanted mode is at least 10−5 lower than the index of the cladding (neffnclad – 10−5); thus the mode is not guided. We fix noc to this value. The mode profiles calculated with the parameter values (noc, roc, nic and ric) are stored and can then be used to compare different EDPs. In the third and fourth steps, we calculate the values of the Z1 and Z2 criteria by sweeping the values of the two-doped ring parameters, which are ρ1, r1,in, r1,out, ρ2, r2,in and r2,out. Finally, in the fifth step, we determine the optimized design by evaluating Ztot.

In the next sections, we use of this method to find the optimal design for two different types of RIP. Firstly, in section 3, we consider a simple step-index fiber with only one refractive index value in the core (ric = 0 µm). Then, in section 4, we consider a raised index region at the center of the fiber core and examine fiber designs with staircase RIP. It is worth noting that the choice to restrict the design of the 6-mode FM-EDFA to a two-ring erbium doping structure is simply a design choice. The method could be applicable to the exploration of a larger parameter space that would encompass a wider variety of designs. The initial parameter values of the two-ring structure in Table 2 were determined by coarse parameter sweeping. A similar approach could be developed for the case of core pumping but it would require redefining the optimization criteria to consider the different overlap of the various pump modes with the doped regions.

Tables Icon

Table 2. EDP parameters sweeping ranges for the step-index RIP design method.

3. Application of the design method for a step-index fiber

As a first example, in this section, we design a two-ring erbium doped 6-mode FM-EDFA with a step-index RIP using the method previously described. The results show the benefits of using the FOM. We follow the steps outlined in Fig. 6.

3.1 Determination of noc

The goal of this step is to adjust the RIP such that only 6 modes are guided. Since the chosen design is a step-index fiber (ric = 0 µm), and considering that we have set roc = 10 µm, we only need to choose the noc value that will result in 6 guides modes (LP01, LP11a, LP11b, LP21a, LP21b and LP02), it is set to noc = 1.449.

3.2 Calculating the mode profiles

The mode profiles are calculated using COMSOL although, for this simple step-index RIP, the analytical expression of the mode profiles could also be used. The resulting mode profiles have been presented earlier in Fig. 2.

3.3 Determining the optimal design using the figure-of-merit Ztot

We now calculate Z1, Z2, and Ztot for a wide variety of EDPs with two doped rings located in the fiber core. Based on previous simulation results and limiting erbium doping concentrations to achievable values, the sweeping range for each parameter are the presented in Table 2.

To examine if the optimization criteria Z1 and Z2, and the FOM Ztot are good predicators for low DMG, we performed full simulations of the resulting fiber designs under the scenario presented in Table 1, i.e. we solved the radially resolved population rate equations and the signal/pump propagation equations through numerical simulations for each design. In Fig. 7, we present a color map of DMGmax as a function of Z1 and Z2, while the relationship between DMGmax and the FOM Ztot is shown in Fig. 8. Therefore, in Fig. 7, each circle represents a different EDP that was generated and evaluated in the optimization process. The DMGmax value on the colormap was capped to a maximum of 2.5 dB to better illustrate the change in the low DMGmax values (some EDP led to DMG up to 8.65 dB). All the 526 data obtained from the simulations are present in Fig. 7, but data with DMG values above 2.5 dB are assigned the same color in order to enhance the discrimination between the lower DMG values on the color scale. Figure 7. shows the importance of evaluating both criteria. To choose the optimal design between all the generated EDP, we now turn to Ztot as the FOM and we obtain the fiber labeled as EDF-B, with the RIP and EDP profiles shown in Fig. 2. As discussed earlier, this design achieved 21.7 dB of minimal gain and less than 0.46 dB of DMG over the C-band for a fiber length of 30.0 m. EDF-B is characterized by Ztot = 0.0012.

 figure: Fig. 7.

Fig. 7. Color map of the DMGmax as a function of both optimization criteria in the case of a step-index refractive index profile (EDF-A and EDF-B are respectively represented by a black triangle outline and a black star outline filled with their appropriate color from the DMG color scale).

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. DMGmax as a function of the FOM Ztot for a step-index refractive index profile. The inset zooms on the low DMGmax region (EDF-A is shown by a green triangle and EDF-B are shown by a green star).

Download Full Size | PDF

4. Results for a staircase index fiber

We now examine if more complex RIP can be used to further reduce the DMG using the proposed optimization technique. Since we are targeting simple designs compatible with standard fiber fabrication techniques, we add only one additional step to the refractive index profile, resulting in the staircase profile shown in Fig. 5. The sweeping ranges of the different optimization parameters of the RIPs and EDPs are presented in Table 3. An optimized fiber, EDF-C, with the RIP and EDP shown in Fig. 9. left, is obtained. The optimization led to a RIP with a ric of 1.75 µm and a nic of 1.453, and an EDP with the first doped ring having an erbium ion concentration of 1.2 × 1025 ions/m3 located at r = 3 µm, and a second doped ring with a higher erbium doping concentration of 2.0 × 1025 ions/m3 located at r = 8 µm. Figure 9. right displays the six mode profiles as well as the total intensity, which shows that the total intensity in both doped rings is very similar in this optimum design.

 figure: Fig. 9.

Fig. 9. The staircase RIP and the two-ring EDP (left), and the simulated mode profiles plotted with the EDP (right), of the optimized design EDF-C.

Download Full Size | PDF

Tables Icon

Table 3. EDP and RIP parameters sweeping ranges for the staircase RIP design method.

This fiber design is now simulated with the full EDFA model described previously under the scenario presented in Table 1, and the results are shown in Fig. 10. In this case, the optimization technique led to a combination of RIP and EDP that achieved 22.6 dB of minimal gain and reached a low DMGmax of 0.18 dB over the C-band. The optimal fiber length was of 29 m and the maximum noise figure was of 4.08 dB. This fiber reached Ztot = 0.0006. It should be noted that we also investigated designs with a depressed index region at the center of the fiber core and it did not lead to improved performance.

 figure: Fig. 10.

Fig. 10. Gain (left) and DMG (right) for the LP01, LP11a, LP11b, LP21a, LP21b and LP02 modes of the few mode EDF design EDF-C, which has a staircase RIP (P­p = 20 W, Ps = -9 dBm/mode, fiber length = 29 m).

Download Full Size | PDF

Figure 11 shows the spectral gain and DMG of fiber EDF-C when the fiber length is chosen for optimum gain flatness instead of the highest minimum gain. With a length of 27.3 m, we obtain a minimum gain of 21.95 dB and a DMGmax of 0.1820dB.

 figure: Fig. 11.

Fig. 11. Gain (left) and DMG (right) for the LP01, LP11a, LP11b, LP21a, LP21b and LP02 modes of the few mode EDF design EDF-C, which has a staircase RIP and where the length has been optimized for gain flatness (Pp = 20 W, Ps = -9 dBm/mode, fiber length = 27.3 m).

Download Full Size | PDF

5. Sensitivity to the operation point

We now examine the earlier claim that the DMGmax variations will be minimum when the operating point of the amplifier changes. This property follows from the similar level of total signal intensity in the doped regions. If we compare the FOM value Ztot for EDF-B and EDF-C (0.0012 vs 0.0006) we see that EDF-C will perform better than EDF-B. However, both designs are robust with regards to operating conditions because Z2 was minimized in both cases. In Fig. 12, we plot the variations of DMGmax over the C-band as a function of fiber length for both designs, which shows almost no DMGmax variations as the fiber length deviates from the optimal fiber lengths represented by the dash-dotted lines. Although the inversion varies along the fiber length, all signal modes interact with erbium doped regions that are similarly inverted, which preserved the low DMGmax property. By a similar argument, the resulting designs also maintain a low DMGmax for a wide range of pump power and signal power. The variations of DMGmax as a function of the pump power injected in the cladding is presented in Fig. 13. DMGmax slightly increases with pump power, but so does the minimum gain, with the result that the ratio DMGmax/Gainmin in linear scale is approximately constant at 0.007197 ± 0.000004 for EDF-B and 0.005915 ± 0.000004 for EDF-C, over a pump power range from 15 to 30 W. The DMGmax as a function of the modal signal input power is shown in Fig. 14. It is seen that the DMGmax decreases as the input signal power per mode increases but, as the minimum gain also decreases, the ratio DMGmax/Gainmin is 0.00672 ± 0.00006 for EDF-B and 0.00555 ± 0.00005 for EDF-C over an input signal mode power range from -25 dBm to -5 dBm. In Fig. 12, 13 and 14, the parameter value corresponding to the optimized design are indicated by the vertical dash-dotted lines, while other parameters are held constant to the values used in the optimization, i.e. 20 W of pump power, -9 dBm/mode of signal power and the optimum fiber length. Figure 13 and Fig. 14 both show the increased robustness of the EDF-C design to changes in the operation point of an amplifier.

 figure: Fig. 12.

Fig. 12. DMGmax (over the C-band) as a function of the fiber length for EDF-B (dashed blue) and EDF-C (solid red). The dash-dotted vertical lines show the optimal length (Ps = 20 W and Ps = -9 dBm/mode). The inset shows the variation in Gainmin.

Download Full Size | PDF

 figure: Fig. 13.

Fig. 13. DMGmax (over the C-band) as a function of the pump power for EDF-B (blue dashed) and EDF-C (red solid). The dash-dotted vertical line indicates the pump power used to optimize the design (Ps = -9 dBm/mode, fiber length = 30 m for EDF-B and fiber length = 29 m for EDF-C). The inset shows the variation in Gainmin.

Download Full Size | PDF

 figure: Fig. 14.

Fig. 14. DMGmax (over the C-band) as a function of the signal input power in each mode for EDF-B (blue dashed) and EDF-C (red solid). The dash-dotted line shows the input signal power in each mode used to optimize the two fiber designs (Pp = 20 W, fiber length = 30 m for EDF-B and fiber length = 29 m for EDF-C). The inset shows the variation in Gainmin.

Download Full Size | PDF

6. Conclusion

We proposed novel FOM and design method based on two optimization criteria that results in low DMG in cladding-pumped few-mode EDFAs. This method has low computational cost, which allows to explore a wide variety of possible erbium-doping profiles. The interest of this method was validated by finding the best designs with two-doped ring regions in six-mode EDFs. Even if a two-doped ring 6-mode FM-EDFA was used as an example in this paper, the method in itself is not restricted to this fiber design. The design of a first EDF, with a simple step-index refractive index profile, showed that the usual criteria based on equal overlap of the fiber mode profile with the erbium doping is not sufficient to identify the optimum designs. By using the proposed FOM, we can identify the fiber designs that will provide the lowest DMG. The best design of a two-ring erbium doped few-mode EDF was achieved with the staircase refractive index profile for which a DMGmax of 0.18 dB and a minimal gain of 22.6 dB were obtained. The novelty of this FOM is that it considers the local saturation of the erbium ion population by the total signal intensity of all modes, which impacts the optimal distribution of the erbium doping in the core. With this method, fiber design optimization can be performed without having to resort to full EDFA simulations. Lastly, the use of this novel FOM also has the advantage of identifying designs that will maintain a low DMG for a very wide range of signal power, pump power and fiber length values. We believe that this method greatly simplifies the design of erbium-doped fibers for cladding-pumped few mode EDFAs.

Funding

Natural Sciences and Engineering Research Council of Canada (CRDPJ 538379-18).

Disclosures

The authors declare no conflicts of interest.

Data availability

All the data presented in this paper are not publicly available but could potentially be shared upon request.

References

1. P. J. Winzer, D. T. Neilson, and A. R. Chraplyvy, “Fiber-optic transmission and networking: the previous 20 and the next 20 years [Invited],” Opt. Express 26(18), 24190–24239 (2018). [CrossRef]  

2. R. J. Essiambre, G. J. Foschini, P. J. Winzer, G. Kramer, and B. Goebel, “Capacity Limits of Optical Fiber Networks,” J. Lightwave Technol. 28(4), 662–701 (2010). [CrossRef]  

3. D. J. Richardson, J. M. Fini, and L. E. Nelson, “Space-division multiplexing in optical fibres,” Nat. Photonics 7(5), 354–362 (2013). [CrossRef]  

4. A. F. Herbster and M. A. Romero, “Few-mode erbium-doped fiber amplifier design method based on the signal-pump overlap integral,” Opt. Eng. 53(9), 096101 (2014). [CrossRef]  

5. P. Sillard, “Next-generation fibers for space-division-multiplexed transmissions,” J. Lightwave Technol. 33(5), 1092–1099 (2015). [CrossRef]  

6. T. Mizuno and Y. Miyamoto, “High-capacity dense space division multiplexing transmission,” Opt. Fiber Technol. 35, 108–117 (2017). [CrossRef]  

7. M. Salsi, C. Koebele, D. Sperti, P. Tran, P. Brindel, H. Mardoyan, S. Bigo, A. Boutin, F. Verluise, P. Sillard, M. Astruc, L. Provost, F. Cerou, and G. Charlet, “Transmission at 2×100Gb/s, over two modes of 40km-long prototype few-mode fiber, using LCOS-based mode multiplexer and demultiplexer,” Optics InfoBase Conference Papers 19(17), PDPB9 (2011). [CrossRef]  

8. Y. Jung, S. U. Alam, D. J. Richardson, S. Ramachandran, and K. S. Abedin, “Multicore and multimode optical amplifiers for space division multiplexing,Optical Fiber Telecommunications VII, pp. 301–333, 2019, doi: 10.1016/B978-0-12-816502-7.00008-7.

9. Y. Jung, S.-U. Alam, and D. J. Richardson, “Optical Amplifiers for Mode Division Multiplexing,” Handbook of Optical Fibers, pp. 849–873, 2019, doi: 10.1007/978-981-10-7087-7_49.

10. A. F. Herbster and M. A. Romero, “Design of a FM-EDFA with Gain Correction for Few-Mode WDM Optical Networks,” 2019 SBFoton International Optics and Photonics Conference (SBFoton IOPC), pp. 1–5, 2019, doi: 10.1109/SBFoton-IOPC.2019.8910227.

11. M. Wada, T. Sakamoto, S. Aozasa, T. Mori, T. Yamamoto, N. Hanzawa, and F. Yamamoto, “L-band 2-LP mode EDFA with low modal dependent gain,” 2015 Optical Fiber Communications Conference and Exhibition (OFC), 2015, doi: 10.1364/OFC.2015.Tu3C.3.

12. N. Bai, E. Ip, Y. Luo, G. D. Peng, T. Wang, and G. Li, “Experimental study on multimode fiber amplifier using modal reconfigurable pump,” Optics InfoBase Conference Papers 2(2), OW1D.3 (2012). [CrossRef]  

13. G. Lopez-Galmiche, Z. Sanjabi Eznaveh, J. E. Antonio-Lopez, A. M. Velazquez Benitez, J. Rodriguez Asomoza, J. J. Sanchez Mondragon, C. Gonnet, P. Sillard, G. Li, A. Schülzgen, C. M. Okonkwo, and R. Amezcua Correa, “Few-mode erbium-doped fiber amplifier with photonic lantern for pump spatial mode control,” Opt. Lett. 41(11), 2588–2591 (2016). [CrossRef]  

14. Y. Jung, E. L. Lim, Q. Kang, T. C. May-Smith, N. H. L. Wong, R. Standish, F. Poletti, J. K. Sahu, S. U. Alam, and D. J. Richardson, “Cladding pumped few-mode EDFA for mode division multiplexed transmission,” Opt. Express 22(23), 29008–29013 (2014). [CrossRef]  

15. Y. Jung, S. U. Alam, Z. Li, A. Dhar, D. Giles, I. P. Giles, J. K. Sahu, F. Poletti, L. Grüner-Nielsen, and D. J. Richardson, “First demonstration and detailed characterization of a multimode amplifier for space division multiplexed transmission systems,” Opt. Express 19(26), B952–B957 (2011). [CrossRef]  

16. M. Wada, T. Sakamoto, S. Aozasa, T. Mori, T. Yamamoto, and K. Nakajima, “Differential modal gain reduction of L-band 5-mode EDFA using EDF with center depressed core index,” J. Lightwave Technol. 35(4), 763–767 (2017). [CrossRef]  

17. M. Yamada, D. Nobuhira, S. Miyagawa, O. Koyama, and H. Ono, “Doping concentration distribution in 2-signal LP-mode ring-core erbium-doped fiber,” Appl. Opt. 56(36), 10040–10045 (2017). [CrossRef]  

18. A. Gaur, G. Kumar, and V. Rastogi, “Modal gain equalization of 18 modes using a single-trench ring-core EDFA,” J. Opt. Soc. Am. B 35(9), 2211–2216 (2018). [CrossRef]  

19. Y. Amma, T. Hosokawa, H. Ono, K. Ichii, K. Takenaga, S. Matsuo, and M. Yamada, “Ring-core multicore few-mode erbium-doped fiber amplifier,” IEEE Photonics Technol. Lett. 29(24), 2163–2166 (2017). [CrossRef]  

20. Q. Kang, E. Lim, Y. Jun, X. Jin, F.P. Payne, S. U. Alam, and D. J. Richardson, “Gain equalization of a six-mode-group ring core multimode EDFA,” 2014 European Conference and Exhibition on Optical Communication, (ECOC), 2014, doi: 10.1109/ECOC.2014.6964017.

21. Y. Jung, Q. Kang, H. Zhou, R. Zhang, S. Chen, H. Wang, Y. Yang, X. Jin, F. P. Payne, S. U. Alam, and D. J. Richardson, “Low-loss 25.3 km few-mode ring core fiber for mode-division multiplexed transmission,” J. Lightwave Technol. 35(8), 1363–1368 (2017). [CrossRef]  

22. Y. Jung, Q. Kang, L. Shen, S. Chen, H. Wang, Y. Yang, K. Shi, B. C. Thomsen, R. A. Correa, Z. S. Eznaveh, J. C. A. Zacarias, J. E. Antonio-Lopez, P. Barua, J. K. Sahu, S. U. Alam, and D. J. Richardson, “Few mode ring-core fibre amplifier for low differential modal gain,” 2017 European Conference and Exhibition on Optical Communication, (ECOC), 2017, doi: 10.1109/ECOC.2017.8345893.

23. Y. Wakayama, K. Igarashi, D. Soma, H. Taga, and T. Tsuritani, “Novel 6-Mode Fibre Amplifier with Large Erbium-Doped Area for Differential Modal Gain Minimization,” ECOC 2016; 42nd European Conference on Optical Communication, 2016.

24. N. K. Fontaine, B. Huang, Z. Sanjabieznaveh, H. Chen, C. Jin, B. Ercan, A. Velázquez-Benetez, S.H. Chang, R. Ryf, A. Schülzgen, J. Carlos Alvarado, P. Sillard, C. Gonnet, J. E. Antonio-Lopez, and R. Amezcua Correa, “Multi-mode Optical Fiber Amplifier Supporting over 10 Spatial Modes,” Optical Fiber Communication Conference Postdeadline Papers, pp. Th5A.4, 2016, doi: 10.1364/OFC.2016.Th5A.4.

25. Q. Kang, E. L. Lim, Y. Jung, J. K. Sahu, F. Poletti, C. Baskiotis, S. U. Alam, and D. J. Richardson, “Accurate modal gain control in a multimode erbium doped fiber amplifier incorporating ring doping and a simple LP01 pump configuration,” Opt. Express 20(19), 20835–20843 (2012). [CrossRef]  

26. Y. Zeng, Y. Fang, Y. Qin, O. Xu, J. Li, and S. Fu, “Rigorous FM-EDF design with an oversized two-layer erbium ion distribution for C-band DMG mitigation,” J. Opt. Soc. Am. B 38(12), F1–F7 (2021). [CrossRef]  

27. C. Jin, B. Ung, Y. Messaddeq, and S. LaRochelle, “Tailored modal gain in a multi-mode erbium-doped fiber amplifier based on engineered ring doping profiles,” Photonics North 8915, 89150A (2013). [CrossRef]  

28. Z. Sanjabi Eznaveh, H. Chen, J.E. Antonio Lopez, J.C. Alvarado Zacarias, B. Huang, A. Amezcua Correa, C. Gonnet, P. Sillard, G. Li, A. Schülzgen, R. Ryf, and R. Amezcua Correa, “Ultra-low DMG multimode EDFA,” 2017 Optical Fiber Communications Conference and Exhibition (OFC), pp. 1–3, 2017.

29. Q. Kang, E.-L. Lim, F. Poletti, Y. Jung, C. Baskiotis, S.-u. Alam, and D. J. Richardson, “Minimizing differential modal gain in cladding-pumped EDFAs supporting four and six mode groups,” Opt. Express 22(18), 21499–21507 (2014). [CrossRef]  

30. K. S. Abedin, T. F. Taunay, M. Fishteyn, D. J. DiGiovanni, V. R. Supradeepa, J. M. Fini, M. F. Yan, B. Zhu, E. M. Monberg, and F. V. Dimarcello, “Cladding-pumped erbium-doped multicore fiber amplifier,” Opt. Express 20(18), 20191–20200 (2012). [CrossRef]  

31. C. Matte-Breton, R.-J. Essiambre, C. Kelly, Y. Messaddeq, and S. LaRochelle, “Multicore Cladding-Pumped Fiber Amplifier With Annular Erbium Doping for Low Gain Compression,” J. Lightwave Technol. 40(6), 1836–1846 (2022). [CrossRef]  

Data availability

All the data presented in this paper are not publicly available but could potentially be shared upon request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (14)

Fig. 1.
Fig. 1. Fiber core divided into K rings (annular regions), each ring k characterized by its erbium ion concentration ${\rho _{Er,(k )}}$ and mode intensity distribution ${\mathrm{\Gamma }_{s,({\textrm{m},{\mathrm{\lambda }_i},k} )}}$.
Fig. 2.
Fig. 2. The mode profile distributions and the total signal intensity of EDF-A and EDF-B across the fiber core. The green solid and dashed lines are the two-ring erbium doping profile of EDF-A and EDF-B respectively. The inset shows the RIP used for both designs.
Fig. 3.
Fig. 3. Gain (left) and DMG (right) over the C-band in EDF-A for the six modes LP01, LP11a, LP11b, LP21a, LP21b and LP02 (P­p = 20 W, Ps = -9 dBm/mode, fiber length = 23 m).
Fig. 4.
Fig. 4. Gain (left) and DMG (right) over the C-band in EDF-B for the six modes LP01, LP11a, LP11b, LP21a, LP21b and LP02 (P­p = 20 W, Ps = -9 dBm/mode, fiber length = 30 m).
Fig. 5.
Fig. 5. EDP (green) and RIP (black) with the optimization parameters used in the simulations.
Fig. 6.
Fig. 6. Schematic of the methodology used to design low-DMG few-mode EDFs.
Fig. 7.
Fig. 7. Color map of the DMGmax as a function of both optimization criteria in the case of a step-index refractive index profile (EDF-A and EDF-B are respectively represented by a black triangle outline and a black star outline filled with their appropriate color from the DMG color scale).
Fig. 8.
Fig. 8. DMGmax as a function of the FOM Ztot for a step-index refractive index profile. The inset zooms on the low DMGmax region (EDF-A is shown by a green triangle and EDF-B are shown by a green star).
Fig. 9.
Fig. 9. The staircase RIP and the two-ring EDP (left), and the simulated mode profiles plotted with the EDP (right), of the optimized design EDF-C.
Fig. 10.
Fig. 10. Gain (left) and DMG (right) for the LP01, LP11a, LP11b, LP21a, LP21b and LP02 modes of the few mode EDF design EDF-C, which has a staircase RIP (P­p = 20 W, Ps = -9 dBm/mode, fiber length = 29 m).
Fig. 11.
Fig. 11. Gain (left) and DMG (right) for the LP01, LP11a, LP11b, LP21a, LP21b and LP02 modes of the few mode EDF design EDF-C, which has a staircase RIP and where the length has been optimized for gain flatness (Pp = 20 W, Ps = -9 dBm/mode, fiber length = 27.3 m).
Fig. 12.
Fig. 12. DMGmax (over the C-band) as a function of the fiber length for EDF-B (dashed blue) and EDF-C (solid red). The dash-dotted vertical lines show the optimal length (Ps = 20 W and Ps = -9 dBm/mode). The inset shows the variation in Gainmin.
Fig. 13.
Fig. 13. DMGmax (over the C-band) as a function of the pump power for EDF-B (blue dashed) and EDF-C (red solid). The dash-dotted vertical line indicates the pump power used to optimize the design (Ps = -9 dBm/mode, fiber length = 30 m for EDF-B and fiber length = 29 m for EDF-C). The inset shows the variation in Gainmin.
Fig. 14.
Fig. 14. DMGmax (over the C-band) as a function of the signal input power in each mode for EDF-B (blue dashed) and EDF-C (red solid). The dash-dotted line shows the input signal power in each mode used to optimize the two fiber designs (Pp = 20 W, fiber length = 30 m for EDF-B and fiber length = 29 m for EDF-C). The inset shows the variation in Gainmin.

Tables (3)

Tables Icon

Table 1. FM-EDFA simulation parameters for evaluation of gain, DMG and noise figure

Tables Icon

Table 2. EDP parameters sweeping ranges for the step-index RIP design method.

Tables Icon

Table 3. EDP and RIP parameters sweeping ranges for the staircase RIP design method.

Equations (4)

Equations on this page are rendered with MathJax. Learn more.

D M G ( λ ) = max m ϵ M ( G ( m i , λ ) G ( m j , λ ) ,
Z 1 = max λ i ( max m M ( k = 1 K Γ s , ( m , λ i , k ) ρ E r , ( k ) ) min m M ( k = 1 K Γ s , ( m , λ i , k ) ρ E r , ( k ) ) ) ,
Z 2 = max k ϵ E D P ( I ( k , ρ 0 ) ) min k ϵ E D P ( I ( k , ρ 0 ) ) ,
Z t o t = ( Z 1 1 ) ( Z 2 1 ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.