Fermionic Molecular Dynamics models a system of fermions by means of a trial many-body state composed of an antisymmetrized product of single-particle states which are localized gaussians in coordinate and momentum space. The parameters specifying them are the analogue to the variables in classical molecular dynamics. The time-dependent variational principle yields the equations of motion which are solved for heavy-ion collisions. These calculations show a great variety of phenomena including explosion, sequential fragmentation and multifragmentation. In this contribution we present an approximation for the expectation value of the two-body interaction, define the ensemble of events and show that large fluctuations arise therein.