probability - an example of a quantum system for which wigner function transitions to negative values

I want to check my understanding of the Wigner transform and try to understand why and how exactly the probabilistic interpretation drops down as the function goes to zero and then to negative values

So, suppose we have a free quantum oscillator with integer eigenvalues of energy and (boson) occupation.

Questions:

• 
  

  Is it possible to start the system in a state that will have a well defined Wigner phase space probability density (with an analog interpretation as a Liouville phase space density), and it will freely evolve into a state for which the Wigner density becomes negative-valued?. In other words, is the positive-definite property of a Wigner density invariant under free evolution?

  
  


• 
  

  is there a simple system (hopefully a simple harmonic oscillator based state) for which the Wigner density transitions from being positive-definite to being negative at certain regions?

  
  


• 
  

  what happens with the probability interpretation in the boundary of the transition between positive-definite and negative at certain regions in general?

  
  

Answer

The Wigner function is the fourier transform in one variable of the density matrix for the single particle $\rho(x,y)$. If you Fourier transform in y, you get the Wigner function $\rho(x,p)$. This is important to understand, because it explains why the Wigner function is at all interesting, and why it obeys simple dynamics. It also shows that it doesn't have a probability interpretation away from the classical limit, because only the on-diagonal density matrix elements are probabilities.

Question 1: is is possible for Wigner's density to be negative after starting positive?

The answer is yes for a general potential, but for the special case of a Harmonic oscillator, the answer is no, because time evolution just rotates phase space. For the other special case of a free particle, a Gaussian wavefunction just spreads into a wider Gaussian, so it is not possible there either. For this reason, you have a hard time getting an example.

It is also true that semiclassically, the motion is along the classical trajectory, with shear according to the change in period with increasing J. So if a semiclassical Wigner matrix is positive near a single non-chaotic trajectory, it won't become negative, near the trajectory, at least not for a long time.

But it is very easy to get any kind of Wigner densities away from the semiclassical regime.

Question 2: What about the "probability density on phase space" interpretation?

This interpretation is faulty. The density matrix is only a probability for on-diagonal elements. So the Fourier transform of the Wigner density in either x or p has a diagonal probability density interpretation at all times, but the x,p position Wigner function is nothing. It's just a complicated encoding of off-diagonal matrix elements that is not particularly special or useful, away from the semi-classical limit.

To see that it is no good, consider a plane-wave state, whose phase space Wigner density is the plane-wave profile times a delta function at a certain value of p. If you scatter any lump into approximate plane waves, the end state of scattering is a bunch of complex stuff that has no interpretation as a phase-space probability.

The proper interpretation of Wigner's phase space density is that it is a Fourier transform of the density matrix, nothing more.