24小时热门版块排行榜    

查看: 1772  |  回复: 6

mjjwhs

新虫 (小有名气)

[求助] gamess 能量分解,读取gaussian初始猜测

尝试了很多次,都失败了。
我的输入文件如下,请大家看看哪里出现了错误,谢谢!
  1  $contrl scftyp=rhf runtyp=eda  icharg= 0 mult=1 $end
2  $system mwords=300 $end
  3  $scf CONV=1.0d-08 DIIS=.true. $end
  4  $basis gbasis=sto ngauss=3 $end
  5  $guess guess=moread  norb=21  $end
  6  $morokm iatm(1)=4  $end
  7  $data
。。。
8  $end
9  $vec
10 $end
11 $vec1
12 $end
13 $vec2
14 $end
错误信息如下:
515  INITIAL GUESS ORBITALS GENERATED BY MOREAD   ROUTINE.
516  *** ERROR IN READMO:   PROBLEM READING ORBITALS!
517  POSSIBLY A DAMAGED OR MANGLED ORBITAL INPUT GROUP?
518  ERROR OCCURED AT ORBITAL=     2 (MODULUS 100=   2),
519           ITS LINE NUMBER=     1
520  DATA READ FROM INPUT WAS ORBITAL=     1 LINE=     3
521  EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Sun Mar  4 20:47:57 2012

这里norb 是输入supermolecule 的轨道数目吗?
需要给出supermolecule 的初始猜测吗?
需要读入虚轨道的系数吗?
再次感谢!
回复此楼

» 猜你喜欢

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

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

mjjwhs

新虫 (小有名气)

没有人回帖啊,自己顶起啊,
真心求教各位啦
2楼2012-03-05 10:46:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

0000372

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
9  $vec
。。。。
10 $end
中间是不是应该有点什么?总感觉少点东西,可不知道少什么
3楼2012-03-05 21:42:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

不喜欢豆芽

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
mjjwhs(金币+3): 有帮助 2012-03-05 22:54:57
ben_ladeng(金币+2): 感谢应助 2012-03-06 00:11:04
不知道能不能解决问题,不过你这两个地方貌似需要改一下:
runtyp=morokuma
6  $morokm iatm(1)=4  $end这里应该加上单体的电荷。
4楼2012-03-05 21:49:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

不喜欢豆芽

木虫 (小有名气)

【答案】应助回帖


mjjwhs(金币+1): 有帮助 2012-03-06 16:17:07
还有可能是你的坐标输入有问题,你看一下你的坐标前面有没有远在的核电核数,没有这个会出错的。
5楼2012-03-05 21:51:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zklkhp

铜虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by 0000372 at 2012-03-05 21:42:20:
9  $vec
。。。。
10 $end
中间是不是应该有点什么?总感觉少点东西,可不知道少什么

$VEC group         (optional, relevant for all SCFTYP's)
                     (required if GUESS=MOREAD)

      This group consists of formatted vectors, as written
onto file PUNCH in a previous run.  It is considered good
form to retain the titling comment cards punched before
the $VEC card, as a reminder to yourself of the origin of
the orbitals.

      For Morokuma decompositions, the names of this group
are $VEC1, $VEC2, ... for each monomer, computed in the
identical orientation as the supermolecule.  For transition
moment or spin-orbit coupling runs, orbitals for states
one and possibly two are $VEC1 and $VEC2.

可惜我的问题看了这个还是没明白。。
6楼2012-05-01 21:42:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

$VEC中的数据是HF/DFT/MCSCF计算产生的波函(轨道因子),保存在PUNCH文件里,需要另外计算。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
7楼2012-05-02 00:30:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mjjwhs 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 7/350 2026-03-17 20:20 by 花125533
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 化学工程321分求调剂 +11 大米饭! 2026-03-15 14/700 2026-03-17 17:11 by ruiyingmiao
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 285化工学硕求调剂(081700) +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 材料与化工304求B区调剂 +7 邱gl 2026-03-11 8/400 2026-03-17 09:36 by 努力学习赚彩礼
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
请填处理意见