  We have developed Monte Carlo wave function simulation schemes to study cold
collisions between magnesium atoms in a strong red-detuned laser field. In
order to address the strong-field problem, we extend the Monte Carlo wave
function framework to include the partial wave structure of the
three-dimensional system. The average heating rate due to radiative collisions
is calculated with two different simulation schemes which are described in
detail. We show that the results of the two methods agree and give estimates
for the radiative collision heating rate for $^{24}$Mg atoms in a
magneto-optical trap based on the $^1$S$_0$-$^1$P$_1$ atomic laser cooling
transition.
