第三章Gaussian軟件的使用new

1、Gaussian98/03的功能 分子的能量和結(jié)構(gòu) 過渡態(tài)的能量和結(jié)構(gòu) 振動頻率 紅外和拉曼光譜（包括預(yù)共振拉曼） 熱化學(xué)性質(zhì) 成鍵和化學(xué)反應(yīng)能量 化學(xué)反應(yīng)路徑 分子軌道 原子電荷 電多極矩 旋光性 振動-轉(zhuǎn)動耦合 非諧性的振動分析和振動-轉(zhuǎn)動耦合 電子親和能和電離勢 極化和超極化率（靜態(tài)的和含頻的） 靜電勢和電子密度 NMR 屏蔽和磁化系數(shù) 自旋-自旋耦合常數(shù)幾何構(gòu)型的構(gòu)建幾何構(gòu)型的構(gòu)建GaussianView例例1：使用使用HF方法，優(yōu)化方法，優(yōu)化 H2O2分子分子%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculatio

2、n0 1 HO 1 0.9O 2 1.4 1 105.0H 3 0.9 2 105.0 1 120.0%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1H 0.000 0.000 0.000 O 0.000 0.900 0.000 O 1.350 1.262 0.000 H 1.464 1.742 -0.752 內(nèi)坐標(biāo)表示內(nèi)坐標(biāo)表示 笛卡爾直角坐標(biāo)表示笛卡爾直角坐標(biāo)表示上述兩個作業(yè)將對上述兩個作業(yè)將對H2O2分子的結(jié)構(gòu)進行完全優(yōu)化，包括所有的鍵長鍵角和二面角。分子的結(jié)構(gòu)進行完全優(yōu)化，包括所有的鍵長鍵角和二面角。

3、通過將內(nèi)坐標(biāo)定義成變量，可以對分子結(jié)構(gòu)進行部分優(yōu)化通過將內(nèi)坐標(biāo)定義成變量，可以對分子結(jié)構(gòu)進行部分優(yōu)化%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1 HO 1 B1O 2 B2 1 A1H 3 B1 2 A1 1 D1B1=0.9 B2=1.4A1=105.0D1=120.0 通過對兩個鍵長和鍵角使用同通過對兩個鍵長和鍵角使用同一變量定義可以控制分子的對一變量定義可以控制分子的對稱性稱性%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculati

4、on0 1 HO 1 B1O 2 B2 1 A1H 3 B1 2 A1 1 D1B1=0.9A1=105.0 變量部分變量部分 D1=120.0B2=1.4 常量部分常量部分通過設(shè)定常量，可以使程序通過設(shè)定常量，可以使程序只對分子結(jié)構(gòu)進行只對分子結(jié)構(gòu)進行部分部分優(yōu)化優(yōu)化例例2： C2H4分子分子 平面形分子平面形分子# HF/STO-3G OPTC2H4 opt0 1CC 1 B1H 1 B2 2 A1H 1 B2 2 A1 3 180.0H 2 B2 1 A1 3 0.0H 2 B2 1 A1 4 0.0B1=1.32B2=1.09A1=120.0 通過通過同一個變量同一個變量控制控制C-H

5、鍵長鍵長 把二面角定義為把二面角定義為180和和0的的常量常量 來控制分子的平面構(gòu)型來控制分子的平面構(gòu)型 例例3： CH3F分子分子 分子點群分子點群C3V # HF/STO-3G OPTCH3F C3v opt0 1CF 1 B1H 1 B2 2 A1H 1 B2 2 A1 3 D1H 1 B2 2 A1 3 D1B1=1.38B2=1.09A1=110.6D1=120.0 通過常量D1來控制C3V對稱性例例4： NH3分子分子 分子點群分子點群C3V # HF/STO-3G OPTNH3 C3v opt0 1NH 1 B1H 1 B1 2 A1H 1 B1 2 A1 3 D1B1=0.99

6、A1=109.0D1=120.0 例例5： 呋喃分子呋喃分子 分子點群分子點群C2v 1.431.341.54 例例6： 乙炔分子乙炔分子 直線形分子直線形分子練習(xí)練習(xí) 寫出分子寫出分子 H2C=C=O分子輸入坐標(biāo)分子輸入坐標(biāo) 其中其中C=C: 1.35 C=O: 1.20 C-H: 1.09 寫出苯分子內(nèi)坐標(biāo)寫出苯分子內(nèi)坐標(biāo) C=C: 1.30 C-H: 1.09Gaussian的文件 Gjf 計算任務(wù)文件 Chk 臨時文件 Rwf 讀寫文件 Out 輸出文件Gaussian程序界面和輸入文件的構(gòu)造程序界面和輸入文件的構(gòu)造%chk=water.chk % Section： 行首以行首以%開始

7、，段開始，段后無空行后無空行%rwf=water.rwf#p hf/6-31g scfcyc=250 scfcon=8 Route Section：行首以行首以# #開始，段后加空行開始，段后加空行 Water ennergy title：作業(yè)的簡要描述，段后加空行作業(yè)的簡要描述，段后加空行 0 1 Molecular Specification：O 分子說明部分，段后通常加空行分子說明部分，段后通常加空行H 1 R1H 1 R1 2 a1R1=1.04a1=104.0 Gaussian作業(yè)的格式作業(yè)的格式% Section(link 0) 定義計算過程中的臨時文件定義計算過程中的臨時文件%c

8、hk=name.chk.chk文件在計算中記錄分子幾何構(gòu)型，分子軌道，力常數(shù)矩陣等信息文件在計算中記錄分子幾何構(gòu)型，分子軌道，力常數(shù)矩陣等信息%rwf=name.rwf.rwf文件主要在作業(yè)重起時使用，當(dāng)計算量比較大時，文件主要在作業(yè)重起時使用，當(dāng)計算量比較大時，.rwf文件通常會非常大，此時需要將文件通常會非常大，此時需要將之分割保存之分割保存 內(nèi)存使用控制內(nèi)存使用控制%mem=n 控制運行過程中使用內(nèi)存的大小控制運行過程中使用內(nèi)存的大小,可以以可以以W或者或者MB，GB為單位為單位default：6000000W48MB 綜合考慮到計算的需要和硬件水平，內(nèi)存并非給得越多越好，最有效率的方法

9、是綜合考慮到計算的需要和硬件水平，內(nèi)存并非給得越多越好，最有效率的方法是根據(jù)作業(yè)類型估算所需要內(nèi)存的大小根據(jù)作業(yè)類型估算所需要內(nèi)存的大小Route Section Route Section以以# 開始，開始，# 控制作業(yè)的輸出控制作業(yè)的輸出#N 正常輸出；默認(rèn)正常輸出；默認(rèn) (沒有計算時間的信息沒有計算時間的信息)#P 輸出更多信息。包括每一執(zhí)行模塊在開始和結(jié)束時與計算機系統(tǒng)有關(guān)的各種輸出更多信息。包括每一執(zhí)行模塊在開始和結(jié)束時與計算機系統(tǒng)有關(guān)的各種 信息信息 ( (包括執(zhí)行時間數(shù)據(jù)，以及包括執(zhí)行時間數(shù)據(jù)，以及SCF計算的收斂信息計算的收斂信息) )#T 精簡輸出：只打印重要的信息和結(jié)果。精

10、簡輸出：只打印重要的信息和結(jié)果。 Route Section主要由方法，基組，任務(wù)類型三部分組成主要由方法，基組，任務(wù)類型三部分組成Gaussian程序能完成的任務(wù)類型：程序能完成的任務(wù)類型：Sp單點能量計算（默認(rèn)任務(wù)類型）單點能量計算（默認(rèn)任務(wù)類型）ADMP&BOMP動力學(xué)計算Opt分子幾何結(jié)構(gòu)優(yōu)化分子幾何結(jié)構(gòu)優(yōu)化Force計算核的受力Freq振動分析振動分析Stable波函數(shù)穩(wěn)定性測試Irc反應(yīng)途徑計算反應(yīng)途徑計算Volume計算分子體積IrcMax在指定反應(yīng)途徑上找能量最大值Density=CheckPoint僅計算布居分析Scan勢能面掃描Guess=Only僅作分子軌道初猜P

11、olar極化率和超極化率計算ReArchive從.chk文件中提取存檔Gaussian 程序中的部分關(guān)鍵詞程序中的部分關(guān)鍵詞Guess（初始軌道猜測） Guess Guess=read 要求從chk文件讀取初始軌道猜測 Guess=mix Guess=alterSCF收斂 SCFCYC=n default=64 收斂圈數(shù)opt 優(yōu)化平衡態(tài)優(yōu)化平衡態(tài) opt=z-matrix 內(nèi)坐標(biāo)，直角坐標(biāo)直接是opt 優(yōu)化過渡態(tài)優(yōu)化過渡態(tài) opt=(ts,z-matrix,noeigentest,readfc) optcyc=nfreq Freq 頻率分析計算常用的量子化學(xué)計算方法量子力學(xué)理論Born-Op

12、penheimer近似非相對論近似單電子近似Hartree-Fock方方 程程Roothaan方方 程程自洽場從頭算SCF-ab initio密度泛函法DFT超 HFLCMTO-X耦合電子對CEPA組 態(tài)相互作用CI微擾處理MP多組態(tài)自洽場MCSCF價電子從頭算EP(VP)模擬從頭算SAMO分子碎片法MF梯度近似GGA浮動球高斯法FSGOAM1C-EHMOEHMOIT-EHMOMCNDOCNDOMINDOINDOMNDONDDOPM3MSW-XDV-XLCAO-X局域密度近 似LDA從頭算法Ab Initio半從頭算法Slater X半經(jīng)驗法Semi-emperical獨立電子對IEPA第一原

13、理計算常用的量子化學(xué)計算方法量子力學(xué)理論Born-Oppenheimer近似非相對論近似單電子近似Hartree-Fock方方 程程Roothaan方方 程程自洽場從頭算SCF-ab initio密度泛函法DFT超 HFLCMTO-X耦合電子對CEPA組 態(tài)相互作用CI微擾處理MP多組態(tài)自洽場MCSCF價電子從頭算EP(VP)模擬從頭算SAMO分子碎片法MF梯度近似GGA浮動球高斯法FSGOAM1C-EHMOEHMOIT-EHMOMCNDOCNDOMINDOINDOMNDONDDOPM3MSW-XDV-XLCAO-X局域密度近 似LDA從頭算法Ab Initio半從頭算法Slater X半經(jīng)驗

14、法Semi-emperical獨立電子對IEPA基組的選擇 Route Section采用自由格式，大小寫不敏感采用自由格式，大小寫不敏感 同一行不同項之間可以使用空格，逗號和同一行不同項之間可以使用空格，逗號和“/”連接；連接； 例：例：#p hf/6-31g scfcyc=230 scfcon=8 #p,hf/6-31g,scfcyc=230,scfcon=8 關(guān)鍵詞可以通過關(guān)鍵詞可以通過 keyword =option,keyword(option), keyword=(option1, option2, .)，keyword( option1, option2, .)指定指定 例：例：

15、#p HF/STO-3G opt #p HF/STO-3G opt=(TS,readfc) Title部分必須輸入，但是程序并不執(zhí)行，起標(biāo)識和說明作用部分必須輸入，但是程序并不執(zhí)行，起標(biāo)識和說明作用Title SectionCharge & Multipl. 輸入分子的電荷和自旋多重度輸入分子的電荷和自旋多重度 例：例： 電荷電荷 多重度多重度2s+1 H2O 0 1 H3O+ 1 1 NO 0 2 電荷多重度部分通常也算作分子說明電荷多重度部分通常也算作分子說明Molecular Specification 分子說明部分主要用來定義原子核相對位置分子說明部分主要用來定義原子核相對位置

16、原子核相對位置可以用笛卡爾坐標(biāo)，內(nèi)坐標(biāo)原子核相對位置可以用笛卡爾坐標(biāo)，內(nèi)坐標(biāo)(Z-matrix)，或者是二者混合，或者是二者混合 表示表示笛卡爾坐標(biāo)是內(nèi)坐標(biāo)的一種特殊形式笛卡爾坐標(biāo)是內(nèi)坐標(biāo)的一種特殊形式分子坐標(biāo)的格式為分子坐標(biāo)的格式為: 元素符號, x, y, z 元素符號 (n) 原子1 鍵長 原子2 鍵角 原子3 二面角 (0,180) 以原子在分子中的序數(shù)表示例例1：使用使用HF方法，優(yōu)化方法，優(yōu)化 H2O2分子分子%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1 HO 1 0.9O 2 1.4 1 1

17、05.0H 3 0.9 2 105.0 1 120.0%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1H 0.000 0.000 0.000 O 0.000 0.900 0.000 O 1.350 1.262 0.000 H 1.464 1.742 -0.752 內(nèi)坐標(biāo)表示內(nèi)坐標(biāo)表示 笛卡爾直角坐標(biāo)表示笛卡爾直角坐標(biāo)表示上述兩個作業(yè)將對上述兩個作業(yè)將對H2O2分子的結(jié)構(gòu)進行完全優(yōu)化，包括所有的鍵長鍵角和二面角。分子的結(jié)構(gòu)進行完全優(yōu)化，包括所有的鍵長鍵角和二面角。通過將內(nèi)坐標(biāo)定義成變量，可以對分子結(jié)構(gòu)進行部分優(yōu)

18、化通過將內(nèi)坐標(biāo)定義成變量，可以對分子結(jié)構(gòu)進行部分優(yōu)化%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1 HO 1 r1O 2 r2 1 a1H 3 r1 2 a1 1 d1r1=0.9 r2=1.4a1=105.0d1=120.0 通過對兩個鍵長和鍵角使用同通過對兩個鍵長和鍵角使用同一變量定義可以控制分子的對一變量定義可以控制分子的對稱性稱性%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1 HO 1 r1O 2 r2 1 a1

19、H 3 r1 2 a1 1 d1r1=0.9a1=105.0 變量部分變量部分 d1=120.0r2=1.4 常量部分常量部分通過設(shè)定常量，可以使程序通過設(shè)定常量，可以使程序只對分子結(jié)構(gòu)進行只對分子結(jié)構(gòu)進行部分部分優(yōu)化優(yōu)化例例2： C2H4分子分子 平面形分子平面形分子# HF/STO-3G OPTC2H4 opt0 1CC 1 r1H 1 r2 2 a1H 1 r2 2 a1 3 180.0H 2 r2 1 a1 3 0.0H 2 r2 1 a1 4 0.0r1=1.32r2=1.09a1=120.0 通過通過同一個變量同一個變量控制控制C-H鍵長鍵長 把二面角定義為把二面角定義為180和和

20、0的的常量常量 來控制分子的平面構(gòu)型來控制分子的平面構(gòu)型 例例3： CH3F分子分子 分子點群分子點群C3V # HF/STO-3G OPTCH3F C3v opt0 1CF 1 r1H 1 r2 2 a1H 1 r2 2 a1 3 -bH 1 r2 2 a1 3 br1=1.38r2=1.09a1=110.6b=120.0 通過常量b來控制C3V對稱性例例4： NH3分子分子 分子點群分子點群C3V # HF/STO-3G OPTNH3 C3v opt0 1NH 1 r1H 1 r1 2 a1H 1 r1 2 a1 3 br1=0.99a1=109.0b=120.0運行界面GaussianO

21、verlay1Overlay9,10,11,99L101L102L122Overlay0L0L001 Gaussian 程序的結(jié)構(gòu)程序的結(jié)構(gòu)link0: 初始化程序，控制初始化程序，控制占位占位link1: 讀入并處理讀入并處理計算執(zhí)行路徑，計算執(zhí)行路徑，建建立要執(zhí)行的立要執(zhí)行的link列表列表link9999: 終止計算終止計算Overlay99L9999Gaussian03 程序界面程序界面Preferences:對對Gaussian程序進行初程序進行初始化設(shè)置始化設(shè)置自定義外部文字編輯器，用來打開自定義外部文字編輯器，用來打開.out文文件件link.exe所在的文件夾所在的文件夾臨時文

22、件存放文件夾臨時文件存放文件夾Gaussian03 圖形工具欄圖形工具欄開始作業(yè)開始作業(yè) 暫停當(dāng)前作業(yè)暫停當(dāng)前作業(yè) 當(dāng)前當(dāng)前l(fā)ink后暫停后暫停 終止當(dāng)前作業(yè)和批處理終止當(dāng)前作業(yè)和批處理 恢復(fù)當(dāng)前作業(yè)恢復(fù)當(dāng)前作業(yè) 在當(dāng)前作業(yè)完成后終止批處理在當(dāng)前作業(yè)完成后終止批處理 終止當(dāng)前作業(yè)終止當(dāng)前作業(yè) 編輯或建立批處理編輯或建立批處理 批處理專用批處理專用打開外部編輯器打開外部編輯器 編輯輸出文件編輯輸出文件-Link1- 指定多步任務(wù)Gaussian作業(yè)輸出示例作業(yè)輸出示例 （H2O SP） Entering Link 1 = C:G03Wl1.exe PID= 1572. Copyright (c)

23、 1988,1990,1992,1993,1995,1998,2003, Gaussian, Inc. All Rights Reserved. This is the Gaussian(R) 03 program. It is based on the the Gaussian(R) 98 system (copyright 1998, Gaussian, Inc.), the Gaussian(R) 94 system (copyright 1995, Gaussian, Inc.), the Gaussian 92(TM) system (copyright 1992, Gaussian

24、, Inc.), the Gaussian 90(TM) system (copyright 1990, Gaussian, Inc.), the Gaussian 88(TM) system (copyright 1988, Gaussian, Inc.), the Gaussian 86(TM) system (copyright 1986, Carnegie Mellon University), and the Gaussian 82(TM) system (copyright 1983, Carnegie Mellon University). Gaussian is a feder

25、ally registered trademark of Gaussian, Inc.RESTRICTED RIGHTS LEGEND Use, duplication or disclosure by the US Government is subject to restrictions as set forth in subparagraph (c)(1)(ii) of the Rights in Technical Data and Computer Software clause at DFARS 252.227-7013. 歷歷史史和和版版權(quán)權(quán)Cite this work as:

26、Gaussian 03, Revision B.02, M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, J. A. Montgomery, Jr., T. Vreven, K. N. Kudin, J. C. Burant, J. M. Millam, S. S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G. A. Petersson, H. Nakats

27、uji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, J. E. Knox, H. P. Hratchian, J. B. Cross, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, P. Y. Ayal

28、a, K. Morokuma, G. A. Voth, P. Salvador, J. J. Dannenberg, V. G. Zakrzewski, S. Dapprich, A. D. Daniels, M. C. Strain, O. Farkas, D. K. Malick, A. D. Rabuck, K. Raghavachari, J. B. Foresman, J. V. Ortiz, Q. Cui, A. G. Baboul, S. Clifford, J. Cioslowski, B. B. Stefanov, G. Liu, A. Liashenko, P. Pisko

29、rz, I. Komaromi, R. L. Martin, D. J. Fox, T. Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, M. Challacombe, P. M. W. Gill, B. Johnson, W. Chen, M. W. Wong, C. Gonzalez, and J. A. Pople, Gaussian, Inc., Pittsburgh PA, 2003.引引用用信信息息* Gaussian 03: x86-Win32-G03RevB.02 16-Apr-2003 08-Apr-2005 * %chk

30、=h2o.chk %rwf=h2o.rwf Default route: MaxDisk=2000MB - #p hf/6-31g - 1/38=1/1; 2/17=6,18=5,40=1/2; 3/5=1,6=6,11=9,16=1, 25=1,30=1/1,2,3; 4/1; 5/5=2,32=1,38=5/2; 6/7=2,8=2,9=2,10=2,28=1/1; 99/5=1,9=1/99; Leave Link 1 at Fri Apr 08 01:21:37 2005, MaxMem= 0 cpu: 1.0 重重新新打打印印作作業(yè)業(yè)輸輸入入和和標(biāo)標(biāo)準(zhǔn)準(zhǔn)途途徑徑作業(yè)執(zhí)行過程作業(yè)執(zhí)行過

31、程1/38=1/1; 2/17=6,18=5,40=1/2;3/5=1,6=6,11=9,16=1,25=1,30=1/1,2,3;4/1;5/5=2,32=1,38=5/2;6/7=2,8=2,9=2,10=2,28=1/1;99/5=1,9=1/99;L101L202L301L302L303L401L502L601L9999L101 讀取標(biāo)題和分子輸入部分； L202 坐標(biāo)重新定向，計算對稱性，檢查變量；L301 產(chǎn)生基組信息； L302 計算重疊積分，動能積分和勢能積分；L303 計算多極積分； L401 形成分子軌道初猜； L502 自洽場迭代求解；L601 布居數(shù)和相關(guān)分析；L999

32、9 完成計算和輸出L101 讀取標(biāo)題和分子輸入部分讀取標(biāo)題和分子輸入部分(Enter C:G03Wl101.exe) - h2o sp - Symbolic Z-matrix: Charge = 0 Multiplicity = 1 O 0. 0. 0. H 0.662 0. 0.468 H -0.662 0. 0.468 Leave Link 101 at Fri Apr 08 01:21:38 2005, MaxMem= 6291456 cpu: 0.0L202 坐標(biāo)重新定向，計算對稱性，檢查變量坐標(biāo)重新定向，計算對稱性，檢查變量(Enter C:G03Wl202.exe) Input o

33、rientation: - Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z - 1 8 0 0.000000 0.000000 0.000000 2 1 0 0.662000 0.000000 0.468000 3 1 0 -0.662000 0.000000 0.468000 - Distance matrix (angstroms): 1 2 3 1 O 0.000000 2 H 0.810721 0.000000 3 H 0.810721 1.324000 0.000000 Stoichiomet

34、ry H2O Framework group C2VC2(O),SGV(H2) Deg. of freedom 2 Full point group C2V NOp 4 Largest Abelian subgroup C2V NOp 4 Largest concise Abelian subgroup C2 NOp 2 Standard orientation: - Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z - 1 8 0 0.000000 0.000000 0.093600 2 1 0 0.0

35、00000 0.662000 -0.374400 3 1 0 0.000000 -0.662000 -0.374400 - Rotational constants (GHZ): 1289.0084222 572.1186153 396.2468433 Leave Link 202 at Fri Apr 08 01:21:40 2005, MaxMem= 6291456 cpu: 1.0L301 產(chǎn)生基組信息產(chǎn)生基組信息(Enter C:G03Wl301.exe) Standard basis: 6-31G (6D, 7F) There are 7 symmetry adapted basis

36、 functions of A1 symmetry. There are 0 symmetry adapted basis functions of A2 symmetry. There are 2 symmetry adapted basis functions of B1 symmetry. There are 4 symmetry adapted basis functions of B2 symmetry. Integral buffers will be 262144 words long. Raffenetti 1 integral format. Two-electron int

37、egral symmetry is turned on. 13 basis functions, 30 primitive gaussians, 13 cartesian basis functions 5 alpha electrons 5 beta electrons nuclear repulsion energy 10.8432720049 Hartrees. IExCor= 0 DFT=F Ex=HF Corr=None ExCW=0 ScaHFX= 1.000000 ScaDFX= 1.000000 1.000000 1.000000 1.000000 IRadAn= 0 IRan

38、Wt= -1 IRanGd= 0 ICorTp=0 NAtoms= 3 NActive= 3 NUniq= 2 SFac= 2.76D+00 NAtFMM= 60 Big=F Leave Link 301 at Fri Apr 08 01:57:33 2005, MaxMem= 6291456 cpu: 1.0L302 計算重疊積分，動能積分和勢能積分計算重疊積分，動能積分和勢能積分L303 計算多極積分計算多極積分(Enter C:G03Wl302.exe) NPDir=0 NMtPBC= 1 NCelOv= 1 NCel= 1 NClECP= 1 NCelD= 1 NCelK= 1 NCe

39、lE2= 1 NClLst= 1 CellRange= 0.0. One-electron integrals computed using PRISM. One-electron integral symmetry used in STVInt NBasis= 13 RedAO= T NBF= 7 0 2 4 NBsUse= 13 1.00D-06 NBFU= 7 0 2 4 Leave Link 302 at Fri Apr 08 01:57:35 2005, MaxMem= 6291456 cpu: 1.0 (Enter C:G03Wl303.exe) DipDrv: MaxL=1. L

40、eave Link 303 at Fri Apr 08 01:57:36 2005, MaxMem= 6291456 cpu: 0.0L401 形成分子軌道初猜形成分子軌道初猜(Enter C:G03Wl401.exe) Harris functional with IExCor= 205 diagonalized for initial guess. ExpMin= 1.61D-01 ExpMax= 5.48D+03 ExpMxC= 8.25D+02 IAcc=1 IRadAn= 1 AccDes= 1.00D-06 HarFok: IExCor= 205 AccDes= 1.00D-06

41、IRadAn= 1 IDoV=1 ScaDFX= 1.000000 1.000000 1.000000 1.000000 Harris En= -76.0428456360852 Initial guess orbital symmetries: Occupied (A1) (A1) (B2) (A1) (B1) Virtual (A1) (B2) (B1) (A1) (B2) (B2) (A1) (A1) The electronic state of the initial guess is 1-A1. Leave Link 401 at Fri Apr 08 01:57:37 2005,

42、 MaxMem= 6291456 cpu: 0.0 L502 自洽場迭代求解自洽場迭代求解 (Enter C:G03Wl502.exe) Warning! Cutoffs for single-point calculations used. Closed shell SCF: Requested convergence on RMS density matrix=1.00D-04 within 128 cycles. Requested convergence on MAX density matrix=1.00D-02. Requested convergence on energy=5.

43、00D-05. No special actions if energy rises. Using DIIS extrapolation, IDIIS= 1040. Integral symmetry usage will be decided dynamically. Keep R1 integrals in memory in canonical form, NReq= 424251. IEnd= 18776 IEndB= 18776 NGot= 6291456 MDV= 6285583 LenX= 6285583 Symmetry not used in FoFDir. MinBra=

44、0 MaxBra= 1 Meth= 1. IRaf= 0 NMat= 1 IRICut= 1 DoRegI=T DoRafI=F ISym2E= 0 JSym2E=0. Cycle 5 Pass 1 IDiag 1: E= -75.9273519590932 Delta-E= -0.000003762089 Rises=F Damp=F DIIS: error= 5.91D-05 at cycle 5 NSaved= 5. NSaved= 5 IEnMin= 5 EnMin= -75.9273519590932 IErMin= 5 ErrMin= 5.91D-05 ErrMax= 5.91D-

45、05 EMaxC= 1.00D-01 BMatC= 4.02D-08 BMatP= 1.89D-06 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: 0.711D-02-0.177D-01 0.103D+00-0.460D+00 0.137D+01 Coeff: 0.711D-02-0.177D-01 0.103D+00-0.460D+00 0.137D+01 Gap= 0.746 Goal= None Shift= 0.000 RMSDP=4.58D-05 MaxDP=1.68D-04 DE=-3.76D-06 OVMax= 2.54D-

46、04 SCF Done: E(RHF) = -75.9273519591 A.U. after 5 cycles Convg = 0.4576D-04 -V/T = 1.9895 S*2 = 0.0000 KE= 7.673560164965D+01 PE=-2.025358608090D+02 EE= 3.902963519532D+01 Leave Link 502 at Fri Apr 08 01:57:41 2005, MaxMem= 6291456 cpu: 2.0 L601 布居數(shù)和相關(guān)分析布居數(shù)和相關(guān)分析(Enter C:G03Wl601.exe) Copying SCF den

47、sities to generalized density rwf, ISCF=0 IROHF=0. * Population analysis using the SCF density. * Orbital symmetries: Occupied (A1) (A1) (B2) (A1) (B1) Virtual (A1) (B2) (B2) (B1) (A1) (A1) (B2) (A1) The electronic state is 1-A1. Alpha occ. eigenvalues - -20.52887 -1.44410 -0.81350 -0.57281 -0.51130

48、 Alpha virt. eigenvalues - 0.23490 0.33109 1.13536 1.16373 1.17135 Alpha virt. eigenvalues - 1.43636 1.50054 1.68694 Condensed to atoms (all electrons): 1 2 3 1 O 8.090714 0.312339 0.312339 2 H 0.312339 0.374540 -0.044575 3 H 0.312339 -0.044575 0.374540 Mulliken atomic charges: 1 1 O -0.715392 2 H 0

49、.357696 3 H 0.357696 Sum of Mulliken charges= 0.00000 Atomic charges with hydrogens summed into heavy atoms: 1 1 O 0.000000 2 H 0.000000 3 H 0.000000 Sum of Mulliken charges= 0.00000Electronic spatial extent (au): = 17.2272 Charge= 0.0000 electrons Dipole moment (field-independent basis, Debye): X=

50、0.0000 Y= 0.0000 Z= -2.4602 Tot= 2.4602 Quadrupole moment (field-independent basis, Debye-Ang): XX= -6.8696 YY= -4.2696 ZZ= -6.1388 XY= 0.0000 XZ= 0.0000 YZ= 0.0000 Traceless Quadrupole moment (field-independent basis, Debye-Ang): XX= -1.1103 YY= 1.4898 ZZ= -0.3795 XY= 0.0000 XZ= 0.0000 YZ= 0.0000 Octapole moment (field-independent basis, Debye-Ang*2): XXX= 0.0000 YYY= 0.0000 ZZZ= -1.6189 XYY= 0.0000 XXY= 0.0000 XXZ= -0.5144 XZZ= 0.0000 YZZ= 0.0000 YYZ= -0.9389 XYZ= 0.0000 Hexadecap