Repulsive interatomic potentials calculated at three levels of theory

The high-energy repulsive interaction between nuclei at distances much smaller than the equilibrium bond length is the key quantity determining the nuclear stopping power and atom scattering in keV and MeV radiation events. This interaction is traditionally modeled within orbital-free density functional theory with frozen atomic electron densities, following the Ziegler-Biersack-Littmark (ZBL) model. In this work, we calculate atom pair specific repulsive interatomic potentials with the ZBL model, and compare them to two kinds of quantum chemical calculations - second-order Møller-Plesset perturbation theory in flexible Gaussian basis sets as well as density functional theory with numerical atomic orbital basis sets - which go well beyond the limitations in the ZBL model, allowing the density to relax in the calculations. We show that the repulsive interatomic potentials predicted by the two quantum chemical models agree within $\sim$ 1% for potential energies above 30 eV, while the ZBL pair-specific potentials and universal ZBL potentials differ much more from either of these calculations. We provide new pair-specific fits of the screening functions in terms of 3 exponentials to the calculations for all pairs $Z_1$-$Z_2$ for $1 \leq Z_i \leq 92$, and show that they agree within $\sim 2$% with the raw data. We use the new potentials to simulate ion implantation depth profiles in single crystalline Si and show very good agreement with experiment. However, we also show that under channeling conditions, the attractive part of the potential can affect the depth profiles. The full data sets of all the calculated interatomic potentials as well as analytic fits to the data are shared as open access.