Gaussian程序使用：分子平衡几何构型优化及分子性质计算.doc

Gaussian程序使用：分子平衡几何构型优化及分子性质计算【实验目的与要求】计算化学，其本质是对分子体系薛定谔方程所代表的化学理论通过大型计算机程序的求解，模拟化学的各种实验研究。作为一门计算化学实验课程，主要目的是从实际操作出发，掌握程序的使用，以便得到预期的结果。对于所涉及的理论和方法，只要求结合程序的演算能够定性予以理解。本实验主要涉及优化分子几何构型的程序输入及计算结果的解读。【实验原理】1. 在结构化学中，曾用“变数分离”方法对于单电子体系(氢原子和类氢离子)的Schrdinger方程进行精确求解。但是对于多电子的分子体系，由于第i个电子与其余电子间的排斥能取决于所有电子的坐标，使这种分离变为不可能。但可以在定核近似下将核的运动分离出去后，在固定的核势场中近似求解多电子体系的能量本征方程。具体做法是，对第i个电子，可以假定一个单电子的分子轨道(单电子近似)，并将它用现成的原子轨道线性展开(LCAO近似)。这时，Schrdinger方程由微分方程变成一个齐次线性的代数方程组。求解该方程组，即求各分子轨道能级及相应的分子轨道展开系数。具体过程是在给定的核坐标下，先猜测一组展开系数(极端情况均为0)，代入方程组得到一组新的系数，再代入方程组求解，周而复始，直到前后两组系数相同，称为“自恰场迭代”。这就是HF自恰场分子轨道方法。2. 量子化学中的基组是用于描述体系波函数的若干具有一定性质的函数。基组是量子化学从头计算的基础，在量子化学中有着非常重要的意义。在量子化学计算中，根据体系的不同，需要选择不同的基组，构成基组的函数越多，基组便越大，计算的精度也越高，计算量也随之增大。要提高量子化学计算精度，必须加大基组的规模，即增加基组中基函数的数量，增大基组规模的一个方法是劈裂原子轨道，也就是使用多个基函数来表示一个原子轨道。劈裂价键基组就是应用上述方法构造的较大型基组，所谓劈裂价键就是将价层电子的原子轨道用两个或以上基函数来表示。常见的劈裂价键基组有3-21G、4-21G、4-31G、6-31G、6-311G等，在这些表示中前一个数字用来表示构成内层电子原子轨道的高斯型函数数目，“-”以后的数字表示构成价层电子原子轨道的高斯型函数数目。如6-31G所代表的基组，每个内层电子轨道是由6个高斯型函数线性组合而成，每个价层电子轨道则会被劈裂成两个基函数，分别由3个和1个高斯型函数线性组合而成。劈裂价键基组能够比STO-NG基组更好地描述体系波函数，同时计算量也比最小基组有显著的上升需要根据研究的体系不同而选择相应的基组进行计算。【实验内容】1. 输入的基本组成部分。如图所示：命令部分：作业类型部分总是必需的，它以第一行的第一列用符号“#”开始，此行的其余部分是自由格式。#p HF/STO-3G为不优化仅仅用此方法计算得出结果。#p HF/STO-3G opt=z-matrix表示对所输入变量进行循环，找到满足最大力、均方根力、最大位移、均方根位移的结果。标题输入部分：这一部分也是必需的，但程序并不执行，它出现在最后结果中，起标示和说明作用，应该是作业的特征信息。该部分也以一个空行结束。分子输入：这一部分第一行是分子的电荷和自旋多重度。例如，闭壳层中性分子输入“0 1”，带电的负离子自由基为“-1 2”。其余的行说明原子核的相对位置，用内坐标输入最方便。所谓内坐标，是用键长、键角和二面角定义原子核的位置。2按以下数据输入关于水分子的计算：（1）H2O#p HF/6-31G(d,p)Water energy标题输入部分0 1分子输入O H 1 0.956H 1 0.956 2 104.5这是H2O 的HF/STO-3G 基态的单点计算(固定核坐标)，未引入变量。若优化这些参数则（2）#p HF/6-31G(d,p) opt=z-matrixWater 6-31G(d,p) structure 0 1OH 1 r1H 1 r2 2 a1r1 = 0.96r2 = 0.97 a1 = 104.5变量说明：这是一个优化过程，优化键长r1、r2和键角a1，并极小化HF/6-31G能量。（3）如果分子说明写成：0 1OH 1 rH 1 r 2 ar=0.96a=104.5两个键长为同一变量，意味优化受到限制。通过变量限制可控制分子对称性(此例为C2v，上例为Cs)。（4) 用相同方法(#p HF/6-31G(d,p) opt=z-matrix)进行乙烯分子C2V构型优化。#p HF/6-31G(d,p) opt=z-matrixC2H4 opt0 1C C 1 r1H 1 r2 2 a1 H 1 r2 2 a1 3 180.0 H 2 r2 1 a1 3 0.0 H 2 r2 1 a1 4 0.0 r1=1.32r2=1.09 a1=120.0（5）用相同方法(#p HF/6-31G(d,p) opt=z-matrix)进行CH3F 分子C3v构型优化。#p HF/6-31G(d,p) opt=z-matrixCH3F (C3v) opt0 1C F 1 r1H 1 r2 2 aH 1 r2 2 a 3 bH 1 r2 2 a 3 -br1=1.38r2=1.09a=110.6b=120.0（5）用相同方法(#p HF/6-31G(d,p) opt=z-matrix)独立计算C2H2分子的结构。#p HF/6-31G(d,p) opt=z-matrixC2H2 opt0 1C C 1 r1X 2 1.0 1 90.0X 1 1.0 2 90.0 3 0.0H 2 r2 3 90.0 1 180.0H 1 r2 4 90.0 2 180.0r1=1.26r2=1.09（6）用相同方法(#p HF/6-31G(d,p) opt=z-matrix)独立计算呋喃分子C4H4O的结构。#p HF/6-31G(d,p) opt=z-matrixC4H4O opt0 1O C 1 r1C 1 r1 2 a1C 2 r2 1 a2 3 0.0C 3 r2 1 a2 2 0.0H 2 r3 1 a3 3 180.0H 3 r3 1 a3 2 180.0H 4 r4 2 a4 1 180.0H 5 r4 3 a4 1 180.0r1=1.40r2=1.35r3=1.08r4=1.09a1=107.0a2=108.0a3=120.0a4=127.0【实验结果】（1）H2O分子的能量计算 (NN)、分子总能(E)以及电子动能(KE)、电子与核的吸引能(PE)、电子排斥能(EE)和FMO能 (即HOMO与LUMO能), 核对各能量之间的关系；分子中的各个原子上的静电荷以及化学键的键电荷，分子的偶极矩。NN (nuclear repulsion energy) =9.20655 HartreesE(RHF) =-76.02323 a.u.KE= 75.81501 a.u.PE=-198.93098 a.u.EE= 37.88619 a.u.HOMO能=-0.49716 a.u. LUMO能=0.21238 a.u.关系：E=NN+KE+PE+EE 各原子上的静电荷(a.u.)：1 O: -0.67321 2 H：0.33661 3 H：0.33661化学键的键电荷(a.u.)：1O与2H：0.309 1O与3H： 0.309分子的偶极矩：2.1842 D（2-7），记录第一次和最后一次的分子总能，偶极矩，原子电荷，分子几何构型。（2）H2O分子的几何构型优化第一次：E=-76.02254 a.u.，偶极矩为2.1911 D，原子电荷(a.u.)分别为1O -0.67632，2H 0.33767，3H 0.33865，分子所属点群为Cs，几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 O 2 2 H 1 0.942758( 1) 3 3 H 1 0.941795( 2) 2 105.241( 3) -最后一次：E =-76.02362 a.u.，偶极矩为2.1476 D，原子电荷(a.u.)分别为1O -0.67069，2H 0.33534，3H 0.33534，分子所属点群为Cs,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 O 2 2 H 1 0.943043( 1) 3 3 H 1 0.943061( 2) 2 105.968( 3) -（3）H2O分子结构的优化第一次：E=-76.02299 a.u.，偶极矩为2.1870 D，原子电荷(a.u.)分别为1O -0.67458，2H 0.33729，3H 0.33729，分子所属点群为C2V,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 O 2 2 H 1 0.942614( 1) 3 3 H 1 0.942614( 2) 2 105.446( 3) -最后一次：E =-76.02362 a.u.，偶极矩为2.1476 D，原子电荷(a.u.)分别为1O -0.67069，2H 0.33534，3H 0.33534，分子所属点群为C2V,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 O 2 2 H 1 0.943053( 1) 3 3 H 1 0.943053( 2) 2 105.968( 3) -(4) 乙烯分子C2V结构的优化第一次：E=-78.03756 a.u.，偶极矩为2.1870 D，原子电荷(a.u.)分别为1C -0.25602，2C -0.25602，3H 0.12801，4H 0.12801, 5H 0.12801, 6H 0.12801，分子所属点群为D2h,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 C 2 2 C 1 1.318788( 1) 3 3 H 1 1.074226( 2) 2 121.425( 6) 4 4 H 1 1.074226( 3) 2 121.425( 7) 3 180.000( 10) 0 5 5 H 2 1.074226( 4) 1 121.425( 8) 3 0.000( 11) 0 6 6 H 2 1.074226( 5) 1 121.425( 9) 4 0.000( 12) 0 -最后一次：E=-78.03884 a.u.，偶极矩为0.0000，原子电荷(a.u.)分别为1C -0.25437，2C -0.25437，3H 0.12718，4H 0.12718, 5H 0.12718, 6H 0.12718，分子所属点群为D2h,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 C 2 2 C 1 1.316479( 1) 3 3 H 1 1.076391( 2) 2 121.744( 6) 4 4 H 1 1.076391( 3) 2 121.744( 7) 3 180.000( 10) 0 5 5 H 2 1.076391( 4) 1 121.744( 8) 3 0.000( 11) 0 6 6 H 2 1.076391( 5) 1 121.744( 9) 4 0.000( 12) 0 -（5）CH3F (C3v) opt第一次：E=-139.03895 a.u.，偶极矩为2.0454 D，原子电荷(a.u.)分别为1C 0.08377，2F -0.40853，3H 0.10825，4H 0.10825, 5H 0.10825,分子所属点群为C3v,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 C 2 2 F 1 1.360832( 1) 3 3 H 1 1.082144( 2) 2 108.721( 5) 4 4 H 1 1.082144( 3) 2 108.721( 6) 3 120.000( 8) 0 5 5 H 1 1.082144( 4) 2 108.721( 7) 3 -120.000( 9) 0 -最后一次：E=-139.03974 a.u.，偶极矩为1.9819 D，原子电荷(a.u.)分别为1C 0.08660，2F -0.40719，3H 0.10686，4H 0.10686, 5H 0.10686,分子所属点群为C3v,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 C 2 2 F 1 1.364820( 1) 3 3 H 1 1.082638( 2) 2 109.185( 5) 4 4 H 1 1.082638( 3) 2 109.185( 6) 3 120.000( 8) 0 5 5 H 1 1.082638( 4) 2 109.185( 7) 3 -120.000( 9) 0 -(6) C2H2分子的结构第一次：E=-76.80898 a.u.，偶极矩为0.0000，原子电荷(a.u.)分别为1C -0.24727，2C -0.24727，3H 0.24727，4H 0.24727,分子所属点群为Dh,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 C 2 2 C 1 1.163942( 1) 3 X 2 1.000000( 2) 1 90.000( 6) 4 X 1 1.000000( 3) 2 90.000( 7) 3 0.000( 10) 0 5 3 H 2 1.053127( 4) 3 90.000( 8) 1 180.000( 11) 0 6 4 H 1 1.053127( 5) 4 90.000( 9) 2 180.000( 12) 0 -最后一次：E=-76.82184 a.u.，偶极矩为0.0000，原子电荷(a.u.)分别为1C -0.23345，2C -0.23345，3H 0.23345，4H 0.23345,分子所属点群为Dh,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 C 2 2 C 1 1.186066( 1) 3 X 2 1.000000( 2) 1 90.000( 6) 4 X 1 1.000000( 3) 2 90.000( 7) 3 0.000( 10) 0 5 3 H 2 1.056877( 4) 3 90.000( 8) 1 180.000( 11) 0 6 4 H 1 1.056877( 5) 4 90.000( 9) 2 180.000( 12) 0 -（7）呋喃C4H4O分子的结构第一次：E=-228.62293 a.u.，偶极矩为1.0915 D，原子电荷(a.u.)分别为1O -0.59964，2C 0.17372，3C 0.17372，4C -0.20870，5C -0.20870，6H 0.17647，7H 0.17647，8H 0.15833，9H 0.15833，分子所属点群为C2v,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 O 2 2 C 1 1.349886( 1) 3 3 C 1 1.349886( 2) 2 105.967( 9) 4 4 C 2 1.337001( 3) 1 111.267( 10) 3 0.000( 16) 0 5 5 C 3 1.337001( 4) 1 111.267( 11) 2 0.000( 17) 0 6 6 H 2 1.065250( 5) 1 115.598( 12) 3 180.000( 18) 0 7 7 H 3 1.065250( 6) 1 115.598( 13) 2 180.000( 19) 0 8 8 H 4 1.070281( 7) 2 126.956( 14) 1 180.000( 20) 0 9 9 H 5 1.070281( 8) 3 126.956( 15) 1 180.000( 21) 0 -最后一次：E=-228.63294 a.u.，偶极矩为0.7714 D，原子电荷(a.u.)分别为1O -0.54337，2C 0.16901，3C 0.16901，4C -0.22309，5C -0.22309，6H 0.16617，7H 0.16617，8H 0.15960，9H 0.15960，分子所属点群为C2v,几何构型如下表所示：CD Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J - 1 1 O 2 2 C 1 1.343671( 1) 3 3 C 1 1.343671( 2) 2 107.146( 9) 4 4 C 2 1.339105( 3) 1 110.794( 10) 3 0.000( 16) 0 5 5 C 3 1.339105( 4) 1 110.794( 11) 2 0.000( 17) 0 6 6 H 2 1.068418( 5) 1 116.292( 12) 3 180.000( 18) 0 7 7 H 3 1.068418( 6) 1 116.292( 13) 2 180.000( 19) 0 8 8 H 4 1.070276( 7) 2 126.704( 14) 1 180.000( 20) 0 9 9 H 5 1.070276( 8) 3 126.704( 15) 1 180.000( 21