24小时热门版块排行榜    

查看: 1462  |  回复: 4

zhilv

金虫 (初入文坛)


[交流] 【求助】molpro试算

请教如何利用molpro中的caspt2计算HXeOH的过渡态,
关键是求一个输入文件,谢谢各位的帮助!
回复此楼

» 收录本帖的淘帖专辑推荐

我的收藏

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
zhilv(金币+1): 2011-02-26 09:00:52
caspt2没有解析hessian,算数值hessian比较慢,而且不能利用对称性。可以用其它方法产生数值或解析hessian,提高效率。例如用hf,dft。如果是这些方法不适用的体系,也可以用小基组、小活性空间的caspt2。

首先在一个低阶方法优化的过渡态产生hessian。输入:
***, generate initial hessian for transition state optimization
memory,1000,m
file,2,hessian.dat

geometry={
...
}

BASIS={
...
}

{rks,b3lyp3;}

{optg,ROOT=2,maxit=1;}
---

hessian默认被存到2号文件,所以要给2号文件命名,以后使用。因为只需要该点的hessian,一步优化就可以(maxit=1)。察看输出文件,找到存放hessian的具体位置(例如,5302.2)。

在同样的初始过渡态结构,用caspt2做优化的时候,不要升级hessian文件(NOHESS)。另外,为了避免和caspt2的新2号文件冲突,上一步的存有hessian的2号文件改为3号,5302.2也改为5303.3。输入:

***, transition state
file,3,hessian.dat;
memory,2000,m

geometry={
...
}

BASIS={
...
}

{HF;
...}

{MULTI;
...
}

{rs2c;
...
}

{optg,ROOT=2;HESSIAN,HESSREC=5302.3;
METHOD,QSD,NOHESS;}
---

优化结束后,人工控制hessian的更新。这需要在新的过渡态重复第一、第二个计算。直到第二个计算中的结构不再更新。一般需要3~5次。

[ Last edited by beefly on 2011-2-25 at 23:08 ]
2楼2011-02-25 23:06:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhilv

金虫 (初入文坛)


引用回帖:
Originally posted by beefly at 2011-02-25 23:06:31:
caspt2没有解析hessian,算数值hessian比较慢,而且不能利用对称性。可以用其它方法产生数值或解析hessian,提高效率。例如用hf,dft。如果是这些方法不适用的体系,也可以用小基组、小活性空间的caspt2。

首先 ...

谢谢beefly!
3楼2011-02-26 09:04:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhilv

金虫 (初入文坛)


引用回帖:
Originally posted by beefly at 2011-02-25 23:06:31:
caspt2没有解析hessian,算数值hessian比较慢,而且不能利用对称性。可以用其它方法产生数值或解析hessian,提高效率。例如用hf,dft。如果是这些方法不适用的体系,也可以用小基组、小活性空间的caspt2。

首先 ...

beelfy,您好!
我按您的方法,去试算了下,总是无法正常结束。
我想我是对您所给的输入方法没理解透彻,所以没有试算成功。
不知道您有没有一份完整的输入文件(CASPT2算过渡态)?有的话能否
发给我一份?我的邮箱lvzhi0619@163.com
4楼2011-02-27 11:16:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

974358979

铁虫 (初入文坛)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
1415360楼: Originally posted by zhilv at 2011-02-27 11:16:01
beelfy,您好!
我按您的方法,去试算了下,总是无法正常结束。
我想我是对您所给的输入方法没理解透彻,所以没有试算成功。
不知道您有没有一份完整的输入文件(CASPT2算过渡态)?有的话能否
发给我一份?我的 ...

您好,有没有CASPT2计算方面的详细解释,我是新手,不知道该怎样用来计算小体系,
5楼2013-03-01 21:18:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhilv 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:13 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 13:02 by ECLTXnJnRSXf
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:28 by ECLTXnJnRSXf
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:00 by HajghKOMynmz
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 10:38 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:30 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:27 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 09:52 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 3/150 2026-11-12 08:26 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:30 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:30 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 05:19 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
信息提示
请填处理意见