HPS2异构体结构与相对稳定性

在MP2/6-311++G(d,p)和QCISD(t)/6-311++G(3df,2p)(单点)水平下计算得到了包括9个异构体和15个过渡态的HPS2势能面.其中,异构体trans-HSPS(E1)的能量最低, 其次是cis-HSPS(E2)和HPS(S)(C2_V,E3),能量分别比trans-HSPS(E1)高3.43和14.17 kJ/mol.计算结果表明,异构体E1,E2,E3和立体的三元环结构的异构体HP(S)S(C_S,E4), 及与E4共存的trans-HPSS(E5)和cis-HPSS(E6)具有一定的稳定性.在QCISD(t)/6-311++G(3df,2p)//MP2/6-311++G(d,p)并包含零点能水平下,PH和S₂反应生成的E6和E5分别越过65.75和71.73 kJ/mol的势垒就可以异构化为具有较高动力学稳定性的产物E4, 计算结果对实验认定的产物是cis-HPSS(E6)的结论进行了修正.     

HAsO~2异构体结构、相对稳定性与体系势能面

在MP2/6-311++G(d,p)和QCISD(T)/6-311++G(3df,2p)(单点)水平下计算得到了包括9个异构体和10个过滤态的HAsO~2体系势能面。在势能面上，异构体cis-HOAsO(E1)的能量是最低的，其次是trans-HOAsO(E2)和HAsO(O)(C~2~V,E3)，能量分别比cis-HOAsO高13.15和192.74kJ/mol。根据体系的势能面，异构体E1，E3及cis-HOOAs(E6),trans-HOOAs(E5)具有一定的动力学稳定性，在实验中应该可以观测到。AsH和O~2反应的第一步产物将会异构化为具有较高动力学稳定性的异构体E3；而OH和AsO反应可直接生成E1。计算结果与HPO~2,HPS~2,HNO~2,HNS~2等价电子相同的分子的势能面进行了比较。     

HSCCS自由基势能面的理论研究

在CCSD(T)/6-311G(d,f)//MP2/6-311G(d,f) ZPE水平下,计算得到含有8个异构体和11个过渡态的HSCCS自由基体系势能面,讨论了异构体的结构与稳定性及异构体之间的异构化过程.结果表明异构体m1的能量最低,除m1以外,异构体m2和m3的能量也比较低,在MP2水平上,过渡态TS1的能量比异构体m2仅高3.9kJ/mol,而在CCSD(T)水平上,TS1的能量比m2低14.6 kJ/mol,这说明异构体m2可以迅速转化为能量更低的m1.异构体m3的能量比异构体m1高49.99 kJ/mol,可以推断,在合适的实验条件下可以观测到异构体m1.     

HAsS2异构体结构与稳定性

在MP2/6-311++G(d,p)和QCISD(t)/6-311++G(3df,2p)(单点)水平下计算得到9个异构体和10个过渡态的HAsS₂体系势能面.异构体cis-HSAsS(E1)的能量最低,其次是trans-HSAsS(E2)、具有AsSS三元环的立体HAs(S)S(C_{s,E3)和HAs(S)S(C2v,E4)结构的异构体,能量分别比cis-HSAsS高1.46,60.78和93.63kJ/mol.根据体系的势能面,异构体E1,E2,E3和E4具有一定的动力学稳定性.AsH和S2第一步反应产物将会异构化为具有较高动力学稳定性的异构体E3,而SH和AsS第一步反应产物将会异构化为E1.计算结果与HNO2,HNS2,HPO2,HPS2和HAsO2等价电子相同的分子的势能面进行了比较.}     

HPO_2异构体结构和相对稳定性

在MP2/6-311++G（d,p）和QCISD（t）/6-311++G（3df,2p）（单点）水平下计算得到了包括8个异构体和12个过渡态的HPO2体系势能面.在势能面上,异构体cis-HOPO（EI）的能量是最低的,其次是trans-HOPO（E2）和HPO（O）（C2v,E3）,能量分别比cis-HOPO高10.99和48.36 kJ/mol.根据体系的势能面,只有异构体E1和E3具有较高的动力学稳定性,在实验中应该可以观测到.PH和O2直接反应生成的cis-HPOO（E5）和trans-HPOO（E6）经过几步势垒较低的异构化过程就可以异构化为具有较高动力学稳定性的产物E1;而OH和PO反应可直接生成E1.计算结果较好地解释了相关实验.     

NC2S+离子的结构和稳定性的理论研究

在密度泛函和从头算理论水平下计算了单重态的NC2S+离子的结构、能量、光谱以及稳定性. 在B3LYP/6-311G(d)水平下, 得到8个异构体, 它们由15个过渡态相连接. 在CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE水平下, 得到能量最低的异构体是直线型的具有1Σ电子态的NCCS+(1)(0.0 kJ/mol), 其次是直线型的异构体CNCS+(2)(54.8 kJ/mol). 两个低能量的异构体1和2及另外一个高能量的直线型异构体CCNS+(3)(323.8 kJ/mol)都具有相当大的动力学稳定性, 这三个异构体在具备一定条件的实验室和星际条件下是可以进行观测的. 分析了这3个异构体的成键性质.     

HPO2异构体结构和相对稳定性

在MP2/6-311++G(d,p)和QCISD(t)/6-311++G(3df,2p)(单点)水平下计算得到了包括8个异构体和12个过渡态的HPO₂体系势能面. 在势能面上, 异构体cis-HOPO(E1)的能量是最低的, 其次是trans-HOPO(E2)和HPO(O)(C_2V, E3), 能量分别比cis-HOPO高10.99和48.36 kJ/mol. 根据体系的势能面, 只有异构体E1和E3具有较高的动力学稳定性, 在实验中应该可以观测到. PH和O₂直接反应生成的cis-HPOO(E5)和trans- HPOO(E6)经过几步势垒较低的异构化过程就可以异构化为具有较高动力学稳定性的产物E1; 而OH和PO反应可直接生成E1. 计算结果较好地解释了相关实验.     

CH_3CN+CH_3→CH_2CN+CH_4氢消除反应的直接动力学研究

采用直接动力学方法,对乙腈与甲基的反应进行了理论研究.在BHandHLYP/6-311G(d,p)和MP2/6-311G(d,p)水平下获得,稳定点的几何结构、振动频率及最小能量路径(MEP),在G3(MP2)和MC-QCISD水平下对能量信息进一步确认.利用正则变分过渡态理论,结合小曲率隧道效应校正(CVT/SCT)方法计算了该反应在220K～2000K的速率常数,与实验值符合得很好.     

高级量子化学从头计算法研究N2和H2O分子间相互作用

在MP2/6-311++G(3d,3p)电子相关校正水平上,对N₂和H₂O分子间可能存在的氢键复合物进行全自由度能量梯度优化,发现了一个接近于直线的弱氢键总能量极小结构(1),进一步在高级电子相关校正的MP4SDTQ和CCSD(T)水平,用6-311++G(3d,3p)基组加上(3s3p2d1f)键函数,用MP4和CCSD(T)计算的结构1的结合能分别为-5.061kJ/mol和-4.715kJ/mol.     

HPOS体系异构体结构与稳定性

在MP2／6－311＋＋G(d,p)和QCISD（T）／6－311＋＋G（3df,2p)（单点）水 平下计算得到了HPOS体系势能面上18个异构体和25个过渡态及解离碎片等驻点，并 分析了这些异构体的结构及异构化过程，讨论了可能的解离方式。在得到的异构体 中，有8个异构体是动力学较稳定的，它们是dis-HOPS,trans-HOPS,trans-HSPO, cis-HSPO,HP(O)S(Cs),trans-HPSO,cis-HPSO和HP（O）S（C1）。这些异构体在实 验中应该可以观测到。理论研究表明，P与S原子较强的超价能力在降低异构体能量 ，提高异构体动力学稳定性方面起到了关键的作用。得到的计算结果与HPO2， HPS2，HNOS等价电子相同的体系进行了比较。     

Structures and stabilities of HPO_2 isomers

The potential energy surface of HPO2 system including eight isomers and twelve transition states is predicated at MP2/6-311++G(d, p) and QCISD(t)/6-311++G(3df,2p)(single-point) levels of theory. On the potential energy surface, cis-HOPO(E1) is found to be thermodynamically and kinetically most stable isomer followed by trans-HOPO(E2) and HPO(O)(C2v, E3) at 10.99 and 48.36 kJ/mol higher, respectively. Based on the potential energy surface, only E1 and E3 are thermodynamically stable isomers, and should be experimentally observable. The products cis-HPOO(E5) and frans-HPOO(E6) in the first-step reaction of HP with O2 can isomerize into isomer E1 that has higher stability. The reaction of OH with PO will directly lead to the formation of isomer E1. The computed results are well consistent with the previous experimental studies.     

Structures and stabilities of HPS_2 isomers

The potential energy surface of HPS2 system containing nine isomers and fifteen transition states is obtained at MP2/6-311++G(d, p) and QCISD(t)/6-311++G(3df, 2p)(single-point) levels. On the potential energy surface, the lowest-lying frans-HSPS(EI) is found to be thermodynami-cally the most stable isomer followed by cis-HSPS(E2) and HP(S)S(C2v, E3) at 3.43 and 14.17 kJ/mol higher, respectively. The computed results show that species E1, E2, E3, stereo HP(S)S(Cs, E4) with PSS three-membered ring, isomers trans-HPSS(E5) and cis-HPSS(E6) which coexist with E4 are kinetically stable isomers. The products E6 and E5 in the reaction of HP with S2 can be isomerized into higher kinetic stable isomer E4 with 65.75 and 71.73 kJ/mol reaction barrier height, respectively. The predicated results may correct the possible inaccurate conclusion in that the product was experimentally assigned as isomer cis-HPSS(E6).     

Structures and stabilities of HPO<Subscript>2</Subscript> isomers

The potential energy surface of HPO₂ system including eight isomers and twelve transition states is predicated at MP2/6-311++G(d, p) and QCISD(t)/6-311++G(3df,2p)(single-point) levels of theory. On the potential energy surface, cis-HOPO(E1) is found to be thermodynamically and kinetically most stable isomer followed by trans-HOPO(E2) and HPO(O)(C_2v, E3) at 10.99 and 48.36 kJ/mol higher, respectively. Based on the potential energy surface, only E1 and E3 are thermodynamically stable isomers, and should be experimentally observable. The products cis-HPOO(E5) and trans-HPOO(E6) in the first-step reaction of HP with O₂ can isomerize into isomer E1 that has higher stability. The reaction of OH with PO will directly lead to the formation of isomer E1. The computed results are well consistent with the previous experimental studies.     

Structures and stabilities of HPS2 isomers

The potential energy surface of HPS² system containing nine isomers and fifteen transition states is obtained at MP2/6-311++G(d, p) and QCISD(t)/6-311++G(3df, 2p)(single-point) levels. On the potential energy surface, the lowest-lying trans-HSPS(E1) is found to be thermodynamically the most stable isomer followed by cis-HSPS(E2) and HP(S)S(C_2v, E3) at 3.43 and 14.17 kJ/mol higher, respectively. The computed results show that species E1, E2, E3, stereo HP(S)S(C_s, E4) with PSS three-membered ring, isomers trans-HPSS(E5) and cis-HPSS(E6) which coexist with E4 are kinetically stable isomers. The products E6 and E5 in the reaction of HP with S² can be isomerized into higher kinetic stable isomer E4 with 65.75 and 71.73 kJ/mol reaction barrier height, respectively. The predicated results may correct the possible inaccurate conclusion in that the product was experimentally assigned as isomer cis-HPSS(E6).     

Ab initio studies of the conformations of ethynyl formate, cyano formate and related ethynyl and cyano compounds

The optimized geometries, relative free energies and related thermodynamic properties, harmonic frequencies, and dipole moments have been calculated at the HF and MP2 levels for ethynyl formate (1a), ethynyl acetate (1b), cyano formate, HCO₂CN (1c), cyano acetate (1d), S-ethynyl thioformate (2a), S-ethynyl thioacetate (2b), S-cyano thioformate (2c), S-cyano thioacetate (2d), N-ethynylformamide (3a), N-ethynylacetamide (3b), N-cyanoformamide (3c), and N-cyanoacetamide (3d) with the gaussian 98 program. For ethynyl formate, the calculation for 25 °C at the MP2/6-311++G(df,pd) level predicts that the Z isomer is more stable by 1.23 kcal/mol. For S-ethynyl thioformate, calculations at the MP2/6-311++G(2d,2p) level predict that the E isomer is favored by 0.71 kcal/mol at 25 °C. The E isomers of N-ethynylformamide and N-ethynylacetamide were found at all levels to be more stable than the Z isomers at 25 °C. For cyano formate and cyano acetate, calculations at the MP2/6-311++G(df,pd) level predict that the Z isomers are more stable at 25 °C by 1.50 and 2.72 kcal/mol, respectively. At this level and temperature, the Z isomers of 2c, 2d, 3c, and 3d are predicted to have free energies of 0.46, −0.07, 1.22, and 2.28 kcal/mol, respectively, relative to the E conformations. Z to E free-energy barriers at 25 °C of 8.63, 10.64, 17.63, 7.39, and 14.03 kcal/mol were calculated for 1a, 2a, 3a, 1c, and 3c at the HF/6-311G(d,p) level, and at the HF/6-311+G(d,p) level, the free-energy barrier for 2c was 7.08 kcal/mol.     

Conformational stability, vibrational assignmenents, barriers to internal rotations and ab initio calculations of 2-aminophenol (d 0 and d3)

The Raman (3700-100 cm(-1)) and infrared (4000-400 cm(-1)) spectra of solid 2-aminophenol (2AP) have been recorded. The internal rotation of both OH and NH2 moieties produce ten conformers with either Cs or C1 symmetry. However, the calculated energies as well as the imaginary vibrational frequencies reduce rotational isomerism to five isomers. The molecular geometry has been optimized without any constraints using RHF, MP2 and B3LYP levels of theory at 6-31G(d), 6-311+G(d) and 6-31++G(d,p) basis sets. All calculations predict 1 (cis; OH is directed towards NH2) to be the most stable conformation except RHF/6-31++G(d,p) basis set. The 1 (cis) isomer is found to be more stable than 8 (trans; OH is away from the NH2 moiety and the NH bonds are out-of-plane) by 1.7 kcal/mol (598 cm(-1)) as obtained from MP2/6-31G(d) calculations. Aided by experimental and theoretical vibrational spectra, cis and trans 2AP are coexist in solution but cis isomer is more likely present in the crystalline state. Aided by MP2 and B3LYP frequency calculations, molecular force fields, simulated vibrational spectra utilizing 6-31G(d) basis set as well as normal coordinate analysis, complete vibrational assignments for HOC6H4NH2 and DOC6H4ND2 have been proposed. Furthermore, we carried out potential surface scan, to determine the barriers to internal rotations of NH2 and OH groups. All results are reported herein and compared with similar molecules when appropriate.     

几种含芯电子相关能修正的G2和G2(ACI)方法

针对冻结芯电子近似,在MP2/6-311G(d,p)级别上对G2、MP2/6-311G(d)和MP2/6-311G(d,p)级别上对G2(QCI)方法进一步考虑了芯电子相关能修正,尝试建立了G2(fu2)、G2(QCI/ful)和G2(QCI/fu2)方法。G2-l test set 的反应能量计算结果表明,这些方法进一步减小了经验修正量;G2(QCI/ful)和G2(QCI/fu2)也比G2(QCI)的总体精度有所提高;但G2(fu2)在G2基础上,总体精度没有改善。G2(fu2)、G2(QCI/ful)和G2(QCI/fu2)计算G2-l test set反应能量的平均绝对偏差分别为5.11、4.74和4.81kJ mol-1,G2和G2(QCI)分别为5.09和4.97kJ mol-1.