Electronic and thermal refractive index changes in ... - OSA Publishing

Electronic and thermal refractive index changes in Ytterbium-doped fiber amplifiers M. S. Kuznetsov,1∗ O. L. Antipov,1,5 A. A. Fotiadi,2,3,4 and P. Mégret2 1. Department of Nonlinear Dynamics and Optics Division, Institute of Applied Physics RAS, 46 Ul’yanov Street, 603950, Nizhny Novgorod, Russia 2. University of Mons, Electromagnetism and Telecommunication Department, 31 Boulevard Dolez, B-7000 Mons, Belgium 3. Ioffe Physico-Technical Institute of the RAS, 26, Polytekhnicheskaya Str., St. Petersburg 194021, Russia 4 Ulyanovsk State University, 42, Lev Tolstoy Str., Ulyanovsk 432970, Russia 5 Nizhny Novogorod State University, 23, Gagarin Ave., Nizhny Novgorod 603950, Russia ∗

[email protected]

Abstract: We develop a theoretical framework to analyze the mechanism of refractive index changes (RIC) in double-clad Yb3+ doped optical fibers under resonant core or clad pumping, and with signal amplification. The model describes and compares thermal and electronic contributions to the phase shifts induced on the amplified signal at 1064 nm and the probe signal at 1550 nm, i.e. located inside and outside of the fiber amplification band, respectively. The ratio between the thermal and electronic phase shifts is evaluated as a function of the pump pulse duration, the gain saturation, the amplified beam power and for a variety of fiber parameters. © 2013 Optical Society of America OCIS codes: (060.2300) Fiber measurements; (060.2400) Fiber properties; (060.3510) Lasers, fiber; (140.3510) Lasers, fiber; (140.6810) Thermal effects.

References 1. M. J. F. Digonnet, R. W. Sadowski, H. J. Shaw, and R. H. Pantell, “Resonantly enhanced nonlinearity in doped fibers for low-power all-optical switching: A review,” Opt. Fiber Technol. 3, 44–64 (1997). 2. E. Desurvire, Erbium-doped Fiber Amplifiers: Principles and Applications (Wiley, 1994). 3. J. W. Arkwright, P. Elango, G. R. Atkins, T. Whitbread, and M. F. Digonnet, “Experimental and theoretical analysis of the resonant nonlinearity in ytterbium-doped fiber,” J. Lightwave Technol. 16, 798–806 (1998). 4. A. Fotiadi, O. Antipov, and P. Mégret, “Dynamics of pump-induced refractive index changes in single-mode Yb-doped optical fibers,” Opt. Express 16, 12658–12663 (2008). 5. D. J. Richardson, J. Nilsson, and W. A. Clarkson, “High power fiber lasers: current status and future perspectives,” J. Opt. Soc. Am. B 27, 63–92 (2010). 6. 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). 7. A. Smith and J. Smith, “Mode instability in high power fiber amplifiers,” Opt. Express 19, 10180–10192 (2011). 8. C. Jauregui, T. Eidam, J. Limpert, and A. Tünnermann, “The impact of modal interference on the beam quality of high-power fiber amplifiers„” Opt. Express 19, 3258–3271 (2011). 9. 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). 10. V. Spirin, C. López-Mercado, D. Kinet, P. Mégret, I. Zolotovskiy, and A. Fotiadi, “Single longitudinal-mode brillouin fiber laser passively stabilized at pump resonance frequency with dynamic population inversion grating,” Laser Phys. Lett. 10, 015102 (2013).

#193649 - $15.00 USD Received 11 Jul 2013; revised 31 Aug 2013; accepted 5 Sep 2013; published 16 Sep 2013 (C) 2013 OSA 23 September 2013 | Vol. 21, No. 19 | DOI:10.1364/OE.21.022374 | OPTICS EXPRESS 22374

11. S. Turitsyn, A. Bednyakova, M. Fedoruk, A. Latkin, A. Fotiadi, A. Kurkov, and E. Sholokhov, “Modeling of cw Yb-doped fiber lasers with highly nonlinear cavity dynamics,” Opt. Express 19, 8394–8405 (2011). 12. H. Bruesselbach, D. C. Jones, M. S. Mangir, M. Minden, and J. Rogers, “Self-organized coherence in fiber laser arrays,” Opt. Lett. 30, 13–15 (2005). 13. T. Fan, “Laser beam combining for high-power, high-radiance sources,” IEEE J. Sel. Top. Quantum Electron. 11, 567–577 (2005). 14. H. Bruesselbach, S. Wang, M. Minden, D. Jones, and M. Mangir, “Power-scalable phase-compensating fiberarray transceiver for laser communications through the atmosphere,” J. Opt. Soc. Am. B 22, 347–353 (2005). 15. A. Fotiadi, N. G. Zakharov, O. Antipov, and P. Mégret, “All-fiber coherent combining of Er-doped amplifiers through refractive index control in Yb-doped fibers,” Opt. Lett. 34, 3574–3576 (2009). 16. A. Fotiadi, O. Antipov, M. Kuznetsov, K. Panajotov, and P. Mégret, “Rate equation for the nonlinear phase shift in Yb-doped optical fibers under resonant diode-laser pumping,” J. Hologr. and Speckle 5, 1–4 (2009). 17. A. Fotiadi, O. Antipov, and P. Mégret, “Resonantly induced refractive index changes in Yb-doped fibers: the origin, properties and application for all-fiber coherent beam combining,” in Frontiers in Guided Wave Optics and Optoelectronics, B. Pal, ed. (INTECH, 2010), chap. 11, pp. 209–234. 18. H. Tünnermann, J. Neumann, D. Kracht, and P. Weßels, “Gain dynamics and refractive index changes in fiber amplifiers: a frequency domain approach,” Opt. Express 20, 13539–13550 (2012). 19. A. Fotiadi, O. Antipov, M. Kuznetsov, and P. Mégret, “Refractive index changes in rare earth-doped optical fibers and their applications in all-fiber coherent beam combinig,” in Coherent Laser Beam Combining, A. Brignon, ed. (John Wiley & Sons, 2013), chap. 7, pp. 193–230. 20. S. Stepanov, A. Fotiadi, and P. Mégret, “Effective recording of dynamic phase gratings in Yb-doped fibers with saturable absorption at 1064 nm,” Opt. Express 15, 8832–8837 (2007). 21. M. Davis, M. Digonnet, and R. Pantell, “Thermal effects in doped fibers,” J. Lightwave Tech. 16, 1013–1023 (1998). 22. D. Brown and H. Hoffman, “Thermal, stress, and thermo-optic effects in high average power double-clad silica fiber lasers,” IEEE J. Quantum Electron. 37, 207–217 (2001). 23. O. Antipov, O. Eremeykin, A. Savikin, V. Vorob’ev, D. Bredikhin, and M. Kuznetsov, “Electronic changes of refractive index in intensively pumped Nd:YAG laser crystals,” IEEE J. Quantum Electron. 39, 910–918 (2003). 24. R. Soulard, R. Moncorgé, A. Zinoviev, K. Petermann, O. Antipov, and A. Brignon, “Nonlinear spectroscopic properties of Yb3+ -doped sesquioxides Lu2 O3 and Sc2 O3 ,” Opt. Express 18, 11173–11180 (2010). 25. V. Gainov and O. Ryabushkin, “Effect of optical pumping on the refractive index and temperature in the core of active fibre,” Quantum Electronics 41, 809–814 (2011). 26. H. Tünnermann, J. Neumann, D. Kracht, and P. Weßels, “All-fiber phase actuator based on an erbium-doped fiber amplifier for coherent beam combining at 1064 nm,” Opt. Lett. 36, 448–450 (2011). 27. M. Bass, E. Van Stryland, D. Williams, and W. Wolfe, Handbook for Optics, 2nd ed. (MGH, 1995). 28. V. Privalko, Handbook for Physical Chemistry of Polymers (Naukova Dumka, Kiev, 1984). 29. M. Melkumov, I. Bufetov, K. Kravtsov, A. Shubin, and E. Dianov, Cross Sections of Absorption and Stimulated Emission of Yb3+ Ions in Silica Fibers Doped with P2 O5 and Al2 O3 (FORC, Moscow, 2004). 30. A. Fotiadi, O. Antipov, I. Bufetov, E. Dianov, and P. Mégret, “Comparative study of pump-induced refractive index changes in aluminum and phosphate silicate Yb-doped fibers,” in Conference on Lasers and ElectroOptics/Quantum Electronics and Laser Science Conference and Photonic Applications Systems Technologies, (Optical Society of America, Washington, D.C., 2009), paper JWA9. 31. J. W. Dawson, M. J. Messerly, R. J. Beach, M. Y. Shverdin, E. A. Stappaerts, A. K. Sridharan, P. H. Pax, J. E. Heebner, C. W. Siders, and C. Barty, “Analysis of the scalability of diffraction-limited fiber lasers and amplifiers to high average power,” Opt. Express 16, 13240–13266 (2008). 32. G. Agrawal, Nonlinear Fiber Optics, 3rd ed. (Academic Press, 2001). 33. H.-G. Unger, Planar Optical Waveguides and Fibres (Oxford University Press, 1977). 34. E. Tugolukov, Solution of Problems of Thermal Conductivity by Finite Integral Transforms: a Tutorial (TSTU, Tambov, 2005). 35. H. Carslow and J. Jaeger, Conduction of Heat in Solids (Science, Moscow, 1964). 36. S. Stepanov and M.P. Sánchez, “Phase population gratings recorded in ytterbium doped fiber at 1064 nm,” in 22nd Congress of the International Commission for Optics: Light for the Development of the World, R. Rodríguez-Vera and R. Díaz-Uribe, eds., Proc. SPIE 8011, 801153 (2011). 37. A. Luikov, Analytical Heat Diffusion Theory (Academic Press, 1968).

1.

Introduction

The understanding of the origin of the refractive index changes (RICs) in intensively-pumped rare-earth doped fibers appears to be very important because the pump- and signal-induced index changes could significantly impact the dynamics of the fiber lasers and amplifiers, the fiber mode structures, the fiber Bragg gratings, and so on [1–11]. Besides, the enhanced nonlinear #193649 - $15.00 USD Received 11 Jul 2013; revised 31 Aug 2013; accepted 5 Sep 2013; published 16 Sep 2013 (C) 2013 OSA 23 September 2013 | Vol. 21, No. 19 | DOI:10.1364/OE.21.022374 | OPTICS EXPRESS 22375

phase shift could be used for optical switching [1, 3], coherent beam combining [12–19] and adaptative interferometry [20]. Two main mechanisms of the RIC have been discussed recently: the well-known thermal index change δ nT caused by the heating owing to the matrix absorption and quantum defects [1,2,5,21,22], and the athermal index change δ nN (also called electronic or "Kramers-Kroning" effect) related to the population change of the active-ion levels with different polarizabilities [1, 3, 4, 16–18]. The latter effect has been previously investigated mainly for crystal lasers [23, 24] and now attracts a great interest for fiber lasers and amplifiers [7, 8, 19, 25, 26]. This paper is devoted to a comparative study of the thermal and electronic RICs in optical fiber doped by Yb3+ -ions leading to phase shifts of the probe beams at 1064 nm and 1550 nm, i.e. inside and outside the fiber amplification band, respectively, under different pumping conditions, different amplified beam powers and for a variety of the fiber parameters. Theoretical calculations are in good agreement with our experimental results reported earlier [4,15,17,19,20]. 2.

Description of the refractive index changes in the optical fibers

For analyzing the RIC and the phase shift of a probe beam in an active fiber, the fiber is described as in Fig. 1 by the infinite composite cylinder consisting of a Yb3+ -doped core (green), a silica-glass pumping cladding (blue) and a plastic cladding (red). The following fiber parameters have been used for calculations: the core radius, r0 is 1.8 µm; the glass radius, r1 is 62.5 µm; the external plastic radius, r2 is 125 µm; Yb-ion doping concentration in the core is 8.56 × 1019 cm−3 ; the fiber length, l is 2 m. Other physical parameters of the fibers are given in Table 1. Only calculation results related to aluminum-silicate (AS) fibers are discussed here. The results for phosphate-silicate (PS) fibers can be obtained with the same formalism using the corresponding data shown in Table 1.

r2 r1

2r0

Figure 1. Scheme of the cross-section of the fiber light guide. r1 is the glass radius, r0 is the core radius activated by Yb3+ -ions, red part denotes the layer of plastic r1 < r < r2 .

The phase shift of the probe beam at λ in the optical fiber of length l associated with the RIC induced by the pumping beam at λ p ∼ 976 nm in the presence or absence of an amplified beam at λA ∼ 1064nm can be described by the following expression [32] as a fixed testing-beam structure is assumed: ´l´∞ ˆ l k0 0 0 ∆n(z, r,t) | G(r) |2 rdrdz ´∞ (1) ∆ϕ = ∆β (z)dz = 2 0 0 | G(r) | rdr

where ∆β = k0 ∆n is the variation of the propagation constant β (z) associated with the RIC, k0 = 2π /λ , the refractive index change ∆n is the sum of the thermal and electronic contributions, ∆n(z, r,t) = δ nT (z, r,t) + δ nN (z, r,t), G(r) is the radial distribution of the probe beam field [33]: J0 (ur/r0 ) if r ≤ r0 J0 (u) K0 (wr/r0 ) if r > r0 G(r) = A K0 (w)

G(r) = A

(2)


Table 1. Parameters values used in calculations

Parameter Glass density, ρ1 , g/cm3 Specific heat capacity of glass, c p1 , cal/(g K) (1 cal= 4.1868 J) Thermal conductivity of glass, K1 , W/(cm K) Plastic density, ρ2 , g/cm3 Specific heat capacity of plastic, c p2 , cal/(g K) Thermal conductivity of plastic, K2 , W/(cm K) Energy allocated to heat from the upper therm 2 F5/2 of Yb3+ ions in glass, cm−1

PS fiber 2.2 0.188

AS fiber 2.2 0.188

Reference, note [27] [27]

0.014 1.19 0.3 0.002 0

0.014 1.19 0.3 0.002 0

Energy allocated to heat from the low therm 2F7/2 of Yb3+ ions in glass, cm−1 Life time of the level 2F5/2 τ1 , ms Emission cross-section at 976 nm, σ21 (ν p ), pm2 Absorption cross-section at 976 nm, σ12 (ν p ), pm2 Emission cross-section at 1.01 µm, σ21 (νL ), pm2 Absorption cross-section at 1.01 µm, σ12 (νL ), pm2 Emission cross-section at 1064 nm, σ21 (νA ), pm2 Absorption cross-section at 1064 nm, σ12 (νA ), pm2 Polarizability difference at 1550 nm, ∆p, cm3 Polarizability difference at 1060 nm, ∆p, cm3 Glass refractive index, n0 Derivative of refractive index of glass with respect to temperature, ∂ n/∂ T , K−1 Thermal expansion coefficient of glass δT , K−1 Poisson’s coefficient of glass, b

740

1210

[5, 27] [27] [28] [28] Pumping in the low level of the therm [29]

1.276 1.46 1.38 0.485 0.068

0.83 2.97 2.69 0.495 0.086

[29] [29] [29] [29] [29]

0.13 0.0016

0.3 0.0046

[29] [29]

0.9 × 10−26 1.28 × 10−26 1.5 1.2 × 10−5

1.2 × 10−26 1.74 × 10−26 1.5 1.2 × 10−5

[16, 19, 30]

5.1 × 10−7 0.164

5.1 × 10−7 0.164

[27] [27]

[27] [31]

where u and w are defined by: u2 = r02 (n20 k02 − β 2 );

w2 = r02 (β 2 − n21k02 );

V 2 = u2 + w2 = k02 r02 [n20 − n21]

(3)

and u and w are linked by the characteristic Eq.: uJ1 (u) wK1 (w) = J0 (u) K0 (w)

(4)

In these relations, Ji is the Bessel function of order i, Ki is the modified Bessel function of order i, A is the normalization factor, ni is the refractive index of the core (i = 0) and the cladding (i = 1). 2.1. Thermal contribution to refractive index changes The thermal contribution to the RIC in the fiber can be defined by the expression:

δ nT (z, r,t) =

∂n δT ∂T

(5)


where ∂ n/∂ T is the thermo-optic coefficient of the silica glass, and δ T (z, r,t) is the temperature variation of the fiber. The temperature evolution in the fiber under pumping and in the presence of amplified signal is described by the following Eq.:

∂T − a2i ∇2 T = Q(z, r,t) ∂t

(6)

where a2i = Ki /(ρi c pi ) is the thermal diffusivity, Ki is the thermal conductivity, ρi and c pi are the density and heat capacity (i = 1 for glass and i = 2 for plastic), Q(z, r,t) is the heat source inside the fiber. In the case of pumping at 976 nm (in the lowest sublevel of 2F5/2 -state) and amplification signal at 1064 nm, the heat source can be described by the sum:

αAnr IA + αPnr IP hνB1 N2 (z, r,t) νBL N2 (z, r,t)IL σ21 (νL ) + + ρ1 c p1 ρ1 c p1 τ1 ρ1 c p1 νL νB3 {[σ21 (νA ) + σ12(νA )] N2 (z, r,t) − σ12 (νA )N0 (z, r,t)} IA + ρ1 c p1 νA

Q(z, r,t) =

(7)

where h is the Planck’s constant, νB1 is the frequency of the nonradiative transition between sublevels of the ground state, N2 (z, r,t) is the population of the excited level, τ1 is the relaxation time of the excited level (0.83 ms and 1.245 ms for AS and PS fibers, respectively [29]), N0 is the total population of excited and unexcited ions, I p is the intensity of the pumping beam, ν p is the pumping frequency, IL is the luminescence intensity, IA is the amplified beam intensity, νA is the amplified beam frequency, σ21 (ν ) is the emission cross-section at frequency ν , σ12 (ν ) is the absorption cross-section at frequency ν , νL is the effective frequency of the luminescence corresponding to λL ∼ 1.01 µm , νB3 = ν p − νA , νBL = ν p − νL , αPnr and αLnr are the absorption coefficients of the pumping and amplified beams in the silica glass host, respectively. The first summand in (7) corresponds to the gray losses in the matrix, and the other three summands describe to thermalization of the pump, luminescence, and amplified beam respectively. In the case of a linear heat transfer from plastic to air, the boundary condition on the outer surface of the plastic cylinder can be described by the Newton’s law: ∂ T (r,t) (8) + η [T (r,t) − Tair ] = 0 K2 ∂ r r2 r2 where η is the heat transfer coefficient from plastic to air. The increase of temperature δ T (z, r,t) of the fiber (inner cladding and core regions) can be found from (6) - (8) and written by the following expression [34]: ˆ r0 ˆ t ∞ 1 K1 µn r µn r ′ J J0 exp −µn2 (t − t ′) Q(z, r′ ,t ′ )dt ′ r′ dr′ (9) δ T (z, r,t) = ∑ 2 0 Z a a a 1 1 0 0 n=1 n 1

where J0 is the Bessel function of the zero order, r is the transverse (radial) coordinate of the fiber, Zn is described by the following expression: K2 K1 r12 2 2 J0 (Ψn,1,1 ) + J1 (Ψn,1,1 ) + 2 0.5r22C22 (µn ) J02 (Ψn,2,2 ) + J12 (Ψn,2,2 ) Zn = 2 2a1 a2

+r22C2 (µn )D2 (µn ) [J0 (Ψn,2,2 )Y0 (Ψn,2,2 ) + J1 (Ψn,2,2 )Y1 (Ψn,2,2 )] +0.5r22 D22 (µn ) Y02 (Ψn,2,2 ) + Y12 (Ψn,2,2 ) − 0.5r12C22 (µn ) J02 (Ψn,1,2 ) + J12 (Ψn,1,2 ) −0.5r12C2 (µn )D2 (µn ) [J0 (Ψn,1,2 )Y0 (Ψn,1,2 ) + J1 (Ψn,1,2 )Y1 (Ψn,1,2 )]


−0.5r12 D22 (µn )

2 Y0 (Ψn,1,2 ) + Y12 (Ψn,1,2 )

(10)

where ψn,i, j = µn ri /a j and µn is the n-th positive root of the Eq.: µn K2 µn K2 C2 (µn ) J0 (Ψn,2,2 ) − J1 (Ψn,2,2 ) + D2 (µn ) Y0 (Ψn,2,2 ) − Y1 (Ψn,2,2 ) = 0 a2 η a2 η

(11)

with: J0 (Ψn,1,1 ) − D2 (µn )Y0 (Ψn,1,2 ) J0 (Ψn,1,2 ) (K1 a2 /K2 a1 )J0 (Ψn,1,2 )J1 (Ψn,1,1 ) − J1(Ψn,1,2 )J0 (Ψn,1,1 ) D2 ( µ n ) = J0 (Ψn,1,2 )Y1 (Ψn,1,2 ) − J1 (Ψn,1,2 )Y0 (Ψn,1,2 ) C2 (µn ) =

(12) (13)

In these relations, J1 , Y0 and Y1 are the Bessel function of the first order, Neumann functions of the zero and first order, respectively. When the fiber is pumped, generation of heat is created inside the core and the cladding. This leads to a modification of the refractive index through the thermo-optic effect. But as the temperature is increased, the fiber length and section will expand and these secondary effects will also lead to phase shifts induced on the probe beam. A careful analysis is thus necessary to compare the relative contributions of these two secondary thermal effects with the thermooptic effect acting directly on the refractive index n. Let us define δ ϕ1 , δ ϕ2 and δ ϕ3 , the phase shift contributions arising respectively from the local thermo-optic coefficient, the transverse thermal expansion, and the longitudinal thermal expansion. The thermal contribution δ ϕ1 associated with the thermo-optic coefficient is described by the formula: ´l´∞ ∂ n 0 0 δ T (z, r,t) | G(r) |2 rdrdz ´∞ (14) δ ϕ1 = k0 2 ∂T 0 | G(r) | rdr To calculate the fiber deformation, it is necessary to solve the problem of thermoelasticity in a quasi-static formulation with a given temperature field in the approximation of plain strain. The effect of radial thermal expansion on the phase shift can be accounted for as follows: increasing the radius of the fiber r0 leads to a change of the mode and the wavenumber β . Then the phase shift δ ϕ2 associated with this expansion is expressed as: ˆ l ˆ r0 ˆ l 2 ∂ β ∂ r0 ∂β δ ϕ2 = Tav (z,t)dz = δT (1 + b)r0 δ T (r, z,t)rdrdz (15) 2 ∂ r0 0 r0 0 0 ∂ r0 ∂ Tav

where δT is the coefficient of thermal expansion of glass, b is the Poisson’s coefficient of glass, ˆ 2 r0 Tav (z,t) = 2 δ T (r, z,t)rdr (16) r0 0 is the average temperature over the fiber core cross-section, s u2 β = k02 n20 − 2 r0

∂β = ∂ r0 [J02 (u)+J12 (u)]r02 2u2

+

J1 (u)K0 (w) 2

r0 2π (n0 − n1)/λ wK1 (w) [K02 (wr1 /r0 )−K12 (wr1 /r0 )]r12 +[K12 (w)−K02 (w)]r02 J12 (u)

(17)

(18)

2w2 K12 (w)


The fiber elongation makes the following contribution to the phase change:

δ ϕ3 = δT β

ˆ

0

l

2 r12

ˆ

r1

δ T (r, z,t)rdrdz

(19)

0

where δT is the thermal expansion coefficient of glass. Estimations show that contributions to the phase shift caused by the elongation δ ϕ3 and the expansion δ ϕ2 of the fiber are much smaller than the contribution δ ϕ1 due to changes in refractive index under the influence of temperature. In fact, ´r ´r δT 0 1 δ T (r, z,t)rdr 0 1 | G(r) |2 rdr δ ϕ2 < δ ϕ3 ≈ δ ϕ1 ∂ n 2 ´ r1 (20) 2 ∂ T r1 0 δ T (r, z,t) | G(r) | rdr

so the ratio δ ϕ3 /δ ϕ1 is smaller by more than 2 orders of magnitude at used parameters, and the contribution associated with the thermo-optic coefficient is the main thermal contribution. 2.2. Electronic (population) contribution to refractive index changes

The athermal (electronic) component of the RIC of doped silica fibers δ nN at the presence of the population change δ N of Yb3+ ions can be described by the following expression [23]:

δ nN =

∂n δN ∂N

(21)

where ∂ n/∂ N = (2π /n0)FL2 ∆p; FL = (n20 + 2)/3 is the factor of local field (the Lorenz factor), n0 is the undisturbed refractive index and ∆p(λ ) is the polarizability difference of medium particles at the probe wavelength expressed as: ∆p(λ ) =

n0 4π 3 FL2

∞ 0

λ2 λ ′2 − λ 2

σ12 (λ ′ ) + σ21(λ ′ ) − σesa (λ ′ ) dλ ′

(22)

ffl ∞ where 0 ()dλ ′ stands for the Cauchy principal part of the integral, σesa (λ ) is the excited-state absorption cross-section at λ . In Yb-doped materials, the well-allowed UV transitions to the 5d-electron shell and the charge-transfer transition are characterized by oscillator forces that are several orders of magnitude higher than the forces of optical transitions inside the 4f-electron shell. As a consequence, in the IR spectrum band, the polarizability difference ∆p(λ ) is expressed from Eq. (22) as a sum of contributions from the near-resonance transitions (between the ground and excited states) and non-resonance UV transitions [17]. The calculated spectrum of the polarizability difference in Yb-doped fiber within the amplification band is presented in Fig. 2 where one can see that within the Yb-fiber amplification band (around 1060 nm) the non-resonant contributions to the polarizability difference are dominating over the resonant parts. This fact has been confirmed in the experiment [20] for aluminum-silicate fibers. System of levels of ytterbium ions in the glass can be regarded as quasi two-level system provided thermodynamic equilibrium of the population distribution over the sublevels set in each multiplet is assumed. But the description of this system as a two-level one, leads to the fact that absorption and stimulated emission cross-sections are significantly different functions of the wavelength. Under this assumption and assumption of uniform distribution of the activation, the average population N¯ 2 (z) of the metastable level, given by: ˆ 2 r0 N2 (z, r)rdr (23) N¯ 2 (z) = 2 r0 0 #193649 - $15.00 USD Received 11 Jul 2013; revised 31 Aug 2013; accepted 5 Sep 2013; published 16 Sep 2013 (C) 2013 OSA 23 September 2013 | Vol. 21, No. 19 | DOI:10.1364/OE.21.022374 | OPTICS EXPRESS 22380

Polarizability difference, 10−27 cm3

50 40 30 20 10 0 −10 −20 −30 −40 840

880

920

960

1000 1040 1080 1120 1160 1200

Wavelength, nm

Figure 2. The difference of Yb3+ -ion polarizabilities in excited and ground states for phosphate silicate fibers (non resonant part in cyan, resonant part in green and total part in blue) and for aluminum-silicate fibers (non resonant part in gray, resonant part in red and total part in orange).

and is computed from the expressions: ∂ N¯ 2

N¯ 2 + ∂t τ1

Pp Score = σ12 (ν p ) (Γ p0 N0 − Γ p2 N¯ 2 ) hν p P P PL p L −σ21 (ν p )Γ p2 N¯ 2 − σ21 (νL )ΓL2 N¯ 2 + σ12 (νL ) (ΓL0 N0 − ΓL2 N¯ 2 ) hν p hνL hνL PA PA −σ21 (νA )ΓA2 N¯ 2 + σ12 (νA ) (ΓA0 N0 − ΓA2N¯ 2 ) hνA hνA

∂ Pp = − [(Γ p0 N0 − Γ p2N¯ 2 )σ12 (ν p ) − Γ p2N¯ 2 σ21 (ν p )] Pp − αPnr Pp ∂z

(24) (25)

∂ PL (26) = ΓL2 N¯ 2 σ21 (νL )PL + N¯ 2 ζ − (ΓL0 N0 − ΓL2 N¯ 2 )σ12 (νL )PL − αLnr PL ∂z ∂ PA (27) = ΓA2 N¯ 2 σ21 (νA )PA − (ΓA0N0 − ΓA2 N¯ 2 )σ12 (νA )PA − αAnr PA ∂z where, Pp , PL and PA are powers of the pump, the luminescence and the amplified beam, respectively, Score is the area of the core, αPnr , αLnr and αAnr are coefficients of nonresonant losses for the pump, luminescence and amplified waves, respectively. ζ is the Langevin source of the luminescence: hνL Score 2r0 2 (28) ζ= l τ1 4π where νL is the effective luminescence frequency. In these expressions, the parameter Γ p,L,A; j , is defined as follow: ´ r0 N j (z, r)I p,L,A (r)rdr r2 ´ Γ p,L,A; j (z) = 0 ´ ∞ 0 (29) 2 0 I p,L,A (r)rdr 0r0 N j (z, r)rdr

if I p,L,A is the intensity of the pump, luminescence and amplified beams respectively. Calculations showed that one can neglect the luminescence power PL in solving this problem, so it is no longer taken into account in the following theoretical estimations.


Smode hν p,A

sat = For the fiber under simulations, the saturation powers defined as Pp,A

[σ12 (ν p,A )+σ21 (ν p,A )]τ1 are 530 mW, 0.41 mW and 8.2 mW for the clad pumping at 976 nm, and core pumping at 976 nm and 1064 nm, respectively. 3.

Comparison of thermal and electronic contributions to the phase shift of a probing beam during pumping without strong amplified beam.

The ratio of the thermal and electronic contributions to RIC was found to depend on the powers of the pump and amplified beams. In this section we consider this ratio in the case of a small amplified beam power, so in the limit PA → 0, without saturation of the population inversion. Investigation of the thermal and electronic contributions to the phase shift of the fiber was carried out using analytical and numerical method. Coefficients µn of Eq. (9) were numerically calculated by solving Eq. (11) and then Eq. (9) was numerically integrated by the predictorcorrector finite-difference method scheme. Calculations showed that, depending on the conditions, thermal or electronic contributions can be the dominant effect. Refractive index change, 10−5

100

10−1

10−2

10−3 10−1

plastic cladding

core 100

101

air 102

103

Transverse coordinate, µm

Figure 3. Distribution of refractive index changes (RIC) at λ = 1550 nm and z0 = 20 cm in the transverse coordinate shown on logarithmic scales. The heat source is assumed to be uniformly distributed in the fiber core (pumping in the clad of 145 mW). Orange curve is the steady-state electronic contribution to RIC, green curve is the thermal contribution to RIC at t = 0.5 s, and red curve is the steady-state thermal contribution to RIC.

Graphics of the distribution of RIC on the distance from the fiber axis are shown on Fig. 3. The orange trace gives the steady-state electronic RIC contribution and one can see a modification of the core refractive index of 9 × 10−7. The thermal RIC contribution is given after 0.5 s in the green trace and in thermal steady-state in the red trace. One can see that the core refractive index change due to the thermal effect is 2 × 10−7 after 0.5 s, so much lower than the electronic contribution and then reaches 8 × 10−6 in steady-state, which is here much higher than the electronic contribution. Fig. 3 clearly illustrates that there should exist a time for which the two contributions are equal: this is the alignment time estimated to 0.757 s as explained below. In dynamics, the ratio of the thermal δ ϕT and the electronic δ ϕN contributions could be characterized through the alignment time tm that is defined as the time needed for these two contribution to become equal: ˆ lˆ 0

0

∞

δ nT (z, r,tm ) | G(r) |2 rdrdz =

ˆ lˆ 0

∞

δ nN (z, r,tm ) | G(r) |2 rdrdz

(30)

0

The analytical estimate of the alignment time of electronic and thermal contributions tm is obtained from (30) under assumptions of rectangularly switching in time heating Q(t), #193649 - $15.00 USD Received 11 Jul 2013; revised 31 Aug 2013; accepted 5 Sep 2013; published 16 Sep 2013 (C) 2013 OSA 23 September 2013 | Vol. 21, No. 19 | DOI:10.1364/OE.21.022374 | OPTICS EXPRESS 22382

K2 /η r2 ≫ 1 and αPnr = αLnr = αAnr = 0. After some computations, tm is found to be: ( ) r2 K2 (1 − p) 4π FL2 ∆pτ1 r2 ηξ (1 + p) tm = − ln 1 − 2a22 η n0 ∂∂Tn r02 hvB1 (1 − p)

(31)

where: r12 ρ1 c p1 −1 p= 2 r2 ρ2 c p2 ´ l ´ r0 N2 (z, r) | G(r) |2 rdrdz r02 ξ = ´∞ 0 0 ´ l ´ r0 2 | G(r) |2 rdr N2 (z, r)rdrdz 0

(32) (33)

0 0

If Q(r) is uniform over the cross-section of the fiber core (pumping to the clad), i.e. it does not depend on r, then relation (33) implies ξ = 0.38 at λ = 1550nm and ξ = 0.69 at λ = 1060 nm. It is worth to notice that (31) is the upper limit of the equalizing time, i.e. the time for which the thermal contribution accounted through the term n = 1 in the expression (9) becomes equal to the electronic contribution to the phase shift. To observe the alignment of thermal and electronic contributions to the phase shift in a finite time, it is necessary that the ratio ε of the electronic steady-state contribution and thermal steady-state contribution to the phase shift is less than one: −1 ∂n 2 1 ln(r2 /r1 ) 1/2 + ln(r1 /r0 ) + ε = 4π FL2 ∆pτ1 ξ n0 r0 hνB1 + ηcr ) the RIC electronic contribution always predominates over the thermal contribution even in steady-state. When one substitutes the parameters mentioned above (included in Eq. (35)), ηcr takes the value 0.003 cal/(cm2 K s) and this value clearly obeys K2 /(ηcr r2 ) ≈ 12.3 ≫ 1. This result suggests that, when the heat transfer coefficient η < ηcr , one can use Eq. (31) to find the alignment time of the contributions. For the calculations presented here, the heat transfer coefficient from plastic to air is assumed to be η = 0.000 118cal/(cm2 K s) [35], if not otherwise stated. Discussing kinetics of the thermal and electronic contributions to the RIC, we have noted that the heat diffusion time of the core tcore ≈ 4r02 /a21 ≈ 16 µs is much less than the Yb-ion life time of the level 2F5/2 , τ1 = 830 µs and the time for which the general temperature balance between the 2 (1−p) fiber and its surrounding is achieved tbal ≈ r2 K2a ≈ 18 s (this relationship is determined by 2η 2

the decrement of the first dominating term in Eq. (9)). Figure 4 shows the temporal behavior of the RIC at 1064 nm on the fiber axis simulated for two cladding pump powers: 145 mW in Fig. 4(a), and 100 W in Fig. 4(b). The comparison of these Figs. highlights a strong increase of the total RIC with increase of the pump power due to income of both the thermal and electronic contributions. However, income of electronic RIC at growing pump power is less than the increment of the thermal component. At high pump power, the population becomes saturated but heat continues to rise. The electronic contribution dominates at the beginning, while the thermal RIC dominates in the steady-state. #193649 - $15.00 USD Received 11 Jul 2013; revised 31 Aug 2013; accepted 5 Sep 2013; published 16 Sep 2013 (C) 2013 OSA 23 September 2013 | Vol. 21, No. 19 | DOI:10.1364/OE.21.022374 | OPTICS EXPRESS 22383

Refractive index change, 10−6

Refractive index change, 10−6

(a) 145 mW pumping

2.7 2.4 2.1 1.8 1.5 1.2 0.9 0.6 0.3 0 0

1

2

3

4

(b) 100 W pumping

20 16 12 8 4 0

5

0

Time, s

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

Time, s

Figure 4. Time dependence of the refractive index change ∆n(t, z0 ) at 1550 nm for z0 = 20 cm in the case of cladding pumping of 145 mW(a) and 100 W (b). The total refractive index change is in black, the electronic refractive index change in orange and the thermal refractive index change in green. Curves are shown for input signal at 1060 nm with powers of 0 mW (solid curves), 10 mW (dashed curves), and 100 mW (dashed-dotted curves).

The increase of the input signal power increases heating and decreases of the population inversion. So, the thermal contribution increases whereas the electronic contribution decreases Figure 5 shows the time behaviors of all the contributions to the phase shift on the probe beam at 1550 nm for the case of cladding pumping [Fig. 5(a)] and core pumping [Fig. 5(b)], and without signal amplification. The orange curve is the electronic contribution to the phase shift, and in each case, after a rapid transient, this electronic phase shift saturates. Green, blue and red curves show respectively the thermal contributions to phase shift arising from the thermo-optic coefficient, the elongation of the fiber, and the lateral expansion of the fiber. As expected, the thermo-optic effect on the phase shift is much larger than the two other thermal effects. The comparative analysis of all the contribution to the phase shift on the probe beam at 1550 nm for the cladding or core pumping (in the absence of an amplified beam) highlights several stages in the RIC dynamics: domination of the electronic component at the beginning, alignment of the electronic and thermal contribution, and domination of the thermal RIC in steady-state. (a) Cladding pumping

103

Phase shift, radian

Phase shift, radian

103 101 10−1 10−3 10−5

10−7 10−5 10−4 10−3 10−2 10−1

Time, s

100

101

102

(b) Core pumping

101 10−1 10−3 10−5 10−7 10−5 10−4 10−3 10−2 10−1

100

101

102

Time, s

Figure 5. Time dependence of the phase shift of the probe beam at 1550 nm without signal at 1060 nm and for a pump power of 145 mW. Orange curve is the electronic contribution to the phase shift, green curve is the thermal contribution to the phase shift due to thermooptic coefficient, blue curve is the thermal contribution to the phase shift associated with the elongation of the fiber, red curve is the thermal contribution to the phase shift associated with the lateral expansion of the fiber, and dashed-black curve is the total phase shift. Axes are plotted on logarithmic scales


It should be noticed that the relative electronic contribution to the phase shift is less than the relative electronic contribution to RIC. This is explained by the fact that (1) the electronic phase shift is proportional to the integral of the product of the RIC and the mode profile of the probe beam, and (2) the thermal RIC is located on whole fiber structure whereas the electronic RIC is located only inside the core of the fiber. As the mode of the probe beam is present outside the core too, this leads to the decrease of the electronic contribution to the phase shift. The expression of the phase shift associated with the electronic part of RIC for the cases of core and clad pumping (without an amplified signal) can be found from Eq. (1) and with the steady-state solutions of (24) to (27): ˆ ˆ 4π FL2 ∆pξ τ1 4π 2FL2 ∆pξ l ∞ N2 (z, r)rdrdz = [Pp (0) − Pp(l)] δ ϕN = (36) 2 λ n0 r0 /2 0 0 λ n0 r02 hν p Score Theoretical estimate of the electronic part of the phase shift given by expression (36) coincides with the value obtained by numerical modeling. It is 0.75π for 145 mW cladding pumping (the pump intensity does not depend on the transverse coordinate) and is 2.87π for similar core pumping. The experimental value of the phase shift in the core pumped Yb-doped fiber was about 3.8π [4,17]. This difference can be explained [16] by taking into account the doping profile which is uniform in the simulation case and was gaussian-like in the experiment but with mean value given in table 1. The alignment time of the thermal and electronic phase shift for clad pumping with power of 145 mW is estimated to be tm = 839 ms from expression (31) and is close to the calculated value of 757 ms (the intersection point of the green and orange curves on Fig. 5(a)). 4.

Comparing of thermal and electronic contributions to the phase shift of a probing beam during pumping in the presence of the amplified beam

The presence of the amplified beam which depopulates the Yb3+ ions excited state can strongly change the ratio of RIC contributions. One can see from the Figs. 4(a)-4(b) that, with increase of the amplified beam power, the thermal contribution to RIC increases, while the electronic contribution decreases. The thermal loading increases due to growing rate of the cycle: the excited states populated by the pump are depopulated by the amplified beam. This process is more efficient for higher pump power levels causing stronger effect of the amplified beam on the thermal and electronic RICs. The total phase shift shown in Fig. 6 also depends on power of the pumping and amplified beams. The electronic component dominates over the thermal one during the transient stage. However, in steady-state, the thermally induced phase shift is higher than the electronic contribution. An increase of the amplified beam power leads to depopulation of the excited level of the Yb3+ ions and a decrease of the electronic contribution. The thermal component continues to increase due to higher rate of the excited state population - depopulation cycle. For a given set of parameters, the alignment time of the electronic and thermal components of the phase shift is less for cladding pumping than for the core pumping. The alignment time of the phase-shift components shown in Fig. 7 depends on the pump and amplified beam powers, the fiber length and the heat transfer coefficient η . The increase of the amplified beam power reduces the population inversion and increases the heat evolution, so the alignment time is reduced. Reduction of the pump power causes a decrease of both excited state population and amplified beam power, leading to a strong reduction of the thermal contribution and a less drastic decrease of the electronic contribution. Therefore, the alignment time of both contributions increases with the reduction of the pump power. Lower values of heat-transfer coefficient η provide higher thermal effect, and so the electronic contribution could no longer dominate over the thermal one, the alignment time #193649 - $15.00 USD Received 11 Jul 2013; revised 31 Aug 2013; accepted 5 Sep 2013; published 16 Sep 2013 (C) 2013 OSA 23 September 2013 | Vol. 21, No. 19 | DOI:10.1364/OE.21.022374 | OPTICS EXPRESS 22385

160 140

Phase shift, radian

120 100 80 60 40 20 0 0

0.1

0.2

0.3

0.4

0.5

Time, s

Figure 6. Time dependence of the phase shift at 1060 nm for cladding pumping of 100 W. The total phase shift is in black, the electronic phase shift in orange and the thermal phase shift in green. Curves are shown for input signal powers of 0 mW (solid curves), 10 mW (dashed curves), and 100 mW (dashed-dotted curves). 2.4

Alignement time, s

2.0 1.6 1.2 0.8 0.4 0 10−14

10−12

10−10

10−8

10−6

10−4

10−2

100

Relative input power of amplified beam

Figure 7. The dependence of the alignment time of total contribution to the the phase shift on the input signal power at 1060 nm for cladding pumping at 145 mW for the curves in orange and blue, and for cladding pumping at 100 W for the curves in red, green and black. The fiber length is equal to 2 m (orange, red and black curves), 20 cm (blue curve), and 10 m (green curve). The heat transfer coefficient η is equal to 0.000 118 cal/(cm2 K s) at curves blue, orange, green, and red, and 0.2388 cal/(cm2 K s) at curve black

increases. The estimation results are in good qualitative agreement with direct measurements of the phase shift induced in Yb-doped fiber by 1064 nm radiation [36]. 5.

The electronic and thermal phase shifts for the amplified pulse train

The repetitively-pulsing amplified signal combined with CW pumping is an operational regime commonly used with real amplification systems. The numerical simulation of the phase shift at the signal wavelength of 1060 nm highlights interesting dynamical behavior with antiphase behavior of the electronic and thermal phase shift contributions during the signal pulse propagation [Fig. 8]. The reduction of the excited state population during the amplified pulses causes negative RIC, while the thermal load induced by the stimulated transitions to the ground-state


leads to the increase of the refractive index. For amplified pulses with duration in nanoseconds - microseconds time scale, i.e. much shorter than the excited-state life time, the decrease of the electronic index change is stronger than the income of the thermal RIC. As a result the total index decreases during the amplified pulses, while the time-averaged refractive index increases due to increase of the thermal contribution. 60 Phase shift, radian

Phase shift, radian

60

40

20

0 0

10

20

30

40

50

40

20

0 200.99 201.01 201.03 201.05 201.07 201.09 201.11

Time, ms

Time, ms

Figure 8. Time dependence of the phase shift (a) in the amplified rectangular pulse train at 1060 nm, and (b) its zoom corresponding to 20th pulse. 10 W cladding pumping is used and the signal input power is 100 mW. The pulse train period is 10 ms, the pulse duration is 20 µs. Curve in black is the total phase shift, curve is green is the thermal component and curve in orange is the electronic component.

6.

Conclusion

We have developed a theoretical model for the description of thermal and electronic mechanisms of refractive index changes in Yb-doped optical fiber. The first or second effects can dominate depending on parameters of the system. At the initial stage of the pumping or pulse amplification, the electronic contribution dominates for any types of Yb-doped fibers. The thermal contribution becomes more significant at CW operation, specially when the amplitudes of the amplified and pumping beams (in saturation conditions) are increased, when fiber core diameter is increased, and when the heat transfer coefficient to air is decreased. The electronic index change predominates over the thermal contribution during amplified signal pulses with duration less than the excited-state life time even in the case of CW pumping. This work was supported by European "ERA-NET RUS INTENT", FEDER-Wallonia "Mediatic" and FP7 IRSES projects, the Interuniversity Attraction Pole program IAP PVII-35 of the Belgian Science Policy and by the Ministry of Education and Science of the Russian Federation and Russian Academy of Sciences. Appendix. Particular solutions of the heating problem 1. To give insight on the thermal behavior, one can simplify the system by considering a fiber without a plastic coating and assume simpler boundary conditions such as constant temperature at the cladding-air interface or no heat transfer from cladding to air (the fiber is perfectly insulated). In that way, we end up with a simple cylinder whose analytical solutions can be found in the literature [37]. We assume that the heat source Q(r,t) is uniformly distributed in the core (r ≤ r0 ) and has a rectangular time profile of duration τ p . When the temperature at the external boundary of the glass is fixed (ideal heat sink at the external boundary), the expression of the temperature on


the axis (r = 0) of the activated fiber could be evaluated from Eqs. (9) - (13) as: r ∞ J1 µn 0 r1 r0 Sn (t, τ p , τ1 ) δ T (t) = 2 τ1 Q ∑ 3 2 r1 µ J ( µ n) n=1 n 1

(37)

where τ1 = r12 /a21 , µn (parameter without dimension) is the n-th positive root of the Eq. J0 (µ ) = 0 and Sn (t, τ p , τ1 ) is given by: i h  1 − exp −µn2 τt if t ≤ τ p 1 i (38) Sn (t, τ p , τ1 ) = h  1 − exp −µ 2 τ p exp −µ 2 t−τ p if t ≥ τ p n τ1 n τ1

For a sufficient long time after the switch-off of the heat source, the increase of temperature δ T goes to zero. Another limiting case occurs when the heat flow at the external boundary of the glass is equal to zero (heat-insulated external glass boundary). The temperature on the axis of the activated fiber (r = 0) is expressed as: r ∞ J1 µn 0 r1 r2 r0 δ T (t) = 2 τ1 Q ∑ 3 2 (39) Sn (t, τ p , τ1 ) + 02 QU(t, τ p ) r1 r1 n=1 µn J0 (µn ) where µn (parameter without dimension) is n-th positive root of the Eq. J1 (µ ) = 0 and U(t, τ p ) is given by: ( t if t ≤ τ p U(t, τ p ) = (40) τ p if t ≥ τ p In this last case, the increase of temperature δ T has a permanent term after the switch-off of the heat source, and this term increases linearly with the heat duration τ p . 2. In general, the population and the heat source depend on the transverse coordinate r, the boundary condition at the external plastic-clad is arbitrary and the steady-state solution on the axis (r = 0,t → ∞) can be expressed from (9) as: ˆ ˆ ˆ ′ 1 r0 1 r ′′ K2 r2 K1 r1 1 r0 ′ ′ ′ + ln r Q(r )dr + 2 r Q(r′′ )dr′′ dr′ (41) δ Tst = + ln r1 K2 r0 a21 0 η r2 a1 0 r ′ 0 and if Q does not depend on the transverse coordinate (cladding pumping for instance), this solution can be simplified to: Qr2 K2 r2 K1 r1 1 + ln + (42) δ Tst = 20 + ln r1 K2 r0 2 η r2 2a1 For a rectangular heat source, the temperature on the core axis can be found from (9) as: ˆ r0 ∞ 1 K1 µn r ′ Q(r′ )r′ dr′ Sn (t, τ p , 1) δ T (t) = ∑ J0 (43) 2 2 a1 n=1 Zn a1 µn 0 where µn (parameter in s−1/2 ) is n-th positive root of the Eq. (11). In the case of the cladding pumping, the heat source Q can be assumed to independent of r, and the solution (43) simplified to: ∞ 1 K1 r0 µn r0 Sn (t, τ p , 1) δ T (t) = ∑ Q J1 (44) 3 a1 n=1 Zn a1 µn
