A fundamentally more straightforward way to model multiple atomic collisions is provided by molecular dynamics (MD) simulations, in which the time evolution of a system of atoms is calculated by solving the equations of motion numerically. Special MD methods have been devised in which the number of interactions and atoms involved in MD simulations have been reduced in order to make them efficient enough for calculating ion ranges. The MD simulations this automatically describe the nuclear stopping power. The electronic stopping power can be readily included in molecular dynamics simulations, either as a frictional force or in a more advanced manner by also following the heating of the electronic systems and coupling the electronic and atomic degrees of freedom.