gaussian结构优化计算二硫键键能与C-S键能复现文献问题

LBH0532 · 发表于 Post on 2023-4-2 22:13:16

各位老师好，我正在参考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

c26.gjf (1.77 KB, 下载次数 Times of downloads: 1)

poc3480-sup-0001-supplementary.docx (1.33 MB, 下载次数 Times of downloads: 3)

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

fatpig · 发表于 Post on 2023-4-3 10:02:49

可以先按文献方法重复对方的数据，通常来讲你用大基组优化后的能量应该是低于对方这种只算一个单点能的...

sobereva · 发表于 Post on 2023-4-3 10:06:12

仔细看此文
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387（http://bbs.keinsci.com/thread-6600-1-1.html）

跟程序无关

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

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

[Gaussian/gview] gaussian结构优化计算二硫键键能与C-S键能复现文献问题

浏览过的版块