| 查看: 1854 | 回复: 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还不是太熟悉，专家们给支个招吧！ |
回复此楼» 猜你喜欢
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有176人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
分子动力学模拟合作需求
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯能量单位转换：优化和过渡态中，所得能量的单位及单位转换
已经有7人回复
- 关于过渡态和能量优化的问题
已经有7人回复
- 过渡态优化，iop(1/8=1),能量曲线依然震荡，我还要hold住吗？
已经有17人回复
pwzhou
铁杆木虫 (正式写手)
- QC强帖: 2
- 应助: 51 (初中生)
- 金币: 7551.4
- 散金: 13
- 红花: 14
- 帖子: 412
- 在线: 316.9小时
- 虫号: 56225
- 注册: 2005-01-07
- 性别: GG
- 专业: 化学动力学
9楼2013-11-15 16:02:38
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