M² Beam Quality Measurement by the Knife-Edge Method

Scope

This article describes the procedure for measuring the M² beam quality parameter of a laser beam using the knife-edge method, conforming to ISO 11146-1 (Test methods for laser beam widths, divergence angles, and beam propagation ratios). Coverage includes the D4σ second-moment beam width definition, the multi-position hyperbolic fit, and the dominant sources of measurement error. Scanning-slit and CCD-based variants of the same fundamental measurement are discussed for context but not in procedural detail.

Definition

The beam quality parameter $M^2$ (also called the beam propagation ratio) compares the divergence of a real beam to that of an ideal Gaussian beam of the same waist size. For a Gaussian beam of waist radius $w_0$ propagating in the $z$ direction, the beam radius evolves as:

$$w(z) \;=\; w_0 \sqrt{1 + \left(\frac{z - z_0}{z_R}\right)^2}, \qquad z_R \;=\; \frac{\pi w_0^2}{\lambda},$$

where $z_0$ is the position of the waist and $z_R$ is the Rayleigh range. The far-field divergence half-angle is $\theta = \lambda / (\pi w_0)$.

For a real (non-Gaussian) beam, the same hyperbolic functional form applies, but the divergence is scaled by $M^2 \geq 1$:

$$w(z) \;=\; w_0 \sqrt{1 + M^4 \left(\frac{z - z_0}{z_R^G}\right)^2}, \qquad z_R^G \;=\; \frac{\pi w_0^2}{\lambda},$$

equivalently expressed in terms of beam waist and divergence:

$$M^2 \;=\; \frac{\pi w_0 \theta}{\lambda}.$$

$M^2 = 1$ corresponds to an ideal diffraction-limited Gaussian beam. Higher values indicate the beam diverges faster than the corresponding ideal Gaussian, requiring tighter optics or larger apertures to deliver the same beam quality at distance.

Typical $M^2$ values:

For asymmetric beams (such as edge-emitter outputs), $M^2$ is reported separately for the two principal axes, $M_x^2$ and $M_y^2$.

Beam width definition

ISO 11146 specifies the second-moment beam width — the D4σ width — as the standard width for $M^2$ extraction:

$$w^2 \;=\; 4 \sigma^2 \;=\; \frac{4 \iint (x - x_c)^2 \, I(x, y) \, dx \, dy}{\iint I(x, y) \, dx \, dy},$$

where $I(x, y)$ is the intensity distribution and $x_c$ is the centroid. For a Gaussian beam, the D4σ width equals the conventional $1/e^2$ intensity diameter exactly. For non-Gaussian beams, the D4σ width remains well-defined and produces an $M^2$ value consistent with the divergence-based definition.

Alternative width definitions — full-width-half-maximum (FWHM), $1/e$, $1/e^2$ — produce different numerical values for non-Gaussian beams and must not be used for $M^2$ extraction without explicit conversion. ISO 11146 is unambiguous: use D4σ.

Knife-edge method principle

A knife-edge is translated transversely across the beam at fixed propagation distance $z$. As the edge moves from fully blocking to fully unblocking the beam, the transmitted power $P(x)$ traces a sigmoid curve. The intensity profile is the derivative of this curve:

$$I(x) \;=\; -\frac{dP}{dx}.$$

For a Gaussian beam, the transmitted power as a function of edge position has the form of an erf function:

$$P(x) \;=\; \frac{P_\text{max}}{2} \left[ 1 - \text{erf}\left(\frac{\sqrt{2}(x - x_c)}{w}\right) \right],$$

where $w$ is the $1/e^2$ intensity radius. Fitting the measured $P(x)$ to this form extracts $w$ directly.

For non-Gaussian beams, the second-moment width can be extracted from the knife-edge data either by numerical differentiation followed by second-moment integration, or by fitting to a generalized lineshape. For most beams of practical interest, the Gaussian erf fit gives a beam width within $\sim 5\%$ of the true D4σ.

The knife-edge method has two important practical advantages over CCD-based methods:

1. No saturation issues for high-power beams; the integrating detector accepts the full beam power minus the masked portion.
2. No wavelength-limited sensor: knife-edge measurements work at wavelengths where CCDs are not available (deep UV, mid-IR, THz).

The corresponding limitation is that the knife-edge method measures one-dimensional integrated profiles. The full two-dimensional intensity distribution is reconstructed from two orthogonal scans (knife edges along $x$ and $y$); for highly non-symmetric beams or beams with complex spatial structure, CCD methods are preferred.

Equipment

The focusing lens transforms the input beam to produce a measurable waist. For a high-quality lens (f-number $> 4$, low aberrations), the lens preserves $M^2$ and the measured waist parameters apply to the input beam after correction by the lens focal length.

For input beams that are already focused (such as fiber outputs), the lens may be unnecessary; the input beam waist is measured directly in the existing focus region.

Procedure

ISO 11146-1 specifies the measurement procedure in detail. The essential elements:

1. Configure the focused beam

Place the focusing lens in the beam path. Identify the approximate waist position $z_0$ by observing the beam diameter through the focus region. The waist position should be accessible to the knife-edge translation stage with at least $\pm 2 z_R$ of travel on either side.

For an input beam of half-divergence $\theta_\text{in}$ and lens focal length $f$, the focused waist has radius approximately

$$w_0 \;\approx\; \frac{\lambda f}{\pi w_\text{in}},$$

where $w_\text{in}$ is the input beam radius at the lens. Choose $f$ so that $w_0$ is in the range $20{-}200~\mu$m for convenient knife-edge measurement.

2. Acquire knife-edge scans at multiple z positions

ISO 11146-1 requires at least 10 measurement positions distributed as follows:

• Half within $\pm z_R$ of the waist (the near-field region)
• Half beyond $\pm 2 z_R$ from the waist (the far-field region)

Even distribution is not required; the asymmetric distribution (more points near the waist where small width variations are visible, more points in the far field where divergence is established) gives the best fit.

At each $z$ position, perform a knife-edge scan:

1. Translate the knife from fully blocking to fully unblocking the beam
2. Record transmitted power $P(x)$ at $\geq 30$ knife positions
3. Subtract any background (with knife fully blocking)

Record both the absolute $z$ position and the resulting transmitted power curve.

3. Repeat for orthogonal axis

Rotate the knife edge by 90° (or use a second knife edge oriented perpendicular to the first). Repeat the scan series for the $y$ axis. The orthogonal axis must be measured at the same $z$ positions to enable proper $M_x^2$ and $M_y^2$ extraction for asymmetric beams.

For circularly symmetric beams, only one orientation is necessary; for any beam from a non-symmetric source (edge-emitter, fiber-coupled with poor mode matching, etc.), both axes are required.

4. Extract beam width at each z

For each knife-edge scan, fit the transmitted power to the erf function:

$$P(x) \;=\; \frac{P_\text{max}}{2} \left[ 1 - \text{erf}\left(\frac{\sqrt{2}(x - x_c)}{w(z)}\right) \right] + P_\text{bg}.$$

The fit has four free parameters: $w(z)$, $x_c$, $P_\text{max}$, and $P_\text{bg}$. The fitted $w(z)$ is the $1/e^2$ intensity radius at this $z$ position.

For non-Gaussian beams, the erf fit may not capture the lineshape well; in this case, numerically differentiate $P(x)$ to obtain $I(x)$ and compute the second moment directly:

$$w(z)^2 \;=\; 4 \cdot \frac{\sum_i (x_i - x_c)^2 I(x_i)}{\sum_i I(x_i)}.$$

5. Fit the hyperbolic propagation

Plot $w(z)^2$ versus $z$ for each axis. The data follows:

$$w(z)^2 \;=\; w_0^2 + \left(\frac{M^2 \lambda}{\pi w_0}\right)^2 (z - z_0)^2.$$

Equivalently, in terms of the asymptotic divergence half-angle $\theta = M^2 \lambda / (\pi w_0)$:

$$w(z)^2 \;=\; w_0^2 + \theta^2 (z - z_0)^2.$$

Perform a least-squares fit on $w(z)^2$ versus $z$, with three free parameters: $w_0$, $z_0$, and $\theta$ (or equivalently $M^2$).

6. Compute M²

From the fit:

$$M^2 \;=\; \frac{\pi w_0 \theta}{\lambda}.$$

Report:

• $M^2$ for both axes ($M_x^2$ and $M_y^2$)
• Waist radius $w_0$ for both axes
• Waist position $z_0$
• Rayleigh range $z_R$
• 95% confidence interval on $M^2$ from the fit
• Number of data points $N$ and range $z_\text{min}$ to $z_\text{max}$ in units of $z_R$

Worked example

A 1064 nm fiber-coupled solid-state laser is focused through an aspheric lens with $f = 50$ mm. The focused beam is scanned over $\pm 4$ mm around the apparent focus position, with 10 knife-edge scans per axis at:

$z = -4.0, -3.0, -2.0, -1.0, -0.4, 0.0, 0.4, 1.0, 2.0, 4.0$ mm

(relative to nominal focus). At each $z$, the erf fit to the knife-edge scan gives:

Fitting $w^2$ versus $z$ to the hyperbolic form gives, for the x-axis:

$$w_{0,x} \;=\; 56.2~\mu\text{m}, \qquad z_{0,x} \;=\; 0.012~\text{mm}, \qquad \theta_x \;=\; 0.071~\text{rad}.$$

The implied $M^2_x = \pi w_0 \theta / \lambda = \pi \cdot 56.2 \times 10^{-6} \cdot 0.071 / (1064 \times 10^{-9}) = 11.8$.

For the y-axis: $w_{0,y} = 54.1~\mu$m, $\theta_y = 0.069$ rad, $M^2_y = 11.0$.

The fitted $z_0$ differs slightly from nominal because the operator's reference position was not exactly at the true waist; the fit recovers the true position.

The resulting $M^2 \approx 11$ is consistent with a multi-mode fiber-coupled solid-state laser. The slight asymmetry ($M^2_x > M^2_y$) reflects the residual asymmetry of the fiber-coupled output mode.

Sources of measurement error

Insufficient points outside the Rayleigh range. The hyperbolic fit's divergence parameter is constrained primarily by far-field data points ($|z - z_0| > z_R$). Measurements that cluster near the waist produce poorly-constrained divergence and large $M^2$ uncertainty. ISO 11146 specifies the asymmetric distribution (half near, half far) for this reason.

Lens aberrations. Spherical aberration in the focusing lens redistributes intensity in the focal region and can artificially inflate $M^2$. Use of a low-aberration aspheric or apochromatic lens, and avoidance of operation at large fractions of the lens NA, mitigate this contribution. Standard guideline: lens $f$-number $> 4$.

Beam clipping. If the beam diameter exceeds the lens aperture at any point in the optical path, the clipped beam has artificially elevated divergence and $M^2$. The clear aperture must accommodate $\geq 3 w$ for negligible clipping.

Background light. Stray light reaching the detector during the knife-edge scan adds a baseline that shifts the apparent beam centroid and broadens the apparent profile. Subtract background measured with knife fully blocking; for very weak beams, use lock-in detection.

Knife-edge diffraction. At small beam sizes (waist comparable to wavelength), diffraction off the knife edge distorts the measured profile. The lower beam-size limit for accurate knife-edge measurements is approximately $w_0 \geq 50 \lambda$, or $\sim 50~\mu$m at 1064 nm.

Detector saturation. During the unblocked portion of the scan, the detector receives the full beam power. Operation in the linear range of the detector must be verified — particularly for pulsed lasers, where peak power may saturate even at low average power.

Beam position drift during the measurement. A 10-position $M^2$ measurement takes minutes; if the beam drifts laterally during this time, the measured profile centroids shift and the apparent width may be inflated. Active beam stabilization or short measurement times mitigate this.

Mistaken use of FWHM or $1/e^2$ width. Reporting $M^2$ using a non-D4σ width is the most common literature error. The numerical $M^2$ value from a Gaussian-fit $1/e^2$ width and from a true D4σ measurement differ for non-Gaussian beams. For diffraction-limited or near-Gaussian beams, the difference is small (a few percent); for high-$M^2$ beams with complex structure, the difference can be 20% or more.

Polarization-dependent reflectivity at high incidence angles. For beams measured at non-normal incidence on lenses or beam splitters, polarization-dependent transmission may distort the measured profile. Use S- or P-polarized input where possible, or characterize the polarization dependence.

Validation

The fit residuals on $w^2$ versus $z$ should be randomly distributed about zero with no systematic trend. Visible curvature in the residuals indicates either lens aberration, beam clipping, or fit failure.

For a diffraction-limited reference source (single-mode HeNe or single-mode fiber output), the measured $M^2$ should fall in the range 1.0–1.1. Values significantly higher indicate measurement error, not source quality.

Repeated measurement on the same source under nominally identical conditions should agree within $\pm 5\%$. Drift larger than this indicates source instability, beam drift, or measurement chain issues.

For asymmetric beams, the geometric mean $\sqrt{M_x^2 M_y^2}$ is often quoted as a single quality metric; report both axes to enable proper system design.

References

For the controlling international standard, see ISO 11146-1:2005, Lasers and laser-related equipment — Test methods for laser beam widths, divergence angles, and beam propagation ratios. For the original theoretical treatment of beam quality and the M² formalism, see Siegman (1990, 1993) on new developments in laser resonators and beam quality. For practical comparison of knife-edge, slit, and CCD methods, see Roundy (1999) and the SPIE proceedings on laser beam profiling.