Gaussian程序使用分子平衡几何构型优化及分子性质计算.docx
《Gaussian程序使用分子平衡几何构型优化及分子性质计算.docx》由会员分享,可在线阅读,更多相关《Gaussian程序使用分子平衡几何构型优化及分子性质计算.docx(15页珍藏版)》请在冰点文库上搜索。
Gaussian程序使用分子平衡几何构型优化及分子性质计算
aussian程序使用:
分子平衡几何构型优化及分子性质计算
Gaussian程序使用:
分子平衡几何构型优化及分子性质计算
【实验目的与要求】
计算化学,其本质是对分子体系薛定谔方程所代表的化学理论通过大型计算机程序的求解,模拟化学的各种实验研究。
作为一门计算化学实验课程,主要目的是从实际操作出发,掌握程序的使用,以便得到预期的结果。
对于所涉及的理论和方法,只要求结合程序的演算能够定性予以理解。
本实验主要涉及优化分子几何构型的程序输入及计算结果的解读。
【实验原理】
1.在结构化学中,曾用“变数分离”方法对于单电子体系(氢原子和类氢离子)的Schrödinger方程进行精确求解。
但是对于多电子的分子体系,由于第i个电子与其余电子间的排斥能取决于所有电子的坐标,使这种分离变为不可能。
但可以在定核近似下将核的运动分离出去后,在固定的核势场中近似求解多电子体系的能量本征方程。
具体做法是,对第i个电子,可以假定一个单电子的分子轨道(单电子近似),并将它用现成的原子轨道线性展开(LCAO近似)。
这时,Schrödinger方程由微分方程变成一个齐次线性的代数方程组。
求解该方程组,即求各分子轨道能级及相应的分子轨道展开系数。
具体过程是在给定的核坐标下,先猜测一组展开系数(极端情况均为0),代入方程组得到一组新的系数,再代入方程组求解,周而复始,直到前后两组系数相同,称为“自恰场迭代”。
这就是HF自恰场分子轨道方法。
2.量子化学中的基组是用于描述体系波函数的若干具有一定性质的函数。
基组是量子化学从头计算的基础,在量子化学中有着非常重要的意义。
在量子化学计算中,根据体系的不同,需要选择不同的基组,构成基组的函数越多,基组便越大,计算的精度也越高,计算量也随之增大。
要提高量子化学计算精度,必须加大基组的规模,即增加基组中基函数的数量,增大基组规模的一个方法是劈裂原子轨道,也就是使用多个基函数来表示一个原子轨道。
劈裂价键基组就是应用上述方法构造的较大型基组,所谓劈裂价键就是将价层电子的原子轨道用两个或以上基函数来表示。
常见的劈裂价键基组有3-21G、4-21G、4-31G、6-31G、6-311G等,在这些表示中前一个数字用来表示构成内层电子原子轨道的高斯型函数数目,“-”以后的数字表示构成价层电子原子轨道的高斯型函数数目。
如6-31G所代表的基组,每个内层电子轨道是由6个高斯型函数线性组合而成,每个价层电子轨道则会被劈裂成两个基函数,分别由3个和1个高斯型函数线性组合而成。
劈裂价键基组能够比STO-NG基组更好地描述体系波函数,同时计算量也比最小基组有显著的上升需要根据研究的体系不同而选择相应的基组进行计算。
【实验内容】
1.输入的基本组成部分。
如图所示:
命令部分:
作业类型部分总是必需的,它以第一行的第一列用符号“#”开始,此行的其余部分是自由格式。
#pHF/STO-3G为不优化仅仅用此方法计算得出结果。
#pHF/STO-3Gopt=z-matrix表示对所输入变量进行循环,找到满足最大力、均方根力、最大位移、均方根位移的结果。
标题输入部分:
这一部分也是必需的,但程序并不执行,它出现在最后结果中,起标示和说明作用,应该是作业的特征信息。
该部分也以一个空行结束。
分子输入:
这一部分第一行是分子的电荷和自旋多重度。
例如,闭壳层中性分子输入“01”,带电的负离子自由基为“-12”。
其余的行说明原子核的相对位置,用内坐标输入最方便。
所谓内坐标,是用键长、键角和二面角定义原子核的位置。
2.按以下数据输入关于水分子的计算:
(1)H2O
#pHF/6-31G(d,p)
Waterenergy••••••••••••••••••••••标题输入部分
01•••••••••••••••••••••••••••••••••分子输入
O
H10.956
H10.9562104.5
这是H2O的HF/STO-3G基态的单点计算(固定核坐标),未引入变量。
若优化这些参数则
(2)
#pHF/6-31G(d,p)opt=z-matrix
Water6-31G(d,p)structure
01
O
H1r1
H1r22a1
r1=0.96
r2=0.97
a1=104.5
变量说明:
这是一个优化过程,优化键长r1、r2和键角a1,并极小化HF/6-31G能量。
(3)如果分子说明写成:
01
O
H1r
H1r2a
r=0.96
a=104.5
两个键长为同一变量,意味优化受到限制。
通过变量限制可控制分子对称性(此例为C2v,上例为Cs)。
(4)用相同方法(#pHF/6-31G(d,p)opt=z-matrix)进行乙烯分子C2V构型优化。
#pHF/6-31G(d,p)opt=z-matrix
C2H4opt
01
C
C1r1
H1r22a1
H1r22a13180.0
H2r21a130.0
H2r21a140.0
r1=1.32
r2=1.09
a1=120.0
(5)用相同方法(#pHF/6-31G(d,p)opt=z-matrix)进行CH3F分子C3v构型优化。
#pHF/6-31G(d,p)opt=z-matrix
CH3F(C3v)opt
01
C
F1r1
H1r22a
H1r22a3b
H1r22a3-b
r1=1.38
r2=1.09
a=110.6
b=120.0
(5)用相同方法(#pHF/6-31G(d,p)opt=z-matrix))独立计算C2H2分子的结构。
#pHF/6-31G(d,p)opt=z-matrix
C2H2opt
01
C
C1r1
X21.0190.0
X11.0290.030.0
H2r2390.01180.0
H1r2490.02180.0
r1=1.26
r2=1.09
(6)用相同方法(#pHF/6-31G(d,p)opt=z-matrix))独立计算呋喃分子C4H4O的结构。
#pHF/6-31G(d,p)opt=z-matrix
C4H4Oopt
01
O
C1r1
C1r12a1
C2r21a230.0
C3r21a220.0
H2r31a33180.0
H3r31a32180.0
H4r42a41180.0
H5r43a41180.0
r1=1.40
r2=1.35
r3=1.08
r4=1.09
a1=107.0
a2=108.0
a3=120.0
a4=127.0
【实验结果】
(1)H2O分子的能量计算
(NN)、分子总能(E)以及电子动能(KE)、电子与核的吸引能(PE)、电子排斥能(EE)和FMO能(即HOMO与LUMO能),核对各能量之间的关系;分子中的各个原子上的静电荷以及化学键的键电荷,分子的偶极矩。
NN(nuclearrepulsionenergy)=9.20655Hartrees
E(RHF)=-76.02323a.u.
KE=75.81501a.u.
PE=-198.93098a.u.
EE=37.88619a.u.
HOMO能=-0.49716a.u.LUMO能=0.21238a.u.
关系:
E=NN+KE+PE+EE
各原子上的静电荷(a.u.):
1O:
-0.673212H:
0.336613H:
0.33661
化学键的键电荷(a.u.):
1O与2H:
0.3091O与3H:
0.309
分子的偶极矩:
2.1842D
(2-7),记录第一次和最后一次的分子总能,偶极矩,原子电荷,分子几何构型。
(2)H2O分子的几何构型优化
第一次:
E=-76.02254a.u.,偶极矩为2.1911D,原子电荷(a.u.)分别为1O-0.67632,2H0.33767,3H0.33865,分子所属点群为Cs,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11O
22H10.942758
(1)
33H10.941795
(2)2105.241(3)
------------------------------------------------------------------------
最后一次:
E=-76.02362a.u.,偶极矩为2.1476D,原子电荷(a.u.)分别为1O-0.67069,2H0.33534,3H0.33534,分子所属点群为Cs,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11O
22H10.943043
(1)
33H10.943061
(2)2105.968(3)
------------------------------------------------------------------------
(3)H2O分子结构的优化
第一次:
E=-76.02299a.u.,偶极矩为2.1870D,原子电荷(a.u.)分别为1O-0.67458,2H0.33729,3H0.33729,分子所属点群为C2V,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11O
22H10.942614
(1)
33H10.942614
(2)2105.446(3)
------------------------------------------------------------------------
最后一次:
E=-76.02362a.u.,偶极矩为2.1476D,原子电荷(a.u.)分别为1O-0.67069,2H0.33534,3H0.33534,分子所属点群为C2V,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11O
22H10.943053
(1)
33H10.943053
(2)2105.968(3)
------------------------------------------------------------------------
(4)乙烯分子C2V结构的优化
第一次:
E=-78.03756a.u.,偶极矩为2.1870D,原子电荷(a.u.)分别为1C-0.25602,2C-0.25602,3H0.12801,4H0.12801,5H0.12801,6H0.12801,分子所属点群为D2h,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11C
22C11.318788
(1)
33H11.074226
(2)2121.425(6)
44H11.074226(3)2121.425(7)3180.000(10)0
55H21.074226(4)1121.425(8)30.000(11)0
66H21.074226(5)1121.425(9)40.000(12)0
------------------------------------------------------------------------
最后一次:
E=-78.03884a.u.,偶极矩为0.0000,原子电荷(a.u.)分别为1C-0.25437,2C-0.25437,3H0.12718,4H0.12718,5H0.12718,6H0.12718,分子所属点群为D2h,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11C
22C11.316479
(1)
33H11.076391
(2)2121.744(6)
44H11.076391(3)2121.744(7)3180.000(10)0
55H21.076391(4)1121.744(8)30.000(11)0
66H21.076391(5)1121.744(9)40.000(12)0
------------------------------------------------------------------------
(5)CH3F(C3v)opt
第一次:
E=-139.03895a.u.,偶极矩为2.0454D,原子电荷(a.u.)分别为1C0.08377,2F-0.40853,3H0.10825,4H0.10825,5H0.10825,分子所属点群为C3v,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11C
22F11.360832
(1)
33H11.082144
(2)2108.721(5)
44H11.082144(3)2108.721(6)3120.000(8)0
55H11.082144(4)2108.721(7)3-120.000(9)0
------------------------------------------------------------------------
最后一次:
E=-139.03974a.u.,偶极矩为1.9819D,原子电荷(a.u.)分别为1C0.08660,2F-0.40719,3H0.10686,4H0.10686,5H0.10686,分子所属点群为C3v,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11C
22F11.364820
(1)
33H11.082638
(2)2109.185(5)
44H11.082638(3)2109.185(6)3120.000(8)0
55H11.082638(4)2109.185(7)3-120.000(9)0
------------------------------------------------------------------------
(6)C2H2分子的结构
第一次:
E=-76.80898a.u.,偶极矩为0.0000,原子电荷(a.u.)分别为1C-0.24727,2C-0.24727,3H0.24727,4H0.24727,分子所属点群为D∞h,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11C
22C11.163942
(1)
3X21.000000
(2)190.000(6)
4X11.000000(3)290.000(7)30.000(10)0
53H21.053127(4)390.000(8)1180.000(11)0
64H11.053127(5)490.000(9)2180.000(12)0
------------------------------------------------------------------------
最后一次:
E=-76.82184a.u.,偶极矩为0.0000,原子电荷(a.u.)分别为1C-0.23345,2C-0.23345,3H0.23345,4H0.23345,分子所属点群为D∞h,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11C
22C11.186066
(1)
3X21.000000
(2)190.000(6)
4X11.000000(3)290.000(7)30.000(10)0
53H21.056877(4)390.000(8)1180.000(11)0
64H11.056877(5)490.000(9)2180.000(12)0
------------------------------------------------------------------------
(7)呋喃C4H4O分子的结构
第一次:
E=-228.62293a.u.,偶极矩为1.0915D,原子电荷(a.u.)分别为1O-0.59964,2C0.17372,3C0.17372,4C-0.20870,5C-0.20870,6H0.17647,7H0.17647,8H0.15833,9H0.15833,分子所属点群为C2v,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11O
22C11.349886
(1)
33C11.349886
(2)2105.967(9)
44C21.337001(3)1111.267(10)30.000(16)0
55C31.337001(4)1111.267(11)20.000(17)0
66H21.065250(5)1115.598(12)3180.000(18)0
77H31.065250(6)1115.598(13)2180.000(19)0
88H41.070281(7)2126.956(14)1180.000(20)0
99H51.070281(8)3126.956(15)1180.000(21)0
------------------------------------------------------------------------
最后一次:
E=-228.63294a.u.,偶极矩为0.7714D,原子电荷(a.u.)分别为1O-0.54337,2C0.16901,3C0.16901,4C-0.22309,5C-0.22309,6H0.16617,7H0.16617,8H0.15960,9H0.15960,分子所属点群为C2v,几何构型如下表所示:
CDCentAtomN1Length/XN2Alpha/YN3Beta/ZJ
------------------------------------------------------------------------
11O
22C11.343671
(1)
33C11.343671
(2)2107.146(9)
44C21.339105(3)1110.794(10)30.000(16)0
55C31.339105(4)1110.794(11)20.000(17)0
66H21.068418(5)1116.292(12)3180.000(18)0
77H31.068418(6)1116.292(13)2180.000(19)0
88H41.070276(7)2126.704(14)1180.000(20)0
99H51.070276(8)3126.704(15)1180.000(21)0
------------------------------------------------------------------------
【思考与讨论】
1、内坐标的定义,特别是平面和非平面二面角的右手规则。
内坐标是描述体系内部运动的坐标。
坐标平面内坐标轴上的点横坐标x与纵坐标y二者至少有一个为0,当且仅当该点是原点时x与y同时为0,所以有xy=0.对于N原子分子,其内坐标共有3N-6(直线型分子为3N-5)个。
常用的内坐标有三种,键长、键角和二面角,即用键长、键角和二面角定义原子核的位置。
平面二面角的右手规则:
用右手,大