[gmx-users] Pressure Question
stephane.teletchea at univ-nantes.fr
Wed Nov 5 20:54:14 CET 2014
Le 04/11/2014 18:00, Johnny Lu a écrit :
> If my NVT simulation of a protein in 30k molecules of water has a pressure
> of 11 bar (error 0.5 bar from g_energy), will the dynamics (not
> distribution of conformations) change enough that the mechanism inferred
> from this simulation be significantly more unreliable than the mechanism
> inferred from a 1 bar simulation? (Will the reviewers cut my paper into
> Thanks again.
Considering only your "NVT" parameters for your simulation,
I would consider it totally irrelevant to talk about "pressure" where
your constrain the volume.
This value or any other one has not really a meaning in this situation,
and I seen
many variations in the pressure value in this microcanonical ensemble
without paying too much
attention on it.
In an "NPT" simulation, then you should be able to find back a normal 1
bar simulation I think.
Do you have any reason to do first an NPT simulation, and "then" an NVT one?
I would personally let the system equilibrate in NVT, then swith to the
more natural NPT,
provided actual code and force fields are now good enough in this ensemble.
Team Protein Design In Silico
UFIP, UMR 6286 CNRS,
UFR Sciences et Techniques,
2, rue de la Houssinière, Bât. 25,
44322 Nantes cedex 03, France
Tél : +33 251 125 636
Fax : +33 251 125 632
http://www.ufip.univ-nantes.fr/ - http://www.steletch.org
More information about the gromacs.org_gmx-users