<div dir="ltr">Hello everyone, <div>I am doing NVT MD calculations using pw.x code. During the simulation, my pressure was constantly oscillating at a particular value of around 50 GPa. Many literature review papers<i> mention this line: "<b>The supercell </b></i><b><i>volume was relaxed in such a way to minimize the system's energy." </i><i>In principle, during the simulations, one could adjust the </i><i>supercell volume to make the external pressure of the inherent </i></b><i><b>structures equal to zero.</b>" I</i> have tried a lot of different things to maintain my pressure at zero, but as soon as I manage to achieve oscillatory zero pressure, some other parameters seem to shift. For example, I ran "md" at 3000K using the rescaling method, then restarted and ran "vc-relax" with fixed ion positions to achieve lower pressure, as shown below.</div><div>        " total   stress  (Ry/bohr**3)                   (kbar)     P=      123.26<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=      116.29<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=       89.51<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=       77.12<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=       60.67<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=       55.80 "<br></div><div><br></div><div>However, in the same vc-relax output<b>,</b><b> temperature</b> moved as follows:</div><div><br></div><div>"    Ekin =     0.00000000 Ry    T =    0.0 K  Etot =   -13138.95318083</div>     Ekin =     0.30599294 Ry    T =  253.6 K  Etot =   -13138.94891345<br>     Ekin =     1.41097581 Ry    T =  714.8 K  Etot =   -13138.96700552<br>     Ekin =     1.65432345 Ry    T =  878.9 K  Etot =   -13139.19284122<br>     Ekin =     1.87157146 Ry    T = 1101.3 K  Etot =   -13139.61462781<br>     Ekin =     1.67673049 Ry    T = 1142.5 K  Etot =   -13140.01208055 " <div><br></div><div>Is this temperature gain due to electron kinetic energy or why does T begin at zero? Is there something I am missing? Additionally, I used "vc-md" to control the temperature and pressure, but after that, I was unable to restart the file</div><div>My exact methodology is:</div><div>1.     restart_mode ='from_scratch',----> MD calculation (NVT) at 1000K with simple rescaling method up to 10 picoseconds</div><div>2.    restart_mode = "restart"-----> <font face="arial, sans-serif">'vc-relax' up to pressure </font>convergence threshold (0.5d0 Kbar) !!! { problem: temperature not maintained shown above} </div><div>My main question is:</div><div>Can I restart the calculation from MD to vc-md/relax and again from vc-md/relax to md ?</div><div>And if I run vc-relax or vc-md calculations, my supercell volume fluctuates slightly, requiring me to switch from an old supercell celldm to a new one in order to maintain pressure zero..</div><div><br></div><div>regards</div><div>Jayraj Anadani</div><div>Ph.D Scholar, SPU </div><div><br></div><div><br></div></div>