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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » CAS过渡态优化的总能量为零？
242/3返回列表上一页123下一页
查看: 2143  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 娱乐 +3 Tide man 2026-08-03 3/150 2026-08-04 09:58 by gltch
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 jRl3mE6ddZGq 2026-08-03 10/500 2026-08-04 09:26 by v1TJC4qJauAB
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化，科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭