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