Wigner functions in optics: describing beams as ray bundles and pulses as particle ensembles

This tutorial gives an overview of the use of the Wigner function as a tool for modeling optical field propagation. Particular emphasis is placed on the spatial propagation of stationary fields, as well as on the propagation of pulses through dispersive media. In the first case, the Wigner function gives a representation of the field that is similar to a radiance or weight distribution for all the rays in the system, since its arguments are both position and direction. In cases in which the field is paraxial and where the system is described by a simple linear relation in the ray regime, the Wigner function is constant under propagation along rays. An equivalent property holds for optical pulse propagation in dispersive media under analogous assumptions. Several properties and applications of the Wigner function in these contexts are discussed, as is its connection with other common phase-space distributions like the ambiguity function, the spectrogram, and the Husimi, P, Q, and Kirkwood–Rihaczek functions. Also discussed are modifications to the definition of the Wigner function that allow extending the property of conservation along paths to a wider range of problems, including nonparaxial field propagation and pulse propagation within general transparent dispersive media.