分子平衡几何构型优化及分子性质计算.docxVIP

实验目的1、学会从实际操作出发，掌握程序的使用，得到正确的数据。2、学习Gaussian程序使用，运用程序进行几种分子模型的构建及优化。3、利用Gaussian程序对分子体系薛定谔方程所代表的化学理论加深理解。 4、掌握构建分子模型的方法和分子几何构型的输入方法，为目标分子设定计算坐标，模拟化学分子，能够正确解读计算结果，采集有用的结果数据。二、实验原理1、Gaussian程序可以作为功能强大的工具，用于研究许多化学领域的课题，例如取代基的影响，化学反应机理，势能曲面和激发能等等。它还可以预言分子和化学反应的许多性质，如，分子能量和结构、过渡态的能量和结构、电子密度分布、热力学性质、振动频率、红外和拉曼光谱、NMR化学位移、极化率和静电势，等等。本实验教材的重点，通过驻点(分子和反应势能面上的极小点和鞍点)的优化和性质计算，进行结构与性质关系的预测和化学反应动力学，包括反应速率和反应机理的预测。2、HF方程自洽场近似：在结构化学中，“变数分离”方法对于单电子体系(氢原子和类氢离子)的Schr?dinger方程进行精确求解。但是对于多电子的分子体系，由于第i个电子与其余电子间的排斥能取决于所有电子的坐标，使这种分离变为不可能。但可以在定核近似下将核的运动分离出去后，在固定的核势场中近似求解多电子体系的能量本征方程。 具体做法是，对第i个电子，可以假定一个单电子的分子轨道(单电子近似)，并将它用现成的原子轨道线性展开(LCAO近似)。这时，Schr?dinger方程由微分方程变成一个齐次线性的代数方程组。求解该方程组，即求各分子轨道能级及相应的分子轨道展开系数。具体过程是在给定的核坐标下，先猜测一组展开系数(极端情况均为0)，代入方程组得到一组新的系数，再代入方程组求解，周而复始，直到前后两组系数相同，称为“自恰场迭代”。这就是HF自恰场分子轨道方法。3、基组用于描述体系波函数的若干具有一定性质的函数。基组是量子化学从头计算的基础，在量子化学中有着非常重要的意义。A: Gauss 函数中心l: 角量子数三、仪器设备与软件计算机、GaussView程序、Gaussian 程序四、实验数据处理（一）、在水分子中(NN)、分子总能(E)以及电子动能(KE)、电子与核的吸引能(PE)、电子排斥能(EE)和FMO能 (即HOMO与LUMO能), 核对各能量之间的关系；分子中的各个原子上的静电荷以及化学键的键电荷，分子的偶极矩。(1)核排斥能(NN)9.2065538966a.u.分子总能(E)-76.0232281985a.u.电子动能(KE)75.81500790115a.u.电子与核的吸引能(PE)-198.9309839994a.u.电子排斥能(EE)37.88619400319a.u.HOMO能-0.49716a.u.LUMO能0.21238a.u.(2)核对各能量之间的关系分子总能分别为-20.56036， -1.34091， -0.70404， -0.56899 ，-0.49716得=-76程序运算出的结果-76.0232281985比较相近，误差计算得0.0008788%。可知该计算方法可行。（3）分子中的各个原子上的静电荷1 O -0.6732132 H 0.3366063 H 0.336606（4）化学键的键电荷=0.308856e×2=0.617712e（5）分子的偶极矩=2.1842(二）、请记录第二题至第七题中第一次和最后一次的分子总能，偶极矩，原子电荷，分子几何构型。（2）H2O分子的几何构型优化1O2H 1 r13H 1r2 2 a1(a.u.)各原子电荷(e)点群分子几何构型偶极矩优化前-76.02254165671O -0.6763172H 0.3376703H 0.338647Csr1=0.960r2=0.970a1=104.52.1911优化后-76.0236151O -0.6706862H 0.3353423H 0.335344Csr1=0.943043r2=0.943061a1=105.9682.1476(3)H2O分子结构的优化1O 2H 1 r1 3H 1 r 2 a1(a.u.)各原子电荷(e)点群键长r(?)键角偶极矩优化前-76.02298508281O -0.6745822H 0.3372912H 0.337291C2vr1=0.960a2=104.52.1870优化后-76.0236151O -0.6706862H 0.3353433H 0.335343C2vr1=0.943053a1=105.9682 2.1476(4)乙烯分子结构的优化1C2