碳管笼烯的稳定性问题

碳笼烯C_(60)及C_(70)等的发现启示人们去探索合成更多的新型封闭型碳烯化合物,可以设想,由六元碳环构成的石墨层平面可以卷成顶和底为p元环、侧面为紧接着顶底的p个五元环、中间k层p个六元环围成的碳管笼烯C_n(n=2p(k+2)),其结构见图1。笼烯的每个C原子都与3个相邻的C原子成键,并有一个与笼面垂直的单电子占据轨道,可以相互作用构成笼面共轭大π键。我们在     

三金属氮化物富勒烯中富勒烯笼和内嵌团簇的尺寸效应

内嵌金属富勒烯是指金属原子或其团簇内嵌于碳笼之中而形成的复合结构.实验和理论研究都表明内嵌金属原子或团簇对碳笼的电荷转移作用在这类分子中是普遍存在的,并且电荷转移作用能够稳定化原本高度活性的碳笼.三金属氮化物富勒烯(trimetallic nitride fullerenes,TNFs)因其新奇的结构和性质吸引了广泛的研究兴趣.研究者已经对Sc3N@C68,Sc3N@C78和Sc3N@C80进行了广泛而深入的研究并阐明了它们的结构;而且,三金属氮化物团簇(M3N,M=Sc,Y,La等)在各种碳中的稳定化机理已经被阐明,即能够嵌入氮化物团簇的碳笼是那些LUMO+3与LUMO+4能级能量差大的碳笼.对于含有更重的内嵌团簇的TNFs,目前已经取得相当的研究进展,然而仍然有很多未解之谜等待揭开.如:为什么烟灰中Sc3N@C2n,Y3N@C2n和La3N@C2n(n=34,39和40)的产量戏剧性地下降?三金属氮化物M3N@C2n中,电荷转移相互作用像单金属内嵌富勒烯M@Cn中的情况一样吗?为了揭开上述未解之谜,我们对M3N@C80进行系统研究的基础上(ChemPhysChem2006,7:1306),通过密度函数理论方法,对M3N@C2n(n=34,39,40)进行了更广泛更深入的研究.结果表明,对于C68-6140,金属原子位于五元环-五元环共边(B55)的前面;对于C78-5,金属原子位于六元环-六元环共边(B66)的前面;对于C80-7,Sc原子位置不固定,而Y、La原子位于六元环(R6)的前面.Sc3N在所有考虑的碳笼内都是平面的,Y3N在C68和C78中是微弱锥化的,然而在C80-Ih中是平面的;相反,La3N在所有碳笼中都是锥化的.总之,内嵌团簇趋向于在碳笼中形成具有最大金属-金属(M-M)距离的结构,以便于张力能得到最大限度的释放.内嵌锥形团簇的那些碳笼与相应的富勒烯笼比较起来,结构上发生了重大的变形,此种情况下相应分子中存在重要的M-C以及M-M排斥力,相应的分子是不稳定的.对优化后的M3N@C2n的结构分析表明,Y-Y/Y-N和La-La/La-N距离都分别小于自由态Y3N、La3N中的M-M/M-N距离.这些结构中,金属原子与碳笼的成键特性或内嵌团簇与碳笼之间的分子轨道相互作用都是次要的.对于较小的碳笼和较大的内嵌团簇形成的体系,首先满足的应该是几何尺寸上的要求.这些结果清楚地表明,内嵌团簇尺寸较大的三金属氮化物富勒烯中存在重大的M-M以及M-C排斥相互作用并随着碳笼的增大而减小.这些计算结果揭示了内嵌团簇和碳笼的尺寸对TNFs的结构和稳定性的影响,能够良好地解释M3N在碳笼中的位置以及M3N@C2n在产量上的差异并能够为这类物质的结构测定、结构-性质关系阐释提供重要的参考价值.     

N8H8环状异构体的结构与稳定性的理论研究

采用密度泛函理论的B3LYP方法在6-311 G**基组水平上对N8H8氮氢环状化合物可能存在的构型进行了几何优化,得到74种稳定异构体,应朋A然键轨道理论NBO和分子中的原子理论AIM分析了这些化合物成键特征和相对稳定性,G3MP2方法汁算了各异构体的能量及生成热.研究结果表明:N原子孤对电子到相邻的氮氮键的超共轭作用是影响氮氮键长变化的丰要因素;N8H8环状异构体的稳定性顺序为:六元环>七元环>八元环,五元环>三元环>四元环,六元环是这些N8H8环状异构体中最稳定的,最不稳定的是四元环,G19是所有环状异构体中能量最低的:M3能量最高,稳定性最差,A7密度最大.     

(BCO)12与(CH)12稳定性的环张力分析

基于密度泛函理论的B3LYP方法, 对具有等瓣相似性的(BCO)12和(CH)12的10种异构体结构的稳定性进行了计算对比研究, 这10种异构体由三元、四元、五元和六元环组成. 环张力分析表明对羰基硼笼体系, 三元环起主要的稳定化作用, 而四元环是张力的主要来源, 对碳氢笼体系, 五元环起主要的稳定化作用. 电子差分密度表明羰基硼笼中的三元环与碳氢笼中的三元环有不同的电子结构, 导致了它们不同的张力表现. 核独立化学位移(NICS)分析表明, 尽管σ芳香性不是稳定性的决定因素, 但对笼的稳定性有一定的影响.     

B24N24团簇的结构与稳定性

采用密度泛函理论,在B3LYP/6-31G^*水平下,对B₂₄N₂₄笼状团簇的12种异构体进行了优化,并对它们的几何构型、化学键性质、振动光谱和稳定性进行了探讨.研究表明:具有S₈对称的含有2个八元环、8个四元环和16个六元环的结构h是B₂₄N₂₄笼状团簇最稳定的异构体,只存在B-N键,而无N-N和B-B键.含有五元环结构的稳定性最低.B-B和N-N键对的数目越多,结构的稳定性越低.12种异构体的稳定性顺序为h＞a＞b＞I＞g＞l＞c＞k＞j＞d＞e＞f.     

(XY)12(X=B,Al;Y=P)团簇的结构与稳定性

采用B3LYP／6031G*方法,对(XY)12(X=B,Al;Y=P)笼状团簇的同分异构体进行优化,筛选出能量最低的构型．讨论它们的几何构型、HOMO-LUMO能隙、生成焓、核独立化学位移(NICS)和自由能．得到(BP)12和(AlP)12团簇的最稳定构型均为具有Th对称性的四、六元环组成的笼,亚稳态结构中含有五元环。     

碱土金属Mg外掺杂硼笼B40Mg6的储氢性能

对外掺杂碱土金属Mg的B₄₀硼笼的储氢性能进行了密度泛函理论研究. 结果表明, B₄₀笼含有2个六元环与4个七元环. 多个Mg原子对B₄₀笼进行外掺杂时不会发生成簇现象, 有利于进一步储氢. Mg原子外掺杂的B₄₀笼对H₂分子的平均吸附能介于物理吸附与化学吸附之间(0.1~0.8 eV). 体系的储氢密度达到7.60%(质量分数). 储氢结构能在常温常压下释放H₂分子, 因此, Mg原子外掺杂的B₄₀笼比Mg合金具有更好的储氢性能.     

叠氮二氢硼多聚体结构和性质的理论研究(英文)

本文采用DFT-B3LYP方法,以不同基组对叠氮二氢硼多聚体(H2BN3)n (n=1-4)进行计算研究.二聚体(H2BN3)2(C2h对称性)中含B2N2平面四元环结构.船式(Cs对称性)和椅式(C3v对称性)三聚体(H2BN3)3的结合能相近(-122 和 -126 kJ·mol-1),其中均含B3N3六元环结构.拥有B4N4八元环结构的四个四聚体的结合能只有稍微差别.与单体相比,簇合物的结构参数变化较大.由ΔG0T可知,298.2 K下单体形成二聚体在热力学上是不利的,而形成三聚体和四聚体是有利的.     

(XY)12 (X═B, Al; Y═P)团簇的结构与稳定性

采用B3LYP/6-31G*方法, 对(XY)₁₂ (X═B, Al; Y═P)笼状团簇的同分异构体进行优化, 筛选出能量最低的构型. 讨论它们的几何构型、HOMO-LUMO能隙、生成焓、核独立化学位移(NICS)和自由能. 得到(BP)₁₂和(AlP)₁₂团簇的最稳定构型均为具有T_h对称性的四、六元环组成的笼, 亚稳态结构中含有五元环.     

N8H8环状异构体的结构与稳定性的理论研究

采用密度泛函理论的B3LYP方法在6-311＋＋G**基组水平上对N₈H₈氮氢环状化合物可能存在的构型进行了几何优化, 得到74种稳定异构体, 应用自然键轨道理论NBO和分子中的原子理论AIM分析了这些化合物成键特征和相对稳定性, G3MP2方法计算了各异构体的能量及生成热. 研究结果表明: N原子孤对电子到相邻的氮氮键的超共轭作用是影响氮氮键长变化的主要因素; N₈H₈环状异构体的稳定性顺序为: 六元环＞七元环＞八元环, 五元环＞三元环＞四元环, 六元环是这些N₈H₈环状异构体中最稳定的, 最不稳定的是四元环, G19是所有环状异构体中能量最低的; M3能量最高, 稳定性最差, A7密度最大.     

Aromaticity of Heterofullerenes C18BxNy (x+y=2) and Their Molecular Ions

The aromaticity of all possible substituted fullerene isomers of C18N2, C18B2, C18BN, and their molecularions which originate from the C20 (Ih) cage were studied by the topological resonance energy (TRE) and the percentage topological resonance energy methods. The relationship between the aromaticity of C18BxNy isomers and the sites where the heteroatoms dope at the C20 (Ih) cage is discussed. Calculation results show that at the neutral and cationic states all the isomers are predicted to be antiaromatic with negative TREs, but their polyvalent anions are predicted to be aromatic with positive TREs. The most stable isomer is formed by heteroatom doping at the 1,11-sites in C18N2, C18B2, and C18BN. Heterofullerenes are more aromatic than C20. The stability order in the neutral states is C18N2>C18BN>C18B2>C20. The stability order in closed-shell is C18B2 8->C20 6->C18BN6->C18N2 4-. This predicts theoretically that their polyvalent anions have high aromaticity.     

[B20H18]n-(n=0,2,4,6)离子异构体的结构、稳定性与电子计数规则

对单个多面体硼烷(Polyhedral boranes)已有较多的理论研究，由多个多面体通过共用一个或多个顶点而构成的稠合型硼烷(Macropolyhedral borones)，具有多种多样的结构类型，并已被大量合成出来，目前对稠合型硼烷结构及成键特性的理论研究尚不充分，它们不能再以简单的closo，nido和arachno分类，Wade规则也不再适于解释其结构，理论上各种电子计数规则已有不少报道，对预言和发现新的分子十分重要。     

星际尘埃HC_(3-x)NC_x(x=0～3)与HC_(5-x)NC_x(x=0～5)分子的DFT研究

运用密度泛函理论（DFT）中的B3LYP方法,计算了星际尘埃HC_3N,HC_5N分子 及它们的同分异构体的分子结构、相对稳定性、偶极矩、转动常数等,对于HC_(3- x)NC_x(x = 0～3),HC_3N分子比其它的异构体要稳定,将计算所得的HC_3N, HC_2NC的转动常数B_0与实际的观测值B_0相比较,发现两者能较好地吻合。比较 HC_(5-x)NC_x(x = 0～5)的各同分异构体的总能量,发现HC_5N的相对稳定性最高 ,但它的偶极矩相对较低。将HC_5N的转动常数B_0的理论计算值与实际观测值进行 比较,其相对标准误差仅0.6%。     

Ab initio studies on isomers of macropolyhedral borane ions [B20H18]n (n = 0, -2, -4)

Various isomers of macropolyhedral borane ions [B20H18]n (n = 0, -2, -4) are investigated by using the density functional theory methods at RB3LYP/6-31+G* and RB3LYP/6-31G* levels to obtain the optimized geometries, harmonic vibrational frequencies, electron structures, and the stability order. The calculated results show that optimized bond lengths are consistent with the available experimental values and the natural populations, taking [a2 -B20H18]4- (4) as an example, are also in agreement with NMR spectra. The calculated vibrational frequencies are all real, so all of these isomers could be stable, among which [a2 -B20H18]2- (3) and [a2 -B20H18]0 (7) are considered for the first time in this paper. On the basis of the contour maps of molecular orbitals, the delocalized characteristic of molecular orbitals and the possible redox mechanism of these ions are also discussed. Moreover, the analysis on counting of skeletal bonding electrons shows that the isomers (1)-(6) obey the electronic requirement predicted by the mno rule, whereas the newly predicted isomer (7) does not match the mno rule.     

A theoretical study of the structures, energetics, stabilities, reactivities, and out-of-plane distortive tendencies of skeletally substituted benzenes (CH)(5)XH and (CH)(4)(XH)(2) (X = B(-), N(+), Al(-), Si, P(+), Ga(-), Ge, and As(+)).

Ab initio molecular orbital theory at Hartree-Fock (HF), post-Hartree-Fock (MP2 and CCSD(T)), and the hybrid density functional theory (B3LYP) calculations were done on the mono-(CH)(5)XH and diskeletally substituted (CH)(4)(XH)(2) benzenes (X = B(-), N(+), Al(-), Si, P(+), Ga(-), Ge, and As(+)). The computed relative energies of the disubstituted isomers show interesting trends. While the ortho-isomer is the most stable for X = Ga(-), Ge, and As(+), meta was found to be the most stable for X = B(-), N(+), Al(-) and Si, and para was found to be the most stable for X = P(+). Various intricate factors that govern the relative stabilities, such as the sum of bond strengths in the twin Kekule forms, rule of topological charge stabilization (TCS), and electrostatic repulsion were critically examined. The sum of bond strengths in the twin Kekule forms was proved to be quite a successful measure in predicting the relative stability orders between ortho- and meta-/para-isomers. The rule of TCS breaks down especially in the presence of overwhelming factors such as the differences in the cumulative bond strengths of the two positional isomers; however, the stability ordering between the para- and meta-isomers is successfully predicted in most cases. The tendency for ring puckering increases a great deal especially when the substituents are from 3rd or 4th row. Extension of the popular inverse relationship between the thermodynamic stability and reactivity was found to be inapplicable for this class of compounds. The computed singlet-triplet energy differences and the chemical hardness (eta) values indicate that the skeletal substitution weakens the pi-strength of the benzenoid system and increases their reactivity.     

The distribution pattern of squares on the stability of F4F6-(BN)n (n = 10–22) polyhedrons

To gain an insight into the structures and stability of F₄F₆-(BN)_n polyhedrons with alternation of B and N atoms, a density functional theory study was performed on all isomers of F₄F₆-(BN)_n polyhedrons with n between 10 and 22. The calculation results demonstrate that the lowest energy isomers do not contain B₄₄ bonds (the bonds shared by two squares) and the energies of those isomers containing B₄₄ bonds increase with the number of B₄₄ bonds linearly, indicating that the energetically favored structures of F₄F₆-(BN)_n polyhedrons satisfy the isolated square rule and square adjacency penalty rule. The structural analysis reveals that the stability is determined by the pyramidalization of B and N atoms at the square–square fusion. The binding energy is fitted to the numbers of edges and a model is proposed for predicting the relative stability of these B–N polyhedral molecules.     

杂富勒烯的电子结构与性质关系的研究──Ⅱ.C_(57)X_3,C_(56)X_4(X=B、N、P)的位置异构体

用 CNDO／2方法在 586微机上计算了 C57X3C56X4（X=B、N、P）的 234个位置异构体的电子结构。在C57X3（X=B、N、P）位置异构体中，C57X3（1，2，9）（X=B、N、P）分别是最稳定的。对于C56X4（X=B、N、P）位置异构体，C56B4（1，2，9，8），C56N4（1，2，9，12）和C56P4（1，2，9，12）分别是最稳定的，但稳定性都比C60差。其氧化或还原性都比 C60好，将它们和 C60比较，与 X相距一个或两个键的 C原子电荷密度增加或减少较多，其余电或亲核反应性增加；X与 X，X与 C原子之间 Wiberg Order都减少较多，其键的强度削弱；邻近杂原子的 C与 C原子之间 Wiberg Order减少或增加很少，其键的强度稍有削弱或增强。     

Aromaticity of All Possible C26N2 Isomers

The aromaticity of all possible heterofullerenes C₂₆N₂ and C₂₈ based on T_d symmetry has been studied by means of the topological resonance energy and percentage topological resonance energy methods. The relationship between the aromaticity of the C₂₆N₂ isomers and the sites where nitrogen atoms dope at the C₂₈ cage has been discussed. The calculation results show that the most stable isomer of C₂₆N₂ derivatives is formed by nitrogen atoms doping at the two tetrahedral vertices. C₂₆N₂ isomers are more stable than C₂₈, but the C₂₆N₂^2? isomers are less stable than C₂₈ ^4?4. The effect of nitrogen substitution on C₂₈ stability was investigated by the topological charge stabilization rule.     

硼/氮掺杂富勒烯C20芳香性的争论

用拓扑共振能(TRE)和百分拓扑共振能(%TRE)方法对硼和氮掺杂富勒烯C20的芳香性进行了研究, 并与核独立化学位移(NICS)方法的结果进行了比较. 研究结果表明, 在C20-2nX2n(X=B, N; n=1,2,3,4)中只有C16N4, C14N6和C12N8具有芳香性, 而其它化合物都具有反芳香性. 这与NICS方法的预测结果不一致. 对NICS判据和TRE方法对C20和C20-2nX2n(B, N; n=1,2,3,4)的芳香性进行判断时所得出的不一致结果及其原因进行了讨论, 认为用NICS判据和2(n+1)2规则对这些化合物芳香性的预测是不可靠的.     

The effects of atom pyramidalization and square distribution on the stability of F4F6‐(BN)n polyhedrons

The structures and stability of F₄F₆‐(BN)_n polyhedrons (n = 20–30) with the alternation of B and N atoms were studied with DFT method. The calculation results reveal that the atoms at square–square fusions with large pyramidalization angles are remarkably extruded out of the surfaces of (BN)_n polyhedrons. The energetically favorable isomers do not contain square–square bonds and the energies of those isomers containing square–square bonds increase with the number of square–square bonds linearly, demonstrating that the energetically favorable structures of F₄F₆‐(BN)_n polyhedrons satisfy the isolated square rule and square adjacency penalty rule. The atom pyramidalization determines the stability of the isomers. The binding energy is fitted to the numbers of vertices formed from different faces and a model is proposed to predict the relative stability of these polyhedral molecules. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011