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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » CAS过渡态优化的总能量为零？
242/3返回列表上一页123下一页
查看: 2232  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 15:47 by xiangy672
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[基金申请] 国社科又开始会评了，不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:10 by l0VvVHGBGRLv
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭