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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » CAS过渡态优化的总能量为零？
242/3返回列表上一页123下一页
查看: 1786  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

pwzhou

铁杆木虫 (正式写手)
引用回帖:
11楼: Originally posted by sujialing at 2013-11-15 17:31:51
我再尝试，可是还没摸索出来，不大会用自然轨道的，你之前用过不，方便给我个小例子吗？...

我有firefly的例子，和gamess应该稍微有点差别，等上班了我再贴出来。

[ 发自小木虫客户端 ]
赞一下 回复此楼
13楼2013-11-15 19:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)
引用回帖:
12楼: Originally posted by sujialing at 2013-11-15 17:33:03
我用CAS的优化好的轨道，在输入文件中写入其HESSIAN阵，应该可以吧，我现在是这样尝试的，不过还没出结果呢。...

估计不行。

[ 发自小木虫客户端 ]
回复此楼
14楼2013-11-15 19:29:28
已阅   回复此楼   关注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的回帖

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

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

木虫 (小有名气)
我很不想说的一个事实就是我采用输入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 的主题更新
242/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭