Cascade energies, and the rate of their occurrence, were used to model the damage that was due to a flux of 14MeV neutrons. The cascades were generated by molecular dynamics simulations, for recoil energies of 5 to 25keV. The numbers of each type of defect, the cluster size-distributions and the positions of the defects within the crystal, were monitored as a function of time. Simulated defect cluster densities, as a function of doses of up to 0.1dpa at room temperature, were compared with experimental data. It was found that the stimulated cluster densities were within a factor of about 2, of the experimental results, over several orders of magnitude of the dose.

Simulation of the Kinetics of Defect Accumulation in Copper under Neutron Irradiation. H.L.Heinisch, B.N.Singh: Journal of Nuclear Materials, 1999, 271-272, 46-51