24小时热门版块排行榜    

查看: 1773  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +8 冬十三 2026-03-15 9/450 2026-03-17 22:37 by 冬十三
[考研] 085601材料工程专硕求调剂 +5 慕寒mio 2026-03-16 5/250 2026-03-17 21:31 by hmn_wj
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +3 小李0854 2026-03-16 3/150 2026-03-17 13:40 by 热情沙漠
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
信息提示
请填处理意见