24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1665  |  回复: 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的回帖

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的回帖
查看全部 7 个回答

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的回帖
信息提示
请填处理意见