I'm working in molecular dynamics with transition metal clusters. The aim is find the melting point, so, it's necessary scale the velocities to increment temperature.
I know that I should use instantaneous temperature derived from the equipartition theorem but my question is about how to increment the temperature scaling the velocities.
Any advice to code the increment of temperature and make a caloric curve (Energy.vs. Temperature)?