计算化学公社

标题: 请教一下关于高斯基组的选择 [打印本页]

作者
Author: 知黑守白 时间: 2022-1-3 11:37
标题: 请教一下关于高斯基组的选择
如题，在优化结构和计算频率的时候使用级别为 B3LYP/(6-31G** & LANL2DZ)

现在要计算单点能，想使用更高级别的混合基组，所以选择了 B3LYP/(6-311G** & def2-TZVPP)

想请教各位大佬，优化用 LANL2DZ 基组单点用 def2-TZVPP这样使用是否合理？

若不合理能否给出合理建议？

另外，计算的是生物蛋白-小分子簇模型体系，含亚铁血红素，有文献中加了弥散，是否有必要？

作者
Author: fffff 时间: 2022-1-3 12:49
新手回答，如有错误麻烦老师们指出，谢谢~
1.是否加DFT-D3色散校正，大概率会有弱相互作用，不加色散校正结果不可靠，不知道该不该加索性都加，加上不会带来坏处。
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）
2.B3LYP算含过渡金属体系不是最好的选择，但是很大的体系就一个Fe的话用B3LYP问题也不是很大。
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
3.Lanl2DZ精度低，可以换Lanl2TZ，
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）
3.计算单点能确实需要更高的计算级别，但是B3LYP不能用来算能量，一般算普通有机体系的能量使用M06-2X，但是你的体系有过渡金属，可以换TPSS0-D4或者PBE0-D4(BJ)
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
4.计算弱相互作用能的话可以加弥散函数。
谈谈弥散函数和“月份”基组
http://sobereva.com/119
总结（个人建议）：
几何优化+频率计算：用B3LYP-D3(BJ)/(6-31G*&Lanl2TZ)或者PBE0-D3(BJ)/(6-31G*&Lanl2TZ)，Fe（第四周期）还不是必须要考虑相对论效应，也可以在ORCA里用RI-PBE0-D3(BJ)/def2-TZVP(-f)
算能量：PBE0-D4/def2-TZVP
算弱相互作用等具体任务类型参考sob老师关于基组和泛函选择的博文。

作者
Author: sobereva 时间: 2022-1-3 21:55
但凡让别人推荐计算级别，一定要把当前情况交代清楚，仔细看
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

假设你算的是亚铁血红素一类的情况，优化建议对过渡金属用SDD或lanl2TZ，此文已经明确说了lanl2DZ已经很过时了
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

至于算单点能，你不说清楚算单点能的目的是什么也明显不可能告诉你该用什么。诸如涉及到计算配位键键能，给配体部分用6-311G**明显太low

带过渡金属的情况强烈不建议用B3LYP，仔细看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
用TPSSh-D3(BJ)、PBE0-D3(BJ)都明显比B3LYP更合适。若是计算有机金属反应能量问题用wB97M-V更好。

没特殊情况没必要带弥散，仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

作者
Author: 知黑守白 时间: 2022-1-4 21:29

sobereva 发表于 2022-1-3 21:55
但凡让别人推荐计算级别，一定要把当前情况交代清楚，仔细看
在网上求助计算化学问题的时候必须把问题描述 ...

好的卢老师！好好研究了一下您发的教程

现在是这样的：

优化的计算级别：B3LYP/(6-31G** & LANL2DZ) （优化已经做完了）

打算进行单点计算的级别：PBE0/def2TZVPP

从计算级别而言，这样选取是否可行呢？

作者
Author: sobereva 时间: 2022-1-5 00:55

知黑守白发表于 2022-1-4 21:29
好的卢老师！好好研究了一下您发的教程

现在是这样的：

单点计算加上D3或D4校正
这种切换泛函的情况很莫名其妙
仔细看
量子化学研究中切换泛函应当注意的问题
http://sobereva.com/415（http://bbs.keinsci.com/thread-9878-1-1.html）

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