高斯软件应用.doc

用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??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??**** Gaussian中偶极矩的单位的换算 用Gaussian进行极化分子的计算时，最后会输出沿着X,Y,Z方向的偶极矩，但是其单位用的是Debye。??典型输出如下：??---------------------------------------------?? Charge= 0.0000 electrons?? Dipole 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 for two charges of 10–10franklin separated by 1 Angstrom, D=10–18Fr cm ≈ 3.335 64 × 10–30C 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( 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, E1.2等表达的是分子轨道的对称性信息。??是群论中经常用到的对称性的表示方法。 ??主轴