G09和G16使用相同参数计算HOMO/LUMO的结果差异

PLwang

各位老师好！

       我最近在尝试复现一篇文章（Impact of the valence band energy alignment at the hole-collecting interface on the photostability of wide band-gap perovskite solar cells - ScienceDirect）中图1D的有机分子HOMO计算结果。该文章在SI中明确，该计算使用的是B3LYP functional at the 6-311G(d,p) level of DFT，G09。我在重复该结果时，使用了相同的函数和基组，但使用的是G16而非G09。使用G16计算出来的结果与文献中给出的结果相差巨大（-4.1296eV,文献中计算的结果为-5.48eV）

      我在社区里查到了一些老师的建议（最好使用同一软件，此处最好使用G09复现）。我目前在尝试安装G09并重新运行该计算程序。但还是非常困惑，不同计算程序可能导致如此巨大的差异吗？

乐平

Gaussian16 的默认收敛标准比 Gaussian09 的高。
比如 Gaussian16 默认 SCF=tight, 默认 INT=ultrafine。所以，你看似是相同的泛函和基组，但是结果不同，这很正常。

sobereva

那叫泛函不叫函数

仔细看此文
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678（http://bbs.keinsci.com/thread-38911-1-1.html）
谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573（http://bbs.keinsci.com/thread-20122-1-1.html）

G09和G16的差异达不到这份儿上，从其它方面找原因。

orangesea

可尝试使用G09Defaults关键词

PLwang

非常非常感谢各位的帮助。目前该问题已解决。sobereva老师说的是对的，造成我和文献中计算巨大差异的原因是建模上存在问题。我建模时，没有从头建模，而是将此处原本的Br原子直接替换成了C原子。替换后没有及时加氢，从而导致体系结构优化结果存在较大的问题，另外也感谢orangesea老师的提醒，这个指令很有帮助。