第三章Gaussian软件的使用new课件.ppt

《第三章Gaussian软件的使用new课件.ppt》由会员分享，可在线阅读，更多相关《第三章Gaussian软件的使用new课件.ppt（62页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Gaussian98/03的功能 分子的能量和结构 过渡态的能量和结构 振动频率 红外和拉曼光谱（包括预共振拉曼） 热化学性质 成键和化学反应能量 化学反应路径 分子轨道 原子电荷 电多极矩 旋光性 振动-转动耦合 非谐性的振动分析和振动-转动耦合 电子亲和能和电离势 极化和超极化率（静态的和含频的） 静电势和电子密度 NMR 屏蔽和磁化系数 自旋-自旋耦合常数几何构型的构建几何构型的构建GaussianView例例1：使用使用HF方法，优化方法，优化 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 内坐标表示内坐标表示 笛卡尔直角坐标表示笛卡尔直角坐标表示上述两个作业将对上述两个作业将对H2O2分子的结构进行完全优化，包括所有的键长键角和二面角。分子的结构进行完全优化，包括所有的键长键角和二面角。

3、通过将内坐标定义成变量，可以对分子结构进行部分优化通过将内坐标定义成变量，可以对分子结构进行部分优化%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 常量部分常量部分通过设定常量，可以使程序通过设定常量，可以使程序只对分子结构进行只对分子结构进行部分部分优化优化例例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的的常量常量 来控制分子的平面构型来控制分子的平面构型 例例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： 乙炔分子乙炔分子 直线形分子直线形分子练习练习 写出分子写出分子 H2C=C=O分子输入坐标分子输入坐标 其中其中C=C: 1.35 C=O: 1.20 C-H: 1.09 写出苯分子内坐标写出苯分子内坐标 C=C: 1.30 C-H: 1.09Gaussian的文件 Gjf 计算任务文件 Chk 临时文件 Rwf 读写文件 Out 输出文件Gaussian程序界面和输入文件的构造程序界面和输入文件的构造%chk=water.chk % Section： 行首以行首以%开始

7、，段开始，段后无空行后无空行%rwf=water.rwf#p hf/6-31g scfcyc=250 scfcon=8 Route Section：行首以行首以# #开始，段后加空行开始，段后加空行 Water ennergy title：作业的简要描述，段后加空行作业的简要描述，段后加空行 0 1 Molecular Specification：O 分子说明部分，段后通常加空行分子说明部分，段后通常加空行H 1 R1H 1 R1 2 a1R1=1.04a1=104.0 Gaussian作业的格式作业的格式% Section(link 0) 定义计算过程中的临时文件定义计算过程中的临时文件%c

8、hk=name.chk.chk文件在计算中记录分子几何构型，分子轨道，力常数矩阵等信息文件在计算中记录分子几何构型，分子轨道，力常数矩阵等信息%rwf=name.rwf.rwf文件主要在作业重起时使用，当计算量比较大时，文件主要在作业重起时使用，当计算量比较大时，.rwf文件通常会非常大，此时需要将文件通常会非常大，此时需要将之分割保存之分割保存 内存使用控制内存使用控制%mem=n 控制运行过程中使用内存的大小控制运行过程中使用内存的大小,可以以可以以W或者或者MB，GB为单位为单位default：6000000W48MB 综合考虑到计算的需要和硬件水平，内存并非给得越多越好，最有效率的方法

9、是综合考虑到计算的需要和硬件水平，内存并非给得越多越好，最有效率的方法是根据作业类型估算所需要内存的大小根据作业类型估算所需要内存的大小Route Section Route Section以以# 开始，开始，# 控制作业的输出控制作业的输出#N 正常输出；默认正常输出；默认 (没有计算时间的信息没有计算时间的信息)#P 输出更多信息。包括每一执行模块在开始和结束时与计算机系统有关的各种输出更多信息。包括每一执行模块在开始和结束时与计算机系统有关的各种 信息信息 ( (包括执行时间数据，以及包括执行时间数据，以及SCF计算的收敛信息计算的收敛信息) )#T 精简输出：只打印重要的信息和结果。精

10、简输出：只打印重要的信息和结果。 Route Section主要由方法，基组，任务类型三部分组成主要由方法，基组，任务类型三部分组成Gaussian程序能完成的任务类型：程序能完成的任务类型：Sp单点能量计算（默认任务类型）单点能量计算（默认任务类型）ADMP&BOMP动力学计算Opt分子几何结构优化分子几何结构优化Force计算核的受力Freq振动分析振动分析Stable波函数稳定性测试Irc反应途径计算反应途径计算Volume计算分子体积IrcMax在指定反应途径上找能量最大值Density=CheckPoint仅计算布居分析Scan势能面扫描Guess=Only仅作分子轨道初猜Polar

11、极化率和超极化率计算ReArchive从.chk文件中提取存档Gaussian 程序中的部分关键词程序中的部分关键词Guess（初始轨道猜测） Guess Guess=read 要求从chk文件读取初始轨道猜测 Guess=mix Guess=alterSCF收敛 SCFCYC=n default=64 收敛圈数opt 优化平衡态优化平衡态 opt=z-matrix 内坐标，直角坐标直接是opt 优化过渡态优化过渡态 opt=(ts,z-matrix,noeigentest,readfc) optcyc=nfreq Freq 频率分析计算常用的量子化学计算方法量子力学理论Born-Oppenh

12、eimer近似非相对论近似单电子近似Hartree-Fock方方 程程Roothaan方方 程程自洽场从头算SCF-ab initio密度泛函法DFT超 HFLCMTO-X耦合电子对CEPA组 态相互作用CI微扰处理MP多组态自洽场MCSCF价电子从头算EP(VP)模拟从头算SAMO分子碎片法MF梯度近似GGA浮动球高斯法FSGOAM1C-EHMOEHMOIT-EHMOMCNDOCNDOMINDOINDOMNDONDDOPM3MSW-XDV-XLCAO-X局域密度近 似LDA从头算法Ab Initio半从头算法Slater X半经验法Semi-emperical独立电子对IEPA第一原理计算常

13、用的量子化学计算方法量子力学理论Born-Oppenheimer近似非相对论近似单电子近似Hartree-Fock方方 程程Roothaan方方 程程自洽场从头算SCF-ab initio密度泛函法DFT超 HFLCMTO-X耦合电子对CEPA组 态相互作用CI微扰处理MP多组态自洽场MCSCF价电子从头算EP(VP)模拟从头算SAMO分子碎片法MF梯度近似GGA浮动球高斯法FSGOAM1C-EHMOEHMOIT-EHMOMCNDOCNDOMINDOINDOMNDONDDOPM3MSW-XDV-XLCAO-X局域密度近 似LDA从头算法Ab Initio半从头算法Slater X半经验法Sem

14、i-emperical独立电子对IEPA基组的选择 Route Section采用自由格式，大小写不敏感采用自由格式，大小写不敏感 同一行不同项之间可以使用空格，逗号和同一行不同项之间可以使用空格，逗号和“/”连接；连接； 例：例：#p hf/6-31g scfcyc=230 scfcon=8 #p,hf/6-31g,scfcyc=230,scfcon=8 关键词可以通过关键词可以通过 keyword =option,keyword(option), keyword=(option1, option2, .)，keyword( option1, option2, .)指定指定 例：例：#p H

15、F/STO-3G opt #p HF/STO-3G opt=(TS,readfc) Title部分必须输入，但是程序并不执行，起标识和说明作用部分必须输入，但是程序并不执行，起标识和说明作用Title SectionCharge & Multipl. 输入分子的电荷和自旋多重度输入分子的电荷和自旋多重度 例：例： 电荷电荷 多重度多重度2s+1 H2O 0 1 H3O+ 1 1 NO 0 2 电荷多重度部分通常也算作分子说明电荷多重度部分通常也算作分子说明Molecular Specification 分子说明部分主要用来定义原子核相对位置分子说明部分主要用来定义原子核相对位置原子核相对位置可

16、以用笛卡尔坐标，内坐标原子核相对位置可以用笛卡尔坐标，内坐标(Z-matrix)，或者是二者混合，或者是二者混合 表示表示笛卡尔坐标是内坐标的一种特殊形式笛卡尔坐标是内坐标的一种特殊形式分子坐标的格式为分子坐标的格式为: 元素符号, x, y, z 元素符号 (n) 原子1 键长 原子2 键角 原子3 二面角 (0,180) 以原子在分子中的序数表示例例1：使用使用HF方法，优化方法，优化 H2O2分子分子%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1 HO 1 0.9O 2 1.4 1 105.0H 3

17、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 内坐标表示内坐标表示 笛卡尔直角坐标表示笛卡尔直角坐标表示上述两个作业将对上述两个作业将对H2O2分子的结构进行完全优化，包括所有的键长键角和二面角。分子的结构进行完全优化，包括所有的键长键角和二面角。通过将内坐标定义成变量，可以对分子结构进行部分优化通过将内坐标定

18、义成变量，可以对分子结构进行部分优化%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 a1H 3 r1 2

19、 a1 1 d1r1=0.9a1=105.0 变量部分变量部分 d1=120.0r2=1.4 常量部分常量部分通过设定常量，可以使程序通过设定常量，可以使程序只对分子结构进行只对分子结构进行部分部分优化优化例例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和和0的的常量常量

20、来控制分子的平面构型来控制分子的平面构型 例例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运行界面GaussianOverlay1O

21、verlay9,10,11,99L101L102L122Overlay0L0L001 Gaussian 程序的结构程序的结构link0: 初始化程序，控制初始化程序，控制占位占位link1: 读入并处理读入并处理计算执行路径，计算执行路径，建建立要执行的立要执行的link列表列表link9999: 终止计算终止计算Overlay99L9999Gaussian03 程序界面程序界面Preferences:对对Gaussian程序进行初程序进行初始化设置始化设置自定义外部文字编辑器，用来打开自定义外部文字编辑器，用来打开.out文文件件link.exe所在的文件夹所在的文件夹临时文件存放文件夹临时

22、文件存放文件夹Gaussian03 图形工具栏图形工具栏开始作业开始作业 暂停当前作业暂停当前作业 当前当前link后暂停后暂停 终止当前作业和批处理终止当前作业和批处理 恢复当前作业恢复当前作业 在当前作业完成后终止批处理在当前作业完成后终止批处理 终止当前作业终止当前作业 编辑或建立批处理编辑或建立批处理 批处理专用批处理专用打开外部编辑器打开外部编辑器 编辑输出文件编辑输出文件-Link1- 指定多步任务Gaussian作业输出示例作业输出示例 （H2O SP） Entering Link 1 = C:G03Wl1.exe PID= 1572. Copyright (c) 1988,19

23、90,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, Inc.),

24、 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 federally reg

25、istered 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. 历历史史和和版版权权Cite this work as: Gaussian

26、 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. Nakatsuji, M.

27、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. Ayala, K. Mo

28、rokuma, 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. Piskorz, I. K

29、omaromi, 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=h2o.chk

30、 %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 重重新新打打印印作作业业输输入入和和标标准准途途径径作业执行过程作业执行过程1/38=1/

31、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 读取标题和分子输入部分； L202 坐标重新定向，计算对称性，检查变量；L301 产生基组信息； L302 计算重叠积分，动能积分和势能积分；L303 计算多极积分； L401 形成分子轨道初猜； L502 自洽场迭代求解；L601 布居数和相关分析；L9999 完成计算和输

32、出L101 读取标题和分子输入部分读取标题和分子输入部分(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 坐标重新定向，计算对称性，检查变量坐标重新定向，计算对称性，检查变量(Enter C:G03Wl202.exe) Input orientati

33、on: - 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 Stoichiometry H2O F

34、ramework 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.000000 0.

35、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 产生基组信息产生基组信息(Enter C:G03Wl301.exe) Standard basis: 6-31G (6D, 7F) There are 7 symmetry adapted basis functio

36、ns 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 integral sy

37、mmetry 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 IRanWt= -1 I

38、RanGd= 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 NCelE2= 1 N

39、ClLst= 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. Leave Lin

40、k 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 IRadAn=

41、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, MaxMem=

42、 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.00D-05.

43、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= 0 MaxBra

44、= 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-05 EMaxC

45、= 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-04 SCF D

46、one: 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 布居数和相关分析布居数和相关分析(Enter C:G03Wl601.exe) Copying SCF densities t

47、o 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 Alpha v

48、irt. 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.357696

49、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= 0.0000 Y

50、= 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