I am working with a BCC single crystal solid. Considering a small simulation box of 35 X 35X 100 unit cells at room temperature and zero pressure. After that I want to preheat my sample to 2100K. For this,
Firstly, I raised the temperature from 300K to 2100K keeping the pressure at 0 using NPT and I achieved it.
And then, ran a NVE ensemble for 60000 steps to equilibrate the system using the given script.
velocity bulk create 300 12345 rot yes dist gaussian
fix equilibration bulk npt temp 300 2100 0.1 iso 0 0 1 drag 10
fix 1 bulk nve
I wanted to know, even after giving a long run of NVE. is it important to minimise the system, by giving a box/relax statement. Or, is NVE taking care of the equilibration of the crystal.