@走走 又一次更新了本文，增加了CCSD(T)/cc-pVTZ和CCSD(T)/Aug-cc-pVTZ下的结果。

更新了本文，取消了冻核近似(对后HF不计算内核电子相关能的近似)，因为不是横向对比，原理上应该会更好一些。未来可能还会补充一下CCSD(T)的结果。

学过基础无机化学的同学们都知道，实验上一氧化碳分子偶极矩很弱，而且是（正指向负）是氧指向碳的，就是说氧这边带微弱正电，碳那边带微弱负电。之前有同学问过，为什么自己计算出来的结果不符合，这里简单计算一下。计算使用最常见的量子化学计算软件Gaussian，首先要清楚，Gaussian输出的偶极矩方向是【负指向正】的，是与常规定义相反的，因此结果必须取相反数，这是有些同学数据处理错误的重要原因。一氧化碳为C∞v点群分子，程序的标准取向会将分子摆在z轴上并让分子的质心在坐标原点。因此只要读输出的偶极矩在z轴上分量即可，另两个分量一定是0。再强调一点，有些同学用计算出来的Mulliken电荷布居来判断偶极矩的大小和方向，这是不妥的，加了弥散函数的Mulliken电荷布居意义非常差，而且Mulliken电荷布居只是一个很随意的经验模型，应该结合电荷密度分布直接计算偶极矩更为妥当。计算采用的方法分别有Hartree-Fock(HF)，Moller-Plesset的二阶微扰(MP2)和不考虑三重激发的四阶微扰(MP4(SDQ))，单激发和双重激发的耦合簇(CCSD)，普通杂化泛函中普遍认为计算偶极矩非常好的PBE0和Gaussian支持的两种双杂化泛函B2PLYP-D3(BJ)和DSD-PBEP86-D3(BJ)。基组分别采用Dunning的3-zeta相关一致性基组(CC-pVTZ)和对其每个角动量加弥散函数的版本(Aug-CC-pVTZ)。首先简单介绍一下，计算量大致是HF≈普通杂化泛函