<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Dear all,<br><br>I would like to study MD with NPT ensemble. I am learning example 18 in QE 4.2.1 which uses cp.x to perform molecular dynamics simulation of SiO2.<br><br>For the vc-cp simulation, input file specify tempw=300 and press=1.0. However, the &nbsp; Averaged Physical Quantities<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; accomulated&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; this run<br>&nbsp;&nbsp; ekinc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00001&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00001 (AU)<br>&nbsp;&nbsp; ekin&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 94.60238&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 94.60238 (AU)<br>&nbsp;&nbsp;
 epot&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; -343.05127&nbsp;&nbsp;&nbsp; -343.05127 (AU)<br>&nbsp;&nbsp; total energy&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; -211.98234&nbsp;&nbsp;&nbsp; -211.98234 (AU)<br>&nbsp;&nbsp; temperature&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.75835&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.75835 (K )<br>&nbsp;&nbsp; enthalpy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; -211.93045&nbsp;&nbsp;&nbsp; -211.93045 (AU)<br>&nbsp;&nbsp; econs&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; -211.93020&nbsp;&nbsp;&nbsp; -211.93020 (AU)<br>&nbsp;&nbsp; pressure&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2.19645&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2.19645 (Gpa)<br>&nbsp;&nbsp; volume&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; 1526.61474&nbsp;&nbsp;&nbsp; 1526.61474 (AU)<br><br>Why the averaged temperature and
 pressure in the output file are far&nbsp; away from the target value set in the input? <br>One more question, what is the purpose of using cp, cp-restart, vc-cp, and vc-cp-restart together? Is this a routine to study MD, i.e. combing these four? In other words, could we just get the vc-cp part for investigating NPT ensemble?<br><br>Best regards<br>Jingyun<br><br><br><br></td></tr></table><br>