直链和单环共轭多烯本征行列式递推公式及本征能级的推导
Derivations of the Recurrence Formula for Eigendeterminants and Eigenenergy Levels of Linear and Monocyclic Conjugated Polyenes
通讯作者:
收稿日期: 2021-10-22 接受日期: 2022-12-7
| 基金资助: |
|
Received: 2021-10-22 Accepted: 2022-12-7
The Hückel molecular orbital method is an approximate model to describe the interactions of π-electrons in a conjugate system. In structural chemistry textbooks, the recursive formulas for the eigendeterminants and eigenenergy levels of conjugated π-electrons systems have been directly provided, without derivation. Based on mathematical induction, this work shows how these recursive formulas are derived for linear and monocyclic conjugated polyenes, and clarifies the characteristics and regularity of the eigenenergy levels. This information is suitable for students who lack mathematical knowledge with regard to topics such as the group theory and graph theory.
Keywords:
本文引用格式
张冬菊, 蔡欣懿, 宋其圣.
Zhang Dongju.
共轭多烯是有机化学中一类重要的有机化合物,其结构特征是全部原子均在同一平面,每个碳原子提供一个垂直于分子平面的p轨道,形成离域共轭体系。本文所指共轭多烯仅指直链和单环共轭多烯。1931年休克尔(Hückel)提出了处理共轭体系π电子运动的近似分子轨道模型,简称HMO法。该方法用于理解共轭分子的结构和性质,图像直观、概念清晰,在有机化学中有广泛应用[1]。
HMO方法求解直链和单环共轭多烯π电子薛定谔方程的核心是对Hückel行列式方程的求解,唐敖庆、江元生等[2, 3]前辈基于分子轨道图形理论推导出了直链共扼多烯和轮烯本征多项式的一般表达式,找到了分子整体与其片段之间的分子图形和本征多项式的对应关系,获得了求解高阶久期行列式方法。由于唐-江方法涉及较多的群论、图论等数学知识,在本科教学中通常仅给出本征行列式的递推公式[4–7],但教学过程中发现学生对公式表示难于理解,期望明白公式的来源。本文在早期工作的基础上[8–11],用数学归纳法系统地对直链和单环共轭多烯本征多项式的递推公式及其三角函数表达式进行了重新推导,校正了一些常见的表达错误,推导过程仅需要初等数学知识,方便学生理解,适宜于本科化学类各专业的理论教学。
直链共轭烯烃和单环共轭烯烃的本征行列式分别用gn(x)和pn(x)表示(下标n表示共轭烯烃的C原子数目),即:
令
下面给出这些本征行列式递推公式和对应本征值方程解的通式。
1 直链共轭烯烃本征行列式的递推公式
用数学归纳法证明如下:
所以,直链共轭烯烃本征行列式的递推公式可归纳为(1)式。
2 直链共轭烯烃本征行列式的三角函数通式
《结构化学》教科书中通常直接给出直链共轭烯烃本征行列式的三角函数通式,即:
该通式可用下列两种方法予以导出。
方法一:数学归纳法
令x = 2cosθ,则:
归纳以上结果,即可得到(2)式。
方法二:无穷级数法
将一变量为y的函数f(y)展开为一无穷级数,使其展开系数为gn,
上式分别乘以−xy和y2得:
上述三式相加,(3) + (4) + (5),得:
由递推公式(1)式可知,(6)式右端从第二项开始,所有系数都等于零,从而得:
将x = 2cosθ代入上式,得:
将
由于n = −1时,sin(n + 1)θ = 0
所以,
将(8)式与(3)式比较,即可得到直链共轭烯烃本征行列式的三角函数通式(2)式。
3 直链共轭烯烃整体与片段的递推关系
证明如下:
还可进一步证明:
进一步可证明:
即:一个共轭直链多烯可以拆分为(9)式所示的任意两个片段的组合。
4 单环共轭烯烃本征行列式的递推公式
证明如下:
类似地,将p5、p6等更高阶本征行列式依次展开、降阶,可得:
本征行列式的递推公式可归纳为(10)式。
5 单环共轭烯烃本征行列式的三角函数通式
将(2)式代入(10)式,得:
即单环共轭烯烃本征行列式的三角函数通式为:
6 本征值方程解的通式
(1) 直链共轭烯烃。
由(2)式,得本征值方程:
解之得:
即:
所以
考虑到直链共轭烯烃的HMO轨道关于x = 0的对称分布(图 1),与(14)式等价的表达式如下:
图1
本征行列式gn具有下列性质(容易用数学归纳法予以证明):
(15) 式表明,如果x是本征值方程的一个根,则−x也是其根,即直链共轭烯烃的能级是成对出现的,对称分布。若n为奇数,则当
(2) 单环共轭烯烃。
由(11)式,得单环共轭烯烃的本征值方程为:
即:
n =偶数:
n =奇数:
将(17)、(18)式合并表达为:
所以,单环共轭烯烃π电子能级的通式为:
考虑到n为偶数的单环共轭烯烃HMO轨道关于x = 0的对称分布(图 2),(20)式也可简化为:
图2
由图 2还可以看出,对于n为奇数的情况,当n = 4m + 1时,成键轨道多于反键轨道,共轭分子对应的负离子更稳定,如C5H5−、C9H9−;当n = 4m + 3时,反键轨道多于成键轨道,共轭分子对应的正离子更稳定,如C3H3+、C3H3+。対于n为偶数的情况,当n = 4m + 2时,n个π电子将填充n/2个成键轨道填满,形成稳定的闭壳层结构,比相同数目的直链共轭烯烃更稳定,具有芳香性;当n = 4m时,有两个电子占据两个非键轨道,且自旋平行,构成环辛四烯的三重态(D4d)。三重态的环辛四烯不是稳定结构,将通过形变消除简并,形成非平面的船式结构。通过对这些本征能级特征和规律的认识,有助于学生加深对“结构决定性能”的认识和理解。
作者注意到,在几本近年出版的结构化学教材中,对于单环共轭烯烃本征行列式递推公式(10)和(11)式,常有错误,多是忽略了式中的(−1)n,致使给出的公式仅适用于n为偶数的情况。本文对这些公式进行了统一校正,使其也适用于n为奇数的单环共轭烯烃。
7 结语
本文主要基于数学归纳法系统讨论了直链及单环共轭体系本征行列式的递推公式及本征能级的通式,校正了一些常见的公式表达错误,全部公式推导仅需初等数学知识,便于学生理解和应用, 可作为结构化学课程教学的补充材料供学生学习使用。
参考文献
/
| 〈 |
|
〉 |

