| 查看: 2037 | 回复: 2 | ||
1466036890铜虫 (正式写手)
|
[求助]
求助利用molpro做CASPT2计算的输入文件 已有1人参与
|
| 求采用molpro软件的CASPT2输入文件，想采用CASPT2计算无能垒反应最低能量路径上各点的单点能，之前只是采用molpro做过显相关CC计算，从来没有做过CASPT2计算，希望虫友们给予帮助。 |
回复此楼» 猜你喜欢
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有102人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
分子动力学模拟合作需求
已经有4人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26695.9
- 红花: 349
- 帖子: 2029
- 在线: 1461.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
1466036890: 金币+20, ★★★很有帮助, 已收到，非常感谢 2023-05-11 10:50:11
感谢参与，应助指数 +1
1466036890: 金币+20, ★★★很有帮助, 已收到，非常感谢 2023-05-11 10:50:11
|
链接：https://pan.baidu.com/s/1yf4zrbVU8ZPnDrkKg-o-mQ?pwd=hgfk 提取码：hgfk --来自百度网盘超级会员V9的分享 可以看看这个例子 之前算过的一个示例。 |
2楼2023-05-10 21:48:47
1234zou
金虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 1941.5
- 红花: 5
- 帖子: 92
- 在线: 31.6小时
- 虫号: 2914121
- 注册: 2014-01-04
- 性别: GG
- 专业: 理论和计算化学
|
还有另一种途径是不学习Molpro各种语法和功能，直接采用MOKIT程序黑箱、自动式的多参考计算，这里以两根O-H键均被拉长的水分子为例，写一个高斯gjf文件 %mem=4GB %nprocshared=2 #p CASPT2/cc-pVDZ mokit{CASPT2_prog=Molpro} 0 1 O -0.23497692 0.90193619 -0.068688 H 1.26502308 0.90193619 -0.068688 H -0.73568721 2.31589843 -0.068688 输入文件使用最常用的gjf格式，想要什么理论方法直接写（哪怕高斯不支持该方法也行，因为并不是仅调用高斯算的），符合书写直觉。中间CASPT2_prog=Molpro表示调用Molpro做CASPT2计算，一样符合直觉。用automr提交，命令如下 automr h2o.gjf >h2o.out 2>&1 仅20秒结束，不用肉眼看轨道、挑选轨道、做轨道局域化、大小基组投影、担心收敛性、调试关键词等各种曲折操作，直接算出来就是CASSCF(4,4)和CASPT2(4,4)，程序会自动确定活性轨道是两根O-H的成键、反键轨道。h2o.out输出文件里有各个方法的电子能量、双自由基特征、未成对电子数等大量信息，直观清楚 E(CASCI) = -75.90802109 a.u. E(CASSCF) = -75.90823063 a.u. ----------------------- Radical index ----------------------- biradical character (2c^2) y0= 0.165 tetraradical character(2c^2) y1= 0.139 Yamaguchi's unpaired electrons (sum_n n(2-n) ): 1.121 Head-Gordon's unpaired electrons(sum_n min(n,(2-n))): 0.607 Head-Gordon's unpaired electrons(sum_n (n(2-n))^2 ): 0.315 ------------------------------------------------------------- Leave subroutine do_cas at Sat Aug 12 14:46:05 2023 Enter subroutine do_mrpt2... CASPT2 based on CASSCF orbitals. Frozen_core = F, CASPT2 using program molpro $molpro -t 1 -n 2 -m 250m h2o_uhf_gvb4_CASSCF_CASPT2.com E(ref) = -75.90823031 a.u. IP-EA shift = 0.25 (default) E(corr) = -0.13847648 a.u. E(CASPT2) = -76.04670678 a.u. Leave subroutine do_mrpt2 at Sat Aug 12 14:46:16 2023 另有fch文件可用GaussView和Multiwfn直接打开（不管你用什么量化程序算，都会给出fch文件，十分通用） CASSCF初始轨道 h2o_uhf_uno_asrot2gvb4_s.fch CASSCF收敛轨道，也是自然轨道 h2o_uhf_gvb4_CASSCF_NO.fch MOKIT程序免费、开源，见https://gitlab.com/jxzou/mokit 中文简介见https://gitlab.com/jxzou/mokit/-/blob/master/README_zh.md |
3楼2023-08-12 14:54:16
