We review phase-space simulation techniques for fermions, showing how a Gaus- sian operator basis leads to exact calculations of the evolution of a many-body quantum system in both real and imaginary time. We give particular application to the Hubbard model and to the problem of molecular dissociation of bosonic molecules into pairs of fermionic atoms.