资源描述:
《GAMESS程序-能量分解方案》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第十讲:GAMESS程序-能量分解方案和DRC计算2008年5月20日K.Morokuma,J.Chem.Phys.,55,1236(1971)W.Chen,M.S.Gordon,J.Phys.Chem.,100,14316(1996).作用能定义ES:ClassicalCoulumbicinteractionsofOCC(A)andOCC(B)EX:exchangeinteraction,theelectrondelocalizationbetweenOCC(A)andOCC(B)PL:polarizationinterac
2、tionbetweenOCC(A)andVIR(A)orOCC(B)andVIR(B)CT:chargetransferbetweenOCC(A)andVIR(B)orOCC(B)andVIR(A)MIX:theremainingtermsMorokuma原方案各能量项EES=E1-E0ECPL=E2-E1EEX=E3-E1ECT=E2+E3-E1-E4分解方案ES的计算:PL的计算EX的计算CT的计算作用能的计算[]轨道不迭代$MOROKMBSSE=.F.MOROKM=.T.RVS=.F.IATM(1)=ICHM(1)CTP
3、SPL=.T.K.Morokuma,J.Chem.Phys.,55,1236(1971)W.Chen,M.S.Gordon,J.Phys.Chem.,100,14316(1996).输入文件和输出文件输入文件输出结果第十一讲:势能面的交叉问题N.Matsunagaetal.JCP,104,7988-7996(1996)势能面相交的一般原则双原子分子只有不同对称性或不同自旋多重度的势能面多原子体系不受以上的限制GAUSSIAN程序的处理#pcas(6,6,slaterdet)/6-31G**opt=conicaloptcyc=1
4、00ch2s-t01ch1r1h1r12a1r1=1.10a1=99.5结果如何计算-GAMESS程序以CH2(3B1)与CH2(1A1)为例:r(CH)a(HCH)Energy3-B-11.10155130.9955-38.96057261-A-11.1386899.9417-38.9415486R=1.11A=115$contrlscftyp=mcscfmult=1ispher=1$end$basisgbasis=ccd$end$MCSCFCISTEP=GUGA$END$DRTgroup=c2vfors=.t.n
5、mcc=1ndoc=3nval=3istsym=1$end$datamethylene...cc-pVDZbasissetCnv2Carbon6.0Hydrogen1.00.00.936-0.596$end结果单重态$contrlscftyp=mcscfmult=3ispher=1$end$basisgbasis=ccd$end$MCSCFCISTEP=GUGA$END$DRTgroup=c2vfors=.t.nmcc=1ndoc=2nalp=2nval=2istsym=3$end$datamethylene
6、...cc-pVDZbasissetCnv2Carbon6.0Hydrogen1.00.00.936-0.596$end三重态ISTSYM选择12345678C1ACiAgAuCsA’A”C2ABC2vA1A2B1B2C2hAgBuBgAuD2AB1B2B3D2hAgB1gB2gB3gAuB1uB2uB3u求交叉点$contrlruntyp=mexispher=1nzvar=3$end$basisgbasis=ccd$end$guessguess=moreadnorb=24$end$zmatizmat(1)=1,1,2
7、,1,1,3,2,2,1,3$end$mexscf1=mcscfmult1=3nmos1=24scf2=mcscfmult2=1nmos2=24nrdmos=3$end$det1istsym=3ncore=1nact=6nels=6$end$det2istsym=1ncore=1nact=6nels=6$end$datamethylene...cc-pVDZbasissetCnv2Carbon6.0Hydrogen1.00.00.936-0.596$end$VEC1$END$VEC2$END结果$contrlscftyp=m
8、cscfmult=3ispher=1nzvar=3$end$basisgbasis=ccd$end$zmatizmat(1)=1,1,2,1,1,3,2,2,1,3$end$MCSCFCISTEP=GUGA$END$guessguess=moreadnorb=24$e