[gmx-users] error when using g_membed tool
Kieu Thu Nguyen
kieuthu2212 at gmail.com
Tue Jan 8 05:29:18 CET 2013
Dear All,
I want to use g_membed tool to embed protein into lipids membrane. When i
type the command
g_membed -f input.tpr -p merged.top -xyinit 0.1 -xyend 1.0 -nxy 1000
the error is
Step 1 Warning: pressure scaling more than 1%, mu: 21.3036 21.3036 1.51453
Segmentation fault (core dumped)
My sample.mdp file for making the input.tpr file is
; TIMESTEP IN MARTINI
integrator = md
tinit = 0.0
dt = 0.02
nsteps = 1000
nstcomm = 1
comm-grps =
nstxout = 100
nstvout = 100
nstfout = 0
nstlog = 100
nstenergy = 100
nstxtcout = 100
xtc_precision =
xtc-grps =
energygrps = Protein
freezegrps = Protein
freezedim = Y Y Y
energygrp_excl = Protein Protein
; NEIGHBOURLIST and MARTINI
nstlist = 10
ns_type = grid
pbc = xyz
rlist = 1.4
; MARTINI and NONBONDED
coulombtype = Shift
rcoulomb_switch = 0.0
rcoulomb = 1.2
epsilon_r = 15
vdw_type = Shift
rvdw_switch = 0.9
rvdw = 1.2
DispCorr = No
; MARTINI and TEMPRATURE/PRESSURE
tcoupl = Berendsen
tc-grps = Protein DSPC Ion_W
tau-t = 2.0 2.0 2.0
ref-t = 300 300 300
pcoupl = Berendsen
pcoupltype = semiisotropic
tau-p = 4.0 4.0
compressibility = 3e-4 3e-4
ref-p = 1.0 1.0
; GENERATE VELOCITIES FOR STARTUP RUN
gen-vel = no
; MARTINI and CONSTRAINTS
constraints = none
constraint_algorithm = Lincs
continuation = no
lincs_order = 4
lincs_warnangle = 30
I inreased tau_p value as the gromacs website guiding but this error still
exists.
What should i do ?
Thanks in advance !
KT
More information about the gromacs.org_gmx-users
mailing list