计算化学公社

标题: gaussian结构优化计算二硫键键能与C-S键能复现文献问题 [打印本页]

作者
Author:
LBH0532    时间: 2023-4-2 22:13
标题: gaussian结构优化计算二硫键键能与C-S键能复现文献问题
各位老师好,我正在参考Density functional theory calculations on S―S bond dissociation energies of disulfides
Yang, Y.-M., Yu, H.-Z., Sun, X.-H., and Dang, Z.-M. (2016) Density functional theory calculations on S―S bond dissociation energies of disulfides. J. Phys. Org. Chem., 29: 6– 13. doi: 10.1002/poc.3480.
这篇文章,考察文中二硫键键能与C-S键能的关系,具体方法是求2倍自由基和分子的自由能之差,文中SI提到c26这个分子根据M06-2X/6-311++G(2df,2p)//M06-2X/6-31G(d)
单点能计算-1938.074444,但是我从几何优化到单点能全部使用gaussian的M06-2X/6-311++G(2df,2p)算得-1938.0701233 (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 1)

(, 下载次数 Times of downloads: 3)

在自由基能量接近的情况下这个小数点后第三位的差值乘以hartree到kcal/mol的627导致的键能差别就很可观了。我现在有三个问题
1.文章中同时使用大小不一的两个基组是否合适?
2.我自己的输入文件是否存在一些问题导致计算出现偏差?
3.换用orca等其他量化软件能否解决这种问题?

作者
Author:
fatpig    时间: 2023-4-3 10:02
可以先按文献方法重复对方的数据,通常来讲你用大基组优化后的能量应该是低于对方这种只算一个单点能的...
作者
Author:
sobereva    时间: 2023-4-3 10:06
仔细看此文
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387http://bbs.keinsci.com/thread-6600-1-1.html

跟程序无关

确保几何结构一样、计算设置精确一样(并且注意看下文注意程序版本的差异)、数据读取无误,并且文章本身的数据正确,才可能保证文章里的数据能重现
谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573http://bbs.keinsci.com/thread-20122-1-1.html





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3