Fw: [gmx-users] Reg:1 particles communicated to PME node 0 are more than 2/3 times the cut-off
jalemkul at vt.edu
Sun Feb 10 21:18:25 CET 2013
On 2/10/13 3:14 PM, Shima Arasteh wrote:
> I get the same error, but the information written in my log file is as :
> Energies (kJ/mol)
> U-B Proper Dih. Improper Dih. CMAP Dih. LJ-14
> 8.26180e+04 5.53024e+04 6.50084e+02 -2.83474e+02 8.50267e+03
> Coulomb-14 LJ (SR) Disper. corr. Coulomb (SR) Coul. recip.
> -1.06730e+05 7.84233e+04 -7.99529e+03 -9.89064e+05 -2.54063e+04
> Potential Kinetic En. Total Energy Temperature Pres. DC (bar)
> -9.03982e+05 3.53794e+05 -5.50188e+05 4.26000e+02 -1.34125e+02
> Pressure (bar) Constr. rmsd
> -1.97048e+04 5.76582e-04
> The temperature is around 400, it seems normal, but the pressure is negative! Is it the problem which causes my system get this fatal error?
An instantaneous value of pressure is meaningless.
> I changed the pressure coupling to Berendsen and decrease the time step to 1fs. How would I be sure that my modifications is beneficial? By checking the pressure in log file passing the time?
Well, for one, the system won't blow up ;)
The assessment of any equilibration is basically the same - do the observables
of interest (temperature, pressure, etc) converge to the desired value?
Justin A. Lemkul, Ph.D.
Department of Biochemistry
jalemkul[at]vt.edu | (540) 231-9080
More information about the gromacs.org_gmx-users