24小时热门版块排行榜    

查看: 1894  |  回复: 1

Estheriii

新虫 (初入文坛)

[交流] 在gromsacs里加电场后就mdrun不起来了 已有1人参与

做的是dmf的模拟里面放有1mol nacl ,下面附上我的md.mdp信息,求大神帮助!!
title                = OPLS Lysozyme MD
; Run parameters
integrator        = md                ; leap-frog integrator
nsteps                = 250000       
dt                = 0.002                ; 2 fs
; Output control
nstxout                = 500                ; save coordinates every 2 ps
nstvout                = 500                ; save velocities every 2 ps
nstxtcout        = 500                ; xtc compressed trajectory output every 2 ps
nstenergy        = 500                ; save energies every 2 ps
nstlog                = 500                ; update log file every 2 ps
; Bond parameters
continuation        = yes                ; Restarting after NPT
constraint_algorithm = lincs        ; holonomic constraints
constraints        = all_bonds                ; all bonds (even heavy atom-H bonds) constrained
lincs_iter        = 1                ; accuracy of LINCS
lincs_order        = 4                ; also related to accuracy
; Neighborsearching
ns_type                = grid                ; search neighboring grid cells
nstlist                = 5                ; 10 fs
rlist                = 0.8                ; short-range neighborlist cutoff (in nm)
rcoulomb        = 0.8                ; short-range electrostatic cutoff (in nm)
rvdw                = 0.8                ; short-range van der Waals cutoff (in nm)
; Electrostatics
coulombtype        = Ewald                ; Particle Mesh Ewald for long-range electrostatics
pme_order        = 4                ; cubic interpolation
fourierspacing        = 0.16                ; grid spacing for FFT
; Temperature coupling is on
tcoupl                = V-rescale        ; modified Berendsen thermostat
tc-grps                = MOL        SOL  NA  CL        ; two coupling groups - more accurate
tau_t                = 0.2        0.2  0.2  0.2        ; time constant, in ps
ref_t                = 298         298  298  298        ; reference temperature, one for each group, in K
; Pressure coupling is on
pcoupl                = Parrinello-Rahman        ; Pressure coupling on in NPT
pcoupltype        = isotropic        ; uniform scaling of box vectors
tau_p                = 2.0                ; time constant, in ps
ref_p                = 1.0                ; reference pressure, in bar
compressibility = 4.5e-5        ; isothermal compressibility of water, bar^-1
; Periodic boundary conditions
pbc                = xyz                ; 3-D PBC
; Dispersion correction
DispCorr        = EnerPres        ; account for cut-off vdW scheme
; Velocity generation
gen_vel                = no                ; Velocity generation is off
nh-chain-length = 1  
E-x             = 1 3e5 0      

出现的错误是
step 0: Water molecule starting at atom 619 can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.

Back Off! I just backed up step0b.pdb to ./#step0b.pdb.4#

Back Off! I just backed up step0c.pdb to ./#step0c.pdb.4#
Wrote pdb files with previous and current coordinates
step 0
WARNING: Listed nonbonded interaction between particles 6 and 2
at distance 3f which is larger than the table limit 3f nm.

This is likely either a 1,4 interaction, or a listed interaction inside
a smaller molecule you are decoupling during a free energy calculation.
Since interactions at distances beyond the table cannot be computed,
they are skipped until they are inside the table limit again. You will
only see this message once, even if it occurs for several interactions.

IMPORTANT: This should not happen in a stable simulation, so there is
probably something wrong with your system. Only change the table-extension
distance in the mdp file if you are really sure that is the reason.



step 1: Water molecule starting at atom 478 can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.

Back Off! I just backed up step1b.pdb to ./#step1b.pdb.1#

Back Off! I just backed up step1c.pdb to ./#step1c.pdb.1#
Wrote pdb files with previous and current coordinates
段错误 (核心已转储)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

japmail

银虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
我也遇见同意的问题了,请问您解决了,怎么解决的
我还遇见,做npt的时候遇见的,nvt,em都没问题啊

WARNING: Listed nonbonded interaction between particles 1 and 13
at distance 2.504 which is larger than the table limit 2.016 nm.

This is likely either a 1,4 interaction, or a listed interaction inside
a smaller molecule you are decoupling during a free energy calculation.
Since interactions at distances beyond the table cannot be computed,
they are skipped until they are inside the table limit again. You will
only see this message once, even if it occurs for several interactions.

IMPORTANT: This should not happen in a stable simulation, so there is
probably something wrong with your system. Only change the table-extension
distance in the mdp file if you are really sure that is the reason.


step 0Warning: Triclinic box is too skewed.
         Box (3x3):
            Box[    0]={-8.19631e-01,  0.00000e+00,  0.00000e+00}
            Box[    1]={-0.00000e+00, -8.19631e-01, -0.00000e+00}
            Box[    2]={-0.00000e+00, -0.00000e+00, -8.19631e-01}
         Can not fix pbc.
Warning: Triclinic box is too skewed.
         Box (3x3):
            Box[    0]={-8.19631e-01,  0.00000e+00,  0.00000e+00}
            Box[    1]={-0.00000e+00, -8.19631e-01, -0.00000e+00}
            Box[    2]={-0.00000e+00, -0.00000e+00, -8.19631e-01}
         Can not fix pbc.
段错误 (核心已转储)
一个高尚而又纯粹的爷
2楼2015-09-04 16:21:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Estheriii 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 27博士 +6 锦鲤送福至 2026-11-10 6/300 2026-11-14 21:23 by xupeii
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-14 15:58 by KXypaGgBBlw7
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-13 11:37 by Equinoxhua
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1Nb0IVlmV6v0 2026-11-11 4/200 2026-11-13 11:18 by Equinoxhua
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:11 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 11:04 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 9/450 2026-11-13 11:02 by Equinoxhua
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:28 by ECLTXnJnRSXf
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:22 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:00 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 C6CRn66RYYXS 2026-11-10 7/350 2026-11-12 10:16 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 10:03 by HajghKOMynmz
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 3/150 2026-11-12 08:26 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 06:03 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:30 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
信息提示
请填处理意见