Drift in the system even with 'fix momentum' for an Fe-Cr system with one defect

You describe your simulation to humans instead of just sharing the input deck. Are you expecting anyone to run your simulation in their head based on how they interpret/remember LAMMPS commands? Well, good luck with the answers!

I can only make two educated guesses:

  1. If a sample drifts, it can be due to poor conservation of energy due to incorrect settings (eg large timestep) or to excess kinetic energy due to the creation of defects.
  2. Having a fix momentum in the production run is hardly a good idea. If the sample still drifts, increase the equilibration run by 1M steps.
  3. Since it is the first post, Please Read This First: Guidelines and Suggestions for posting LAMMPS questions. It will avoid so much pain.

Well, I made three in the end.

1 Like