黄金规则用于N-3＋N3体系电子转移的研究

基于B3P86/6 311＋G优化的N3与N－3分子几何,确定了N3＋N－3基态电子转移体系 的六种不同的耦合机理,及各种形式耦合络合物的几何性质、活化能、稳定化能、耦合矩阵 元和态密度,并利用黄金规则计算了电子转移速率,讨论了各耦合方式对电子转移速率的影响 .     

黄金规则用于N_3~- N_3体系电子转移的研究

基于 B3P86/6 311＋ G优化的 N3与 N－ 3分子几何 ,确定了 N3＋ N－ 3基态电子转移体系的六种不同的耦合机理 ,及各种形式耦合络合物的几何性质、活化能、稳定化能、耦合矩阵元和态密度 ,并利用黄金规则计算了电子转移速率 ,讨论了各耦合方式对电子转移速率的影响 .     

电子转移反应O~2+O~2^.^-→O~2^.^-+O~2的从头算研究

利用abinitio方法,在UHF,UMP2及不同基组3-21G,6-31G^*,6-311+G^*和UMP2(full)/6-311+G^*水平上,研究了O~2/O~2^.^-自交换电子转移反应。优化了电子转移前后反应物和产物的结构,研究了体系能量的变化,计算了自交换电子转移反应的内重组能。对UHF方法和UMP2方法的计算结果进行了比较,并与实验结果进行了对照。结果表明UHF方法由于没有考虑组态相互作用,计算结果存在较大偏差,UMP2(full)/6-311+G^*水平上的计算结果与实验值吻合较好。在UMP2(full)/6-311+G^*水平上计算了气相自交换电子转移反应速率常数。在优化了电子转移复合物结构的基础上考虑了溶剂效应的影响,计算了水溶液中的溶剂重组能。研究结果表明O~2/O~2^.^-体系电子转移反应的活化能主要来源于溶剂重组能的贡献。最后计算了该反应在水溶液中的反应速率常数。理论计算结果与实验值吻合得很好。     

电子转移反应O2+O2·-→O2·-+O2的从头算研究

利用abinitio方法,在UHF,UMP2及不同基组3-21G,6-31G^*,6-311+G^*和UMP2(full)/6-311+G^*水平上,研究了O~2/O~2^.^-自交换电子转移反应。优化了电子转移前后反应物和产物的结构,研究了体系能量的变化,计算了自交换电子转移反应的内重组能。对UHF方法和UMP2方法的计算结果进行了比较,并与实验结果进行了对照。结果表明UHF方法由于没有考虑组态相互作用,计算结果存在较大偏差,UMP2(full)/6-311+G^*水平上的计算结果与实验值吻合较好。在UMP2(full)/6-311+G^*水平上计算了气相自交换电子转移反应速率常数。在优化了电子转移复合物结构的基础上考虑了溶剂效应的影响,计算了水溶液中的溶剂重组能。研究结果表明O~2/O~2^.^-体系电子转移反应的活化能主要来源于溶剂重组能的贡献。最后计算了该反应在水溶液中的反应速率常数。理论计算结果与实验值吻合得很好。     

1, 2-二硒方酸气相酸性和芳香性的理论研究

在RHF/6-311G^*^*水平优化得到1,2-二硒方酸(3,4-二羟基-3-环丁烯-1,2-二硒酮)三种平面构象异构体的平衡几何构型。进一步用MP2(full)/6-311G^*//RHF/6-311G^*^*方法计算三种异构体的单点能量,发现ZZ型异构体是能量最低构象,且ZZ和ZE型能量非常接近。用优化的最稳定构象ZZ型异构体在RHF/6-311G^*^*//RHF/6-311G^*^*,RHF/6-311+G^*^*//RHF/6-311+G^*^*,MP2(full)/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其气相酸性[ΔGⅲ~(~2~9~8~K~)]和同键反应芳香性稳定化能(HASE)。用基团加和法(groupincrementapproach)在RHF/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其磁化率增量(Λ)。计算结果指出标题化合物的键长发生了平均化,同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定。     

铁原子与氮分子间的相互作用——单侧双配位构型

用密度泛函理论的B3LYP/6-311+G(d)方法对单侧双配位FeN2体系(简记为S-FeN2)不同自旋多重度的稳定态、范德华力作用态和过渡态的多个电子态的几何结构、电子结构、能量和振动频率进行了计算比较研究. 结果表明, S-FeN2体系三种自旋态间, Fe—N 距离R1和N—N 距离R2值均比较接近; 能量最低的是15B2态, 相近态有15B1、13B1和13B2, 彼此能差约25 kJ·mol-1. 三重态电子结构复杂, 单重态能量普遍偏高; 基组态Fe原子与N2间存在强的σ-π电子对排斥而无有效轨道重叠和电子转移, 其它组态4s13d7、4s13d64p1和3d74p1, Fe 和N2间发生σ(sd)-π和π-π*轨道重叠作用, 有少量电子转移, 体系呈现一定的离子性特征, 活化N2键长基本不超过120 pm. Fe 原子的电子单或双重被激发到由N2反键轨道为主要成分的分子轨道上时, 能使N2活化到单键程度甚至解离.     

ClO/ClO~-耦合体系结构性质及反应通道研究

利用从头算法和密度泛函理论对ClO/ClO~-体系进行了研究。首先应用密度泛 函理论的六种方法(B3LYP,BLYP,B3P86,BP86,BHLYP,LSDA)和从头算理论的CCSD方法 在6-3+G~*,6-311+G~*及6-311+G(3df)基组水平上对单体进行了优化,选出最适合 该体系的方法和基组B3LYP/6-311+G(3df);然后在B3LYP/6-311+G(3df)水平上计算 了沿各种反应通道的离解能,并且进行了校正。找出了存在的中间体及中间体异构 化的过渡态,并进行了IRC路径解析。对各稳定体进行了频率分析和成键分析。结 果表明,单体ClO和单体ClO~-结合为稳定的中间体后,其离解方式主要是向着 ClOO+Cl~-和ClOO~-+Cl两种方式进行,两种离解方式的离解能分别为-33.39和82. 88kJ/mol,并且前者是一个离解性电子转移过程,经历一个电子转移过渡态。     

H+CH3NO2H2+CH2NO2反应途径和变分速率常数计算研究

采用MP2(FULL)/6-311G**从头算方法, 优化了H+CH3NO2H2+ CH2NO2反应的过渡态结构, 得出该反应的正逆反应的活化位垒分别是82.73和57.14 kJ*mol-1. 沿IRC分析指出该反应是一个H-H键生成和C-H键断裂的协同反应, 而且在反应途径上存在一个引导反应进行的振动模式, 这一反应模式引导反应进行的区间在-0.7～0.2( amu)1/2*a0之间; 在1 000～1 400 K温度范围内, 运用变分过渡态理论(CVT), 计算了该反应的速率常数, 计算结果与实验相一致.     

溶液中M（H2O）6^2＋／3＋（M=V,Cr,Mn,Fe,Co）自交换电子转移体系内层重组能和活

通过建立电子转移过程的活化模型和重组模型,提出了用量子化学从头算方法研究电子转移过程内层重组能和活化能的新方法。在UMP2／6－311G水平上获得了5对过渡金属水合离子体系M（H2O）6^2 /3 (M=V,Cr,Mn,Fe,Co)自交换反应的内层重组能和活化能,获得了与Marcus电子转移理论相一致的结果。     

C02+OH2/C03+OH2电子转移反应性的理论研究

提出了研究Co2+OH2/Co3+OH2反应体系电子转移反应性的接触距离依赖关系的分析方案和ab initio计算的应用方法,并讨论验证了此方案及其相应模型的可行性,分析了有关动力学量的接触距离依赖关系.详细的结果表明,用精确PES法得出的活化能与用非谐振子势得出的活化能吻合较好,它们明显优于谐势模型.对分布函数随接触距离从1.20～0.35 nm改变而从10-2变到10-5.偶合矩阵元随接触距离的增加呈指数性降低.有效电子偶合要求接触距离＜0.75 nm.在0.50～0.75 nm范围内,相应的电子发射系数值在1.0～10-6之间.电子因子使得定域ET速率也指数性的随接触距离的增加而降低,而对分布函数对总电子转移速率的贡献与电子因子的贡献则相反.球平均ET速率随接触距离的变化呈抛物线变化,并在接触距离为0.5 nm时有最大值.此最大值与总观测ET速率非常接近.对于此偶合体系,气态时ET速率是106 L@mol-1@s-1.进一步来说,实验上难于确定此类水合体系尤其是未饱和中间组分的电子结构和PES,abinitio算法在讨论其ET反应性方面能起到一个有效的辅助作用.     

Quantum-chemical Study on the Structures and Properties of Uracil-BX3 （X = F, Cl） Complexes

1 INTRODUCTION The intermolecular interaction of bases in DNA or RNA is of immense interest and significance to che- mists and biologists alike. The interactions of these bases with metal cations, solvent molecules and other small molecules or ions would affect the struc- ture and biological properties or recognition process,which has been investigated widely[1～8]. Boron contained compounds are electron deficient com- pounds and have been extensively used as catalysts in chemical react…     

类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应. 在B3LYP/6- 311＋G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C＝SiLiBr与HF, H2O或NH3发生插入反应的机理相同. QCISD/6-311＋＋G(d,p)//B3LYP/6-311＋G(d,p)计算的三个反应的势垒分别为148.62, 164.42和165.07 kJ•mol－1, 反应热分别为－69.63, －43.02和－28.27 kJ•mol－1. 相同条件下发生插入反应时, 反应活性都是H—F＞H—OH＞H—NH2.     

离子簇合物N+2Hen(n=1,2,3)的ab initio研究

利用Ｇａｕｓｓｉａｎ９２量子化学从头算程序，选择ＵＭＰ２／６－３１１Ｇ方法，对Ｎ＾＋２Ｈｅ，Ｎ＾＋２Ｈｅ２和Ｎ＾＋２Ｈｅ３等离子簇合物进行了几何优化，并通过频率分析确认了体系的稳定构型。     

NH_2(~2B_1)与CH_4反应机理的理论研究

在密度泛函理论 B3LYP/6 -31 1 G*水平下 ,研究了 NH2 与 CH4的反应机理 .通过振动频率和内禀反应坐标 ( IRC)分析 ,对反应过渡态进行了确认 .在 QCISD( T) /6 -31 1 G*水平下进行了单点能计算 ,并进行了零点能校正 ,结果表明 ,反应 NH2 + CH4NH3 + CH3 是主要的反应通道 .     

NH_(2-3)~(0-1 )离解能等的高级ab initio计算与评价

FeiQi等人[1]最近对H2N－H、H2N －H的离解能以及NH3、NH2的电离能，用真空紫外光电离实验进行了重新测定，并同时得到了H2N－H 的离解能．他们还在QCISD（T）／6－311＋G（3df，ZP）／／MPZ（full）／6刁IG（d）水平（GZ理论的参考水平）上对这些数据及其它相关分子的某些性质     

类锗烯H2C=GeLiCl与RH(R=F,OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了不饱和类锗烯H2C=GeLiCl与RH(R=F, OH, NH2)的插入反应. 在B3LYP/6-311+G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C=GeLiCl与HF、H2O 或NH3发生插入反应的机理相同. QCISD/6-311++G(d,p)//B3LYP/6-311+G(d,p)计算的三个反应的势垒分别为173.53、194.48和209.05 kJ·mol-1, 反应热分别为60.18、72.93和75.34 kJ·mol-1. 相同条件下发生插入反应时, 反应活性顺序都是H—F>H—OH>H—NH2.     

Mechanistic insights into triterpene synthesis from quantum mechanical calculations. Detection of systematic errors in B3LYP cyclization energies

Most quantum mechanical studies of triterpene synthesis have been done on small models. We calculated mPW1PW91/6-311+G(2d,p)//B3LYP/6-31G* energies for many C30H51O+ intermediates to establish the first comprehensive energy profiles for the cationic cyclization of oxidosqualene to lanosterol, lupeol, and hopen-3beta-ol. Differences among these 3 profiles were attributed to ring strain, steric effects, and proton affinity. Modest activation energy barriers and the ample exothermicity of most annulations indicated that the cationic intermediates rarely need enzymatic stabilization. The course of reaction is guided by hyperconjugation of the carbocationic 2p orbital with parallel C-C and C-H bonds. Hyperconjugation for cations with a horizontal 2p orbital (in the plane of the ABCD ring system) leads to annulation and ring expansion. If the 2p orbital becomes vertical, hyperconjugation fosters 1,2-methyl and hydride shifts. Transition states leading to rings D and E were bridged cyclopropane/carbonium ions, which allow ring expansion/annulation to bypass formation of undesirable anti-Markovnikov cations. Similar bridged species are also involved in many cation rearrangements. Our calculations revealed systematic errors in DFT cyclization energies. A spectacular example was the B3LYP/6-311+G(2d,p)//B3LYP/6-31G* prediction of endothermicity for the strongly exothermic cyclization of squalene to hopene. DFT cyclization energies for the 6-311+G(2d,p) basis set ranged from reasonable accuracy (mPW1PW91, TPSSh with 25% HF exchange) to underestimation (B3LYP, HCTH, TPSS, O3LYP) or overestimation (MP2, MPW1K, PBE1PBE). Despite minor inaccuracies, B3LYP/6-31G* geometries usually gave credible mPW1PW91 single-point energies. Nevertheless, DFT energies should be used cautiously until broadly reliable methods are established.     

Theoretical investigation of adenine–BX3 (X=F,Cl) complex

Geometries and binding energies are predicted at B3LYP/6-311+G* level for the adenine–BX₃ (X=F,Cl) systems and four conformers with no imaginary frequencies have been obtained for both adenine–BF₃ and adenine–BCl₃, respectively, and single energy calculations using much larger basis sets (6-311+G(2df,p)) and aug-cc-pVDZ were carried out as well. The most stable conformer is BF₃ or BCl₃ connected to N₃ of adenine and with the stabilization energy of 22.55 or 20.59 kcal/mol at B3LYP/6-311+G* level (BSSE corrected). The analyses for the combining interaction between BX₃ and adenine with natural bond orbital method (NBO) and the atom-in-molecules theory (AIM) have been performed. The results indicate that all the conformers were formed with σ–p type interactions between adenine and BX₃, in which pyridine-type nitrogen or nitrogen atom of amino group offers its lone pair electron to the empty p orbital of boron atom and the concomitances of charge transference from adenine to BX₃ were occurred. Frequency analysis suggested that the stretching vibration of BX₃ underwent a red shift in complexes. Adenine–BF₃ complex was more stable than adenine–BCl₃ although the distance of B–N is shorter in the later.     

Density functional theory studies on the dissociation energies of metallic salts: relationship between lattice and dissociation energies

The formation and physicochemical properties of polymer electrolytes strongly depend on the lattice energy of metal salts. An indirect but efficient way to estimate the lattice energy through the relationship between the heterolytic bond dissociation and lattice energies is proposed in this work. The heterolytic bond dissociation energies for alkali metal compounds were calculated theoretically using the Density Functional Theory (DFT) of B3LYP level with 6‐311+G(d,p) and 6‐311+G(2df,p) basis sets. For transition metal compounds, the same method was employed except for using the effective core potential (ECP) of LANL2DZ and SDD on transition metals for 6‐311+G(d,p) and 6‐311+G(2df,p) calculations, respectively. The dissociation energies calculated by 6‐311+G(2df,p) basis set combined with SDD basis set were better correlated with the experimental values with average error of ca. ±1.0% than those by 6‐311+G* combined with the LANL2DZ basis set. The relationship between dissociation and lattice energies was found to be fairly linear (r>0.98). Thus, this method can be used to estimate the lattice energy of an unknown ionic compound with reasonably high accuracy. We also found that the dissociation energies of transition metal salts were relatively larger than those of alkaline metal salts for comparable ionic radii. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 827–834, 2001     

Theoretical prediction of the contact distance dependence of the electron transfer reactivity of the ClO/ClO− coupling system

On the basis of the structures and properties of the ClO/ClO^? system obtained at the density functional theory (DFT) (UB3LYP) level, employing the 6‐311+G(3df) standard basis set, the electron transfer reactivity of this system is investigated. The results indicate that there are five possible stable coupling complexes that correspond to the generous minima on the global potential energy surfaces (PES). The most stable coupling complex is planar EC4, in which there is a O? O linkage with two trans‐Cl atoms. Their stabilization energies are calculated to be 20.57 (EC1: C₁), 20.54 (EC2: C₂, ²B), 20.69 (EC3: C₁), 20.70 (EC4: C_s, ²A′), and 20.69 (EC.5: C_2h, ²B_u) kcal/mol at the B3LYP/6‐311+G(3df) level; with the correction of the basis set superposition error (BSSE), the stability order of these encounter complexes is EC4 > EC.5 > EC3 > EC1 > EC2. Based on the five encounter complexes, five coupling modes are designed for the study of the electron transfer reactivity of this system. The dissociation energy curves at the activated states and the corresponding activation energies of these five coupling modes are obtained and are compared at the B3LYP/6‐311+G(3df) and MP2/6‐311+G* levels. The inapplicability of DFT methods has also been discussed in this article in predicting the energy curves, especially with a long contact distance, in which DFT methods give the abnormal behavior for the dissociations of the complexes caused by the “inverse symmetry breaking” problem. On the basis of the golden rule of the time‐dependent perturbation theory, the electron transfer reactivity and the contact distance dependence of the various electron transfer kinetics parameters (e.g., activation energy, coupling matrix element) have been analyzed at the UMP2(full)/6‐311+G* level. The electron transfer can take place over a range of contact distances, but the most effective coupling distance corresponds to only a small range. The coupling orientation analyses also indicate that the most favorable coupling mode to the electron transfer does not always correspond to the most stable encounter complex mechanism. Some highly energetic coupling modes are more favorable for the electron transfer. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005