用Gaussian计算CO的偶极矩方向

用Gaussian计算CO的偶极矩方向

如果用Gaussian简单做一下CO的计算,会发现g03得到的电荷分布是氧上带负电荷。但是实际情况是氧上面应该是正电荷，C原子上带负电荷.

输入文件如下:

-----------------------------------------------------

%chk=CO6311G.chk

%mem=6MW

%nproc=1

# hf/6-311++g(3df,3pd) geom=connectivity pop=full

CO opt 6-311G++(3df,3pd)

0 1

C

O 1 B1

B1 1.10255600

1 2 3.0

2

------------------------------------------------------

这是个经典的问题，一般的量化计算都会给出不正确的结果。也就是氧原子上带负电。O的电负性是3.44，C的是2.55，按常理，O要电子的能力强，所以带负电。但其实这是不正确的。

碳原子上带负电荷是从CO分子成键的特点分析得到的，一般认为CO中的一个π是由O原子提供电子对形成的配位键，基于这一点，尽管O的电负性比C大，但C原子因为接受了O的电子对而显负电。这个概念主要用来说明在羰基配合物中为什么配位原子是C而不是O。

根据文献启发，发现使用DFT，如下基组，可以算出Mulliken分析中，C带负电荷，但是用NBO或者AIM，都还是C带正电荷。

C 0

6-311+G(d)

****

O 0

6-311+G

****

O 0

D 1 1.00 0.000000000000

0.0100000D+00 0.1000000000D+01

****