[gmx-users] pressure scaling more than 1%
Sebastian Waltz
sebastian.waltz at physik.uni-freiburg.de
Fri Jun 11 10:13:37 CEST 2010
Hi all,
my simulation runs well for the first 700000 steps with a
step size of 0.001ps and when I get the warning that the
pressure scaling is more than 1% every 10000 step or so.
So, is my tau_p of 0.2 to small or is the equillibration
not done well enough.
I am simulating a small peptide in CHCl3 with an handmade
CHCl3 model.
title = ttt
cpp = /lib/cpp
include = -I../top
constraints = none
;define = -DFLEX_SPC
;define = -DPOSRES ; for possition
restraints
integrator = md
emtol = 100.0
emstep = 0.005
dt = 0.001 ; ps !
nsteps = 2000000 ; total 2 ns
nstcomm = 1
nstxout = 0
nstvout = 0
nstfout = 10
nstlog = 0
nstenergy = 10
nstxtcout = 10
xtc_grps = Protein+PTS+AIC
nstlist = 5
ns_type = grid
rlist = 1.2
coulombtype = PME
rcoulomb = 1.2
rvdw = 1.2
fourierspacing = 0.12
pme_order = 4
optimize_fft = yes
Tcoupl = berendsen
tc-grps = Protein+PTS+AIC CHCl3
tau_t = 0.1 0.1
ref_t = 300 300
energygrps = Protein+PTS+AIC CHCl3
Pcoupl = berendsen
tau_p = 0.2
compressibility = 5.4e-5
ref_p = 1.0
gen_vel = yes
gen_temp = 300
gen_seed = 913529
Thanks
More information about the gromacs.org_gmx-users
mailing list