24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kevin8871

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 28.5
帖子: 8
在线: 4.1小时
虫号: 1075839

[交流] 【求助】用NAMD 做NPT模拟遇到下面的问题

本人用NAMD做NPT模拟，遇到了下面的问题
FATAL ERROR: Periodic cell has become too small for original patch grid!
Possible solutions are to restart from a recent checkpoint,
increase margin, or disable useFlexibleCell for liquid simulation.
FATAL ERROR: Periodic cell has become too small for original patch grid!
Possible solutions are to restart from a recent checkpoint,
increase margin, or disable useFlexibleCell for liquid simulation.
这是为什么？下面是我的conf file
structure       3HT_20.psf
coordinates    3HT_20_centered.pdb

set outputname    NPT file for 3 HT

#continuing a run
set inputname    nimi
binVelocities    3HT_20_mini.restart.vel
bincoordinates    3HT_20_mini.restart.coor
extendedsystem    3HT_20_mini.restart.xsc
firsttimestep    6000
numsteps          2506000

#############################################################
## SIMULATION PARAMETERS                                  ##
#############################################################

# Input

paraTypeCharmm       on
parameters          ht_monomer.par
# Force-Field Parameters
exclude          scaled1-4
1-4scaling       1.0
cutoff             12.
switching          on
switchdist       10.
pairlistdist       13.5

# Integrator Parameters
timestep          2.0  ;# 2fs/step
rigidBonds       all  ;# needed for 2fs steps
nonbondedFreq    1
fullElectFrequency  2
stepspercycle    10

# Constant Temperature Control
langevin          on ;# do langevin dynamics
langevinDamping    5    ;# damping coefficient (gamma) of 5/ps
langevinTemp       300
langevinHydrogen off ;# don't couple langevin bath to hydrogens

#PBC
wrapall    on
cellBasisVector1 31.312999725341797 0.0 0.0
cellBasisVector2    0.0 32.558000564575195 0.0
cellBasisVector3    0.0 0 45.066001892089844
cellOrigin       -3.0093722216406604e-6 1.4230939996195957e-5 2.136936927854549e-6

wrapAll          on

# PME (for full-system periodic electrostatics)
PME                yes
PMEGridSizeX       42
PMEGridSizeY       43
PMEGridSizeZ       55

# Constant Pressure Control (variable volume)
useGroupPressure    yes ;# needed for rigidBonds
useFlexibleCell    no  ;# no for water box yes for membrane
useConstantArea    no  ;# no for waterbox, maybe for membrane

langevinPiston       on
langevinPistonTarget  1.01325 ;#  in bar -> 1 atm
langevinPistonPeriod  100.
langevinPistonDecay 50.
langevinPistonTemp 300

# Output
outputName       $outputname

restartfreq       2500    ;# 2500steps = every 5ps
dcdfreq          1000    ;# 1000steps = every 2ps
xstFreq          1000
outputEnergies    500
outputPressure    500

#############################################################
## EXTRA PARAMETERS                                     ##
#############################################################

#############################################################
## EXECUTION SCRIPT                                     ##
#############################################################

# Minimization

run 2500000 ;# 5ns

回复此楼

» 猜你喜欢

科研求助已经有3人回复
面上没中，邀请各位路过的虫友分析一下分数已经有15人回复
博士去二本高校当辅导员，如何调整心态？已经有10人回复
面上有专家说实验设备不是我们单位的已经有6人回复
27细胞生物学申博已经有5人回复
江苏省自然基金什么时候出结果已经有6人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有11人回复
闲聊已经有3人回复
两块石头已经有10人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有15人回复

» 本主题相关价值贴推荐，对您同样有帮助:

关于超级计算机的安装权限问题：不能用sudo命令，该如何是好？已经有18人回复
word参考文献-尾注引用，出了点问题，求高手解答已经有3人回复
签证时用身份证的问题已经有12人回复
1,3-二羰基化合物单碘代的问题请做过类似反应的指教下已经有5人回复
【求助】在NAMD运行中使用CHARMM格式的力场文件，但力场参数来自AMBER力场，会出问题� 已经有4人回复
【求助】做NAMD的1K4C.pdb时候出现的错误已经有30人回复
【求助】Gromacs和NAMD比Amber快很多吗？如何选择？已经有7人回复
【专家会诊】bay 讨论namd,gmx处理水溶液,蛋白等MD 模拟;讨论水及水通道已经有50人回复
【求助】NAMD和Gromacs哪个简单？已经有5人回复
【求助】namd tutorial遇到问题已经有10人回复
【求助】NAMD安装问题已经有19人回复
【在线答疑】关于NAMD/VMD 已经有210人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-11-28 15:14:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wally8962

木虫 (著名写手)

应助: 9 (幼儿园)
金币: 5366.7
帖子: 1062
在线: 471.1小时
虫号: 600024

★ ★ ★
zh1987hs(金币+3):thank you 2010-12-03 17:01:47
kevin8871(金币+2): 2010-12-05 09:27:13

langevinPistonPeriod 100.
langevinPistonDecay 50.

太小。改大。

ug建议的是200、100.实际上可能需要更大。

赞一下(1人)

回复此楼

4楼2010-12-03 16:52:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 11 个回答

leewz922

金虫 (小有名气)

应助: 6 (幼儿园)
金币: 1141.8
帖子: 101
在线: 89.8小时
虫号: 609327

★ ★ ★
zh1987hs(金币+3):谢谢 2010-12-02 22:11:55

把下面这几行去掉，试试
cellBasisVector1 31.312999725341797 0.0 0.0
cellBasisVector2    0.0 32.558000564575195 0.0
cellBasisVector3    0.0 0 45.066001892089844
cellOrigin       -3.0093722216406604e-6 1.4230939996195957e-5 2.136936927854549e-6

赞一下(1人)

回复此楼

2楼2010-12-02 20:55:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖