XeO$_{3}$和XeOF$_{2}$与三种不同杂化方式N原子的Xe$cdots$N非共价作用的理论研究

本文采用密度泛函理论和二级微扰理论方法研究了XeO$_3$和XeOF$_2$与一系列不同杂化含氮给体配合物之间的相互作用. 含氮给体与XeO$_3$和XeOF$_2$之间的相互作用能在6.5$sim$19.9 kcal/mol之间. 相互作用强弱顺序为sp$^3$$>$sp$^2$$>$sp,且XeO$_3$比XeOF$_2$高. 约化密度梯度图证明,对于sp$^2$和sp$^3$杂化的给体,空间位阻效应在相互作用中起着一定的影响. 能量成分分析发现,Xe$cdots$N作用主要的吸引部分是静电相互作用. 在XeO$_3$配合物中,极化能的占比大于色散,而在XeOF$_2$配合物中则相反. 除表面静电势($V_mathrm{s}$)外,其他5个相互作用参数之间的线性相关系数均较大.     

共轭体系分子轨道能的共振同系规律

同系线性规律已有较多研究.最近戴乾圜等用共振论阐释了同系线性规则,提出了共振同系因子,并用式(1)计算了电子吸收光谱数据的线性相关性.     

键能的分子轨道理论研究Ⅱ.键能与Mulliken布居对键强度的判断的比较

用键能EAB和Mulliken布居对化学键强度的判别进行了分析比较.结果表明,键能判据比Mulliken布居判据所得结论更符合实际情况.作为衡量原子间化学键强度的尺度,不仅应考虑原子轨道间的布居因素,还应考虑分子轨道(或原子轨道)的能量因素.     

键能的分子轨道理论研究2: 键能与Mulliken布居对键强度的 判断的比较

用键能E~A~B和Mulliken布居对化学键强度的判别进行了分析比较。结果表明,键能判据比Mulliken布居判据所得结论更符合实际情况。作为衡量原子间化学键强度的尺度,不仅应考虑原子轨道间的布居因素,还应考虑分子轨道(或原子轨道)的能量因素。     

键能的分子轨道理论研究 1: 理论公式

从LCAO-MO出发, 给出了一个计算键能的近似方法, 即EAB(i)-∑∑CaiSabCbiεi为第i个占据分子轨道(MO)中的一对电子对A-B键键能的贡献。对所有分子轨道求和即为该键的键能: EAB=∑EAB(i)。按该方法, 不仅可以计算各种不同分子中每两个相键连原子间的键能, 还可以从MO及AO角度分析每一具体键, 如σ, π, δ键的键能以及各AO对键能的贡献。该方法虽有别于求键焓和平衡离解能De, 但计算结果和De的实验值甚相符合。通过对键能的分析研究, 能较好地揭示原子间的相互作用关系及化学键的强弱, 从而可进一步探讨化学反应活性, 反应速率等化学性质。     

键能的分子轨道理论研究2: 键能与Mulliken布居对键强度的 判断的比较

用键能E~A~B和Mulliken布居对化学键强度的判别进行了分析比较。结果表明,键能判据比Mulliken布居判据所得结论更符合实际情况。作为衡量原子间化学键强度的尺度,不仅应考虑原子轨道间的布居因素,还应考虑分子轨道(或原子轨道)的能量因素。     

键能的分子轨道理论研究 1: 理论公式

从LCAO-MO出发, 给出了一个计算键能的近似方法, 即EAB(i)-∑∑CaiSabCbiεi为第i个占据分子轨道(MO)中的一对电子对A-B键键能的贡献。对所有分子轨道求和即为该键的键能: EAB=∑EAB(i)。按该方法, 不仅可以计算各种不同分子中每两个相键连原子间的键能, 还可以从MO及AO角度分析每一具体键, 如σ, π, δ键的键能以及各AO对键能的贡献。该方法虽有别于求键焓和平衡离解能De, 但计算结果和De的实验值甚相符合。通过对键能的分析研究, 能较好地揭示原子间的相互作用关系及化学键的强弱, 从而可进一步探讨化学反应活性, 反应速率等化学性质。     

杂苯C5H5X（X=N,P,As,Sb,Bi）芳香性的核独立化学位移（NICS）与异构体稳定化能（ISE）研究

应用量子化学方法,通过核独立化学位移(NICS)和异构体稳定化能(ISE)的计算,研究了苯及第五主族元素取代杂苯分子C5H5X(X=N,P,As,Sb,Bi)的芳香性与稳定性.局域轨道定位函数局部最大值的计算结果表明,分子中C—X成键强度与实验稳定性顺序一致.从头算与密度泛函理论对分子的化学位移计算结果各异,计算值与实验值相关分析表明,Hartree-Fock方法对所研究体系的NICS比密度泛函理论具有更好的相关性.在分子环平面上方0.8～0.9处的NICS是芳香性判据的最佳选择,由自然定域分子轨道分解NICS最大处的zz张量值,结果显示π键对分子的芳香性起主要贡献.异构体稳定化能与NICS(max)的zz张量及π键(NICS(max)zzπ)均有很好的相关性,可以表征杂苯分子C5H5X全局芳香性,其顺序为:苯>吡啶>磷杂苯>砷杂苯>锑杂苯>铋杂苯.特别地,对这类分子π轨道的研究发现不包含X原子的π轨道将产生异常大的π键芳香性,这一现象可为分子磁性设计提供理论指导.     

dSH扭曲程序的建立以及在两种基组(NBO, LFMO)下对扭曲分子的能量分解

介绍了d_SH 扭曲法产生的背景, 为d_SH扭曲法研究Mills-Nixon效应提供了一个自动化程序和方法. 通过在程序中引入随机函数模拟并代替了人工调节构型参数的过程, 提高了效率并拓宽了程序的应用范围. 利用此程序算得三元苯并分子C₆B₃H₃ 和C₆O₃的一系列的扭曲结构, 并利用自然键轨道(NBO)和定域片断分子轨道(LFMO)两种基组分别进行了NBO能量分析和Morokuma SCF能量分解, 并以动态的方式分析了扭曲过程中垂直共振能及其各分量的变化趋势, 比较了不同基组和不同分析方法下的扭曲的驱动力和阻力. 由分析知, 由于NBO基组的非绝对定域性和NBO能量分析方法的一次性对角化直接导致NBO基组及其能量分析方法在d_SH扭曲能量分析中的不合理.     

氨基酸侧链与氧化鸟嘌呤碱基对体系的氢键作用

应用量子化学方法,分别在气相和水溶液中对氨基酸侧链与氧化鸟嘌呤碱基对(8-oxo-G∶C)形成的三体复合物的氢键键能、几何结构、电荷分布及二阶稳定化能进行了研究.结果表明,水溶液的存在削弱了复合物中的氢键强度,电荷分布变化明显,水溶液中形成氢键位点的电荷变化量约为气相中的10倍,而几何结构变化不明显、对于酶与DNA之间的相互作用的研究需在水溶液中进行.水溶液对带电三体复合物中8-oxo-G∶C与氨基酸侧链间的氢键有较大影响,键能平均减小了69.23 k J/mol,不带电复合物仅减小了3.60k J/mol.水溶液中三体复合物中8-oxo-G∶C间的氢键受侧链的影响不大,且与侧链带电与否无关,带电复合物和不带电复合物的氢键强度分别减小了24.57和30.05 k J/mol,且二阶稳定化能越大,其对应的氢键键长越短.     

dSH扭曲程序的建立以及在两种基组(NBO, LFMO)下对扭曲分子的能量分解

介绍了d_SH 扭曲法产生的背景, 为d_SH扭曲法研究Mills-Nixon效应提供了一个自动化程序和方法. 通过在程序中引入随机函数模拟并代替了人工调节构型参数的过程, 提高了效率并拓宽了程序的应用范围. 利用此程序算得三元苯并分子C₆B₃H₃ 和C₆O₃的一系列的扭曲结构, 并利用自然键轨道(NBO)和定域片断分子轨道(LFMO)两种基组分别进行了NBO能量分析和Morokuma SCF能量分解, 并以动态的方式分析了扭曲过程中垂直共振能及其各分量的变化趋势, 比较了不同基组和不同分析方法下的扭曲的驱动力和阻力. 由分析知, 由于NBO基组的非绝对定域性和NBO能量分析方法的一次性对角化直接导致NBO基组及其能量分析方法在d_SH扭曲能量分析中的不合理.     

Assessing stereoelectronic effects in dipolar cycloadditions yielding fused thiazolopyridone rings

We report a combined experimental and computational study on the cycloadditions of bicyclic 1,3-thiazolium-4-olates, derived from thiazolidin-2-thiones, with asymmetrically-substituted acetylenes. These results provide further mechanistic insights into the above dipolar cycloadditions and enable an unequivocal characterization by NMR spectroscopy of regiochemical patterns as previous derivatives had substituents at both C-2 (in the dipole) and C-6 (in products). Accordingly, new dihydrothiazolopyrid-2-ones have been obtained from a thioisomünchnone lacking substitution at C-2. With the aim of assessing the steric hindrance as well as the facial stereoselection induced by a bulky group on the Si face (relative to C-7a) of the mesoionic heterocycle, a chiral thioisomünchnone has also been obtained along with the resulting optically active thiazolopyridones. A computational study of these particular cycloadditions, largely based on an NBO analysis, allowed us to evaluate the influence of substituents on intermolecular steric repulsions, charge transfers, as well as solvent effects.     

Theoretical study of activation Fe–O bond of FeO by CO in the gas phase

The entire reaction mechanism for the gas phase CO–CO₂ conversion by FeO⁺ is discussed by means of the density functional theory and the intrinsic reaction coordinate approach. The calculated results have strongly indicated that the reaction of is a spin-forbidden reaction between the quartet and the sextet potential energy surfaces (PESs). There is a crossing point between the quartet and the sextet potential energy surfaces which may play a significant role in this reaction, by which the activation energy can be decreased from −15.1 to −56.4 kJ mol⁻¹ at the reaction system.     

Theoretical investigation of C–HM interactions in organometallic complexes: A natural bond orbital (NBO) study

The nature of C–HM agostic interactions in model metal complexes [M²⁺(CH₂CH₃)(PH₃)_nCl] (where M = Sc, Ti, V, Ti, Cr, Mn, Fe, Co, Ni, Cu, Zn; n = 1, 2, 3, 4) was studied with the natural bond orbital analysis (NBO) approach using density functional theory (DFT) optimized geometries at the B3LYP/6-31G(d,p) level of theory. The effect of nature of metal, coordination number, oxidation state and ligand field effects on the agostic interaction is examined. A set of 20 crystal structures of organometallic complexes taken from the Cambridge Structural Database (CSD) was studied computationally employing AIM theory and NBO analysis, and the applicability of these methods was critically accessed in demarcating the two types of interaction.     

分子内结构环境对解离能的影响

为了理解化学键的这一结构效应, 本文对具有相同化学键而分子内结构环境不同的系列分子进行了计算研究, 讨论了化学键结构环境对解离能的影响.     

[EPy][BF4]和[EPy][PF6]离子对的气相及液相阴阳离子相互作用研究

采用从头算HF/6-31G和密度泛函理论B3LYP/6-31+G(d,p)方法, 对乙基吡啶四氟硼酸盐([EPy][BF4])和乙基吡啶六氟磷酸盐([EPy][PF6])的离子对进行了结构优化和频率分析, 并利用自洽反应场(SCRF)的导体极化连续模型(CPCM)考察了离子对液态下的结构及相互作用, 得到了两种离子对的红外光谱及气相、液相下最稳定结构. 由两种离子对的几何参数可知, 阴阳离子通过氢键相互作用, 两种离子液体的红外光谱特征值与实验值比较吻合. 应用自然键轨道(NBO)理论分析了吡啶阳离子及离子对中的原子电荷分布和电荷转移情况, 结果表明两种离子对中阴阳离子间存在静电相互作用和氢键作用. 通过几何参数、相互作用能及NBO分析研究发现, 液相下由于周围电荷的中和作用, 离子对中阴阳离子的相互作用明显降低.     

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

A series of 1-D zinc(II) coordination polymers with 3-D supramolecular networks: synthesis,structural investigation,and NBO analysis

Abstract

N,N′-Bis(pyridin-4-ylmethylene)naphthalene-1,5-diamine (L) acts as a bipyridine analogue linker ligand towards {Zn₇(μ₄-O)₂(OAc)₁₀}, {Zn₂(NCS)₂(OAc)₂}, and {Zn(N₃)₂} nodes and allows construction of three new 1-D coordination polymers, the linear chain [Zn₇(μ₄-O)₂(OAc)₁₀(L)]_n (1), [Zn(NCS)(OAc)(L)]_n (2) in ladder-type geometry and the zigzag chain [Zn(N₃)₂(L)]_n (3). Structural characterization reveals that in 1 acetate anionic ligands connect seven Zn(II) ions through the bridging coordination modes μ₃-η¹,η² and μ₂-η¹,η¹. The resulting heptanuclear node is located on an inversion center and therefore consists of four crystallographically distinct cations; their coordination spheres correspond to distorted octahedra or tetrahedra. The Zn(II) ions in polymer 2 exhibit distorted trigonal bipyramidal {ZnN₃O₂} coordination; μ₂-η¹,η¹ coordinated acetate and terminal thiocyanate ligands lead to inversion-symmetric [Zn₂(NCS)₂(OAc)₂] secondary building units (SBU), which are further linked by the N,N′-bipyridine analogue L. Terminal coordination of two anionic azide ligands and the bridging bipyridine L result in coordination polymer 3, in which the cations adopt distorted tetrahedral {ZnN₄} coordination. In all crystalline solids 1–3, adjacent 1-D chains interact through π–π stacking and non-classical (C???H···O, C???H···π) hydrogen bonds, leading to 3-D supramolecular architectures. Differences in their 3-D arrangement are due to variations in the anionic co-ligands, subtle conformational differences in the semi-rigid linker and the variable coordination sphere about the zinc cations. Thermogravimetric investigations indicate differences in both thermal stability and decomposition mode. Natural bond orbital (NBO) analysis provides a convenient basis for investigating the intramolecular bonding interactions and delocalization effects in these molecular systems. Finally, solids 1–3 exhibit intense luminescence at room temperature.     

Kinetic Study and NBO Analysis of the Dehydrogenation Mechanism of Five‐membered Ring Heterocyclic 2,5‐Dihydro‐[furan,thiophene, selenophene

The theoretical study of the dehydrogenation of 2,5‐dihydro‐[furan ( 1 ), thiophene ( 2 ), and selenophene ( 3 )] was carried out using ab initio molecular orbital (MO) and density functional theory (DFT) methods at the B3LYP/6‐311G**//B3LYP/6‐311G** and MP2/6‐311G**//B3LYP/6‐311G** levels of theory. Among the used methods in this study, the obtained results show that B3LYP/6‐311G** method is in good agreement with the available experimental values. Based on the optimized ground state geometries using B3LYP/6‐311G** method, the natural bond orbital (NBO) analysis of donor‐acceptor (bond‐antibond) interactions revealed that the stabilization energies associated with the electronic delocalization from non‐bonding lone‐pair orbitals [LP(e)_X3] to δ*_{C(1)  H(2)} antibonding orbital, decrease from compounds 1 to 3 . The LP(e)_X3→δ*_{C(1)  H(2)} resonance energies for compounds 1 – 3 are 23.37, 16.05 and 12.46 kJ/mol, respectively. Also, the LP(e)_X3→δ*_{C(1)  H(2)} delocalizations could fairly explain the decrease of occupancies of LP(e)_X3 non‐bonding orbitals in ring of compounds 1 – 3 ( 3 > 2 > 1 ). The electronic delocalization from LP(e)_X3 non‐bonding orbitals to δ*_{C(1)  H(2)} antibonding orbital increases the ground state structure stability, Therefore, the decrease of LP(e)_X3→δ*_{C(1)  H(2)} delocalizations could fairly explain the kinetic of the dehydrogenation reactions of compounds 1 – 3 (k ₁ >k ₂ >k ₃ ). Also, the donor‐acceptor interactions, as obtained from NBO analysis, revealed that the (_C(4)C(7)→δ*_{C(1)  H(2)} resonance energies decrease from compounds 1 to 3 . Further, the results showed that the energy gaps between (_C(4)C(7) bonding and δ*_{C(1)  H(2)} antibonding orbitals decrease from compounds 1 to 3 . The results suggest also that in compounds 1 – 3 , the hydrogen eliminations are controlled by LP(e)→δ* resonance energies. Analysis of bond order, natural bond orbital charges, bond indexes, synchronicity parameters, and IRC calculations indicate that these reactions are occurring through a concerted and synchronous six‐membered cyclic transition state type of mechanism.     

Orbital Energy Order and the Through Bond Interaction in Homonuclear Diatomic Molecules

The homonuclear diatomic molecules are the simplest systems having both the σ framework and the lone pair orbitals n_a and _b for investigating their through space and through bond interaction. The striking orbital energy order n_g～ n_a+ n_b > n_n ～ n_a - n_b has been accounted for by the through bond interaction. However, when the p-content in the lone pair orbitals n_a and n_b decreases, one may have the reverse orbital energy order: n_g < n_g. A reverse orbital energy order has been found in F₂ and Cl₂, whose n_a and n_b are almost pure s-type atomic orbitals. The reverse order also occurs in molecule N₂ when the internuclear distance is larger man 1.5 Å. It is also found that the detail through space and through bond interaction and the eventual orbital energy order for n_g and n_u can be accounted for by the Fock operator within the localized molecular orbital space.