Abstract A realistic characterization of the interaction of ions and neutral particles with device walls is important for any edge plasma calculation. Present reflection models vary in detail and computational efficiency. This paper presents a data set for the distribution of the reflection coefficient, R N , over reflected energy, polar, and azimuthal angles, as functions of incident polar angle and energy. These results have been computed using a vectorized version of the TRIM Monte Carlo code. The data are stored using an algorithm for reducing the data into three one-dimensional distributions, resulting in a realistic reflection model which can be used very efficiently in plasma edge calculations.