高斯软件应用.doc

用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 B1 B1 1.102556001 2 3.02-这是个经典的问题，一般的量化计算都会给出不正确的结果。也就是氧原子上带负电。O的电负性是3.44，C的是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.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 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( 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等表达的是分子轨道的对称性信息。是群论中经常用到的对称性的表示方法。 主轴对称，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) 的大作中提到: 】 : 我听一个师兄讲G