ABSTRACTWe present results of molecular dynamics computer simulation studies of the threshold energy for point defect production in silicon. We employ computational cells with 8000 atoms at ambient temperature of 10 K that interact via the Stillinger-Weber potential. Our simulations address the orientation dependence of the defect production threshold as well as the structure and stability of the resulting vacancy-interstitial pairs. Near the <111> directions, a vacancy-tetrahedral-interstitial pair is produced for 25 eV recoils. However, at 30 eV recoil energy, the resulting interstitial is found to be the <110> split dumbbell configuration. This Frenkel pair configuration is lower in energy than the former by 1.2 eV. Moreover, upon warming of the sample from 10 K the tetrahedral interstitial converts to a <110> split before finally recombining with the vacancy. Along <100> directions, a vacancy-<110> split interstitial configuration is found at the threshold energy of 22 eV. Near <110> directions, a wide variety of closed replacement chains are found to occur for recoil energies up to 45 eV. At 45 eV, the low energy vacancy-<110> split configuration is found. At 300 K, the results are similar. We provide details on the atomic structure and relaxations near these defects as well as on their mobilities.