计算化学公社

 找回密码 Forget password
 注册 Register
Views: 5169|回复 Reply: 0

[其它] 将Gaussian的ONIOM产生的fch文件给Multiwfn做波函数分析要注意的问题

[复制链接 Copy URL]

4万

帖子

99

威望

4万

eV
积分
89975

管理员

公社社长+计算化学玩家

发表于 Post on 2019-11-18 23:53:06 | 显示全部楼层 Show all |阅读模式 Reading model
2022-Jan-13注:下文有些说法不对,见http://bbs.keinsci.com/thread-27408-1-1.html里面我发的截图

老有些人用Gaussian做ONIOM任务产生的chk/fch文件直接给Multiwfn作为波函数分析用的输入文件,如果没搞清楚情况,结果可能是完全错的。这里专门澄清一下情况:

Gaussian做ONIOM(H:L)形式的两层ONIOM任务时,总共做三个计算,最后一个计算是L方法算整个体系,因此最终得到的chk文件里存的信息是"L方法算的整个体系”。如果你用的L层是分子力场,那显然chk里根本没有波函数信息,Multiwfn根本就没法读入。如果你用的L层是半经验方法,由于Multiwfn不支持半经验的波函数,因此对电子结构的分析结果完全没法用,十分不靠谱,但那些只依赖于结构的分析倒是仍可以正常使用。如果你用的L层是DFT,那么chk里的波函数就是DFT结合L层的基组对整个体系算出来的分子轨道,此时Multiwfn可以进行分析,但显然对于高层部分原子的分析只是在L层对应的比较便宜的基组下实现的。

如果你就是想对高层区域在你定义的High级别下做波函数分析,那么应当把高层区域抠出来(如果边界有共价键被打断的话则恰当地加氢饱和),然后算个单点,之后给Multiwfn做波函数分析。如果你用的低层是分子力场,而且想让这些MM区域原子的电荷对高层的影响体现出来,相当于让电子嵌入也体现在波函数分析中,那么单独算高层时应当用charge关键词设置背景电荷,电荷数值就是低层原子的电荷(用CHELPG、MK等对静电势重现性较好的原子电荷为宜),电荷的位置就是低层原子所在的位置。

用ONIOM时,out=wfn功能是没法用的(至少对于G09、G16 A.03都是如此),会导致Gaussian最后在导出wfn文件时出错。实际上一般情况也没必要用wfn文件作为Multiwfn的输入文件,通常用fch就行了。

评分 Rate

参与人数
Participants 3
eV +15 收起 理由
Reason
cnu-hjp + 5
段时间 + 5 谢谢分享
暖空 + 5 谢谢老师,很有帮助!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2023-2-7 03:57 , Processed in 0.191397 second(s), 22 queries .

快速回复 返回顶部 返回列表 Return to list