计算化学公社

标题: 键能计算误差过大 [打印本页]

作者
Author: xxtnenu 时间: 2024-5-29 09:15
标题: 键能计算误差过大
请问各位老师，学生使用b3lyp/6-311G(d,p) em=gd3bj进行结构优化，m062x/ma-tzvp计算单点能。使用Shermo程序计算其在60摄氏度下的键解离焓。计算后发现醚键的解离竟高达200多kcal/mol，不知哪个环节出现了问题。

是否将结构优化也使用m062x/ma-tzvp？

作者
Author: KAIMISITERUI 时间: 2024-5-29 15:07
你给的输出文件里直接用b3lyp/6-311G(d,p) em=gd3bj级别下的数据算的键解离焓是合适的，大概率是数据搞错了，你的SP文件里又在新的级别下重新优化了整体的结构，是不是这里搞错了？

作者
Author: KAIMISITERUI 时间: 2024-5-29 15:08
你给的输出文件里直接用b3lyp/6-311G(d,p) em=gd3bj级别下的数据算的键解离焓是合适的，大概率是数据搞错了，你的SP文件里又在新的级别下重新优化了整体的结构，是不是这里搞错了？

作者
Author: xxtnenu 时间: 2024-5-30 10:36

KAIMISITERUI 发表于 2024-5-29 15:08
你给的输出文件里直接用b3lyp/6-311G(d,p) em=gd3bj级别下的数据算的键解离焓是合适的，大概率是数据搞错了 ...

谢谢老师回复，是数据计算搞错了。

作者
Author: xxtnenu 时间: 2024-5-30 10:54

KAIMISITERUI 发表于 2024-5-29 15:08
你给的输出文件里直接用b3lyp/6-311G(d,p) em=gd3bj级别下的数据算的键解离焓是合适的，大概率是数据搞错了 ...

请问老师，学生刚刚又看到另一篇文献计算同一个化合物使用的是M062xdef2tzvp,计算了298.15K下的键解离能，同样位置的C-F的BDE文章计算为113.6Kcal/mol，我进行了构象搜索在333.15K下计算结果是125.5Kcal/mol，相同位置的C-O的BDE它是79.6Kcal/mol，我的是97.05Kcal/mol。这个级别也会有这么大差别吗？学生是否需要也提升结构优化方法基组的级别呢，谢谢老师。

作者
Author: KAIMISITERUI 时间: 2024-5-30 11:21

xxtnenu 发表于 2024-5-30 10:54
请问老师，学生刚刚又看到另一篇文献计算同一个化合物使用的是M062xdef2tzvp,计算了298.15K下的键解离能 ...

建议你用文献中的结构先复现一下他的结果，然后在相同的条件下换成你现在使用的级别，这样才能确定是否是计算级别的影响

作者
Author: xxtnenu 时间: 2024-5-30 11:56

KAIMISITERUI 发表于 2024-5-30 11:21
建议你用文献中的结构先复现一下他的结果，然后在相同的条件下换成你现在使用的级别，这样才能确定是否是 ...

谢谢老师，我尝试一下。

作者
Author: xxtnenu 时间: 2024-5-31 09:42

KAIMISITERUI 发表于 2024-5-30 11:21
建议你用文献中的结构先复现一下他的结果，然后在相同的条件下换成你现在使用的级别，这样才能确定是否是 ...

老师您好，由于大部分文献都没有上传计算结构的坐标，因此我从其他文献找到结构并摆放近似的位置，以m062x/def2-tzvp进行结构优化和计算单点能，但得到的结果仍有较大差距。

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