计算化学公社

标题: 小分子表面静电势全部为正求助 [打印本页]

作者
Author: VincenzoChen 时间: 2021-9-17 19:57
标题: 小分子表面静电势全部为正求助
各位老师好，我想计算一个小分子的表面静电势，使用! PWPB95 D3 def2-QZVPP def2/J def2-QZVPP/C RIJCOSX tightSCF noautostart miniprint nopop进行了单点计算之后按照教程用Multiwfn计算了小分子的表面静电势，用的都是默认的基本参数，最后得到的结果发现小分子表面的静电势全部为正，这种情况正常吗？

作者
Author: wzkchem5 时间: 2021-9-17 19:58
这个小分子带净电荷吗？如果带，正常，否则不正常

作者
Author: VincenzoChen 时间: 2021-9-17 19:59

wzkchem5 发表于 2021-9-17 19:58
这个小分子带净电荷吗？如果带，正常，否则不正常

小分子带一个正点

作者
Author: VincenzoChen 时间: 2021-9-17 20:00

wzkchem5 发表于 2021-9-17 19:58
这个小分子带净电荷吗？如果带，正常，否则不正常

好的，谢谢老师

作者
Author: snljty 时间: 2021-9-17 20:19
本帖最后由 snljty 于 2021-9-17 20:28 编辑

另外如果你直接用的这个单点的gbw文件，而输入文件没写%mp2 Density relaxed end，那得到的只是SCF级别密度而不是双杂化密度，精度可能还不如普通泛函。我相信用这个密度做静电势计算不是你的本意。不过orca 至少4.x版本还不支持这个泛函的解析导数，应该还算不了后HF密度。一般来说也没有必要用双杂化泛函的计算画静电势。

作者
Author: sobereva 时间: 2021-9-18 06:44
你当关键词下用的是PWPB95双杂化泛函的中间步骤的杂化泛函的波函数做的分析，算静电势的质量还不如常用的比如B3LYP杂化泛函。

根据我的全面测试，普通泛函里wB97M-V的静电势质量是最好的之一。如果要更好点（但改进已经不明显了），可以用DSD-PBEP86-D3(BJ)波函数对应的静电势，但记得必须让程序产生双杂化级别的波函数，需要明确要求把双杂化的自然轨道存到molden文件里

作者
Author: fhh2626 时间: 2021-9-18 09:52

VincenzoChen 发表于 2021-9-17 19:59
小分子带一个正点

小分子带正电表面静电势当然是正的...

欢迎光临计算化学公社 (http://bbs.keinsci.com/)