How to calculate Rayleigh range and waist of a Gaussian laser beam lets you determine the propagation geometry of a fundamental TEM₀₀ laser beam in the paraxial regime. The formulas were derived in 1966 by Hermann Kogelnik and Tingye Li at Bell Labs (Appl. Opt. 5:1550) and codified by Anthony Siegman in his book "Lasers" (1986, Stanford). The Rayleigh range is z_R = π·w₀²/λ (or with beam quality M²: z_R = π·w₀²/(λ·M²)), the axial distance over which the beam cross-section doubles. The waist radius evolves as w(z) = w₀·√(1 + (z/z_R)²), and the far-field divergence is θ = λ/(π·w₀) (half-angle). It is the foundation for laser focusing in industrial cutting (Trumpf, IPG, Coherent), confocal microscopy (Leica, Zeiss, Nikon, Olympus), optical traps (Nobel Prize 2018, Arthur Ashkin), free-space optical communication (NASA Lunar Laser Communication, Starlink ISL), and the design of 3-4 km arm cavities at LIGO Hanford-Livingston (Caltech-MIT) and Advanced Virgo Cascina (INFN-EGO).

1. Fundamental Gaussian-beam formulas

A TEM₀₀ Gaussian beam propagating along z has intensity I(r,z) = I₀·(w₀/w(z))²·exp(-2r²/w²(z)). Key parameters: w₀ (waist radius, defined at 1/e² intensity, 86.5% power), z_R (Rayleigh range), w(z) = w₀·√(1 + (z/z_R)²), R(z) = z·(1 + (z_R/z)²) (wavefront curvature), ζ(z) = arctan(z/z_R) (Gouy phase, π total over propagation). For real beams with M² > 1, formulas scale by M². The BPP (Beam Parameter Product) = w₀·θ = M²·λ/π is propagation-invariant and defines the laser's "brightness."

2. Parameters for typical lasers

• HeNe 632.8 nm, w₀ = 0.5 mm, M² = 1.1 → z_R = 1.12 m, θ = 0.40 mrad. Interferometry, alignment.
• Nd:YAG 1064 nm, w₀ = 100 µm focused, M² = 1.3 → z_R = 22.7 mm, θ = 4.4 mrad. Industrial metal cutting Trumpf, IPG.
• CO₂ 10.6 µm, w₀ = 50 µm focused, M² = 1.2 → z_R = 0.62 mm, θ = 56.2 mrad.
• ArF excimer 193 nm, w₀ = 1 mm, M² = 50 → z_R = 32 cm, θ = 3.1 mrad. ASML DUV lithography.
• Telecom 1550 nm singlemode silica fiber, w₀ ≈ 5.2 µm mode-field, divergence in air ≈ 95 mrad — basis of AT&T, Verizon networks.
• Ti:Sapphire 800 nm CPA fs, w₀ = 50 µm, M² = 1.1 → z_R = 9.8 mm. Ultrashort pulse lasers at LCLS, ELI Beamlines.
• LIGO 1064 nm, w₀ ≈ 5 cm at end mirror over 4 km arm, z_R ≈ 7.4 km: designed to minimize diffraction loss.

3. Limits and validity regimes

Classical formulas are valid in the paraxial regime (θ < 30°), with linear polarization and single-mode TEM₀₀. For high-NA focusing (NA > 0.5), the Richards-Wolf formalism (1959) including vector and evanescent effects is required — fundamental for STED microscopy (Hell, Moerner, Betzig, Nobel 2014). For multimode beams (M² > 1), w(z) and θ scale by √M². For non-Gaussian beams (super-Gaussian, top-hat, Bessel-Gauss, vortex OAM beams) specialized treatments apply. For nonlinear propagation (self-focusing, Kerr lens), the beam collapses or modulates — relevant at PW lasers like LCLS at SLAC, NIF Livermore, Apollon France, ELI-Beamlines Czech Republic. For slit and aperture diffraction, w₀ sets the Abbe resolution limit.

4. Worked numerical examples

Example 1 — HeNe 632.8 nm focused at w₀ = 100 µm. z_R = π·(10⁻⁴)²/(632.8×10⁻⁹) = 49.6 mm. Divergence θ = 632.8e-9/(π·10⁻⁴) = 2.01 mrad. At z = 1 m from waist, w(1m) = 100·√(1 + (1/0.0496)²) = 2.02 mm.

Example 2 — Trumpf Nd:YAG 4 kW for stainless cutting, w₀ = 50 µm with M² = 1.3. z_R = π·(50×10⁻⁶)²/(1064×10⁻⁹·1.3) = 5.68 mm. Peak intensity I₀ = 2P/(π·w₀²) = 1.02×10¹¹ W/m² = 1 GW/cm² — sufficient for metal ablation.

Example 3 — Advanced LIGO, 1064 nm laser on 4 km arm end mirror. w₀ ≈ 5 cm at beam splitter, z_R = π·(5×10⁻²)²/1064×10⁻⁹ = 7,378 m. Over 4 km, w(4 km) = 5·√(1 + (4000/7378)²) ≈ 5.9 cm — manageable on 34 cm diameter test masses, fundamental for gravitational-wave detection sensitivity (GW150914 event 2015, Nobel 2017).

5. Practical applications

• Industrial lasers: Trumpf (Germany), IPG (USA), Coherent, Lumibird, Lockheed Martin design cutting, welding, marking heads using Gaussian beam optics
• Advanced microscopy: confocal (Leica, Zeiss LSM, Nikon, Olympus FV4000), STED, two-photon, light-sheet at MIT, Stanford, Caltech, MPI Munich
• Optical traps: tweezers (Ashkin Nobel 2018) use focused Gaussian beams to manipulate cells and nanoparticles
• Optical telecommunications: AT&T, Verizon, Comcast operate single-mode fibers with Gaussian mode-field
• Gravitational-wave detection: Advanced LIGO Hanford-Livingston (Caltech-MIT), Advanced Virgo Cascina (INFN-EGO), KAGRA Japan use Fabry-Perot cavities with Gaussian beams over 3-4 km arms
• Free-space optical communication: NASA LCRD/LLCD lunar laser, Starlink ISL inter-satellite, Mynaric LCT

6. Comparison with other formalisms

The Richards-Wolf formalism (1959) generalizes for high NA and includes vector polarization effects. Hermite-Gauss (TEM_mn) describes higher-order modes with nodes. Laguerre-Gauss (LG_pl) includes vortex beams with orbital angular momentum (lh̄). Bessel-Gauss describes non-diffractive beams with rotational Bessel profile. For ultrashort pulses, space-time coupling and chromatic dispersion enter. Codes Zemax OpticStudio, GLAD, COMSOL Multiphysics, Lumerical FDTD, ASAP implement all formalisms.

7. Special cases and advanced regimes

For laser cavities, Kogelnik-Li formulas determine the stable TEM₀₀ mode via |g₁·g₂|<1 (g_i = 1-L/R_i). For Fabry-Perot cavities at LIGO, finesse F = π·√(R_1·R_2)/(1-R_1·R_2) reaches 450 with 99.995% mirror reflectivity. For STED microscopy (Hell-Moerner-Betzig, Nobel 2014), the depletion vortex beam achieves sub-diffraction resolution ~30 nm. For multipath kernels and Bessel non-diffractive beams (Durnin 1987), classical formulas don't apply. OAM vortex beams (Allen 1992) carry orbital angular momentum, used in quantum communication. Adaptive optics at Keck, VLT-ESO, Subaru telescopes correct the beam phase to recover the diffraction-limited Gaussian profile.

Updated July 2026