分子连接性指数的种类

《分子连接性指数的种类》由会员分享，可在线阅读，更多相关《分子连接性指数的种类（12页珍藏版）》请在装配图网上搜索。

1、分子连接性指数的种类A1.陈艳，冯长君.部分有机污染物气相色谱保留指数的QSRR研究J.计算机与 应用化学，200724 (10)： 14041408.理论与方法：距离调和指数T的定义和计算：C用隐氢图来表示化学结构，隐氢图中的顶点和边分别对应分子结构中的基团 (或非氢原子)和化学键。用距离矩阵D描述分子中基团的连接性，设隐氢图中有 n个顶点，则距离矩阵D=d. 乂，dij为分子结构中基团i和基团j之间最短路径 1J n 沟 ij 上的化学键数(也称拓扑距离)。距离矩阵中每一行的矩阵元dij之和以S.表示， 由此构成S矩阵。1J1为了区分杂原子和多重键，对拓扑距离进行修正，用相对键长取代拓扑距

2、离。 相对键长定义如下:规定 C-C 键的键长为 1.0，以便和相邻碳原子间的拓扑距离为 1 一致，其它各类相邻顶点间的拓扑距离则定义为相应的键长与C-C键键长的相 对值，例如:C-0键的相对键长为1.43/1.54=0.9286;C=O键的相对键长为 1.22/1.54=0.7922，依次类推。以乙酸甲醋为例，其隐氢图及D矩阵和S矩阵为：| D 7422C 匚一一 o C-0. 0000LOOOO.000D.92861.7922勺珈1.92862. 85721.00000. 00000. 79220+ 92861. 85721.92860. 92860. 00001.72080. 92862

3、. 85722. 64941. 85720. 00000. 92b6-1.79220. 79221. 72082. 6494O 0000-7. 5780n定义S矩阵中矩阵元(Si)的调和平均数的平方根为距离调和指数，用Tc表示。即:5. 50668, 29240.50.5=2.5057丿式中“工”是对分子中每个原子求和，例如乙酸甲醋的Tc值为:11/7.5780 + 1/4.5780 + 1/5.5066 + 1/8.2924 + 1/6.9546D1.堵锡华，陈艳，冯长君.连接性指数与多氯化萘理化性质的相关研究J.科中 技人学学报(自然科学版)，2006, 36 (11)： 124128.理

4、论与方法：考察文献所列75种多氯化萘(PCNs)分了的结构与其正辛醇/空气分配系数1gkoA和过 冷液体蒸汽压lg (pL/Pa)的变化规律，发现PCNs的这些性质与氯原了数日及在环上的位 置密切相关，为此，本文在文献方法基础上，提出一种新的原子点价值： = (4 - h )(X -2.25-2.25 )i i i C式中：0匚为第i个非氢原了点价，Xi为该原了的电负性(Pauling标度)；XC为碳原了的电 负性值；h.为与该原了直接连接的氢原了数。式中的(X - 2.25 )/(X - 2.25 )是对杂原了与i i C 碳原了的电负性值首先进行校正，在减去2.25 的基础上(通过反复试验

5、基础上得出的最佳 值)，得到杂原了相对于碳原了的电负性值，这一相对值能反映出连接在碳原了上的杂原了 的电了构象，也能蕴涵分了中原了的化学环境，对于萘环上连接的氯原了点价值，考虑到不 同位置的氯原了对分了的贡献有所不同，因此对式进行修正：=(4 h )(X 2.25 )(1 + 0.05 N- 2.25 )i i i i C式中N.表小蔡环上连接的各个氯原子的位置数。根据式可以计算得到连接在18位的氯 i原子的0Cl值分别为：0C11 = 12. 740, 0C12 = 13. 347, 0占13. 953, 0严.560, 0金15. 167， 0 =15. 773， 0= 16. 380，

6、0= 16. 987。Cl6Cl7Cl8EP5. 彭国文,肖方竹,聂长明.烷氧基氯硅烷类化合物气相色谱保留指数与结构关系研究J.南华大学学报（自然科学版），2007,21 （3）： 69.理论与方法：对于由n个顶点构成的分子图，其距离矩阵是一个nXn对称矩阵，本文提出的一组拓扑指数来源于距离矩阵，定义为：（、m与i相距m个边的顶点D = M ij 0其他定点其中 n 为分子隐氢图中的原子数。 由于化合物的性质因成键原子的支化度不同而不同。因此本文考虑到在距离 矩阵D中第一列增加各原子支化度参数Bi，Bi按下式定义：Bi=Zi-hi+1i i i式中Zi表示与分子中原子i的核外价电子数，hi表示

7、与原子i相连的氢原子数。 考虑到烷氧基氯硅烷化合物类的中心原子si的特殊性，其支化度定义为Bsi=Zsi-hi+3si si i 同时我们考虑用分子中原子 i 的平衡电负性对分子图的顶点着色,在距离矩阵 D 的第二列再增加各原子的平衡电负性的开方的倒数得：d 11仇e2血2.-J.1心其中e.= x 0.5, x.表示原子i的平衡电负性。然后对Q作如下处理：iimiM=Q X Q /mQm其中q /为Qm的转置矩阵，则新的拓扑指数NP定义为：匕mmNP =ln（久）m i，m其中久山为矩阵M的特征值之和。由于随着原子边间距的增加，原子的相互作用 越来越小，一般取m =1, 2, 3的一组拓扑指

8、数NP. NP2、NP3即能较好的表征 分子的结构信息，因此式中m取1，2, 3。例如甲烷氧基三氯硅烷的隐氢图及 增广矩阵为： 52S0田）7()00050636 40100 =70591 6)01020637 8()0100S0SSO 7)00仃-S0580 7()00080SS0 7)020250636 4()002仃a =70内6200020637 8()2()()80，泓）12()2()-S0倾7202020刚73()()()0()()()()0 637 8()()可得到甲烷氧基三氯硅烷的NP、NP2和NP3值分别为5.724 1，5.864 1，5.858 4。我们用上述方法对烷氧基

9、氯硅烷类化合物进行了计算。S1. 舒元梯. 氯代苯生物活性的分子拓扑研究 J. 达县师范高等专科学校学报 2006, 16 (5)： 4649.xph 0-125 lg+J)T3)3.98式(3)中，3.为第i原子的点价，xP为i原子的pauling电负性值，hi为i 原子直接键连的成o键的氢原子数目，o.为i原子成o键的电子数，n为i原子 ip 在形成离域n键时所提供的Pz电子数，3.98为氟元素的Pauling电负性。例:对氯甲苯(4-氯甲苯)(P-ClC6H4CH3)-CH3碳原子,=Cv, -CH=碳原子,-Cl原子的原点价分别为：x p h 0.125 1g (n + 1 )3(-C

10、H3) = P P i塩 p + 1 =2.554-30- 125-lgl /3.98=1.827833.983(=C)=2.553-0-lg2 /3.98=1.72923(-CH=)=2.553-1-lg2 /3.98=1.08853(-Cl)=3.161-0-lg3 /3.98=0.41511.王克强，杜仲祥. 硅烷摩尔体积与其分子结构定量关系的拓扑化学研究 J. 洛 阳师范学院学报， 2001， 2： 5658.理论与方法：摩尔体积与分了结构密切相关，其不仅与分了中基团的种类和数目有关，而且还与分了 中基团的连接性有关，为揭示摩尔体积与分了结构之间的定量关系，我们用图论方法对分了 结构进

11、行表征。根据基团贡献法对基团划分方式，我们将硅烷基团分为8种：“CH3-”、 “-CH2-”、“=CH-”、“=C=”、“SiH3-”、“-SiH2-”、“=SiH-”、“=Si=”，其中前 4 种为烷烃基团， 后4种为硅烷基团。分了中基团的连接性主要影响分了的形状，基团间的连接性不同，分了 形状也就有差别。根据文献的研究结果，距离矩阵D能够反映丰富的分了结构信息，我们 用距离矩阵D表征分了中基团的连接性，用染色矩阵C表征分了中基团的特性，用形状因 了矩阵S表征基团对分了形状的影响。设烷烃或烯烃的分了图G中有n个顶点，则距离矩 阵D= (d.)，则矩阵元d.表示分了图G中顶点(基团)i和顶点j

12、之间最短径上的边数; nxnij染色矩阵C= (c.)】，矩阵元C.是区别分了图G中不同顶点的染色因了，对不同的顶点，i 1xni我们用不同的染色因子加以区分；形状因子矩阵S=(s.)1，矩阵元s.是反映基团i对分子 i 1xni形状的影响程度。染色因子c和形状因子用计算机优化得到(具体方法见下文)。对于硅烷 分了中的8种基团，本文采用的染色因了为:“CH3-”,c = 3.1926、“-CH2-”,c = 1.8029、“=CH-”, c = 0.4606、 “=C=”， c=-0.9457;、 “SiH3-” ， c = 3.2699、 “-S.H2-” ， c = 1.9921、 “=S

13、iH-” c= 0. 9240、 “=Si=”， c= ，为简化起见，我们根据化合物的类别，确定形状因子，以减少形状 因子的数日，对于烷烃，本文取形状因了 s= 1.0000；硅烷形状因子1. 8359。通过对矩阵D、 C和S进行数学处理，可得到表征分子结构的拓扑指数W = CD E0P = SD E33式中，矩阵E=(1 1 11)T, Dk(k=0, 1, 2, n)为D的k阶矩阵，可通过算符Q作 用于矩阵D得到，具体方法见文献。既然分子结构可用拓扑指数W和P3表征，那么摩尔体积V必然为W和卩3的函数，3m3注意到摩尔体积V具有加和性，我们将摩尔体积V分解为基团的贡献和结构影响两部分, m

14、m即：V = f( W， P ) = f( W) + f( P3)m3式中第一项是在忽略分子结构影响条件下基团对摩尔体积的贡献；第二项是结构对摩尔体积 的影响。在系统分析实验数据的基础上，我们发现摩尔体积随着W的增大而增大，随着p3 的增大而减小。1.谢爱理，余训民.关于氯代苯在盐水溶液中的活度系数的研究J.上饶师范学 院学报，2000, 20 (6)： 4447 二 | n.l / (叭陥 |(4)式(4)中n*为盐中阳离子的最高有效主量子数，nX为盐中阴离子的最高主 量子数，Rmx为相联结的阴阳离子的有效离子半径之和，单位为纳米，vm和vx分别是M和X的价电子数。用式(4)计算复杂阴离子盐

15、时，nX和VX取阴离子 中各原子的平均值，例如计算K2SO4的Ym值时，对于so 2-离子，方法为：4nx=(ns+4no)/(l+4)=(3+4x2)/5=2.2Vp=(Vs+4Vo)/(l+4)=(6+4x6)/5=6X所以 K2SO4 的 Ym 值为：YM=(2x3.3227)/(2.2x(0.133+0.149)x6=4.37301.余训爽，余训民氯代芳烃化合物水溶解性和分配系数的QSAR研究J.科 技通报，2006，22 (1)： 49.1)m(Z - h )+ E (L x L ) Y i iiijin2(L +1)ii式(1)中n、m、L、Z.分别为分子隐氢结构图中顶点原子i的最

16、高主量子 数、成键的电子数、杂化方式和价电子数，h.为与i原子直接连接的氢原子数， L.为与原子i成键的相邻原子j的杂化方式。L.或的取值通过逼近法求得，对 于 sp3 杂化：则 L.=h.+1，L.=4(3- h.)/N.；对于 sp2 杂化：L.=Z.+ h.，L. =4(Z.- h.)/ Ni iji ji i i ji i j对于sp杂化：Li=Zi+ h.，L. =4(Z.-2 h.)/ N，N为顶点j原子在元素周期表中所处 的族数。2.余训爽，李克华，余训民.多氯代二苯并咲喃的lgKOW与结构信息连接性指数 的相关性研究J.安徽大学学报(自然科学版)，2007，31 (1)： 78

17、81.m(Z - h )+ S (p + h ) x (p + h )(2) T i i ii ij j( 2 )in2 ( p + h )iii式(2)中，n.、m.和Z.分别为分子图中顶点原子i的最高主量子数、成键的电 子数和价电子数；h.或h,为与i或j原子直接连接的氢原子数；p.或p,为顶点原子i或j成键时p轨道在杂化轨道中所占的比例，如在sp2化中p = 2，对于配位 3原子如Cl其p =1。(a)(b) (c)(d)图 2.1 子图类型根据分子图论学知识，上述子图分别为：(a)路径(m=0，1，2)；(b)簇(m=3)； (c)路径/簇(m=4)； (d)链(m=4)。分别用p、c

18、、pc和ch标在拓扑指数的右下角以示区别。例如：一氯代苯(C6H5C1)分子的隐氢图如下图2.2：Cl原子的点价鮎T=mi(Zj- hi)+Z(LixLj)/ni2(Li+1),其中成键的电子数m=1, 价电子数Z=7，与Cl原子直接连接的氢原子数h=0，则L=h+1=1(Cl原子为 iii isp3杂化)；午=0,顶点j原子在元素周期表中所处的族数N.=4,与Cl原子成键的3图 2.2 一氯代苯分子的隐氢图相邻原子j的杂化方式=4(彳勺)/件=4(4-0)/4=4,故： %t=1 x (7-0)+lx4/32(l+l)=0.6111 同理可求：1T=2T=3T=1.75； NT=j6T=i.

19、9583；叹丁=2.3429。123465一氯代苯分子连接性指数0Tp, 1Tp, 3TC, 4Tpc，6Tch的计算公式分别为： 叫=丫(罕)j-0-5=3x(1.75)-0-5+2x(1.9583)-0-5+(2.3429)-0.5+(0.6111)-0-5=5.6295 1T=Z(iTjT) -o.5=2x(1.75x1.75)-o.5+2x(1.75x1.9583)-o.5+2x(1.9583x2.3429)-o.5 +(0.6111)-0.5=4.02393Tc=(2.3429x1.9583x1.9583x0.6111)-0.5=0.4268 4Tpc=2x(2.3429x1.958

20、3x1.9583x0.6111x1.75)-0.5=0.6452 6Tch=(2.3429x1.9583x1.75x1.75x1.75X 1.9583)-0.5=0.14413. 余训爽.新的分子价连接性指数及其应用J.长江大学学报(自科版，2005, 2 (4)： 133135.理论与方法：在分析无机物分子的结构特点及组成原子的结构特征的基础上,提出了一个计算无机分子中原子特征值3 S的新方法：5 s = Im Z +m / N M 2 n2i i i i j i式中，Z.、m.和n.分别为无机物分子中原子i的价层电子数、成键电子数和价层最高主量子i i i数；N.为与无机分子中原子i相键合

21、的原子j的族数。可见，3s随n.的增大而减小，随乙jii.和m.的增大而增大，随N.的增大而减小。.j4. 杨道武，余训民. 新定义的价连接性指数与有机物理化活性的相关性研究 J. 有机化学， 2004， 5(24)：525536.理论与方法：实验结果表明，化合物的理化活性主要依赖分子的大小、形状和支化度等内在结构。分子的大小、形状和支化度在定量构效关系研究中一般用分子中的含碳数NC、拓扑指数T来Cp表征。拓扑指数T应用最广泛的是Randic的支化度指数mx及Kier等在mx基础上改进的分p子价连接性指数mxv。虽然mxv比mx使用范围更广，但5 v对主族元素中同一周期元素取相i 同值，因而导

22、致同碳数分子中出现简并现象，使得分子的结构选择性差，又因未考虑顶点原 子i的成键情况，如醛、酯、酸中羰基氧和醚基氧的5 v取值都为6,叔胺与腈中的氮的5 v均 ON 为5等，这显然是不合理的，因为分子结构不同。前不久杨锋等对5 v的取值进行了改进，i定义的价点价为5 -，但仍然不能区分同一氧化态不同杂化方式的原子价点价的计算问题，i如醛、醚中氧的5 -取值都为3。解决问题的关键仍是价点价5 v的取值问题，考虑到分子的 Oi性质主要由组成各原子的价层电子结构和成键环境等特性所决定，研究发现这种特性可用原 子的四个量子数和所在周期表中的族数来描述。事实上，角量子数反映了原子轨道的成键能 力，可用杂

23、化方式L或L来表征，磁量子数和自旋量子数可用成键电子数M.、孤电子对数ijip.来替代，原子的大小可用最高主量子数n.和族数N.表征，结合作者研究成果，在考虑顶 i i i点原子结构特征与相邻键合原子成键环境的情况下，对顶点原子价点价5 y （以示区别5 ,5 vi i i 和5 - ）重新定义如下：5 Y = In x Z x（Z h ）+（L x L ）/t（M h ）（L + p ）（n - 1）式中n，n，L； M； p分别为分子隐氢结构图中顶点原子i的最高主量子数、 i i i i i连接的顶点原子数、杂化方式、成键的电子数、未成键的孤电子对数，L.为与顶 点原子i成键的相邻顶点j原

24、子的杂化方式。九或L.的取值通过逼近法求得，对 于 sp3 杂化：则 L.=h.+1, L.=4（3- h.）/N.；对于 sp2 杂化：L =Z+ h，L. =4（Z- h.）/N.；i iji ji i i j i i j对于sp杂化：Li=Zi+ hi，L. =4（Zi-2 hi）/ N.，N为顶点j原子在元素周期表中所处 的族数。以价连接矩阵为基础，将Kier等改进的价连接性指数稍加变形，5 y构建 的一种能表征含多重键、杂原子化合物的新结构信息价连接性指数m3Y定义如 下：m3 Y=En（5 y ）0.55 y和m3 Y的物理意义与创新之处为：（1）m3 Y构建的基础是5 y，而5

25、y定义的基础是考虑了 i原子的价层电子结构特征（如 n.，N.，Z.，p. 等内因）和成键的化学环境（如hi， n.，M.，L.或 L. . .i . . .等外因），因而丫，13 Y蕴涵了分子的大小、形状、支化度和成键情况等内在 的结构要素。分子中氢原子越多则5 y越小，其m3 Y也就越小；反之hi越少则5 y 越大，其m3 Y也越大；对相同碳数NC的有机化合物，5 y越大则m3 Y也越大。 m3Y的大小不同，反映出化合物理化活性的差异性。因此6 Yi包含了丰富的原 子结构信息和成键性质，可以期待与有机化合物的理化活性具有较高的结构-性 质/活性的相关性。（2）m3 Y反映了分子的大小。含碳

26、原子数NC的多，成键数越多，分子体积 越大，则相应的m3Y就越大，因此m3Y具有加和性性质；m3Y反映了分子的 形状及支化度，同碳原子数NC的分子，分子的形状和支化度不同，相应的m3 Y 也不相同。（3）m3Y反映了分子的成键性质状况，具有高度的结构选择性和区分性。 对于同碳数NC结构相同的一元或二元取代基的链烃、芳烃，取代基位置虽不同， 但nx v（m=0，1）的值相同，而m3 Y （m=0，1）值则不同。如3-甲基庚烷与4-甲基 庚烷的 0x v 均为 6.4058，lx v 均为 3.8081,而 03 y 分别为 14.7930，14.7840,13 y 分别为27.2276, 27.

27、3111；如1,3-二氯苯和1,4-二氯苯的0x v和0X -均为5.5772,1X v和 1X -均为 3.8646，而 03Y分别为 17.6074, 17.5209,y分别为 36.9315, 36.9067。1.周莲.一个新的拓扑指数用于有机化合物的QSPR和QSAR研究J.西南民族 大学学报自然科学版，2004，30(6)： 727731.理论与方法：分子拓扑学就是在研究问题时仅考虑分子中原子间以何种形式的键相连，而 不涉及具体键长、键角等参数，用拓扑方法可得到分子的拓扑图(图 1)，若令 G=V, E是一个隐氢图，其中F是顶点集，E是边集。则邻接矩阵A=a., lj nxn 此处n

28、为顶点数，ajj为矩阵A的矩阵元。对于任意两个相邻顶点，产1，其他 情况下a. =0。距离矩阵D=d.，dij为拓扑距离，即沿顶点i和j之间的最短路ljlj nx n lj径所包含的边的数目。将矩阵A的任一行或列加和得顶点度v.，相应的矩阵为顶 i点度矩阵F。同理。将矩阵D的任一行或列加和得距离和s.，相应的矩阵为加和 i矩阵S。由F，S两个矩阵的矩阵元，我们定义一个新的拓扑指数ZH为： 鬲=(N/工叮订“式(1)中N为顶点数例如，2 一甲基丁烷的隐氢图如图1所示。邻接矩阵A和顶点度知阵F分别为:A =距离矩阵D及加和A、则 Z=5/(1/8+3/5+2/6+1/9+1/8)1/览1.9?5H

29、矩阵F，S均包含了丰富的分子结构37勾信息，特别是矩阵F与分子图具有1/1对应的关系，ZH指数结合了矩阵F和S的顶点度V.及距离和s.,因此可以期待 HiiZH指数具有较高的结构一性质括性相关性及结构选择性。Wiener在提出W指数的同时引出路径数的概念，所谓路径数（P“）是指两 原子间间隔若干个化学单键（即相距的边数州）的数目Dij就是距离矩阵中的矩 阵元。例如 A 分子的距离矩阵中， dij=1 的个数为 12， dij=2 的个数为 18， dij=3 的个数为12。由于州对每个原子均计算了两次，因此，路径数（Pn）应为其个 数的1/2： P1 =6，P2 =9，P3 =6。对于A分子的州=4个数为零，故P4 =0。据文 献报道，P2由Godell和Scantlebury于1964年提出，用以反映分子支化度;P3由 Wiener提出，称为极化数，用以表征分子的形状;P4是由我国学者王化云引入的.