<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   Averaged Physical Quantities<br>                      accomulated      this run<br>   ekinc         :        0.00001       0.00001 (AU)<br>   ekin          :       94.60238      94.60238 (AU)<br>  
 epot          :     -343.05127    -343.05127 (AU)<br>   total energy  :     -211.98234    -211.98234 (AU)<br>   temperature   :        0.75835       0.75835 (K )<br>   enthalpy      :     -211.93045    -211.93045 (AU)<br>   econs         :     -211.93020    -211.93020 (AU)<br>   pressure      :       -2.19645      -2.19645 (Gpa)<br>   volume        :     1526.61474    1526.61474 (AU)<br><br>Why the averaged temperature and
 pressure in the output file are far  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>