24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1664  |  回复: 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 的主题更新
信息提示
请填处理意见