Theoretical analysis of the mode coupling induced by heat of large-pitch micro-structured fibers*
Zhang Hai-Taoa)†, Chen Danb), Hao Jiea), Yan Pinga), Gong Ma-Lia)‡
Center for Photonics and Electronics, State Key Laboratory of Tribology, Department of Precision Instruments, Tsinghua University, Beijing 100084, China
Southwest Institute of Technical Physics, Chengdu 610041, China

Corresponding author. E-mail: zhanghaitao@mail.tsinghua.edu.cn

Corresponding author. E-mail: gongml@mail.tsinghua.edu.cn

Project supported by the National Natural Science Foundation of China (Grant No. 61475081)

Abstract

In this paper, a theoretical model to analyze the mode coupling induced by heat, when the fiber amplifier works at high power configuration, is proposed. The model mainly takes into consideration the mode field change due to the thermally induced refractive index change and the coupling between modes. A method to predict the largest average output power of fiber is also proposed according to the mode coupling theory. The largest average output power of a large pitch fiber with a core diameter of 190 μm and an available pulse energy of 100 mJ is predicted to be 540 W, which is the highest in large mode field fibers.

Keyword: 42.55.–f; 42.55.Wd; fiber amplifier; large pitch fiber; mode coupling
1. Introduction

Increasing the fiber core diameter is a prime method to reduce the nonlinear effect of a pulsed fiber amplifier. By employing micro-structured fiber (MSF) with non-continuous interface to increase the difference in loss between modes, the core diameter can be enlarged further to 100  μ m. It is demonstrated that near diffraction limited beam quality laser output with a pulse width of 20 ns and pulse energy of 22 mJ is achieved by employing a large pitch fiber (LPF) with 135-μ m core diameter.[1] It seems that the LPF is one of the most promising approaches for core diameters beyond 100  μ m.[24] However, MSFs, including large pitch fibers, suffer from mode instabilities when laser powers exceed a certain value due to losses of confinement compared with that in continuous interface, resulting in a dramatic decrease in beam quality.[58]

According to the 135-μ m MSF fiber structure, we for the first time fully investigate the output laser beam quality change trend of the discontinuous interface of MSF with core diameter varying from 51  μ m to 190  μ m in laser amplification. The mode competition caused by the mode gain and loss is considered. What is more, we focus on the influence of average output power on beam quality of the nearly 200-μ m fiber. With the increase of the fiber core size, the mode loss difference decreases. Taking the beam quality factor M2 < 1.5 as the upper limit value, the core diameter of MSF can be expanded to nearly 200  μ m. This is the biggest fiber core for the high beam quality laser output. Its biggest advantage is to effectively reduce the nonlinear effects in the pulse amplification, increasing the available pulse energy, which can reach 80  mJ/m. When M2 < 2, the average power can reach 540  W. This theoretical laser output power level has broken all of the available work intervals of the pulse fiber laser at mJ level (the pulse energy of more than 20  mJ has an average single pulse power lower than 150  W[2]).

The calculated results in this paper provide the reference for the laser available work interval of currently very large mode area MSFs. For further lowering the pulse nonlinear effect and obtaining MSF with a core diameter of more than 200  μ m, and also with a high quality beam output laser average power above 500  W, it must make a bigger breakthrough on the structure. But for pulsed fiber amplifiers requiring high pulse energy, the available pulse energy is as high as 100  mJ. So the situation where high energy and high average power cannot be realized at the same time should be broken.[9]

A super large mode area and large pitch fiber suitable for fiber drawing with diameter as large as 190-μ m has been designed, which is, to our knowledge, the largest core-diameter fiber able to maintain near diffraction limit beam quality. In this paper, in order to analyze the average power output ability of the 190-μ m Yb-doped large pitch fiber amplifier, we take into consideration the thermally induced mode coupling in mode competition. We arrive at the conclusion that the largest average output power of this amplifier can reach 540 W and its available pulse energy can reach 100  mJ when working in the near diffraction-limited configuration (M2 < 2.0).

2. Theoretical model

The large pitch fiber used in this paper and the characteristic mode distribution are shown in Fig.  1. The fiber has 18 air-holes that are hexagonally arranged in two circles with the central hole missing: the core diameter (dcore) is 190  μ m, the diameter of air hole cladding (d) is 19  μ m, and the air-hole pitch (Γ ) is 104.5  μ m. The fundamental mode (FM) has a mode diameter of 180  μ m, and a propagation loss of 1 dB/m. The first higher order mode (HOM) has a mode diameter of 273  μ m and a propagation loss of 19  dB/m, which is 18  dB/m higher than that of the FM. As to this fiber, the available pulse energy, which is defined as Eavailable = hν (ntrpntrs)Adoped, with ntrp, s = N0σ ap, s(σ ap, s + σ ep, s) , [10] can be 80  mJ/m, making it perfectly suitable for applications in high average power, high peak power-pulsed amplifiers.

Fig.  1. End-face of the large pitch fiber used and its several major modes.

A large core-diameter fiber amplifier with no interface can achieve high peak power output but the output beam quality would suffer a decrease when the output average power exceeds a certain value. The physical mechanism behind this is not quite clear, except for the fact that the influence of nonlinear effect is ruled out because the threshold has nothing to do with the peak power and is only influenced by the average power.[5] There are many explanations of this phenomenon.[1114] One explanation is that it is due to the formation of induced long-period gratings.[12] Here the interference pattern of the FM between an HOM is mapped into the population inversion and, therefore, locally modifies the refractive index of the core, resulting in coupling between FM and HOMs. But due to the short length of MSF in general application, the influence of induced long-period gratings is not significant. According to the power scaling estimation, [15]the thermally induced mode deformation is one of the most important factors to restrict the possible output power of MSF, especially for large mode area (LMA) fibers. The heat produced by the fiber changes the refractive index, causing the fiber modes to shift. Therefore, the newly generated modes couple with the original modes, resulting in the decrease of beam quality.[16] For this reason, the thermal effect introduces an independent limitation to avoid multimode output. Plenty of research has focused on the beam quality deterioration caused by thermally induced mode coupling, [13, 16, 17] but they have failed to combine thermally induced mode coupling with mode competition.

In this paper, we substitute the thermally induced mode coupling into mode competition for a 190-μ m core-diameter fiber amplifier at high power levels, acquiring the relationship between the average output power and the beam quality. The influence of thermally induced mode coupling on laser beam quality is also given. A schematic figure of the laser amplifier using the 190-μ m core-diameter fiber is shown in Fig.  2.

Fig.  2. Schematic diagram of the fiber amplifier.

Suppose that all of the modes can propagate along the fiber, then the population of the fiber core, the signal power, and the pump power of each mode will be given by the following space and time independent steady-state rate equations:[18]

where τ is the spontaneous lifetime of the upper lasing level; h is the Planck constant; ν s and ν p are the laser and pump frequencies, respectively; N1 and are N2 the population densities of the lower and upper lasing levels; N0 is the doping concentration distribution; and are the pump powers in the forward and backward directions, respectively; and are the signal powers of the i-th transverse mode in the forward and backward directions, respectively; σ α p (σ α s) and σ ep (σ es) are the pump absorption (emission) and signal absorption (emission) cross sections, respectively; α p and α si are the laser and pump loss factors of the i-th mode, respectively; φ p and φ si are the normalized intensity distributions of the pump and i-th mode; and, dij is the power coupling coefficient between the i-th and j-th modes.

The thermal effect appears when the output power is high enough to exceed the threshold, resulting in a change of the refractive index. According to Ref.  [16], the thermally induced refractive index change in MSF can be described as

where a and b are the radii of core and inner-clad, respectively. The heart density Q(z) is expressed as ; η is the quantum defect, which is described as η = 1 − λ p/λ s; α is the absorption coefficient of pump light, which is expressed as α = σ apntrpAdoped/Apump, with ntrp = N0σ ap/(σ ap + σ ep), Adoped and Apump being the areas of active region and pump cladding; β = dn/dT is the thermal– optic coefficient; and, H is the convection heat transfer coefficient. In order to investigate the complex thermal mechanics in LPF in a quick and simple way, besides utilizing the simplified average thermally induced index parameters, we replace the thermal conductivity κ by κ = (Asilicaκ silica + Aairκ air)/(Asilica + Aair) in the summary form, with Asilica and Aair being the areas of the silica material and air holes in the inner-cladding, respectively.

The thermally induced index change affects the mode profiles, which would change the final output beam quality of the fiber laser system. In an amplifier, the eigenmodes change along the fiber as the thermal perturbation caused by the residual pump power and the waste heat change. The mode coupling effect is a function of fiber parameters, pump direction, pump power, seed power, etc. The mode coupling coefficient between two modes can be described by the overlap integral of their normalized field distributions, and is defined as[19]

where Ei and Ej are the normalized field distributions of the i-th mode and j-th mode, respectively; and, A denotes the fiber cross section area.

Fig.  3. Mode coupling in heating process.

The heat-induced index change relates to pump power, while pump laser power differs at different positions of the gain fiber; therefore, the eigenmode distribution and the coupling coefficients change constantly along the fiber length direction. As shown in Fig.  3, taking the coupling between LP01 and LP21 for example, we describe the mode instability procedure. Dividing the fiber with length L into N segments, each with a length of Δ L, then according to Eqs.  (4) and (5) and COMSOL software, we can derive the and of the first segment Δ L1 after the index change has been induced by initial pump power. The coupling coefficient can be calculated according to Eq.  (6) by and and the seed laser modes LP01 and LP21. Substituting the obtained values into rate equations  (1)– (3), we can then derive the pump power and the signal power signal1 after propagating a length of Δ L1. We can then calculate and , and can be calculated by , , , and . As the seed laser power of the second segment, signal1 should be substituted into the second segment Δ L2 and introduce mode competition. The laser mode and power after mode competition of N segments can be derived this way. The equivalent output mode field diameter (MFD) and beam quality (M2) can then be calculated. In this paper, we only give the influence of the heat instability on output beam quality by adopting a relatively simple method, the results need to be supported and improved by further theoretical and experimental researches.

3. Results and discussion

The fiber laser’ s output beam quality and equivalent mode field diameter under a certain pump power can be calculated according to the mode coupling analysis. Then, by the limit of beam quality, we obtain the utmost pump power and the corresponding maximum output laser power. The mode coupling tends to happen between lower order modes, especially with FM. The calculation results show that coupling between LP21 and LP01 takes place most easily. So only the coupling between these two modes is taken into consideration in the following calculation. The calculation parameters are shown in Table  1. A forward pump is used in mode coupling because the pump light can be absorbed almost completely in this case; therefore, little residual heat is introduced, resulting in FM output without influence of heat.

Table 1. Parameter values used in the calculation of index change.

Figure  4 shows the LP01 and LP21 distributions after the heat has induced index change under different pump powers. We can obtain that, as the pump power is higher, the laser modes have a stronger tendency to shrink into the fiber core. This is because higher pump power results in higher residual heat and larger index change. On the other hand, the change in the fiber core is bigger than in cladding, resulting in the enlargement of equivalent numerical aperture (NA).

Fig.  4. Simulated laser modes after index change under different pump powers.

Fig.  5. Simulated results of the M2 and MFD of output laser under different pump powers.

Figure  5 gives the simulated results of the M2 and MFD of output laser under different initial pump powers, from which we can see that the MFD increases and the output beam quality decreases with pump power increasing. The situation deteriorates under the higher power operation. This is due to the mode instability accelerating under the higher pump power, then the mode coupling coefficient is higher under high power operation, resulting in LP21 accounting for a higher proportion in the output laser.

Taking M2 = 2.0 as the upper limit of the output beam quality, the threshold pump power is obtained to be 780  W. The simulated results of M2 and MFD of the signal at different positions of the fiber when pump power is Ppump= 780  W are shown in Fig.  6. The LP21 is excited due to mode coupling, which accounts for a higher proportion in the total output power; therefore, the M2 increases at the beginning. In the middle segments of the fiber we see the decrease of M2. This is due to the decrease of mode coupling coefficient and the shrink in proportion of LP21 in the output power with the decrease of pump power. In the end segment of the fiber the M2 increases again because the heat decreases with the decrease of pump power and the laser modes tends to spread out to the cladding as heat decreases. In this case, M2 takes a dominant position in laser modes at different positions of the fiber. The MFD increases along the fiber, also owing to the influence of pump power decreasing on laser mode distribution. The output laser power in fiber end is 540  W, with an available pulse energy of almost 100  mJ, which is suitable for pulse amplifier with a repetition rate of more than 5  kHz.

Fig.  6. Simulated results of M2 and MFD of the signal laser at different positions under Ppump = 780  W.

4. Conclusions

In this paper, we propose a method to analyze the mode coupling due to heat-induced index change. The method is based on the mode competition model, which is similar to the actual functioning mechanism of the gain fiber. We also predict the maximum average output power according to the mode coupling theory. For the large pitch fiber with a core diameter of 190  μ m, we analyze the influence of mode coupling on the MFD and beam quality of output laser. The result shows that, with the increase of the initial pump power, the fiber mode instability is enhanced, the MFD of output laser increases, and the beam quality declines. The beam quality is the main factor restricting the promotion of average power for large mode area fiber. In the case of M2 ≤ 2, the maximum average power reaches 540  W and the available pulse energy is 80  mJ/m, which is suitable for high energy, high average power pulse amplifier with a repetition rate of more than 5  kHz.

Reference
1 Jansen F, Stutzki F, Liem A, Jauregui C, Limpert J and Tünnermann A 2012 Advanced Solid-State PhotonicsOptical Society of America, AT1A. 4 [Cited within:1]
2 Stutzki F, Jansen F, Liem A, Jauregui C, Limpert J and Tünnermann A 2012 Opt. Lett. 37 1073 DOI:10.1364/OL.37.001073 [Cited within:2] [JCR: 3.385]
3 Eidam T, Rothhardt J, Stutzki F, Jansen F, Hädrich S, Carstens H, Jauregui C, Limpert J and Tünnermann A 2011 Opt. Express 19 255 DOI:10.1364/OE.19.000255 [Cited within:1] [JCR: 3.546]
4 Stutzki F, Jansen F, Jauregui C, Limpert J and Tünnermann A 2011 Opt. Lett. 38 97 [Cited within:1] [JCR: 3.385]
5 Eidam T, Wirth C, Jauregui C, Stutzki F, Jansen F, Otto H J and Tünnermann A 2011 Opt. Express 19 13219 [Cited within:2] [JCR: 3.546]
6 Wirth C, Schreiber T, Rekas M, Tsybin I, Peschel T, Eberhardt R and Tünnermann A 2011SPIE LASE International Society for Optics and Photonics75801H-75801H-6 [Cited within:1]
7 Eidam T, Hädrich S, Jansen F, Stutzki F, Rothhardt J, Carstens H and Tünnermann A 2011 Opt. Express 19 8656 DOI:10.1364/OE.19.008656 [Cited within:1] [JCR: 3.546]
8 Stutzki F, Jansen F, Eidam T, Steinmetz A, Jauregui C, Limpert J and Tünnermann A 2011 Opt. Lett. 36 689 DOI:10.1364/OL.36.000689 [Cited within:1] [JCR: 3.385]
9 Wan P, Yang L M and Liu J 2013 Opt. Express 21 29854 DOI:10.1364/OE.21.029854 [Cited within:1] [JCR: 3.546]
10 Limpert J, RÖser F, Schimpf D N, Seise E, Eidam T, Hädrich S and Rothhardt J 2009 IEEE J. Sel. Top. Quantum Electron. 15 159 DOI:10.1109/JSTQE.2008.2010244 [Cited within:1] [JCR: 4.078]
11 Jiang Z and Marciante J R 2008 J. Opt. Soc. Am. B 25 247 DOI:10.1364/JOSAB.25.000247 [Cited within:1] [JCR: 2.21]
12 Jauregui C, Eidam T, Limpert J and Tünnermann A 2011 Opt. Express 19 3258 DOI:10.1364/OE.19.003258 [Cited within:1] [JCR: 3.546]
13 Smith A V and Smith J J 2011 Opt. Express 19 10181 [Cited within:1] [JCR: 3.546]
14 Digonnet M J F, Sadowski R W, Shaw H J and Pantell R H 1997 Opt. Fiber Technol. 3 44 DOI:10.1006/ofte.1997.0189 [Cited within:1] [JCR: 1.187]
15 Dawson J W, Messerly M J, Beach R J, Shverdin M Y, Stappaerts E A, Sridharan A K, Pax P H, Heebner J E, Siders C W and Barty C P J 2008 Opt. Express 16 13240 DOI:10.1364/OE.16.013240 [Cited within:1] [JCR: 3.546]
16 Brown D C and Hoffman H J 2005 Thermal IEEE J. Quantum Electron. 37 207 [Cited within:3]
17 Ke W, Wang X, Bao X and Shu X 2013 Opt. Express 21 14272 DOI:10.1364/OE.21.014272 [Cited within:1] [JCR: 3.546]
18 Gong M, Yuan Y, Li C, Yan P, Zhang H and Liao S 2007 Opt. Express 15 3236 DOI:10.1364/OE.15.003236 [Cited within:1] [JCR: 3.546]
19 Okamoto K 2006 Fundamentals of Optical Waveguides2nd edn San Diego Academic [Cited within:1]