衡量芳香性的方法以及在Multiwfn中的计算 - 图文

式中R'_r和R''_r代表实际分子中的pi键键长,n1和n2分别代表kekule式对应结构的单键和双键的数目,k是力常数。因此HOSE描述了分子从实际几何结构变化到Kekule式描述的电子定域化的几何结构过程中所需的能量的负值,或者说是源自电子离域效应产生的实际结构相对于kekule结构的稳定化能。

关于计算离域化能其它的一些细节讨论和历史上的方法可参见Chem.Rev.,105,3773(2005)。

5.2 芳香稳定化能(Aromatic stabilization energy, ASE)

ASE和离域化能都是基于能量来分析芳香性的标准,体现电子离域导致的能量变化,但ASE的计算方法与之截然不同。ASE计算过程是构建一个开环的等键化学反应。反应能如果算出来是正值,即吸热,说明体系是芳香性;如果是负,即放热,说明是反芳香性。例如苯的ASE可按如下方式计算,能量是B3LYP/6-31G*下得到的

一般认为ASE>5kcal/mol就可以算是有芳香性了。计算ASE时构建的等键反应有任意性,不同的构建方式结果会有所不同,如果构建得不合适,甚至结果会不靠谱,详见

Computational Chemistry-Introduction to the Theory and Applications of Molecular and Quantum Mechanics 2ed(Errol Lewars)一书307页的讨论。

ASE最大的弱点就是计算极其麻烦。对于一个复杂的大体系,或者诸如金属芳香性体系,构建一个合理、有意义的等键反应几乎不可能。构建方式的不唯一性也导致不容易对不同文献的结果进行对比。所以我不认为如今使用ASE讨论芳香性有多大价值,除了对于简单小体系。

5.3 磁共振能(Magnetic resonance energies, MRE)

MRE是一种计算ASE和环电流的方法,详见JACS,128,2873(2006)、JPCA,111,8873(2007)。

MSE基于图论和Huckel理论,只要知道分子拓扑关系即可进行计算,计算时只需要一个参数即Huckel共振积分。此方法只能用于pi芳香性研究。

这个方法计算时要考虑体系全部可能的环路。比如菲不是只包含三个六元环路,而是包括三个六元环路、两个10元环路和一个14元环路。对于每个环路,比如第i个环路,可以计算相应的A_i值,它为正值和负值说明这个环路的感应环电流起抗磁(电流方向符合左手定则)和顺磁(电流方向与左手定则相反)作用。A_i也被称为环路共振能(CRE),将当前体系的所有A_i加起来,得到的就是当前体系的MRE,可视为是离域化能。

根据I_i=4.5*I_0*A_i*S_i/S_0可以计算第i个环路的环电流,其中I_0是苯的环电流,S_i是第i个环的面积,S_0是苯环的面积。将所有环路的环电流累加到一起,就能绘制出整体感应电流图。这和AICD那些方法得到的不一样,那些方法可以得到空间中每个具体的点感应电流密度,而MRE的电流图给出的只是每个键上电流大小和方向信息,如下所示。

MRE方法的优点是不依赖于量化计算,不依赖于具体几何结构,只根据拓扑关系就能迅速得到离域化能和感应电流分布。而且不像大多数芳香性判据那样只能算每个小环的芳香性(比如对菲只能算三个六元环),而是可以研究体系中任意环路的稳定化能和环电流。此方法的缺点也很多,即只能用于平面pi体系、精度仅是HMO级别、不能用于非平衡构型或非基态或开壳层、体现不了键的拉伸产生的影响、不能用于诸如金属簇等体系等等。MRE最适合的体系就是多环芳烃(可以含杂原子),特别是尺寸很大的情况。实现MRE的程序据笔者所知目前没有公开的,不过估计也不会很难写。

6 基于电子离域性的方法

本节介绍的方法都比较直接地与电子离域性相联系。

6.1 多中心键级

多中心键级定义为

其中P是单电子密度矩阵,S是重叠矩阵,a,b,c...代表基函数序号,A,B,C...代表环上的原子序号,并且这些原子是按照拓扑关系依次相邻的。Patrick Bultinck等人的文章鼓吹应当把加和号内的项的各种置换组合都考虑进去,实际上根本没这个必要,分析结果和不考虑置换时一样,而计算量还增大了很多。

多中心键级计算速度很快,理论依据强,是寡人强烈推荐使用的衡量芳香性的方法。它的数值越大说明芳香性越强。它也可以用在开壳层,只要把密度矩阵换成alpha和beta电子的密度矩阵分别计算,然后加和即可。由于对称性的原因,sigma和pi电子的多中心键级是可加和的,因此多中心键级可以分别研究sigma和pi芳香性。多中心键级的一个遗憾是没法用来分析反芳香性,反芳香性的体系只是多中心键级很小而已,故和芳香性很弱的体系没法区分开。另外,虽然多中心键级计算量不大,对于六元环瞬间就能出结果,但是计算量却随着环内原子数以阶乘方式增加,当环很大的时候,如9中心时,计算量就挺大了。

需要注意的是计算多中心键级时切不可使用弥散函数,否则结果会缺乏意义。虽然也有诸如多中心离域化指数(Theor.Chem.Acc.,105,292(2001))等类似的指标能够用来衡量芳香性,可以减小多中心键级的基组依赖性,尤其是能够解决不适合弥散函数的问题,但计算明显更费时。

多中心键级可以在Multiwfn当中计算,目前最多支持六中心,下个版本支持到10中心。启动Multiwfn后输入Gaussian的fch文件的路径(不能用其它类型文件作为输入),然后选9,再选2,输入环上的原子序号(必须在环上依次相邻),如3,7,10,9,8,4,然后六中心键级数值会立刻输出(输入几个原子序号就计算几中心键级)。注意顺时针和逆时针输入环上的原子序号结果会有微小差异,原理上都对,取谁都可以,也可以取平均值。

如果想计算仅由pi或仅由sigma电子贡献的多中心键级,则应先进入主功能100里的功能22。分子在哪个平面上就选哪个平面,然后程序会自动探测出哪些MO是pi轨道。如果接下来想计算sigma部分,就选1将pi轨道的占据数都设为0;如果想计算pi部分,就选2把除了pi轨道以外的MO占据数都设为0。之后再照常计算多中心键级,所用的密度矩阵就是扣掉了sigma或pi贡献的密度矩阵,也即得到了pi或sigma的多中心键级。

(PS:设定轨道占据数也可以用主功能6的选项26,由此可以手动将指定的MO的占据数设为0来排除它的贡献。但是想排除sigma或pi轨道的贡献,用户就得自己去观察轨道图形来找出pi轨道,这比较费事,远不如直接用主功能100里的功能22那么方便。)

6.2 ELF-sigma/pi和LOL-sigma/pi

电子定域化函数(Electron localization function, ELF)是Becke提出的衡量电子定域性的实空间函数,目前有极为广泛的应用。此函数的在《电子定域性的图形分析》

(http://hi.http://www.35331.cn//sobereva/item/4a8aeae2291d8d0f8c3ea8d8)中从图形的角度进行了介绍,在 物理化学学报,27,2786(2011) 中作者对ELF的函数形式和物理意义做了详细的

讨论,因此关于ELF在本文就不再具体介绍了。可以认为,被数值越高的ELF等值面包围的空间,电子越容易在这个空间内离域,同时越难离域到这个空间之外。

在JCP,120,1670(2004)中,作者提出了ELF-sigma和ELF-pi用于研究sigma和pi芳香性。前者就是在计算ELF时只考虑pi以外的轨道,后者就是计算ELF时只考虑pi轨道。注意哪怕是对平面体系,ELF也不能分解为ELF-sigma和ELF-pi的加和,但使用ELF-sigma和ELF-pi研究问题确实有实际价值,所以也就无视这个理论上的小缺陷了。

ELF等值面的数值(isovalue)取得越大,等值面内的空间(称为域)就会越来越小。如果域内含有多个ELF极大点(叫可约域),则随着isovalue的增加,这个连通的域就会逐渐分解为多个小的孤立的域,最终每个域内只包含一个ELF极大点(叫不可约域)。不同isovalue时苯的ELF-pi和ELF-sigma等值面如下图所示

可见,对于ELF-pi等值面,当isovalue=0.91时,在苯环上方、下方对应离域pi电子的连通的可约域恰好要分解为六个碳原子上的不可约域,故0.91就是讨论芳香性时候苯的ELF-pi值。红箭头所示的位置被称为二分点(bifurcation point)。不同体系ELF-pi值不同,pi芳香性越强的体系ELF-pi可约域分解得越晚,即ELF-pi值越大。多数体系并非像苯这样是每个键都对称的,因此ELF-pi可约域在环上的某些键上分解得早,在有的键上分解得晚,ELF-pi应当取最先出现分解时的值,而最先分解处也对应了环共轭最弱的位置。

ELF-sigma函数在构成共价键的原子间总有极大点,对应一块ELF-sigma值较高的区域。衡量芳香性所取的ELF-sigma值指的是环上的这些相邻的高ELF-sigma区域之间的二分点,对于苯,就是上图红箭头所示的位置。由于是在isovalue=0.71时恰好在此处从连通的可约

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4