In this chapter, we will develop methods for mapping operator equations to equivalent c-number equations. This results in a continuous phase-space representation of a many-body quantum system, using phase-space distributions instead of density matrices. In order to do this, we will first find c-number representations of operators. We have already seen one such representation, the Glauber–Sudarshan P-representation for the density matrix. We now introduce several more such representations. The main focus will be on the truncated Wigner representation, valid at large photon number, and the positive P-representation, which uses a double-dimensional phase space and exists as a positive probability for all quantum density matrices.

Phase-space techniques have a great advantage over conventional matrix-type solutions to the Schrödinger equation, in that they do not have an exponential growth in complexity with mode and particle number. Instead, the equations that describe the dynamics of these c-number representations of the density matrix are Fokker–Planck equations. These have equivalent stochastic differential equations, which behave as c-number analogs of Heisenberg-picture equations of motion.

This means that problems that would be essentially impossible to solve using conventional number-state representations can be transformed into readily soluble differential equations. In many cases, no additional approximations, such as perturbation theory or factorization assumptions, are needed.