Fiber_Optics_Physics_Technology
.pdf
42 |
Chapter 3. Treatment with Wave Optics |
Figure 3.7: Observed intensity distribution of all modes existing at V = 6. For reproduction purposes the gray scale of the original photograph has been reduced to a binary black and white. From [145] with kind permission.
becomes much smaller than the core radius and we approach the multimode case. This confirms our previous heuristic assumption.
The LP01 mode exists all the way down to arbitrarily small V , i.e., at any arbitrarily long wavelength. None of the other modes has this property. However, the existence of the fundamental mode down to V = 0 must not be taken literally: We have used the approximation that the fiber cladding is infinitely wide. However, at some very long wavelength the field will penetrate far enough into the cladding to reach the outside surface. There is always a practical limit for the longest wavelength supported in the fiber, often dictated by wavelengthdependent losses.
3.11Number of Modes
In order to find the total number of modes at any given V , we have to note that there are degeneracies. For example, any mode can exist in two distinct,
3.12. A Remark on Microwave Waveguides |
43 |
mutually orthogonal polarization states which are identical in any other respect as long as we stick to our approximation that fibers are circularly symmetric. (We will later look at small deviations from this symmetry.) Then, all LP0p modes are actually pairs (doublets) of modes.
For m = 0, we have to additionally note that the azimuthal dependence of the solution can be written either with sin or with cos; again, these are two mutually orthogonal variants. Taking both this and the polarization degeneracy into account, LPmp with m = 0 are actually quartets of modes. Let us consider, as an example, the situation at V = 6: From Fig. 3.3 we see six nominal modes, two of which are pairs, and four, quartets. The total number thus is 20. Asymptotically one can approximate
V → ∞ total number of modes = V 2/2 (step index fiber).
For gradient index fibers, there is a similar approximation with V 2/4.
Strictly speaking, the modes do not have precisely linear polarization. This is due to the fact – neglected in our approximation – that the fiber is not a perfectly homogenous material but has a step in its refractive index. This leads to distortions of the field which produces some deviation of the modes, mostly for higher-order modes. We can live with that because we are mostly interested in the fundamental mode.
3.12A Remark on Microwave Waveguides
The reader may or may not be aware that discrete modes also exist in microwave waveguides. Microwave waveguides are metal pipes with conducting walls. This enforces a node of the electrical field on the boundary. In contrast, optical fibers are weakly guiding conduits. Therefore we could use here an approximation which is not valid in microwave guides whereas there one finds di erent types of modes and uses a di erent terminology [37]. Many of the modes derived here are linear combinations of metallic waveguide modes; the following table presents the correspondence:
LP modes |
Microwave guide modes |
LP01 |
HE11 |
LP11 |
HE21, EH01 |
LP21 |
HE31, EH11 |
LP02 |
HE12 |
LP31 |
HE41, EH21 |
Nevertheless, similarities do exist between optical fibers and microwave waveguides. The biggest di erence may be that they always have a minimum frequency even for the fundamental mode; below, no mode is supported at all. This can be traced directly to the conducting walls. In Sect. 4.5.2, we will present a particular case in which a fiber with special structure has a well-defined finite lower cuto even for the fundamental mode, too.
3.13Energy Transport
We have calculated the modal structure of fibers under the assumption of circular symmetry. Waveguiding arises from the guiding of modes by the refractive
44 |
Chapter 3. Treatment with Wave Optics |
index structure. As soon as a fiber is bent, the circular symmetry is broken. And, as may be expected, then additional energy loss arises; bending losses are discussed in Chap. 5.2. Here we can already note this much:
The fiber core may be decisive for waveguiding, but it would be an oversimplification that the light power is guided in the core exclusively. We have seen that the field amplitude decays radially like an exponential function; this implies that there is always a certain fraction of power well outside the core. Di erent modes have di erent geometric field distributions; the fraction outside the core must therefore also be di erent for di erent modes.
The energy transport out of (or into) a certain volume element is described by the Poynting vector
S = E × H. |
(3.65) |
The direction of propagation is perpendicular to the plane containing E and H. Disregarding anisotropic materials, all three vectors are pairwise orthogonal. The reader is reminded that by convention, the direction of polarization designates the direction of the oscillation of E (historically there was once a convention to refer to H, but that has long been obsolete).
In most cases, one describes the energy transport of a wave by giving its irradiance or intensity. By this, one means the temporal average of the instantaneous intensity
Iinst = E(t) H(t).
I gives power per unit area and thus has units of W/m2. Obviously I is the temporal average of the Poynting vector:
I = |S| .
In the special case that the wave is harmonic, the relations between rms and |
||||||||||
ˆ |
√ |
|
ˆ |
|
√ |
|
|
|
|
|
|
2 He . Then we find the intensity as |
|||||||||
peak values are E = |
|
2 Ee and H = |
|
|
|
|||||
|
|
|
I = |
1 |
|
ˆ ˆ |
|
|||
|
|
|
|
2 |
E H. |
|||||
Using the relation |
|
|
ˆ |
|
|
|
|
|
ˆ |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
E = Z0 H, |
|
||||||
one can also write |
|
|
|
|
1 |
|
|
|||
|
|
|
I = |
|
ˆ2 |
|
||||
|
|
|
2Z0 |
E |
. |
|||||
When the wave propagates in nonmagnetic matter with real index n, this becomes
ˆ2
E .
Power P is the integral of intensity over the cross-sectional area s:
P = |
I ds. |
|
s |
After a calculation which we will not show in detail here (see [50, 120]), one finds the energy fraction in the core (PK) and in the cladding (PM) for the fundamental mode:
3.13. Energy Transport |
45 |
P /P K
1.0
|
|
|
|
51 |
|
|
|
|
41 |
|
|
|
31 |
|
0.5 |
|
21 |
12 |
22 |
|
11 |
03 |
||
01 |
|
02 |
|
|
0.0 |
|
|
|
|
0 |
2 |
4 |
6 |
8 |
V
Figure 3.8: Frequency dependence of the partition of power between core and cladding. The relative power inside the core is shown as a function of V for all modes up to V = 8. Modes with m = 0, 1 are guided when only a tiny fraction of power is in the core (at their cuto these curves begin at zero). It is true for all modes that for large V the fraction approaches unity. After [50] with kind permission.
PK
P
PM
P
= |
|
w |
2 1 + |
|
J0(u) |
2 |
||||||
|
|
|
||||||||||
|
V |
|
|
|
|
J1(u) |
|
|||||
|
|
|
|
w |
|
2 |
|
|
J0(u) |
2 |
||
= |
1 − |
|
|
|
|
1 + |
|
|
||||
V |
|
|
J1(u) |
|||||||||
For very large V , the mode is strongly concentrated in the core. As V decreases, the field begins to spread into the cladding. At the mode cuto , practically all the power is in the cladding. This is precisely what causes the loss of waveguiding. For modes with m ≥ 2, waveguiding is lost already when the power fraction outside the core exceeds 1/m. Figure 3.8 illustrates the situation.
We have seen that fields of di erent modes have di erent cladding penetrations. For any given mode, the degree of spreading into the cladding depends on wavelength. This observation is also relevant for the fiber’s dispersion, which will be treated in Chap. 4.
Chapter 4
Chromatic Dispersion
A light signal propagating in an optical fiber is subject to a variety of ways in which it can get distorted. Many of these are based on di erent propagation velocities for di erent parts of the signal. After such distortion, there is a risk that the signal arrives at the receiver in such a mangled form that it may be impossible to correctly decipher it.
We already encountered one such mechanism: modal dispersion in multimode fibers. Now we will address such distortions as they arise in single-mode fibers.
At the center of explanation is the fact that the refractive index of the fiber glass, just like that of any other material, depends on wavelength (or frequency). No light signal is ever truly monochromatic; rather, it contains Fourier components from a certain spectral interval. In other words, a light pulse of finite duration by necessity has a nonzero spectral width. Di erent frequency components, however, will propagate with di erent velocity. This gives rise to di erential transit time and thus to signal distortion called delay distortion.
The wavelength dependence of the refractive index is behind three di erent contributions to delay distortion. They are collectively called chromatic dispersion. Individually, they are
Material dispersion. Dm: This contribution arises directly from the wavelength dependence of the index. Material dispersion is not specific to fibers but can be found in any bulk glass. It is independent of geometry and (given the material) depends solely on wavelength.
Waveguide dispersion. Dw: In the particular geometry of optical fibers, there is a modification to the di erential propagation time. The reader is reminded that we saw in Chap. 3.13 that the signal power is partitioned between core and cladding; the splitting ratio depends on the wavelength. On the other hand, core and cladding indices are slightly di erent. As the wavelength is varied, we have a crossover from mostly core index to mostly cladding index. The result is a contribution to the wavelength dependence of the e ective index.
Profile dispersion. Dp: Strictly speaking, the index di erence between core and cladding itself, and thus Δ, is also wavelength-dependent. (Core index and
F. Mitschke, Fiber Optics, DOI 10.1007/978-3-642-03703-0 4, |
47 |
c Springer-Verlag Berlin Heidelberg 2009
48 |
Chapter 4. Chromatic Dispersion |
cladding index do not vary “in parallel.”) This gives rise to another correction which, however, is often much smaller than material and waveguide dispersion.
Another reason for dispersive distortions in single-mode fibers is related to the state of polarization of the light. As mentioned above, each mode can be decomposed into two mutually orthogonal parts. An ideal fiber has perfect circular symmetry; then both polarization states (polarization modes) propagate with identical velocities. However, real-world glass always has at least some residual birefringence; this implies a slightly di erent e ective index for both states. One can argue that the term “single-mode fiber” is a misnomer: Even when it is true that only a single geometric field amplitude distribution (LP01) can propagate, it still consists of two polarization modes. This is why in real fibers there is polarization mode dispersion. Typically it is a small contribution; we will discuss it further below.
To characterize dispersion, one normally specifies the size of the e ect per distance. For both modal and polarization mode dispersion, one can write the dispersion parameter
D = |
1 |
|
δτ , |
(4.1) |
|
L |
|||||
|
|
|
|||
where δτ designates the di erence of propagation time after distance L. Units of ps/km are commonly used. For chromatic dispersion, including material, waveguide, and profile dispersion, the following specification is common:
D = |
1 |
|
dτ |
. |
(4.2) |
|
|
||||
|
L dλ |
|
|||
Here, D contains the three parts just mentioned:
D = Dm + Dw + Dp . |
(4.3) |
This dispersion parameter indicates the propagation time di erence per distance and per wavelength di erence; therefore, units of ps/(nm km) are commonly used.
4.1Material Dispersion
For any glass, the refractive index varies with wavelength. This gives rise to chromatic defects of lenses and the color-separating capability of prisms. For historical reasons, it became common practice (see, e.g., Schott glass catalog [11]) to characterize types of glasses by giving their indices at three wavelengths:
nD, the refractive index at wavelength 589.30 nm (yellow, Fraunhofer’s D line of sodium),
nF at wavelength 486.13 nm (blue–green, Fraunhofer’s F line of hydrogen), and
nC at wavelength 656.27 nm (red, Fraunhofer’s C line of hydrogen).
4.1. Material Dispersion |
49 |
This choice was motivated by the availability of narrowband emission lamps at these wavelengths. As a further characterization, often the Abbe number is given, defined by
vD = |
nD − 1 |
. |
(4.4) |
|
nF − nC |
|
|
Obviously this is a metric of the wavelength dependence of the index near the central (yellow) wavelength.
Glass catalogs often specify only nD and vD. It is instructive, however, to survey the variation of the index over a wide spectral range. This is schematically shown in Fig. 4.1. There are absorption bands due to electronic transitions in the ultraviolet at wavelengths on the order of 100 nm and molecular vibrations in the infrared around 10 μm. In the vast interval in between, including all of the visible and near-infrared, pure silica glass does not exhibit any resonances. This is the reason, of course, why it appears “crystal-clear” to the eye. The position of the absorption resonances is reflected in the refractive index.
Figure 4.1: Schematic plot of refractive index vs. wavelength. It is dominated by resonances which occur both in the ultraviolet and the infrared.
Pure fused silica (SiO2) has a refractive index of nD = 1.456, decreasing slightly toward longer wavelengths. As long as one stays clear of the resonances, one can empirically describe the wavelength dependence with interpolation formulas. One of the best known such formulas is Sellmeier’s equation
m |
A |
λ2 |
|
|||
|
|
|||||
n2(ω) = 1 + |
|
|
j |
|
, |
(4.5) |
|
2 |
|
2 |
|||
j=1 |
λ |
|
− λj |
|
||
but there are also variants to this. Coe cients Aj denote the resonance strengths and λj the pertaining wavelengths. These coe cients are tabulated in the literature (for various glass types, see [11], and for fibers with various doping materials and concentrations, see e.g., [22]). In most cases three terms in the sum are considered su cient; sometimes, five. Let us emphasize again that a Sellmeier curve is an empirical fit: the coe cients λj must not be construed to indicate the resonance wavelengths in a literal sense.
50 |
Chapter 4. Chromatic Dispersion |
4.1.1Treatment with Derivatives to Wavelength
We now turn to the actually observed propagation times and the scatter thereof. Consider a plane monochromatic wave with angular frequency ω and wave number β. It is well known that it propagates with phase velocity
vph = ω/β , |
(4.6) |
whereas the propagation of a signal, like a wave packet, is governed by the group velocity
vgr = dω/dβ . |
(4.7) |
|||||||||
The group propagation time then is |
|
|
|
|
|
|
|
|
|
|
τ = |
|
L |
= L |
dβ |
(4.8) |
|||||
vgr |
dω |
|||||||||
|
|
|
|
|
|
|||||
= |
L |
dβ |
|
dλ |
. |
(4.9) |
||||
|
|
|||||||||
|
|
|
dλ dω |
|
||||||
Since β = nk0 = 2πn/λ, we can rearrange the first fraction in the last line:
|
dβ |
= 2π |
d |
|
n |
= |
2π |
λ |
dn |
|
− |
n . |
(4.10) |
||||||||||||||
|
dλ |
|
|
λ |
|
λ2 |
|
dλ |
|
||||||||||||||||||
|
|
|
dλ |
|
|
|
|
|
|
|
|||||||||||||||||
The second fraction can be rearranged using λ = 2πc/ω to yield |
|
||||||||||||||||||||||||||
|
|
|
dλ |
= 2πc |
|
d |
|
1 |
|
= |
|
|
|
2πc |
. |
(4.11) |
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
dω |
|
dω |
ω |
|
− ω2 |
|
|
||||||||||||||||||
Combining, we obtain |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
τ = L |
2π |
|
2πc |
|
n |
− |
λ |
dn |
|
|
(4.12) |
|||||||||||||
|
|
|
λ2 |
|
dλ |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
ω2 |
|
|
|
|
|
|
|
|
|
|
||||||||
or, since λω = 2πc,
τ = |
L |
n − λ |
dn |
(4.13) |
||
|
|
|
. |
|||
c |
dλ |
|||||
We have now expressed the wavelength dependence of the group propagation time as a function of the easily measured quantities n and λ. We emphasize that we here find the bracketed expression taking the role of the usual index n which only appears in the first term of that expression. Therefore we call the bracketed expression the group index:
dn
ngr = n − λ dλ (4.14)
Figure 4.1 shows that throughout the visible and near infrared, n decreases with increasing wavelength; thus, in this range ngr > n. Figure 4.2 shows a comparison of both indices for fused silica.
The scatter of arrival times at the receiver is obtained from
δτ = |
dτ |
δλ . |
(4.15) |
|
dλ |
||||
|
|
|
4.1. Material Dispersion |
51 |
1,48
ngr
n
1,45
0,5 |
1,0 |
1,5 |
2,0 |
|
|
λ / μm |
|
Figure 4.2: Refractive index and group index as a function of wavelength, calculated from a three-term Sellmeier’s equation.
Herein, the derivative is
dτ |
= |
L |
|
d |
n |
− |
λ |
dn |
|
dλ |
|
|
|
||||||
|
c dλ |
|
dλ |
||||||
Ld2n
=− c λ dλ2 .
The contribution of the material dispersion Dm to the dispersion coe cient ,
|
Dm = |
1 |
|
|
dτ |
, |
(4.16) |
|||
|
|
|
|
|||||||
|
|
L dλ |
|
|||||||
is then |
|
|
|
|
|
|
|
|
|
|
|
|
λ d2n |
|
|||||||
|
Dm = − |
|
|
|
. |
(4.17) |
||||
|
c |
dλ2 |
||||||||
In some cases the signal may occupy a quite broad spectral band. Then one must take into account that D varies with wavelength: Dispersion would not be described to su cient accuracy by D alone. In such case one can additionally specify the dispersion slope:
Sm = |
dDm |
. |
(4.18) |
|
|||
|
dλ |
|
|
4.1.2Treatment with Derivatives to Frequency
An alternative terminology to describe dispersion arises when one takes derivatives not with respect to wavelength but with respect to frequency. One starts from a series expansion of the propagation constant β
β(ω) = n(ω) ωc = β0 + β1(ω − ω0) + 12 β2(ω − ω0)2 + ··· (4.19)
with
β |
= |
dmβ |
|
. |
(4.20) |
|
|
m |
|||||
m |
|
dω |
|
|
|
|
|
|
|
ω=ω0 |
|
|
|
|
|
|
|
|
|
|
52 |
Chapter 4. Chromatic Dispersion |
Let us assess the meaning of βm:
|
|
β0 = β(ω = ω0) = kn |
|
(4.21) |
||||||||||||
where n is the ordinary (phase) index. |
|
|
|
|
|
|
|
|||||||||
β |
= |
|
|
dβ |
|
|
|
|
|
|
|
|
|
|||
1 |
|
|
|
dω |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
ω=ω0 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
d |
|
|
|
ω |
|
|
|
|
|
||
|
= |
|
|
|
n(ω) |
|
|
|
|
|||||||
|
dω |
c |
|
|
|
|||||||||||
|
|
|
|
|
ω=ω0 |
|
|
|
||||||||
|
|
|
1 |
n(ω |
|
|
|
dn(ω) |
|
|
||||||
|
= |
|
) + ω |
. |
||||||||||||
|
|
|
|
|||||||||||||
|
|
|
c |
0 |
|
|
|
|
dω |
|
|
|||||
|
|
|
|
|
|
|
|
ω=ω0 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
By inserting ω = 2πc/λ, this can be written as |
|
|
|
|||||||||||||
|
β1 = |
1 n(λ0) |
|
λ dn(λ) |
|
. |
(4.22) |
|||||||||
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
c |
|
|
dλ |
λ=λ0 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
The bracketed expression is the group index ngr as introduced in Eq. (4.14), which implies that
|
|
β1L = |
|
L |
ngr = τ |
(4.23) |
|||||||||||
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
||||
holds. Then |
|
|
|
|
|
|
1 |
|
|
|
ngr |
|
|
|
|
||
|
|
β1 |
= |
|
|
= |
. |
|
|
(4.24) |
|||||||
|
|
|
|
|
|||||||||||||
|
|
|
|
|
vgr |
c |
|
||||||||||
β1 is of the order of β1 ≈ 5 ns/m. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
For β2 we find |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
β2 |
= |
dβ1 |
|
= |
1 |
2 |
dn |
+ ω |
d2n |
. |
(4.25) |
||||||
dω |
c |
dω |
|
||||||||||||||
|
|
|
|
|
|
|
|
dω2 |
|
||||||||
This quantity is called group velocity dispersion parameter or GVD parameter. It is commonly given in units of ps2/km. The GVD parameter is preferred by theorists over the dispersion coe cient D which is widely used by technicians. To convert, one uses
Dm = |
dβ1 |
= |
dβ1 |
|
dω |
= β2 |
d |
|
2πc |
|
(4.26) |
dλ |
dω dλ |
dλ |
|
||||||||
|
|
|
|
λ |
|
||||||
and thus
Dm = −β2 |
2πc |
= − |
ω |
β2 . |
(4.27) |
λ2 |
λ |
There are some cases in which higher-order dispersion terms become relevant. Then one can also specify the third-order dispersion (TOD) β3. To convert between dispersion slope S and β3, one can use
(2πc)2 |
|
4πc |
|
|
||
Sm = |
|
β3 + |
|
β2 . |
(4.28) |
|
λ4 |
λ3 |
|||||
|
|
|
|
|||
