| 查看: 2220 | 回复: 23 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
sujialing木虫 (小有名气)
|
[求助]
CAS过渡态优化的总能量为零？
|
||||
|
我采用GAMESS的CAS进行过渡态的优化，正常结束，但是我的反应物和过渡态的TOTAL ENERGY都是0，没有取代基的时候是能算对的。我没弄明白这是怎么回事儿，请大家帮忙给解决一下吧！ 不知道有没有人用CAS做过振动分析，遇到过这样的情况没，出现下面的错误提示： TOO MANY ITERATIONS IN MCCPCG *** MAXIT= 300 EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Tue Nov 12 15:01:26 2013 我做振动分析的结构和轨道都是上一步优化好的，确实是过渡态。输入文件如下： $CONTRL SCFTYP=MCSCF RUNTYP=HESSIAN NZVAR=51 MAXIT=100 COORD=ZMT MULT=2 ISPHER=1 $END $SYSTEM TIMLIM=525600 MEMORY=500000000 $END $BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 DIFFSP=.FALSE. $END $MCSCF MAXIT=100 $END 尝试过加大MAXIT不过没有效果，还是出现上面的错误提示。 刚刚接触Gamess还不是太熟悉，专家们给支个招吧！ |
回复此楼» 猜你喜欢
中南大学化学化工学院招收电池储能、理论计算方向博士生
已经有16人回复
-
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
物理化学论文润色/翻译怎么收费?
已经有146人回复
-
关于合成MIL-125的求助
已经有0人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13336.5
- 散金: 27307
- 红花: 301
- 帖子: 4683
- 在线: 928.6小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★ ★ ★
gkf高: 金币+4 2013-11-15 23:08:20
gkf高: 金币+4 2013-11-15 23:08:20
|
过渡态优化不要求精确的HESSIAN，只要提供一个大致的虚频方向就行。像Molpro就可以用CASPT2解析梯度+HF、CASSCF级别的解析HESSIAN或DFT、MP2的数值HESSIAN。为了保险，也可以每隔几步就人工更新一次HESSIAN。在得到过渡态以后，更新HESSIAN重算一次。 GAMESS有个类似的做法：$STATPT HESS=CALC IHREP=0 $END，但是必须在第一步用MCQDPT计算一次数值HESSIAN。 用$STATPT HESS=READ IHREP=0 $END可以读入HESSIAN，但是手册里没有提到这种用法，可能无效。建议先用小基组对一个小分子做测试。 |
17楼2013-11-15 21:08:02
pwzhou
铁杆木虫 (正式写手)
- QC强帖: 2
- 应助: 51 (初中生)
- 金币: 7551.4
- 散金: 13
- 红花: 14
- 帖子: 412
- 在线: 316.9小时
- 虫号: 56225
- 注册: 2005-01-07
- 性别: GG
- 专业: 化学动力学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
sujialing: 金币+3 2013-11-14 14:26:01
gmy1990: 金币+3 2013-11-14 16:14:27
sujialing: 金币+7, ★★★很有帮助 2013-11-15 17:33:36
感谢参与，应助指数 +1
sujialing: 金币+3 2013-11-14 14:26:01
gmy1990: 金币+3 2013-11-14 16:14:27
sujialing: 金币+7, ★★★很有帮助 2013-11-15 17:33:36
|
你这是CASSCF没收敛，所以Gamess最后给你报一个0的能量，实际上在之前的迭代中是有计算出的没有收敛的能量的。 CASSCF如果300步还不收敛，估计是你的轨道没选对。 下面这个是你除坐标之外完整的输入文件吗？如果是的话，这样的输入文件恐怕是不行的，都没看到你用的活化空间是多大？另外算CASSCF必然是要用事先算好的轨道，并在必要条件下（事实上是大多数条件下）进行轨道交换的，这个在你的输入文件中都没体现。建议你好好看看CASSCF如何计算并且好好研究一下Gamess中关于CASSCF的说明。 $CONTRL SCFTYP=MCSCF RUNTYP=HESSIAN NZVAR=51 MAXIT=100 COORD=ZMT MULT=2 ISPHER=1 $END $SYSTEM TIMLIM=525600 MEMORY=500000000 $END $BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 DIFFSP=.FALSE. $END $MCSCF MAXIT=100 $END |
2楼2013-11-14 13:44:55
sujialing
木虫 (小有名气)
- 应助: 10 (幼儿园)
- 金币: 1541.1
- 散金: 41
- 红花: 1
- 帖子: 257
- 在线: 671.9小时
- 虫号: 1218664
- 注册: 2011-03-02
- 性别: MM
- 专业: 理论和计算化学
3楼2013-11-14 14:32:44
pwzhou
铁杆木虫 (正式写手)
- QC强帖: 2
- 应助: 51 (初中生)
- 金币: 7551.4
- 散金: 13
- 红花: 14
- 帖子: 412
- 在线: 316.9小时
- 虫号: 56225
- 注册: 2005-01-07
- 性别: GG
- 专业: 化学动力学
4楼2013-11-14 14:41:26
20