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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【求助】用NAMD 做NPT模拟遇到下面的问题
111/1返回列表
查看: 2370  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂 +6 吃吃吃才有意义 2026-03-19 6/300 2026-03-20 10:47 by 尽舜尧1
[考研] 296求调剂 +3 www_q 2026-03-18 7/350 2026-03-20 10:39 by barlinike
[考研] 295复试调剂 +5 简木ChuFront 2026-03-19 5/250 2026-03-20 10:38 by 无际的草原
[考研] 265求调剂 +6 梁梁校校 2026-03-17 6/300 2026-03-20 08:59 by ZHANG0tao
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +8 步川酷紫123 2026-03-13 8/400 2026-03-20 08:47 by xingguangj
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 10/500 2026-03-19 16:05 by 余麟余
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7 小材化本科 2026-03-18 7/350 2026-03-19 10:46 by Linda Hu
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
信息提示
关闭
请填处理意见
关闭