We present a coupled-mode model of transverse mode instability in high-power fiber amplifiers, which takes the effect of gain saturation into account. The model provides simple semi-analytical formulas for the mode instability threshold, which are valid also for highly saturated amplifiers. The model is compared to recently published detailed numerical simulations of mode instability, and we find reasonably good agreement with our simplified coupled-mode model.
© 2014 Optical Society of America
The recent discovery of transverse mode instability (TMI) in high-power fiber amplifiers  has motivated experimental [2–4] and theoretical [5–13] efforts to characterize and explain this phenomenon. TMI manifests itself as an instability in which the signal power fluctuates on a ms timescale between the fundamental mode (FM) and the first few higher-order modes (HOMs) when the average signal power is increased above a certain threshold. The experimental and theoretical research into this phenomenon has provided strong evidence that the underlying physical mechanism responsible for TMI is a form of stimulated thermal Rayleigh scattering (STRS) , caused by quantum defect heating of the active fiber. Essentially, the presence of light in the FM and HOM of a multimode active fiber results in a mode beating pattern in the signal intensity, which leads to a corresponding spatially and temporally oscillating thermal profile. Since the refractive index of silica is temperature dependent, a thermally induced long period grating (LPG) is generated, which has the correct period to couple the two modes. Due to the fact that this thermal nonlinearity is extremely slow, with a typical timescale of several ms, a phase-lag between this LPG and the mode beating pattern will exist if the latter changes due to noise or transients in the signal, which is necessary to induce strong coupling between the modes .
Several numerical models of TMI have been based on variants of the beam propagation method (BPM) [5–10]. While such numerical simulations can potentially include all relevant physical effects, the need to include the temporal dynamics, either by performing the simulation in the time domain or by including several frequency components in the frequency domain, leads to a numerically very intensive problem. Indeed, these numerical codes have so far employed various approximations, such as limiting the number of frequency components in the simulation [6, 7], expanding the electric field in a small number of transverse modes [8, 10] or limiting the simulation time in time domain algorithms .
In previous publications, we have proposed an alternative model of TMI in terms of a coupled-mode model [11, 12]. While this model has the advantage of numerical efficiency, and even allows for semi-analytical solutions, it does rely on a number of approximations, the most limiting of which is the neglect of spatiotemporal oscillations in the gain coefficient due to gain saturation. This approximation leads to an overestimate of the strength of the thermally induced LPG, which in turn leads to an underestimate of the TMI threshold power. In this paper, we show how this important effect can be included to lowest order in a coupled-mode formulation of TMI in fiber amplifiers.
2. Coupled-mode model
The change in temperature relative to the cooling fluid ΔT obeys the heat equation8], where it was found that TMI is not strongly dependent on the longitudinal heat diffusion. In fibers with a mode beat length comparable to the core radius, this approximation may break down and longitudinal heat diffusion would need to be taken into account. Assuming a steady-state solution to the rate equations of the gain medium, the heat load is given by Eq. (6) that the presence of light in the HOM leads to an oscillating heat load Q, which results in a thermally induced LPG. However, it is also clear from Eq. (7) that the strength of this grating is reduced if the gain is highly saturated, since the maxima in the mode beating intensity pattern coincide with minima in the gain coefficient under the assumption of fixed g0 and Isat. In the extreme case where I >> Isat, we see that the oscillations in the heat load disappear, and no LPG is induced.
Assuming that the fiber amplifier is operating below the TMI threshold, where the excitation of the HOM is small compared to the FM, we can expand Eq. (6) in powers of Ĩ/(I0 + Isat) which yields10] to show that TMI does not occur under steady-state conditions, as also pointed out in .
As in  we can now derive coupled-mode equations, which are valid also for non-steady-state conditions, for the mode amplitudes pi by solving Eq. (4) in the frequency domain using an appropriate Green’s function12]. The resulting coupled-mode equations in terms of mode power Pi = |pi|2 are Eqns. (15), in which spatial hole burning is only approximately taken into account. It is given by
The nonlinear gain coefficient χ is given byEq. (19) that the nonlinear gain coefficient χ is reduced as the signal power becomes comparable to the saturation power. This is due to the above mentioned reduction in the strength of the thermally induced LPG. While it is possible to calculate the integral in Eq. (19) for each z in a numerical solution of the coupled-mode equations, we shall show that a more efficient approach is to approximate the inner integral in Eq. (19) so that the integral only need to be calculated once. We assume that the mode function of the FM ψ1 in the denominator of the integrand can be approximated to be constant over the doped region, such that Eq. (19). We hereby obtain the approximate expression 11, 12] Eqns. (22) are similar to the ones derived in , except for the fact that the nonlinear gain coefficient is now z-dependent and has a reduced value due to the effect of gain saturation.
In order to check the approximation used in obtaining the approximate expression for A given in Eq. (21), we compare the value of χ calculated using the exact expression in Eq. (19) and the approximate expression given in Eq. (21) for a typical double-clad SIF and a range of values of P1/Psat. Two SIFs, called Fiber A and Fiber B, with different radii of the doped region are considered to investigate the effect of a reduced doped region. The relevant fiber parameters are given in Table 1 and the result is shown in Fig. 1. It is clear that the nonlinear gain coefficient decreases significantly as the signal power becomes comparable to the saturation power, which is due to the additional weakening of the thermally induced LPG caused by gain saturation. The approximate expression for χ, calculated using Eq. (21), is seen to provide a reasonable approximation to the exact value. The approximation is seen to be better for a fiber with a reduced doped region of the core, which can be understood by considering the integral over the doped region in Eq. (19). For a smaller doped region, the approximation of a nearly constant FM profile over the doped area improves.
3. Solution of the coupled-mode equations
In order to solve the coupled-mode equations given by Eqns. (22) we make the approximation that the light in the FM is undepleted by TMI. Since the expansion in Eq. (13) is only valid for a weak excitation of the HOM, our model applies only to operation below or near the TMI threshold, and hence the approximation of an undepleted FM should be valid. In addition, we introduce a propagation equation for the pump light. We thus consider the following system of 3 coupled ODEs12] that a likely source for seeding the TMI is intensity noise in the input signal. In this case the input power spectral density is Eqns. (26a) and (26b) numerically using a standard numerical ODE solver followed by numerical quadrature. Likewise, the integral in Eq. (25) required for calculating χ0 is evaluated for a range of frequencies using standard numerical quadrature techniques.
We calculated the HOM content ξ(L) as a function of output signal power for the fiber given in Table 1. The result is shown in Fig. 2 for both co- and counter-pumped operation, and is also compared to the result obtained when the effect of the gain oscillations are neglected as in [11,12]. The input signal power was taken to be 10 W and the input HOM content ξ(0) = 0.01, which should correspond to nearly optimal coupling of the signal into the fiber core. The relative intensity noise was assumed to be −100 dBc/Hz, which we believe is representative of realistic experimental conditions. Defining the TMI threshold power Pth as the output signal power at which ξ(L) = 0.1, we find Pth = 723 W for counter-pumped operation and Pth = 672 W for co-pumped operation. Compared to the corresponding value of Pth = 260 W obtained when gain oscillations are neglected, it is clear that this effect is very important in this case. Interestingly, there is little difference in TMI threshold between the co- and counter-pumped cases. When gain saturation is low and the effect of the associated gain oscillations can be ignored, this is expected, as shown in . Apparently, including the effects of gain saturation in the model approximately preserves the symmetry between co- and counter-pumping.
We have also performed the same calculation for Fiber B. The result is shown in Fig. 3. The TMI thresholds for co- and counter-pumped operation are 1047 W and 1286 W, respectively, which is substantially higher than for Fiber A. This is due to the combination of the reduced overlap between the modes with the doped region and the longer fiber length. From the dashed curve in Fig. 3 we see that the TMI threshold, calculated under the assumption of negligible gain oscillations, is 331 W for Fiber B. This increase of 27% compared to Fiber A is due to the reduced mode overlap with the doped region. However, the increase in TMI threshold when taking gain oscillations into account is much larger, 78% for counter-pumped and 56% for co-pumped. This increase in TMI threshold thus cannot solely be explained by the reduced mode overlap, but is also caused by the higher gain saturation. We have verified this conclusion by simulating a shorter version of Fiber B, which has the same length as Fiber A. In this case, the TMI threshold for co-pumped operation drops to 767 W compared to 1047 W for the original Fiber B. This is due to the fact that the gain is much less saturated in this case, which in turn implies that the gain oscillations are weaker and the thermally induced LPG is stronger. However, this shorter version of Fiber B has a very poor efficiency due to the lower pump absorption, and is thus not considered a realistic case.
4. Comparison to BPM calculations by Smith and Smith
In a recent paper , Smith and Smith investigated the effect of gain saturation on TMI numerically using a BPM code. In this section we compare our calculation of the nonlinear gain coefficient χ with a calculation of the same quantity obtained by the BPM presented in Fig. 5 in . The fiber parameters for the two SIFs used for the comparison are given in Table 2. The first fiber, here called Fiber C, is a short fiber with a small pump cladding, which provides high pump intensity and hence lower gain saturation, while the second fiber, called Fiber D, is a longer fiber with a much larger pump cladding, which therefore operates with much higher gain saturation. We have calculated the nonlinear gain coefficient χ using the approximate expression given in Eq. (21) as a function of z for a fixed frequency offset Δν = 1.1 kHz and the specific pump powers given in . The result for Fiber C is shown in Fig. 4, while the result for Fiber D is shown in Fig. 5. Compared to Fig. 5 in  we find a reasonably good agreement with our results, but with some interesting differences. For Fiber C, we find an initially lower value of χ for the counter-pumped case, which rises to a slightly lower value than in . For the co-pumped case, χ falls off more quickly with z compared to  and also reaches a lower value. The same is true for the co-pumped case for Fiber D, while the counter-pumped case shows an initial value of χ slightly higher than in , which then falls off with z to a slightly lower value. Interestingly, the counter-pumped case in the lower plot of Fig. 5 in  shows a slightly increasing nonlinear gain χ with z, which according to our analysis would seem to indicate a decreasing gain saturation with z. However from Fig. 2 in , it appears that the gain saturation is in fact increasing with z. A possible explanation of this apparent discrepancy could be thermal lensing, since the larger heat load toward the output end of the fiber will cause the overlap between the HOM and FM to increase, which will increase the nonlinear gain.
To test the accuracy of the approximate calculation of χ for Fiber C and D, we again compare the exact and approximate expressions for a fixed frequency offset Δν = 1.1 kHz and for a range of values of P1/Psat. The result is shown in Fig. 6 and shows good agreement between the exact and approximate calculations even when gain saturation is substantial. The discrepancy between the calculation in  and our model is thus unlikely to be due to this approximation, but is rather due to other approximations, such as the fact that we do not include thermal lensing, and the fact that we treat the gain oscillations to first order only.
In summary we have shown how the effects of gain saturation can be included approximately in a coupled-mode model of TMI. The resulting model is numerically very efficient and easy to implement, and can be used to obtain quantitative estimates of the TMI threshold for a given fiber amplifier under conditions where gain saturation is substantial. Applying our model to particular SIF designs (Fiber A & B), we found that the estimate for the TMI threshold increased by approximately a factor of 3 compared to a model which does not include the effect of gain saturation. Only in the case of a short, rod-type fiber (Fiber C) did we find a nonlinear gain comparable to the prediction of the model when neglecting gain saturation effects. We also found that the benefit of reducing the doped area of the core is significantly greater than predicted by the simpler model, due to the increased gain saturation. Including gain saturation in the model breaks the symmetry between co- and counter-pumping, with the latter having the higher TMI threshold.
Comparing our model to BPM simulations presented in , we found reasonably good agreement between the nonlinear gain coefficient calculated with our model and the BPM results, although our model results in slightly lower values for the nonlinear gain coefficient. This is likely due to the various approximations employed by our model, most importantly the neglect of thermal lensing and the fact that gain saturation is only included to lowest order of approximation. We thus expect our model to be able to provide good quantitative estimates of the TMI thresholds of fiber amplifier designs under conditions where gain saturation is significant.
References and links
1. T. Eidam, C. Wirth, C. jauregui, F. Stutzki, F. Jansen, H.-J. Otto, O. Schmidt, T. Schreiber, J. Limpert, and A. Tünnermann, “Experimental observations of the threshold-like onset of mode instabilities in high power fiber amplifiers,” Opt. Express 19, 13218–13224 (2011). [CrossRef] [PubMed]
2. F. Stutzki, H.-J. Otto, F. Jansen, C. Gaida, C. Jauregui, J. Limpert, and A. Tünnermann, “High-speed modal decomposition of mode instabilities in high-power fiber lasers,” Opt. Lett. 36, 4572–4574 (2011). [CrossRef] [PubMed]
3. H.-J. Otto, F. Stutzki, F. Jansen, T. Eidam, C. Jauregui, J. Limpert, and A. Tünnermann, “Temporal dynamics of mode instabilities in high-power fiber lasers and amplifiers,” Opt. Express 20, 15710–15722 (2012). [CrossRef] [PubMed]
4. M. M. Johansen, M. Laurila, M. D. Maack, D. Noordegraaf, C. Jakobsen, T. T. Alkeskjold, and J. Lægsgaard, “Frequency resolved transverse mode instability in rod fiber amplifiers,” Opt. Express 21, 21847–21856 (2013). [CrossRef] [PubMed]
10. S. Naderi, I. Dajani, T. Madden, and C. Robin, “Investigations of modal instabilities in fiber amplifiers through detailed numerical simulations,” Opt. Express 21, 16111–16129 (2013). [CrossRef] [PubMed]