[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