24小时热门版块排行榜    

查看: 1404  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhilv

金虫 (初入文坛)


[交流] 【求助】molpro试算

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

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

我的收藏

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注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的回帖
查看全部 5 个回答
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的回帖

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的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
信息提示
请填处理意见