24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhilv

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 481.5
帖子: 40
在线: 19.7小时
虫号: 644776

[交流] 【求助】molpro试算

请教如何利用molpro中的caspt2计算HXeOH的过渡态，
关键是求一个输入文件，谢谢各位的帮助！

回复此楼

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

查看全部散金贴

1楼 2011-02-25 20:52:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

974358979

铁虫 (初入文坛)

应助: 3 (幼儿园)
金币: 535.4
帖子: 22
在线: 169小时
虫号: 2075703

★
小木虫: 金币+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的回帖

查看全部 5 个回答

beefly

专家顾问 (职业作家)

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

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 481.5
帖子: 40
在线: 19.7小时
虫号: 644776

引用回帖:

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的回帖