# Pulse Propagation in Presence of Polarization Mode Dispersion and Chromatic Dispersion in Single Mode Fibers

Authors: Hassan Abid Yasser; Nizar Salim Shnan
Journal: International Journal of Optics
 (2013)

---

## Abstract

The presence of (first and second orders) polarization mode dispersion (PMD), chromatic dispersion, and initial chirp makes effects on the propagated pulses in single mode fiber. Nowadays, there is not an accurate mathematical formula that describes the pulse shape in the presence of these effects. In this work, a theoretical study is introduced to derive a generalized formula. This formula is exactly approached to mathematical relations used in their special cases. The presence of second-order PMD (SOPMD) will not affect the orthogonality property between the principal states of polarization. The simulation results explain that the interaction of the SOPMD components with the conventional effects (chromatic dispersion and chirp) will cause a broadening/narrowing and shape distortion. This changes depend on the specified values of SOPMD components as well as the present conventional parameters.


---

## Body

## 1. Introduction

As a pulse propagates through a light-wave transmission system with a polarization mode dispersion (PMD), the pulse is spilt into a fast and slow one and therefore becomes broadened. This kind of PMD is commonly known as first-order PMD (FOPMD). Under FOPMD, a pulse at the input of a fiber can be decomposed into two pulses with orthogonal states of polarization (SOP). Both pulses will arrive at the output of the fiber undistorted and polarized along different SOPs, the output SOPs being orthogonal [1]. Both the PSPs and the differential group delay (DGD) are assumed to be frequency independent when only first-order PMD is being considered [2].PMD also has higher-order components. Second-order PMD (SOPMD) is formed due to dependence of DGD on wavelength, because of that PMD is random variable. Statistical nature of PMD makes compensation very difficult [3]. Statistical properties of first- and second-order PMD have been reported in [4–6]. In these approaches, mean DGD plays a central role. Further, analytical calculations of pulse broadening and PMD compensation have also been reported [7, 8]. However, these analytical results make use of autocorrelation of the PMD vector [4].Also, it is shown that FOPMD and SOPMD vectors are statistically dependent and tend to be mutually perpendicular. It is important to calculate the probability that a certain value of SOPMD happens, when the value of DGD is known in a transmission line [8, 9]. The frequency dependence of the DGD causes polarization-dependent chromatic dispersion (PCD), resulting in polarization-dependent pulse compression and broadening. While PCD is a simple mechanism, it is a relatively minor component of the SOPMD vector since the frequency dependence of PMD vector has a statistical tendency to point away from the PSP. The dominant impairment mechanism is due to the depolarization, which causes distortion of the pulse shape. Statistical analyses and numerical simulations indicate a strong interaction between CD and the SOPMD in general [10–13].Chirp in the transmitted signal has a significant impact on second-order impairments. Impairments were found to be dependent on chirp and CD and to be highly dependent on launch SOP. Impairments due to SOPMD and higher-order PMD occur for larger signal bandwidth, particularly when these PMD components combine with chromatic dispersion or signal chirp [1].The chromatic dispersion and SOPMD can have a significant influence on pulse shape. Whereas FOPMD deals with two pulses delayed relative to each other, pulse distortions due to SOPMD can be visualized in terms of the interplay of six different pulse replicas. These mechanisms result in pulse overshoots and the generation of satellite pulses [2, 14].Section1 presents introduction of CD, FOPMD, SOPMD, and chirp effects. Section 2 discusses theoretical treatments prevalent in the area. Section 3 describes the suggested theory to find the reconstructed pulses under the effects of these effects. Section 4 gives special results, while Section 5 presents the results and discussion. Finally in Section 6, conclusions of the work done have been drawn.


## 2. Statement of the Problem

Consider a generic input fieldA→in(T)=Ain(T)|s〉, where Ain(T) represents the field complex amplitude and |s〉 represents the input Jones polarization vector, which depends on the azimuth θ and the ellipticity η of the input SOP. Mathematically, the SOP is defined as [15]
(1)|s〉=[ab]=[cosθcosη-isinθsinηcosθcosη+icosθsinη].The SOP and the pulse shape of the propagated pulses are affected by the FOPMD. The SOPMD induces other changes at the pulse properties. However, the SOPMD contains the PCD that may be taken as an additional CD. The origin CD is influenced by the chirp. These effects will be studied simultaneously in this paper to explain an accurate expressions of the reconstructed pulses.The frequency response of a single mode fiber with CD and PMD has the form [2, 14]
(2)H(w)=TPMD(w)exp[iβ2Lw22],
where β2 is the fiber’s GVD parameter, L is the fiber length, and β2 is the carrier frequency. The matrix TPMD(w) is responsible for PMD effects, which is a 2×2 unitary matrix that is defined as [16]
(3)TPMD=exp[-iw2τ→(w)·σ→],
where τ→(w) is the PMD vector, |τ→(w)| is the DGD between the PSPs, and σ→=(σ1,σ2,σ3) is the vector of Pauli spin matrices which can be defined as [1]
(4)σ0=[1001],σ1=[100-1],σ2=[0110],σ3=[0-ii0].
Here τ→·σ→=τ1σ1+τ2σ2+τ3σ3 is a 2×2 matrix. Now, the output field is given by
(5)A→out(w)=exp[iβ2Lw22]Ain(w)TPMD(w)|s〉=Aout(w)|t(w)〉,
where Ain(w) is the Fourier transform of Ain(T) and |t(w)〉=TPMD(w)|s〉 is the output SOP. Note that all the fiber birefringence effects are included in the frequency-dependent vector |t(w)〉, while the CD effects and pulse characteristics are accounted for inside the output field Aout(w).In the three-dimensional Poincare representation, PMD is described by the dispersion vectorτ→=Δττ^, where Δτ=|τ→| and τ^ is the direction of one of the two orthogonal PSP’s of TPMD(w). In a first-order approximation both Δτ and τ^ are frequency independent; that is, τ→=Δτ0τ^0, where Δτ0 is the DGD between two input signals polarized along the two PSP’s of the fiber. The Taylor series of τ→ vector up to second order is defined as [12, 13]
(6)τ→=Δτ0τ^0+(Δτ′τ^0+Δτ0τ^w)w,
where τ^0=τ^|w=w0, Δτ0=Δτ|w=w0, Δτ′=∂Δτ/∂w|w=w0, τ^w=∂τ^/∂w|w=w0, w0 is the central frequency, and w is the respective deviation to the central frequency w0.It is important to note thatτ^w is not unit vector and the vectors τ^0 and τ^w are orthogonal in Stokes space. The value of Δτ′ is responsible for PCD, while the factor τ^w represents the depolarization.


## 3. Theoretical Representation

In this section, we will explain the effects of PMD, chromatic dispersion, and initial chirp on the degree of polarization, state of polarization, and pulse shape. Substituting (6) into (3) and using the notation τ^w=2km^, one may obtain
(7)TPMD=exp[-i2(φ1τ^0+φ2m^)·σ→],
where φ1=w(Δτ0+Δτ′w) and φ2=w2Δτ02k. Under the assumption that the combined effect of PMD may be represented as a concatenation of FOPMD and SOPMD processes, (7) may be rewritten as follows:
(8)TPMD=(C1σ0-iτ^0·σ→S1)(C2σ0-im^·σ→S2),
where we have used the abbreviations C1, C2, S1, and S2 to represent cos(φ1/2), cos(φ2/2), sin(φ1/2), and sin(φ2/2), respectively. Using the property τ^0⊥m^ and the identity (a^·σ→)(b^·σ→)=a^·b^+ia^×b^·σ→, (8) will be
(9)TPMD=C1C2σ0-i(S1C2τ^0+S2C1m^+S1S2n^)·σ→,
where n^=τ^0×m^. Note that the vectors τ^0, m^, and n^ constituted the orthogonal basis of Stokes space. However, to advance analysis, let
(10)Q→=[S1C2τ1+S2C1m1+S1S2n1S1C2τ2+S2C1m2+S1S2n2S1C2τ3+S2C1m3+S1S2n3],
such that (9) will be
(11)TPMD=C1C2σ0-iQ→·σ→.
Using (11), it may be found that
(12)TPMD†TPMD=C12C22σ0+|Q→|2σ0.
Also, using (10), it is easy to get
(13)|Q→|2=S12C22+S22C12+S12S22.
Substituting (13) into (12) will make TPMD†TPMD=σ0. The eigenvalues of the matrix TPMD are
(14)λ1,2=C1C2∓C12C22-1.
In presence of PMD, we have C12C22<1, such that the eigenvalues are different and complex and their product is λ1λ2=1. Using these facts and the analysis presented in [16], we concluded a very important conclusion, which is the principle states of polarization (PSPs) vectors stay orthogonal but the orientation of antiparallel PSP’s vector may be changed. That is, the orthogonality property will not be affected by SOPMD.Return back to (11), which can be written as
(15)TPMD=[C1C2-iq1-i(q2-iq3)-i(q2+iq3)C1C2+iq1].
Using (11) and the definition of |s〉 from (1), the output SOP may be
(16)|t(w)〉=TPMD|s〉=[a(C1C2-iq1)-ib(q2-iq3)-ia(q2+iq3)+a(C1C2+iq1)].
Using (10) into (16) yields
(17)|t(w)〉=C1C2[ab]+S1C2[θ1θ4]+S2C1[θ2θ5]+S1S2[θ3θ6],
where
(18)[θ1θ4]=-iτ^0·σ→|s〉,[θ2θ5]=-im^·σ→|s〉,[θ3θ6]=-in^·σ→|s〉.In order to explain the effects of PMD on the signal DOP, it is interesting to rewrite (16) as follows:
(19)|t(w)〉=[txty]=[u1a+u2b-u2*a+u1*b],
where u1=C1C2-iq1 and u2=-i(q2-iq3).The corresponding Stokes vector is(20)t→=[t1t2t3]=[|tx|2-|ty|2txty*+tytx*i(txty*-tytx*)].
Since different frequency components of the input broadband signal will undergo different polarizations evaluation when they propagate through the fiber, the signal DOP will degrade. Consequently, the weighted average of the Stokes vector is necessary, which is calculated by multiplication the normalized spectral intensity function I=|Ain(w)|2 and integrating over w. The DOP will be [14]
(21)DOP=|∫-∞∞t1Idw|2+|∫-∞∞t2Idw|2+|∫-∞∞t3Idw|2∫-∞∞Idw.
Note that the parameters t1, t2, and t3 are functions of Δτ0, Δτ′, and k. That is, the signal DOP is a function of FOPMD and SOPMD. The integrations in (21) are very hard to solve analytically because the parameters t1, t2, and t3 are complicated functions of frequency. However, determining the analytic form of the signal DOP is beyond this paper.Now, the trigonometric functions in (17) will be simplified as follows:(22a)C1C2=14[ψ1+ψ1*+ψ2+ψ2*],(22b)S1C2=-i4[ψ1-ψ1*+ψ2-ψ2*],(22c)S2C1=-i4[ψ1-ψ1*-ψ2+ψ2*],(22d)S1S2=14[-ψ1-ψ1*+ψ2+ψ2*],where
(23)ψq=exp[i2(Δτ0w+αqw2)],q=1,2,α1=Δτ′+2kΔτ0, and α2=Δτ′-2kΔτ0. Substituting (18) and (22a)–(22d) into (17) will get
(24)|t(w)〉=14[ψ1|s1〉+ψ1*|s2〉+ψ2|s3〉+ψ2*|s4〉],
where(25a)|s1〉=[I+(in^-(τ^+m^))·σ→]|s〉,(25b)|s2〉=[I+(in^+(τ^+m^))·σ→]|s〉,(25c)|s3〉=[I+(-in^-(τ^-m^))·σ→]|s〉,(25d)|s4〉=[I+(-in^+(τ^-m^))·σ→]|s〉.Assuming that the input pulse is a normalized Gaussian pulse, this is defined as [1]
(26)A(T)=exp(-1+iC2T2T02),
where C and T0 are the initial chirp and pulse width. Substituting (24) and the Fourier transform of (26) into (5) yields
(27)A→out(L,w)=14[F1+|s1〉+F1-|s2〉...+F2+|s3〉+F2-|s4〉],
where
(28)Fq±=2πT021+iCexp[±iwΔτ02-(T02+αq±C-iαq±)w22(1+iC)].
Also(29a)α-1=β2L+α1=β2L+(Δτ′+2kΔτ0),(29b)α-2=β2L-α1=β2L-(Δτ′+2kΔτ0),(29c)α-3=β2L+α2=β2L+(Δτ′-2kΔτ0),(29d)α-4=β2L-α2=β2L-(Δτ′-2kΔτ0).The inverse Fourier transform on (27) will obtain
(30)A→out(L,T)=14[A1|s1〉+A2|s2〉+A3|s3〉+A4|s4〉],
where(31a)Aq=Mqexp(iφq)exp[-(T+(-1)q+1Δτ0/2)22Tq2],(31b)Tq=T01+2α-qCT02+α-q2(C2+1)T04,(31c)φq=12[tan-1(α-qT02+α-qC)-β],(31d)Mq=T0Tq,q=1,2,3,4,(31e)β=(T+(-1)q+1Δτ0/2)2[CT02+α-q(C2+1)](T02+α-qC)2+α-q2.Here, Tq, φq, and Mq represent the reconstructed pulse width, generated phase, and amplitude, respectively. Using the vectors in (25a)–(25d) into (30) and conducting a routine simplification, the output power equation takes the following formula:
(32)Pout(L,T)=14[|A1|2+|A2|2+|A3|2+|A4|2]+τ^0·s^2Re(A2A4*-A1A3*)+m^·s^4[|A2|2+|A3|2-|A1|2-|A4|2]-n^·s^2Im(A2A4*-A1A3*),
where Re and Im represent real and imaginary parts, respectively. Equation (32) represents the most important achievement in this paper as it contains the effects of CD, chirp, FOPMD, and SOPMD. It is clear from this equation that the second term results from the FOPMD effects and the third term results because of the SOPMD effects, while the fourth term results from the interaction between FOPMD and SOPMD. The first term has nothing to do with the SOP or PSPs, so it represents the CD effects in the absence of other effects. By the adoption of Table 1, (32) will be reduced to special formulas according to the following effects that accompany the propagated pulse.Table 1
Parameters specifications for different effects.
Case
Effects
Parameters specification
1
CD only
α
-
1
=
α
-
2
=
α
-
3
=
α
-
4
=
β
2
L
,
T
1
=
T
2
=
T
3
=
T
4
,
τ
^
0
=
m
^
=
n
^
=
0
,
Δ
τ
0
=
Δ
τ
′
=
2
k
=
0
2
FOPMD and CD
α
-
1
=
α
-
2
=
α
-
3
=
α
-
4
=
β
2
L
,
T
1
=
T
2
=
T
3
=
T
4
,
τ
^
0
≠
0
,
m
^
=
n
^
=
0
3
SOPMD and CD
α
-
1
=
-
α
-
2
=
Δ
τ
′
+
2
k
Δ
τ
0 α-3=-α-4=Δτ′-2kΔτ0 T1≠T2≠T3≠T4 If C≠0; T1=T2; T3=T4 If C=0, τ^0=n^=0,m^≠0
4
Orthogonal SOPMD and CD,2k≠0, Δτ′=0
α
-
1
=
α
-
4
=
β
2
L
+
2
k
Δ
τ
0
,
α
-
2
=
α
-
3
=
β
2
L
-
2
k
Δ
τ
0
,
T
1
=
T
4
;
T
2
=
T
3
τ
^
0
=
n
^
=
0
,
m
^
≠
0
5
Parallel  SOPMD and CD,2k=0, Δτ′≠0
α
-
1
=
α
-
3
=
β
2
L
+
Δ
τ
′
,
α
-
2
=
α
-
4
=
β
2
L
-
Δ
τ
′
,
T
1
=
T
3
;
T
2
=
T
4
τ
^
0
=
n
^
=
0
,
m
^
≠
0
6
SOPMD and FOPMD(general)
α
-
1
=
-
α
-
2
=
Δ
τ
′
+
2
k
Δ
τ
0
, α-3=-α-4=Δτ′-2kΔτ0 If C≠0 then  T1≠T2≠T3≠T4 If C=0 then T1=T2; T3=T4, τ^0⊥m^⊥n^≠0


## 4. Simplified Analytic Forms

In this section, we will address the mathematical analysis for special cases resulting from (32) to highlight theoretically the individual effects of the phenomena under study as well as to check the validity of the evidence contained in the previous section.(a) CD Only. In this simple case, (32) will be reduced to
(33)Pout(L,T)=T0T1exp-T2T12,
where T1 is computed from (31b) and Table 1. Here there is no PSPs because the CD has nothing to do with SOP, but this case is causing broadening or compressing of the output pulse depending on β2C sign. Note that this formula is well known in the scientific researches.(b) FOPMD and CD. In this case, (32) will be
(34)Pout(L,T)=T0T1[R-exp-T+2T12+R+exp-T-2T12],
where R∓=(1∓τ^0·s^)/2, T1 is as in the previous case, and T±=T±Δτ0/2. We note here that the PSP is exactly ±τ^0, and the CD will be broadened or compressed each of the two orthogonal components. Neglecting the CD effects will make the orthogonal components not suffer from any change expect the shifting to right or left by Δτ0/2. However, the resulted pulse will be broadened by Δτ0.(c) SOPMD and CD. Assuming that the SOPMD will be compensated, then
(35)Pout(L,T)=R-[|A1|2+|A4|2]+R+[|A2|2+|A3|2].
Here R∓=(1∓m^·s^)/2. In this case, the PSPs are ±m^, where there are two orthogonal components and each one contains a sum of two pulses that are shifted to right and left.(d) Orthogonal SOPMD and CD. Suppose that Δτ′=0; thus
(36)Pout(L,T)=T0[R-[exp-T+2T12+exp-T-2T12]/T1….+R+[exp-T+2T22+exp-T-2T22]/T2],
where R∓ is as in the previous case and T1,2 are computed from (31b) and Table 1. The PSPs are to be ±m^ here, where there are two orthogonal components and each one has different amplitudes and widths. This means that the orthogonal SOPMD will create the new PSPs, which are ±m^.(e) Parallel SOPMD and CD. Suppose that 2k=0; thus
(37)Pout(L,T)=T02[1T1exp-T+2T12+1T2exp-T-2T22],
where T1,2 are computed from (31b) and Table 1. There is no relationship here with the SOP or the PSPs. This confirms that the changes will be to the pulse width, where each component will be broadened or compressed, and the result is equivalent to a polarization-dependent CD. This result is physically similar to that of case (a).(f) FOPMD and SOPMD. This case involved general where all terms of (32) will be used in creating the resulting pulse. In this case, PSPs will not be ±τ^0, but they are created from both ±τ^0 and ±m^; that is, the presence of SOPMD will cause a change in the fiber optic PSPs and this is what is called depolarization.


## 5. Results and Discussion

Equation (32) has been used for having simulated the behavior of propagated pulse that the input pulse has been considered as a normalized Gaussian one (one can see (23)). The orthogonal vectors τ^0, m^, and n^ are randomly generated. The direction of s^ vector for these orthogonal ones is important such that it will determine the output pulse shape. If s^ is parallel to τ^0, the effects of SOPMD will not appear, while the case that s^ is parallel to m^ will prevent the effects of FOPMD. But as s^ is parallel to n^, the individual effects of FOPMD and SOPMD will be discarded save the last term of (32) which represents the interaction between two kinds of PMD. Since these vectors are random ones, the pulse shape differ from simulation process to another one. To determine the actual behavior, it should calculate the statistical average. However, this leads to not clarifying individual effects for parameters in research. So, the pulse’s behavior is proposed by using a specific simulation to point out the effects of the parameters β2, C, Δτw, and k. We have not focused on the investigation of Δτ0 as FOPMD is a well-covered in the scientific literatures. Some parameters are considered constant within a simulation process as L=50km, T0=10ps, and Δτ0=5ps, and the others parameters β2, Δτw, C, and k are changed within simulation. All figures are labeled by colors: blue, red, and green stand for the cases k=0,5,10ps, respectively.Figure1 illustrates the output pulse at Δτw=0 and β2=0 with various values of k and C. If C=0, the effects are only concerned to attenuation that increases by increasing k. This is in a good agreement with applied physics because k stands for normal component of SOPMD. The effects of k are shown in the form of broadening/compression and distortion in output pulses by increasing C. The calculated variables in Figure 2 are the same as those in Figure 1 except that β2=1ps2/km. As is seen that pulses are tended to broaden because the chromatic dispersion. The amount of broadening will raise by increasing the chirp. At higher value of k, the SOPMD tends to balance chromatic dispersion effect. So, the output pulses are shown less broadening at a higher value of k. As the pulse with high values of k is affected by the parallel component of SOPMD and chromatic dispersion, then it tends to be distorted.Figure 1
The resulted pulses withΔtw=0ps2 and β2=0ps2/km for different values of chirp and k.Figure 2
The resulted pulses withΔtw=0ps2 and β2=1ps2/km for different values of chirp and k.The calculated variables in Figure3 are the same as those in Figure 1 except that Δτw=30ps2. As is seen the action of the two components of SOPMD is given to rise broadening/compression and distortion which increase their rates with increasing C. For high values of k the pulses tend for broadening and vice versa. This is corresponding to applications as long the components of SOPMD are orthogonal. The broadening is attributed to the interaction among C, Δτw, and k. The calculated variables in Figure 4 are the same as those in Figure 3 except that β2=1ps2/km. Note that the effects become in contrast to these in Figure 3 because the presence of chromatic dispersion is about to rebalance among those parameters mentioned above so as to be projected with existing high values of C.Figure 3
The resulted pulses withΔtw=30ps2 and β2=0ps2/km for different values of chirp and k.Figure 4
The resulted pulses withΔtw=30ps2 and β2=1ps2/km for different values of chirp and k.The results are implicitly similar with those in Figure 1 in [12], where the power pulses are exceeded “1” as PSPs are rotated at a high value of C. Overshoot in power may happen by increasing the SOPMD effects, especially, with higher values of C.


## 6. Conclusions

The above equations and results illustrated the following conclusions. The short pulses are more affected by the initial chirp and vice versa. The SOPMD will add a new complexity that is presented by PCD and depolarization. The interaction between the FOPMD and SOPMD will cause a more degradation on the resulted pulses. The compensation of FOPMD will make the pulse under the affection of SOPMD only, which may destroy the pulse because of the presence ofΔτw and k. That is, the omitting or compensation of FOPMD may set serious variations. The increasing amount of k will make a more distortion on the propagated signal, while the variation of Δτw presents a new type of CD.The resulted pulses are affected, too, byβ2C, β2Δτ′, and CΔτ′. The presence of CD only broadens/narrows the resulted pulses. The mixing of CD and PMD pictures will alter the known behavior, where the values of Δτw and k are very important which enhanced/complicated the certain problem. However, the problem is very small for C=0, while the variations of pulse shape will be large for C≠0.



---

*Source: 345154-2013-12-11.xml*