请选择 进入手机版 | 继续访问电脑版

计算化学公社

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

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

[复制链接]

2万

帖子

25

威望

2万

eV
积分
54442

管理员

公社社长

发表于 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就行了。

评分

参与人数 1eV +5 收起 理由
暖空 + 5 谢谢老师,很有帮助!

查看全部评分

第12届北京科音初级量子化学培训班将于六月举办!是零基础和初学者速成量子化学计算的非常难得和宝贵的机会,欢迎参加!详情和报名方式见http://www.keinsci.com/workshop/12KEQC
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计约6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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, 2020-5-31 18:25 , Processed in 0.149237 second(s), 27 queries .

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