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

24小时热门版块排行榜    

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

sujialing

木虫 (小有名气)
我很不想说的一个事实就是我采用输入HESSIAN阵，让它自动读取然后计算的方法，在历时三天两夜的计算后，失败了
赞一下 回复此楼
越努力，越幸运！
19楼2013-11-18 09:47:50
已阅   回复此楼   关注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。在得到 ...

我再问您一问题啊，高斯的过渡态的优化，我会做，也明白一些道理，但是到GAMESS上就有点晕了，您提到的“STATPT HESS=CALC IHREP=0 $END，但是必须在第一步用MCQDPT计算一次数值HESSIAN” 是直接用MCQDPT计算RUNTYP=HESSIAN？但是这是做振动分析的时候，还是我理解错误了，您在帮我解答一下吧！谢谢啦！
第二种方法通过我的实践证明确实是无效的，算不出来。
赞一下 回复此楼
越努力，越幸运！
20楼2013-11-18 09:53:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
21楼: Originally posted by beefly at 2013-11-18 12:44:28
用STATPT HESS=CALC IHREP=0 $END优化过渡态，gamess自动计算HESSIAN...

恩，谢谢啦。我也是这么算的，我按照您给的提示第一步先用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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)
引用回帖:
23楼: Originally posted by beefly at 2013-11-18 22:58:41
没看到输出文件

昨天太晚了，就没告诉你，我应该是内存控制不好，我把内存改了改，现在在算呢，一会儿看看是不是朝着预期走的，太麻烦你了
赞一下 回复此楼
越努力，越幸运！
24楼2013-11-19 08:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sujialing 的主题更新
152/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭