How to calculate chromatic dispersion with the Sellmeier equation

Compute refractive index n(λ) of optical glasses, silica, sapphire, Si, Ge with the Sellmeier equation (1872): n²(λ) = 1 + ΣBᵢλ²/(λ²-Cᵢ). Abbe number, group index, GVD for fiber optics, EUV lithography, ultrafast lasers — base of Schott, Corning, Carl Zeiss SMT, ASML, IBM, Intel.

Calculation parameters
Was this calculator useful?

How to calculate chromatic dispersion with the Sellmeier equation lets you precisely determine the refractive index n of an optical material as a function of wavelength λ. The formula was derived in 1872 by Wilhelm Sellmeier (University of Göttingen) as an empirical improvement over Cauchy (1836) and Lorentz, based on the classical theory of resonant oscillators. The standard three-term form is n²(λ) = 1 + B₁λ²/(λ²-C₁) + B₂λ²/(λ²-C₂) + B₃λ²/(λ²-C₃), with λ in micrometers (µm). For Schott BK7 glass, B₁=1.03961212, B₂=0.231792344, B₃=1.01046945, C₁=0.00600069867, C₂=0.0200179144, C₃=103.560653. It underpins the design of camera lenses, fiber optics, EUV lithography systems, advanced microscopy, and lasers — used by Carl Zeiss SMT, Schott Mainz, Corning, OFS, ASML, IBM, Intel, and academic labs at MIT, Stanford, Caltech, Bell Labs.

1. The Sellmeier equation and Abbe number

The Sellmeier equation describes anomalous and normal dispersion of a dielectric material via sums of resonant terms: each Cᵢ is the squared resonance wavelength (in µm²) and Bᵢ the oscillator strength. The Abbe number ν_d = (n_d - 1)/(n_F - n_C) — where n_d is at 587.6 nm (Fraunhofer D), n_F at 486.1 nm (F) and n_C at 656.3 nm (C) — quantifies dispersion: ν_d > 50 indicates "crown" glasses (BK7 ν_d=64.17, low dispersion), ν_d < 50 indicates "flint" glasses (SF11 ν_d=25.76, high dispersion). Group velocity dispersion (GVD) ∝ d²n/dλ² is critical for ultrafast laser pulse compression (CPA, Mourou-Strickland Nobel 2018) and for fiber-optic communications.

2. Sellmeier coefficients for typical materials

  • Fused silica SiO₂ (Malitson 1965) — B₁=0.6961663, B₂=0.4079426, B₃=0.8974794, C₁=0.0684043, C₂=0.1162414, C₃=9.896161 µm². Range 200-3500 nm. Basis of Corning telecom fiber.
  • BK7 borosilicate (Schott) — Carl Zeiss SMT standard glass, range 350 nm-2000 nm.
  • SF11 heavy flint (Schott) — high dispersion for achromats, ν_d=25.76.
  • Sapphire Al₂O₃ (Tropf 1997) — high hardness, range 300 nm-5000 nm.
  • Silicon Si (Tatian 1984) — B₁=10.6684293, B₂=0.003043475, B₃=1.54133408. Range 1100-7000 nm, basis of integrated photonics.
  • Germanium Ge — IR transparency 2-15 µm, FLIR thermal imaging.
  • ZnSe — CO₂ laser windows (10.6 µm), range 600 nm-21 µm.
  • CVD diamond — extreme IR optics, range 230 nm-100 µm, n=2.42 at 590 nm.

3. Validity range and corrections

The Sellmeier equation is valid for isotropic transparent materials in the spectral range of applicability (typically 200 nm - 5 µm for silica, 400 nm - 2 µm for BK7). Near absorption resonances (UV for most glasses, IR for molecular bonds), Sellmeier diverges and the full Lorentz model with absorption is needed. For birefringent crystals (calcite, quartz, BBO, KTP), separate Sellmeier equations are needed for ordinary n_o and extraordinary n_e rays. For doped glasses or laser crystals (Nd:YAG, Yb:KGW, Ti:sapphire), Sellmeier requires correction for active-ion absorptions. Group dispersion and group velocity dispersion (GVD) require dn/dλ and d²n/dλ², critical for ultrashort laser pulses at LCLS Stanford, European XFEL Hamburg, ELI Beamlines. For temperature variations, dn/dT (~10⁻⁶ K⁻¹ typical) is needed. Codes Zemax, Code V, FRED, Lumerical FDTD implement all corrections.

4. Worked numerical examples

Example 1 — Fused silica at λ=1550 nm (telecom). λ=1.55 µm; λ²=2.4025. n²=1 + 0.6961663·2.4025/(2.4025-0.0684043²) + 0.4079426·2.4025/(2.4025-0.1162414²) + 0.8974794·2.4025/(2.4025-9.896161). Result: n_silica(1550) ≈ 1.4440 — matches Corning SMF-28 single-mode fiber tables.

Example 2 — BK7 at λ=587.56 nm (Na D line). Standard Schott coefficients yield n_d ≈ 1.5168 — IUPAC reference for standard glass. Abbe number ν_d = (1.5168-1)/(1.5224-1.5143) ≈ 63.8 (close to tabulated 64.17).

Example 3 — Silicon at λ=1310 nm (telecom O-band). Tatian 1984 coefficients yield n_Si(1310) ≈ 3.5037. For silicon photonics waveguides at IBM, Intel, IMEC Leuven, CEA-Leti this is the baseline for designing couplers, MZI modulators, and Ge/Si photodetectors.

5. Practical applications of Sellmeier

  • Optical lens design: Carl Zeiss SMT, Schott Mainz, Hoya, Nikon, Canon, Pentax use Sellmeier for every element of achromatic, apochromatic, microscope, and EUV lenses
  • Telecom fiber optics: AT&T, Verizon, Comcast, Lumen, Vodafone backbone Internet at 1310/1550 nm with chromatic dispersion calculated via Sellmeier on doped silica
  • Integrated photonics: silicon waveguides at 1310/1550 nm designed via Sellmeier at CEA-Leti, IMEC Leuven, IBM Almaden, Intel Hillsboro
  • DUV/EUV lithography: ASML EUV scanners use Sellmeier for projection optics design (Mo/Si multilayer mirrors)
  • Ultrafast lasers: LCLS Stanford, European XFEL, BELLA Berkeley, ELI Beamlines use Sellmeier for dispersion management with gratings, prisms, dispersive mirrors in fs and attosecond lasers
  • Confocal microscopy: Leica, Olympus, Nikon, Zeiss immersion lenses use Sellmeier for chromatic aberration correction

6. Comparison with other dispersion formulas

The Cauchy equation (1836) n(λ) = A + B/λ² + C/λ⁴ is a simplified approximation valid only in the visible, replaced by Sellmeier for extended range. The Conrady equation n(λ) = n₀ + A/λ + B/λ^3.5 is occasionally used for optical glasses. The Schott equation n²(λ) = a₀ + a₁λ² + a₂λ⁻² + a₃λ⁻⁴ + a₄λ⁻⁶ + a₅λ⁻⁸ is internally used in Schott Glass Catalog. The Herzberger formula for IR glasses is a 6-parameter variant. For laser crystals, NIST and Brookhaven Standard use semi-empirical formulas. Group dispersion requires d/dλ of Sellmeier. All algorithms are implemented in Zemax OpticStudio, Code V, FRED, Lumerical FDTD for industrial optical design.

7. Special cases and advanced materials

For birefringent crystals (calcite CaCO₃, crystalline quartz, BBO, KTP, LiNbO₃), separate Sellmeier equations are used for ordinary n_o and extraordinary n_e rays — fundamental for polarizers, λ/4 plates, and second-harmonic generation (SHG) at LCLS, European XFEL, ELI Beamlines. For nonlinear phase-matching crystals, the propagation angle is computed from Sellmeier to optimize conversion efficiency. For photonic crystal fibers (PCF) and hollow-core fibers, dispersive properties combine material Sellmeier with geometric structure, developed at OFS Specialty Photonics, Corning, Heriot-Watt, IPHT Jena. For doped materials (Yb:KGW, Nd:YAG, Ti:sapphire), Sellmeier requires correction for laser-active ions. Spatial dispersion in metamaterials and topological photonics requires k-dependent formalism beyond classical Sellmeier.

Updated July 2026

Frequently Asked Questions

What is the Sellmeier equation and why is it important?

The Sellmeier equation (1872) describes the refractive index n(λ) of an optical material via sums of resonant oscillators: n²(λ) = 1 + ΣBᵢλ²/(λ²-Cᵢ). Wilhelm Sellmeier derived it at the University of Göttingen as an improvement over Cauchy (1836). It is the basis of lens, fiber-optic, and integrated-photonics design used by Carl Zeiss, Schott, Hoya, Corning.

How do you compute n(λ) with Sellmeier?

n²(λ) = 1 + B₁λ²/(λ²-C₁) + B₂λ²/(λ²-C₂) + B₃λ²/(λ²-C₃), with λ in µm. For BK7 (B₁=1.0396, B₂=0.2318, B₃=1.0105, C₁=0.00600, C₂=0.02002, C₃=103.56), at λ=0.5876 µm we get n≈1.5168. For fused silica at λ=1.55 µm: n≈1.444 (telecom standard). Coefficients tabulated in Schott, Hoya, and Ohara catalogs.

What is the Abbe number?

The Abbe number ν_d = (n_d-1)/(n_F-n_C) quantifies a glass's dispersion: n_d at 587.6 nm (Fraunhofer D), n_F at 486.1 nm (F), n_C at 656.3 nm (C). Crown glasses have ν_d > 50 (BK7=64.17, low dispersion); flint glasses ν_d < 50 (SF11=25.76, high dispersion). Combining crown and flint yields achromatic doublets for microscopy and photography.

What are the Sellmeier coefficients for fused silica?

For fused silica SiO₂ (Malitson 1965): B₁=0.6961663, B₂=0.4079426, B₃=0.8974794, C₁=0.0684043², C₂=0.1162414², C₃=9.896161². Valid 200-3500 nm. n_silica(589 nm)=1.458, n(1310 nm)=1.447, n(1550 nm)=1.444. Basis for Corning, OFS, Prysmian fiber design and Carl Zeiss SMT UV fused-silica lenses.

What are the limits of Sellmeier?

Valid for isotropic transparent materials in the applicability range (~200 nm-5 µm for silica). Near absorption resonances (UV or IR), Sellmeier diverges and full Lorentz model with absorption is needed. Birefringent crystals (calcite, BBO) need separate Sellmeier for ordinary/extraordinary rays. Doped glasses (Nd:YAG) need correction. Ultrashort pulses need GVD = d²n/dλ².

How does Sellmeier apply to telecom fiber optics?

Singlemode silica fibers operate at 1310 nm (O-band, zero dispersion) and 1550 nm (C-band, minimum 0.2 dB/km attenuation). Sellmeier computes n(λ) and chromatic dispersion D = -λ/c · d²n/dλ², which at 1550 nm is ~17 ps/(nm·km) for undoped silica. Compensation with DCF (Dispersion Compensating Fiber) and DSF (Dispersion-Shifted Fiber) by Corning and OFS is designed via Sellmeier.

How does Sellmeier apply to integrated photonics?

Crystalline silicon (Si SOI) waveguides at 1310/1550 nm have n_core=3.50, n_cladding=1.44 (silica). Sellmeier provides n(λ) for designing MZI modulators, directional couplers, ring resonators, WDM demultiplexers. Companies IBM, Intel, ST Microelectronics, IMEC Leuven, CEA-Leti develop Si photonic chips for data centers, biosensors, and quantum computing.

Comments (0)

No comments yet. Be the first!