乙炔分子碳碳三键的解离能该怎么计算呢？

shuiningzhu

各位老师，乙炔分子，碳碳三键的解离能该怎么计算呢？
我用b3lyp/6-31g(d)优化，M062x/def2TZVP计算，与实验值差距太大。而其他的烯烃则很接近。



乙炔分子：
%chk=acetylene-opt.chk
# opt freq  b3lyp/6-31g(d)

Title Card Required

0 1
C                 -3.29758700    1.63538872    0.00000000
C                 -2.10258700    1.63538872    0.00000000
H                 -4.35858700    1.63538872    0.00000000
H                 -1.04158700    1.63538872    0.00000000

--link1--
%oldchk=acetylene-opt.chk
%chk=acetylene-sp.chk

# M062x/def2TZVP freq geom=allcheck guess=read


解离片段：
%chk=radical-opt.chk
# opt freq ub3lyp/6-31g(d)

Title Card Required

0 4
C                 -3.29308700    1.63538872    0.00000000
H                 -4.36308700    1.63538872    0.00000000

--link1--
%oldchk=radical-opt.chk
# uM062x/def2TZVP freq geom=allcheck guess=read

zjxitcc

计算严重不合理。至少存在如下问题：
（1）当前计算级别 对计算键能 过于拉胯。体系很小，所以结构优化建议(U)M062X/def2TZVP起步。高精度单点能建议(U)CCSD(T)/cc-pVTZ起步。
（2）切换计算级别后，忘记了结构优化、直接进行freq计算，计算得到的热力学量无意义。freq需与opt在同一计算机别下进行。
（3）贴输入文件时最好采用代码块展示，这样空行数目基本不会有问题。

shuiningzhu

zjxitcc 发表于 2023-4-2 09:55
计算严重不合理。至少存在如下问题：
（1）当前计算级别 对计算键能 过于拉胯。体系很小，所以结构优化建 ...

感谢指正。