第10届北京科音分子动力学与GROMACS培训班将于11月11~14日于北京举行,是系统、正确、迅速上手分子动力学计算、掌握GROMACS使用的最好机会!报名已经开始,详情见此链接,欢迎参加

计算化学公社

 找回密码
 现在注册!
查看: 2617|回复: 0

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

[复制链接]

3万

帖子

99

威望

3万

eV
积分
73712

管理员

公社社长+计算化学玩家

发表于 2019-11-18 23:53:06 | 显示全部楼层 |阅读模式
老有些人用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就行了。

评分

参与人数 2eV +10 收起 理由
段时间 + 5 谢谢分享
暖空 + 5 谢谢老师,很有帮助!

查看全部评分

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395(已满),2号:466017436(已满),3号:764390338(可加)。合计8000人,讨论范畴相同
思想家公社的门口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!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|小黑屋|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )

GMT+8, 2021-10-21 20:19 , Processed in 0.210327 second(s), 25 queries .

快速回复 返回顶部 返回列表