24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2220  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 +16 2261744733 2026-04-08 16/800 2026-04-09 07:58 by Sammy2
[考研] 一志愿0703化学招61最终排名62化学求调剂 +13 招61排名62 2026-04-07 15/750 2026-04-08 23:06 by 招61排名62
[考研] 308求调剂 +16 倘若起风了呢 2026-04-05 16/800 2026-04-08 23:03 by 猪会飞
[考研] 专硕0854初试考材科基,求调剂 +7 3220548044 2026-04-06 10/500 2026-04-08 21:59 by hypershenger
[考研] 材料专硕(0856) 339分求调剂 +16 哈哈哈鹅哈哈哈 2026-04-05 16/800 2026-04-08 16:02 by luoyongfeng
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸! 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 生物与医药求调剂 +7 heguanhua 2026-04-05 8/400 2026-04-06 18:41 by macy2011
[考研] 085405软件工程301分求调剂,专硕可跨专业,四六级已过 +3 静静想想 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[考研] 材料工程310专硕调剂 +14 捞捞我…. 2026-04-04 15/750 2026-04-06 14:18 by lqwchd
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 生物与医药调剂 +4 十七sa 2026-04-05 4/200 2026-04-05 20:05 by lys0704
[考研] 329求调剂 +17 miaodesi 2026-04-02 20/1000 2026-04-05 18:33 by 蓝云思雨
[考研] 288求调剂 一志愿哈工大 材料与化工 +13 洛神哥哥 2026-04-03 13/650 2026-04-05 17:27 by zzx2138
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 材料调剂 +11 吴棂颖! 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
信息提示
请填处理意见