We show that polarization modulated and polarization multiplexed transmission may be significantly impaired by the polarization scattering induced by intra-channel cross-phase modulation and four-wave mixing. In polarization multiplexed transmission, channel interleaving may be used to mitigate the effect when two-pulse collisions are dominant.
© 2011 OSA
It is well known that intra-channel cross-phase modulation (XPM) may produce significant polarization scattering in soliton transmission . Intra-channel XPM is however a significant process in highly dispersed pulse transmission, so that we may expect a similar polarization scattering to occur in polarization multiplexed or polarization modulated highly dispersed pulse transmission even if the transmission is single channel. This effect may impair high baud rate single channel transmission like it impaired multi-channel soliton transmission. Using a perturbation approach , we analyze in this Letter the effect of intra-channel collisions in polarization modulated or polarization multiplexed systems, and show that they may indeed induce signal depolarization similar to that caused by collisions between wavelength division multiplexed solitons . Simple analytical expressions capturing the essence of the phenomenon and quantitatively predicting the observed amount of depolarization are derived. We finally show that a relative delay of the two multiplexed channels by half of a symbol duration (channel interleaving) reduces the scattering by more than one order of magnitude.
2. First order perturbation theory
For moderate amounts of polarization mode dispersion (PMD), the evolution of the vector electric field E of a polarized optical field in a nonlinear dispersive fiber, averaged over the fast polarization evolution, is well described by the Manakov equation [1, 3]Eq. (2) is extended to terms where the four-wave mixing interaction is resonant with the pulse E0, which are those with Tk = Tj + Tl hence k = j + l. In a three dimensional space (the Stokes space) spanned by the unit vectors (e⃗1, e⃗2, e⃗3), let us now introduce the Pauli spin vector as σ⃗ = (|x〉〈x| – |y〉〈y|)e⃗1 + (|y〉〈x| + |x〉〈y|)e⃗2 + i(|y〉〈x| – |x〉〈y|)e⃗3, where |x〉 and |y〉 are unit polarization vectors parallel to the x and y axis respectively. Let us then define En = En|s〉 where |s〉 is the unit polarization state, and the Stokes vectors of the field is s⃗ = 〈s|σ⃗|s〉. After multiplying by both sides of Eq. (2), adding the hermitian conjugate, and integrating over t we obtain Equation (5) may be solved using a perturbation approach, by entering at right hand side for En, n = 0, k,l, j and |sn〉, their expressions in the absence of nonlinearity. To this purpose, let us assume that the input pulses have a well-defined polarization, uniform across the pulses, that is En(z = 0,t) = En(0,t)|s〉, where the unit polarization vector |s〉 is independent of t. A purely dispersive propagation preserves the uniformity of the polarization hence, within the perturbation approach, 〈s0|σ⃗|sl〉〈sj+l|sj〉 at right hand side of Eq. (5) may be assumed as time independent quantities, and pulled out of the integral.
It is interesting to consider the limit case of two-pulse interaction, in which two pulses of the four at right hand side of Eq. (5) overlap, that is the particular case j = 0. Using that and i[σ⃗(σ⃗ ·s⃗l) – (σ⃗·s⃗l)σ⃗]/2 = s⃗l × σ⃗ , we obtain from Eq. (5)(5) and (6) become
The case of two-pulse interaction described by Eqs. (6) and (8) corresponds to pure XPM. In this case, the component parallel to s⃗0 at the right hand side of Eqs. (6) and (8) is zero, hence dlnU0/dz = 0. This means that two-pulse collisions does not produce exchange of energy between the two pulses but only a change in their polarization. In addition, there is no polarization scattering if the pulses are either co-polarized (s⃗0 parallel to s⃗l) or cross-polarized (s⃗0 anti-parallel to s⃗l). Note once again the strong analogy between the two-pulse case and the polarization scattering caused by inter-channel XPM in wavelength division multiplexed and polarization multiplexed soliton transmission .
Four pulse interactions described by Eqs. (5) and (7) with j ≠ 0, on the other hand, correspond to non-degenerate intra-channel four-wave mixing. In this case, intra-pulse power transfer may take place. This is reflected by the fact that 〈s0|σ⃗|sl〉 may have a component parallel to s⃗0. Furthermore, polarization scattering may also occur between three co-polarized pulses and one cross-polarized. To see this, let us assume without loss of generality that s⃗0 = e⃗2. In this case, in Eqs. (5) or (7) the rate of power change dlnU0(z)/dz is proportional to |〈s0|σ⃗|sl〉 ·e⃗2|2 = (1 +s⃗l ·s⃗0)/2, whereas the rate of depolarization ds⃗0/dz have components along e⃗1 and e⃗3 proportional to |〈s0|σ⃗|sl〉 ·e⃗1,3|2 = (1 – s⃗l · s⃗0)/2. Therefore, when s⃗l ·s⃗0 = −1, hence when the pulses at l and 0 are orthogonally polarized, ds⃗0/dz may have non-zero depolarizing components parallel to e⃗1 and e⃗3 if the pulse at j is co-polarized with the pulse at j + l so that |〈sj+l|sj〉| = 1.
3. Numerical example: polarization multiplexed QPSK signal
Let us now evaluate the variance of the polarization spread in a polarization-multiplexed quadrature phase shift keying (QPSK) signal, where all pulses have the same amplitude Ak = A and the phases of the pulses on each polarization are 0, π/2, π and 3π/2, assuming synchronous pulse streams on the two polarizations. To evaluate the effect of two-pulse collisions, we notice that in this case the possible pulse polarizations are the two linear polarizations s⃗ = (0,±1,0) and the two circular polarizations s⃗ = (0,0,±1). Let us assume that the pulse at T0 = 0 has the linear polarization s⃗0 = (0, 1,0). The interacting pulses will be with equal probability linearly polarized, in which case s⃗0(0) × s⃗l(0) = 0, or circularly polarized, in which case s⃗0(0) × s⃗l(0) = (±1,0,0). Assuming a random message, the polarization spread Δs⃗0 =s⃗0(z) – s⃗0(0) is on average zero, 〈Δs⃗0〉 = 0, whereas its variance along e⃗1 isEq. (10).
The estimate of the polarization scattering caused by multiple pulse collisions requires the knowledge of the polarization and the absolute phase of the 16 transmitted pulses. They are in our example ±|±s2〉, ±|s±〉, ±i| ±s2〉, and ±i|s±〉, where , . Assume that |s0〉 = | + s2〉. We have 〈s2|σ⃗|s2〉 = e⃗2, 〈s2|σ⃗| – s2〉 = e⃗1 + ie⃗3, 〈s2|σ⃗|s+〉 = (1 – i)e⃗1/2+ (1+ i)(e⃗2 + e⃗3)/2, 〈s2|σ⃗|s−〉 = (1+ i)e⃗1/2+ (1 – i)(e⃗2 – e⃗3)/2, and 〈s2|s±〉 = (1 ± i)/2 〈 − s2|s±〉 = (1 ∓ i)/2, 〈s2| – s2〉 = 〈s+|s−〉 = 0.
When the interaction is with a pulse l circularly polarized or co-polarized with the pulse in 0, the four wave mixing interaction has a component parallel to s⃗0, which produces a change of the pulse energy U0. This means that only the intra-channel two-pulse interaction is a real soliton-like collision, where there is no power transfer between the two pulses but only a polarization change. Four-pulse interaction (not unexpectedly) produces a polarization rotation of the pulse at 0 together with a phase dependent and polarization dependent power variation.
Notice that the average of 〈sj+l|sj〉 over all pulse is zero, so that the average polarization displacement is also zero. Squaring the expression for Δs⃗0(z) · e⃗1,3, we obtain a four-fold sum of terms proportional to product of the real part of a1 = 〈s0|σ⃗|sl〉〈sj+l|sj〉 and a2 = 〈s0|σ⃗|sl′〉〈sj′+l′|s′j〉. Let us now use the expansion . The first term is zero for the phase symmetry of the channel. The second term is the real part of 〈s0|σ⃗|sl〉〈sj+l|sj〉 · 〈sj′|sj′+l′〉〈sl′|σ⃗|s0〉, which for isotropic polarization distribution has non-zero average only if j = j′ and l = l′. Using now that 〈|sj+l|sj〉|2〉 = 1/2 and 〈|〈s0|σ1,3|sl〉|2〉 = 1/2, and the statistical independence of the two terms, we obtain 〈Re[2iG〈s0|σ1,3|sl〉〈sj+l|sj〉] Re [2iG〈s0|σ1,3|sl′〉〈sj′+l′|sj′〉]〉 = δj,j′δl,l′|G|2/2. Notice that, being the polarization of the pulse at 0 parallel to e⃗2, depolarization of this pulse is caused by the fluctuations of the first and third component of the Stokes vector only.
The contribution of the four pulse collision to the polarization spread along e⃗1 and e⃗3 is the sum of the average squared of the individual elements of the sum, that isEq. (7).
In Fig. 1 we show the depolarization of the pulse Stokes vectors on the Poincaré sphere in a transmission systems employing two polarization multiplexed QPSK signals at 50 Gbaud/s (200 Gbit/s aggregate rate). The pulses have a 5 ps full-width at half maximum Gaussian profile.
The system is made of 7 spans of zs = 100 km each hence total length L = 700 km, with precompensation −β″z*, full compensation after each span, and post-compensation at the receiver −β″(zs – z*) so that 100% of the total fiber dispersion is compensated. Figure 2 has been obtained with no pre-compensation, z* = 0. The fiber dispersion is 16 ps/nm/km, loss is 0.25 dB/km, and nonlinear coefficient γ = 1.3 W−1 km−1. The average input power is (left to right) 3, 7 and 13 dBm. The polarization rotation caused by intra-channel collisions is evident.
Let us now make a quantitative comparison between theory and simulations, for the same system analyzed in Fig. 1. In Fig. 2a we show the standard deviation of the first component of the Stokes vector, 〈[Δs0(L) · e⃗1]2〉1/2, for 7 dBm input power vs. the pre-compensated length of fiber. The red dot-dashed line is the contribution of the two-pulse collisions alone, the purple dashes line is the contribution of the four pulse collisions alone, the blue solid line is the combination of the two, so it is the observed value of the standard deviation of the polarization jitter along e⃗1. Along e⃗3, the contribution of the two-pulse collisions to polarization jitter is zero, and only the contribution of the four pulse collisions is effective, and its variance is equal to the variance of the four pulse collisions along e⃗1. Consequently, the purple dashed curve represents the observed standard deviation of the fluctuations of the third component of s⃗0, that is 〈[ Δs⃗0(L) · e⃗3]2〉1/2. The same expression may be obtained for the standard deviation of the term parallel to e⃗2, hence being s⃗0 parallel to e⃗2, to the fluctuations of the pulse amplitude [the second term at left hand side of Eq. (7)]. Scattered symbols represent the results directly obtained from simulations, with the exception of the symbols over the two-pulse collision curve, which have been inferred using . The imperfect agreement at low pre-compensation is probably caused by the small depolarization of the pulses induced by the nonlinear interaction, more effective for small pulse pre-dispersion. The simulations were performed using the Manakov equation, yet the full equations including random mode coupling gave in our numerical example identical results if the fiber PMD did not exceed 0.05 ps/km1/2, and results within 10% if the fiber PMD did not exceed 0.1 ps/km1/2.
A possible way to significantly reduce the polarization scattering is polarization interleaving. If the two orthogonally polarized channels are offset by half of a symbol time, the interaction is between pulses either parallel or orthogonally polarized and two-pulse collision induced polarization rotation is to first order absent. This approach is expected to be effective only when two-pulse scattering is dominant over multiple-pulse scattering, that is for moderate amounts of pre-compensation, because four-pulse interactions produce polarization scattering also when all pulses are either co- or cross-polarized. To verify the effect of channel interleaving, we show in Fig. 2b by solid lines and in a logarithmic scale the standard deviation of the fluctuations of the first component of the Stokes vector, 〈[Δs⃗0(L) · e⃗1]2〉2–pulse vs. the dispersion pre-compensation z* obtained using Eq. (10), with no interleaving. Points show the results of simulations. The three curves refer to an average input power of 3 (lower solid curve, triangles), 7 (intermediate solid curve, squares) and 13 (upper solid curve, circles) dBm. In the same figures, dashed black lines connect triangles, squares and circles representing the result of simulations with pulse interleaving. Again, triangles, squares and circles refer to 3, 7 and 13 dBm input power. The curves show a reduction of the polarization jitter by more than one order of magnitude, down to a level compatible with that expected from the sole contribution of four-pulse collisions. The reduction of depolarization obtained with pulse interleaving relies on the validity of the Manakov equation, which neglects PMD, hence the benefit of this technique may become less significant as PMD increases.
We have shown that intra-channel collisions may introduce significant polarization scattering in systems in which non orthogonal polarizations are transmitted. This phenomenon is the intra-channel analogous of the inter-channel polarization scattering in wavelength division multiplexed transmission. A first order perturbation theory including both two-pulse collision and four pulse collisions and quantitatively predicting the expected amount of depolarization has been presented. Finally, we have predicted and verified by computer simulations that channel interleaving reduces polarization jitter by more than one order of magnitude compared to the case of synchronous channels.
The work at Fondazione Ugo Bordoni was carried out in the framework of the Project “ATENA” supported by the Italian Ministry of Economic Development.
References and links
2. A. Mecozzi, M. Tabacchiera, F. Matera, and M. Settembre, “Intra-channel nonlinearity in differentially phase-modulated transmission,” Opt. Expr. 19, 3990–3995 (2011). [CrossRef]
3. D. Marcuse, C. R. Menyuk, and P. K. A. Wai“Application of the Manakov-PMD equation to studies of signal propagation in optical fibers with randomly varying birefringence,” J. Lightwave Technol. 15, 1735–1746 (1997). [CrossRef]