24小时热门版块排行榜    

查看: 1786  |  回复: 23

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 的主题更新
信息提示
请填处理意见