The crystal lattice of cubic scandium fluorine (ScF₃) exhibits negative thermal expansion (NTE) over a wide range of temperatures from 10 K to 1100 K. Here the NTE effect in ScF₃ is studied using atomistic simulations based on empirical and ab initio molecular dynamics (AIMD) in the isobaric-isothermal (NpT) ensemble. The temperature dependence of the average lattice constant, the Sc–F–Sc bond angle distribution and the radial distribution functions were obtained. Crossover from the NTE to positive thermal expansion occurring at about 1100 K is reproduced by AIMD simulations in agreement with the known experiment data. At the same time, empirical MD model fails to reproduce the NTE behaviour and suggests an expansion of the ScF₃ lattice with increasing temperature. However, both MD models predict strong anisotropy of fluorine atom thermal vibration amplitude, being larger in the direction orthogonal to the Sc–F–Sc atom chain.