TOPICS
Search

Wigner Distribution


The Wigner distribution of a square-integrable complex function f is the time-frequency function

 W_f(x,xi)=int_(-infty)^inftyf(x+t/2)f(x-t/2)^_e^(-2piixit)dt.
(1)

Its marginal integrals, interpreted weakly when necessary, recover the position and frequency energy densities,

int_(-infty)^inftyW_f(x,xi)dxi=|f(x)|^2
(2)
int_(-infty)^inftyW_f(x,xi)dx=|f^^(xi)|^2,
(3)

where f^^ is the Fourier transform with kernel e^(-2piixix). If ||f||_2=1, both marginals are probability density functions. Under this normalization, the real-valued Wigner distribution is a quasiprobability distribution, since it need not be nonnegative and therefore need not be a joint probability density function.

For a normalized wave function psi, the momentum convention used in quantum mechanics is

 w_psi(x,p)=1/(2pih)W_psi(x,p/(2pih))=1/(2pih)int_(-infty)^inftypsi(x+y/2)psi(x-y/2)^_e^(-ipy/h)dy,
(4)

where p is momentum and h>0 is the reduced Planck constant (Case 2008, p. 938). The factor 1/(2pih) ensures that the integral of w_psi over phase space is 1.

For the quantum harmonic oscillator with Hamiltonian

 H(x,p)=(p^2)/(2m)+(momega^2x^2)/2,
(5)

where m>0 and omega>0, the normalized stationary wave functions are scaled Hermite functions. Their Wigner functions are

 w_n(x,p)=((-1)^n)/(pih)e^(-r^2)L_n(2r^2),
(6)

where n=0, 1, 2, ..., L_n is a Laguerre polynomial, l=sqrt(h/(momega)), and r^2=x^2/l^2+l^2p^2/h^2 (Curtright et al. 1998). The ground state n=0 gives a positive Gaussian function, whereas the first excited state has w_1(0,0)=-1/(pih). A coherent state, obtained by displacing the ground state to position x_0 and momentum p_0, has the Gaussian function (Case 2008, pp. 940-942)

 w_(x_0,p_0)(x,p)=1/(pih)exp[-((x-x_0)^2)/(l^2)-(l^2(p-p_0)^2)/(h^2)].
(7)

For a time-independent Hamiltonian H=p^2/(2m)+V(x) with V a polynomial of degree at most 2, the exact quantum evolution of the Wigner function w_t obeys the classical transport equation

 (partialw_t)/(partialt)=(partialH)/(partialx)(partialw_t)/(partialp)-(partialH)/(partialp)(partialw_t)/(partialx)={H,w_t},
(8)

where {H,w_t} is the Poisson bracket. The quantum corrections vanish because they contain derivatives of V of order at least 3 (Case 2008, pp. 941-942). If phi_t is the flow determined by Hamilton's equations, then

 w_t(x,p)=w_0(phi_(-t)(x,p)).
(9)

For the harmonic oscillator this is a rigid clockwise rotation in the scaled coordinates (x/l,lp/h) with angular speed omega, for every initial state. In particular, a coherent state's Gaussian function retains its shape while its center follows the classical orbit.


See also

Fourier Transform, Fractional Fourier Transform, Probability Density Function

Explore with Wolfram|Alpha

WolframAlpha

More things to try:

References

Case, W. B. "Wigner Functions and Weyl Transforms for Pedestrians." Am. J. Phys. 76, 937-946, 2008. https://doi.org/10.1119/1.2957889.Cohen, L. "Generalized Phase-Space Distribution Functions." J. Math. Phys. 7, 781-786, 1966. https://doi.org/10.1063/1.1931206.Curtright, T.; Fairlie, D.; and Zachos, C. "Features of Time-Independent Wigner Functions." Phys. Rev. D 58, 025002, 1998. https://doi.org/10.1103/PhysRevD.58.025002.Wigner, E. "On the Quantum Correction for Thermodynamic Equilibrium." Phys. Rev. 40, 749-759, 1932. https://doi.org/10.1103/PhysRev.40.749.

Cite this as:

Weisstein, Eric W. "Wigner Distribution." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/WignerDistribution.html

Subject classifications