Dear Justin:<br><br>Thanks for your response. Here is the complete my .mdp file:<br><br>title               =  pdm<br>cpp                 =  /lib/cpp<br>constraints         =  all_bonds<br>integrator          =  md<br>dt                  =  0.004   ; ps !<br>
nsteps              =  2500000   ; total 10ns.<br>nstcomm             =  1<br>nstxout             =  50000<br>nstvout             =  50000<br>nstfout             =  0<br>nstlog              =  5000<br>nstenergy           =  5000<br>
nstxtcout           =  25000<br>nstlist             =  10<br>ns_type             =  grid<br>pbc                 =  xyz<br>coulombtype         =  PME<br>rlist               =  1.4<br>rcoulomb            =  1.4<br>rvdw                =  1.4<br>
fourierspacing      =  0.20<br>pme_order           =  4<br>ewald_rtol          =  1e-5<br>; Berendsen temperature coupling is on in one groups<br>Tcoupl              =  berendsen<br>tc_grps             =  PDM<br>tau_t               =  0.1<br>
ref_t               =  300<br>; Energy monitoring<br>energygrps          =  PDM<br>; Isotropic pressure coupling is now on<br>Pcoupl              =  berendsen<br>pcoupltype          =  isotropic<br>;pc-grps            =  PDM<br>
tau_p               =  1.0<br>ref_p               =  1.0<br>compressibility     =  4.5e-5<br><br>; Generate velocites is off at 300 K.<br>gen_vel             =  yes<br>gen_temp            =  300.0<br>gen_seed            =  100000<br>
<br><br><br><div class="gmail_quote">On Thu, May 14, 2009 at 4:06 PM, Justin A. Lemkul <span dir="ltr">&lt;<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div><div></div><div class="h5"><br>
<br>
Yanmei Song wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear All:<br>
<br>
I have question about the pressure coupling. I have done a 10ns simulation with 19800 atoms for 120 large molecules using the following pressure coupling.<br>
<br>
Tcoupl              =  berendsen<br>
tc_grps             =  PDM<br>
tau_t               =  0.1<br>
ref_t               =  300<br>
<br>
Pcoupl              =  berendsen<br>
pcoupltype          =  isotropic<br>
;pc-grps            =  PDM<br>
tau_p               =  1.0<br>
ref_p               =  1.0<br>
compressibility     =  4.5e-5<br>
<br>
Then I did g_energy for the last 3ns and got the results:<br>
<br>
Energy                      Average       RMSD     Fluct.      Drift  Tot-Drift<br>
-------------------------------------------------------------------------------<br>
Potential                    -98061          0          0    0.61668    1850.04<br>
Temperature                 303.561    109.602    109.602 0.000181791   0.545372<br>
Pressure (bar)               4.4884    11110.8    11109.8  -0.169835   -509.506<br>
<br>
For such a long run the pressure drift is still too much and seem hasn&#39;t approached 1bar. Does it mean the system hasn&#39;t reach equilibrium yet. I did a similar system by using the same method. it just take 2 or 3ns to reach the equilibrium. and the pressure is around 1.01after the run. The only difference is the cutoff changing from 1.2 to 1.4.  Does the cufoff of 1.4 is too large to make the system running slower. Or the pressure coupling method is not working well.  Anyone can give me any suggestions?<br>

<br>
<br>
</blockquote>
<br></div></div>
I think it will depend on the interplay of other parameters as well.  Posting a complete .mdp file may be more helpful.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">
-- <br>
Yanmei Song<br>
Ph.D. Candidate<br>
Department of Chemical Engineering<br>
Arizona State University<br>
<br>
<br></div>
------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</blockquote>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
_______________________________________________<br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Yanmei Song<br>Ph.D. Candidate<br>Department of Chemical Engineering<br>Arizona State University<br>