Hi, Justin,
Thanks for your rapid reply. Let me decribes the question in more details.
First ,two gaseous mathane systems and one liquid water system are equilibrated
at 300 bar and 250 K in advance. Then, the gas systems are put at the both
sides of liquid. The combined system is performed for 3ns with the
Parrinello-Rahman barostat, which produces a average pressure of 327 bar. I am
strange to the value and several tests are run, mentioned previously. All these
test runs take the end gro file of the 3ns-run as the beginning. When the
berendsen barostat is used, the pressure is OK. When the Parrinello-Rahman
barostat is used, the questionary pressure occurs again. Although I will try to
make a long run for equilibrium, I think the equilibrium has reached because
the ~330 bar is always obtained no matter how long simulations are sampled,
from 0.2 ns to 3 ns. The version of Gromacs I used is 4.0.4. By the way, how
long did you equilibrate your membrane system and how many atoms in it? Thanks
again.
Guang-Jun
> Date: Sat, 28 Feb 2009 07:43:29 -0500
> From: jalem...@vt.edu
> To: gmx-users@gromacs.org
> Subject: Re: [gmx-users] Cannot obtain the expected pressure using the
> Parrinello-Rahman coupling in a gas-liquid system
>
>
>
> GuoGuangjun wrote:
> > In fact, I once run the gas-liquid interface for 3 ns, the pressure is
> > still as high as 327 bar.
>
> There are two possibilities in my mind:
>
> 1. 3 ns is still not long enough.
> 2. Your gas phase is causing the problem. If the pressure is right for water
> alone, then the problem likely comes from the inclusion of your gas phase and
> whatever parameters you're using for whatever species you have in this phase.
> Examine your parameters and consider how they were derived and whether or not
> they are adequate.
>
> Also, what version of Gromacs are you using? There were some changes made to
> the Parrinello-Rahman barostat between 3.3.x and 4.0.x.
>
> -Justin
>
> > Guang-Jun
> >
> >
> > > Date: Sat, 28 Feb 2009 07:19:30 -0500
> > > From: jalem...@vt.edu
> > > To: gmx-users@gromacs.org
> > > Subject: Re: [gmx-users] Cannot obtain the expected pressure using
> > the Parrinello-Rahman coupling in a gas-liquid system
> > >
> > >
> > >
> > > GuoGuangjun wrote:
> > > >
> > > > Hi, All,
> > > > I am simulating an interface system composed of methane gas and liquid
> > > > water in the NPT ensemble. The temperature is well controlled by using
> > > > either the Berendsen or the Nose-Hoover T-couplings. As for pressure,
> > > > the Berendsen P-coupling also works well. However, the
> > Parrinello-Rahman
> > > > P-coupling always produces a higher pressure than the expected. For
> > > > example, I set ref_p=300 bar in the mdp file but obtain an average
> > > > pressure as ~330 bar (over 0.5 ns). I test the tau_p in the range
> > 0.1 ~
> > > > 10 ps, the case does not change. When I test the sin! ! gle-phase
> > system,
> > > > i.e., pure water, the Parrinello-Rahman P-coupling can indeed produce
> > > > the expected pressure with the deviation less than 2 bar. What is the
> > > > matter in the two-phase system?
> > > >
> > >
> > > When I've used the Parrinello-Rahman barostat for membrane systems,
> > it works
> > > just fine. Perhaps 500 ps is too short of a timeframe to arrive at
> > the average.
> > > Heterogeneous systems often take longer to equilibrate.
> > >
> > > -Justin
> > >
> > > > Thanks a lot.
> > > > Guang-Jun
> > > >
> > > >
> > ------------------------------------------------------------------------
> > > > 立刻下载 MSN 保护盾,保障 MSN 安全稳定! 现在就下载!
> > > > <http://im.live.cn/safe/>
> > > >
> > > >
> > > >
> > ------------------------------------------------------------------------
> > > >
> > > > _______________________________________________
> > > >! ! ; gmx-users mailing list gmx-users@gromacs.org
> > > > http://www.gromacs.org/mailman/listinfo/gmx-users
> > > > Please search the archive at http://www.gromacs.org/search before
> > posting!
> > > > Please don't post (un)subscribe requests to the list. Use the
> > > > www interface or send it to gmx-users-requ...@gromacs.org.
> > > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
> > >
> > > --
> > > ========================================
> > >
> > > Justin A. Lemkul
> > > Graduate Research Assistant
> > > Department of Biochemistry
> > > Virginia Tech
> > > Blacksburg, VA
> > > jalemkul[at]vt.edu | (540) 231-9080
> > > http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
> > >
> > > ========================================
> > > _______________________________________________
> > > gmx-users mailing list gmx-users@gromacs.org
> > > http://www.gromacs.org/mailman/listinfo/gmx-users
> > > Please! ! search the archive at http://www.gromacs.org/search before
> > posting!
> > > Please don't post (un)subscribe requests to the list. Use the
> > > www interface or send it to gmx-users-requ...@gromacs.org.
> > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
> >
> > ------------------------------------------------------------------------
> > 你和你的爱人都去过哪里?快来微软地图一起绘制爱的地图吧! 立即查看!
> > <http://ditu.live.com/?form=MICHAJ >
>
> --
> ========================================
>
> Justin A. Lemkul
> Graduate Research Assistant
> Department of Biochemistry
> Virginia Tech
> Blacksburg, VA
> jalemkul[at]vt.edu | (540) 231-9080
> http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
>
> ========================================
> _______________________________________________
> gmx-users mailing list gmx-users@gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-requ...@gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
_________________________________________________________________
上Windows Live 中国首页,下载最新版 MSN!
http://im.live.cn/
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php