应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【求助】用NAMD 做NPT模拟遇到下面的问题
111/1返回列表
查看: 2349  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kevin8871

新虫 (初入文坛)

[交流] 【求助】用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
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-11-28 15:14:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leewz922

金虫 (小有名气)
★ ★ ★
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的回帖

kevin8871

新虫 (初入文坛)
我已经去掉了，我先做了NVT，然后用NVT的做NPT，但是我做到后面的时候出错了，我查看了下，我的box的边长变成了16 16 25， 但是我的cutoff 是12，我觉得边长应该大于cutoff的两倍，所以最少应该是24，请问你有什么办法修改这个嘛，我现在是不知道怎么弄了。
赞一下 回复此楼
3楼2010-12-03 14:30:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)
★ ★ ★
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的回帖

leewz922

金虫 (小有名气)
★ ★ ★
zh1987hs(金币+3):谢谢 2010-12-03 22:49:09
kevin8871(金币+2): 2010-12-05 09:27:20
PME 参数设置应该被 2 、3、 5整除 且稍大于box边长

改成这样试试
PMEGridSizeX        40
PMEGridSizeY        40
PMEGridSizeZ        50
赞一下(1人) 回复此楼
5楼2010-12-03 20:07:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin8871

新虫 (初入文坛)
★
zh1987hs(金币+1):呵呵~鼓励交流 2010-12-05 10:46:52
我老板的建议是不要PME，而用PMEGridSpacing 1.0代替，我准备按照wally8962  的建议试一下，谢谢大家的建议
赞一下(1人) 回复此楼
6楼2010-12-05 09:26:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin8871

新虫 (初入文坛)
按照wally的建议改了下,多跑了很多,但是到最后同一个问题又出现了,请问应该怎么解决?
赞一下 回复此楼
7楼2010-12-05 13:36:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)
引用回帖:
Originally posted by kevin8871 at 2010-12-05 13:36:13:
按照wally的建议改了下,多跑了很多,但是到最后同一个问题又出现了,请问应该怎么解决?

继续改大。
有的时候需要改到1000、500
赞一下 回复此楼
8楼2010-12-05 14:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin8871

新虫 (初入文坛)
我是用NVT的结果做NPT,是NVT的conf file也要改到200 100是吧
赞一下 回复此楼
9楼2010-12-05 15:09:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin8871

新虫 (初入文坛)
噢,看错了,NVT的不用,我刚把NPT的改成了300 和150,跑了两步就死了,我看了xst文件第一布的时候就小于24了,这是为什么
赞一下 回复此楼
10楼2010-12-05 15:11:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ruthxu

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by kevin8871 at 2010-12-05 15:11:22
噢,看错了,NVT的不用,我刚把NPT的改成了300 和150,跑了两步就死了,我看了xst文件第一布的时候就小于24了,这是为什么

你好，我在用namd模拟npt下的磷脂膜，但是也遇到了Periodic cell has become too small for original patch grid的问题，看到你发过帖子，想问一下，你是怎么解决的？非常感谢!
赞一下 回复此楼
11楼2014-08-05 10:53:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kevin8871 的主题更新
111/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 4/200 2026-03-16 21:52 by Losir
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[教师之家] 焦虑 +7 水冰月月野兔 2026-03-13 9/450 2026-03-16 10:00 by Quakerbird
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 一志愿华中农业大学071010，总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子！
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭