An efficient combination of the binary collision approximation with classical dynamics for stable low energy radiation damage calculations

A new technique is developed in order to estimate stable low energy radiation damage. It is based on a combination of the binary collision approximation, used to simulate the collisional phase, and quenching classical dynamics with periodic boundary conditions, used to simulate the crystal relaxation. The technique makes the simulation fast enough to collect statistics over sets as well as the energy dependence of vacancy production. The use of energy thresholds and recombination volumes in the binary collision approximation is discussed in light of the results. (orig.)