实验64二氟二氯甲烷分子结构的从头算法 宏观物体的运动状态用位置、运动速度和加速度来描述，其运动规律遵循牛顿三定律。但对 于质量很小的微观粒子，如分子、原子、电子等，其运动有两个不同于宏观物体的特点：量子化 和波粒二象性，这些粒子既具有波动性又具有粒子性，其能量和动量是不连续的，作跳跃式增 减，经典力学无法反映微观粒子的这些特点。在20世纪初，为了描述微观粒子的这些运动行为， 物理学家们建立了新的科学一一量子力学。在量子力学的理论框架下，微观粒子的状态用波函数 来描述，波函数的统计解释是波粒二象性的一个体现，波函数随时间的变化即粒子的运动规律满 足薛定谔方程。量子力学基本思想的正确性已由各种情况下用此方程得出的结论与实验结果相符 而得到了验证。 分子结构理论的主要任务是根据原子核和电子相互作用及其运动的规律，运用量子力学原 理，探求分子的化学组成、结构和性质等化学规律。分子轨道理论是一种广泛应用的分子结构理 论，从头计算方法即是在分子轨道理论基础上发展起来的，由于在理论上的严格性和计算结果的 可靠性，使它在各种量子化学计算方法中居于主导地位。它不仅是理论化学研究的一种必不可少 的理论工具，并且广泛地应用于化学、固体物理学、生物学、医药学等领域。 【实验目的】 1.了解分子轨道理论和从头算法的基础知识。 2.学习分子构型的建立和优化 3.掌握分子结构计算的基本方法和相关软件的使用。 【预习要求】 1.什么是从头算法？ 2.什么是Hlatree-Fork方法和密度泛函理论？其基本思想是什么？ 3.什么是基组？ 4.二氟二氯甲烷分子的电子组态如何？哪些是价壳层？哪些是内价壳层？哪些是内层轨道？ 【实验原理】 一、从头算法 分子是一个多粒子体系，运动十分复杂。除了电子运动外，还包括分子作为整体的平动和转 动以及分子内各原子核之间的相对振动。同时还存在电子的自旋和轨道耦合，电子自旋与核自旋 之间以及核与核之间的自旋相互作用等。理论计算上要同时考虑这些方面显然是不实际的，必须 分清主次，作合理的近似。电子的质量较核轻得多，而且电子的运动速度比核的速度大得多，引 入Borm-Oppenheimeri近似，把核运动与电子运动分开考虑，进而引入单电子近似，认为每个电子都 在各核和其他电子组成的平均作用势场中独立地运动，做了这样的近似后，电子的运动状态用单 电子函数，米描述，这种单电子函数就称为分子轨道。从头算法就是在这些近似的基础上，仅仅 利用普朗克常数、申子质量和电量三个基本物理常数及元素的原子序数，不借助于任何经验参 数，计算体系全部电子的分子积分，近似求解薛定谔方程。以下介绍从头算法的Hatree-Fork方法和
实验 6-4 二氟二氯甲烷分子结构的从头算法 宏观物体的运动状态用位置、运动速度和加速度来描述，其运动规律遵循牛顿三定律。但对 于质量很小的微观粒子，如分子、原子、电子等，其运动有两个不同于宏观物体的特点：量子化 和波粒二象性，这些粒子既具有波动性又具有粒子性，其能量和动量是不连续的，作跳跃式增 减，经典力学无法反映微观粒子的这些特点。在 20 世纪初，为了描述微观粒子的这些运动行为， 物理学家们建立了新的科学——量子力学。在量子力学的理论框架下，微观粒子的状态用波函数 来描述，波函数的统计解释是波粒二象性的一个体现，波函数随时间的变化即粒子的运动规律满 足薛定谔方程。量子力学基本思想的正确性已由各种情况下用此方程得出的结论与实验结果相符 而得到了验证。 分子结构理论的主要任务是根据原子核和电子相互作用及其运动的规律，运用量子力学原 理，探求分子的化学组成、结构和性质等化学规律。分子轨道理论是一种广泛应用的分子结构理 论，从头计算方法即是在分子轨道理论基础上发展起来的，由于在理论上的严格性和计算结果的 可靠性，使它在各种量子化学计算方法中居于主导地位。它不仅是理论化学研究的一种必不可少 的理论工具，并且广泛地应用于化学、固体物理学、生物学、医药学等领域。 【实验目的】 1．了解分子轨道理论和从头算法的基础知识。 2．学习分子构型的建立和优化。 3．掌握分子结构计算的基本方法和相关软件的使用。 【预习要求】 1． 什么是从头算法？ 2． 什么是Hatree-Fork方法和密度泛函理论？其基本思想是什么？ 3． 什么是基组？ 4． 二氟二氯甲烷分子的电子组态如何？哪些是价壳层？哪些是内价壳层？哪些是内层轨道？ 【实验原理】 一、从头算法 分子是一个多粒子体系，运动十分复杂。除了电子运动外，还包括分子作为整体的平动和转 动以及分子内各原子核之间的相对振动。同时还存在电子的自旋和轨道耦合，电子自旋与核自旋 之间以及核与核之间的自旋相互作用等。理论计算上要同时考虑这些方面显然是不实际的，必须 分清主次，作合理的近似。电子的质量较核轻得多，而且电子的运动速度比核的速度大得多，引 入Born-Oppenheimer近似，把核运动与电子运动分开考虑，进而引入单电子近似，认为每个电子都 在各核和其他电子组成的平均作用势场中独立地运动，做了这样的近似后，电子的运动状态用单 电子函数 i 来描述，这种单电子函数就称为分子轨道。从头算法就是在这些近似的基础上，仅仅 利用普朗克常数、电子质量和电量三个基本物理常数及元素的原子序数，不借助于任何经验参 数，计算体系全部电子的分子积分，近似求解薛定谔方程。以下介绍从头算法的Hatree-Fork方法和
密度泛函理论。 二、Hatree-Fork方法(HF) 考虑到电子是费米子，分子体系波函数必须满足交换反对称要求，则闭壳层分子体系（含2 个电子)波函数用Slater行列式表示 Ψ= 2∑-lyPw(2.y(2n-lw4(2my (6-4-1) 其中a和B分别表示自旋向上和自旋向下的电子，P是置换算符，P为置换次数，"，a和W,B为分 子的单电子轨道波函数？，和电子自旋波函数的乘积。 在做了上述近似后，由N个原子核和2个电子组成的闭壳层分子体系的哈密顿量可以表示为 (这里取原子单位，此时电荷的原子单位为1个电子所带电量，能量的原子单位为272V) H=H+H, (6-4-2) H=-∑+∑21=24 (6-4.3) 4=-∑22 (6-44) %-羽 (6-45) 其中，H,称为单电子算符，是所有电子动能及核库仑势能的和：H,称为双电子算符，是分子内 所有电子两两之间的排斥势能之和：h,代表了分子内单个电子的动能和电子与原子核之间的库仑 势能之和。 分子体系能量的期望值为 E=(平|HI平) (64-6) 将(64-2)、(643)、(644)和(6-45)式代入上式可得 B=229-22,-K,) (6-4-7) 其中 ，-ae22Io0p②） (6-4-8) K,-(a0ne合o2e0 (6-4.9)
密度泛函理论。 二、Hatree-Fork方法（HF） 考虑到电子是费米子，分子体系波函数必须满足交换反对称要求，则闭壳层分子体系（含 2n 个电子）波函数用Slater行列式表示 1 1 1 ( 1) { (1) (2) (2 1) (2 )} (2 )! p n n p P n n n （6-4-1） 其中 和 分别表示自旋向上和自旋向下的电子， P 是置换算符， p 为置换次数, i 和 i 为分 子的单电子轨道波函数 i 和电子自旋波函数的乘积。 在做了上述近似后，由N个原子核和 2n个电子组成的闭壳层分子体系的哈密顿量可以表示为 （这里取原子单位，此时电荷的原子单位为 1 个电子所带电量，能量的原子单位为 27.2eV） H H H 1 2 （6-4-2） 2 2 1 1 1 [ ] 2 n i i i i i Z H h r （6-4-3） 1 2 2 i i i Z h r （6-4-4） 2 1 i j ij H r （6-4-5） 其中， H1 称为单电子算符，是所有电子动能及核库仑势能的和； H2 称为双电子算符，是分子内 所有电子两两之间的排斥势能之和； i h 代表了分子内单个电子的动能和电子与原子核之间的库仑 势能之和。 分子体系能量的期望值为 E H | | （6-4-6） 将（6-4-2）、（6-4-3）、（6-4-4）和（6-4-5）式代入上式可得 1 2 (2 ) n i ij ij i i i j E J K （6-4-7） 其中， 12 1 (1) (2) | | (1) (2) ij i j i j J r （6-4-8） 12 1 (1) (2) | | (2) (1) Kij i j i j r （6-4-9）
6,=(@,0lh01g)+∑(2，-K) (6-4-10) 应用变分原理并利用正交归一化条件(,》=，对(647)式作变分，调整一组{9， 使能量最小，经过简化后得到 FD=Φ (6-4-11) 其中F、E和分别为F(Hartrec-.Fock哈密顿算符)、G,和g,的矩阵形式，且 F=h+∑(2J-K,) F0,=80 (6-4-12) J,0=e212K9,=((2112e0 (64-13 i, (6-4-12)式称为Hartree-Fock方程，方程的解是一组最优化的分子轨道，称为Hartree-Fock分子轨 道。Hartree-Fock方程直接求解十分困难，有效的方法是将Hartree-.Fock分子轨道用原子轨道的线性 组合来近似（缩写为LCAO-MO近似） (6-4-14) 其中，是一组原子轨道，cm是组合系数。将(6414)式代入(6-47)一(64-10)式，经过 与前面类似的变分后可得 ∑(Enm-e,Sm)n=0(i=l2,n:4,v=l2,m） (6-4-15) 其中 Fm=uhIv)+∑∑Pw2(uwla-(ol2vp S=(uv). (urio）=∬x0x.0xi2)z.2drda .2zd, Pn=2∑c (6-4-15)式称为Hartree-Fock-Roothaan(HFR)方程，其矩阵形式为： FC-SC8 (6-4-16)
(1) (1) | (1) (2 ) i i i ij ij j h J K （6-4-10） 应用变分原理并利用正交归一化条件 i j ij ，对（6-4-7）式作变分，调整一组{ i }， 使能量最小，经过简化后得到 F ε （6-4-11） 其中 F 、 ε 和 分别为 F （Hartree-Fock哈密顿算符）、 i 和 i 的矩阵形式，且 (2 ) i j j j F h J K F i i i （6-4-12） 12 1 (1) (2) | | (2) j j j J r 12 1 (1) (1) (2) | | (2) (1) Kj j j i j r （6-4-13） （6-4-12）式称为Hartree-Fock方程，方程的解是一组最优化的分子轨道，称为Hartree-Fock分子轨 道。Hartree-Fock方程直接求解十分困难，有效的方法是将Hartree-Fock分子轨道用原子轨道的线性 组合来近似（缩写为LCAO-MO近似） m i i c （6-4-14） 其中， 是一组原子轨道， i c 是组合系数。将（6-4-14）式代入（6-4-7）—（6-4-10）式, 经过 与前面类似的变分后可得 ( ) , ) 0 ( 1, 2, , 1, 2, , F S c i i i n m ； （6-4-15） 其中， F h P | (1) | {2 } S | ， 1 2 12 1 2 12 1 (1) (1) (2) (2) 1 (1) (1) (2) (2) v v d d r d d r 1 2 n j j j P c c （6-4-15）式称为Hartree-Fock-Roothaan（HFR）方程，其矩阵形式为： FC=SCε （6-4-16）
通过求解上式就可以得到分子的单电子波函数及轨道能量。首先要假设一套系数C,代入上式计算 得到一组本征能量值￡，再带入(6416)式得到新的C矩阵，重复此过程，直到相继两次的变化 在误差范围内，则认为已经达到自洽收敛，求解完成。 二、密度泛函理论(0FT) 1.DFT方法介绍 HF方法中，仅考虑了电子之间的自旋交换效应，没有考虑到电子之间的关联作用。而密度泛 函方法则通过引入交换关联泛函同时考虑到了这两种相互作用，在计算很多分子尤其是大分子的 电子结构方面显示了它的优越性。 1964年，Hohenberg和Kohn提出，分子的基态能量和波函数只由电子概率密度p(r)决定，分 子基态能量 E。=Ep】(这里[P表示是pr)的泛函关系) (6-4-17) l96s年，Kohn和Sham提出了分子的单电子的Kohn-Sham方程 12(e()dr+V-Ip(()(r)(64-18) 其中等号右端的e,为轨道电子的能量本征值，鸟为Kohn-Sham单电子轨道波函数：等号左端第一 项为轨道电子动能，第二项为轨道电子与核之间的库仑势能，第三项为轨道电子与其他电子之间 的排斥势能，第四项代表了电子之间的交换一关联势，其表达式为 '[pr】=6 EIp(r)]/o (6-4-19) 其中E【(r】称为交换一关联能泛函，它分为交换泛函和关联泛函两项，即 E[P]=E,Ie]+E[P] (6-4-20) E[P]和Ep的精确形式很难得到，各种近似的E,[p和E[P]在计算中被广泛使用。 分子的电子密度p(r)可由单电子轨道，(r)确定 pr)=∑14(r)f (6-4-21) 由此，Kohn和Sham进一步给出了含有n个电子和m个原子核的分子的基态精确能量 E=2m1v1叭-子 F-RP (64-22) *号rear+E.pl 其中第一项为所有的电子动能，第二项为所有电子与核的库仑势能，第三项为所有电子之间的排 斥势能，第四项为电子之间的交换一关联能。 在计算过程中，先给一个初始值P,可求得选定的E[p],代入(64-19)式可得到'c,再 由Kohn-Sham方程(6-4-8)得到一套Kohn-Sham轨道{}和能量本征值{e,},将{}代入(6-4 21)式得到改进的P,再重新开始.，如此反复循环直到收敛为止
通过求解上式就可以得到分子的单电子波函数及轨道能量。首先要假设一套系数C，代入上式计算 得到一组本征能量值 ε ，再带入（6-4-16）式得到新的C矩阵，重复此过程，直到相继两次的变化 在误差范围内，则认为已经达到自洽收敛，求解完成。 二、密度泛函理论（DFT） 1.DFT方法介绍 HF方法中，仅考虑了电子之间的自旋交换效应，没有考虑到电子之间的关联作用。而密度泛 函方法则通过引入交换关联泛函同时考虑到了这两种相互作用，在计算很多分子尤其是大分子的 电子结构方面显示了它的优越性。 1964 年，Hohenberg和Kohn提出，分子的基态能量和波函数只由电子概率密度 ( ) r 决定，分 子基态能量 0 0 [ ] （这里 [ ] 表示是 ( ) r 的泛函关系） （6-4-17） 1965 年，Kohn和Sham提出了分子的单电子的Kohn-Sham方程 2 1 1 ( ) ( ) ( [ ( )]) ( ) ( ) 2 2 xc i i i Z d V ' ' ' r r r r r r r R r r （6-4-18） 其中等号右端的 i 为轨道电子的能量本征值， i 为Kohn-Sham单电子轨道波函数；等号左端第一 项为轨道电子动能，第二项为轨道电子与核之间的库仑势能，第三项为轨道电子与其他电子之间 的排斥势能，第四项代表了电子之间的交换—关联势，其表达式为 [ ( )] [ ( )]/ V E xc xc r r （6-4-19） 其中 [ ( )] E xc r 称为交换—关联能泛函，它分为交换泛函和关联泛函两项，即 [ ] [ ] [ ] E E E xc x c （6-4-20） [ ] E x 和 [ ] E c 的精确形式很难得到，各种近似的 [ ] E x 和 [ ] E c 在计算中被广泛使用。 分子的电子密度 ( ) r 可由单电子轨道 ( ) i r 确定 2 1 ( ) | ( ) | n i i r r （6-4-21） 由此，Kohn和Sham进一步给出了含有n个电子和m个原子核的分子的基态精确能量 2 0 1 1 E ( ) | | ( ) ( ) 2 1 ( ) ( ) [ ] 2 n m i i i xc Z d d d E ' ' ' r r r r r R r r r r r r （6-4-22） 其中第一项为所有的电子动能，第二项为所有电子与核的库仑势能，第三项为所有电子之间的排 斥势能，第四项为电子之间的交换—关联能。 在计算过程中，先给一个初始值 ，可求得选定的 [ ] E xc ，代入（6-4-19）式可得到 Vxc ，再 由Kohn-Sham方程（6-4-18）得到一套Kohn-Sham轨道 { } i 和能量本征值 { }i ，将 { } i 代入（6-4- 21）式得到改进的 ，再重新开始 .，如此反复循环直到收敛为止
2.交换和关联泛丽 下面简要介绍一下常用的交换和关联泛函。 (1)Slater和Becke交换泛函 Slater定域自旋交换泛函（简写为S)是Slater基于均匀电子气模型用统计方法导出的电子交换 泛函，其形式为 E=-a2 ySrdr且a-号 (6-4-23) 均匀电子气模型认为无限个电子在连续均匀分布的正电荷气氛中运动，整个体系呈电中性而且在 空间无限延伸。考虑到分子中的正电荷不是均匀分布的，它仅位于各原子核处，P应随空间位置 不同而迅速变化，所以在(6-4-23)式为P(r)的积分。若考虑到分子体系并不是无限大，而且 电子密度也不均匀，有效的方法是让交换泛函包括电子密度的非局域梯度项。如常用的Bcke交换 泛函（简写为B) E.-E.6xsinh (6-4-24) 其中,b为可调参数 (2)LYP关联泛函 LYP关联泛函是Lce,Yag和Par在1988年提出的关联泛函，其中包括了电子密度的非局域梯 度修正项。LYP关联泛函的形式比较复杂，这里不再给出，可查找相关文献资料。 在进行密度泛函理论的计算时可以选用不同的交换泛函与关联泛函搭配进行分子轨道的理论 计算。 (3)B3LYP泛函 在杂化泛函中，交换泛函是HF交换项、定域交换项和梯度修正交换项的线性组合，再加上定 域和（或）非定域梯度修正的关联泛函。B3LYP杂化泛函是Becke于1993年提出的杂化泛函，是 一种常用的杂化泛函，其具体形式为 EmL=ES+Co(EIF-ES)+CAE+EN+C(ELWP-EIWN3)(6-4-25) 其中E俨表示HF交换项，E3表示Vosko,Wlk和Nusair在1980年提出的关联泛函Ⅲ。(6-4 25)式中含有三个参数，所以又被称为Becke三参数泛函。三个参数值分别为：c0=0.20, c.=0.72和c.=0.81。 三、基组 前面提到（式(6414）)，分子的单电子轨道波函数用原子轨道的线性组合来表示，在计算 时，一般可以用事先定义好的一组单电子函数{之}来代替原子轨道，么就称为基函数。在计算 过程中，基函数的选取十分重要。基函数应能较真实地反映分子系统的波函数，同时应减少计算 困难，有效地节省计算时间。下面介绍常用的高斯型基函数(GTO)。 原子中电子轨道由径向部分R(r)和角向部分Y(，p)的乘积表示
2. 交换和关联泛函 下面简要介绍一下常用的交换和关联泛函。 （1）Slater和Becke交换泛函 Slater定域自旋交换泛函（简写为S）是Slater 基于均匀电子气模型用统计方法导出的电子交换 泛函，其形式为 3 3 1/3 4/3 ( ) ( ) 2 4 S E d x r 且 2 3 （6-4-23） 均匀电子气模型认为无限个电子在连续均匀分布的正电荷气氛中运动，整个体系呈电中性而且在 空间无限延伸。考虑到分子中的正电荷不是均匀分布的，它仅位于各原子核处， 应随空间位置 不同而迅速变化，所以在（6-4-23）式为 4/ 3 ( ) r 的积分。若考虑到分子体系并不是无限大，而且 电子密度也不均匀，有效的方法是让交换泛函包括电子密度的非局域梯度项。如常用的Becke交换 泛函（简写为B） 4/3 2 1 [ ] 1 6 sinh B S x x x E E b d bx x （6-4-24） 其中 4/3 x | | ，b 为可调参数。 （2）LYP关联泛函 LYP关联泛函是Lee，Yang 和Parr在 1988 年提出的关联泛函，其中包括了电子密度的非局域梯 度修正项。LYP关联泛函的形式比较复杂，这里不再给出，可查找相关文献资料。 在进行密度泛函理论的计算时可以选用不同的交换泛函与关联泛函搭配进行分子轨道的理论 计算。 （3）B3LYP泛函 在杂化泛函中，交换泛函是HF交换项、定域交换项和梯度修正交换项的线性组合，再加上定 域和（或）非定域梯度修正的关联泛函。B3LYP杂化泛函是Becke于 1993 年提出的杂化泛函，是 一种常用的杂化泛函，其具体形式为 3 3 3 0 ( ) ( ) B LYP S HF S B VWN LYP VWN E E c E E c E E c E E xc x x x x x c c c c （6-4-25） 其中 HF E x 表示HF交换项， VWN3 E c 表示Vosko，Wilk 和Nusair在 1980 年提出的关联泛函III。（6-4- 25）式中含有三个参数，所以又被称为Becke三参数泛函。三个参数值分别为： 0 c 0.20 ， 0.72 x c 和 0.81 c c 。 三、基组 前面提到（式（6-4-14）），分子的单电子轨道波函数用原子轨道的线性组合来表示，在计算 时，一般可以用事先定义好的一组单电子函数 来代替原子轨道， 就称为基函数。在计算 过程中，基函数的选取十分重要。基函数应能较真实地反映分子系统的波函数，同时应减少计算 困难，有效地节省计算时间。下面介绍常用的高斯型基函数（GTO）。 原子中电子轨道由径向部分 R r( ) 和角向部分 Y( , ) 的乘积表示
(r.0.p)=R(r)Y.(0.) (6-4-26) 高斯型基函数的径向部分为 R(a,r)=cx"y"=exp(-ar2) (6-4-27) 式中是高斯指数，c是归一化常数。角向部分同原子角向函数。高斯型函数的缺点是r很大或很 小时都不能真实地描述原子的轨道波函数。解决的办法是选用缩并的高斯型基函数，即选择几个 高撕斯型基函数，用固定的线性组合系数重新组合一个新的函数（通常称为STO-G),这就是基 组。下面介绍几种常用的基组，在计算程序中己经存储了这些基组。 (1)6-31G基组 每个内层轨道用个ST0-6G拟合，而每个价轨道用个STO3G和一个ST0.1G拟合。这是考 虑到在分子中，原子内层电子轨道形变较小，而价电子轨道形变较大，所以把内层轨道和价轨道 分开处理。 (2)6-311G基组 每个内层轨道用一个ST0-6G拟合，而每个价轨道用一个ST0-3G,一个STO-1G和另一个ST0- 1G拟合。 (3)6-311++G基组 在6311G基组的基础上对每个非氢原子加上一组s型和一组即型的高斯函数。对每个氢原子加 上一组s型的高斯函数。这是考虑到分子的某些电子轨道，在大处仍有可观的电子密度，加扩散基 函数可以较好地反映这一点。 (4)6-311++G(dD) 在6-311+G基组的基础上对非氢原子加上一组型的高斯函数，对氢原子加上一组p型的高斯 函数。加极化高斯函数的目的是为了更好地反映分子内电子所受电场与在原子环境中的不同。 (5)6-311++G(3df3pd) 在6-311+G基组的基础上对非氢原子加上三组d型和一组f型的高斯函数，对氢原子加上三组p 型和一组d型的高斯函数。 (6)ce-pVTZ和cc-pVQ☑ Dunning提出的关联一致基组，其中已为各原子增加了极化基函数。cc-pVTZ和cc-pVQZ基组 的具体形式见表641。 表6-4-1.cc-DVTZ和cc-pVOZ基组 原子 ce-pVTZ ce-pVQZ H和He 38,2p,1d 4s,3p,2d1f B-Ne 45,3p,2d,1f 5s,4p,3d,2f,1g Al-Ar 5s,4p.2d.If 6s,5p,3d2￡1g (7)AUG-ce-pVTZHAUG-cc-pVQZ 分别对上面两个基组添加扩散函数。AUG-cc-pVTZ在cc-pVTZ的基础上为氢原子添加一组s 型，一组p型和一组d型高斯函数，为B-Ne和A-Ar添加一组s型，一组p型，一组d型和一组型高斯 函数。AUG-cc-pVQZ在cc-pVQZ的基础上为氢原子添加一组s型，一组p型，一组d型和一组型扩高
( , , ) ( ) ( , ) nlm nl lm r R r Y （6-4-26） 高斯型基函数的径向部分为 2 ( , ) c exp( ) n m l R r x y z r （6-4-27） 式中 是高斯指数，c是归一化常数。角向部分同原子角向函数。高斯型函数的缺点是r很大或r很 小时都不能真实地描述原子的轨道波函数。解决的办法是选用缩并的高斯型基函数，即选择几个 高斯型基函数，用固定的线性组合系数重新组合一个新的函数（通常称为STO-nG），这就是基 组。下面介绍几种常用的基组，在计算程序中已经存储了这些基组。 （1） 6-31G 基组 每个内层轨道用一个STO-6G拟合，而每个价轨道用一个STO-3G和一个STO-1G拟合。这是考 虑到在分子中，原子内层电子轨道形变较小，而价电子轨道形变较大，所以把内层轨道和价轨道 分开处理。 （2） 6-311G 基组 每个内层轨道用一个STO-6G拟合，而每个价轨道用一个STO-3G，一个STO-1G和另一个STO- 1G拟合。 （3） 6-311++G基组 在 6-311G基组的基础上对每个非氢原子加上一组s型和一组p型的高斯函数。对每个氢原子加 上一组s型的高斯函数。这是考虑到分子的某些电子轨道，在大r处仍有可观的电子密度，加扩散基 函数可以较好地反映这一点。 （4） 6-311++G（d,p） 在 6-311++G基组的基础上对非氢原子加上一组d型的高斯函数，对氢原子加上一组p型的高斯 函数。加极化高斯函数的目的是为了更好地反映分子内电子所受电场与在原子环境中的不同。 （5） 6-311++G（3df,3pd） 在 6-311++G基组的基础上对非氢原子加上三组d型和一组f型的高斯函数，对氢原子加上三组p 型和一组d型的高斯函数。 （6） cc-pVTZ和cc-pVQZ Dunning 提出的关联一致基组，其中已为各原子增加了极化基函数。cc-pVTZ和cc-pVQZ基组 的具体形式见表 6-4-1。 表 6-4-1. cc-pVTZ和cc-pVQZ基组 （7） AUG-cc-pVTZ和AUG-cc-pVQZ 分别对上面两个基组添加扩散函数。AUG-cc-pVTZ在cc-pVTZ的基础上为氢原子添加一组s 型，一组p型和一组d型高斯函数，为B-Ne和Al-Ar添加一组s型，一组p型，一组d型和一组f型高斯 函数。AUG-cc-pVQZ在cc-pVQZ的基础上为氢原子添加一组s型，一组p型，一组d型和一组f型扩高 原子 cc-pVTZ cc-pVQZ H和He B-Ne Al-Ar 3s, 2p, 1d 4s, 3p, 2d, 1f 5s, 4p, 2d, 1f 4s, 3p, 2d, 1f 5s, 4p, 3d, 2f, 1g 6s, 5p, 3d, 2f, 1g
斯函数，为B-N和A-Ar添加一组s型，一组p型，一组s型，一组型和一组g型高斯函数。 表6.42.CF,C1,分子的波函 基组 原始高斯函药 总个数 缩并后的高斯函数 总个数 3 2 04 UG-cc-DVOZ 2221p243030g 6s15024d30E30e 71824d30.30g1 四、CFC,分子的电子组态和轨道电离能的计算 CF,C2分子的电子组态为： (core)(1a,b,产2a,广b, inner valence (3a,F(2b(4a2b,)a2(5a(3b,(3b,)(6a(4b,)2(2a,24b2 outer valence 表643中是光电子谱实验和前人用0VGF方法选用6-311+G(,p)基组计算得到的外层轨道电 离能。从表中可见理论计算的轨道顺序与实验顺序有差异，可见理论计算仍有待进一步提高，但 理论计算的电离能与实验值已经很接近，这说明理论计算结果是可靠的。 表64-3.前人光电子谱及理论计算得到的CF,C分子外价壳层轨道的垂直电离能（单位V) 光电 31 b119 +G(d.p) 21284 5a b204 .62 a122.4 3a23.35) 本实验中利用HF和DFT-B3LYP两种理论方法选择6-311G、6-311+G、6-311+G(dp)、cc pVTZ和AUG-cc-pVTZ对CFCh分子外层电子轨道和电离能进行计算，并利用DFT-BLYP/6-31I+G 计算电离能较低的前六个轨道的电子密度分布，计算程序利用Chem3D和Gaussian09程序包。 【实验器材】 计算机一台。 【实验内容】
斯函数，为B-Ne和Al-Ar添加一组s型，一组p型，一组s型，一组f型和一组g型高斯函数。 表 6-4-2. CF2Cl2分子的波函数 基组 原始高斯函数 总个数 缩并后的高斯函数 总个数 6-31G 6-311G 6-311++G 6-311++G (d,p) 6-311++G (3df,3pd) AUG-cc-pVTZ AUG-cc-pvQZ C, F (10s,12p) (11s,15p) (12s,18p) (12s,18p,6d) (12s,18p,18d,10f) (19s,18p,18d,20f) (22s,21p,24d,30f,30g) Cl (16s,30p) (13s,30p) (14s,33p) (14s,33p,6d) (14s,33p,18d,10f) (42s,51p,18d,20f) (43s,60p,24d,30f,30g) 158 164 184 214 324 457 755 C, F Cl (3s,6p) (4s,9p) (4s,9p) (6s,15p) (5s,12p) (7s,18p) (5s,12p,6d) (7s,18p,6d) (5s,12p,18d,10f) (7s,18p,18d,10f) (5s,12p,18d,20f) (6s,15p,18d,20f) (6s,15p,24d,30f,30g) (7s,18p,24d,30f,30g) 53 81 101 131 241 283 533 四、CF2Cl2分子的电子组态和轨道电离能的计算 CF2Cl2分子的电子组态为： 26 2 2 2 2 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 1 2 1 1 2 1 1 2 1 1 2 2 (core) (1a ) (1b ) (2a ) (1b ) inner valence (3a ) (2b ) (4a ) (2b ) (1a ) (5a ) (3b ) (3b ) (6a ) (4b ) (2a ) (4b ) outer valence 表 6-4-3 中是光电子谱实验和前人用 OVGF 方法选用 6-311+G(d, p)基组计算得到的外层轨道电 离能。从表中可见理论计算的轨道顺序与实验顺序有差异，可见理论计算仍有待进一步提高，但 理论计算的电离能与实验值已经很接近，这说明理论计算结果是可靠的。 表 6-4-3. 前人光电子谱及理论计算得到的CF2Cl2分子外价壳层轨道的垂直电离能（单位为eV） 光电子谱实验值 OVGF/6-311+G(d,p) 4b2(12.26) 2a2 (12.53) 4b1 (13.11) 6a1(13.45) 3b2(14.36) 3b1(15.9) 5a1 (16.3) 1a2 (16.9) 2b1(19.3) 4a1(19.3) 2b2(20.4) 3a1(22.4) 4b2(11.96) 2a2(12.33) 4b1(12.84) 6a1(13.15) 3b2(14.18) 5a1(16.29) 3b1(16.41) 1a2(17.21) 2b2(19.44) 4a1(19.51) 2b1(20.62) 3a1(23.35) 本实验中利用HF和DFT-B3LYP两种理论方法选择 6-311G、6-311++G、6-311++G（d,p）、ccpVTZ和AUG-cc-pVTZ对CF2Cl2 分子外层电子轨道和电离能进行计算，并利用DFT-BLYP/6-311++G 计算电离能较低的前六个轨道的电子密度分布，计算程序利用Chem3D和Gaussian09 程序包。 【实验器材】 计算机一台。 【实验内容】
1.基础性实验内容 (1)利用Chem3D程序绘制CFCl2分子构型并做初步能量优化： (2)利用Gaussian09程序进行结构优化。选用HF方法，基组选用6-311+G和6-311+G (dp),比较两个基组的优化结果。 (3)利用GaussianO9程序计算分子轨道和电离能。选用HF和DFT-BLYP两种方法，基组选用 631IG、631+G、6-311+G(dp)、ce-pVTZ和AUG-c-pVTZ基组，列表比较计算结果，并与 前人和理论计算作比较。 (4)利用Chem3D程序绘制能量最低的六个轨道的电子密度分布，计算方法和基组选用DFT BLYP6-311+G,根据电子密度分布分析分子中原子间的成键特性。 2.设计性实验内容 选择几种分子，利用Gaussian0(9程序选用不同方法和基组计算它们的红外光谱(IR光谱)，查 找前人的R光谱测量结果，将计算结果与前人研究结果进行对比，分析不同方法和基组在计算R 光谱方面的差异。 【注意事项】 (1)计算前认真阅读Chem3D和GaussianO9的相关操作说明。 (2)不同方法和基组计算得到的轨道顺序可能不同，要注意计算结果中的轨道序数与轨道标 识应一一对应。 (3)由于计算方法和基组选用较多，因此计算内容和相关文0件应建立对应的文件夹和文件， 并妥善命名和保存。 【思考题】 1,什么是HF方法和密度泛函理论？ 2.如何建立二氟二氯甲烷的分子模型 3.如何对分子进行结构优化？ 4.如何计算分子的轨道电离能？ 5.如何计算分子的轨道电子密度？从轨道电子密度能得到分子轨道的什么信息？ 【参考文献】 Gaussian0(9说明书 [2]Chem3D说明书. B)陈正隆，孙小强.量子化学.南京大学出版社，2007。 4新瑞发，龙梅.高斯量子化学计算与应用.内蒙古科学技术出版社，2009. 5)曾炜，孙延辉，李红.Chem0ice2008实用教程.化学工业出版社，2009 【6]张培新，陈建华，魏群，等.掺杂材料分子模拟与计算.科学出版社，2012
1. 基础性实验内容 （1）利用Chem3D程序绘制CF2Cl2分子构型并做初步能量优化。 （2）利用Gaussian09 程序进行结构优化。选用HF方法，基组选用 6-311++G和 6-311++G （d,p），比较两个基组的优化结果。 （3）利用Gaussian09 程序计算分子轨道和电离能。选用HF和DFT-BLYP两种方法，基组选用 6-311G、6-311++G、6-311++G（d,p）、cc-pVTZ和AUG-cc-pVTZ 基组，列表比较计算结果，并与 前人和理论计算作比较。 （4）利用Chem3D程序绘制能量最低的六个轨道的电子密度分布，计算方法和基组选用DFTBLYP/6-311++G，根据电子密度分布分析分子中原子间的成键特性。 2. 设计性实验内容 选择几种分子，利用Gaussian09 程序选用不同方法和基组计算它们的红外光谱（IR光谱），查 找前人的IR光谱测量结果，将计算结果与前人研究结果进行对比，分析不同方法和基组在计算IR 光谱方面的差异。 【注意事项】 （1）计算前认真阅读Chem3D和Gaussian09 的相关操作说明。 （2）不同方法和基组计算得到的轨道顺序可能不同，要注意计算结果中的轨道序数与轨道标 识应一一对应。 （3）由于计算方法和基组选用较多，因此计算内容和相关文 0 件应建立对应的文件夹和文件， 并妥善命名和保存。 【思考题】 1．什么是HF方法和密度泛函理论？ 2．如何建立二氟二氯甲烷的分子模型？ 3．如何对分子进行结构优化？ 4．如何计算分子的轨道电离能？ 5．如何计算分子的轨道电子密度？从轨道电子密度能得到分子轨道的什么信息？ 【参考文献】 [1] Gaussian09 说明书． [2] Chem3D说明书． [3] 陈正隆，孙小强．量子化学．南京大学出版社，2007． [4] 靳瑞发，龙梅．高斯量子化学计算与应用. 内蒙古科学技术出版社，2009. [5] 曾炜，孙延辉，李红. ChemOffice 2008 实用教程. 化学工业出版社，2009. [6] 张培新，陈建华，魏群，等. 掺杂材料分子模拟与计算. 科学出版社，2012