1、用 Gaussian 计算 CO 的偶极矩方向如果用 Gaussian 简单做一下 CO 的计算,会发现 g03 得到的电荷分布是氧上带负电荷。但是实际情况是氧上面应该是正电荷，C 原子上带负电荷 .输入文件如下:-%chk=CO6311G.chk%mem=6MW%nproc=1# hf/6-311+g(3df,3pd) geom=connectivity pop=fullCO opt 6-311G+(3df,3pd)0 1CO 1 B1B1 1.102556001 2 3.02-这是个经典的问题，一般的量化计算都会给出不正确的结果。也就是氧原子上带负电。O 的电负性是3.44，C 的是 2.
2、55，按常理，O 要电子的能力强，所以带负电。但其实这是不正确的。碳原子上带负电荷是从 CO 分子成键的特点分析得到的，一般认为 CO 中的一个 是由 O 原子提供电子对形成的配位键，基于这一点，尽管 O 的电负性比 C 大，但 C 原子因为接受了 O 的电子对而显负电。这个概念主要用来说明在羰基配合物中为什么配位原子是 C 而不是 O。 根据文献启发，发现使用 DFT，如下基组，可以算出 Mulliken 分析中，C 带负电荷，但是用 NBO 或者AIM，都还是 C 带正电荷。C 06-311+G(d)*O 06-311+G*O 0D 1 1.00 0.0000000000000.01000
3、00D+00 0.1000000000D+01*Gaussian 中偶极矩的单位的换算用 Gaussian 进行极化分子的计算时，最后会输出沿着 X,Y,Z 方向的偶极矩，但是其单位用的是 Debye。典型输出如下：-Charge= 0.0000 electronsDipole moment (Debye):X= 0.1018 Y= -0.0475 Z= -38.8636 Tot= 38.8638-关于 Debye 的单位定义如下:Non-SI unit of electric dipole moment. It is equal to the electric dipolemoment fo
4、r two charges of 1010franklin separated by 1 Angstrom, D=1018Fr cm 3.335 64 1030C m 或者说：1a.u.=0.97174E22V/m2=2.542Debye 其中 1 a.u.偶极矩定义为1 au = two point electron unit charges seperated by 1 bohr (0.529 angstrom) = 2.54 Debye 在 Gaussian 的计算输出结果中，会有如下的信息-Orbital occupation is: 40 A1( 1) alpha and 39 A1
5、( 1) beta 6 A2( 1,) alpha and 6 A2( 1) beta, and 0.00 alpha elec in 1 orbitals plus 0.00 beta elec in 1 orbitals 36 E1( 2,) alpha and 36 E1( 2) beta, and 0.00 alpha elec in 0 orbitals plus 0.50 beta elec in 1 orbitals 29 E2( 2) alpha and 29 E2( 2) beta Total number electrons: 352.0 -这里的 A1.1, A2.1,
6、E1.2 等表达的是分子轨道的对称性信息。是群论中经常用到的对称性的表示方法。 主轴对称，A，反对称，B 下标 1，垂直于主轴的 C2 对称，反对称，则用 2 括号中 1， （可能是上标） ，表示垂直于主轴的晶面对称，2，反对称 E 是二维的，A，B 是一维的对称不可约表示如何可视化 Gaussian 的分子轨道是利用 gaussian 的 cugegen 的程序。利用计算产生的 CHK 文件，然后 用 formchk 将 chk 转化为.fch 的文件。利用 chem3d 打开，或者 gaussview 打开， 在 chem3d 里面有一个 toolbar 的 view 栏，选择 molecular orbital，即可。 当然可以根据不同需要，如颜色等的设置。 在 gaussview 里面;有一个 results 的 toolbar，选择 surface 即可。 【 在 pandafu (panda) 的大作中提到: 】 : 我听一个师兄讲 Gaussian98 具有以下功能，清高手指点：Draw the LUMO and HOMO as : a three-dimension picture. : 请问 Gaussian98 能做到吗?
