Gaussian beam ($I(r,z)=2P{e}^{-2r^2/w(z{)}^2}/(\pi w(z{)}^2)$)

A fundamental Gaussian beam is an ideal beam of light whose transverse wave amplitude has a Gaussian function profile. Here consider a circularly symmetric beam in vacuum, with propagation along $z$, two transverse coordinates $x,y$, and radius $r=\sqrt{x^2+y^2}$ from the optical axis. With the beam waist at $z=0$, total power $P>0$, and vacuum wavelength $\lambda$, its time-averaged light intensity, meaning irradiance in $\mathrm{W}{\mathrm{m}}^{-2}$, isHere $w_0$ is the minimum beam radius, and $z_R$ is the Rayleigh range. Far from the beam waist, $w(z)\sim \theta |z|$ with beam divergence half-angle $\theta =\lambda /(\pi w_0)$. These are the fundamental circular Gaussian beam laws within the scalar paraxial approximation, appropriate when $w_0\gg \lambda$ and the beam divergence is small. The underlying envelope solves the parabolic wave equation exactly; this does not make it an exact vector solution of Maxwell equations. ^[ref].

The radius convention matters. Writing $U$ for the complex wave amplitude envelope, the same $w(z)$ is the wave amplitude $1/e$ radius and the light intensity $1/e^2$ radius:At $r=w(z)$ these ratios are $e^{-1}\approx 0.367879$ and $e^{-2}\approx 0.135335$, respectively. The light intensity $1/e$ radius is instead $w(z)/\sqrt{2}$. None of these radii is a sharp edge: the ideal Gaussian function has nonzero tails at every finite radius. Kogelnik and Li explicitly use the wave amplitude $1/e$ radius in their fundamental-mode treatment. ^[ref].

For a numerical example, take vacuum wavelength $\lambda =632.8\mathrm{nm}$ and beam waist radius $w_0=0.500\mathrm{mm}$. Substitution gives Rayleigh range $z_R\approx 1.24114754\mathrm{m}$ and far-field beam divergence half-angle $\theta \approx 0.402852992\mathrm{mrad}$. At $z=z_R$, the beam radius is $w=\sqrt{2}{w}_0\approx 0.707107\mathrm{mm}$ and the on-axis irradiance is half its waist value. At the same wavelength, halving the beam waist radius to $0.250\mathrm{mm}$ quarters the Rayleigh range to $0.310286885\mathrm{m}$ and doubles the beam divergence half-angle to $0.805705984\mathrm{mrad}$. At fixed power, that smaller beam waist also quadruples the peak irradiance.

The contours in the left panel are $r=\pm w(z)$ in a plane through the optical axis; they mark an light intensity level, not an outer boundary. The right panel compares the exact model radii with the far-field asymptotes $w=\theta z$. The horizontal and vertical axes use different length units, so the drawn slopes are not physical angles.

The normalization expresses conserved power during free propagation. In transverse polar coordinates, an annulus has area $2\pi rdr$, soA proof outline is the substitution $s=2r^2/w(z{)}^2$ in this integral. Spreading reduces peak irradiance as $w^{-2}$ while increasing transverse area as $w^2$, leaving power unchanged. In particular, the disk of radius $w(z)$ contains $1-e^{-2}\approx 0.8646647168$, or $86.4665\%$, of the total power; the local light intensity at its rim is only $13.5335\%$ of the on-axis light intensity. These are different quantities.

The left panel compares normalized wave amplitude magnitude and irradiance; the right panel integrates the transverse area to show enclosed power. For an ideal absorbing circular aperture of radius $a$ centered on the optical axis, the same fraction gives the power passing the aperture at that plane. An aperture of radius $w$ clips about $13.5335\%$ of the incident power; radius $2w$ clips only $e^{-8}\approx 0.000335463$, or $0.0335463\%$. After clipping, extra diffraction means the transmitted field no longer obeys the unmodified ideal Gaussian beam law. Every finite aperture cuts some ideal tail. Kogelnik and Li distinguish their large-aperture beam analysis from their finite-aperture diffraction analysis. ^[ref].

The useful tradeoff is that a small bright spot cannot remain narrow indefinitely. An optical lens can create a smaller beam waist, but diffraction then makes the beam spread faster. The particularly simple relation $w_0\theta =\lambda /\pi$ ties localization directly to beam divergence. A proof outline for the propagation law starts with $U(x,y,0)=U_0{e}^{-(x^2+y^2)/w_0^2}$. The parabolic wave equation separates into two transverse factors, and a Gaussian integral in each coordinate givesHere $q$ is dimensionless and $k$ is the wavenumber; the carrier convention is $e^{i(kz-\omega t)}$, with angular frequency $\omega$. Taking the squared complex modulus gives the stated $w(z)$ and its $w^{-2}$ light intensity scaling. Each transverse factor contributes $q^{-1/2}$, explaining the $q^{-1}$ prefactor here and the different prefactor in the preserved Gaussian beam with one transverse coordinate.

Historically, Gaussian beam theory developed alongside early laser resonator theory and lens-guided transmission. In their 1966 review, Kogelnik and Li credit Boyd and Gordon, and Goubau and Schwering, with early beam treatments; they summarize earlier work rather than claim the whole construction as a new discovery. Their paper brings beam propagation, focusing and finite-aperture diffraction into a common practical account. ^[ref].

Fourier localization of a Gaussian beam expresses the tradeoff between a narrow transverse field and a broad distribution of transverse spatial frequencies. At a waist, takeTwo Gaussian integrals giveReducing $w_0$ therefore broadens the angular spectrum method representation rather than merely shrinking a picture in the same set of propagation directions.^[ref]

Widths must state their weighting. Normalize $|U{|}^2$ and $|\tilde{U}{|}^2$ to define transverse position and transverse-wavenumber distributions. Their one-coordinate standard deviations areThese are optical field-width statements under the stated Fourier convention, not an inference that each ray has simultaneously measured position and angle. Using amplitude widths instead would give different numerical factors.

For a strongly paraxial spectrum, ${\theta }_x\simeq k_x/k$ and ${\sigma }_{{\theta }_x}=1/(k{w}_0)$. The usual far-field Gaussian $1/e^2$ light intensity half-angle is $2{\sigma }_{{\theta }_x}=\lambda /(\pi w_0)$. For original $\lambda =632.8$ nm and $w_0=0.50$ mm, ${\sigma }_x=0.25$ mm and ${\sigma }_{{\theta }_x}\simeq 0.2014$ mrad; the conventional divergence half-angle is about 0.4029 mrad. Halving the waist doubles both angular measures.

An illustrative apparatus places a camera in a plane conjugate to the waist and another measurement at a lens's back focal plane, where transverse spatial frequency maps to position in the paraxial model. Comparing widths tests localization and angular spread, provided magnification, focal length and detector response are accounted for. Finite aperture clipping changes the field and produces extra spectral structure, so the exact Gaussian width product cannot simply be assigned to every clipped beam.

The same spectrum also bounds the paraxial approximation: $k{w}_0\gg 1$ keeps most amplitude at transverse wavenumbers much smaller than $k$. Strong focusing can require vector polarization constraints and a nonparaxial propagation model.

The Gouy phase describes the additional axial phase associated with a focused wave relative to a plane wave. For a fundamental circular Gaussian beam with beam waist at $z=0$ and Rayleigh range $z_R>0$, defineWith the carrier convention $e^{i(kz-\omega t)}$, the on-axis field phase is $kz-\omega t-\psi (z)$ up to a constant. Thus the signed additional phase is $-\psi$ in this convention; defining the positive function $\psi$ does not reverse the physical prediction.^[ref]

The complex field prefactor gives a short derivation. The normalized Gaussian propagation factor isIts modulus describes the on-axis amplitude reduction; its argument is the additional phase. An intensity measurement alone loses this argument, so measuring the beam's width and brightness does not directly measure the Gouy phase.

As $z$ runs from far before to far after the waist, $\psi$ runs continuously from $-\pi /2$ to $+\pi /2$. The total additional phase change has magnitude $\pi$, but is not a discontinuous jump at the waist. For stipulated $w_0=100\mu \mathrm{m}$ and $\lambda =633\mathrm{nm}$, $z_R\approx 49.6\mathrm{mm}$. Between $z=-z_R$ and $z=+z_R$, the additional phase changes by $-\pi /2$; the far-field-to-far-field limit gives $-\pi$.

The number of confined transverse directions matters. A Gaussian beam with one transverse coordinate has prefactor $(1+iz/z_R{)}^{-1/2}$ and half this phase contribution. Simin Feng and Herbert G. Winful's primary analysis relates the effect to transverse confinement and the angular spectrum representation: a localized beam includes transverse wave-vector components, affecting axial propagation.^[ref]

An illustrative Gouy phase interferometry experiment overlaps the focused beam with a coherent reference beam and records the interference phase at different axial planes. Removing the ordinary optical-path contribution and accounting for wavefront curvature isolates the axial phase change. This is a proposed measurement outline, not an experimental result reported by the cited theoretical paper. In an optical cavity, related phase contributions help determine the resonance condition and distinguish transverse modes.

Gouy phase interferometry compares a focused field with a coherent reference to recover the additional Gouy phase. Beam brightness alone discards the phase. An interferometer converts a relative field phase into intensity modulation, but the optical paths and reference wavefront curvature must be modeled.

A published 2023 experiment by Lyubomir Stoyanov and colleagues used a Mach-Zehnder interferometer to study Bessel-Gaussian beams. Its setup included a 532-nanometre laser, beam expansion, two spatial light modulators, focusing lenses and a CCD camera on a rail to follow propagation after the focus. Its Figure 3 also describes a single-lens interferometer arrangement for stronger focusing.^[ref] Structured modes need their own phase model; that experiment should not be silently interpreted as a fundamental-Gaussian measurement.

A separate original protocol illustrates four-step phase retrieval. At one detector point, apply controlled signal-arm phases $\delta =0,\pi /2,\pi ,3\pi /2$ while the beam intensities remain fixed. Write
This phase-shifting interferometry formula removes the unknown common intensity background $C$ and positive modulation amplitude $A$. It returns phase modulo $2\pi$; phase unwrapping requires continuity information between sampled positions. If the interference visibility vanishes, the phase is not observable by this formula.

For a synthetic coherent Gaussian/reference pair, use reference intensity 1 and signal intensity $1/(1+{\zeta }^2)$, with $\zeta =z/z_R$. Then $C=1+1/(1+{\zeta }^2)$, $A=2/\sqrt{1+{\zeta }^2}$ and $\varphi =-\arctan \zeta$ after common path phases are removed. At $\zeta =1$, the four intensities are $(2.5,2.5,0.5,0.5)$, giving $\hat{\varphi }=-\pi /4$. At $\zeta =-1$, the first intensity is still $2.5$, but the quadratures reverse and the retrieved phase is $+\pi /4$.

Curvature bias in Gouy-phase measurements is important away from the axis: the signal adds $k{r}^2/(2R(z))$ to its spatial phase. Fit or remove that term when using multiple pixels, retain calibrated phase steps and account for drift between frames. The recovered phase combines all coherent optical-path contributions until those controls isolate the Gouy phase.

Curvature bias in Gouy-phase measurements occurs when a spatially varying wavefront-curvature wave phase is mistaken for the axial Gouy phase in an interferometric comparison. For a fundamental Gaussian beam with time convention $e^{-i\omega t}$, the signal wave phase containsA curved reference contributes its own quadratic wave phase. After subtracting common axial path terms, the measured relative wave phase isThus a phase-retrieval algorithm can work exactly while the subsequent physical interpretation is biased. The Gaussian confinement wave phase and its arctangent dependence are derived in the primary theoretical treatment.^[ref]

For an original example with a planar reference and stable path wave phase, let $R_s(z)=z[1+(z_R/z{)}^2]$ and $z_R=k{w}_0^2/2$. At $z=z_R$ and fixed detector radius $r=w_0$, the curvature term is 0.5 radians. The true axial contribution is $-\pi /4$, whereas that off-axis pixel gives $-\pi /4+0.5\simeq -0.2854$ radians. At $z=-z_R$ the curvature term reverses. Comparing those two axial planes at this fixed radius gives a wave phase change $-\pi /2+1$, differing from the axial Gouy change by one radian.

The effect is observable even where the signal is appreciable: at $z=z_R$, $r=w_0$ is inside the plane's Gaussian radius $\sqrt{2}{w}_0$. Avoiding only dim pixels does not remove curvature bias.

A useful measurement fits the spatial quadratic wave phase, identifies the beam center and reference curvature, and extrapolates the relative wave phase to the axis. Tilt, decentering and path drift require their own terms; averaging several pixel phases is not automatically equivalent to measuring the on-axis wave phase. Phase unwrapping can be required before fitting a wave phase profile.

The primary 2023 Bessel-Gaussian experiment documents focusing optics and a rail-mounted camera in its interferometer.^[ref] It motivates tracking both propagation position and wavefront geometry, but its structured-mode observations must not be silently reclassified as measurements of the illustrative fundamental-Gaussian model above.

For a circular Gaussian beam with intensity proportional to $e^{-2r^2/w(z{)}^2}$, the displayed $w(z)$ is the intensity $1/e^2$ radius. Its minimum $w_0$ is the beam waist radius, while the larger radius at a cavity mirror is not itself a waist. The Rayleigh range is $z_R=\pi w_0^2/\lambda$.

The complex beam parameter combines a Gaussian beam's wavefront curvature $R_b$ and spot radius $w$: $1/q=1/R_b-i\lambda /(\pi w^2)$. With $q=z+i{z}_R$, the beam waist is $q=i{z}_R$. A real paraxial ray transfer matrix maps it by $q^{\prime }=(Aq+B)/(Cq+D)$. Positive imaginary part selects the physical finite-width branch in this convention.

A real paraxial ray matrix transforms a Gaussian beam parameter by the displayed fractional-linear map. For a lossless round trip with determinant one, a fixed point with positive imaginary part generically requires $|A+D|<2$. Scalar matrices $\pm I$ are degenerate exceptions: their fractional-linear map fixes every physical beam parameter even at the boundary. Finite apertures and resonance phases impose additional physical constraints.

The beam waist is the narrowest cross-section of a freely propagating Gaussian beam. Its radius $w_0$ uses the wave amplitude $1/e$, or light intensity $1/e^2$, convention; its diameter is $2w_0$. Specifying a “spot size” without saying radius or diameter can introduce a factor-of-two error. Within the scalar paraxial approximation, the wavefront at the beam waist is planar. The minimum radius sets both the Rayleigh range and the far-field beam divergence.

The Rayleigh range of a fundamental Gaussian beam in vacuum is the positive axial distance from the beam waist at which its radius grows from $w_0$ to $\sqrt{2}{w}_0$. Thus its contour area doubles and its peak irradiance halves. Within $|z|\le z_R$, the radius stays at most $\sqrt{2}{w}_0$; for $|z|\gg z_R$, beam divergence dominates. At fixed vacuum wavelength, the Rayleigh range scales quadratically with the beam waist radius. The formula here uses the vacuum wavelength, not an unlabelled wavelength in a material.

For the free parabolic wave equation $E_x=i{E}_{zz}/(2k)$, input $E(0,z)=e^{-A{z}^2}$ with $\mathrm{Re}A>0$ propagates asThe branch is continuous from $q(0{)}^{-1/2}=1$. One-dimensional transverse Fresnel propagation or two Gaussian integrals prove the formula. With $A=D^{-2}+ik/(2F)$, the squared envelope is $|q{|}^{-1}\exp [-2z^2/(D^2|q{|}^2)]$. The amplitude power is $-1/2$ because there is only one transverse coordinate. A two-transverse-coordinate Gaussian beam instead has power $-1$ when both coordinates share the same $q$.