应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » CAS过渡态优化的总能量为零？
52/1返回列表
查看: 1772  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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还不是太熟悉，专家们给支个招吧！
回复此楼

» 猜你喜欢
越努力，越幸运！
1楼 2013-11-13 15:05:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
★ ★ ★ ★
gkf高: 金币+4 2013-11-15 23:08:20
引用回帖:
12楼: Originally posted by sujialing at 2013-11-15 17:33:03
我用CAS的优化好的轨道，在输入文件中写入其HESSIAN阵，应该可以吧，我现在是这样尝试的，不过还没出结果呢。...

过渡态优化不要求精确的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，但是手册里没有提到这种用法，可能无效。建议先用小基组对一个小分子做测试。
赞一下 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
17楼2013-11-15 21:08:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

pwzhou

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
2楼: Originally posted by pwzhou at 2013-11-14 13:44:55
你这是CASSCF没收敛，所以Gamess最后给你报一个0的能量，实际上在之前的迭代中是有计算出的没有收敛的能量的。 CASSCF如果300步还不收敛，估计是你的轨道没选对。

下面这个是你除坐标之外完整的输入文件吗？如果 ...

上面这个不是全部的输入文件，你说的那些我都做了，只是没列出来。我的轨道什么的没有选错，就是我要找的轨道，现在存在的问题是，我用
6-31g*优化过渡态是有能量显示的，但是我用6-31+g*优化就出现上面的情况（我知道基组之间不互认，）针对上面它不收敛的问题，有没有什么好的建议？
赞一下 回复此楼
越努力，越幸运！
3楼2013-11-14 14:32:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)
引用回帖:
3楼: Originally posted by sujialing at 2013-11-14 14:32:44
上面这个不是全部的输入文件，你说的那些我都做了，只是没列出来。我的轨道什么的没有选错，就是我要找的轨道，现在存在的问题是，我用
6-31g*优化过渡态是有能量显示的，但是我用6-31+g*优化就出现上面的情况（我 ...

没什么好的建议，CASSCF不收敛的唯一建议就是回去看一下你的轨道。可以详细说一下你的初始轨道猜测是怎么计算得到的吗？
赞一下 回复此楼
4楼2013-11-14 14:41:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭