《Journal of Molecular Graphics and Modelling》：COMPLETE SET OF ISODESMIC REACTIONS FOR THE ACCURATE THERMOCHEMISTRY OF THREE- AND FOUR-ATOM CNO-COMPOUNDS
编辑推荐：
完整等键反应集（CS IDR）的方法学被开发并应用于一个包含154个有机小分子的测试集，包括三原子和四原子线性和环状含CNO结构，以评估它们的标准生成焓（ΔfH°）。CS IDR方法使用十种量子化学方法进行了测试：七种DFT近似、DLPNO-CCSD(T)以及
完整等键反应集（CS IDR）的方法学被开发并应用于一个包含154个有机小分子的测试集，包括三原子和四原子线性和环状含CNO结构，以评估它们的标准生成焓（ΔfH°）。CS IDR方法使用十种量子化学方法进行了测试：七种DFT近似、DLPNO-CCSD(T)以及两种复合方法G4和W1BD。研究发现，DFT方法在估算ΔfH°方面不能提供高精度。就计算质量而言，从头算方法排序为G4 > W1BD > DLPNO-CCSD(T)，其精度分别由平均绝对偏差（MAD）1.29、1.52和2.24 kJ mol-1表征。对计算结果的分析揭示，复合方法具有一个系统误差，其大小与参考反应的焓成比例。消除复合方法的系统误差并改善计算与实验的一致性可以通过使用根据方程ΔfH° = 0.5×(ΔfH°(G4) + ΔfH°(W1BD)) – 0.78 kJ mol-1的平均估计来实现。
量子化学计算在热化学研究中占据关键地位，尤其对于有机化合物标准生成焓（ΔfH°）的精确估算。等键反应（isodesmic reactions）和同键反应（homodesmotic reactions）是经典的计算策略，其中完整同键反应集（Complete Set of Homodesmotic Reactions, CS HDR）方法通过比较结构相似片段来消除电子相关误差，已成功应用于多种有机体系，但其适用性受限于分子内无显著π电子重叠的饱和结构。对于高度不饱和或共轭化合物（如多烯、芳烃等），同键反应的局部电子分布假设不再成立，导致计算精度下降。同时，现有热化学数据库（如ATcT和NIST）中参考ΔfH°值数量有限，尤其缺乏小型CNO化合物的可靠数据，这限制了等键或同键方法的应用范围。为此，研究人员开发了完整等键反应集（Complete Set of Isodesmic Reactions, CS IDR）方法，将参考反应相似度从同键水平降低至等键水平，并利用键分离反应（bond-separation reaction）构造完整集合，从而对三原子和四原子CNO化合物（共154个分子）的ΔfH°进行系统评估。研究对比了七种密度泛函理论（DFT）近似、DLPNO-CCSD(T)以及两种复合方法（G4和W1BD），旨在检验CS IDR在不同量子化学层次下的准确性，并识别复合方法中的系统误差。结果表明，DFT方法在ΔfH°估算中精度不足，而复合方法虽然在绝对偏差上表现较好，但存在与参考反应焓呈正比的系统误差。通过采用平均估计方程ΔfH° = 0.5×(ΔfH°(G4) + ΔfH°(W1BD)) – 0.78 kJ mol
-1，可有效消除该误差，使计算值与实验值高度吻合，其平均绝对偏差（Mean Absolute Deviation, MAD）降至约1.29 kJ mol
-1。这一工作发表于《Journal of Molecular Graphics and Modelling》，为热化学数据匮乏的小型CNO化合物提供了可靠的ΔfH°估算途径，同时完善了等键反应方法学的理论基础。
**Methods and calculations**
研究人员采用的关键技术方法包括：基于键分离反应构造完整等键反应集（CS IDR）的图形化算法，将目标分子分解为更小的结构片段并确保原子和价环境平衡；利用量子化学软件Gaussian 09和Orca 6.0.1进行全优化和频率计算，确认势能面极小点；使用七个DFT泛函、DLPNO-CCSD(T)以及复合方法G4和W1BD计算电子能量；通过MAD评估方法精度，并分析系统误差与参考反应焓的线性关系，提出校正方程。样本为154个三原子和四原子线性和环状CNO化合物，其中107个ΔfH°值来自ATcT数据库1.220版，用于对比验证。
**Results and discussion**
在“Results and discussion”部分，研究人员首先指出CS HDR方法不适用于高度不饱和或共轭化合物，因而采用CS IDR变体。通过计算全部154个化合物的ΔfH°，并与ATcT参考值对比，得出以下结论：DFT方法的MAD普遍较大，精度不足；在从头算方法中，计算质量排序为G4 > W1BD > DLPNO-CCSD(T)，对应的MAD分别为1.29、1.52和2.24 kJ mol
-1。进一步分析发现，G4和W1BD复合方法存在系统误差，其大小与参考反应的焓值呈正比。通过采用平均估计方程ΔfH° = 0.5×(ΔfH°(G4) + ΔfH°(W1BD)) – 0.78 kJ mol
-1，可以显著消除该误差，使MAD进一步降低，计算值与实验值的一致性显著改善。这一发现表明，虽然CS IDR在等键水平上降低了反应相似度，但结合适当的系统误差校正，仍能获得接近复合方法最高水平的精度。
**Conclusion**
一个“理想的”热中性参考反应不需要量子化学支持，正如一个“理想的”量子化学方法，不受模型近似引起的误差影响，对参考反应类型不敏感。然而在现实世界中，用户必须做出合理选择，不仅要在计算资源与所需精度之间权衡，还要选择合适的参考过程。对于高度不饱和或共轭化合物，只能被迫降低参考反应的等键水平，但通过CS IDR和系统误差校正，仍可获得满意的热化学数据。
