24小时热门版块排行榜     先丰纳米 买纳米材料就点这里

查看: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[分析] 新买的色谱柱用了几天就出现小肩峰 +13 张吕大夫 2018-04-16 20/1000 2018-04-21 09:00 by 良辰1314521
[找工作] 在国家专利局里面,博士和硕士有什么不同? +4 喵喵喵博士 2018-04-20 7/350 2018-04-21 07:22 by yumanxi052
[文学芳草园] 随笔?听雨 +3 Z&$@ 2018-04-18 33/1650 2018-04-20 22:28 by myrtle
[公派出国] 关于联合培养,国家筛选申请人主要看什么? +4 csutzk 2018-04-19 4/200 2018-04-20 19:57 by 果果陶陶
[硕博家园] 北京中科院过程工程研究所工程实验室招聘 +3 tianjitz 2018-04-20 7/350 2018-04-20 09:36 by tianjitz
[教师之家] 考高中教师资格证,专业课没过,只有科一科二过了能报小学的面试吗 +3 临溪村姑 2018-04-19 8/400 2018-04-20 08:34 by 临溪村姑
[公务员考试] 全额事业编 +3 烧饼黄又黄 2018-04-16 3/150 2018-04-20 08:07 by 通缉要饭
[考研] 兰州生物制药研究所非全值得去吗?纠结 +6 lucky,big 2018-04-19 8/400 2018-04-19 16:30 by bluehaiku
[硕博家园] 放养模式研一 该如何努力 +7 Stella- 2018-04-19 7/350 2018-04-19 15:55 by 萌萌哒小兔子
[考博] 今日考博复试,自我感觉遭透了…… +7 xuhongming28 2018-04-18 12/600 2018-04-19 10:56 by zkaiqing
[有机交流] 还原羰基,用的三氯化铝,一般用什么溶剂合适! +3 小白鼠110 2018-04-18 3/150 2018-04-19 07:29 by tlllw
[文学芳草园] 时间赠人阅历,世事尽可原谅 +7 孤鸿飘逸 2018-04-18 7/350 2018-04-19 00:14 by 我是本笨蛋
[找工作] 急求博士讲师在苏科大的待遇 +13 georari 2018-04-17 36/1800 2018-04-18 17:18 by andee
[考研] 考不考数学 +4 奇妙能力L 2018-04-18 8/400 2018-04-18 13:34 by 奇妙能力L
[硕博家园] 博士开过题了,可以换题目吗 +7 小沙弥fly 2018-04-17 15/750 2018-04-18 09:14 by erimac
[硕博家园] !看到我!第一次见博导需要准备些神马~ +3 janni_k 2018-04-16 6/300 2018-04-17 21:47 by keithgong
[硕博家园] 中国科学院大学雁栖湖校区附近租房一年 +4 零的开始 2018-04-16 10/500 2018-04-17 17:35 by 五月樱舞001
[考研] 二战 +10 染指流年~ 2018-04-16 23/1150 2018-04-17 15:36 by 染指流年~
[公派出国] 联培的value值变化的意义,大家讨论下? +7 water6459531 2018-04-16 12/600 2018-04-17 09:34 by 为什么这么
[版块工场] 【早起签到贴】2018年4月17日(q q 48626260 进群申请:小木虫——您个人ID) +64 8475 2018-04-17 66/3300 2018-04-17 08:15 by 心碎念想fly
信息提示
请填处理意见