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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 【求助】molpro试算
51/1返回列表
查看: 1173  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhilv

金虫 (初入文坛)

[交流] 【求助】molpro试算

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

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

我的收藏

» 猜你喜欢

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

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

查看全部散金贴
1楼 2011-02-25 20:52:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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

金虫 (初入文坛)
引用回帖:
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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭