Gaussian计算软件的使用.ppt

Gaussian计算软件的使用 优化时常用的关键词 %chk=td-fang-2-2 #P B3LYP/6-31G* POPT(tight, maxcycle=200,restart) SCF(direct,maxcycle=200) GFINPUT IOP(6/7=3) TEST td at 12.26 0 1 ….. 频率或者NMR计算 #P B3LYP/6-31G* FREQ GFINPUT IOP(6/7=3) TEST #P B3LYP/6-31G* NMR GFINPUT IOP(6/7=3) TEST GEN基组和PSEUDO #p B3LYP/GEN PSEUDO=READ GFinput POPT TEST au-au(PMe3)2 05-06-2003 0 1 au 0.0 0.0 0.0 p 0.0 0.0 1.9 auxx3 1.380000 pxx4 2.440000 cp5 1.950000 P 0 6-31+g(d) **** Au 0 CEP-121G **** Au 0 CEP-121G CIS光谱计算 CIS 单点 #P CIS(singlets,Nstates=30, root=1)/6-31G* TEST 或者 #P HF/6-31G* CIS(singlets,Nstates=30, root=1) TEST CIS 优化 #P CIS(singlets,Nstates=30, root=1)/6-31G* POPT FREQ TEST CIS相关关键词 Singlets 只算单重态 Triplets 只算三重态 50-50 单重和三重态各占一半 Root=N 研究CIS计算的第几个激发态Direct 加快计算速度，减小硬盘使用空间 Restart 从Chk文件中重新开始 CIS 计算 Read 从Chk文件中读 CIS 分析激发态电荷情况 用 Density=current 或者Density=CIS 接CIS优化 # CIS(restart)/3-21g* guess=read geom=check POPT TEST 或者 # CIS/3-21g* POPT=restart TEST TD光谱计算 #P B3LYP/6-31G* TD(singlets,Nstates=30, root=1) TEST 如果中间断电可以用下面来接， #P B3lyp/3-21g* TD=restart guess=read geom=check TEST TD相关关键词 Singlets 只算单重态 Triplets 只算三重态 50-50 单重和三重态各占一半 Root=N 研究TD计算的第几个激发态Direct 加快计算速度，减小硬盘使用空间 Read 从Chk文件中读 TD （G03有，好像G98 A11.2也有） SOS Do sum-over states polarizabilities, etc. Gaussian中NBO计算 # B3LYP/6-31G(d,p) Pop=(NBO, NPA, NBORead) Example of NBO bond orders 0 1 C 0.000000 0.665676 0.000000 H 0.919278 1.237739 0.000000 H -0.919239 1.237787 0.000000 C 0.000000 -0.665676 0.000000 H -0.919278 -1.237739 0.000000 H 0.919239 -1.237787 0.000000 $nbo bndidx file … $end 消除自旋污染 SCF Done: E(UHF) = -111.945340085 A.U. after 19 cycles Convg = 0.3946D-08 -V/T = 2.0071 S**2 = 2.0318 S**2=如果没有自旋污染，它会等于s(s+1)，这里s等于1/2的未成对电子数倍 在Gaussian中，选项“iop(5/14=2)” 告诉程序用湮灭后的波函数产生布局数分析,来降低自旋污染。 Gaussian中的不收敛问题的解决 1)由于体系有很多能量相近的能级，导致计算不收敛。 如果计算