Hi,
You will have to come up yourself with a reasonable guess.
If you have a hard crystal, I guess it would be close to the diagonal values.
Note that with PR coupling tau_p has to be large enough, something like 10 ps.
(note that in Gromacs 4.0 a bug has been corrected which scaled tau_p by 16.6
in version 3.3 and before)
But for equilibration I would start with Berendsen.
Berk
Date: Thu, 19 Mar 2009 04:30:44 -0700
From: lastexile...@yahoo.de
Subject: RE: [gmx-users] compressibility tensor components, pressure
coupling anisotropic PR, triclinic systems
To: gmx-users@gromacs.org
Hello,
I do have a crystal system. You refer to an initial guess. So if I understand
correctly, it is a trial and error, if I do not have any data, this is what I
can make out of. The part of the shear stress you are reffering I think that it
is from experimental values?
Now you say that : compressibility (real) = tau_p*compressibility (gromacs) ?
If this is the case I do not think that the units coincide.
(ps/bar) =? (1/bar)
I think that what you mean has to do with equation 3.36 in the manual, where
with these components the W is calculated.
Forgive my ignorance in the issue.
Thank you,
Nikos
--- Berk Hess <g...@hotmail.com> schrieb am Do, 19.3.2009:
Von: Berk Hess <g...@hotmail.com>
Betreff: RE: [gmx-users] compressibility tensor components, pressure coupling
anisotropic PR, triclinic systems
An: "Discussion list for GROMACS users" <gmx-users@gromacs.org>
Datum: Donnerstag, 19. März 2009, 10:23
Hi,
If you have a liquid system (no off-diagonal elasticity), you should not use
full anistropic pressure coupling.
If you have a solid system (for instance a crystal), you will have an elastic
shear stress response
and you can determine and use off-diagonal compressibility values.
You need the compressibility values to set the time scale of the pressure
coupling.
The compressibility does not affect any thermodynamic quantity, only the
dynamics.
The compressibility always occurs as a product with 1/tau_p.
If you want to determine them from a simulation, you need a reasonable initial
guess
for the compressibility and you can use a large tau_p to be on the safe side
with the actual coupling time.
Berk
Date: Wed, 18 Mar 2009 16:02:38 -0700
From: lastexile...@yahoo.de
To: gmx-users@gromacs.org
Subject: [gmx-users] compressibility tensor components, pressure coupling
anisotropic
PR, triclinic systems
Hello,
Searching first the gromacs mailing list I could not find an answer to the
problem I face.
I would like to know the vaules I have to give to the mdp file where it asks
for compressibility.
I have to conduct an NPT simulation using barostat Parrinello-Rahman. My
pressure coupling should be anisotropic. The values I have given up to now are:
tau_p = 1.0
compressibility = 4.5e-5 4.5e-5 4.5e-5 4.5e-5 4.5e-5 4.5e-5
ref_p = 1.0 1.0 1.0 0.0 0.0 0.0
when I try to calculate the elastic constants I get reasonable values for the
diagonal elements, yet not
good for the off-diagonal. The system I have is a monoclinic one, so giving
zeros to xy,xz,yz components is not an option.
An experimental linear compressibility value I have found out is very close to
the value above, so this is why I think I get good values in the diagonal
elements.
I have found out many different questions, including one about monoclinic
systems. In one of them someone was wondering why he had to put as input the
compressibility values, while he is supposed to calculate them from the
simulation. As far as the monoclinic system is concerned the question had to to
with the components in gromacs derived from the linear experimental values...
There was no answer to that.
I have tried some different values myself, for the off - diagonal components
and I get either a crash of my simulation of different results for the elastic
constants. Is there an algorithm which clearly can point me to the right
direction on how to
derive
the right values?
Thank you in advance,
Nikos
What can you do with the new Windows Live? Find out
_______________________________________________
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
_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/
_______________________________________________
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
