[gmx-users] Meaningless pressure values: NVT and NVE ensembles

[Machtens, Jan-Philipp]

Hi,
although this is probably not the solution to your problem, you should update to Gromacs 5.0.2
(since 5.0.1 has a bug which could affect the quantities you calculate --> http://redmine.gromacs.org/issues/1603 ).
Cheers, Jan-Philipp
________________________________________
From: gromacs.org_gmx-users-bounces at maillist.sys.kth.se [gromacs.org_gmx-users-bounces at maillist.sys.kth.se] on behalf of Francesco Mambretti [francesco.mambretti at studenti.unimi.it]
Sent: Wednesday, November 05, 2014 12:14 PM
To: gmx-users at gromacs.org
Subject: [gmx-users] Meaningless pressure values: NVT and NVE ensembles

Dear all, I am using Gromacs 5.0.1 to simulate a simple atomic system of 108 solid LJ argon particles in the void, pbc applied.
Using g_energy, I find - both after NVT runs and NVE runs - ok values for energies and temperatures, but absurd values for pressures (something like 10^3 bar or -10^3 bar). I read lots about it on the web, included an old post on http://article.gmane.org/gmane.science.biology.gromacs.user/6696/match=negative+pressure; what I am not still able to understand is if output values for pressure computed by Gromacs are meaningful.
I am using the right input parameters (density, temperature and so on) - with other simple MD programs they work well.
Moreover, if I understand well, a negative pressure means that the system wants to collapse, but this doesn't make sense for my system, if it is already at the right density!
So, should I trust pressure output values from g_energy? Does Gromacs compute pressure in a way not compatible with NVE and NVT?
If not, how can I check the pressure of my microcanonical or canonical system?

I am applying long-range dispersion corrections, I generate velocities from the temperature; timestep = 1 fs, no pressure coupling, potential has a cutoff and it is shifted to 0 in correspondence of cutoff radius.

Thanks,


Francesco Mambretti
--
Gromacs Users mailing list

* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.


------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------
Forschungszentrum Juelich GmbH
52425 Juelich
Sitz der Gesellschaft: Juelich
Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498
Vorsitzender des Aufsichtsrats: MinDir Dr. Karl Eugen Huthmacher
Geschaeftsfuehrung: Prof. Dr.-Ing. Wolfgang Marquardt (Vorsitzender),
Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt,
Prof. Dr. Sebastian M. Schmidt
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------