Within QMD simulations we demonstrate the effect of virial corrections on heavy ion reactions. Unlike in standard codes, the binary collisions are treated as non-local so that the contribution of the collision flux to the reaction dynamics is covered. A comparison with standard QMD simulations shows that the virial corrections lead to a broader proton distribution bringing theoretical spectra closer towards experimental values. Complementary BUU simulations reveal that the non-locality enhances the collision rate in the early stage of the reaction. It suggests that the broader distribution appears due to an enhanced pre-equilibrium emission of particles
