MD simulations start by defining an initial set of positions and velocities for the particles in a system. Using Newton's laws of motion, the future positions and velocities of these particles are calculated iteratively over small time steps. The forces between particles are typically derived from potential energy functions, like the Lennard-Jones potential or electrostatic interactions.