<br><br>----- Original Message -----<br>From: Vitaly Chaban &lt;vvchaban@gmail.com&gt;<br>Date: Wednesday, November 3, 2010 6:10<br>Subject: [gmx-users] pressure of the liquid in NPT and NVT<br>To: gmx-users@gromacs.org<br><br>&gt; Hi all,<br>&gt; <br>&gt; Please suggest me what is the reason for the below situation.<br>&gt; <br>&gt; The 1728 TIP4P waters system was simulated in NPT under 1bar<br>&gt; maintained with Parrinello-Rahman barostat at 300K. Some results are<br>&gt; below:<br>&gt; <br>&gt; Energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Average&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RMSD&nbsp;&nbsp;&nbsp;&nbsp; Fluct.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Drift&nbsp; Tot-Drift<br>&gt; -----------------------------------------------------------------<br>&gt; --------------<br>&gt; Total <br>&gt; Energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -58121.6&nbsp;&nbsp;&nbsp; 344.637&nbsp;&nbsp;&nbsp; 344.476&nbsp;&nbsp; 0.182549&nbsp;&nbsp;&nbsp; 36.5101<br>&gt; Temperature&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 299.913&nbsp;&nbsp;&nbsp; 4.11684&nbsp;&nbsp;&nbsp;&nbsp; 4.1138 0.00274337&nbsp;&nbsp; 0.548679<br>&gt; Pressure <br>&gt; (bar)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.39811&nbsp;&nbsp;&nbsp; 326.506&nbsp;&nbsp;&nbsp; 317.633&nbsp;&nbsp; -1.30948&nbsp;&nbsp;&nbsp;&nbsp; -261.9<br>&gt; Box-<br>&gt; X&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.74431 0.00844567 0.00843154 -8.45584e-06<br>&gt; -0.00169118<br>&gt; Density <br>&gt; (SI)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 984.768&nbsp;&nbsp;&nbsp; 6.66431&nbsp;&nbsp;&nbsp; 6.65307&nbsp; 0.0067005&nbsp;&nbsp;&nbsp; 1.34011<br>&gt; <br>&gt; ==========================================================<br>&gt; <br>&gt; The density is stable and the system looks like equilibrated one.<br>&gt; <br>&gt; The output GRO was then renamed to conf.gro and the average Box-<br>&gt; X from<br>&gt; the previous run was set as the cubic box side length. <br><br>Unless the box was equilibrated at NPT at the start of the NPT run, or the run was very long, the average value over the whole run is not all that likely to be one from the equilibrium range. Take the last frame.<br><br>&gt; The barostat<br>&gt; was switched off with no other changes in MDP. The same <br>&gt; properties for<br>&gt; NVT are below:<br>&gt; <br>&gt; Energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Average&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RMSD&nbsp;&nbsp;&nbsp;&nbsp; Fluct.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Drift&nbsp; Tot-Drift<br>&gt; -----------------------------------------------------------------<br>&gt; --------------<br>&gt; Total <br>&gt; Energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -58075.8&nbsp;&nbsp;&nbsp; 309.575&nbsp;&nbsp;&nbsp; 308.343&nbsp; -0.477939&nbsp;&nbsp; -95.5887<br>&gt; Temperature&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 299.916&nbsp;&nbsp;&nbsp; 4.08901&nbsp;&nbsp;&nbsp; 4.08678 -0.00233712&nbsp; -0.467428<br>&gt; Pressure <br>&gt; (bar)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -18.5469&nbsp;&nbsp;&nbsp; 286.376&nbsp;&nbsp;&nbsp; 285.574&nbsp;&nbsp;&nbsp; 0.37081&nbsp;&nbsp;&nbsp; 74.1628<br>&gt; <br>&gt; <br>&gt; So the pressure in NPT is ~1.4 and in NVT is ~(-19). How to <br>&gt; explain so<br>&gt; big difference in average values? Is it possible to get a value of<br>&gt; around 1 in the NVT?<br><br>We had a thread about this last week, and numerous more in the past. These observations are strongly dependent on the length of time over which you observe the average pressure, and the size of the system. Note that your fluctuations dwarf your average... See http://www.gromacs.org/Documentation/Terminology/Pressure_Coupling<br><br>Mark