*[Follow-up post: In what sense is the Wigner function a quasiprobability distribution?]*

I’ve never liked how people introduce the Wigner function (aka the Wigner quasi-probability distribution). Usually, they just write down a definition like

(1)

and say that it’s the “closest phase-space representation” of a quantum state. One immediately wonders: What’s with the weird factor of , and what the heck is ? Usually, the only justification given for the probability interpretation is that integrating over one of the variables recovers the probability distribution for the other (if it were measured):

(2)

where is just the density matrix in the momentum basis. But of course, that doesn’t really tell us why we should think of , as having anything to do with the (rough) value of conditional on a (rough) value of .

Well now I have a much better idea of what the Wigner function actually is and how to interpret it. We start by writing it down in sane variables (and suppress ):

(3)

So the first step in the interpretation is to consider the function

(4)

which appears in the integrand. This is just the (position-space) density matrix in rotated coordinates and . There is a strong sense in which the off-diagonal terms of the density matrix represent the quantum coherence of the state between different positions, so indexes how far this coherence extends; large values of indicate large spatial coherence. On the other hand, indexes how far down the diagonal of the density matrix we move; it’s the average position of the two points between which the off-diagonal terms of the density matrix measures coherence. (See the figure below.)

So then we see that the Wigner function is just the Fourier transform of with respect to the coherence variable .

(5)

We can express this more compactly as

(6)

where denotes the Fourier transform from to and denotes the relevant coordinate rotation from to . This interpretation has several of advantages.

First, all the numerical constants in the formula for the Wigner function can be immediately traced back to the Fourier transform. (Different Fourier transform conventions will of course have different constants.)

Second, transforming the density matrix into a Wigner function is a manifestly reversible operation, since both the Fourier transform and the rotation are reversible.^{a } So the Wigner function is just another lens with which we can look at the state of a quantum system. We could just as well always work with Wigner functions when analyzing the dynamics of a quantum system, although the Hamiltonian evolution turns out to be a little uglier to express.

Third, we can also see how the Wigner function is related to the momentum-space representation of the density matrix by “filling in” all the possible applications of Fourier transforms and coordinate rotations. (See the next figure.) This lets us immediately and transparently derive that

(7)

where is the relevant coordinate rotation of .

The (doubly) Fourier transformed Wigner function^{b }

(8)

(aka the Moyal characteristic function, aka the Wigner characteristic function, especially when the arguments are swapped) also has a nice relationship to the Husimi Q function . Here, is just a condensed notation for phase space coordinates and is the (unsqueezed) coherent state centered at . It’s well known that the Husimi Q function is just the Wigner function convolved with a Gaussian:

(9)

where and . (Up to normalization, the Wigner representation of the coherent state at is a Gaussian: .)

But of course, by the convolution theorem, a convolution in real space is just multiplication in Fourier space:

(10)

where

(11)

In fact, the Fourier transformed Wigner function has more meaning than you might think. Structure of the Wigner function inside areas smaller than are known to be associated with the type of long-scale coherence that is typically destroyed by environmental decoherence.^{c } But the Fourier transformed Wigner function inverts this; the sub-Planckian structure now is on large scales, while the supra-Planckian structures (which gives the classical phase space probability distribution) is compressed to small scales. When calculating the Husimi Q function, we see that the Gaussian just suppresses the large-scale structure (i.e. the coherence) while leaving the small-scale structure untouched. (Of course, we could probably have read this off directly from the convolution.)

*[Edited 2014-6-10]*

### Footnotes

(↵ returns to text)

- Remember that, aside from being a complex-valued function on , the density matrix must also be positive. I haven’t bothered to figure out how that requirement translate into restrictions on the Wigner function beyond the simple fact that it’s a real-valued function on . The unit-trace requirement on is equivalent to the requirement that integrating over its domain equals unity.↵
- I’m going to stop trying to bother with normalization.↵
- Furthermore, this structure is responsible for superpositions being useful in quantum enhanced measurements.↵

Oh.. man! That is one of the *BEST* explanations I have seen so passionately explained on the internet on that too a serious subject! Thanks a bunch! Keep writing. 🙂