The Marlowe code was used, with a minimum cut-off energy, to model a cascade that was then treated by using molecular dynamics methods. Thus, a 1keV cascade was simulated by using the Marlowe code and a cut-off energy of 200eV. This was then followed by a 100fs molecular dynamics simulation. The final configuration was compared with that obtained by using molecular dynamics alone. It was found that binary collisions tended to over-estimate the range and cascade volume. Because differences arose in the low-energy part of the trajectories, the performance of the linked simulation improved greatly at higher cut-off values.

M.Jaraiz, G.H.Gilmer, D.M.Stock, T.Diaz de la Rubia: Nuclear Instruments and Methods in Physics Research B, 1995, 102[1-4], 180-2