[gmx-users] pressure scaling
dinter at rc.m-kagaku.co.jp
Tue Aug 3 07:51:24 CEST 2004
I used a tau_p = 10, then I tried different values for compressibility:
1st:4.5e-5 (initial crash), 14.5e-5, 0.5e-5, 0.5e-6, 0.5e-4 and 4.5e-3. But
there isn't many changes, the cash occured allmost at the same step: 303
(4.5e-5),291 (14.5e-5),301 (0.5e-5),301 (0.5e-6),299(0.5e-4) and 293
Ok so if I undestood well the pressure scaling error occurs when the box
enter in resonance, but normaly with low compressibility it must not occur
or at least at different steps (more difficult to change the box) especialy
with low pressure.
So I would like to know which are the parameter which drive this error:
there is the time step (I took:0.0005ps), the tau_p parameter (I took 10ps)
and the compressibility. And what else?
I still not well undestand what does represent tau_p is my draw is a correct
understanding? (see atached file).
Thank for any help.
From: gmx-users-bounces at gromacs.org
[mailto:gmx-users-bounces at gromacs.org]On Behalf Of Berk Hess
Sent: lundi 2 août 2004 18:18
To: gmx-users at gromacs.org
Subject: Re: [gmx-users] pressure scaling
A smaller tau_p is not the solution.
If anything, it should be larger.
What matters in Berendsen pressure coupling is the ratio:
So even if tau_p seems reasonable, the scaling may go wrong
due to a too large compressibility value.
I don't know which value you have used and which value you
should use for your system.
> Try to use lower tau_p (maybe 3ps?).
>On Aug 2, 2004, at 9:31 AM, nicolas Dinter wrote:
>>Dear Gromacs users,
>>I am performing simulation on diphenyl carbonate and I get the presure
>>I used a small pression about 1 bar, tau_p = 10 ps and my time scale is
>>The NVT simulation runs well, to solve this problem for hight pressure is
>>easier but low pressure knowing that this error is related with tau_p and
>>timestep I would like to know how to fix it.
Talk with your online friends with MSN Messenger http://messenger.msn.nl/
gmx-users mailing list
gmx-users at gromacs.org
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request at gromacs.org.
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 3864 bytes
Desc: not available
More information about the gromacs.org_gmx-users