不饱和类烯H2PBLiF的构型及异构化反应

用ＨＦ／６－３１Ｇ＊＊解析梯度方法研究了无机不饱和类烯Ｈ２ＰＢＬｉＦ的结构,共得到３个平衡构型和２个异构化反应的过渡态构型,动力学分析表明,其中两种平衡构型是它们存在和参加化学反应的基本构型,分析了各平衡构型的结构特点及稳定性,给出了各构型的Ｍｕｅｌｌｉｋｅｎ集居数,并简单讨论了２种基本构型的化学活性。     

锂氟类硅烯与乙烯加成反应的理论研究

用从头计算方法研究了锂氟类硅烯与乙烯的加成反应Ｈ_２ＳｉＬｉＦ＋Ｃ_２Ｈ_４→Ｈ_２ＳｉＣ_２Ｈ４＋ＬｉＦ．该反应的过渡态和类卡宾与乙烯的反应相似。反应前后的能量差经零点能校正后仅为－２．４ｋＪ／ｍｏｌ（ＭＰ２／６－３１Ｇ￣＊／／６－３１Ｇ￣＊）．本文分析了孩反应的热力学和动力学性质，计算了反应热力学函数的变化、平衡常数、Ａ因子以及速率常数。     

卤代类卡宾结构和反应活性的理论研究（Ⅱ）：类卡宾CHClFLi的 …

用从头算分子轨道理论，研究了不同卤代数卡宾ＣＨＣｌＦＬｉ的结构和稳定性，优化得到４种平衡构型和相应过渡态地各构型的特别进行了分析，并与相同卤代类卡宾做了比较。     

超价化合物NLi4^n＋和OLi4^n＋电子结构和性质

以６－３１Ｇ＾＊基组利用ＨＦ、ＭＰ２和ＤＦＴ方法优化了超价化合物ＮＬｉ４＾ｎ＋和ＯＬｉ４＾ｎ－的几何构型。研究结果表明，ＭＰ２和ＤＦＴ法计算出的ＯＬｉ４分子解离出Ｌｉ和Ｌｉ２的反应能与已有的实验值吻合，对于ＮＬｉ４分子，得到其解离出Ｌｉ和Ｌｉ２的反应能分别为１９１．７８和５１５．３７ｋＪ／ｍｏｌ（ＭＰ２值）。并预测了ＯＬｉ４＾ｎ＋和ＮＬｉ４＾ｎ＋分子的基振动频率。     

烷基对取代锂氟类硅烯R2SiLiF的构型和热稳定性的影响

利用从头计算法研究了二甲基锂氟类硅烯和二乙基锂氟类硅烯的各异构体的构型和能量.烷基的引入增加了Si的σ电子对Li的授与而降低了F对Si空p轨道的电子回授,使类硅烯的p-络合物构型热稳定性下降,而σ-络合物热稳定性增加.在上述两种作用的共同影响下,三元环构型的热稳定性只有很小的变化.     

过渡金属和硅烯配体相互作用的理论研究（I）：MSiH^＋2的 …

以ＨＦ／６－３１１＋Ｇ＾＊基组研究了过渡金属硅烯离子ＭＳｉＨ＾＋２（Ｍ从Ｓｃ至Ｃｕ）的构型、成键特征以及Ｍ－Ｓｉ键解离能。具有共平面构型的ＭＳｉＨ＾＋２,其Ｍ＝Ｓｉ键带有明显的双键特征,Ｍ＝Ｓｉ键解离能从Ｓｃ至Ｃｕ呈现周期性变化,与Ｍ的金属离子激发能有近似的线性关系。     

环丙基硅烯的理论研究环丙基硅烯C3H5SiH的重排反应及其机理

本文用限制的Ｈａｒｔｒｅｅ－Ｆｏｃｋ解析梯度方法在３－２１Ｇ和６－３１Ｇ＾＊水平上对环丙基硅烯的重排反应及其机理进行了从头算研究。以６－３１Ｇ＾＊优化构型作了二级微扰计算，并计算了各构型的频率。在此基础上得到了重排反应的热焓ΔＨ，自由能ΔＧ和平衡常数Ｋ，用Ｅｙｒｉｎｇ过渡态理论计算了反应的速度常数ｋ（Ｔ），应用Ｗｏｏｄｗａｒｄ－Ｈｏｆｆｍａｎｎ规则讨论了环丙基硅烯重排反应过程中端基的旋转机理。结果     

反式共轭类碳烯CH2=CH-CH=C:LiF的构型及异构化反应

用MP2／6－31G^**解析梯度方法研究了反式共轭类碳烯CH2＝CH－CH＝C：LiF的结构。得到3个平衡构型和2个异构化反应的过渡态。计算结果表明，反式共轭类碳烯有两种比较稳定的平衡构型。分析了各平衡构型的结构特点及稳定性，给出了各种平衡构型及过渡态结构的Mulliken集居数，并讨论了两种稳定构型的化学活性。     

无机不饱和类烯H~2NBLiF的构型及异构化反应的理论研究

用HF/D95^*^*解析梯度方法研究了无机不饱和类烯H~2NBLiF的结构,共得到4个平衡构型和3个异构化反应的过渡态构型。动力学分析表明,其中两种平衡构型是它们存在和参加化学反应的基本构型。分析了各平衡构型的结构特点及稳定性,给出了各构型的Mulliken集居数,并简单讨论了两种基本构型的化学活性。     

溶液中类锗烯H2GeClMgCl的结构与异构化反应

用密度泛函理论(DFF)和 QCISD (quadratic configuration interaction with single and double excitations)方法研究了类锗烯 H2CeCMgCl在气相和五种溶剂中的构型与异构化反应.结果表明,类锗烯 H2GeClMgCl有三种平衡构型.其中p-配合物型构型能量最低,是其存在的主要构型.讨论了溶剂效应对结构、能量与异构化反应的影响.计算模拟了最稳定构型的红外光谱.     

Solvent effects on one-bond B-Li coupling constants in boryllithium compounds

EOM-CCSD 11B-7Li coupling constants and B chemical shifts have been computed for Li-diazaborole and its complexes with one H2O or FLi molecule. B-Li coupling constants for a model compound H(2)BLi and its complexes with up to 4 H2O or FLi molecules have also been obtained in an attempt to resolve discrepancies between the computed values of these properties for isolated Li-diazaborole and experimentally determined values for boryllithium in a THF solution. The presence of solvent molecules increases the ion-pair character of the B-Li bond, with the result that 1J(B-Li) decreases systematically as the basicity and the number of solvent molecules increases. In the presence of even a single solvent molecule, the boron chemical shift for Li-diazaborole increases, and approaches the experimental value. The computed results emphasize the role of the solvent in determining these NMR properties.     

(HAlNH)~n(n=1～6)簇的结构、红外光谱和热力学性 质的研究

用自洽场理论(HF)和密度泛函理论(DFT)的B3LYP方法,在6-31G水平上研究了HAlNH的低聚物(HAlNH)~n(n=1～6)簇的几何构型、电子结构、红外光谱和化学热力学性质,并比较了(HAlNH)~n和(ClAlNH)~n两种低聚物对应结构中化学键强弱,分析了引起(AlN)~n骨架结构发生变化的原因。结果表明,(HAlNH)~n簇的基态结构为C~s(n=1),D~2~h(n=2),D~3~h(n=3),T~d(n=4),C~s(n=5)和D~3~d(n=6)对称点群。HAlNH基态结构中,Al-N键是三重键。在D~2~h(n=2)和D~3~h(n=3)结构中,所有Al-N键均为二重键。在T~d(n=4)和D~3~d(n=6)中,Al-N键为正常单键,而在C~s(n=5)结构中含有三种Al-N键：单键、双键和混合键。振动频率计算表明,结构a～f均为基态稳定结构。热力学计算给出的稳定性顺序为：f>d>e>c>b>a。     

Isomeric solid enols on ring- and amide-carbonyls of substituted 2-carbanilido-1,3-indandiones

Both isomeric enols on ring carbonyl (5b) and on amide carbonyl (6b) derived from N-p-methoxyphenyl-2-carbamido-1,3-indandione (4b) were isolated, and their X-ray structures were determined. X-ray diffraction of the N-o,p-dimethoxy analogue indicated a disorder ascribed to the presence of a 6:4 mixture of 5c and 6c. Calculation (B3LYP/6-31+G*) gave good agreement with observed geometries. The calculated energies indicated that enols 6 are more stable by <1 kcal/mol than enols 5 and much more stable than amides 4.     

2-Deoxy-beta-D-erythro-pentofuranose: hydroxymethyl group conformation and substituent effects on molecular structure, ring geometry, and NMR spin-spin coupling constants from quantum chemical calculations.

The effect of hydroxymethyl conformation (gg, gt, and tg rotamers about the C4-C5 bond) on the conformational energies and structural parameters (bond lengths, bond angles, bond torsions) of the 10 envelope forms of the biologically relevant aldopentofuranose, 2-deoxy-beta-D-erythro-pentofuranose (2-deoxy-D-ribofuranose) 2, has been investigated by ab initio molecular orbital calculations at the HF/6-31G level of theory. C4-C5 bond rotation induces significant changes in the conformational energy profile of 2 (2gt and 2tg exhibit one global energy minimum, whereas 2gg exhibits two nearly equivalent energy minima), and structural changes, especially those in bond lengths, are consistent with predictions based on previously reported vicinal, 1,3- and 1,4-oxygen lone pair effects. HF/6-31G-optimized envelope geometries of 2gg were re-optimized using density functional theory (DFT, B3LYP/6-31G), and the resulting structures were used in DFT calculations of NMR spin-spin coupling constants involving 13C (i.e., J(CH) and J(CC) over one, two, and three bonds) in 2gg according to methods described previously. The computed J-couplings were compared to those reported previously in 2gt to assess the effect of C4-C5 bond rotation on scalar couplings within the furanose ring and hydroxymethyl side chain. The results confirm prior predictions of correlations between 2J(CH), 3J(CH), 2J(CC) and 3J(CC), and ring conformation, and verify the usefulness of a concerted application of these couplings (both their magnitudes and signs) in assigning preferred ring and C4-C5 bond conformations in aldopentofuranosyl rings. The new calculated J-couplings in 2gg have particular relevance to related J-couplings in DNA (and RNA indirectly), where the gg rotamer, rather than the gt rotamer, is observed in most native structures. The effects of two additional structural perturbations on 2 were also studied, namely, deoxygenation at C5 (yielding 2,5-dideoxy-beta-D-erythro-pentofuranose 4) and methyl glycosidation at O1 (yielding methyl 2-deoxy-beta-D-erythro-pentofuranoside 5) at the HF/6-31G level. The conformational energy profile of 4 resembles that found for 2gt, not 2gg, indicating that 4 is an inappropriate structural mimic of the furanose ring in DNA. Glycosidation failed to induce differential stabilization of ring conformations containing an axial C1-O1 bond (anomeric effect), contrary to experimental data. The latter discrepancy indicates that either the magnitude of this differential stabilization depends on ring configuration or that solvent effects, which are neglected in these calculations, play a role in promoting this stabilization.     

Probing furanose ring conformation by gas-phase computational methods: energy profile and structural parameters in methyl beta-D-arabinofuranoside as a function of ring conformation

The potential energy surface of methyl beta-D-arabinofuranoside (3) has been studied by ab initio molecular orbital (HF/6-31G) and density functional theory (B3LYP/6-31G) calculations via minimization of the 10 possible envelope conformers. The partial potential energy surface identified that the global minimum and lowest energy northern conformer was E(2). In the HF calculations, (2)E was the most stable southern conformer, while the density functional theory methods identified (4)E as the local minimum in this hemisphere. Additional calculations at higher levels of theory showed that the B3LYP-derived energies of many of the envelope conformers of 3 are dependent upon the basis set used. It has also been demonstrated that B3LYP/6-31+G//B3LYP/6-31G single point energies are essentially the same as those obtained from full geometry optimizations at the B3LYP/6-31+G level. The northern and southern minima of the B3LYP/6-31+G surface are, respectively, the E(2) and (2)E conformers. The B3LYP/6-31G geometries were used to study the relationship between ring conformation and various structural parameters including bond angles, dihedral angles, bond lengths, and interatomic distances.     

Comparison of Some Computational Methods for Geometry Optimisation of Phosphorus Acid Derivatives

Several economical methods for geometry optimization, that should be applicable to larger molecules, have been evaluated for 19 phosphorus acid derivatives. MP2/cc-pVDZ geometry optimizations are used as reference points and the geometries obtained from the other methods are evaluated with respect to deviations in bond lengths and angles, from the reference geometries. The geometry optimization methods are also compared to the much used B3LYP/6-31G(d) method. Single point energies obtained by subsequent EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) calculations on the respective equilibrium geometries are also reported relative to the energies obtained from the reference geometries. The geometries from HF/MIDI! optimizations were closer to those of the references than the geometries of the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! optimizations. The EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) single point energies obtained from the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! geometries gave a mean absolute deviation (MAD) from that of the reference geometries of 1.4-3.9 kcal mol m 1 . The HF/MIDI! geometries, however, gave EDF1/6-31+G(d) and B3LYP/6-31+G(d,p) energies with a MAD of only about 0.5 and 0.55 kcal mol m 1 respectively from the energies obtained with the reference geometries. Thus, use of HF/MIDI! for geometry optimization of phosphorus acids is a method that gives geometries of near-MP2 quality, resulting in a fair accuracy of energies in subsequent single point calculations, at a much lower computational cost other methods that give similar accuracies.     

A DFT study on the double bond migration of butene catalyzed by ionic pair of 1-ethyl-3-methyl-imidazolium fluoride

The double bond migration of butene catalyzed by 1-ethyl-3-methyl-imidazolium fluoride (EmimF) has been studied using quantum chemical method. The geometries of reactant, transition state and product for the isomerization have been optimized by density functional theory (DFT) at the B3PW91/6-31G(d,p), 6-311++G(d,p) and aug-cc-PVDZ levels. The computed results show that the 4-H atom on imidazole ring of EmimF has a good catalytic activity to the double bond migration of butene and the catalytic reaction of 1-butene to 2-butene is a synergetic and elementary process. The apparent activation energy of isomerization is about 197 kJ/mol.     

Methyl and trifluoromethyl vinyl sulfoxides. Steric and electronic structure

Potential functions of internal rotation about the $C_{sp^2 } $ -S bond in H₂C=CHSOCY₃ (Y = H or F) were determined and stationary points were identified by vibrational analysis at the MP2(full)/6-31+G(d), B3PW91/6-31+G(d), and B3PW91/6-311+G(3df,p) levels. Energetically favorable conformations were established, and rotation barriers and molecular geometry parameters were evaluated. Wave functions [MP2(full)/6-31+G(d)] were analyzed by the NBO method. Energies and dipole moments of bond and lone-pair orbitals, principal types of donor-acceptor interactions, bond orders, and atomic charges were determined.     

A theoretical study of C4H5N isomers: Comparison of cyano and isocyano as substituents

Ab initio molecular orbital calculations using a 3-21G basis set have been used to optimize geometries for pyrrole, CH₃(X)CCH₂, CH₃(H)CCHX (both cis and trans), c-C₃H₅X, and CH₂CHCH₂X, where X is CN and NC. In all the alkenyl derivatives methyl groups are found to adopt the conformation in which the methyl hydrogen eclipses the double bond. 6-31G*∥3-21G level calculations show the alkenyl cyanides to be of similar energy to pyrrole, but the isocyanides are ～20 kcal mol^?1 higher in energy. For both substituents the cyclopropyl derivatives are higher in energy by ～10 kcal mol^?1. At the 6-31G* level ring strain is 27.7 kcal mol^?1 for the cyanide and 30.6 kcal mol^?1 for the isocyanide. Data on the relative energies of RCN and RNC are compared when R is (i) a saturated hydrocarbon, (ii) an unsaturated hydrocarbon, (iii) an α-carbenium ion, (iv) an allyl cation, and (v) an α-carbanion.     

Ab initio calculation of inner-sphere reorganization energies of arenediazonium ion couples

The geometries of a series of substituted arenediazonium cations (p-NO2, p-CN, p-Cl, p-F, p-H, m-CH3, p-CH3, p-OH, p-OCH3, p-NH2) and the corresponding diazenyl radicals were optimized at the HF/6-31G, MP2/6-31G, B3LYP/6-31G, B3LYP/TZP, B3PW91/TZP, and CASSCF/6-31G levels of theory. Inner-sphere reorganization energies for the single electron-transfer reaction between the species were computed from the optimized geometries according to the NCG method and compared to experimental values determined by Doyle et al. All levels of theory predicted a CNN bond angle of 180 degrees in the cation. A bent neutral diazenyl radical was predicted at all levels of theory excepting B3LYP/TZP and B3PW91/TZP for the p-Cl-substituted compound. Inner-sphere reorganization energies determined at the HF, MP2, and CASSCF levels of theory correlated poorly with both experimental results and calculated geometries. Density functional methods correlated best with the experimental values, with B3LYP/6-31G yielding the most promising results, although the ROHF/6-31G survey also showed some promise. B3LYP/6-31G calculations correctly predicted the order of the inner-sphere reorganization energies for the series, excluding the halogen-substituted compounds, with values ranging from 42.8 kcal x mol(-1) for the p-NO2-substituted species to 55.1 kcal x mol(-1) for NH2. The magnitudes of these energies were lower than the experimental by a factor of 2. For the specific cases examined, the closed-shell cation geometries showed the expected geometry about the CNN bond, with variations in the CN and NN bond lengths correlating with the electron-donating/withdrawing capacity of the substituent. As predicted by Doyle et al., a large geometry change was observed upon reduction. The neutral diazenyl radicals showed a nominal CNN bond angle of 120 degrees and variations in the CN and NN bond lengths also correlated with the electron-donating/withdrawing capacity of the substituent. Changes in theta(CNN) and r(CN) both correlated well with calculated lambda(inner). The key parameters influencing inner-sphere reorganization energy were the CN and NN bond lengths and the CNN bond angle. This influence is explained qualitatively via resonance models produced from NRT analysis and is related to the amount of CN double bond character. Based on these observations, B3LYP/6-31G calculations are clearly the most amenable for calculating inner-sphere reorganization energies for the single electron-transfer reaction between cation/neutral arenediazonium ion couples.