计算化学公社
标题:
分子表面静电势计算级别选用问题
[打印本页]
作者Author:
葱葱cong
时间:
2025-3-16 12:20
标题:
分子表面静电势计算级别选用问题
各位老师好,小白想请教一下,体系共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重新优化吗?
作者Author:
Stardust0831
时间:
2025-3-16 16:36
本帖最后由 Stardust0831 于 2025-3-16 16:43 编辑
结构优化可以用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里考虑频率矫正因子算自由能就行。
作者Author:
葱葱cong
时间:
2025-3-17 11:57
好的感谢老师指点
作者Author:
sobereva
时间:
2025-3-17 12:13
1 能
2 如果之后用更好的级别得到算静电势的波函数,这个级别优化就够了
3 wB97M-V不推荐用来优化,没额外好处,Gaussian里也不支持,ORCA里也没解析Hessian
作者Author:
葱葱cong
时间:
2025-3-17 13:41
sobereva 发表于 2025-3-17 12:13
1 能
2 如果之后用更好的级别得到算静电势的波函数,这个级别优化就够了
3 wB97M-V不推荐用来优化,没额 ...
感谢sob老师解答疑惑
欢迎光临 计算化学公社 (http://bbs.keinsci.com/)
Powered by Discuz! X3.3