[gmx-users] pressure and density

Alexander Alexander alexanderwien2k at gmail.com
Sat Jun 25 17:27:04 CEST 2016


Dear gromacs user,

I have a system containing a solid surface(18000 atoms) and on top of it a
heptapeptide solvated in water(SOL  9604), after a good minimisation and
NVT equilibration and before production(which will be in NVT ensemble), I
did NPT equilibration (1 bar) for 200 ps and continued to 1 ns as well, the
average pressure results are coming below.
-----------------------
pcoupl                  = Parrinello-Rahman
Pcoupltype          = isotropic             ;;;;semiisotropic
compressibility    = 4.46e-5             ;;;;0.0   4.46e-5
ref_p                      = 1.0                        ;;;1.0    1.0
tau_p                     = 2.0
refcoord-scaling   = com
------------------------

Statistics over 100001 steps [ 0.0000 through 200.0000 ps ], 1 data sets
All statistics are over 1001 points

Energy                      Average   Err.Est.       RMSD  Tot-Drift
-------------------------------------------------------------------------------
Pressure                   -782.581         39    925.396   -138.729  (bar)


Statistics over 500001 steps [ 0.0000 through 1000.0000 ps ], 1 data sets
All statistics are over 5001 points

Energy                      Average   Err.Est.       RMSD  Tot-Drift
-------------------------------------------------------------------------------
Pressure                   -746.168         18    921.982    3.41813  (bar)


My first question as always is that if these amount of pressure and their
RMSD together are satisfying the reference pressure of 1 bar that I
considered?

The second question: The water density I got is around 920 kg/m^3 in the
middle of water region and far away from the solid surface, while I was
expecting to get something around 997 kg/m^3 in the region far from the
solid surface. And none of the tools like gmx inset-molecule, gmx solvate
... , does not let me add extra water molecule to the system, then, how I
can increase the density of water in that region?

Thanks in advance for your comments.

Regards,
Alex


More information about the gromacs.org_gmx-users mailing list