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

sobereva

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