24小时热门版块排行榜    

查看: 2171  |  回复: 7
本帖产生 1 个 模拟EPI ,点击这里进行查看

plafalines

捐助贵宾 (小有名气)

[求助] 用tinker的amoeba09力场模拟~~总是提示参数未定义~~~怎么办啊~热切求助

用tinker的amoeba09模拟一个NaCl水溶液体系 总是提示:

Random Number Generator Initialized with SEED :      123456789

MECHANIC  --  Some Required Potential Energy Parameters are Undefined

TINKER is Unable to Continue; Terminating the Current Calculation

xyz文件和key文件见附件

求助各位高手

[ Last edited by plafalines on 2011-5-3 at 22:55 ]
回复此楼

» 猜你喜欢

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

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

piaoma

木虫 (小有名气)

【答案】应助回帖

★ ★
plafalines(金币+2): 2011-05-04 22:50:18
御剑江湖(金币+2): 谢谢 2011-05-05 12:11:01
据我所知,amoeba是ponder开发的针对水的力场,考虑极化的;这个力场里定义了Na和Cl吗?
2楼2011-05-04 22:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

plafalines

捐助贵宾 (小有名气)

引用回帖:
Originally posted by plafalines at 2011-05-03 13:40:17:
用tinker的amoeba09模拟一个NaCl水溶液体系 总是提示:

Random Number Generator Initialized with SEED :      123456789

MECHANIC  --  Some Required Potential Energy Parameters are Undefined

T ...

定义了啊
3楼2011-05-04 22:47:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fcenf

木虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
御剑江湖(金币+5, 模拟EPI+1): 感谢无私回答 2011-07-14 19:57:56
引用回帖:
Originally posted by plafalines at 2011-05-03 13:40:17:
用tinker的amoeba09模拟一个NaCl水溶液体系 总是提示:

Random Number Generator Initialized with SEED :      123456789

MECHANIC  --  Some Required Potential Energy Parameters are Undefined

T ...

你的key文件里只有力场参数,没有动力学计算需要的参数,当然会提示有些参数没有定义啦~
这里有我测试用的key文件。你用时把parameters这行后面的改成你自己装的Tinker中amoeba09力场参数文件的位置。其他参数你改成适合你体系的就行了。


parameters       /opt/tinker/params/amoeba09
verbose
save-velocity
archive


polar-eps        0.01

cutoff           9.0
neighbor-list
list-buffer      0.5

ewald
ewald-cutoff     9.0
pme-grid         19 19 19
pme-order        5
fft-package      fftw

a-axis           19.8647
integrate        verlet
thermostat       berendsen
tau-temperature  0.2
barostat         berendsen
tau-pressure     1.0
compress         0.000046
http://home.ustc.edu.cn/~fucf
4楼2011-07-14 00:00:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)


★ ★ ★ ★ ★
ghcacj(金币+5): 谢谢 2011-07-23 09:24:51
你这key文件定义了PME,但是4.2版手册67页上说:
EWALD SUMMATION METHODS
TINKER contains a versions of the Ewald summation technique for inclusion of long
range electrostatic interactions via periodic boundaries. The particle mesh Ewald (PME)
method is available for simple charge-charge potentials, while regular Ewald is provided for
polarizable atomic multipole interactions.
即TINKER中PME方法不用于极化力场的静电计算,它只适于固定电荷模型。
5楼2011-07-23 03:43:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fcenf

木虫 (初入文坛)

★ ★ ★ ★ ★
zh1987hs(金币+5): 鼓励讨论 2011-07-24 13:21:24
引用回帖:
Originally posted by ChemiAndy at 2011-07-23 03:43:35:
你这key文件定义了PME,但是4.2版手册67页上说:
EWALD SUMMATION METHODS
TINKER contains a versions of the Ewald summation technique for inclusion of long
range electrostatic interactions via perio ...

这个设置PME参数,计算时程序也会根据体系实际的力场形式选择计算方法的,如果对于多偶极势采用Ewald算法,那设置PME参数的参数应该会自动忽略掉~
另外我用的tinker版本是5.0的。我看了下源码,mpole的Ewald算法,其虚空间也是用PME算法处理了,不知道是不是新版本改了~可以看下源码里的 empole1.f 和 pmepol.f
http://home.ustc.edu.cn/~fucf
6楼2011-07-23 20:22:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fcenf

木虫 (初入文坛)

★ ★ ★
御剑江湖(金币+3): 谢谢 2011-08-06 09:37:42
引用回帖:
Originally posted by fcenf at 2011-07-23 20:22:55:
这个设置PME参数,计算时程序也会根据体系实际的力场形式选择计算方法的,如果对于多偶极势采用Ewald算法,那设置PME参数的参数应该会自动忽略掉~
另外我用的tinker版本是5.0的。我看了下源码,mpole的Ewald ...

上面的写错了,源码文件应该是 empole.f, induce.f 和 pmepol.f ~~
http://home.ustc.edu.cn/~fucf
7楼2011-07-23 20:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)


恩,你是对的,PME用于从所有点电荷计算当前盒子的周期势场;当有极化的时候,就用Ewald计算多极矩的那部分作用;
8楼2011-08-05 05:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 plafalines 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +5 yangfz 2026-03-17 5/250 2026-03-17 17:01 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 工科材料085601 279求调剂 +3 困于星晨 2026-03-17 3/150 2026-03-17 14:08 by ms629
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
信息提示
请填处理意见