改变方法基组能量竟差了10kcal！

liyuanhe211 · 发表于 Post on 2017-1-20 15:37:18

fhh2626 发表于 2017-1-20 15:13
我是做经典动力学的，一般就用量化来拟合一下经典或者可极化力场，通常都用后HF方法，弥散函数的影响还是 ...

做力场用6-31G(d)我认为更多的是因为弥散耗时且增加精度的效率很低（不实惠），而并非主要因为加了弥散效果反而变差了。
我只是没遇到也没听说过DFT里会出现这样的情况，同时觉得DFT中若出现这样的情况不好理解。

liyuanhe211 · 发表于 Post on 2017-1-20 22:12:57

实验上对应过程的温度大概多少。

agent99 · 发表于 Post on 2017-1-21 01:42:39

wbn 发表于 2017-1-20 09:35
为什么能说10个kcal的差别并不奇怪？除非其中一种方法理论上明显有缺陷，否则两种计算方法结果不能重复， ...

那你这么想：换了个泛函，差了3kcal，是不是觉得还好？换个基组差出3kcal来也很正常吧，再换个色散校正又差出几个kcal，加起来就10kcal了，正好没有cancel out而已

agent99 · 发表于 Post on 2017-1-21 01:51:16

小巨人发表于 2017-1-20 09:02
但是其他过程差的不多，也就顶多5kcal；由于这一步的差异，具体该选哪一步的催化出现了问题，所以您认为 ...

明尼苏达泛函最好用ultrafine grid，不知你是否用了

sobereva · 发表于 Post on 2017-1-21 10:50:15

M062X和B3LYP算热力学数据本来就相差很大，况且基组还不同，更是雪上加霜。
如果是有机体系，直接用M062X/6-311+G** D3的结果就完了，甭管B3LYP/6-31++G** D3BJ的结果。
如果你还心怀疑虑，在ORCA里用较好的双杂化泛函PWPB95-D3结合较大基组def2-TZVP或def2-QZVP算单点对比一下就知道了。

至于弥散，对当前体系加了有好处，绝对不会使结果恶化。氢的弥散是多余的，浪费计算量也没什么意义。

这是一个非常简单的例子，显见不同级别结果差个10kcal/mol司空见惯。

sobereva · 发表于 Post on 2017-1-21 10:56:10

fhh2626 发表于 2017-1-20 15:13
我是做经典动力学的，一般就用量化来拟合一下经典或者可极化力场，通常都用后HF方法，弥散函数的影响还是 ...

弥散函数过度而极化函数过少或没有，只会对后HF造成不良影响，对于DFT、HF这样单电子近似的方法没有这个问题。另外，LZ当前用的基组即便对于后HF还不至于弥散函数过度。

那些早年的力场用6-31G*这种基组完全是因为计算量所限，现在如果再用这种基组的势能面拟合参数就要被吐槽了。比如AMBER03，对二肽模型体系计算骨架二面角参数已经用MP2/cc-pVTZ了。6-31G*给后HF用是不可容忍的。一般来说，做后HF，要想得到不错结果，最起码得def2-TZVP。

也有一些搞力场的研究者对弥散函数乱用，基本还是造成无意义的耗时，不会对结果有什么不良影响。GLYCAM力场就是个不恰当使用基组的例子，拟合参数参考的是B3LYP/6-31++G(2d,2p)势能面，其实还不如用B3LYP/TZVP，算糖体系用3-zeta远比2-zeta加一堆弥散好得多，时间还能节约。

fhh2626 · 发表于 Post on 2017-1-21 15:50:34

sobereva 发表于 2017-1-21 10:56
弥散函数过度而极化函数过少或没有，只会对后HF造成不良影响，对于DFT、HF这样单电子近似的方法没有这 ...

受教了。感谢sob的回答

小巨人 · 发表于 Post on 2017-1-21 22:49:18

liyuanhe211 发表于 2017-1-20 22:12
实验上对应过程的温度大概多少。

40-70℃左右，还有请问下一般可以自发发生的反应，算出来的能量大约多少呢

小巨人 · 发表于 Post on 2017-1-21 22:53:59

sobereva 发表于 2017-1-21 10:50
M062X和B3LYP算热力学数据本来就相差很大，况且基组还不同，更是雪上加霜。
如果是有机体系，直接用M062X/ ...

谢谢sob老师的解答

liyuanhe211 · 发表于 Post on 2017-1-21 23:01:37

小巨人发表于 2017-1-21 22:49
40-70℃左右，还有请问下一般可以自发发生的反应，算出来的能量大约多少呢

你的势垒即使按计算所得能垒较小的的M06-2X来看也有40kcal/mol，根本不可能在70°C发生（合理的时间尺度下），我认为应考虑假想的反应机理或所建模型有误。

horald · 发表于 Post on 2017-1-26 19:57:31

本帖最后由 horald 于 2017-1-27 05:44 编辑

就计算来讲，不同的方法给出不一致的结果是很可能的。假设有两个人做不同的计算，只要你能说得有理, 都可以发文章。纯粹的计算也就没多大意思。所以有实验数据可以对比参考会比较好。
现在既然你有两组不同的结果，弄清楚为什么还是比较重要的。b3lyp-d 和m06-2x(d) 不能说哪个就不好，你计算的条件也不一样。如你不能反驳b3lyp-d 的结果，就选m06-2x的结果，也可能你选错了，与试验不符。我觉得你可以做些对比性的计算来帮助你确定最有可能的路径: 改变basis sets, reactants ..... 单个的计算数据意义也不大。找出多个数据的共同点（or trend) 会更有说服力。

kevin · 发表于 Post on 2017-2-8 16:53:03

楼主，我也遇到了和你相同的问题，m062x比b3lyp（基组一致）要高7kcal/mol之多，m062x的数据高到我都不敢用，但是又都说m062x算能垒要好一些，还真是挺烦人的，请问你最终怎么处理这个问题呢？

sobereva · 发表于 Post on 2017-2-8 17:10:43

kevin 发表于 2017-2-8 16:53
楼主，我也遇到了和你相同的问题，m062x比b3lyp（基组一致）要高7kcal/mol之多，m062x的数据高到我都不敢用 ...

用双杂化来进一步验证。
我算过的体系M062X算势垒的结果明显比B3LYP更接近于双杂化泛函

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 改变方法基组能量竟差了10kcal！

re

浏览过的版块