24小时热门版块排行榜    

查看: 1977  |  回复: 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
本帖内容被屏蔽

6楼2009-09-10 21:02:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hjsjs 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 明天E口面上会评 +4 布布和一二 2026-07-12 5/250 2026-07-13 10:44 by zzybjr
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:35 by emlgkvoijw
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:25 by emlgkvoijw
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:15 by emlgkvoijw
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 3/150 2026-07-13 09:56 by emlgkvoijw
[硕博家园] 大龄残疾硕士的一点执念 +17 3776737374 2026-07-06 27/1350 2026-07-13 06:05 by 鲁弥abc
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 00:24 by v2i90kmrcu
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 00:15 by v2i90kmrcu
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +10 sadpencil 2026-07-12 15/750 2026-07-12 21:50 by sadpencil
[基金申请] 祈祷青基必中 50+6 hadengry 2026-07-09 16/800 2026-07-12 20:31 by 登山虫虫
[基金申请] 生命口会评 +11 wuke100666 2026-07-07 12/600 2026-07-11 19:30 by hhs666
[考博] 27届辽宁大学应届毕业生申博 +3 可乐微风大海 2026-07-09 3/150 2026-07-10 21:23 by 贪玩的研究僧
[基金申请] 关于如何从代码看上不上会 +17 且听虎啸 2026-07-09 23/1150 2026-07-10 19:51 by sdfapple719
[有机交流] chemdraw 20+5 ?慕鱼? 2026-07-08 6/300 2026-07-10 15:37 by 紫枫星域
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[基金申请] 信息学部会评时间 +3 xuel2011 2026-07-06 3/150 2026-07-09 19:05 by mollyzhang_2003
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[基金申请] 化学口会评?应用基础研究和基础研究,希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联,生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
信息提示
请填处理意见