Group Velocity and Group Index

Clock a short pulse through 100 m of BK7 glass and compare with the obvious estimate t=nL/ct=nL/c: the pulse arrives 8.2 ns late. The obvious estimate used the wrong index. We take from aq-em-phase-group-velocity the general machinery of wave packets — the two-frequency beat whose envelope moves at Δω/Δk\Delta\omega/\Delta k, the continuum statement vg=dω/dkv_g=\mathrm{d}\omega/\mathrm{d}k (the tangent to the dispersion curve, where the phase velocity of aq-optics-phase-velocity is the chord), the group index, and the superluminal caveat near absorption bands. What this article owns is the material instance that optics actually uses: substituting k=n(ω) ω/ck=n(\omega)\,\omega/c turns the tangent slope into vg=c/ngv_g=c/n_g with ng=n+ω dn/dω=n−λ dn/dλn_g=n+\omega\,\mathrm{d}n/\mathrm{d}\omega=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda — a group index computable from the same Sellmeier data that tabulates nn. At the BK7 d-line: nd=1.51680n_d=1.51680 but dn/dλ=−0.04180 μm−1\mathrm{d}n/\mathrm{d}\lambda=-0.04180\ \mathrm{\mu m^{-1}}, so ng=1.5414n_g=1.5414 — the envelope of a pulse crawls at vg=c/ng=1.945×108v_g=c/n_g=1.945\times10^8 m/s while its crests do 1.9765×1081.9765\times10^8 m/s. In every transparency window dn/dλ<0\mathrm{d}n/\mathrm{d}\lambda<0 (normal dispersion), so ng>nn_g>n universally there: pulses always run late relative to the crest-speed estimate, by 81.9 ps per metre in BK7. That gap is the timing budget of every fibre link, the calibration constant of optical coherence tomography, and the carrier-envelope slip of every femtosecond laser.

Key takeaways
  • A signal is never monochromatic: it is a wave packet, and its envelope — which carries the energy and the information — moves at the group velocity vg=dω/dkv_g=\mathrm{d}\omega/\mathrm{d}k, the tangent slope of the dispersion curve at the carrier (recalled from aq-em-phase-group-velocity).
  • In a medium, k=nω/ck=n\omega/c gives vg=c/ngv_g=c/n_g with the group index ng=n+ω dn/dω=n−λ dn/dλn_g=n+\omega\,\mathrm{d}n/\mathrm{d}\omega=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda (λ\lambda = vacuum wavelength).
  • Normal dispersion means dn/dλ<0\mathrm{d}n/\mathrm{d}\lambda<0, so ng>nn_g>n: the envelope always lags the crests in a transparency window.
  • BK7 at 587.6 nm: nd=1.51680n_d=1.51680, dn/dλ=−0.04180 μm−1\mathrm{d}n/\mathrm{d}\lambda=-0.04180\ \mathrm{\mu m^{-1}} ⇒ ng=1.5414n_g=1.5414, vg=1.945×108v_g=1.945\times10^8 m/s — an 81.9 ps/m delay beyond the naive c/nc/n estimate.
  • Pulse arrival times, fibre latency, and interferometric coherence positions are set by ngn_g; refraction angles and phase matching are set by nn — two different datasheet columns.
  • Inside anomalous-dispersion bands vgv_g can exceed cc or turn negative; there it stops describing transport, and the signal front still moves at exactly cc (owner's caveat, recalled by name).
  • Key idea: the carrier moves at the chord ω/k=c/n\omega/k=c/n, the envelope at the tangent dω/dk=c/ng\mathrm{d}\omega/\mathrm{d}k=c/n_g with ng=n−λ dn/dλn_g=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda — and since dn/dλ<0\mathrm{d}n/\mathrm{d}\lambda<0 in every transparency window, pulses are always slower than crests: time flight with ngn_g, never with nn.
The core insight

A pulse's carrier crests move at the phase velocity c/nc/n, but the envelope — the thing that arrives, carries energy, and can be timed — moves at the group velocity vg=dω/dk=c/ngv_g=\mathrm{d}\omega/\mathrm{d}k=c/n_g, where ng=n+ω dn/dω=n−λ dn/dλn_g=n+\omega\,\mathrm{d}n/\mathrm{d}\omega=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda. Normal dispersion makes ng>nn_g>n, so in glass the envelope always lags: BK7's nd=1.5168n_d=1.5168 hides a group index of 1.54141.5414.

Where this lives today

The group index is a line item on real datasheets. Fibre-optic links budget latency as ngL/cn_g L/c — 4.88 µs per kilometre of silica at 1550 nm — and high-frequency trading and data-centre time synchronization live or die by that number; the next-order coefficient, GVD, is the pulse-broadening budget owned by aq-optics-group-velocity-dispersion-pulse-broadening. Optical coherence tomography and white-light interferometry position their fringe packets at the group-delay match — retinal thickness maps are group-index measurements — and time-of-flight refractometers likewise read ngn_g. And in every femtosecond oscillator vp≠vgv_p\neq v_g slips the carrier under the envelope each round trip — the carrier-envelope offset that frequency combs lock. At the exotic end, dispersion engineering makes vgv_g a free parameter. Electromagnetically induced transparency manufactures a hyper-steep normal-dispersion slope inside a transparency window, reaching ng∼107n_g\sim10^7 — 17 m/s in Hau's 1999 sodium experiment — and "stopped light" quantum memories park the pulse in atomic coherence. Gain-assisted anomalous dispersion delivers the opposite regime (vg=−c/310v_g=-c/310 in the Wang experiment) as a precision test bed for the information-velocity theorems of aq-em-phase-group-velocity.

The physical picture

Add two sinusoids of slightly different frequencies and watch the sum: a fast carrier oscillation inside a slow beat envelope. The carrier crests move at ω/k\omega/k — the phase velocity. The beat maxima sit where the two components agree in phase, and keeping them in agreement as time runs requires matching phase changes, not phases: the envelope therefore moves at the ratio of differences Δω/Δk\Delta\omega/\Delta k, which in the limit is the derivative dω/dk\mathrm{d}\omega/\mathrm{d}k. That two-line argument — owned in general by aq-em-phase-group-velocity and re-run in Derivation A — contains the entire concept: chord versus tangent on the dispersion diagram. A real pulse is a continuum of such pairs, and to first order every feature of its envelope translates rigidly at the common tangent slope vgv_g.

Figure
The punchline — a two-wave beat in BK7 glass at two instants: the envelope (dashed) advances at v_g = c/n_g = 1.945×10⁸ m/s while the carrier crests (solid; one crest marked blue) advance faster at v_p = c/n = 1.9765×10⁸ m/s, because n_g = 1.5414 exceeds n_d = 1.5168
Beat in BK7: envelope at v_g = 1.945×10⁸ m/s, crests at v_p = 1.9765×10⁸ m/s (n_g = 1.5414, n_d = 1.5168) t = 0 t = Δt: the envelope peak has moved 40 px, the marked crest 65 px — crests gain on the envelope envelope: signal and energy, v_g = c/n_g crests slide through the envelope speed gap exaggerated — in BK7, v_p/v_g = n_g/n = 1.016

Where does the material enter? A pulse's spectrum straddles a band of frequencies, and in a dispersive medium each component travels at its own crest speed c/n(ω)c/n(\omega). Blue components lag red ones (normal dispersion), so the point of constructive interference — the envelope peak — drifts backward relative to the crests: envelope slower than carrier, ng>nn_g>n. The group index simply packages the two ways a component can be slow: the nn term (crests slow) plus the ω dn/dω\omega\,\mathrm{d}n/\mathrm{d}\omega term (neighbouring components de-phase). The second figure shows both indices across the visible for BK7: ngn_g rides above nn everywhere in the window, the gap widening toward the blue where dispersion steepens.

Figure
BK7 across the visible, computed from the Schott Sellmeier equation: the group index n_g = n − λ·dn/dλ (upper curve) rides above the phase index n (lower curve); at the d-line, n = 1.5168 and n_g = 1.5414 — a gap of 0.0246 that makes every pulse 81.9 ps/m late relative to the naive c/n estimate
BK7: group index n_g above phase index n across the whole visible window 1.52 1.54 1.56 1.58 400 500 600 700 vacuum wavelength λ (nm) index n_d = 1.5168 n_g = 1.5414 → v_g = 1.945×10⁸ m/s gap n_g − n = 0.0246 n_g(λ) — what a stopwatch sees n(λ) — what a protractor sees normal dispersion: dn/dλ negative everywhere here ⇒ n_g exceeds n across the window

Think of ngn_g as "the index a stopwatch sees" and nn as "the index a protractor sees". Refraction angles, fringe positions and phase matching read nn; arrival times, latency budgets and coherence positions read ngn_g. They never coincide in a real material — a frequency-independent index would violate the dispersion relations that close this chapter (aq-optics-kramers-kronig-relations).

The core equations

The definition first, recalled by name from aq-em-phase-group-velocity, which owns the general wavepacket construction. Expanding the dispersion relation ω(k)\omega(k) about the carrier k0k_0, the envelope of any narrow-band packet translates at the tangent slope — evaluated at the carrier.

vg=dωdk∣k0(1)v_g=\left.\frac{\mathrm{d}\omega}{\mathrm{d}k}\right|_{k_0} \qquad (1)

where: vgv_g — group velocity, the translation speed of the envelope (and, in transparent media, of energy and information); k0k_0, ω0\omega_0 — carrier wavenumber and frequency. The material form follows by differentiating the in-medium dispersion relation k=n(ω) ω/ck=n(\omega)\,\omega/c (from aq-optics-plane-waves-isotropic-media). The derivative of the product defines the group index, the single most-used number in pulse optics.

ng≡cvg=c dkdω=n+ωdndω(2)n_g\equiv\frac{c}{v_g}=c\,\frac{\mathrm{d}k}{\mathrm{d}\omega}=n+\omega\frac{\mathrm{d}n}{\mathrm{d}\omega} \qquad (2)

where: ngn_g — group index; n(ω)n(\omega) — phase index; ω dn/dω\omega\,\mathrm{d}n/\mathrm{d}\omega — the dispersive correction: positive in normal dispersion (dn/dω>0\mathrm{d}n/\mathrm{d}\omega>0), so ng>nn_g>n there. Where dispersion vanishes, ng=nn_g=n and envelope rides with crests. Tabulated glass data come as n(λ)n(\lambda) with λ\lambda the vacuum wavelength, so the practical form trades the frequency derivative for a wavelength derivative via ω=2πc/λ\omega=2\pi c/\lambda (Derivation B).

ng=n−λdndλ(3)n_g=n-\lambda\frac{\mathrm{d}n}{\mathrm{d}\lambda} \qquad (3)

where: λ\lambda — vacuum wavelength; dn/dλ\mathrm{d}n/\mathrm{d}\lambda — the Sellmeier slope, negative throughout every transparency window, which is why the correction −λ dn/dλ-\lambda\,\mathrm{d}n/\mathrm{d}\lambda is positive. BK7 at λ=587.6\lambda=587.6 nm: n=1.51680n=1.51680, dn/dλ=−0.04180 μm−1\mathrm{d}n/\mathrm{d}\lambda=-0.04180\ \mathrm{\mu m^{-1}} ⇒ ng=1.51680+0.58756×0.04180=1.5414n_g=1.51680+0.58756\times0.04180=1.5414. What the group index buys is a clock. The envelope of a pulse crossing a length LL of medium arrives after the group delay — the formula behind every fibre latency budget and time-of-flight refractometer.

τ=Lvg=ng Lc(4)\tau=\frac{L}{v_g}=\frac{n_g\,L}{c} \qquad (4)

where: τ\tau — envelope transit time; LL — geometric path. Per metre of BK7: ng/c=5.141n_g/c=5.141 ns/m versus the naive n/c=5.060n/c=5.060 ns/m — the pulse runs 81.9 ps/m late. Using nn where ngn_g belongs is the classic timing bug. One order deeper sits the variation of the group velocity across the pulse's own spectrum: the group-velocity dispersion β2=d2k/dω2=d(1/vg)/dω\beta_2=\mathrm{d}^2k/\mathrm{d}\omega^2=\mathrm{d}(1/v_g)/\mathrm{d}\omega, which broadens and chirps pulses rather than merely delaying them. It is named here for the map and owned by aq-optics-group-velocity-dispersion-pulse-broadening.

k(ω)≈k0+ω−ω0vg+12 β2 (ω−ω0)2(5)k(\omega)\approx k_0+\frac{\omega-\omega_0}{v_g}+\tfrac12\,\beta_2\,(\omega-\omega_0)^2 \qquad (5)

where: the first-order term translates the envelope at vgv_g (this article); the β2\beta_2 term spreads it (the GVD article). Truncating at first order is legitimate when the pulse is narrow-band and the medium transparent.

Derivation AThe envelope of a two-frequency beat moves at dω/dk

The goal is the group velocity from first principles in its cleanest setting: superpose two monochromatic waves straddling a carrier, factor the sum into carrier × envelope, read off the two speeds, and pass to the limit. The construction is recalled from aq-em-phase-group-velocity, which owns the general wavepacket and stationary-phase machinery; the optical extension — each wavenumber supplied by a refractive index, k=n(ω) ω/ck=n(\omega)\,\omega/c — is what turns the tangent slope into c/ngc/n_g and makes the result computable from glass catalogues.

Step 1. Superpose equal amplitudes at (ω±Δω, k±Δk)(\omega\pm\Delta\omega,\ k\pm\Delta k), the wavenumbers fixed by the medium through k(ω)=n(ω) ω/ck(\omega)=n(\omega)\,\omega/c:

E=E0cos⁡[(k+Δk)z−(ω+Δω)t]+E0cos⁡[(k−Δk)z−(ω−Δω)t].E=E_0\cos\big[(k+\Delta k)z-(\omega+\Delta\omega)t\big]+E_0\cos\big[(k-\Delta k)z-(\omega-\Delta\omega)t\big].

Why: two components are the minimal spectrum with any structure — the simplest object that can carry a timeable feature, since one component alone is featureless (the lesson of the phase-velocity article).

Step 2. Apply the sum-to-product identity:

E=2E0 cos⁡(Δk z−Δω t) cos⁡(kz−ωt).E=2E_0\,\cos(\Delta k\,z-\Delta\omega\,t)\,\cos(kz-\omega t).

Why: the sum factors exactly — a fast carrier cos⁡(kz−ωt)\cos(kz-\omega t) multiplied by a slow envelope cos⁡(Δk z−Δω t)\cos(\Delta k\,z-\Delta\omega\,t); no approximation has been made yet.

Step 3. Read off the two speeds by the constant-phase argument applied to each factor: carrier crests at vp=ω/kv_p=\omega/k, envelope maxima where Δk z−Δω t=const\Delta k\,z-\Delta\omega\,t=\text{const}, i.e. moving at venv=Δω/Δkv_{\text{env}}=\Delta\omega/\Delta k. Why: the beat maximum sits where the two components agree in phase; keeping the agreement as time advances matches phase changes, so the envelope speed is a difference quotient — the finite-difference tangent — while the carrier speed is the chord.

Step 4. Shrink the detuning: vg=lim⁡Δk→0Δω/Δk=dω/dkv_g=\lim_{\Delta k\to0}\Delta\omega/\Delta k=\mathrm{d}\omega/\mathrm{d}k, evaluated at the carrier. A general narrow-band pulse is a Fourier integral of such pairs, and to first order in the expansion of k(ω)k(\omega) every pair's envelope moves at this same slope — the packet translates rigidly at vgv_g (the continuum upgrade is owned by aq-em-phase-group-velocity). With k=nω/ck=n\omega/c, the slope is dk/dω=(n+ω dn/dω)/c\mathrm{d}k/\mathrm{d}\omega=(n+\omega\,\mathrm{d}n/\mathrm{d}\omega)/c, so vg=c/ngv_g=c/n_g. Why: the derivative is the same for every nearby pair only because the expansion was truncated at first order — the second-order spread of tangent slopes is precisely GVD, deferred by name to the pulse-broadening article. *Verified: the sum-to-product factorization and the envelope speed Δω/Δk\Delta\omega/\Delta k symbolically (sympy), plus vg=c/ngv_g=c/n_g with BK7 numbers, in verify_aq-optics-group-velocity-group-index.py.*

Go deeper · Derivation B · 1 minThe group index in frequency and wavelength forms
Derivation BThe group index in frequency and wavelength forms

The goal is the working formula pair ng=n+ω dn/dω=n−λ dn/dλn_g=n+\omega\,\mathrm{d}n/\mathrm{d}\omega=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda, derived by differentiating the dispersion relation and then converting to the vacuum-wavelength variable in which every glass catalogue is written — ending with the real BK7 evaluation that anchors this article. The two forms are one identity connected by the chain rule; sign discipline (which variable, which slope) is where all practical mistakes live, so each sign is tracked explicitly.

Step 1. Differentiate k(ω)=n(ω) ω/ck(\omega)=n(\omega)\,\omega/c by the product rule:

dkdω=1c(n+ωdndω)≡ngc,vg=(dkdω)−1=cng.\frac{\mathrm{d}k}{\mathrm{d}\omega}=\frac{1}{c}\left(n+\omega\frac{\mathrm{d}n}{\mathrm{d}\omega}\right)\equiv\frac{n_g}{c},\qquad v_g=\left(\frac{\mathrm{d}k}{\mathrm{d}\omega}\right)^{-1}=\frac{c}{n_g}.

Why: the tangent slope of ω(k)\omega(k) is the reciprocal of dk/dω\mathrm{d}k/\mathrm{d}\omega; the bracket is defined as the group index because it plays for envelopes exactly the role nn plays for crests.

Step 2. Change variable to vacuum wavelength with ω=2πc/λ\omega=2\pi c/\lambda, so dω/dλ=−2πc/λ2\mathrm{d}\omega/\mathrm{d}\lambda=-2\pi c/\lambda^2 and hence dn/dω=(dn/dλ)(dλ/dω)=−λ22πcdndλ\mathrm{d}n/\mathrm{d}\omega=\big(\mathrm{d}n/\mathrm{d}\lambda\big)\big(\mathrm{d}\lambda/\mathrm{d}\omega\big)=-\frac{\lambda^2}{2\pi c}\frac{\mathrm{d}n}{\mathrm{d}\lambda}. Why: catalogues tabulate n(λ)n(\lambda), not n(ω)n(\omega); the minus sign appears because frequency and wavelength run in opposite directions, and it is the sign that makes normal dispersion come out with ng>nn_g>n in both variables.

Step 3. Substitute into the bracket: ω dndω=2πcλ⋅(−λ22πcdndλ)=−λdndλ\omega\,\dfrac{\mathrm{d}n}{\mathrm{d}\omega}=\dfrac{2\pi c}{\lambda}\cdot\left(-\dfrac{\lambda^2}{2\pi c}\dfrac{\mathrm{d}n}{\mathrm{d}\lambda}\right)=-\lambda\dfrac{\mathrm{d}n}{\mathrm{d}\lambda}, giving ng=n−λ dn/dλn_g=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda. Why: the 2πc2\pi c cancels identically — the wavelength form is exact, not approximate, and both forms are the same number for the same wave.

Step 4. Evaluate for BK7 at the d-line from the Schott Sellmeier fit: n(0.58756 μm)=1.51680n(0.58756\ \mathrm{\mu m})=1.51680 and dn/dλ=−0.04180 μm−1\mathrm{d}n/\mathrm{d}\lambda=-0.04180\ \mathrm{\mu m^{-1}}, so

ng=1.51680−0.58756×(−0.04180)=1.51680+0.02456=1.5414,vg=c1.5414=1.945×108 m/s.n_g=1.51680-0.58756\times(-0.04180)=1.51680+0.02456=1.5414,\qquad v_g=\frac{c}{1.5414}=1.945\times10^8\ \mathrm{m/s}.

Why: the correction term +0.0246+0.0246 is three times the F-to-C index spread of the same glass (0.00806) — a group-index effect larger than the dispersion optical designers already track to five decimals. *Verified: the chain-rule identity n+ω dn/dω=n−λ dn/dλn+\omega\,\mathrm{d}n/\mathrm{d}\omega=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda symbolically (sympy), the analytic Sellmeier derivative against a finite difference, and ng=1.5414n_g=1.5414, vg=1.945×108v_g=1.945\times10^8 m/s in verify_aq-optics-group-velocity-group-index.py.*

Limiting cases

CaseResultMeaning
Non-dispersive mediumdn/dω=0\mathrm{d}n/\mathrm{d}\omega=0 ⇒ ng=nn_g=n, vg=vpv_g=v_pvacuum (exactly); envelope rides with crests
Normal dispersion (glass window)dn/dλ<0\mathrm{d}n/\mathrm{d}\lambda<0 ⇒ ng>nn_g>n, vg<vpv_g<v_pBK7: 1.5414 vs 1.5168 — pulses run 81.9 ps/m late
Plasma / hyperbolic familyvpvg=c2v_p v_g=c^2 ⇒ ng=1/n>1n_g=1/n>1crests above cc, envelope below — pinned product
Anomalous band (on a resonance)ng<1n_g<1, ≤0\le0 possible ⇒ vg>cv_g>c or negativeenvelope reshapes; not a transport speed there
Steep normal dispersion (EIT)ng∼107n_g\sim10^7 ⇒ vg∼30v_g\sim30 m/sslow light: bicycle-speed pulses in cold vapours
High frequency ω→∞\omega\to\inftyn→1n\to1, ng→1n_g\to1every medium turns vacuum-like; fronts at exactly cc

Worked examples

1 · The group index of BK7 at the d-line (confidence builder)

From the Schott Sellmeier data for BK7 — nd=1.51680n_d=1.51680 at λ=587.56\lambda=587.56 nm, dn/dλ=−0.04180 μm−1\mathrm{d}n/\mathrm{d}\lambda=-0.04180\ \mathrm{\mu m^{-1}} — find the group index and group velocity, and compare vgv_g with vpv_p.

Step 1. Wavelength form of the group index (equation 3), with λ\lambda in µm to match the slope's units: ng=n−λ dn/dλ=1.51680−0.58756×(−0.04180)n_g=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda=1.51680-0.58756\times(-0.04180).

Step 2. The correction is 0.58756×0.04180=0.024560.58756\times0.04180=0.02456, positive because the slope is negative: ng=1.51680+0.02456=1.5414n_g=1.51680+0.02456=\mathbf{1.5414}.

Step 3. Speeds: vg=c/ng=299 792 458/1.5414=1.945×108 m/sv_g=c/n_g=299\,792\,458/1.5414=\mathbf{1.945\times10^8\ m/s}, against vp=c/nd=1.9765×108v_p=c/n_d=1.9765\times10^8 m/s. The envelope moves 1.6% slower than the crests (vp/vg=ng/n=1.016v_p/v_g=n_g/n=1.016): each crest is born at the rear of a pulse, slides forward through it, and dies at the front. n_g = 1.5414, v_g = 1.945×10⁸ m/s.

2 · Timing a pulse through 100 m of glass — the naive index fails

A pulsed rangefinder fires through L=100L=100 m of BK7 (round trip through a monolithic delay block). Compute the transit time predicted by the naive estimate τ=nL/c\tau=nL/c and by the correct τ=ngL/c\tau=n_g L/c, and the resulting timing error.

Step 1. Naive (phase-index) estimate: τp=ndL/c=1.51680×100/(2.99792458×108)=505.95 ns\tau_p=n_d L/c=1.51680\times100/(2.99792458\times10^8)=\mathbf{505.95\ ns}.

Step 2. Correct (group-index) transit: τg=ngL/c=1.54136×100/(2.99792458×108)=514.14 ns\tau_g=n_g L/c=1.54136\times100/(2.99792458\times10^8)=\mathbf{514.14\ ns}.

Step 3. The pulse is late by Δτ=τg−τp=(ng−n)L/c=0.02456×100/c=8.19 ns\Delta\tau=\tau_g-\tau_p=(n_g-n)L/c=0.02456\times100/c=\mathbf{8.19\ ns} — i.e. 81.9 ps for every metre of glass. Read as distance, a rangefinder calibrated with nn would misplace its target by c Δτ/(2n)≈0.8c\,\Delta\tau/(2n)\approx0.8 m over this path. The same accounting at telecom wavelengths (silica: n=1.4440n=1.4440, ng=1.4626n_g=1.4626 at 1550 nm) sets the latency of every fibre span at 4.879 µs/km — always from ngn_g. Δτ = 8.19 ns over 100 m; group delay 5.141 ns per metre of BK7.

Physics of the portrait

The group velocity is the youngest of the classical wave speeds — it had to wait for physicists to admit that "the" velocity of light in a medium was ill-posed, and it was born not in optics but in the mathematics of water waves.

Origins · 1 mingroup velocity and how this was discovered
GV
group velocity
the envelope speed v_g = dω/dk = c/n_g — the speed a pulse and its information actually travel · electromagnetic optics of media

The group velocity entered physics sideways. William Rowan Hamilton, reworking optics and mechanics into one variational structure, noted in 1839 that in a dispersive medium the velocity of a group of waves differs from the velocity of its individual crests — an observation that sat unused for decades (Stokes revived it as an examination question in 1876). It was Rayleigh's Theory of Sound (1877) that made vg=dω/dkv_g=\mathrm{d}\omega/\mathrm{d}k standard equipment, tying it to the energy transport of a wave train and supplying the practical forms used ever since; the deep-water case — crests at twice the group speed, visibly overtaking their own packet — made the abstraction undeniable. Optics inherited the concept wholesale and gave it its sternest test: near absorption lines the group velocity exceeds cc or turns negative, and Sommerfeld and Brillouin's 1914 dissection of a signal entering a Lorentz medium established what survives there (nothing kinematic — the front arrives at exactly cc). The modern era turned vgv_g from a diagnosis into a design variable: fibre engineers budget latency with the group index, ultrafast laboratories manage the carrier-envelope slip vp≠vgv_p\neq v_g cycle by cycle, and electromagnetically induced transparency has parked ngn_g at 10710^7 — light at bicycle speed — en route to stopping it altogether.

History
1839
— Hamilton distinguishes group from phase propagation in dispersive media
1877
— Rayleigh's Theory of Sound systematizes v_g = dω/dk and links it to energy transport
1914
— Sommerfeld and Brillouin dissect signal arrival in a Lorentz medium: precursors, signal velocity, front at exactly c
1999
— Hau slows light to 17 m/s in a cold sodium cloud via electromagnetically induced transparency
2000
— Wang, Kuzmich and Dogariu measure v_g = −c/310 in gain-assisted caesium; the front still arrives at c

Going further

Texts. Hecht, Optics, 5th ed., §7.2.2 (Group Velocity) · Foot, Atomic Physics, §7.4 · Loudon, The Quantum Theory of Light, 3rd ed., Ch. 1 · Saleh & Teich, Fundamentals of Photonics, 2nd ed., §5.6 · Demtröder, Laser Spectroscopy, Vol. 1, Ch. 3 · Jackson, Classical Electrodynamics, 3rd ed., §7.8 (Superposition and Group Velocity)

Provenance. Definitions vp=ω/kv_p=\omega/k, vg=dω/dkv_g=\mathrm{d}\omega/\mathrm{d}k and the group-index identities ng=n+ω dn/dω=n−λ dn/dλn_g=n+\omega\,\mathrm{d}n/\mathrm{d}\omega=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda follow directly from k=nω/ck=n\omega/c (Hecht §7.2.2; Saleh & Teich §5.6). A verifier should reproduce the two-wave beat derivation and recompute the glass example from standard Sellmeier coefficients; the anomalous-dispersion inequalities and the front-velocity result (vfront=cv_{\text{front}}=c) are standard (Loudon Ch. 1; Jackson §§7.8, 7.11), recalled here from aq-em-phase-group-velocity rather than re-derived. Addendum for this edition: the previous edition's fused-silica anchor (n=1.458n=1.458, dn/dλ≈−1.7×10−5 nm−1\mathrm{d}n/\mathrm{d}\lambda\approx-1.7\times10^{-5}\ \mathrm{nm^{-1}}, ng≈1.468n_g\approx1.468, 4.89 ns/m) is not reproducible from the Malitson Sellmeier fit, which gives dn/dλ(589 nm)=−3.50×10−5 nm−1\mathrm{d}n/\mathrm{d}\lambda(589\ \mathrm{nm})=-3.50\times10^{-5}\ \mathrm{nm^{-1}} and ng=1.4790n_g=1.4790 (delay 4.933 ns/m) — the slope was understated by half. This edition anchors on BK7 instead (nd=1.51680n_d=1.51680, dn/dλ=−0.04180 μm−1\mathrm{d}n/\mathrm{d}\lambda=-0.04180\ \mathrm{\mu m^{-1}}, ng=1.5414n_g=1.5414, vg=1.945×108v_g=1.945\times10^8 m/s), with both glasses recomputed analytically and by finite difference in verify_aq-optics-group-velocity-group-index.py.

Builds on: Phase Velocity · Plane Waves in Homogeneous Isotropic Media · The Refractive Index · Dispersion: Normal and Anomalous Leads to: Sellmeier and Cauchy Dispersion Equations · Kramers–Kronig Relations · The Lorentz Oscillator Model · Absorption Coefficient and Beer–Lambert Law

Why is the envelope speed a derivative rather than the ratio ω/k?

Because an envelope is an agreement between components, and agreements move at the rate their differences dictate. The beat maximum of two superposed waves sits where their phases coincide; a moment later each phase has advanced by its own ω dt\omega\,\mathrm{d}t, and the coincidence point must shift by whatever dz\mathrm{d}z keeps the phase difference zero: Δk dz=Δω dt\Delta k\,\mathrm{d}z=\Delta\omega\,\mathrm{d}t, i.e. speed Δω/Δk\Delta\omega/\Delta k. The ratio ω/k\omega/k answers a different question — how fast one component's own phase pattern slides.

In n_g = n − λ·dn/dλ, do I use the vacuum wavelength or the wavelength inside the glass?

The vacuum wavelength — always. The formula was derived by substituting ω=2πc/λ\omega=2\pi c/\lambda, which is the vacuum relation, and Sellmeier coefficients are fitted with λ\lambda meaning vacuum (or air) wavelength. Feeding in the medium wavelength λ0/n\lambda_0/n double-counts the index and corrupts the correction term by a factor of order nn. The confusion arises because Rayleigh's other form, vg=vp−λmed dvp/dλmedv_g=v_p-\lambda_{\text{med}}\,\mathrm{d}v_p/\mathrm{d}\lambda_{\text{med}}, genuinely uses the in-medium wavelength — a different formula in different variables. Safe practice: work from catalogue n(λvac)n(\lambda_{\text{vac}}) data and convert nothing.

Can the group index be measured directly, or only computed from Sellmeier fits?

Directly, and routinely. A white-light or low-coherence interferometer locates the position where the envelope of the fringe packet is centred — the group-delay match — so inserting a sample of thickness LL shifts the packet by (ng−1)L(n_g-1)L, giving ngn_g without any differentiation of data. Time-of-flight measurement of a picosecond pulse through a known length does the same cruder and faster. Comparing the interferometric ngn_g with the refractometric nn for the same sample is a classic laboratory exercise: for BK7 at the d-line you expect 1.54141.5414 versus 1.51681.5168 — an extra 24.6 µm of optical path per millimetre of glass, about 42 fringes, which the apparatus resolves easily.

Why is the group index larger than the phase index in every transparency window, for every material?

Because transparency windows sit between resonances, and there the index always climbs with frequency. The Lorentz-oscillator response (owned by aq-optics-lorentz-oscillator-model) makes dn/dω>0\mathrm{d}n/\mathrm{d}\omega>0 wherever absorption is negligible — normal dispersion is not a habit of glass but a theorem about lossless response functions, formalized by the Kramers–Kronig relations that close this chapter. Positive dn/dω\mathrm{d}n/\mathrm{d}\omega feeds straight into ng=n+ω dn/dω>nn_g=n+\omega\,\mathrm{d}n/\mathrm{d}\omega>n. The inequality fails only where transparency fails: inside an absorption band the slope turns negative and ngn_g dives below nn — even below zero.

If v_g can exceed c near a resonance, what actually limits how fast a message crosses the medium?

The front velocity — and it is always exactly cc. A genuine message must contain a surprise — a switch-on, a bit flip — and such a feature has spectrum out to arbitrarily high frequency, where every medium's electrons are too sluggish to respond: n→1n\to1 and the leading edge propagates as in vacuum. Sommerfeld and Brillouin made this exact in 1914 (recalled from aq-em-phase-group-velocity): ahead of z=ctz=ct the field is identically zero, whatever ngn_g is doing. Measured superluminal and negative group velocities — Wang's −c/310-c/310 — are real envelope reshaping of smooth, band-limited pulses whose "information" already resided in their leading tails; put a sharp front on the pulse and it arrives at cc, every time.

How do I get n_g from Sellmeier coefficients in practice?

Differentiate analytically, never numerically from tabulated points. With n2(λ)=1+∑jBjλ2/(λ2−Cj)n^2(\lambda)=1+\sum_j B_j\lambda^2/(\lambda^2-C_j), one differentiation gives dn/dλ=−λn∑jBjCj(λ2−Cj)2\mathrm{d}n/\mathrm{d}\lambda=-\frac{\lambda}{n}\sum_j\frac{B_j C_j}{(\lambda^2-C_j)^2} — manifestly negative between resonances, which re-proves ng>nn_g>n there. Then ng=n−λ dn/dλn_g=n-\lambda\,\mathrm{d}n/\mathrm{d}\lambda with λ\lambda in the same units as Cj\sqrt{C_j} (µm for glass catalogues). For BK7 this yields −0.04180 μm−1-0.04180\ \mathrm{\mu m^{-1}} and ng=1.5414n_g=1.5414 at the d-line; for fused silica, −0.0350 μm−1-0.0350\ \mathrm{\mu m^{-1}} and ng=1.4790n_g=1.4790 at 589 nm — note the widely copied value 1.468 comes from an understated slope and does not survive the differentiation.

Before this

Where this leads