24小时热门版块排行榜    

查看: 3339  |  回复: 10

574384607

银虫 (小有名气)

[交流] 【求助】namd 出错Periodic cell has become too small for original patch grid! 已有3人参与

体系在水盒子里做模拟,运行过程中,
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.
好像是周期性水盒子变小了,这样的话把cellBasisVector调大些预计可以处理掉吗?
还是别的地方出了问题哪?
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

superdirac

木虫 (正式写手)

★ ★
ghcacj(金币+2): 谢谢 2011-03-23 12:50:33
wrapAll                 on
后面再加一行。
margin                  3
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
8楼2011-03-22 15:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

574384607

银虫 (小有名气)

引用回帖:
Originally posted by 574384607 at 2010-08-10 15:20:23:
体系在水盒子里做模拟,运行过程中,
FATAL ERROR: Periodic cell has become too small for original patch grid!
Possible solutions are to restart from a recent checkpoint,
increase margin, or disabl ...

不过我已经把盒子边长设的时候每个边调大了2A了
2楼2010-08-10 15:24:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyanfei

木虫 (小有名气)

574384607(金币+1): 2011-03-26 07:35:24

遇到这个问题了,还在研究中……
以虚养心,以德养身,以仁义养天下万物,以道养天下万世
3楼2011-02-22 13:11:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)

★ ★
ghcacj(金币+1): 谢谢 2011-02-24 13:34:35
ghcacj(金币+1): 谢谢 2011-03-22 12:11:43
574384607(金币+1): 2011-03-26 07:35:28

导致这个问题的原因很多
把config贴出来吧
4楼2011-02-22 16:16:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crispsoft

金虫 (小有名气)

574384607(金币+1): 2011-03-26 07:35:37
引用回帖:
Originally posted by wally8962 at 2011-02-22 16:16:33:
导致这个问题的原因很多
把config贴出来吧

您好,我的计算过程中也出现了类似问题(FATAL ERROR: Periodic cell has become too small for original patch grid!),conf文件如下,请指教:structure          ../common/protein_ionized.psf
coordinates        ../common/protein_ionized.pdb

set temperature    310
set outputname     eq1

firsttimestep      0


# Input
paraTypeCharmm            on
parameters          ../common/par_all27_prot_lipid.inp
temperature         $temperature


# Force-Field Parameters
exclude             scaled1-4
1-4scaling          1.0
cutoff              12.0
switching           on
switchdist          10.0
pairlistdist        13.5


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

# Constant Temperature Control
langevin            on    ;#
langevinDamping     5     ;#
langevinTemp        $temperature
langevinHydrogen    off    ;#

# Periodic Boundary Conditions
cellBasisVector1    60.0    0.   0.
cellBasisVector2     0.   60.0   0.
cellBasisVector3     0.    0   50.0  
cellOrigin          15.01   5.22  -9.89

wrapAll             on

# PME (for full-system periodic electrostatics)
PME                 yes
PMEGridSpacing      1.0

#manual grid definition
PMEGridSizeX        72
PMEGridSizeY        72
PMEGridSizeZ        64  


# Constant Pressure Control (variable volume)
useGroupPressure      yes ;# needed for rigidBonds
useFlexibleCell       no
useConstantArea       no

langevinPiston        on
langevinPistonTarget  1.01325 ;#  in bar -> 1 atm
langevinPistonPeriod  200.0
langevinPistonDecay   100.0
langevinPistonTemp    $temperature

#ionized.pdb是被我固定了CA 碳的文件
if {1} {
constraints        on
consref            ../common/ionized.pdb
conskfile          ../common/ionized.pdb
conskcol           B
constraintscaling  10.0
}

# Output
outputName          $outputname

restartfreq         1000     
dcdfreq             1000
xstFreq             1000
outputEnergies      1000
outputPressure      1000

# Minimization
minimize            1000
reinitvels          $temperature

run 10000000
5楼2011-03-22 09:37:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)

★ ★
ghcacj(金币+2): 谢谢 2011-03-22 12:12:10
574384607(金币+1): 2011-03-26 07:35:42
langevinPistonPeriod  200.0
langevinPistonDecay   100.0
先把这个改大,这个值太小了,至少400、200,再看看情况
6楼2011-03-22 10:12:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)

★ ★ ★
ghcacj(金币+3): 谢谢 2011-03-22 12:12:18
574384607(金币+1): 2011-03-26 07:35:47
这个算是约束,不算固定(fix)
正确的步骤需要若干步动力学平衡,逐步减小约束K值,直到0时进入正式模拟。这是原则,目标是尽量使结构以你期望的构型开始计算。运用之妙,存乎一心。
7楼2011-03-22 11:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

476992939

铜虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖
楼主这个问题是怎么解决的?
9楼2011-10-22 15:04:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
zh1987hs(金币+2): 谢谢 2011-11-18 08:23:49
increase margin, or disable useFlexibleCell
10楼2011-11-18 08:20:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 574384607 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 70/3500 2026-08-16 08:29 by 医学老男孩
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见