<div dir="ltr">Hello, QE community,<div>I am simulating 128 metallic atoms at 3000K under the NVT ensemble approach with a velocity-rescaling thermostat in PWscf MD (BOMD-AIMD). After 2.2 picoseconds, I checked some parameters. The temperature, internal energy, and pressure all equilibrate with time (ps), but when I plot internal energy (E=F+TS) vs. temperature (K), the energy value circulates around 3000K, which I set.</div><div>Is it normal, or am I wrong?</div><div>Does this non-zero pressure affect the total energy?</div><div>Should I perform vc-relax calculation by fixing the cell volume to decrease the pressure?</div><div><ul><li>find the attached graph at plot.png</li></ul></div><img src="cid:ii_lo8o5ec90" alt="plot.png" width="473" height="355"><br><div><br></div><div>Thanks</div><div><br></div><div>regards</div><div>Jayraj Anadani</div><div>SPU Ph.D. Scholar</div></div>