| 查看: 1811 | 回复: 2 | ||
1466036890铜虫 (正式写手)
|
[求助]
求助利用molpro做CASPT2计算的输入文件 已有1人参与
|
| 求采用molpro软件的CASPT2输入文件，想采用CASPT2计算无能垒反应最低能量路径上各点的单点能，之前只是采用molpro做过显相关CC计算，从来没有做过CASPT2计算，希望虫友们给予帮助。 |
回复此楼» 猜你喜欢
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
物理化学论文润色/翻译怎么收费?
已经有205人回复
-
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有27人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有45人回复
-
英国贝尔法斯特女王大学招收 2026 年 CSC/BBSRC phd 机器学习 催化剂设计方向
已经有0人回复
-
法国里昂大学招生CSC 资助博士生
已经有0人回复
-
法国里昂大学接收CSC 资助访问学者
已经有0人回复
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 586 (博士)
- 贵宾: 0.01
- 金币: 26316.1
- 红花: 348
- 帖子: 2021
- 在线: 1440.5小时
- 虫号: 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