We show by numerical modeling that saturation of the population inversion reduces the stimulated thermal Rayleigh gain relative to the laser gain in large mode area fiber amplifiers. We show how to exploit this effect to raise mode instability thresholds by a substantial factor. We also demonstrate that when suppression of stimulated Brillouin scattering and the population saturation effect are both taken into account, counter-pumped amplifiers have higher mode instability thresholds than co-pumped amplifiers for fully Yb3+ doped cores, and confined doping can further raise the thresholds.
© 2013 OSA
In earlier papers we described a stimulated thermal Rayleigh scattering process (STRS) that can account for observed modal instability in large mode area fiber amplifiers, and we described in detail our numerical model of this process [1–5]. Our model computes the mode coupling gain and the laser gain for fiber amplifiers with practically achievable step index profiles operating at realistic pump and signal powers. It was shown that the quantum defect heating associated with laser amplification, in conjunction with a frequency offset between the strong fundamental mode and the weak parasitic higher order mode, leads to a high exponential gain for the parasitic mode. Consequently, mode instability has a sharp power threshold, above which the output beam quality is severely degraded.
The essence of the STRS process responsible for mode instability is that laser gain necessarily deposits quantum defect heat in the core of the amplifier fiber, and the asymmetric heating produced by the asymmetric signal irradiance profile due to interference between modes LP01 and LP11 leads to an asymmetric thermal lens that couples light between those two modes. An additional requirement is that there be a phase shift between the temperature grating and the signal irradiance grating which, in our STRS model, is created by the time lag between an irradiance grating traveling along the fiber and the temperature grating that it creates.
Our model imposes a steady-periodic condition on the gain grating because we assume the frequency offset between modes has a narrow linewidth. Alternative STRS models have been presented by Hansen et al.[6, 7] and by Dong . They used similar approximations, including the steady-periodic assumption. Another model by Ward et al. is based on STRS, but without the steady-periodic assumption. Jauregui et al. also model temperature gratings in fiber amplifiers.
The primary difference in the physics of the various models is that the models of Hansen et al. and Dong assume the profile of the quantum defect heating matches the profile of the signal light, while our model and those of Ward et al. and Jauregui et al. compute the heat profile based on the local population inversion plus either the local pump or signal irradiance. In the model of Ward et al. the increase of signal power in a mode due to laser gain is found from the overlap of the local gain g(x,y,z) with the field of that mode. In our BPM model laser gain increases the total signal field locally and is then apportioned among the modes automatically by diffraction in the presence of the core index step. Either the local signal field growth or the local pump irradiance loss is used to compute the quantum defect heating.
Both our model and that of Hansen et al. use the same steady-periodic assumption. Using this treatment they can both include multiple frequencies for the beat note between the coupled modes to address bandwidth issues. However, we find that in computing thresholds using a single frequency at the gain peak produces nearly identical thresholds to using multiple frequencies contained in the gain band. Since added frequencies slow the computation we use a single frequency in all the computations presented in this paper. The seed power is 10−16 W to approximate quantum noise. Only step index fibers will be discussed.
Hansen et al. and Dong showed that using their assumed heat profile the mode coupling gain is related in a simple way to the laser gain. This would imply that the mode instability threshold is largely determined by the net laser gain, and can be adjusted only by changing the modal profiles and their overlap with the Yb3+ doping profile. In fact, they show that by using a low value of the V parameter or by confining the Yb3+ doping to the central portion of the core, the threshold is raised. However, we show in this report that the existence of depletion of the upper state Yb3+ population breaks the simple connection between laser gain and mode coupling gain, making it possible to design fibers with substantially higher mode coupling thresholds than predicted by Hansen et al. and by Dong. An amplifier that is designed to reach a specific level of laser amplification can be designed to avoid the mode instability. Specifically, by designing the amplifier so it has a higher ratio of pump cladding diameter to core diameter, the STRS threshold can be raised. Of course this requires a longer fiber and a smaller core, a combination that is problematic for SBS suppression. In practice the amplifier must achieve a suitable balance between STRS and SBS suppression. We will also discuss how to achieve this balance.
2. Transverse hole burning
In computing the upper state population fraction nu we use the steady state expressionTable 1. We model a fiber with a 50 μm diameter step index core with a typical numerical aperture of 0.054. The Yb3+ doping density is also typical of high power amplifiers.
Figures 1 and 2 show upper state population profiles, nu(x, y) computed at three points along co-pumped and counter-pumped fiber amplifiers with a cladding diameter of 400 μm. The signal light is all in LP01 here. Near the input end of the co-pumped fiber the pump is strong and the signal is weak so the upper state population is weakly depleted. The shape of the population depletion closely matches the signal irradiance profile here. The undepleted population in the region with Is = 0 is . Further along this fiber, near the crossing of the signal and pump powers, the saturation is strong at the mode center but moderate near the core boundary, while near the output end saturation is strong across the entire core. If the same fiber is counter-pumped the signal and pump powers are approximately equal along the full length of the fiber, so the saturation resembles that of the middle diagram of the co-pumped fiber where the signal and pump powers are equal. The degree of saturation would be reduced in a fiber with a smaller cladding, and increased in one with a larger cladding.
The quantum defect heating is proportional to the pump absorption, which is determined from the value of nu(x, y) and the pump irradiance. The latter is assumed uniform across the pump cladding. We compute the heat deposition rate Q usingFigs. 1 and 2 and multiplying them by the doping profile NYb(x, y).
The portion of the heat profile that is responsible for mode coupling gain (STRS gain) is the antisymmetric part created by the antisymmetric part of the signal irradiance. The population saturation strongly influences the shape of this part of the heating. In Figs. 3 and 4 we show cuts through the oscillatory heat profiles at the same three z locations as in Figs. 1 and 2. These heat profiles create the oscillating asymmetric part of the temperature profile that couples modes LP01 and LP11 in the STRS process. From an examination of these heat profiles it will come as no surprise when we show in the next section that strong saturation causes a significant reduction in STRS gain.
3. Mode coupling gain
3.1. Without hole burning
Both Dong and Hansen et al. showed that if the heat profile matches the signal irradiance profile over the portion of the core that is (uniformly) doped, the gain of mode LP11 satisfies
Under the assumption that the shape of the heat profile matches the irradiance profile, χ is independent of z and the modal powers. It depends on the spatial overlap of the two modes with one another and with the dopant profile. Hansen et al. and Dong both showed that confined doping reduces the value of χ, as does reducing the V parameter below 5 or so.
We used our model to verify the value of χ for one case analyzed by Hansen et al. (Fig. 4 of Ref. ). In order to avoid population saturation we used a pump power of 20 kW, dcore = 40 μm, dclad = 250 μm, and an LP01 seed power of 50 W. Otherwise the fiber parameters were those from Table 1. Our frequency of maximum gain was equal to that of Hansen et al., and our computed a value for χ was within 5% of that of Hansen et al. Of course, our fiber was extremely inefficient, using 20 kW to produce approximately 400 W of signal, and depleting the pump by only 2.5%. However, this exercise does demonstrate good agreement between our model and those of Hansen et al. and Dong in the limit of low saturation.
3.2. With hole burning
The model comparison just described illustrates the lack of realism of the assumption of equal shapes for the heat and irradiance profiles. For the heat profile to match the light profile the pump irradiance must be much stronger than the signal irradiance. However, amplifiers that efficiently convert pump power necessarily experience strong population saturation. To emphasize this contrast between our model and one that assumes matching heat and irradiance profiles and thus constant χ values, we define a new shape factor χ′ asEq. (3), gs is the computed signal laser gain, and gstrs is the STRS gain. We do not distinguish between the nearly equal laser gains for the two modes in this high V fiber with fully doped core. Variables Ps, gs, and gstrs are all z dependent computed values. The denominator is proportional to the total deposited heat. Figure 5 shows plots of χ′ for co-pumped and counter-pumped amplifiers.
The upper plot is for a 100 μm diameter cladding; the lower plot is for a 400 μm diameter cladding. Saturation effects are stronger with the larger cladding because of the reduced pump irradiances. In both fibers the value of χ′ near z = 0 for the co-pumped case is nearly equal to the χ we compute in the limit of zero saturation, as expected because of the low degree population saturation there. As z increases the population saturation strengthens so the value of χ′ falls, implying the mode coupling gain is reduced relative to the laser gain. This gain reduction raises the STRS threshold relative to that of Hansen et al. and Dong. In the small cladding fiber the reduction does not occur until half way along the fiber while in the large cladding fiber it occurs much sooner.
As expected from the nearly constant ratio of pump to signal for the counter-pumped fibers, the value of χ′ is nearly constant along the fiber. However, its value is still reduced to a level seen in the co-pumped fiber near the crossing point near z = 0.5 m for the small clad fiber and near z = 1.4 m for the large clad fiber. It is clear from a comparison of the values of χ′ that the STRS threshold should be higher for larger cladding sizes, and the threshold for an efficient amplifier should exceed the Hansen et al. and Dong threshold.
Another way to view the same information is to plot gcomp, gs, and the total heat rather than χ′. This is done in Fig. 6 for the same pair of fibers, co-pumped, with dclad = 100 μm (upper plot) and 400 μm (lower plot).
The expected trend of higher thresholds for larger cladding sizes is illustrated in Fig. 7. The computed thresholds, defined as 1% of the signal power in LP11, are seen to rise with increasing cladding diameters and the resulting increasing degree of population saturation. Co-and counter-pumped fibers are found to have similar thresholds in this study, a finding that appears to be quite general based on our experience modeling many different designs. In the fiber design used in Fig. 7, saturation appears to raise thresholds more significantly for confined doping than for full doping, probably because saturation is strong across the entire doped zone. For comparison, we also show in the figure the threshold predicted without saturation. More details for the model runs included in Fig. 7 are listed in Tables 2 and 5.
In the following tables Δν is the frequency of LP01 minus the frequency of LP11.
5. Stimulated Brillouin scattering (SBS) suppression
Designing a narrow bandwidth fiber amplifier suitable for beam combining applications requires balancing SBS and STRS. Here we present a simplified version of this balancing act. We define an effective SBS gain length at the STRS threshold by11] is exceeded if 11]. The γ parameter is an SBS gain reduction factor that reduces the effective SBS gain from its nominal value of gB = 5 × 10−11 m/W for silica. The linewidth for silica is approximately 50 MHz and the Stokes shift is approximately 16 GHz for λ = 1032 nm. The value of γ can be increased by broadening the signal linewidth by phase modulation  above 50 MHz, for example, or by introducing a temperature gradient  or a strain gradient [13, 14] along the fiber length to vary the Stokes shift. Another approach is to vary the acoustic velocity, and thus the Brillouin shift, across the core region [15, 16]. Temperature and strain shift the SBS frequency by known amounts, so from the signal linewidth and the temperature and strain gradient one can estimate γ using an SBS model. To avoid SBS the value of γ must satisfy
In Fig. 8 we plot the STRS threshold powers versus the quantity PthresLeff for our example fiber with dcore = 50 μm. Selecting an arbitrary value for PthresLeff on the horizontal axis (or equivalently the value of γ), the values of the STRS threshold can be read from the computed curves. From the four curves it appears that for any value of γ the counter-pumped fully doped amplifier provides the highest STRS threshold power, the co-pumped fully doped amplifier offers the lowest STRS threshold power, while the confined doping, co-pumped amplifiers give intermediate threshold powers.
6. Scaling for other core/cladding sizes
To first order changing the core diameter while keeping the ratio dcore/dclad fixed, and keeping NA = 0.054, the threshold powers do not change. The frequency offset changes proportional to 1/Aeff but the degree of saturation, and thus the value of χ′ is unaltered.
However, for dcore < 30 μm, the value of V is reduced to less than 5, so the mode overlap with the core becomes noticeably weaker, which tends to raise the STRS threshold [6, 8]. At the other extreme, for dcore > 80 μm, thermal lensing becomes significant, and this constriction of the modes may alter the STRS threshold.
As dcore changes the area Aeff changes and this changes the value of γ required to suppress SBS for a given value of PthresLeff according to Eq. (8). Similar powers are possible, while avoiding both SBS and STRS, but the value of γ necessary to avoid SBS scales approximately as . If there is a constraint on γ imposed by beam combining requirements, for example, this implies a lower limit on the core size.
7. Scaling for other LP11 starting powers
The thresholds computed in this study are based on an input signal power of 10−16 W in LP11, with a frequency shift to the STRS gain maximum, or threshold minimum. This input power level corresponds approximately to the quantum noise limit, and serves as a standard point of comparison among models. Actual starting powers will be considerably higher than this if there is amplitude modulation on the pump or signal seed with a modulation frequency near the optimum STRS frequency shift [3, 4]. Thermal noise is probably 2–3 factors of 10 higher than this as well . More realistic thresholds would be based on input powers of perhaps 10−10 W, but the exact level will require measurements of the modulation properties of the pump and seed. The reference thresholds presented here (Pref) can be used to estimate the threshold for a higher input power (Pstart) using
The thresholds computed in this study are also based on negligible signal loss due to absorbing impurities in the glass or due to photodarkening. The presence of these processes adds heat with a profile that contributes to STRS gain, and so reduces the thresholds .
We showed that transverse hole burning, or saturation of the population inversion, in a Yb3+ doped fiber can strongly influence the STRS or mode coupling gain, and this effect can be exploited to substantially raise mode instability thresholds. The degree of saturation can be manipulated by changing the ratio of pump irradiance to signal irradiance through choice of core and cladding sizes, injected seed power, and doping profile. Generally we find that otherwise identical co- and counter-pumped configurations have nearly equal thresholds even though the z dependence of saturation is quite different in the two.
The benefits of high saturation persist when SBS suppression is considered, even though high saturation implies longer fibers with higher PthresLeff values. The benefits of confined doping also survive the requirement of SBS suppression.
Fibers with high saturation have the added advantage of minimizing photodarkening, assuming photodarkening increases with higher upper state population density, as experiments seem to indicate .
We have shown that saturation has a strong effect on thresholds for step index fibers. Other, more complex refractive index profiles might be designed to raise the instability threshold by changing the shapes of the interacting modes. However, saturation can be expected to have a strong influence in those as well.
References and links
4. A.V. Smith and J.J. Smith, “Frequency dependence of mode coupling gain in Yb doped fiber amplifiers due to stimulated thermal Rayleigh scattering,” arXiv:1301.4277 [physics.optics](January 18, 2013).
5. A.V. Smith and J.J. Smith, “Maximizing the mode instability threshold of a fiber amplifier,” arXiv:1301.3489 [physics.optics](January 16, 2013).
10. C. Jauregui, T. Eidam, H.-J. Otto, F. Stutzki, F. Jansen, J. Limpert, and A. Tünnermann, “Temperature-induced index gratings and their impact on mode instabilities in high-power fiber laser systems,” Opt. Express 20, 440–451 (2012) [CrossRef] [PubMed] .
11. G.P. Agrawal, Nonlinear Fiber Optics(Academic, 1995).
12. M. Hildebrandt, S. Büsche, P. Wessels, M. Frede, and D. Kracht, “Brillouin scattering spectra in high-power single-frequency ytterbium doped fiber amplifiers,” Opt. Express 16, 15970–15979 (2008) [CrossRef] [PubMed] .
13. T. Horiguchi, T. Kurashima, and M. Tateda, “Tensile strain dependence of Brillouin frequency shift in silica optical fiber,” IEEE Phot. Tech. Lett. 1, 107–108 (1989) [CrossRef] .
14. J.E. Rothenberg, “Suppression of stimulated Brillouin scattering in single-frequency multi-kilowatt fiber amplifiers,” Proc. SPIE 6873,68730O (2008) [CrossRef] .
15. M.-J. Li, X. Chen, J. Wang, S. Gray, A. Liu, J.A. Demeritt, A.B. Ruffin, A.M. Crowley, D.T. Walton, and L.A. Zenteno, “Al/Ge co-doped large mode area fiber with high SBS threshold,” Opt. Express 15, 8290–8299 (2007) [CrossRef] [PubMed] .
17. A.V. Smith and J.J. Smith, “Spontaneous Rayleigh seeding of stimulated Rayleigh scattering in high power fiber amplifiers,” Opt. Express. (submitted).
18. M.J. Söderlund, J.J. Montiel i Ponsoda, S.K.T. Tammela, K. Ylä-Jarkko, A. Salokatve, and S. Honkanen, “Mode-induced transverse photodarkening loss variations in large-mode-area ytterbium doped silica fibers,” Opt. Express 16, 10633–10640 (2008) [CrossRef] [PubMed] .