计算化学公社

标题: 乙炔分子碳碳三键的解离能该怎么计算呢？ [打印本页]

作者
Author: shuiningzhu 时间: 2023-4-2 09:31
标题: 乙炔分子碳碳三键的解离能该怎么计算呢？
各位老师，乙炔分子，碳碳三键的解离能该怎么计算呢？
我用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

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

作者
Author: shuiningzhu 时间: 2023-4-2 10:16

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

感谢指正。

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