[gmx-users] still problem of pr.mdp

sudheer babu sudheer.pbm07 at gmail.com
Tue Feb 26 11:18:29 CET 2008


Hi gmxusers,
i have used below parameters those what you suggested ,

title               =   popc restrained
define              =  -DPOSRES
constraints         =  all-bonds
integrator          =  md
dt                  =  0.002    ; ps !
nsteps              =  10000    ; total 20 ps.
nstcomm             =  1
nstxout             =  50
nstvout             =  1000
nstfout             =  0
nstlog              =  10
nstenergy           =  10
nstlist             =  10
ns_type             =  grid
rlist               =  0.9
coulombtype         =  PME
rcoulomb            =  0.9
rvdw                =  1.4
pbc                 =  xyz
; Berendsen temperature coupling is on in two groups
Tcoupl              =  Berendsen
tc-grps             =  POPC   Protein   SOL_Cl
tau_t               =  0.1     0.1      0.1
ref_t               =  310     310      310
; Anisotropic pressure coupling is now on
Pcoupl              =  parrinello-rahman
pcoupltype          =  anisotropic
tau_p               =  1.0      1.0     1.0     0     0     0
compressibility     =  4.5e-5   4.5e-5  4.5e-5  0       0       0
ref_p               =  1.0      1.0     1.0     0       0       0
; Energy monitoring
energygrps          =  POPC  Protein  SOL_Cl
; Generate velocites is on at 300 K.
gen_vel             =  yes
gen_temp            =  310
gen_seed            =  173529
~

~

 but still i am getting "tilting structure of POPC" except water and protein
stucture( fine) . I am using protein restrain force 1000 and my system box
size is 7nm.
Pls help me.......
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20080226/028ec1ff/attachment.html>


More information about the gromacs.org_gmx-users mailing list