计算机在化工中的应用-上机实验五-Gaussian.docVIP

计算机化学-Gaussian实验 练习1利用Gaussian03进行甲醛单点能，并查看分子轨道。 输入文件为： %chk=HCHO.chk # RHF/6-31G(d) Pop=Full Test Formaldehyde Single Point 0,1 C 0 0.0 0.0 O 0.0 1.22 0.0 H 0.94 -0.54 0.0 H -0.94 -0.54 0.0 采用GaussView打开Scratch目录中的HCHO.chk文件，打开菜单的“Results”——?“Surfaces”，逐个查看各轨道的形状。 比较：不同的输入(优化结构) # B3LYP/6-31G OPT Formaldehyde 0 1 C1 O2 1 A H3 1 B 2 C H4 1 B 2 C 3 D A = 1.28 B = 1.10 C = 121.0 D = 180.0 比较：不同的输入(振动分析) # B3LYP/6-31G Freq Example: aldehyde 0 1 C1 O2 1 A H3 1 B 2 C H4 1 B 2 C 3 D A = 1.28 B = 1.10 C = 121.0 D = 180.0 比较：甲醛的激发态优化 在Gaussian 中可以进行激发态的优化和频率分析。步骤是，首先进行能量计算得到激发态， 然后在此结构基础上进行结构优化和频率分析。 %Chk=c2ho_es #T RCIS/6-31+G(D) Test Formaldehyde (C2V) Vertical Excited States 0,1 C O, 1, AB H,1, AH, 2, HAB H,1, AH, 2, HAB, 3, 180. AB=1.18429 AH=1.09169 HAB=122.13658 --Link1-- %Chk=c2ho_es %NoSave # RCIS(Root=1,Read)/6-31+G(D) Guess=Read Opt Freq Test Geom=AllCheck 第一个步骤计算最低的三个激发态，第二步骤是利用其结果进行结构优化和频率分析. 在本例的计算中，得到的激发态出现了虚频，其结构显示碳原子要离开原子平面。 由于原来的分子结构中，所有原子都在同一平面上，所以，也得到一个共平面的激发态结构。 这样，就要设法产生一个正确的几何结构。(激发态的结构优化目前没有包含内坐标冗余， 所以在初始结构中定义一个平面结构，得到的激发态结构就也是共平面的)。 比较： 甲醛激发态优化 ，确定甲醛第一激发态的结构，比较红外光谱。 #T RCIS(Root=1)/6-31+G(D) Opt Freq Density=Current Test Formaldehyde n--pi* Cs Optimization 0,1 C O,1,R2 X,1,1.,2,A3 H,1,R4,3,A4,2,90. H,1,R4,3,A4,2,-90. R2=1R4=1A3=145. A4=60. 这里使用了虚原子，使得对于分子的描述更加简单，这个分子也不再是共平面的结构。计算采用校正因子0.8929，结果如下： 计算值 校正值 实验值 平面外弯曲 495 442 683 CH2 摇摆 978 873 898 CH2 剪式 1426 1273 1290 CO 伸缩 1647 1471 1173 对称CH伸缩 3200 2857 2847 反对称CH伸缩 3295 2942 2968 基本与实验值相符 练习2 利用Gaussian03进行水分子结构优化与红外分析 输入文件为： %chk=H2O #P B3LYP/6-311G OPT FREQ Geometry optimization and the frequency calculation of H2O 0,1 O H 1 10 H 1 1.0 2 105.4 练习3 利用GaussView画出苯环分子，并计算单点能。 练习4 1,2-二氯-1,2-二氟乙烷分子几何构型输入法及气态能量计算 实验目的： 掌握分子几何构型的输入方法。 掌握能量计算的方法及对其结果z进行分析。 基本原理： ????? 单点能计算是指对给定几何构型的分子的能量以及性质进行计算,由于分子的几何构型是固定不变的,只是一个点,所以叫单点能计算。单点能的计算可以在不同理论等级,采用不同基组进行。使用单