PMD fundamentals: Polarization mode dispersion in optical fibers

1. Introduction

In the more than 15 years since the introduction of the early concepts (1, 2), the fundamentals of polarization mode dispersion (PMD) in optical fibers have become an important body of knowledge basic for the design of high-capacity optical communication systems. PMD effects are linear electromagnetic propagation phenomena occurring in so-called “single-mode” fibers. Despite their name, these fibers support two modes of propagation distinguished by their polarization. Because of optical birefringence in the fiber, the two modes travel with different group velocities, and the random change of this birefringence along the fiber length results in random coupling between the modes. With current practical transmission technology the resulting PMD phenomena lead to pulse distortion and system impairments that limit the transmission capacity of the fiber. Excellent reviews are available (3, 4), covering the practical aspects and applications of PMD concepts to fiber transmission systems and the effects of PMD on nonlinear fiber transmission (5). In this review we aim to complement these surveys and to collect and synthesize the fundamental concepts and theory of PMD, interweaving and linking the principal laws and key formulas that appear scattered in various places in the literature. We will explore the connection between frequency-domain and time-domain analyses and the isomorphic relation between the three-dimensional (3-D) view using real-valued 3-D Stokes vectors and the two-dimensional (2-D) view using complex-valued 2-D Jones vectors. Isomorphic pairings of operators such as these have been widely used elsewhere in physics such as in mechanics (6), in quantum mechanics (7), and even in the unification of quantum theory and general relativity (8). We borrow this methodology for our purposes.

As a preparation, we will first examine the description of polarization of light in the 3-D space of Stokes vectors (Stokes space) and the 2-D space of Jones vectors (Jones space, defined by the transverse coordinates x and y in the laboratory). Pauli spin matrices and spin vectors are the key to connect these two spaces and will be discussed in Appendix A, together with the necessary spin vector algebra. Because propagation through the fiber rotates the Stokes vectors of the light, we devote the following section to the rotation matrices in both Stokes and Jones space—i.e., the Lie groups SU(2) and SO(3) (9). Next, we examine the PMD vectors and principal states of polarization (PSPs) used to describe light propagation in randomly birefringent fibers and the associated differential group delay (DGD). We do this for a variety of representations and establish their connections. Several laws of infinitesimal rotation have been very valuable for visualizing PMD effects in Stokes space. Their connection to PMD vectors will be discussed. Finally, we review the various rules for PMD vector concatenation valuable for analyzing PMD of concatenated pieces of fiber. These rules appear in sum, integral, and differential form.

2. Notation

The unified view to which we aspire requires a unified notation. We attempt to keep our notation simple and transparent while linking to the notation already established as much as possible. The following is an abbreviated listing of our notation:

x, y, z: fiber coordinates; z is the direction of propagation; x and y are the transverse coordinates—i.e., those of Jones space.

e ^{j(ω₀ t−βz)}: continuous wave traveling in the z direction; j is the imaginary unit, ω₀ the angular carrier frequency, t time, and β the propagation constant.

E, Ẽ: electric field vectors; Ẽ(ω) is the Fourier transform of the complex transverse (x, y) electric field vector E(t) and has a complex amplitude e such that The vector of the real electric field is Re(E e ^{jω₀ t}).

ω: deviation from the angular carrier frequency ω₀ of the light. The optical frequency is (ω₀ + ω).

|s〉: 2-D complex Jones (column) ket vector, |s〉 = (_sy ^s_x). The bra 〈s| indicates the corresponding complex conjugate row vector—i.e., 〈s| = (s ^* _x, s ^* _y). The bra–ket notation is used to distinguish Jones vectors from Stokes vectors. Our Jones vectors are all of unit magnitude—i.e., 〈s|s〉 = s ^* _x s _x + s ^* _y s _y = 1.

ŝ: 3-D Stokes vector of unit length indicating the polarization of the field, and corresponding to |s〉. The component of ŝ are the Stokes parameters. By this definition, s ₃ = 1 for right-circular polarized light (s _y = js _x) conforming with the traditional optics definition (10). However, left-circular definitions are also used in the literature (11). We always use the same letter symbols for corresponding Jones and Stokes vectors. Note that a common phase shift of both components of |s〉 does not change ŝ.

I: 2×2 or 3×3 identity matrix. The distinction should be clear from context.

T: 2×2 unitary transmission matrix in Jones space. Relates output to input via We use the symbols s and t when necessary for clarity to designate respective input and output quantities, as illustrated in Scheme S1.

U: 2×2 Jones matrix, with det(U) = 1. Related to T by where φ₀ is the common phase.

R: 3×3 rotation matrix in Stokes space isomorphic to U. Relates output to input via σ₁, σ₂, σ₃: 2×2 Pauli spin matrices, for our purposes defined as σ⃗: Pauli spin vector in Stokes space, σ⃗ = (σ₁, σ₂, σ₃).

β⃗⋅σ⃗: 2×2 matrix in Jones space, β⃗⋅σ⃗ = β₁σ₁ + β₂σ₂ + β₃σ₃.

β⃗: 3-D birefringence vector in Stokes space describing local fiber properties.

τ⃗: our output PMD vector in Stokes space. Its length τ is the DGD, and its direction is that of the Stokes vector of the slow principal state. In a commonly used notation introduced by Poole (2, 3), the output PMD vector is labeled Ω⃗; its length, also the DGD, is labeled Δτ. The two PMD vectors τ⃗ and Ω⃗ are related by inversion through the s ₃ axis, for reasons discussed in Appendix B in the supplemental data at www.pnas.org.

p̂, q̂, r̂: unit Stokes vectors; p̂ and q̂ are sometimes used to describe the polarization of the slow and fast principal states, respectively, while r̂ is used for a rotation axis.

Subscript ω: indicates differentiation—i.e., ds/dω = s _ω.

3. Jones and Stokes Vectors

This preparatory section deals with the representation of polarization in Jones and Stokes spaces and the connection between the two. Throughout this paper we assume that there is no fiber nonlinearity and no polarization-dependent loss, and that the usual loss term of the fiber has been factored out so that we can deal with unitary transmission matrices, T and U, and 3-D rotation matrices, R, where the dagger denotes the Hermitian conjugate. Note that R has real elements, hence R ^† is simply the transpose of R. The input and output field components Ẽ _s(ω) and Ẽ _t(ω) at frequency ω are related by After dropping the field magnitudes and using the unit Jones vectors (12) to describe the state of polarization this becomes where |s〉 and |t〉 include phase and polarization information.

The Pauli spin matrices allow us to write the components s _i of the Stokes vector corresponding to |s〉 in the compact form (11) With the spin vector σ⃗ (9, 13) the Stokes vector is, simply, In the following we discuss several useful connections between the Jones and Stokes vectors, particularly the projection operator, dot products, and vector superpositions.

As a preparation, recall that the scalar product of two Jones vectors |p〉 and |q〉 is and that their dyadic operator is From Eqs. 3.6 and 3.7 one can see that where Tr stands for the trace operator.

Dot Products.

An immediate application of the projector rule 3.9 is the establishment of a connection between the dot (scalar) products of two Jones vectors |p〉 and |q〉 and of their associated Stokes vectors p̂ and q̂. We use 3.9 with s = p and multiply by |q〉 from the right and 〈q| from the left, obtaining A well-known special case of this connection is a pair of orthogonal Jones vectors with 〈q|p〉 = 0. Their Stokes vectors are antiparallel—i.e., q̂ = −p̂. In the following we will label such pairs of orthogonal Jones vectors, with antiparallel Stokes vectors, as |p〉 and |p ₋〉, for example. Any such pair forms a complete orthogonal set in Jones space, since we have where the second (completeness) relation follows from 3.9.

Let us now look at a pair of input Jones vectors |p _s〉, |q _s〉 and the transmitted vectors |p _t〉, |q _t〉 at the fiber output. They are related by Eq. 3.3. Using Eqs. 3.1 and 3.3, it is easy to show that the input and output dot products are preserved—i.e., When 3.13 is inserted in 3.11, we see that the magnitude of Stokes dot products is also preserved on transmission: i.e., transmission through the fiber is represented by a rotation of the Stokes vectors. This rotation property is, of course, a direct consequence of the orthogonality 3.1 of the matrix R.

4. Rotational Matrix Expressions

We have noted above that the change of the polarization of light on transmission through a fiber can be described as the rotation of its Stokes vector. There are matrix forms that highlight these rotational properties—i.e. the rotation axis r̂ and the rotation angle ϕ. These matrices are basic for an understanding of PMD fundamentals, and they have been used for the measurement of PMD vectors in the laboratory (14). We devote this section to a discussion of several expressions for such matrices in both Jones and Stokes space. First we establish a general connection between the matrices U and R of the two spaces, corresponding to the special unitary and special orthogonal Lie groups SU(2) and SO(3) (9).

Rotational Forms in Stokes Space.

Here we use the relation 4.2 to convert the rotational form 4.4 of the Jones matrix U into the isomorphic form for the matrix R in Stokes space. We insert 4.4 into the right-hand side of 4.2 and obtain having made use of the spin vector rules A.3, A.4, A.7, and A.13 to simplify the resulting expressions. Comparison of 4.6 with 4.2 yields the rotational form of R where the 3-D dyadic r̂r̂ is the projection operator and r̂× is the crossproduct operator One can see from 4.7 that for any Stokes vector ŝ, Rŝ represents a right-handed rotation of ŝ through an angle ϕ about the direction r̂. The Muller Matrix method for measuring PMD vectors uses this form for extracting rotation axis and angle from measured data for R (14).

The elegant compact expression for U given in Eq. 4.5 has an isomorphic counterpart in Stokes space (15) in Here the argument of the exponential is a 3×3 matrix operator. To prove the equivalence of this form with 4.7 one uses the power series expansion of the exponential, and the identities to collect the terms appearing in this expansion.

Elementary Rotations.

Elementary rotations (6) in Stokes space are those that rotate the Stokes vectors around the 1, 2, and 3 axes of the Poincaré sphere. The matrices describing those are special cases of 4.4 and 4.7. We call them U ₁, U ₂, U ₃ and R ₁, R ₂, R ₃. In 4.4, r̂⋅σ⃗ reduces to σ_i, yielding while in 4.7 only r _i = 1 is different from zero. The elements of U _i and R _i are listed in Table 1.

Note that U ₁ and R ₁ describe the rotation caused by a birefringent phase plate with the slow principal axis aligned with the x-axis in Jones space. U ₂ and R ₂ correspond to a phase plate set at 45° angle in Jones space. U ₃ and R ₃ describe a rotation by ϕ/2 around the z-axis.

5. PMD Vectors

PMD phenomena in optical fibers typically used for communications occur because of the presence of birefringence in the fiber. This birefringence changes randomly along the fiber length (3, 4). It stems from asymmetries in the fiber stress and geometry, such as elliptical cross sections, microbends, or microtwists. Often such a fiber is visualized or modeled as a sequence of random birefringent sections whose birefringence axes and magnitude change randomly with z (along the fiber). There are different manifestations of PMD depending on the view taken. In the frequency domain view one sees, for a fixed input polarization, a change with frequency ω of the output polarization. In the time domain one observes a mean time delay of a pulse traversing the fiber which is a function of the polarization of the input pulse. The two phenomena are intimately connected.

There exist special orthogonal pairs of polarization at the input and the output of the fiber called the PSPs. Light launched in a PSP does not change polarization at the output to first order in ω. These PSPs have group delays, τ_g, which are the maximum and minimum mean time delays of the time domain view. The difference between these two delays is called the DGD. Typical mean values of the DGD are 1 to 50 ps for a 500-km long fiber, depending on fiber type. The PMD vector τ⃗ describes both the PSPs and the DGD in the fiber. It is a Stokes vector pointing in the direction of the slow PSP with a length equal to the DGD.

Some insight into the PMD problem can be had simply by contemplating a piece of polarization-maintaining fiber. Its PSPs are the polarizations along the principal axes of birefringence of the fiber. In this case the two axes can be treated separately, and in general have different phase shifts φ and different group delays dφ/dω. One can see that the different values of φ(ω) will also produce changes in the output state of polarization as a function of frequency unless the input is launched in one of the PSPs. It may seem surprising that PSPs occur in a fiber exhibiting random birefringence as a function of z. However, they do, as discovered by Poole and Wagner (2). The DGD grows roughly as the square root of the length of fiber, as is characteristic of a random walk problem. This section is devoted to the discussion of the PMD vectors.

We will discuss four mathematical expressions for the PMD vector, one based on Jones matrix eigenanalysis, one on σ-expansion, one on time domain moments, and one on Müller matrices. These expressions correspond to four different yet closely related views of PMD and form the basis for different experimental approaches, such as time-domain and frequency-domain PMD measurements.

Jones Matrix Eigenvector Analysis.

The analysis of Poole and Wagner (2) is based on the Jones matrix U and the frequency-domain view, identifying the group delay of a narrow-band pulse propagating through the fiber with the frequency derivative of the phase of the field. As is the case for polarization-maintaining fiber, they found two output polarization states t̂ that do not change to first order in frequency. These states have different group delays, and also exhibit the maximum and minimum mean time delays for narrow-band pulses of arbitrary input polarization. Note that our PMD vector τ⃗ is defined for right-circular Stokes space, while the widely used PMD vector Ω⃗ of Poole et al. (2, 3) is defined for left-circular Stokes space. The connection between τ⃗ and Ω⃗ is detailed in Appendix B in the supplemental data at www.pnas.org. We have chosen the symbol τ⃗ as appropriate for a quantity with the dimension of time.

According to Eqs. 2.3 and 3.3 we have a transmission equation relating the input and output Jones vectors |s〉 and |t〉. As pulses are described by wave packets with a finite frequency band, we need to consider the frequency dependence of |t〉. We assume fixed input polarization and phase—i.e., |s〉_ω = 0 (hence ŝ_ω = 0), as is appropriate for a simple pulse entering the fiber at time zero. By differentiating 5.1 and eliminating |s〉, we obtain for the change of the output Jones vector Eq. 5.2 tells us that for most input polarizations, the output polarization will change with frequency in first order. The frequency derivative of the common phase φ₀ identifies a mean group delay τ₀ common to all polarizations Poole et al. (2, 3), in effect, noted that if |t〉 is either of the two orthogonal eigenstates of the operator jU _ω U ^†, then t̂_ω = 0, and so |t〉_ω should be expressible in the form where φ is the phase of |t〉. We show below that the operator jU _ω U ^† is Hermitian and that its trace is zero. Thus its eigenvalues are real and add to zero. We designate them τ/2 and −τ/2. Comparing 5.2 with 5.4, identifying dφ/dω in 5.4 with a group delay τ_g, and using 5.3, we get the two values of group delay The eigenstate of jU _ω U ^† associated with the larger value of group delay we designate |p〉, so that the orthogonal eigenstate having the smaller value of group delay becomes |p ₋〉. The slow PSP, |p〉, thus satisfies the Jones matrix eigenvector equation According to 5.5 the DGD is τ. Since the determinant of a Hermitian matrix is the product of its eigenvalues, we have det(jU _ω U ^†) = −τ²/4. From this, noting that det(U) = 1, we can extract an expression for τ, namely The PMD vector τ⃗ is the Stokes vector of |p〉 multiplied with the DGD

Pauli Spin Matrix Expansion.

We now pause to examine the operator jU _ω U ^†. Differentiation of the unitary condition 3.1 yields the relation Since the Hermitian conjugate of the product of two operators is the product of the Hermitian conjugates of the individual operators taken in the opposite order, we see that the matrix product jU _ω U ^† is Hermitian. We also have det(U) = 1. Consider the expansion Since the determinant of the product of two matrices is the product of their determinants, and det(U) = 1, one can see that det(U(ω + dω)) = 1 also only if the trace of U _ω U ^† is zero. This is the argument necessary to arrive at 5.5. Given these properties of jU _ω U ^†, and its eigenvalues ±τ/2, it follows that it can be written in the form where the last expression is in the spin vector form. To show that 5.11 and 5.8 identify the same vector τ⃗, we substitute 5.11 into the eigenvalue equation 5.6, yielding In view of 3.10, we observe that τ⃗ = τp̂, as in 5.8, is the solution of 5.12. The σ-expansion of the Jones matrix product of 5.11 is, thus, a very practical tool, because it yields the PMD Stokes vector components as the coefficients of the expansion.

In Section 4 we derived the spin vector expansion of a Jones matrix U in terms of the rotational variables r̂ and ϕ. By substituting 4.4 into 5.11 we can gain an expression for τ⃗ in terms of these variables and their frequency derivatives. The result, using the spin vector rules A.5 and A.6 to reduce the intermediate factors, and the identity r̂_ω⋅r̂ = 0, is When U is given in the Caley/Klein form we obtain the DGD from 5.7 as (2) and for the components τ_i

Moments and Mean Signal Delay.

Moments are widely used to characterize pulse transmission, and they have also been applied to PMD phenomena (15–18). Here we describe their use for determining and defining PMD vectors. One should note that moments are very appropriate for the description of signal delays in systems involving two (or more) modes of propagation such as those involving PMD effects, whereas concepts such as group velocity and group delay have a strict definition for an individual mode only. The PMD phenomenon can split an input pulse into two or more pulses at the fiber output, leading to polarization-dependent pulse shapes. To maintain our awareness of this complexity we use the word “signal” for the output rather than “pulse.”

To prepare for the discussion of moments we represent a signal at fiber position z by the complex field vector E(z, t) in the time domain and by its Fourier transform Ẽ(z, ω) in the frequency domain. The tilde distinguishes between them. where all integrals extend from −∞ to +∞. Let the fields be normalized so that the energy in the pulse is where the right equation is Parseval's theorem.

The mean time t _g(z) at which a signal passes location z is defined by its center of gravity or first moment W ₁ Using the convolution theorem, the first moment can be written in either the time domain or the frequency domain We now use the transmission and Jones matrices to express the output field Ẽ _t in terms of the input Ẽ _s, Ẽ _t = T Ẽ _s, and use the derivatives, to substitute in 5.22 and 5.23. Defining the mean signal delay τ_g as the difference between the mean arrival times at output (t) and input (s), one gets where τ₀ = dφ₀/dω as before. Substituting the σ-expansion 5.18 for jU ^† U _ω = ½τ⃗_s⋅σ⃗, the mean signal delay can be written in the form We proceed by representing the input field Ẽ _s by its complex amplitude e, its unit Jones vector |s〉, and the corresponding Stokes vector ŝ (see Eq. 3.5), where the signal energy becomes W = ∫dωe*e. Inserting this above, we obtain (15, 18) for the mean signal delay where τ⃗_s is the input PMD vector as before. This expresses τ_g as the spectral mean (denoted by the 〈 〉 bracket) of the spectral density e*e weighing τ₀(ω) and the dot product between input PMD vector τ⃗_s(ω) and the input Stokes vector. This expression for the mean signal delay is valid for any input pulse shape and spectral variation of its polarization.

Specializing to frequency-independent input polarization ŝ and a narrow pulse spectrum such that τ₀ and τ⃗_s can be regarded constant over the signal band, 5.29 simplifies to (16) This agrees with 5.5, which was derived for the special case of input polarizations aligned with the PSPs. Eq. 5.30 gives the mean signal delay for any input polarization. Inputs at the PSPs are seen to lead to the maximum and minimum delays τ_g. In fact, the mean signal delay τ_g can be simply interpreted as the power-weighted average of the two PSP delays. To see this, consider the representation of the input polarization as the superposition 3.14, |s〉 = a|p〉 + b|p ₋〉, of the PSPs. These delays are (τ₀ ± τ/2) and the power-weighted average delay is We can now use the dot-product rule 3.11 to show the equivalence of 5.31 and 5.30. The simple form of 5.30 allows a determination of the PMD vector τ⃗_s in the time domain (19).

The higher moments W _n = ∫dt t ⁿ E ^† E can be used to determine further detail on the properties of transmitted pulses. The second moment provides information about rms pulse spreading. Proceeding along similar lines as above and the substitution of 5.24, we obtain From 5.25 it follows that where we have used spin-vector rule A.8. To simplify 5.33 further we use 5.28, assume frequency-independent input polarization |s〉, and characterize the complex input amplitude e(ω) by its real amplitude a(ω) and phase δ(ω) The result is For the special case of δ_ω = 0, this form agrees with the result of Karlsson (15). This special case assumes freedom from frequency chirp in the optical pulse and a symmetric pulse shape, or, more precisely, that E*(−t) = E(t).

Müller Matrix Expression for τ⃗.

A fourth way of determining the PMD vector uses Stokes space and the Müller (rotation) matrix R isomorphic to U. The background for this will be given in Section 6. For completeness we show the result here. It is where × denotes the cross-product operator as in 4.8.

This method relies on the determination of the rotation matrix R and its derivative in the frequency domain.

6. Laws of Infinitesimal Rotation

The laws of rotation assume a particularly simple form when the rotations are infinitesimally small—e.g., for a small change in fiber length or frequency. These are well known in mechanics (6) and are widely used in the related phenomena of PMD, coupled modes of propagation, and two-waveguide systems (3, 20–24). Here these laws allow simple geometrical interpretations in Stokes space (on the Poincaré sphere). We shall discuss infinitesimal laws for the birefringence vector, the PMD vector, and the dynamical PMD equation.

Birefringence Vector.

Consider the change of the polarization |s〉(z) of light at fiber location z due to a small length addition dz of fiber. This change is influenced by the fiber's local birefringence characterized by its effective relative dielectric tensor ɛ(z)—i.e., a cross-sectional average of the fiber characteristics for the fiber mode of interest. The change is governed by the wave equation for a spectral component of the effective transverse field vector Ẽ(z) of the mode, where k ₀ = 2π/λ₀ is the propagation constant of free space, and λ₀ is the free-space wavelength. To proceed, we use a σ-expansion of the ɛ-tensor of the form where β₀ is the common propagation constant. The coefficients β_i of the expansion are the components of the local birefringence vector β⃗(z) in Stokes space. This vector has the character of a propagation constant and has been a useful tool in describing birefringence in PMD phenomena (3, 22, 25) (see Appendix B in supplemental data at www.pnas.org for the connection of β⃗ to the vector W⃗ used in part of the literature). In addition to this expansion we use an adiabatic approximation assuming that the polarization, |s〉(z), and ɛ(z) all vary slowly with z, and by setting where |s〉 includes a slowly varying phase. We continue by inserting 6.2 and 6.3 into 6.1. We drop the d ²|s〉/dz ² term in accordance with the adiabatic assumption and find the adiabatic wave equation for the Jones vector |s〉(z), The right-hand side of 6.2 shows the form of the fiber's dielectric tensor as expressed by the components β_i of the birefringence vector β⃗. It is apparent that optical activity (circular birefringence) is included in this description. For the special case of linear birefringence aligned with the fiber's x-axis we have β₂ = β₃ = 0 and β₁ = Δnk ₀, where Δn is the differential refractive index. The corresponding ordinary and extraordinary indices are There is a close relationship between the birefringence vector β⃗ and the concept of a local normal mode used in the theory of optical waveguides for the analysis of guides whose characteristics change with length z, such as our ɛ(z). At a given location z ₀, the local normal modes are defined as the modes of a uniform guide with the characteristics of location z ₀—i.e., a guide with uniform ɛ(z) = ɛ(z ₀) in our case. Inserting this assumption into 6.1, we obtain the field of the two local normal modes Ẽ(z) = e ^{−jβ_M z}|s _M〉 whose polarization |s _M〉 does not change along the length. Because of 3.10 we know that |s _M〉 must be an eigenvector of β⃗⋅σ⃗. The polarization of the local normal modes is, thus, described by the Stokes vectors ±β⃗. The propagation constant β_M of the modes depends on their polarization in the form assuming that β₀ ≫ β.

We now return to the case of interest where β⃗(z) changes along the fiber length z. To describe the change of polarization with z we differentiate the Stokes vector ŝ = 〈s|σ⃗|s〉 and obtain Finally, we combine 6.4 and 6.7, use the spin vector rules A.6, A.7, and A.13, and obtain the law of infinitesimal rotation for birefringence Integration of 6.8 yields a formal expression for the Müller matrix R of the fiber as in 2.5.

PMD Vector τ⃗.

Here we consider the change of polarization at the fiber output due to a small change in frequency ω. The polarization at the input is held constant. We start with Eq. 5.2 written as for the change of the output polarization |s〉_ω, where the subscript ω indicates differentiation. The differential of the corresponding Stokes vector is Finally, we combine 6.9 and 6.10, using the σ-expansion 5.11 for the product U _ω U ^†, apply spin vector rule A.6, A.7, and obtain the law of infinitesimal rotation The geometrical interpretation of this simple law is a rotation of the output Stokes vector on the Poincaré sphere as ω changes. The rotation axis is the PSP p̂ and the rotation rate is the DGD τ.

The infinitesimal rotation law 6.11 allows us to express the PMD vector in terms of the rotation (Müller) matrix R that relates the input and output Stokes vectors ŝ(0) and ŝ(z) by ŝ(z) = Rŝ(0). We differentiate this relation [while keeping the input ŝ(0) fixed] and obtain Comparison of 6.11 and 6.12 yields the operator relationship which we have already listed in 5.37.

7. PMD Vector Concatenation Rules

The PMD vector concatenation rules (13, 16, 22, 23, 26) are a powerful set of simple tools that allow the determination of the PMD vector of an assembly of concatenated fiber sections when the PMD vectors of the individual sections are known. Among their uses is the analysis of the evolution of the PMD vector with fiber length (26), statistical PMD modeling (23), PMD simulation, and the design of multisection PMD compensators. They appear in a variety of related forms, including the sum, differential, and integral formulations for both the first-order and second-order PMD vectors. As a preparation for the discussion of these formulations, we shall first review the relations between the output and input PMD vectors of a single fiber section.

Transformation of PMD Vectors.

Consider a fiber section with rotation matrix R and input and output Stokes vectors ŝ and t̂ as shown in Scheme S2. The corresponding PMD vectors at input and output are τ⃗_s and τ⃗. They are related by as discussed earlier (5.19). The expression for the output vector τ⃗ in terms of R was already given in 6.13; it is As cross-product operators transform like matrices from output to input we obtain the expression for the input PMD vector as When higher-order PMD effects are considered, τ⃗ is usually called the “first-order” PMD vector, its frequency derivative τ⃗_ω is the “second-order” PMD vector, the second derivative τ⃗_ωω is the “third-order” PMD vector, etc. The input/output transformation for the second-order PMD vectors τ⃗_sω and τ⃗_ω is derived by differentiating 7.1 and substituting 7.2, giving with the result This shows that the input and output second-order PMD vectors transform the same way as the first-order vectors of the section. For the third-order PMD vectors one gets the somewhat more complicated relationship

PMD of Two Concatenated Sections.

Here, we consider two concatenated fiber sections with rotation matrices R ₁ and R ₂ as shown in Scheme S3. The output PMD vectors of each of the individual sections are τ⃗₁ and τ⃗₂.

Our goal is to determine the output PMD vector of the combined assembly. The rotation matrix R of the two-section assembly is the matrix product Combining this with 7.2, we find for the PMD vector τ⃗ of the assembly that Now we recall that R ₁ R ₁ ^† = I, apply 7.2 for the individual sections as well as a transformation of matrix operators as in 5.18 to simplify 7.7 into the form This is the basic concatenation rule. It can be generalized to multiple sections as well as differentially small sections. It can also be transformed to the input of the fiber, or, for that matter to any desired fiber cross section. The rule is very similar to that for impedances of a transmission line: to get the PMD vector of an assembly, transform the PMD vectors of each individual section to a common reference cross section and take the sum of all those vectors.

The concatenation rule for the second-order PMD vector τ⃗_ω follows from 7.9 by differentiation and substitutions similar to those above. One obtains

The Differential Concatenation Rule.

Here, we apply the two-section concatenation rule to a long piece of fiber with output PMD vector τ⃗ and a differentially small fiber addition of length Δz as shown in Scheme S4.

We want to answer the question: How does τ⃗(z) change due to the differential addition of length Δz? The answer to this question is already provided by the dynamical PMD equation 6.15. Our discussion will show that the latter is identical to the differential form of the concatenation rule. It will establish a link between this rule and the laws of infinitesimal rotation treated in Section 6. For this purpose we need to determine the rotation matrix R _Δ and the PMD vector τ⃗_Δ of the differential addition Δz. We use the local birefringence vector β⃗ of the addition and the law of infinitesimal rotation (6.8) and get for the change Δŝ of the Stokes vector ŝ entering the element In terms of R _Δ, the output Stokes vector ŝ + Δŝ can be expressed as Comparison of 7.11 and 7.12 yields the desired expression for the rotation matrix of the differential addition To obtain its PMD vector τ⃗_Δ we differentiate 7.13, apply 7.2, and drop terms of second order in Δz with the result The PMD concatenation rule 7.9 for the addition of the differential fiber element is Now we can insert 7.13 and 7.14, divide by Δz, and get which is the dynamical PMD equation 6.15 on transition to infinitesimal Δz.

The differential concatenation rule for the second-order PMD vector is obtained along similar lines or by differentiation of 7.16. It is

Concatenation Rules for Many Sections.

Consider now a fiber consisting of m individual sections, which is sketched in Scheme S5.

The rotation matrices R _n and the output PMD vectors τ⃗_n of each section are considered known. The PMD vectors τ⃗ and τ⃗_ω of the combined assembly can be determined by repeated application of the two-section rules 7.9 and 7.10. The resulting expressions simplify if we define the rotation matrix of the last m − n + 1 sections as where R(m, m) = R _m and R(m, m + 1) = I, and the output PMD vector of the first n sections as τ⃗(n).

Using these definitions, the sum rules for the first- and second-order PMD vectors of the assembly are and