Activated polarization pulling and de-correlation of signal and pump states of polarization in a fiber Raman amplifier

We report on a theoretical study of activated polarization pulling and de-correlation of signal and pump states of polarization based on an advanced vector model of a fiber Raman amplifier accounting for random birefringence and two-scale fiber spinning. As a result, we have found that it is possible to provide de-correlation and simultaneously suppress PDG and PMD to 1.2 dB and 0.035 ps/km 1/2 respectively. ©2011 Optical Society of America OCIS codes: (060.2320) Fiber optics amplifiers and oscillators; (060.4370) Nonlinear optics, fibers; (190.5650) Raman effect. References and links 1. C. D. Poole and R. E. Wagner, “Phenomenological approach to polarization dispersion in long single-mode fibers,” Electron. Lett. 22(19), 1029–1030 (1986). 2. P. K. A. Wai and C. R. Menyuk, “Polarization mode dispersion, decorrelation and diffusion in optical fibers with randomly varying birefringence,” J. Lightwave Technol. 14(2), 148–157 (1996). 3. R. H. Stolen, “Polarization effects in fiber Raman and Brillouin lasers,” IEEE J. Quantum Electron. 15(10), 1157–1160 (1979). 4. Q. Lin and G. P. Agrawal, “Vector theory of stimulated Raman scattering and its application to fiber-based Raman amplifiers,” J. Opt. Soc. Am. B 20(8), 1616–1631 (2003). 5. S. Sergeyev, S. Popov, and A. T. Friberg, “Spun fiber Raman amplifiers with reduced polarization impairments,” Opt. Express 16(19), 14380–14389 (2008). 6. S. Sergeyev, S. Popov, and A. T. Friberg, “Virtually isotropic transmission media with fiber Raman amplifier,” IEEE J. Quantum Electron. 46(10), 1492–1497 (2010). 7. S. Sergeyev and S. Popov, “Two-section fiber optic Raman polarizer for high-speed transmission systems,” in the 13th International Conference on Transparent Optical Networks (June 26–30, 2011) Stockholm, Sweden, Th.A6.7. 8. M. J. Li and D. A. Nolan, “Fiber spin-profile designs for producing fibers with low polarization mode dispersion,” Opt. Lett. 23(21), 1659–1661 (1998). 9. A. Galtarossa, L. Palmieri, A. Pizzinat, B. S. Marks, and C. R. Menyuk, “An analytical formula for the mean differential group delay of randomly-birefringent spun fibers,” J. Lightwave Technol. 21(7), 1635–1643 (2003). 10. M. Martinelli, M. Cirigliano, M. Ferrario, L. Marazzi, and P. Martelli, “Evidence of Raman-induced polarization pulling,” Opt. Express 17(2), 947–955 (2009). 11. V. V. Kozlov, J. Nuño, J. D. Ania-Castañón, and S. Wabnitz, “Theory of fiber optic Raman polarizers,” Opt. Lett. 35(23), 3970–3972 (2010). 12. L. Ursini, M. Santagiustina, and L. Palmieri, “Raman nonlinear polarization pulling in the pump depleted regime in randomly birefringent fibers,” IEEE Photon. Technol. Lett. 23(4), 254–256 (2011). 13. N. J. Muga, M. F. S. Ferreira, and A. N. Pinto, “Broadband polarization pulling using Raman amplification,” Opt. Express 19(19), 18707–18712 (2011). 14. J. E. Heebner, R. S. Bennink, R. W. Boyd, and R. A. Fisher, “Conversion of unpolarized light to polarized light with greater than 50% efficiency by photorefractive two-beam coupling,” Opt. Lett. 25(4), 257–259 (2000). 15. J. Fatome, S. Pitois, P. Morin, and G. Millot, “Observation of light-by-light polarization control and stabilization in optical fibre for telecommunication applications,” Opt. Express 18(15), 15311–15317 (2010). 16. A. Zadok, E. Zilka, A. Eyal, L. Thévenaz, and M. Tur, “Vector analysis of stimulated Brillouin scattering amplification in standard single-mode fibers,” Opt. Express 16(26), 21692–21707 (2008). 17. F. Marino, M. Giudici, S. Barland, and S. Balle, “Experimental evidence of stochastic resonance in an excitable optical system,” Phys. Rev. Lett. 88(4), 040601 (2002). 18. B. Lindner, J. García-Ojalvo, A. Neimand, and L. Schimansky-Geier, “Effects of noise in excitable systems,” Phys. Rep. 392(6), 321–424 (2004). 19. M. I. Dykman, B. Golding, L. I. McCann, V. N. Smelyanskiy, D. G. Luchinsky, R. Mannella, and P. V. E. McClintock, “Activated escape of periodically driven systems,” Chaos 11(3), 587–594 (2001). #153455 $15.00 USD Received 31 Aug 2011; revised 21 Sep 2011; accepted 23 Sep 2011; published 14 Nov 2011 (C) 2011 OSA 21 November 2011 / Vol. 19, No. 24 / OPTICS EXPRESS 24268 20. J. D. Ania-Castañón, V. Karalekas, P. Harper, and S. K. Turitsyn, “Simultaneous spatial and spectral transparency in ultralong fiber lasers,” Phys. Rev. Lett. 101(12), 123903 (2008). 21. S. A. Babin, V. Karalekas, E. V. Podivilov, V. K. Mezentsev, P. Harper, J. D. Ania-Castanon, and S. K. Turitsyn, “Turbulent broadening of optical spectra in ultralong Raman fiber lasers,” Phys. Rev. A 77(3), 033803 (2008).


Introduction
Polarization impairments, viz.polarization mode dispersion (PMD) and polarization dependent gain (PDG), are among major factors limiting the progress in further increase of transmission rates and overall capacity of the next generation of optical networks based on distributed fiber Raman amplification.PMD leads to the pulse broadening caused by varying group velocities for the pulses with different states of polarization (SOPs) [1,2] whereas PDG leads to an uncontrollable output gain variation caused by its dependence on the input signal SOP [3][4][5][6][7].By the traditional approach, spinning the fiber periodically, it is possible to reduce PMD to below 0.04 ps/km 1/2 [8,9], but this is accompanied with a simultaneous increase in Raman PDG [3][4][5][6][7].Sergeyev et al. have shown that it is possible to mitigate both PDG and PMD simultaneously by using a two-section fiber ('two-section approach') in which the first section has no spin and the second one is periodically spun [5][6][7].As follows from [7], this approach demonstrated polarization pulling of the signal SOP to the pump SOP and so can be considered as a two-section Raman polarizer which is the one of the techniques of all-optical polarization control [10][11][12][13][14][15][16].Unlike the one-section fiber Raman polarizer demonstrated by Martinelli et al. [10], Kozlov et al. [11], L. Ursini et al. [12] and N. J. Muga et al. [13], a twosection fiber provides almost the same polarization pulling for all input signal states of polarization and so results in very small polarization dependent gain [7].In view of increasing interest in polarization multiplexing of signals along with Raman amplification in fiber optic communication, polarization pulling approaches cannot be used straightforwardly and so have to be modified.In this paper, we have suggested a new approach based on application of fiber spinning with two different scales, viz.fast spinning with a period of 2.7 m and slow spinning with a period of more than 200 m for 10 km of single mode fiber.In view of the similarity of a two-scale periodically spun fiber Raman amplifier (SFRA) to a periodically driven excitable system (ES), we have found resonance-like de-correlation of pump and signal SOPs and polarization pulling equivalent to the Stochastic Resonance (SR) and Dynamic Localization (DL) in ES [17][18][19].As a result, we have demonstrated that it is possible to provide decorrelation of pump and signal SOPs and simultaneously suppress PDG and PMD to the 1.2 dB and 0.035 ps/km 1/2 respectively.The practical implementation of the obtained results can be in the design of a fiber Raman amplifier-based quasi-isotropic spatially and spectrally transparent media for high-speed fiber optic communication [20,21].

Model of a fiber Raman amplifier with random birefringence and arbitrary spin profile
The evolution of signal and pump states of polarization in a single mode fiber is typically described in terms of the unit vectors ˆˆˆ( , , ) p p p = p , respectively, pointing to positions on the Poincaré sphere (Fig. 1).Due to birefringence, these vectors rotate on the Poincaré sphere around the birefringence vector W i = (2b i cosθ, 2b i sinθ,0) T in the same direction, but at different rates b s and b p which are birefringence strengths (b i = π/L bi where L bi is the beat length) at signal λ s and pump λ p wavelengths.We assume that the birefringence strength 2b i is fixed and the orientation angle θ is driven by a white-noise process (fixedmodulus model [2] where <…> means averaging over the birefringence fluctuations along the fiber, δ(z) is a Dirac delta-function, and σ 2 = 2/L c (L c is the birefringence correlation length).
We neglect herein the fiber twist and, therefore, the birefringence vector for the spun fiber takes the form of W i = R 3 [2A(z)]W i,un , where A(z) is the spin profile, and R 3 (γ) represents rotation in the equatorial plane by angle γ around the z-axis on the Poincaré sphere [9] ( ) Stimulated Raman scattering leads to amplification of the signal wave and changes its direction and so evolution of the signal S and pump P SOPs can be found as follows: Here s 0 is the part of the signal amplitude related to the pump-signal SOPs interaction, ( ) where <…> means averaging over the birefringence fluctuations along the fiber.If the input pump and signal SOPs are parallel Raman gain takes the maximum value and if the SOPs are orthogonal than Raman gain takes the minimum one [4][5][6][7].The difference in gains (polarization dependent gain, PDG) can be found as follows [4][5][6][7] ( ) ( ) ( ) To characterize de-correlation of pump and signal SOPs and polarization pulling in terms of the Stochastic Resonance phenomenon we introduce parameter R which is similar to the signal-to-noise ratio in the models of SR [18]: To calculate s 0 (L), we use a vector model of a fiber Raman amplifier for forward pump with allowance for fiber spin profile and random birefringence and neglecting pump depletion [4]: Here 0 ŝ = s s , and W p (NL) , W s (NL) describes the nonlinear SOP evolution caused by self-and cross-phase modulation (SPM and XPM): , where n 2 is the nonlinear Kerr coefficient and A eff is the effective core area of the fiber.We choose the reference frame in the Stokes space in such a way that the local birefringence vector is ( ) , 2 , 0, 0 , i.e. oriented along the X-axis on the Poincaré sphere.This can be accomplished by a suitable transformation ( ) , and ( ) ( ) , provided the variable 0 s and scalar product x = ⋅ p s are invariant under the rotation.In addition, we consider cases which correspond to the maximum and minimum of Raman gain, viz.
( ) . Applying an averaging procedure considered in [5,6] to Eqs. (7), we find the system of equations: Here In view of ε 3 << ε 4 , ε 5 (rotation of the signal SOP with respect to the pump SOP is much slower than the SOPs rotation with respect the local birefringence vector i W ɶ ), we can average over the fast rotations.It can be done by the transformations ( ) Using this transformation, Eq. ( 7) and averaging over the fast oscillations, we find the following equation for 1 1 h f ( ) If we choose parameters for Raman amplifier as, g = 2.3 dB W −1 km −1 , P in = 5W, L = 10 km, L c = 110 m, we find that 1 0.15 1 c L L ε = << .Thus, we can neglect the first term on the right-hand side of Eq. ( 10).As a result, ( ) ( ) ( ) In addition, this procedure leads to and so self-and cross-phase modulation (SPM and XPM) can be neglected [4][5][6][7].A similar result has been obtained in [10,12] by direct modeling of stochastic Eqs. ( 7), viz. it was obtained that for the pump power P in < 10 W and PMD parameter D p >0.01 ps/km 1/2 SPM and XPM has no contribution to polarization pulling.Using transformation (9) and averaging over the birefringence fluctuations and fast oscillations we find the following equations which we use to calculate parameter R from Eq. ( 6): ( To calculate the spin induced reduction factor (SIRF) for the case of spun fiber we use the standard model of PMD [9] ( ) , , where ( ) ( ) are the mean-square differential group delays (DGD) for two orthogonal SOPs in the case of long-length spun fiber and the same fiber without spin, respectively [9].As follows from [5,6,9], after averaging over fluctuations caused by random birefringence, equations for SIRF take the following form: Finally, the PMD parameter for the spun fiber D ps can be found as follows [9] 2 .
In addition, we use fast periodic spinning along with slow-amplitude and phase modulation: .
Here A 0 is amplitude of the fiber spinning in rad, k lf and k hf are low and high frequencies of the fiber spinning.Practical realization of amplitude-and phase-modulated fiber spinning can be achieved based on the modern technology of direct fiber spinning instead of spinning the preform [8].
To quantify polarization pulling and SOPs de-correlation in terms of SR phenomena, we calculate parameters equivalent to the parameters used in the models of excitable systems [17][18][19], viz.Kramers length <L k > and intrawell relaxation length <L R >.In the absence of low-frequency modulation (k lf = 0), Eqs. ( 8) can be simplified for the case of k hf >>L c , L b with the help of averaging over fast periodic fiber spinning [9].As a result, we have Here ( ) ( ) . Linear stability analysis of Eqs.
(17) near the states 0 0 ˆ, x y results in the following eigenvalues: ( ) Thus, relaxation length can be defined as ( ) Im 0 Λ ≠ then the system oscillates around the states 0 0 ˆ, x y and so the averaged length at which the system is close to the state can be defined as ( )

Results and discussion
We have found 0 0 ˆ, x y , Re{Λ 1,2 }, Im{Λ 1,2 } numerically from Eqs. ( 18), (19).We find also PDG, R and PMD parameter D p numerically from Eqs. ( 5), ( 6), ( 8) and ( 12)-( 15) for twoscale fiber spinning (16).We use parameters typical for a distributed fiber Raman amplifier: L = 10 km, α s = 0.2 dB/km, λ p = 1460 nm, λ s = 1550 nm, g = 2.3 dBW As a result of nonlinearity in the form of Raman amplification and presence of random birefringence and periodic fiber spinning, PDG has maxima and minima as a function of correlation length L c and spinning frequency k lf .For slow amplitude modulation of the fiber spinning (Fig. 3 (a)), the maximum PDG corresponds to the case when L c , k lf →0.As follows from the results of [8,9], fast periodic fiber spinning results in PMD mitigation and so leads to polarization pulling [7][8][9][10][11][12][13].However, with increased frequency of amplitude modulation PDG reaches a minimum of 1.5 dB for L c = 105m and k lf = 6 km −1 and for L c = 55 m and k lf = 30 km −1 PDG reaches the maximum of 6.9 dB (Fig. 3 (a)).In the case of slow phase modulation, PDG reaches a minimum of 1.2 dB for L c = 110 m and k lf = 3 km −1 and a maximum of 14.5 dB for L c = 110 m and k lf = 22 km −1 (Fig. 3 (b)).
The evolution of angle Φ between pump and signal SOPs for the cases of minima and maxima is shown in Figs. 4. As follows from Figs. 3 (a, b) and Figs. 4, the maximum of PDG corresponds to the polarization pulling with cos 1 Φ → , and the minimum of PDG corresponds to the case of SOPs de-correlation for which pump and signal SOPs are freely rotating and so cos 0 Φ → .In addition, fiber spinning with low-frequency phase modulation provides better polarization pulling as compared to the case of amplitude modulation (thick lines in Figs. 4).(b, d)).

R, dB
To characterize polarization pulling and SOPs de-correlation in terms of SR phenomena, we present in Fig. 5 (a-d) the results for PDG, parameter R and PMD parameter D p as a function of modulation frequency k l,f and correlation length L c .For low-frequency amplitude modulation, PDG has a minimum and parameter R has a maximum at frequency of k lf,max = 6 km −1 which coincides with the resonance frequency for excitable systems calculated from Eq. (19).A maximum and minimum at this frequency is evidence of Stochastic Resonance where random birefringence fluctuations are synchronized with an external modulation [18].In view of the maximum in R coinciding with the minimum in PDG, SR for this case corresponds to de-correlation of the pump and input SOPs.For the off-resonance conditions (low and high frequencies of spinning), the random birefringence fluctuations have no synchronization with the external modulation and so there is no activation of de-correlation of pump and signal SOPs.For the case of phase modulated fiber spinning, PDG and parameter R have maxima and minima at frequencies different from SR frequencies (thin lines in Figs. 5 (a, c)).It is quite similar to the case of Dynamic Localization (DL) considered by Dykman et al., where the right and the left escape rates become different to each other depending on the phase of non-adiabatic aperiodic modulation [19].Thus, DL phenomena take the form of either polarization pulling or de-correlated SOPs depending on the phase of phase modulated fiber spinning.It is likely that increased noise power σ 2 (σ 2 = 1/L c ) leads to increased asymmetry of probabilities and so to the further increased PDG and parameter R (solid lines in Figs. 5 (a,  c)).
PDG has a minimum as a function of noise power σ 2 = 1/L c for amplitude and phase modulation (Fig. 5 (b)).When the birefringence fluctuations are large (small correlation lengths) and low (large correlation lengths) polarization pulling is dominating and PDG is high (Fig. 5 (b)).If fluctuations are large, escape events happen at very small intervals and so at larger intervals the fiber looks isotropic with small changes in pump and signal SOPs along the fiber.The case of small birefringence fluctuations corresponds to the polarization maintaining (PM) fiber where pump and signal SOPs initially oriented along the fast or slow axes preserve their orientations along the fiber.For an optimal level of noise above the activation threshold, random hops are synchronized with modulation frequency and so SR takes the form of de-correlation of SOPs (minimum in PDG) [17,18].As follows from Fig. 3(a), PDG is almost independent of modulation frequency k lf for amplitude modulation in the limit of high frequency modulation.It is likely that if the fiber spinning oscillates so fast that the system feels mainly the averaged spinning which can be small and not enough to activate de-correlation of SOPs [18].As a result, with increasing modulation frequency polarization pulling occurs instead of SOPs de-correlation (thick dotted line in Fig. 5 (b)).As mentioned before, low-frequency phase-modulation of fiber spinning can result in a DL phenomenon in the form of either polarization pulling or de-correlated SOPs depending on the phase.This results herein in a minimum in PDG as a function of correlation length (thin lines in Figs. 5 (b)).Parameter R can have maxima and minima as a function of a noise power σ 2 = 1/L c and modulation frequency k lf and so it is likely that our case corresponds to more complex cases of SR and DL.PMD parameter has maxima and minima as function of k l,f and L c −1 which can be an indication of SR and DL presence for polarization properties of fiber described by Eqs.(14).Thus, application of fast (2.7 m period) and slow (200-600 m period) amplitude or phasemodulated fiber spinning can result in activated polarization pulling or de-correlation of pump and signal SOPs.De-correlation is accompanied with simultaneous suppression of PDG and PMD to the 1.2 dB and 0.035 ps/km 1/2 respectively.

Conclusion
In conclusion, we have developed an advanced vector model of a fiber Raman amplifier accounting for fiber spinning and random birefringence.Based on this model and the PMD model previously developed in [9], we characterize a Raman amplifier in terms of PDG, gain fluctuations and PMD.Unlike previous approaches, we demonstrated that it is possible simultaneously to de-correlate pump and signal SOPs and to suppress PDG and PMD.This new technique is based on application of fiber spinning with two different scales: fast spinning with period of 2.7 m and slow spinning with period of 200 -6000 m for 10 km of single mode fiber.In view of the similarity of a two-scale periodically spun fiber Raman amplifier to a periodically driven excitable system, we have found resonance-like decorrelation of pump and signal SOPs and polarization pulling equivalent to the Stochastic Resonance and Dynamic Localization in excitable systems [17][18][19].As a result, we demonstrate that it is possible to provide de-correlation of pump and signal SOPs and simultaneously suppress PDG and PMD to the 1.2 dB and 0.035 ps/km 1/2 respectively.The practical implementation of the obtained results can be in the design of fiber Raman amplifierbased quasi-isotropic spatially and spectrally transparent media for high-speed fiber optic communication [20,21].

Fig. 1 .
Fig.1.Evolution of the pump p and the signal ŝ states of polarization on the Poincaré sphere, as well as the fluctuations of the local birefringence vector Wi = (2bicosθ, 2bisinθ,0) T .Vectors ŝ and p rotate around the local axis (W) at rates bp and bs, while vector (W) rotates randomly in the equatorial plane at the rate σ = Lc -1/2 (Lc is the correlation length), Φ is an angle between pump and signal SOPs.