大概这样：
1，做某羰基化合物的opt。或者scan，找到优势结构后用B3LYP/6-311+G（d）算总能E1
2，用opt结构，删除O原子，用单点能算所得三线态卡宾和三线态O原子，这二者总能E2
3，羰基键能=E2-E1
@sugar0763，作为简单算法（比如定性得比较集中C=O趋势）可行否？
我正在从IR的C=O伸缩振动进行对比，但是算能量的基组我实在用不了太大的。机器太烂