Classical field theory with fermions

Classical field theory simulations have been essential for our understanding of non-equilibrium phenomena in particle physics. In this talk we discuss the possible extension of the bosonic classical field theory simulations to include fermions. In principle we use the inhomogeneous mean field approximation as introduced by Aarts and Smit. But in practice we turn from their deterministic technique to a stochastic approach. We represent the fermion field as an ensemble of pairs of spinor fields, dubbed male and female. These c-number fields solve the classical Dirac equation. Our improved algorithm enables the extension of the originally 1+1 dimensional analyses and is suitable for large-scale inhomogeneous settings, like defect networks.