| Wigner quasi-probability distribution | |
|---|---|
| Name | Wigner quasi-probability distribution |
| Field | Quantum mechanics |
| Introduced | 1932 |
| Introduced by | Eugene Wigner |
| Related | Weyl transform, Moyal bracket, phase space |
Wigner quasi-probability distribution
The Wigner quasi-probability distribution (commonly called the Wigner function) is a representation of a quantum state in phase space that assigns real values analogous to a probability density but allows negative regions. It provides a bridge between classical mechanics and quantum mechanics by enabling phase-space calculations and semiclassical analysis. The Wigner function is widely used in theoretical and experimental studies of quantum optics, quantum information, and quantum statistical mechanics.
The Wigner function was introduced by Eugene Wigner in 1932 in a study of quantum corrections to classical statistical mechanics. It grew out of earlier ideas by Hermann Weyl on a phase-space correspondence (the Weyl transform) and was later connected to work by José Enrique Moyal who developed the statistical formulation known as Moyal bracket mechanics. Historically, the Wigner distribution provided a way to discuss quantum interference and decoherence using phase-space intuition familiar from Liouville's theorem in classical dynamics. It has since become a standard tool in analyses performed at institutions such as Harvard University, University of Cambridge, and laboratories including Bell Labs and MIT Lincoln Laboratory where quantum optics and quantum information experiments are developed.
For a single continuous degree of freedom with coordinate x and momentum p, the Wigner function W(x,p) for a density operator ρ is defined by the Weyl transform: W(x,p) = (1/πħ) ∫_{-∞}^{∞} ⟨x + y | ρ | x - y⟩ e^{-2ipy/ħ}\, dy. This formula links the density matrix in the position basis with a real-valued phase-space function. Equivalently, for a pure state with wavefunction ψ(x), one writes W(x,p) = (1/πħ) ∫ ψ*(x + y) ψ(x - y) e^{-2ipy/ħ}\, dy. The Wigner function obeys marginal properties: integrating over p yields the position probability density |ψ(x)|^2, and integrating over x gives the momentum-space density. The Wigner formalism generalizes to multimode systems, spin coherent states, and discrete-variable systems via finite-dimensional Wigner functions used in quantum computing and quantum error correction studies.
The Wigner function is real and normalized but not a true probability distribution because it may take negative values, a signature of nonclassicality. Its moments reproduce expectation values of symmetrically ordered operators (Weyl ordering). Time evolution under a Hamiltonian H is given by the Moyal equation, which reduces to the classical Liouville equation in the ħ → 0 limit, with quantum corrections expressed via higher-order derivatives (the Moyal series). Negativity of the Wigner function is often used as a quantitative indicator of quantum interference and is closely linked to resources in quantum computation and quantum metrology where negative quasi-probabilities can enable tasks impossible classically.
Common examples include the Wigner functions of the harmonic oscillator eigenstates, coherent states, squeezed states, and Fock states in quantum optics. Coherent states yield Gaussian positive Wigner functions, while Fock states exhibit oscillatory patterns with negative regions. Applications span analysis of decoherence in open quantum systems, phase-space visualization of entanglement in bipartite systems, and evaluation of semiclassical approximations in chemical physics (e.g., molecular dynamics). Experimental reconstruction of Wigner functions is routinely performed in cavity quantum electrodynamics experiments at facilities such as IBM Research and Caltech laboratories, where state tomography confirms nonclassical features for quantum technologies.
The Wigner function is one member of a family of phase-space quasiprobability distributions, including the Glauber–Sudarshan P representation and the Husimi Q function used in quantum optics and statistical descriptions. The P function can be highly singular while the Q function is a smoothed, positive-definite convolution of the Wigner function. These representations are related by Gaussian convolutions or operator orderings (normal, anti-normal, and symmetric). In quantum information, discrete Wigner functions are compared with stabilizer formalism and contextuality criteria (e.g., links to the Kochen–Specker theorem and Bell's theorem) for assessing classical simulability of quantum circuits.
Numerical evaluation of Wigner functions uses fast Fourier transforms for the defining integrals and grid-based schemes for multidimensional phase spaces. Semiclassical propagation methods, such as the van Vleck–Gutzwiller propagator and Wigner–Moyal trajectory approaches, approximate dynamics for large systems. Experimentally, quantum state tomography reconstructs the Wigner function from homodyne measurements with inverse Radon transform techniques (quantum homodyne tomography), or via direct parity measurements in cavity QED. Practical implementations in superconducting qubits and trapped ions use tailored measurement protocols to infer Wigner negativities as benchmarks for quantum control.
Generalizations include multi-component and spin Wigner functions, discrete Wigner distributions for finite Hilbert spaces (useful in quantum computing and quantum error correction), and relativistic Wigner functions in quantum field theory for phase-space transport. Other extensions incorporate open-system effects via Lindblad operators to study dissipative evolution in phase space, and noncommutative geometry adaptations for systems with noncanonical commutation relations. Theoretical developments continue to connect Wigner methods with resource theories of nonclassicality, contextuality, and advances in quantum simulation platforms.