| 查看: 1903 | 回复: 2 | ||
1466036890铜虫 (正式写手)
|
[求助]
求助利用molpro做CASPT2计算的输入文件 已有1人参与
|
| 求采用molpro软件的CASPT2输入文件，想采用CASPT2计算无能垒反应最低能量路径上各点的单点能，之前只是采用molpro做过显相关CC计算，从来没有做过CASPT2计算，希望虫友们给予帮助。 |
回复此楼» 猜你喜欢
PbS量子点紫外吸收
已经有0人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有283人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
七嗪类物质合成求助
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26421.1
- 红花: 349
- 帖子: 2026
- 在线: 1442.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
