Chromatic and Modal Dispersion in Fibers

A fiber carries information as a train of light pulses, and the whole game is keeping those pulses distinct after hundreds of kilometres of glass. Each pulse spreads in time, and once neighbours overlap the receiver can no longer tell a one from a zero — intersymbol interference sets a hard ceiling on the data rate. This broadening is dispersion, and because it piles up with distance the true figure of merit is not bit rate alone but the bit-rate–distance product, quoted in Gb s−1⋅km\mathrm{Gb\,s^{-1}\cdot km}. Two distinct mechanisms do the spreading, and it takes results from earlier in this chapter to separate them. From the guided-mode picture (aq-optics-waveguide-modes-wave-equation) we inherit that a fiber supports a discrete set of modes, each with its own propagation constant β\beta. Modal dispersion is the spread in arrival time between those modes; it lives only in multimode fibers and vanishes the instant single-mode operation (aq-optics-single-multimode-operation) leaves just one. Chromatic dispersion survives even then: any real source emits a band of wavelengths, a single mode's group velocity depends on wavelength, and so the colours of one pulse drift apart. It splits again — a material part from the glass index n(λ)n(\lambda) and a waveguide part from how tightly the mode is confined. Everything is organized by one function: the mode propagation constant β(ω)\beta(\omega). Its first derivative is the group delay per length; its second derivative β2\beta_2 — the group-velocity dispersion — is the pulse-broadening engine, repackaged as a coefficient DD in ps nm−1 km−1\mathrm{ps\,nm^{-1}\,km^{-1}} that you multiply by linewidth and length. In standard silica the material and waveguide parts cancel near 1.31 μm1.31\,\mu\mathrm{m}, while loss bottoms out near 1.55 μm1.55\,\mu\mathrm{m} — a mismatch that has shaped fiber design for forty years.

Key takeaways
  • Dispersion broadens optical pulses as they travel, forcing neighbouring bits to overlap and capping the bit-rate–distance product B⋅LB\cdot L.
  • Modal dispersion comes from guided modes having different group velocities; it exists only in multimode fiber and is killed by single-mode operation (V<2.405V<2.405).
  • Chromatic dispersion comes from one mode's group velocity depending on wavelength together with a source of finite linewidth; it splits into material and waveguide parts.
  • The coefficient D=−2πcλ2β2D=-\tfrac{2\pi c}{\lambda^2}\beta_2 crosses zero near 1.31 μm1.31\,\mu\mathrm{m}, but the loss minimum sits at 1.55 μm1.55\,\mu\mathrm{m}, where D≈+17 ps nm−1km−1D\approx+17\,\mathrm{ps\,nm^{-1}km^{-1}} (β2≈−21.7 ps2 km−1\beta_2\approx-21.7\,\mathrm{ps^2\,km^{-1}}).
  • Chromatic broadening is στ=∣D∣ L σλ\sigma_\tau=|D|\,L\,\sigma_\lambda; step-index modal spread is ≈n1ΔL/c\approx n_1\Delta L/c (tens of ns km−1\mathrm{ns\,km^{-1}}), which graded-index profiles cut by ∼Δ/2\sim\Delta/2.
  • Key idea: pulse broadening is governed by how the group delay β1=dβ/dω\beta_1=\mathrm{d}\beta/\mathrm{d}\omega varies — across modes (modal dispersion) or across the source spectrum (chromatic dispersion, set by β2\beta_2). Single-mode operation removes the first; wavelength and profile design minimize the second.
The core insight

Pulse broadening is governed by how the group delay β1=dβ/dω\beta_1=\mathrm{d}\beta/\mathrm{d}\omega varies: across modes (modal dispersion) or across the source spectrum (chromatic dispersion, set by β2\beta_2). Single-mode operation kills the former; the choice of wavelength and index profile minimizes the latter.

Where this lives today

The zero-dispersion window at 1.31 μm1.31\,\mu\mathrm{m} and the loss minimum at 1.55 μm1.55\,\mu\mathrm{m} do not coincide, and that mismatch is the axis around which optical telecommunications is engineered. Modern long-haul systems run at 1.55 μm1.55\,\mu\mathrm{m} to exploit the ∼0.2 dB km−1\sim0.2\,\mathrm{dB\,km^{-1}} loss and erbium-doped fiber amplifiers, then fight the residual D≈17 ps nm−1km−1D\approx17\,\mathrm{ps\,nm^{-1}km^{-1}} with dispersion-compensating fiber, chirped fiber Bragg gratings, or — in today's coherent transponders — digital signal processing that inverts the accumulated β2L\beta_2 L phase after detection. Because that phase is deterministic, a DSP chip can undo hundreds of kilometres of chromatic dispersion, which is why single-channel rates have climbed past 400 Gb s−1400\,\mathrm{Gb\,s^{-1}} on fiber whose DD has not changed since the 1980s. Elsewhere the same physics is a feature: in the anomalous regime above λ0\lambda_0, β2<0\beta_2<0 lets dispersion cancel the Kerr self-phase-modulation chirp, forming solitons and — in photonic-crystal fiber and microresonators — the supercontinua and frequency combs that now anchor optical atomic clocks, while data centres run graded-index multimode fiber with cheap VCSELs for tens of gigabits over a few hundred metres.

The physical picture

Start with the ray picture of a fat multimode step-index core. Different rays take different paths: one runs straight down the axis, another zig-zags at the steepest angle total internal reflection allows. Both travel at the same speed c/n1c/n_1, but the zig-zag ray covers a longer path, so it arrives late. Because the extra path grows with length, the spread grows linearly with LL. That is modal dispersion, set purely by geometry: the fractional index step Δ\Delta between core and cladding.

Figure
Modal dispersion in a step-index multimode fiber: an axial ray and a ray near the critical angle both travel at speed c/n₁ but cover different path lengths, so one launched pulse arrives smeared by Δτ ≈ L n₁ Δ / c ≈ 49 ns per km for n₁ = 1.48 and Δ = 0.01
Same speed c/n₁, different path length → the steep ray arrives late cladding n₂ cladding n₂ core n₁ one pulse in axial (fast) steep (slow) Δτ ≈ 49 ns/km t Δτ = L n₁ Δ / c · n₁ = 1.48, Δ = 0.01

Now silence every mode but one. The spreading does not stop, because a pulse is a wave packet built from a band of frequencies, and the glass index — plus the way the field spills between core and cladding — changes with wavelength. Each colour has its own group velocity, and the red and blue edges pull apart. This is chromatic dispersion, unavoidable in a real single-mode fiber, and its two parts fight in silica: material dispersion swings from negative to positive through the near-infrared while waveguide dispersion stays negative, so their sum crosses zero near 1.31 μm1.31\,\mu\mathrm{m} — moving that crossing is what dispersion-shifted fibers do.

Figure
Chromatic dispersion of standard single-mode silica fiber: material D_M and waveguide D_W add to a total D that passes through zero near λ₀ ≈ 1.31 μm and reaches D ≈ +17 ps nm⁻¹ km⁻¹ at the 1.55 μm loss minimum
The punchline: D = D_M + D_W, zero near λ₀ ≈ 1.31 μm, +17 at 1.55 μm +20 +10 0 −10 −20 1.1 1.2 1.3 1.4 1.5 1.6 1.7 vacuum wavelength λ (μm) D (ps nm⁻¹ km⁻¹) D_M D_W D = D_M + D_W +17 at 1.55 μm λ₀ ≈ 1.31 μm

Think of β2\beta_2 as a chirp generator: it makes the instantaneous frequency sweep across the stretched pulse — blue crowded to one edge, red to the other. That chirp is deterministic and reversible, which is why a short length of opposite-sign fiber can squeeze the pulse back. Chromatic dispersion is a distortion you can undo — the fact that underwrites dispersion compensation, and the bridge to solitons, where β2\beta_2 is balanced against the Kerr nonlinearity instead of cancelled.

The core equations

The starting object is the mode's propagation constant β(ω)\beta(\omega), the axial phase shift per unit length. Because a pulse occupies a narrow band around a carrier ω0\omega_0, we Taylor-expand β\beta there; the coefficients are the dispersion parameters that follow.

β(ω)=β0+β1(ω−ω0)+12β2(ω−ω0)2+⋯(1)\beta(\omega)=\beta_0+\beta_1(\omega-\omega_0)+\tfrac{1}{2}\beta_2(\omega-\omega_0)^2+\cdots \qquad (1)

where: β(ω)\beta(\omega) — axial propagation constant of the guided mode; ω0\omega_0 — carrier angular frequency; β1=dβ/dω\beta_1=\mathrm{d}\beta/\mathrm{d}\omega — inverse group velocity 1/vg1/v_g; β2=d2β/dω2\beta_2=\mathrm{d}^2\beta/\mathrm{d}\omega^2 — group-velocity dispersion (GVD). The linear term β1\beta_1 sets when the pulse arrives; the quadratic β2\beta_2 sets how much it spreads. The first-order coefficient is the reciprocal group velocity. It packages the familiar group index — the same ngn_g used elsewhere for pulse travel time.

vg=1β1=cng,ng=n−λdndλ(2)v_g=\frac{1}{\beta_1}=\frac{c}{n_g},\qquad n_g=n-\lambda\frac{\mathrm{d}n}{\mathrm{d}\lambda} \qquad (2)

where: vgv_g — group velocity (speed of the pulse envelope); ngn_g — group index; nn — phase (material) index; λ\lambda — vacuum wavelength. A pulse's delay per length is β1=ng/c\beta_1=n_g/c, and dispersion is just the variation of this delay with wavelength. Engineers rarely quote β2\beta_2 directly; they quote the dispersion coefficient DD, the change in group delay per length per unit wavelength. It follows from β2\beta_2 by the chain rule and, for the material part alone, reduces to the curvature of the index.

D=−2πcλ2 β2=−λcd2ndλ2(3)D=-\frac{2\pi c}{\lambda^2}\,\beta_2=-\frac{\lambda}{c}\frac{\mathrm{d}^2 n}{\mathrm{d}\lambda^2} \qquad (3)

where: DD — chromatic dispersion parameter, ps nm−1 km−1\mathrm{ps\,nm^{-1}\,km^{-1}}; β2\beta_2 — GVD parameter, ps2 km−1\mathrm{ps^2\,km^{-1}}; cc — vacuum light speed. The two forms are the same physics with opposite signs (positive DD, anomalous, means β2<0\beta_2<0); the rightmost shows material dispersion is literally the curvature of the Sellmeier index. With DD in hand, pulse broadening is a one-line estimate: multiply the coefficient by length and source linewidth. It is the number a link designer reaches for first.

στ=∣D∣ L σλ(4)\sigma_\tau=|D|\,L\,\sigma_\lambda \qquad (4)

where: στ\sigma_\tau — RMS temporal broadening; LL — fiber length; σλ\sigma_\lambda — RMS spectral width of the source. Linear in length and linewidth, so a narrow-line laser and a short span each help in direct proportion. The total chromatic DD is the sum of two pieces: a material term fixed by the glass and a waveguide term fixed by the geometry. The waveguide part enters through the normalized parameters VV and bb.

D=DM+DW,DW=−n2Δcλ Vd2(Vb)dV2(5)D=D_M+D_W,\qquad D_W=-\frac{n_2\Delta}{c\lambda}\,V\frac{\mathrm{d}^2(Vb)}{\mathrm{d}V^2} \qquad (5)

where: DMD_M — material dispersion from d2n/dλ2\mathrm{d}^2 n/\mathrm{d}\lambda^2; DWD_W — waveguide dispersion; bb — normalized propagation constant; Δ\Delta — fractional index difference. Since V d2(Vb)/dV2>0V\,\mathrm{d}^2(Vb)/\mathrm{d}V^2>0 for the fundamental mode, DWD_W is always negative and grows with confinement — the knob that shifts λ0\lambda_0. Modal dispersion has no λ\lambda in it — it is pure geometry. Comparing the axial ray with the critical-angle ray gives the delay spread directly.

Δτstep=L n1c Δ=Lc n1(n1−n2)n2(6)\Delta\tau_{\text{step}}=\frac{L\,n_1}{c}\,\Delta=\frac{L}{c}\,\frac{n_1(n_1-n_2)}{n_2} \qquad (6)

where: n1n_1 — core index; n2n_2 — cladding index; Δ=(n1−n2)/n1\Delta=(n_1-n_2)/n_1 — fractional index difference; LL — length. Linear in LL and Δ\Delta, typically tens of ns km−1\mathrm{ns\,km^{-1}} — orders of magnitude above chromatic spread, so multimode links are short. Whatever the mechanism, broadening becomes a ceiling on the bit-rate–distance product. Bits must be spaced by a few times στ\sigma_\tau, which fixes the maximum B⋅LB\cdot L.

B L≲14 ∣D∣ σλ(7)B\,L\lesssim\frac{1}{4\,|D|\,\sigma_\lambda} \qquad (7)

where: BB — bit rate; LL — distance; DD — dispersion coefficient; σλ\sigma_\lambda — source linewidth. The product B⋅LB\cdot L, not BB alone, is conserved — doubling the reach halves the rate. Finally, modal dispersion is switched off entirely by forcing single-mode operation. The cutoff condition, inherited from the waveguide-modes analysis, caps the VV-number.

V=2πaλn12−n22<2.405(8)V=\frac{2\pi a}{\lambda}\sqrt{n_1^2-n_2^2}<2.405 \qquad (8)

where: VV — normalized frequency (V-number); aa — core radius; 2.4052.405 — first zero of J0J_0, the cutoff of the LP11\mathrm{LP}_{11} mode. Below it only the fundamental LP01\mathrm{LP}_{01} propagates, so modal dispersion is exactly zero and only chromatic remains.

Derivation AFrom the propagation constant to the dispersion coefficient

The aim is to connect the three ways people talk about chromatic dispersion — the GVD parameter β2\beta_2, the engineer's coefficient DD, and the material curvature d2n/dλ2\mathrm{d}^2 n/\mathrm{d}\lambda^2 — and show they are one quantity in three costumes. We start from the group delay per length, differentiate with respect to wavelength (sources are specified by linewidth in nm), and carry the plane-wave index through to the curvature term.

Step 1. The group delay per length is τ=β1\tau=\beta_1, so the coefficient is D=dβ1/dλ=(dω/dλ) β2D=\mathrm{d}\beta_1/\mathrm{d}\lambda=(\mathrm{d}\omega/\mathrm{d}\lambda)\,\beta_2, using dβ1/dω=β2\mathrm{d}\beta_1/\mathrm{d}\omega=\beta_2. Why: dispersion is the wavelength-dependence of arrival time, and the chain rule routes it through β2\beta_2, the quantity the wave equation hands us.

Step 2. From ω=2πc/λ\omega=2\pi c/\lambda, dω/dλ=−2πc/λ2\mathrm{d}\omega/\mathrm{d}\lambda=-2\pi c/\lambda^2, so D=−(2πc/λ2) β2D=-(2\pi c/\lambda^2)\,\beta_2. Why: this is exactly why DD and β2\beta_2 carry opposite signs — the minus comes from ω\omega and λ\lambda running opposite.

Step 3. For the material part, use β=nω/c\beta=n\omega/c for a plane wave in the glass, giving β1=(1/c)(n−λ dn/dλ)=ng/c\beta_1=(1/c)(n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda)=n_g/c. Why: the group index re-appears exactly as in equation (2), tying the fiber delay to the material group index.

Step 4. Differentiate β1\beta_1 with respect to λ\lambda: the nn and single-derivative terms cancel by the product rule, leaving DM=−(λ/c) d2n/dλ2D_M=-(\lambda/c)\,\mathrm{d}^2 n/\mathrm{d}\lambda^2. Why: only the curvature survives, so material dispersion vanishes where the index has an inflection — near 1.276 μm1.276\,\mu\mathrm{m} in silica. *As a check: at λ=1.55 μm\lambda=1.55\,\mu\mathrm{m}, D=17 ps nm−1km−1D=17\,\mathrm{ps\,nm^{-1}km^{-1}} gives β2=−Dλ2/(2πc)=−21.7 ps2 km−1\beta_2=-D\lambda^2/(2\pi c)=-21.7\,\mathrm{ps^2\,km^{-1}}, the standard SMF value.*

Go deeper · Derivation B · 1 minModal delay spread from the ray picture
Derivation BModal delay spread from the ray picture

The goal is to pin down the size of modal dispersion in a step-index multimode fiber using nothing but ray geometry, and to show it is set entirely by the index step — no wavelength anywhere. We compare the two extreme rays the core supports: the axial ray that runs straight, and the ray that skims at the critical angle. The difference in their transit times is the delay spread that smears a pulse.

Step 1. The axial ray travels the length LL at speed c/n1c/n_1, so its transit time is tmin⁡=Ln1/ct_{\min}=Ln_1/c. Why: this is the fastest route — straight down the axis — setting the lower edge of the arrival window.

Step 2. The steepest guided ray is incident at the critical angle θc\theta_c with sin⁡θc=n2/n1\sin\theta_c=n_2/n_1; its zig-zag path is longer by 1/sin⁡θc1/\sin\theta_c, giving tmax⁡=(Ln1/c)(n1/n2)t_{\max}=(Ln_1/c)(n_1/n_2). Why: it covers the same axial distance along a slant, and 1/sin⁡θc1/\sin\theta_c is how much longer that slant is.

Step 3. The spread is Δτ=tmax⁡−tmin⁡=(Ln1/c)(n1/n2−1)\Delta\tau=t_{\max}-t_{\min}=(Ln_1/c)(n_1/n_2-1), which for weak guidance Δ≪1\Delta\ll1 becomes Δτ≈Ln1Δ/c\Delta\tau\approx Ln_1\Delta/c. Why: with n1≈n2n_1\approx n_2, n1/n2−1≈Δn_1/n_2-1\approx\Delta, so the spread carries no λ\lambda — modal dispersion is purely geometric. *As a check: n1=1.48n_1=1.48, Δ=0.01\Delta=0.01 gives Δτ/L=n1Δ/c=49 ns km−1\Delta\tau/L=n_1\Delta/c=49\,\mathrm{ns\,km^{-1}} — hundreds of times the chromatic spread, hence multimode reach is short.*

Go deeper · Derivation C · 1 minGaussian pulse broadening set by β₂ L
Derivation CGaussian pulse broadening set by β₂ L

Here we prove the headline result that group-velocity dispersion stretches a pulse by an amount controlled by the single combination β2L\beta_2 L, and that the process is lossless and reversible. In the slowly varying envelope approximation, in a frame moving at the group velocity, the envelope obeys a Schrödinger-like equation in which β2\beta_2 plays the role of an inverse mass. Propagating an unchirped Gaussian lets us read the output width off directly and identify the dispersion length.

Step 1. Keeping terms through β2\beta_2 and moving to retarded time T=t−β1zT=t-\beta_1 z, the envelope satisfies ∂A/∂z=−i2β2 ∂2A/∂T2\partial A/\partial z=-\tfrac{i}{2}\beta_2\,\partial^2 A/\partial T^2. Why: subtracting the group delay β1z\beta_1 z leaves β2\beta_2 as the sole agent of distortion — a diffusion equation with an imaginary coefficient.

Step 2. In frequency this integrates trivially: A~(z,ω)=A~(0,ω)exp⁡(i2β2ω2z)\tilde A(z,\omega)=\tilde A(0,\omega)\exp(\tfrac{i}{2}\beta_2\omega^2 z), a pure quadratic phase. Why: the amplitudes are untouched, so dispersion only re-times the spectral components — which is why it is perfectly reversible.

Step 3. Feed in an unchirped Gaussian A(0,T)=exp⁡(−T2/2T02)A(0,T)=\exp(-T^2/2T_0^2); the integral back to time gives a complex width T02−iβ2zT_0^2-i\beta_2 z. Why: the imaginary part is the induced chirp, and its growth with zz stretches the pulse.

Step 4. The real half-width is T1(z)=T01+(β2z/T02)2T_1(z)=T_0\sqrt{1+(\beta_2 z/T_0^2)^2}, negligible until z∼LD=T02/∣β2∣z\sim L_D=T_0^2/|\beta_2| and linear beyond. Why: since LD∝T02L_D\propto T_0^2, short pulses disperse fastest, so high bit rates suffer most. *As a check: for z≫LDz\gg L_D, T1→∣β2∣z/T0T_1\to|\beta_2|z/T_0, reproducing στ=∣D∣Lσλ\sigma_\tau=|D|L\sigma_\lambda once σλ\sigma_\lambda and T0T_0 are linked at the transform limit.*

Limiting cases

Which mechanism dominates depends on the fiber type, the source, and the wavelength:

CaseResultMeaning
Multimode step-indexΔτmodal≈n1Δc L\Delta\tau_{\text{modal}}\approx\dfrac{n_1\Delta}{c}\,LModal dispersion rules; linear in LL and Δ\Delta, tens of ns km−1\mathrm{ns\,km^{-1}}, limiting reach to short links.
Multimode graded-indexΔτmodal≈n1Δ22c L\Delta\tau_{\text{modal}}\approx\dfrac{n_1\Delta^2}{2c}\,LAn optimized profile equalizes mode delays, shrinking the spread by ∼Δ/2\sim\Delta/2 versus step-index.
Single-mode fiberστ=∣D∣ L σλ\sigma_\tau=|D|\,L\,\sigma_\lambdaNo modal dispersion; chromatic DD and linewidth set the limit, residual is PMD ∝L\propto\sqrt{L}.
Zero-dispersion wavelengthD(λ0)=0⇒στ∝∣S∣ L σλ2D(\lambda_0)=0\Rightarrow\sigma_\tau\propto|S|\,L\,\sigma_\lambda^2At λ0≈1.31 μm\lambda_0\approx1.31\,\mu\mathrm{m} first-order DD vanishes; the slope S=dD/dλS=\mathrm{d}D/\mathrm{d}\lambda sets the residual.

Worked examples

1 · Modal budget of a step-index multimode link (confidence builder)

A step-index multimode fiber has core index n1=1.48n_1=1.48 and cladding n2=1.465n_2=1.465. Find the modal delay spread per kilometre and the modal-limited bit-rate–distance product.

Step 1. Fractional index difference: Δ=(n1−n2)/n1=(1.48−1.465)/1.48=0.01014\Delta=(n_1-n_2)/n_1=(1.48-1.465)/1.48=\mathbf{0.01014} — the single geometric number that sets modal dispersion.

Step 2. Delay spread per length from equation (6): Δτ/L=n1Δ/c=(1.48)(0.01014)/(2.998×108)=50 ns km−1\Delta\tau/L=n_1\Delta/c=(1.48)(0.01014)/(2.998\times10^8)=\mathbf{50\ ns\,km^{-1}}. Over 1 km1\,\mathrm{km} the fastest and slowest rays arrive 50 ns50\,\mathrm{ns} apart — enormous next to chromatic spreads in ps.

Step 3. With bits spaced by a few Δτ\Delta\tau, taking B≲1/(4 Δτ)B\lesssim1/(4\,\Delta\tau) gives B⋅L≲1/(4×50×10−9 s km−1)=5×106 (bit/s)⋅kmB\cdot L\lesssim1/(4\times50\times10^{-9}\ \mathrm{s\,km^{-1}})=5\times10^{6}\ \mathrm{(bit/s)\cdot km}. So Δ=0.0101\Delta=0.0101, Δτ/L≈50 ns km−1\Delta\tau/L\approx50\,\mathrm{ns\,km^{-1}}, and B⋅L≲5 (Mb/s)⋅kmB\cdot L\lesssim5\ \mathrm{(Mb/s)\cdot km} — a 1 km link caps near 5 Mb/s, which is why step-index multimode is confined to short, low-rate runs.

2 · Chromatic broadening on a single-mode long-haul span

Standard single-mode fiber at 1.55 μm1.55\,\mu\mathrm{m} has D=17 ps nm−1km−1D=17\,\mathrm{ps\,nm^{-1}km^{-1}}. A DFB laser has σλ=0.1 nm\sigma_\lambda=0.1\,\mathrm{nm}. Find β2\beta_2, the broadening over L=100 kmL=100\,\mathrm{km}, and the bit-rate ceiling.

Step 1. Convert to GVD: β2=−Dλ2/(2πc)=−(17×10−6 s m−2)(1.55×10−6 m)2/(2π×2.998×108)=−21.7 ps2 km−1\beta_2=-D\lambda^2/(2\pi c)=-(17\times10^{-6}\ \mathrm{s\,m^{-2}})(1.55\times10^{-6}\ \mathrm{m})^2/(2\pi\times2.998\times10^8)=\mathbf{-21.7\ ps^2\,km^{-1}}; the negative sign marks the anomalous regime above λ0\lambda_0.

Step 2. Broadening from equation (4): στ=∣D∣ L σλ=(17)(100)(0.1)=170 ps\sigma_\tau=|D|\,L\,\sigma_\lambda=(17)(100)(0.1)=\mathbf{170\ ps} — every kilometre adds 1.7 ps1.7\,\mathrm{ps}, accumulating linearly.

Step 3. Bit-rate ceiling: B≲1/(4στ)=1.47 Gb s−1B\lesssim1/(4\sigma_\tau)=\mathbf{1.47\ Gb\,s^{-1}}, equivalently B⋅L≲1/(4∣D∣σλ)≈147 (Gb/s)⋅kmB\cdot L\lesssim1/(4|D|\sigma_\lambda)\approx147\ \mathrm{(Gb/s)\cdot km}. So β2=−21.7 ps2 km−1\beta_2=-21.7\,\mathrm{ps^2\,km^{-1}}, στ=170 ps\sigma_\tau=170\,\mathrm{ps}, and B≲1.5 Gb s−1B\lesssim1.5\,\mathrm{Gb\,s^{-1}} over 100 km without compensation — thousands of times better than the multimode link, yet far below the rates modern systems reach with compensation.

3 · How far before a pulse doubles?

A 50 ps50\,\mathrm{ps} (half-width T0T_0) transform-limited pulse enters the same fiber, β2=−21.7 ps2 km−1\beta_2=-21.7\,\mathrm{ps^2\,km^{-1}}. Over what distance does it broaden appreciably?

Step 1. The natural scale is the dispersion length LD=T02/∣β2∣=(50)2/21.7=115 kmL_D=T_0^2/|\beta_2|=(50)^2/21.7=\mathbf{115\ km}; below it the pulse is essentially unbroadened.

Step 2. At z=LDz=L_D the width grows to T1=2 T0=71 psT_1=\sqrt{2}\,T_0=\mathbf{71\ ps}, a 41% stretch. A 50 ps pulse survives ∼100 km\sim100\,\mathrm{km} before broadening by half; since LD∝T02L_D\propto T_0^2, a 10 ps pulse would broaden as much in only ∼4.6 km\sim4.6\,\mathrm{km} — the quadratic penalty on high bit rates.

Physics of the portrait

The theory that gives the b(V)b(V) curves behind the coefficient DWD_W came out of Bell Labs at the moment low-loss fiber first became real — the work of Detlef Gloge.

Origins · 1 minDetlef Gloge and how this was discovered
DG
Detlef Gloge
Bell Labs, Crawford Hill · weakly guiding fibers, 1971 · fiber waveguide theory and dispersion

When Corning drew the first low-loss fiber in 1970, the theory to design real systems was still missing: the exact vector modes of a fiber are a thicket of Bessel functions with no transparent picture of dispersion. Detlef Gloge, at Bell Laboratories' Crawford Hill site in Holmdel, cut through it in 1971 with the weakly guiding approximation. Recognizing that a communications fiber has n1≈n2n_1\approx n_2 (a fractional step Δ\Delta of order one percent), he showed the modes collapse into nearly-degenerate linearly polarized LP\mathrm{LP} groups, each described by a single universal normalized propagation constant b(V)b(V). His companion paper that year, "Dispersion in Weakly Guiding Fibers," turned those curves into the working formula for waveguide dispersion, DW∝V d2(Vb)/dV2D_W\propto V\,\mathrm{d}^2(Vb)/\mathrm{d}V^2 — the very term that lets a designer push the zero-dispersion wavelength around. Two years later, with Enrique Marcatili, he built the multimode theory of graded-core fibers, showing that a near-parabolic profile equalizes the mode delays and slashes modal dispersion by the factor Δ/2\Delta/2 that makes graded-index fiber usable. His normalized-parameter language — VV, bb, LP\mathrm{LP} modes — is still the vocabulary every fiber-optics text speaks.

History
1966
— Kao and Hockham argue that purified glass fiber can carry telecommunications, setting the loss target
1971
— Gloge's weakly-guiding analysis yields the LP modes and the b(V) curves behind waveguide dispersion
1973
— Gloge and Marcatili show a near-parabolic graded index equalizes mode delays, cutting modal dispersion
1980
— Mollenauer, Stolen and Gordon observe optical solitons, group-velocity dispersion balancing the Kerr nonlinearity
1982
— dispersion-shifted fiber moves the zero-dispersion wavelength onto the 1.55 μm loss minimum

Going further

Texts. Saleh & Teich, Fundamentals of Photonics, chs. 9–10 · G. P. Agrawal, Fiber-Optic Communication Systems and Nonlinear Fiber Optics (the standard source for DD, β2\beta_2, pulse-broadening) · Hecht, Optics, fiber and dispersion chapters · Loudon, The Quantum Theory of Light (group velocity in media) · Demtröder, Laser Spectroscopy (source spectral width). Key relations and anchor numbers: D=−(2πc/λ2)β2=−(λ/c) d2n/dλ2D=-(2\pi c/\lambda^2)\beta_2=-(\lambda/c)\,\mathrm{d}^2 n/\mathrm{d}\lambda^2; at 1.55 μm1.55\,\mu\mathrm{m}, D=17 ps nm−1km−1⇒β2=−21.7 ps2 km−1D=17\,\mathrm{ps\,nm^{-1}km^{-1}}\Rightarrow\beta_2=-21.7\,\mathrm{ps^2\,km^{-1}}; step-index modal spread Δτ≈Ln1Δ/c=50 ns km−1\Delta\tau\approx Ln_1\Delta/c=50\,\mathrm{ns\,km^{-1}} for n1=1.48n_1=1.48, Δ=0.0101\Delta=0.0101; Gaussian broadening T1=T01+(β2L/T02)2T_1=T_0\sqrt{1+(\beta_2 L/T_0^2)^2} with LD=115 kmL_D=115\,\mathrm{km} for T0=50 psT_0=50\,\mathrm{ps}. Standard values: λ0≈1.31 μm\lambda_0\approx1.31\,\mu\mathrm{m}, loss ≈0.2 dB km−1\approx0.2\,\mathrm{dB\,km^{-1}}, D≈+17 ps nm−1km−1D\approx+17\,\mathrm{ps\,nm^{-1}km^{-1}} at 1.55 μm1.55\,\mu\mathrm{m}, cutoff V=2.405V=2.405.

Builds on: Total Internal Reflection and Optical Guiding · Waveguide Modes and the Wave Equation · Single-Mode and Multimode Operation · Step-Index and Graded-Index Optical Fibers. Leads to: Dispersion Management and Compensation · Polarization-Mode Dispersion and Birefringent Fibers · Nonlinear Optical Fibers and Effects · Fiber Attenuation and Loss Mechanisms.

Ask the tutor: Want to walk through how a coherent receiver's DSP inverts the β2L\beta_2 L phase to undo 100 km of chromatic dispersion in software?

Why is modal dispersion so much bigger than chromatic dispersion?

Because they scale on completely different quantities. Modal dispersion is set by the fractional index step Δ∼0.01\Delta\sim0.01, and Δτ/L=n1Δ/c\Delta\tau/L=n_1\Delta/c reaches tens of ns km−1\mathrm{ns\,km^{-1}}. Chromatic dispersion is the product of a small coefficient DD and a small linewidth σλ\sigma_\lambda, landing in the ps km−1\mathrm{ps\,km^{-1}} range — a thousand times smaller. That is why the industry moved to single-mode fiber for anything long: eliminating modal dispersion alone buys three orders of magnitude of bit-rate–distance product.

If chromatic dispersion just stretches a pulse, why can't I always undo it?

For the linear, first-order part you very nearly can — that is the beauty of β2\beta_2. It imposes a deterministic quadratic phase, so a length of opposite-sign fiber, a chirped grating, or a DSP filter inverts it exactly, which is how coherent systems clean up hundreds of kilometres of accumulated dispersion. What resists undoing is everything beyond first order: the dispersion slope SS that varies across a WDM band, polarization-mode dispersion (random and time-varying), and nonlinear mixing with the Kerr effect — the residuals that ultimately cap a real system.

What actually happens at the zero-dispersion wavelength — is the pulse perfectly preserved?

No. Setting D(λ0)=0D(\lambda_0)=0 cancels only the first-order term. Two things remain. First, the dispersion slope S=dD/dλS=\mathrm{d}D/\mathrm{d}\lambda is nonzero, so different parts of the spectrum see slightly different DD, giving residual broadening ∝∣S∣ L σλ2\propto|S|\,L\,\sigma_\lambda^2. Second, polarization-mode dispersion, indifferent to λ0\lambda_0, still spreads the pulse as L\sqrt{L}. And with dispersion near zero, nonlinear four-wave mixing between WDM channels phase-matches strongly, so dispersion-shifted fiber can be worse for dense multiplexing.

How does a graded-index profile equalize the mode delays?

In a step-index core all rays travel at c/n1c/n_1, so the long-path high-order modes arrive late. A graded index turns nn down away from the axis, so the wide-swinging rays spend their detours in lower-index glass where the local speed c/n(r)c/n(r) is higher. With a near-parabolic profile the speed-up almost exactly compensates the extra path, and all modes arrive together. Gloge and Marcatili showed the optimum shrinks the residual spread from n1Δ/cn_1\Delta/c to about n1Δ2/2cn_1\Delta^2/2c, a factor of Δ/2\Delta/2 — which is why graded-index multimode fiber carries data-centre traffic.

Why quote a bit-rate–distance product instead of just a bit rate?

Because dispersion broadening grows with length, the maximum bit rate always trades off against how far you push it. From στ=∣D∣Lσλ\sigma_\tau=|D|L\sigma_\lambda and B≲1/(4στ)B\lesssim1/(4\sigma_\tau), the length cancels into a constant: B⋅L≲1/(4∣D∣σλ)B\cdot L\lesssim1/(4|D|\sigma_\lambda). Doubling the distance forces you to halve the rate, so the product B⋅LB\cdot L is the invariant that characterizes a fiber-plus-source combination — which is why an engineer's first question about a fiber is not "how fast?" but "how fast times how far?".

Does single-mode fiber literally carry only one mode?

One spatial mode, yes — the fundamental LP01\mathrm{LP}_{01} — provided V<2.405V<2.405, which the core radius and index step are chosen to guarantee. But "one mode" hides a subtlety: that mode has two orthogonal polarization states, and in a real fiber tiny asymmetries (bends, stress, non-circular cores) make them travel at slightly different speeds. The result is polarization-mode dispersion, a residual spreading that scales as L\sqrt{L} rather than LL and cannot be designed away like modal dispersion — the whisper of mode-related spreading that survives single-mode operation.

Before this

Where this leads