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.     

C24和B12N12团簇的结构与稳定性

本文采用量子化学密度泛函理论的B3LYP/6-31G*方法,对C24和B12N12团簇的12种异构体进行了优化,并对它们的几何构型、振动频率、核独立化学位移(NICS)和结合能进行了理论探讨, 比较了C24和B12N12团簇结构的稳定性。研究表明:C24团簇的最稳定几何构型为类石墨结构d,B12N12团簇的最稳定结构为4/6笼状结构g。C24异构体的稳定性大小顺序为d > b > f > c > a > e。B12N12团簇异构体稳定性大小顺序为a > f> c> d > e >b。     

3,4-二硫方酸芳香性和气相酸性的理论研究

在RHF/6-311G^**,RHF/6-311+G^**和B3LYP/6-311+G^**水平优化得到3,4-二硫方酸(3,4-二巯基-3-环丁烯-1,2-二酮)三种平面构象异构体的平衡几何构型.用MP2(Full)/6-311G^**//RHF/6-311G^**方法计算单点能量,发现ZZ型异构体是能量最低构象,且ZZ和ZE型能量非常接近.用优化的最稳定构象ZZ型异构体在RHF/6-311G^**//RHF/6-311G^**,RHF/6-311G^**//RHF/6-311G^**,MP2(Full)/6-311G^**//RHF/6-311G^**和B3LYP/6-311G^**//B3LYP/6-311G^**水平计算其气相酸性(ΔG⁰)和同键反应芳香性稳定化能(HASE).用基团加和法(Group Increment Approach)在RHF/6-311G^**//RHF/6-311G^**和B3LYP/6-311G^**//B3LYP/6-311G^**水平计算其磁化率增量(Λ).计算结果表明,标题化合物的同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定.     

直链小碳簇Cn(n=3～6)微正则解离速率常数的RRKM理论计算(Ⅰ)

用从头计算法研究小碳簇的解离通道及其动力学.以MP2/6-31G^*精度优化直链C₃,C₄,C₅,C₆及其过渡态的结构,并对它们进行振动分析.在此基础上计算了各解离通道的能垒,并根据RRKM理论估算各个通道的微正则解离速率.计算结果表明,这些小碳簇的解离主要表现为均裂方式.     

B2Be2簇的结构与成键性质的研究

用HF／6-31Gabinitio法对B2Be2簇9种异构体27个电子态的结构进行了全构型优化,再用大基组二次组态相互作用QCISD（T）／6-311G方法进行单点计算，9种异构体的稳定性顺序是：h＞i＞g＞e＞f＞c＞a＞d＞b.通过结构转变中的Walsh图、能量间隙、键数参数（BNP）和键电荷的研究，揭示了B2Be2的成键特征.     

直链小碳簇Cn(n=7～10)微正则解离速率常数的RRKM理论计算(Ⅱ)

用从头算法研究直链小碳簇C₇,C₈,C₉和C₁₀的解离通道及其动力学.在MP2/6-31G^*精度上优化了这些碳簇及其过渡态的结构,并对它们进行了振动分析.计算了各解离通道的活化能,并根据RRKM理论估算了各个通道的微正则解离速率,计算结果说明它们的主要解离通道为裂解C₃碎片的方式,这与实验所观察到的小碳簇的解离方式完全一致.     

乙炔基自由基C2H与氧气反应的密度泛函理论研究

应用量子化学从头算和密度泛函理论(DFT)对C₂H自由基和O₂的反应进行了研究.在B3LYP/6-311G^**水平上优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型,并计算出它们的振动频率和零点振动能(ZPVE).各物种的总能量由CCSD(T)/6-311G^**//B3LYP/6-311G^**给出,并对能量进行了零点能校正.计算结果表明,反应物中自由基C₂H中的边端C进攻O₂形成了中间体1 (HCCOO),中间体1是一个加合产物.由中间体1经过不同的反应通道可以生成不同的产物P₁ (HCO+CO), P₂ (HCCO+O), P₃(CO₂+CH), P₄ (C₂O+OH)和P₅ (2CO+H).反应通道之间存在着竞争机制.其中P₁, P₂是主要产物,其次还有一定比例的P5生成,而产物P₃, P₄的生成几率较低.各条反应通道化学反应热的计算与实验吻合较好.     

AlmN和AlmN+(m=2～9)团簇结构与稳定性的量子化学研究

用密度泛函理论(DFT)的B3LYP方法,在6-311G^*水平上对Al_mN和Al_mN⁺(m=2～9)团簇的几何构型、电子结构、振动频率等性质进行了理论研究,给出了以Al_m团簇作为设计Al_mN类结构的母体,考虑在不同位置上结合N原子的结构,可以较快找到Al_mN类团簇基态结构的一种方法.通过对基态结构的离解能和能量二次差分讨论,得到m为奇数的Al_mN团簇比m为偶数的稳定.对基态结构的绝热电离能讨论结果表明,只存在Al-N键的Al2N和Al3N团簇较稳定.     

NCS自由基与NO反应动力学的理论研究

用量子化学密度泛函理论B3LYP/6-31+G^*和高级电子相关校正的偶合簇[CCSD(T)/6-311+G^*]方法,对NCS自由基与NO反应的机理和动力学进行了理论研究,得到了体系的势能面信息和可能的反应机理.计算了反应的热力学参数及反应能垒.采用传统过渡态理论计算了各反应通道的速率常数.研究结果表明,NCS自由基与NO反应中存在4个反应通道,产物分别为OCS+N₂,CS+N₂O,ONS+CN和ONCNS.从能量变化和反应速率两方面考虑,NCS+NOOCS+N₂应为主反应通道.     

Si2Br6的分子振动光谱的理论研究

用量化从头算方法(HF/6-31G^*)和密度泛函方法(B3LYP/6-31G^*)以6-31G标准基组加一个极化函数,对Si₂Br₆分子的平衡几何构型和振动频率分别进行优化和计算,优化的结果与实验结果吻合得较好.按照Pulay的建议对HF/6-31G^*水平上所计算的谐性力场进行标度(标度因子取0.9).用HF/6-31G^*SQM力场所计算的基频预测值和实验值的平均误差为9.4cm^-1,最大误差为23.6cm^-1;用B₃LYP/6-31G^*未标度力场所计算的基频预测值和实验值的平均误差为8.6cm^-1,最大误差为16.6cm^-1;用该密度泛函方法所计算的基频预测值比用HF/6-31G^*的标度后的SQM力场所计算的基频预测值和实验值(除Si-Si键扭转振动基频之外的11条振动基频)吻合得更好.HF/6-31G^*和B3LYP/6-31G^*计算给出Si-Si键扭转振动基频的预测值分别为14cm^-1和9cm^-1.     

A theoretical study on three conformational structures of 2,6‐distyrylpyridine

Ab initio methods at the levels HF/cc‐pVDZ, HF/6‐31G(d,p), MP2/cc‐pVDZ, and MP2/6‐31G(d,p), as well as methods based on density functional theory (DFT) employing the hybrid functional B3LYP with the basis sets cc‐pVDZ and 6‐31G(d,p), have been applied to study the conformers of 2,6‐distyrylpyridine. Bond distances, bond angles, and dihedral angles have been calculated at the B3LYP level. The calculated values were in good agreement with those measured by X‐ray diffraction analysis of 2,6‐distyrylpyridine. The values calculated using the Hartree‐Fock method and second‐order perturbation theory (MP2) were inconsistent. The optimized lowest‐energy geometries were calculated from the reported X‐ray structural data by the B3LYP/cc‐pVDZ method. Three conformations, A, B, and C, were proposed for 2,6‐distyrylpyridine. Calculations at the three levels of theory indicated that conformation A was the most stable structure, with conformations C and B being higher in energy by 1.10 and 2.57 kcal/mol, respectively, using the same method and basis function. The same trend in the relative energies of the three possible conformations was observed at the two levels of theory and with the different basis sets employed. The reported X‐ray data were utilized to optimize total molecular energy of conformation A at the different calculation levels. The bond lengths, bond angles, and dihedral angles were then obtained from the optimized geometries by ab initio methods and by applying DFT using the two basis functions cc‐pVDZ and 6‐31G(d,p). The values were analyzed and compared. The calculated total energies, the relative energies of the molecular orbitals, the gap between them, and the dipole moment for each conformational structure proposed for 2,6‐distyrylpyridine are also reported. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

Study of conformational and optical rotation for the alaninamide

Six stationary points of alaninamide have been located on the potential surface energy (PES) at the B3LYP/6‐311++G(2d,2p) level of theory both in the gas phase and in aqueous solution. In the aqueous solution, to take the water solvent effect into account, the polarizable continuum model (PCM) method has been used. Accurate geometric structures and their relative stabilities have been investigated. The results show that the intramolecular hydrogen bond plays a very important role in stabilizing the global minimum of the alaninamide. Moreover, the consistent result in relative energy using high‐level computations, including the MP2 and MP3 methods with the same basis set [6‐311++G(2d,2p)], indicates that the B3LYP/6‐311++G(d,p) level may be applied to the analogue system. More importantly, the optical rotation of the optimized conformers (both in the gas phase and in aqueous solution) of alaninamide have been calculated using the density functional theory (DFT) and Hartree–Fock (HF) method at various basis sets (6‐31+G*, 6‐311++G(d,p), 6‐311++G(2d,2p) and aug‐cc‐pvdz). The results show that the selection of the computation method and the basis set in calculation has great influence on the results of the optical rotations. The reliability of the HF method is less than that of DFT, and selecting the basis set of 6‐311++G(2d,2p) and aug‐cc‐pvDZ produces relative reliable results. Analysis of the computational results of the structure parameters and the optical rotations yields the conclusion that just the helixes in molecules caused the chiral molecules to be optical active. The Boltzmann equilibrium distributions for the six conformers (both in the gas phase and in the aqueous solution) are also carried out. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A steric deuterium isotope effect in 1,1,3,3-tetramethylcyclohexane

The equilibrium isotope effect (EIE) for the interconversion of the two chair isotopomers of 1-trideutero-1,3,3-trimethylcyclohexane was predicted using geometry and vibrational force constants derived from electronic structure theory at HF, B3LYP, and MP2 levels as input for the program THERMISTP. Agreement between theory and previously reported NMR results is very good (experimental K(eq) = 1.042 +/- 0.001 vs K(eq) = 1.0409 at MP2/6-311G* level, K(eq) = 1.0503 at HF/6-311G*, and K(eq) = 1.0417 at B3LYP/6-311G* level, all at 17 degrees C). In order to investigate the origin of this isotope effect, the calculated EIEs for the monodeuterated isotopomers were analyzed. It has been shown that the hydrogen atom on an axial methyl group which is unusually close to its counterpart on the other axial methyl is responsible for the large (steric) isotope effect in the compound studied.     

Ab initio and DFT energetics of silylenic X–CNSi (X=H, F, Cl, and Br)

Singlet–triplet energy splitting for 24 silylenic reactive intermediates, X–CNSi (where X=H, F, Cl and Br), are compared and contrasted at 11 levels of theory: B1LYP/6-31++G**, B3LYP/6-31++G**, B1LYP/6-311++G**, B3LYP/6-311++G**, MP3/6-31G*, MP3/6-311++G**, MP2/6-31+G**, MP2/6-311++G**, MP4 (SDTQ)/6-311++G**, QCISD(T)/6-311++G** and CCSD(T)/6-311++G**. Each X-substituted silylenic species may either be singlet (s) or triplet (t), with one of the following three structures: 3-X-2-aza-1-silacyclopropenylidene (1_s-X, 1_t-X); [(X-imino)methylene]silylene (2_s-X, 2_t-X); and X-cyanosilylene (3_s-X, 3_t-X). For all X–CNSi species studied, orders of singlet–triplet energy separations (ΔE_s-t,X), appear as a function of electro-negativity (F>Cl>Br>H). For the six H–CNSi isomers (X=H), stability order is: 3_s-H>1_s-H>2_t-H>3_t-H>2_s-H>1_t-H. Likewise, stability order for the six isomers with X=F, is: 3_s-F>3_t-F>1_s-F>1_t-F>2_s-F>2_t-F. For X=Cl, the order of stability is: 3_s-Cl>1_s-Cl>3_t-Cl>2_t-Cl>1_t-Cl>2_t-Cl. Finally, the order of stability for six isomers of Br–CNSi is: 3_s-Br>3_t-Br>1_s-Br>2_s-Br>2_t-Br>1_t-Br. The lowest energy minimum, among all 24 species scrutinized, appears to be the singlet acyclic 3_s-X. Triplet silylene 2_t-H is suggested to be more stable than its corresponding 2_s-H at MP3, MP2 and DFT levels of theory. Comparisons between relative stabilities; multiplicities and geometrical parameters of 1–3 are discussed.     

Investigation on Structure and Bonding of As_4S_4 Isomers

Mineral medicine, especially those containing heavy metals, is one of the characteristics of traditional Chinese medicine. A famous mineral medicine, realgar, containing heavy metal arsenic with a chemical formula of As_4S_4, has the function of detoxification, killing bacteria and viruses, and eliminating dampness and phlegm. Different As_4S_4 isomers are likely to have different drug effects and pharmacological actions. Therefore, it is of great scientific significance to find more stable As_4S_4 isomers. In view of this, ab initio molecular orbital theory and density functional theory(DFT) have been used to study ten isomers of As_4S_4 at the B3LYP/6-31 G*, B3LYP/6-311+G*, B3LYP/6-311+G(3 df, 2 p) and MP2/(6-311+G*, LanL2 MB) levels of theory. In addition to the two isomers having been studied previously, eight new isomers were investigated in the present paper. All the ten As_4S_4 isomers were proved to be true local minima on their potential energy surfaces. The calculated NICS values and molecular orbital analyses showed that, the D_(2d) symmetric As_4S_4, isomer 1, may be s-aromatic. The study proves that ten As_4S_4 isomers are stable thermodynamically, and are highly desirable for the future theoretical study of realgar.     

2,3,7,8-四氯苯并二英分子结构与热力学性质的理论研究

用密度泛函理论(DFT)和从头算(ab initio)方法,在B3LYP/6-31G、 B3LYP/6-31G*、 B3LYP/6-311G*和MP2/6-31G*水平上全优化计算了2,3,7,8-四氯苯并二英(2,3,7,8-TCDD)的几何构型、电子结构和振动频率,并用校正后的频率计算了298～1500 K的标准热力学函数,同时用半经验的PM3 SCF-MO进行了同样的计算,计算结果与实验值及文献值较好地吻合.     

Comparative study on formamide–water complex

The hydrogen bonding of 1:1 complexes formed between formamide and water molecule have been investigated systematically using Hartree–Fock (HF), hybrid density functional theory (B3LYP), and post‐Hartree–Fock (MP2 and CCSD(T)) methods with range of basis sets 6‐31G(d), cc‐pVXZ (X = D, T, Q) and aug‐cc‐pVYZ (Y = D, T). Three stable structures are considered on the potential energy surface of formamide and water system. The optimized geometric parameters and interaction energies for various isomers at different levels are estimated. The IR frequencies, intensities, and frequency shifts are reported. This study shows that B3LYP/aug‐cc‐pVDZ method gives better performance for formamide‐water complexes. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010.     

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

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