Wigner function

Wigner function
Name	Wigner function
Introduced	1932
Inventor	Eugene Wigner
Field	Quantum mechanics
Applications	Quantum optics, Quantum information science

Wigner function

The Wigner function is a quasi-probability distribution function on phase space that represents quantum states in terms of position and momentum. Introduced by Eugene Wigner in 1932, it provides a bridge between classical mechanics and quantum mechanics by encoding quantum coherence and interference in a phase-space picture widely used in quantum optics and quantum information science.

Definition and physical interpretation

The Wigner function is defined for a quantum state (pure or mixed) as a real-valued function on phase space that yields the correct expectation values of symmetrically ordered operators. It gives a phase-space representation analogous to the classical Liouville theorem distribution, but unlike a true probability density it can take negative values, reflecting nonclassical effects such as quantum interference and entanglement. The Wigner function is central to operational descriptions of measurements in systems studied by groups such as those at Bell Labs historically and modern laboratories like Max Planck Institute for Quantum Optics and MIT Lincoln Laboratory.

Mathematical formulation

For a one-dimensional pure state with wavefunction ψ(x), the Wigner function W(x,p) is defined by the Fourier transform W(x,p)=\frac{1}{\pi\hbar}\int_{-\infty}^{\infty}\psi^*(x+y)\psi(x-y)e^{2ipy/\hbar}\,dy. For a mixed state described by density operator ρ, it generalizes to W(x,p)=\frac{1}{\pi\hbar}\int\langle x+y|\rho|x-y\rangle e^{2ipy/\hbar}\,dy. This construction is equivalent to the Weyl transform of ρ and is closely connected to the Weyl correspondence and the Moyal product formalism introduced by José Enrique Moyal. The inverse mapping recovers ρ from W via the Weyl quantization procedure, establishing a one-to-one correspondence between operators on Hilbert space and functions on phase space.

Properties and marginal distributions

The Wigner function has several key properties: it is real and normalized (∬W(x,p)dx dp = 1), and its marginals yield the correct position and momentum probability densities, \int W(x,p)\,dp = |\psi(x)|^2, \int W(x,p)\,dx = |\tilde\psi(p)|^2. It reproduces expectation values of symmetrically ordered operators and transforms covariantly under canonical transformations represented by the symplectic group (e.g., phase-space rotations via the Fourier transform). For Gaussian states (e.g., coherent states of the quantum harmonic oscillator), W(x,p) is positive and identical to a classical Gaussian distribution; for non-Gaussian states it can develop oscillatory regions. The purity Tr(ρ^2) can be expressed as (2πħ)∬W^2(x,p)dx dp, linking Wigner negativity to mixedness.

Quantum-classical correspondence and negativity

The Wigner function is a principal tool for studying the quantum–classical transition and semiclassical approximations such as the WKB approximation and Van Vleck propagator. In the classical limit ħ→0, Wigner functions for suitable states converge to classical phase-space distributions governed by the Liouville equation, while quantum evolution is given by the Moyal equation (an ħ-deformation of the classical Poisson bracket). Negative regions of Wigner functions are indicators of nonclassicality; their presence has been connected to quantum computational speed-up and contextuality in results by researchers studying the resource theory of nonclassicality and works referencing Scott Aaronson and Gottesman–Knill theorem contexts.

Applications in quantum optics and quantum information

In quantum optics, the Wigner function is routinely used to analyze coherent states, squeezed states, and single-photon states measured by techniques such as homodyne detection and quantum state tomography. Experimental reconstruction of Wigner functions has been achieved in setups at institutions like Caltech, Harvard University, and Ecole Normale Supérieure using optical homodyne tomography and in cavity QED and trapped-ion experiments. In quantum information, Wigner representations classify continuous-variable protocols (teleportation, error-correcting codes) and diagnose entanglement and nonclassical resources for universal quantum computing, including work on bosonic codes (e.g., GKP code by Daniel Gottesman, Alexei Kitaev and others). The function also appears in quantum chemistry and condensed-matter physics for semiclassical transport and Berry-phase related phenomena.

Computational methods and phase-space techniques

Numerical evaluation and propagation of Wigner functions employ sampling, pseudo-spectral methods, and grid-based solvers of the Moyal equation. Techniques developed in computational physics include truncated Wigner approximations, positive-P and Husimi Q representations, and stochastic approaches used in quantum optics and ultracold-atom simulations at laboratories like JILA and Copenhagen University groups. Connection to signal processing yields time–frequency analysis tools (e.g., the Wigner–Ville distribution) and fast algorithms leveraging the Fast Fourier Transform. For high-dimensional systems, tensor network methods and Monte Carlo sampling are active research directions to overcome the exponential scaling of phase-space resolution.

Category:Quantum mechanics Category:Quantum optics Category:Phase space methods