| 查看: 2221 | 回复: 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还不是太熟悉，专家们给支个招吧！ |
回复此楼» 猜你喜欢
关于合成MIL-125的求助
已经有0人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有119人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
sujialing
木虫 (小有名气)
- 应助: 10 (幼儿园)
- 金币: 1541.1
- 散金: 41
- 红花: 1
- 帖子: 257
- 在线: 671.9小时
- 虫号: 1218664
- 注册: 2011-03-02
- 性别: MM
- 专业: 理论和计算化学
|
恩，谢谢啦。我也是这么算的，我按照您给的提示第一步先用MCQDPT优化过渡态，但是它只是读取一部分轨道然后任务就终止了，我没分析出是什么原因，输出文件中也没有什么错误提示，如果您有时间的话，帮我看一下吧，谢谢啦！ 我的输入文件$CONTRL SCFTYP=MCSCF RUNTYP=SADPOINT MPLEVL=2 NZVAR=51 MAXIT=100 COORD=ZMT MULT=2 ISPHER=1 NUMGRD=.T. $END $SYSTEM TIMLIM=525600 MEMORY=5000000000 $END $BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 $END $MCSCF MAXIT=1000 $END $STATPT HESS=CALC IHREP=0 $END $DATA 分子描述 $END $GUESS GUESS=MOREAD NORB=158 norder=0 $END $DET NCORE=26 NACT=3 NELS=3 NSTATE=1 $END $VEC $END 我的输出文件：http://pan.baidu.com/disk/home |
22楼2013-11-18 14:08:42
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