福井函数计算求助

nakamori

我最近在计算福井函数过程中，发现存在两种方法，1、单点能计算输出.wfn文件计算福井函数。2、使用.fchk文件得到.cube文件计算福井函数。请问这两种方法结果是否有差异？此外，对于过渡金属配合物二茂铁，使用PBE0/def2svp结构优化后，使用PBE0/def2tzvp的单点能计算fchk文件计算福井函数是否可行？感谢各位老师的解答

邱诺狼

可参考社长写的blog操作：
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484

个人意见（假设你用的是Gaussian）
1.这两种方法理论上不应该有差异，有也是数值误差。这些操作在Multiwfn上就能自动处理。
2.泛函无问题，但建议加上DFT-D3(BJ)色散校正。还有不是做完优化就一了百了了，还要检验是否有虚频、波函数是否稳定。
波函数分析用Def2-TZVP基组做单点是浪费了，可以改成和优化时一样的Def2-SVP，详见这个blog：
谈谈量子化学中基组的选择
http://sobereva.com/336

sobereva

484博文的做法是计算福井函数最简单的做法，没有之一。自己按照Multiwfn手册4.5.4节手动产生wfn文件并用Multiwfn绘制密度差的功能得到福井函数亦可，但步骤更麻烦，而且还没法同时顺带得到各种概念密度泛函框架下定义的量。

分子内色散作用对二茂铁的构象不会产生影响，因此优化不用加D3校正，但如果以后要算柔性配体，色散作用可能明显影响构象，则加D3更适合，参考《谈谈“计算时是否需要加DFT-D3色散校正？”》（http://sobereva.com/413）。但这种情况大多还需要同时用molclus（http://www.keinsci.com/research/molclus.html）做构象搜索，以确保计算福井函数用的是最稳定结构。
福井函数的计算对基组并不敏感，若以后要算较大体系，用到6-311G*或def-TZVP那种较便宜的3-zeta基组就够了。def2-SVP也凑合。

nakamori

sobereva 发表于 2025-11-20 00:19
484博文的做法是计算福井函数最简单的做法，没有之一。自己按照Multiwfn手册4.5.4节手动产生wfn文件并用Mul ...

感谢sob老师的解答。对于您在http://sobereva.com/402中提到的使用VMD对福井函数进行可视化，如果想要可视化亲核进攻，是否在计算cube文件和使用VMD的过程中用N+1替换原文的N，N替换原文的N-1即

sobereva

nakamori 发表于 2025-11-23 14:07
感谢sob老师的解答。对于您在http://sobereva.com/402中提到的使用VMD对福井函数进行可视化，如果想要可 ...

用我前面484博文说的通过Multiwfn直接产生f+格点数据，要么Multiwfn里直接显示等值面，要么导出cube文件放VMD里显示
另外，下文的做法很有用，准确性往往比f+更好
使用Multiwfn通过局部电子附着能(LEAE)考察亲核反应位点、难易及弱相互作用
http://sobereva.com/676（http://bbs.keinsci.com/thread-38416-1-1.html）