Wigner Functions and Separability for Finite Systems
  A discussion of discrete Wigner functions in phase space related to mutually unbiased bases is presented. This approach requires mathematical assumptions which limits it to systems with density matrices defined on complex Hilbert spaces of dimension p^n where p is a prime number. With this limitation it is possible to define a phase space and Wigner functions in close analogy to the continuous case. That is, we use a phase space that is a direct sum of n two-dimensional vector spaces each containing p^2 points. This is in contrast to the more usual choice of a two-dimensional phase space containing p^(2n) points. A useful aspect of this approach is that we can relate complete separability of density matrices and their Wigner functions in a natural way. We discuss this in detail for bipartite systems and present the generalization to arbitrary numbers of subsystems when p is odd. Special attention is required for two qubits (p=2) and our technique fails to establish the separability property for more than two qubits. 
Arthur O. Pittenger and Morton H. Rubin