[gmx-users] The Cut-off for coulombtype heat up the water system?
Yanmei Song
ysong30 at asu.edu
Thu Jun 18 02:13:40 CEST 2009
Dear Users:
I was running a system by non-equilibrium MD using a plain Cut-off for the
electrostatics:
title = water
cpp = /lib/cpp
constraints = all_bonds
integrator = md
dt = 0.004 ; ps !
nsteps = 500000 ; total 8ns.
nstcomm = 1
nstxout = 50000
nstvout = 50000
nstfout = 0
nstlog = 5000
nstenergy = 5000
nstxtcout = 25000
nstlist = 10
ns_type = grid
pbc = xyz
coulombtype = Cut-off
rlist = 0.8
rcoulomb = 0.9
rvdw = 0.8
fourierspacing = 0.12
ewald_rtol = 1e-5
;nemd
cos_acceleration = 0.005
; Berendsen temperature coupling is on in two groups
Tcoupl = berendsen
tc_grps = SOL
tau_t = 0.1
ref_t = 300
; Energy monitoring
energygrps = SOL
; Isotropic pressure coupling is now on
Pcoupl = berendsen
pcoupltype = isotropic
;pc-grps = SOL
tau_p = 1.0
ref_p = 1.0
compressibility = 4.5e-5
; Generate velocites is off at 300 K.
gen_vel = yes
gen_temp = 300.0
gen_seed = 100000
I can not figure out where I did wrong, the temperature of the system is
303.54 after 5ns run ( the temperature turns to 303 in 500ps). Thanks for
the help in advance!
--
Yanmei Song
Ph.D. Candidate
Department of Chemical Engineering
Arizona State University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090617/faf5bbcf/attachment.html>
More information about the gromacs.org_gmx-users
mailing list