分子表面静电势计算级别选用问题

葱葱cong

各位老师好，小白想请教一下，体系共46个原子（前四周期），我先前计算热力学数据用#P B3LYP/6-31G* em=gd3bj opt freq scrf(SMD,solvent=water) nosymm关键词进行几何优化振动分析后，随后使用#P M062X/Def2TZVPP em=gd3 guess=read scrf(SMD,solvent=water) SP nosymm geom=allcheck关键词得到高级别单点能。后续想用Multiwfn计算该分子的静电势，请问各位老师：
（1）能直接用高级别单点得到的波函数直接计算吗？
（2）先前几何优化的B3LYP/6-31G*级别有点低，会不会对静电势结果有影响？
（3）sob老师的博文简谈量子化学计算中DFT泛函的选择 - 思想家公社的门口：量子化学·分子模拟·二次元以及谈谈量子化学中基组的选择 - 思想家公社的门口：量子化学·分子模拟·二次元提到计算静电势用wB97M-V普通泛函、基组用6-311G**或def-TZVP就足够，有必要更换wB97M-V/def-TZVP重新优化吗？

Stardust0831

结构优化可以用B3LYP/6-31G** em=gd3bj。单点其实没必要def2-TZVPP，这个基组有很多高角动量基函数，是给双杂化泛函、post-HF、多参考计算这类的情况用的，正常M06-2X的话，用def2-TZVP就已经足够，想要更高的精度也可以直接def2-QZVP。wB97M-V据我所知、高斯还未内置；如果你想用的话，需要换成orca来算。ESP的话，B3LYP-D3(BJ)/def2-TZVP其实也足够了。
可以试试这个策略：高斯B3LYP-D3(BJ)/6-31G**结构优化，orca里wB97M-V/def2-TZVP（追求更好结果时用def2-QZVP）算单点，然后用orca算出的波函数来算ESP。其实可以全用orca，但是高斯的结构优化更为稳健，只考虑计算性价比的话可以用高斯来结构优化；orca里可以用RI加速，所以计算高精度单点的效率很好、同时有很多高斯不支持的新的高级别理论方法，所以单点可以上orca。如果需要自由能的话就正常用orca的单点和高斯的freq，在shermo里考虑频率矫正因子算自由能就行。

葱葱cong

好的感谢老师指点

sobereva

1 能
2 如果之后用更好的级别得到算静电势的波函数，这个级别优化就够了
3 wB97M-V不推荐用来优化，没额外好处，Gaussian里也不支持，ORCA里也没解析Hessian

葱葱cong

sobereva 发表于 2025-3-17 12:13
1 能
2 如果之后用更好的级别得到算静电势的波函数，这个级别优化就够了
3 wB97M-V不推荐用来优化，没额 ...

感谢sob老师解答疑惑