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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » CAS过渡态优化的总能量为零？
151/2返回列表12下一页
查看: 1770  |  回复: 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的回帖

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的回帖

sujialing

木虫 (小有名气)
引用回帖:
4楼: Originally posted by pwzhou at 2013-11-14 14:41:26
没什么好的建议，CASSCF不收敛的唯一建议就是回去看一下你的轨道。可以详细说一下你的初始轨道猜测是怎么计算得到的吗？...

我用的是HF优化，得到轨道，然后从中挑出我要的轨道，然后改用CAS做单点，完成后检查活化空间，在优化
赞一下 回复此楼
越努力，越幸运！
6楼2013-11-14 15:01:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
5楼: Originally posted by pwzhou at 2013-11-14 14:43:03
通常如果轨道猜测好的话，收敛也就10几圈的事情。...

是啊，我算的没有取代基的时候，两步就收敛了，郁闷
赞一下 回复此楼
越努力，越幸运！
7楼2013-11-14 15:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
你用过MCQDPt优化过渡态吗？
赞一下 回复此楼
越努力，越幸运！
8楼2013-11-15 14:38:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
10楼: Originally posted by pwzhou at 2013-11-15 16:03:27
建议你用MP2的自然轨道试一下。...

我再尝试，可是还没摸索出来，不大会用自然轨道的，你之前用过不，方便给我个小例子吗？
赞一下 回复此楼
越努力，越幸运！
11楼2013-11-15 17:31:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
9楼: Originally posted by pwzhou at 2013-11-15 16:02:38
MCQDPT这个没有解析梯度吧，应该优化不了过渡态吧。...

我用CAS的优化好的轨道，在输入文件中写入其HESSIAN阵，应该可以吧，我现在是这样尝试的，不过还没出结果呢。
赞一下 回复此楼
越努力，越幸运！
12楼2013-11-15 17:33:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
13楼: Originally posted by pwzhou at 2013-11-15 19:29:12
我有firefly的例子，和gamess应该稍微有点差别，等上班了我再贴出来。
...

好的！谢谢啦！我在它的官网上怎么没看到具体的例子，只是有一些文件介绍。
赞一下 回复此楼
越努力，越幸运！
15楼2013-11-15 19:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
14楼: Originally posted by pwzhou at 2013-11-15 19:29:28
估计不行。
...

不要吧，那岂不是我努力这么久的东西要完全放弃了？？到现在也没出结果，也看不出能不能算出来，机器也不给力啊！唉！
赞一下 回复此楼
越努力，越幸运！
16楼2013-11-15 19:46:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
17楼: Originally posted by beefly at 2013-11-15 21:08:02
过渡态优化不要求精确的HESSIAN，只要提供一个大致的虚频方向就行。像Molpro就可以用CASPT2解析梯度+HF、CASSCF级别的解析HESSIAN或DFT、MP2的数值HESSIAN。为了保险，也可以每隔几步就人工更新一次HESSIAN。在得到 ...

好的，我试试，谢谢啦！
赞一下 回复此楼
越努力，越幸运！
18楼2013-11-16 10:41:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sujialing 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭