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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【求助】用NAMD 做NPT模拟遇到下面的问题
111/1返回列表
查看: 2348  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 326求调剂 +3 上岸的小葡 2026-03-15 4/200 2026-03-15 18:50 by 无际的草原
[基金申请] 国自科面上基金字体 +4 iwuli 2026-03-12 5/250 2026-03-15 17:07 by 风云无泪
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 337一志愿华南理工材料求调剂（有希望2吗？） +3 mysdl 2026-03-09 3/150 2026-03-14 02:53 by JourneyLucky
[考研] 293求调剂 +5 上班不着吉 2026-03-09 5/250 2026-03-14 02:37 by JourneyLucky
[考研] 一志愿郑大070303，338分，求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 一志愿安徽大学材料工程专硕313分，求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 298求调剂 +3 Vv呀！ 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
信息提示
关闭
请填处理意见
关闭