查看: 250  |  回复: 5
当前主题已经存档。

hjsjs

铜虫 (正式写手)

[交流] 【求助】gaussian ONIOM 求助

Hi everyone,

I am trying to run an ONIOM calculation with B3LYP/6-31+G(d,p) as a high layer
and Amber as my low layer.  The error at the end of the file says I am missing
atomic parameters (see below for the last part of the calc file).  What are
these atomic parameters and where would I obtain them to get my calculation to
work?

Thanking You in Advance,
Sarah Whittleton
Mount Allison University

Th calcaultion dies almost immediatley and ends with:

I=  105 IAn= 15 Valence= 4.
          JB=  1 J=  106 IAn=  8 IBT= 1 Dist= 1.80D+00
          JB=  2 J=  107 IAn=  8 IBT= 2 Dist= 1.48D+00
          JB=  3 J=  171 IAn=  8 IBT= 1 Dist= 1.80D+00
          JB=  4 J=  104 IAn=  8 IBT= 2 Dist= 1.58D+00
Missing atomic parameters for atom   207 IAtTyp=    20000000
Missing atomic parameters.
Error termination via Lnk1e in /var/local/g03/l101.exe at Fri Nov  7 15:18:04
2003.
Job cpu time:  0 days  0 hours  0 minutes  0.6 seconds.
File lengths (MBytes):  RWF=      7 Int=      0 D2E=      0 Chk=      1
Scr=      1

我的问题与这个人的相同,这个人在某个论坛上发的求助贴,但是,没有人回应,我又碰到相同的问题了  哎, gjf文件是gaussview产生的,应该不会有问题吧,不知道问题出在哪里?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

新虫 (著名写手)


hjsjs(金币+1,VIP+0): 9-10 17:40
难道是分层没分好?猜的
2楼2009-09-10 09:14:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)

★ ★ ★ ★
yjcmwgk(金币+2,VIP+0):呵呵,我同意的的说法——你猜的对! 9-10 12:27
hjsjs(金币+1,VIP+0): 9-10 17:40
yjcmwgk(金币+1,VIP+0): 9-10 19:02
MM方法一般没有通用性(指元素原子),不同的MM方法是针对特定体系制定的,也许你的计算体系包含了 Amber 不支持的元素原子(比如金属),乱猜的。
3楼2009-09-10 09:30:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hjsjs

铜虫 (正式写手)

★ ★ ★ ★
yjcmwgk(金币+4,VIP+0):感谢您的辛劳 9-10 19:02
我查了不支持的原子是氢原子,就是DNA上的磷酸的-OH的氢原子,我又查了MM的原始论文,会不会是atom type有问题,比如氢原子的atom type有,H,HW,HO,HS,HA等等
4楼2009-09-10 17:49:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)

★ ★ ★ ★ ★ ★
loovfnd(金币+2,VIP+0):感谢答案! 9-10 20:40
hjsjs(金币+4,VIP+0): 9-11 23:26
Gview产生的也会有问题的。
你好好看看你这个原子,根据AMBER手册,把它的参数手动写到.gjf文件里边去就得了。

记住AMBER力场的参数在Gaussian里边是不全的。通常要自己加。
增加的方法有三种,以自己给的参数优先、默认参数优先、只读gjf文件里的参数
。方法在Gaussian手册里。
5楼2009-09-10 18:57:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ferlich

金虫 (小有名气)

★ ★ ★ ★
yjcmwgk(金币+1,VIP+0): 9-10 21:49
hjsjs(金币+3,VIP+0): 9-11 23:26
引用回帖:
Originally posted by hjsjs at 2009-9-10 17:49:
我查了不支持的原子是氢原子,就是DNA上的磷酸的-OH的氢原子,我又查了MM的原始论文,会不会是atom type有问题,比如氢原子的atom type有,H,HW,HO,HS,HA等等

从amber自带的gaff.dat或parm*.dat中找找也许会有相应的原子参数
6楼2009-09-10 21:02:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hjsjs 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 985高校硕士每月会发五六千工资!!!!? +93 freies 2018-12-14 116/5800 2018-12-16 01:54 by jie740i
[找工作] 面板厂辞职后就等于失业了吗? +8 W;something 2018-12-15 8/400 2018-12-15 23:57 by precious2012
[硕博家园] +3 1576197172 2018-12-15 3/150 2018-12-15 22:14 by zs110631
[硕博家园] 导师希望把论文一作让给其他老师,是否赞成 +74 kobeone 2018-12-14 96/4800 2018-12-15 20:26 by herokiller
[硕博家园] 求一组数据之间是否有显著差异 +9 爱睡懒觉的宝 2018-12-11 12/600 2018-12-15 10:53 by 张凡云
[公派出国] CSC给的1200欧元在德国生活够吗? +25 JasmineJLU 2018-12-10 25/1250 2018-12-14 17:23 by 了了工作室
[考博] 求过来人指点 +3 考博党2 2018-12-13 6/300 2018-12-14 15:48 by 考博党2
[硕博家园] 搞什么研究, +8 泡泡飞走 2018-12-13 15/750 2018-12-14 14:02 by HYJ黄艳杰
[高分子] 丙烯酸增稠剂爆聚 5+3 随风zz 2018-12-10 5/250 2018-12-14 10:47 by li1023591626
[论文投稿] 2018年中科院期刊分区手机公众号啥时候能更新查询 +7 arui6123479 2018-12-11 10/500 2018-12-14 09:04 by wslwydgv
[论文投稿] 文章被拒了 +15 26℃的白开水 2018-12-13 37/1850 2018-12-13 22:35 by 逝水无痕1990
[硕博家园] 氧化石墨烯 +8 xwwwwww 2018-12-13 8/400 2018-12-13 21:45 by 面朝大海6
[电化学] 参比电极电位不准,请问各位大神 我该怎么办 +4 3153849727 2018-12-13 8/400 2018-12-13 16:17 by huixiong0627
[找工作] 凝聚态物理就业好迷茫 +7 xiaomogu666 2018-12-09 7/350 2018-12-13 15:37 by QQ550303332
[论文投稿] 文章发表的时候不算sci的杂志,这两年变成sci了,能不能把文章算成sci杂志发表里面 +7 wangpeijun 2018-12-09 8/400 2018-12-13 15:32 by wslwydgv
[职场人生] 一般高校教辅岗位,有什么好处? +6 jeff800_800 2018-12-10 10/500 2018-12-13 14:50 by bluice6
[硕博家园] 实验啊 +8 用户胡云云 2018-12-13 14/700 2018-12-13 13:57 by 用户胡云云
[找工作] 材料工程硕士毕业求职 +5 lucifer011 2018-12-10 7/350 2018-12-13 10:26 by 梦想成真2014
[论文投稿] 2018中科院分区各学科一二区期刊数量占比,综合期刊比例特低。。。 +8 supernova1 2018-12-11 11/550 2018-12-12 20:41 by 15201683221
[论文投稿] 期刊升区了 +6 ZXTbryant 2018-12-12 10/500 2018-12-12 11:46 by ZXTbryant
信息提示
请填处理意见