AbstractSpecific features of the angular distributions of accelerated neutral nitrogen atoms at the grazing angles of incidence on the Al(001) crystal surface have been investigated by the computer simulation method. The N–Al pair interaction potential is approximated by the three-parameter Morse potential with the energydependent coefficients. The angular distributions of scattered atoms have been simulated taking into account the interaction between atoms and several atomic layers in the lattice and the atomic displacement during thermal oscillations. The parameters of the pair potential of accelerated neutral nitrogen atoms in the energy range from 10 to 70 keV have been determined according to the best agreement between the calculated dependence of the rainbow scattering angle on the energy of particles incident on the crystal surface and the available experimental data.