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.     

A combined density functional theory and coupled cluster method investigation of the structural properties and stabilities of radical CH(2)CP and its isomers

The doublet potential energy surface of radical system [C(2), H(2), P] is investigated at the UB3LYP/6-311++G(d,p) and UCCSD(T)/6-311++G(2df,2p) (single-point) levels. Eight chainlike and three-membered ring structures are located as energy minima connected by 10 interconversion transition states. At the final UCCSD(T)/6-311++G(2df,2p)//UB3LYP/6-311++G(d,p) level with zero-point vibrational energy correction, species CH(2)CP is found to be thermodynamically the most stable isomer followed by HCCPH, H-cCPC-H, cPCC-H(H), H-cCCP-H, cis-CC(H)PH, trans-CC(H)PH, and CCPH(2) at 11.01, 12.57, 40.07, 43.63, 50.25, 56.82, and 65.36 kcal/mol, respectively. The computed results indicate that the chainlike isomers CH(2)CP and HCCPH and cyclic radical H-cCPC-H possess considerable kinetic stability at extra low pressures and temperatures. Interestingly, radical CCPH(2), whose energy is the highest in all predicted CH(2)CP isomers, can be also regarded as a kinetically stable species with the smallest isomerization barrier of 22.26 kcal/mol at extra low pressures and temperatures. Therefore, considering higher kinetic stability, in addition to the microwave spectroscopy characterized isomer CH(2)CP in previous experiments, the species HCCPH, H-cCPC-H, and CCPH(2) should be considered as excellent candidates for possible experimental observation. Furthermore, the structural nature of stable radical isomers is discussed based on bonding characteristics, single electron spin distribution, and comparison with their analogues.     

Predicting Potential Stable Isomers on the Singlet Surface of the [H,P,C,S] System by the MP2 and QCISD(T) Methods

A detailed singlet potential energy surface of [H,P, C,S] system is investigated by means of the MP2 and QCISD(T) methods. Eight isomers are located on the potential energy surface, and at the final QCISD(T)/6-311++G (3df,2p)//MP2/6-311++G(d,p) level with zero-point energy correction, the chainlike isomer HPCS is found to be kinetically and thermodynamically the most stable species followed by the chainlike HSCP, planar three-membered ring HC(S)P, chainlike HCPS, and stereo three-membered ring HP(C)S, which are predicted to be also kinetically stable isomers and should be experimentally observable provided that accurate experimental conditions are available. The dissociation processes from the kinetically and thermodynamically most stable species HPCS to the low-lying molecular dissociation fragments are not more favorable in energy than the isomerization process from HPCS to HSCP. Therefore, the experimental observation for potential isomer HSCP with C ≡ P triple bond is possible by means of photoisomerization technology using HPCS as precursor.     

Global investigation on the potential energy surface of CH3CN: application of the scaled hypersphere search method

An extensive quantum chemical study of the potential energy surface (PES) for all possible isomerization and dissociation reactions of CH3CN is reported at the DFT (B3LYP/6-311++G(d,p)) and CCSD(T)/ cc-pVTZ//B3LYP/6-311++G(d,p) levels of theory. The pathways around the equilibrium structures can be discovered by the scaled hypersphere search (SHS) method, which enables us to make a global analysis of the potential energy surface for a given chemical composition in combination with a downhill-walk algorithm. Seventeen equilibrium structures and 59 interconversion transition states have been found on the singlet PES. The four lowest lying isomers with thermodynamic stability are also kinetically stable with the lowest conversion barriers of 49.69-101.53 kcal/mol at the CCSD(T)/cc-pVTZ//B3LYP/6-311++G(d,p) level, whereas three-membered-ring isomers c-CH2NCH, c-CH2CNH, and c-CHNHCH can be considered as metastable intermediates which can further convert into the low-lying chain-like isomers and higher lying acyclic isomers with the lowest conversion energies of 21.70-59.99 kcal/mol. Thirteen available dissociation channels depending on the different initial isomers have been identified. A prediction can be made for the possible mechanism explaining the migration of a hydrogen atom in competition with the CC bond dissociation. Several new energetically accessible pathways are found to be responsible for the migration of the hydrogen atom. The present results demonstrate that the SHS method is an efficient and powerful technique for global mapping of reaction pathways on PESs.     

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.     

Theoretical study on the mechanism of the gas-phase reaction of diborane(3) anion with carbon disulfide

The complex potential energy surface of the gas-phase reaction of HB(H)BH- with CS2 to give three low-lying products [B2H3S]- + CS, [BH2CS]- + HBS, and [BH3CS] + BS-, involving nine [B2H3CS2]- isomers and 12 transition states, has been investigated at the CCSD(T)/6-311++G(d,p)/B3LYP/6-311++G(d,p) level. Our calculations are in harmony with the recent experimental and theoretical results, and reveal some new bonding and kinetic features of this reaction system. Our theoretical results may help the further identification of the products [BH2CS]- + HBS and [BH3CS] + BS- and may provide useful information on the chemical behaviors of other electron-deficient boron hydride anions.     

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).     

Structures and Stability of HNS_2 Isomers

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｎｉｔｒｏｇｅｎ ｓｕｌｆｕｒｃｏｎｔａｉｎｉｎｇｃｏｍｐｏｕｎｄｓｈａｖｅａｔ ｔｒａｃｔｅｄｍｕｃｈａｔｔｅｎｔｉｏｎｉｎｔｈｅｆｉｅｌｄｓｏｆｐｏｌｙｍｅｒｉｃｍａｔｅｒｉａｌｓａｎｄａｔｍｏｓｐｈｅｒｅｃｈｅｍｉｓｔｒｙ .1 4 Ａｎｄｓｏｍｅａｖａｉｌａｂｌｅｅｘｐｅｒｉ ｍｅｎｔａｌａｎｄｔｈｅｏｒｅｔｉｃａｌｉｎｆｏｒｍａｔｉｏｎｗａｓｒｅｐｏｒｔｅｄａｂｏｕｔｔｈｅｓｅｃｏｍｐｏｕｎｄｓｗｉｔｈＨ ｅｌｅｍｅｎｔｏｒｏｔｈｅｒｆｕｎｃｔｉｏｎａｌｇｒｏｕｐｓ…     

Theoretical Studies on Structures and Stabilities of C4H2+ Isomers

The structures, energies, stabilities and spectroscopies of doublet C₄H₂⁺ cations were explored at the DFT/B3LYP/6-311G(d,p), CCSD(T)/6-311+G(2df,2pd)(single-point), and G3B3 levels. Ten minimum isomers including the chainlike, three-member-ring, and four-member-ring structures are interconverted by means of 15 interconversion transition states. The potential energy surface was investigated. At the CCSD(T)/6-311+G(2df,2pd) and G3B3 levels, the global minimum isomer was found to be a linear HCCCCH. The structures of the stable isomer and its relevant transition state are further optimized at the QCISD/6-311G(d,p) level. The bonding nature and structure of isomer HCCCCH were analyzed.     

Theoretical prediction regarding structural and thermodynamical characteristics of stable CH3PO2 isomers and unimolecular decomposition mechanisms of species CH3P(O)2, CH3O?PO,and CH2P(O)OH

The detailed isomerization and dissociation reaction potential energy profile of the CH₃PO₂ system was established at the UCCSD(T)/6‐311++G(3df,2p)//UB3LYP/6‐311++G(d,p) level of theory. Seventy minimum isomers were located and connected by 93 optimized interconversion transition states. Furthermore, 32 isomers with high kinetic stability were predicted to be possible candidates for further experimental detection. The bonding nature of the suggested stable isomers was analyzed while their molecular properties including heats of formation, adiabatic ionization potentials, and adiabatic electronic affinities were calculated at the G2, G2(MP2), G3, and CBS‐Q levels. Based on the isomerization and dissociation potential energy surface, possible unimolecular decomposition mechanisms and pathways of the low‐lying molecules CH₃P(?O)₂, CH₃O? P?O, and CH₂?P(?O)OH were discussed. Furthermore, the transition state theory rate constants of the primary unimolecular dissociation channels were also calculated. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Theoretical investigation on the protonation reactions and products of the stable [N,C,C,S] isomers

A theoretical study on the protonation system of [N,C,C,S], [H,N,C,C,S]+, was performed at the B3LYP/6-311++G(d,p) and CCSD(T)/6-311++G(2df,2p) (single point) levels of theory. On the doublet [H,N,C,C,S]+ surface, 24 species were located as energy minima and 10 of them were considered as kinetically stable species. The species HNCCS+ with 2A' state and a shallow W-shaped skeleton was predicted to be the global minimum and kinetically the most stable species, being in good agreement with previous experimental findings. Furthermore, the protonation reactions of the stable [N,C,C,S] isomers were investigated in detail. The calculation results indicated that the [N,C,C,S] isomers may be significantly stabilized upon protonation. Finally, the possible covalent structures of the [H,N,C,C,S]+ isomers with considerable stability were briefly discussed.     

Theoretical studies on dynamics and thermochemistry of the reactions CF(3)CHCl(2)+Cl-->CF(3)CCl(2)+HCl and CF(3)CHFCl+Cl-->CF(3)CFCl+HCl

A dual-level direct dynamics study has been carried out for the two hydrogen abstraction reactions CF(3)CHCl(2)+Cl and CF(3)CHFCl+Cl. The geometries and frequencies of the stationary points are optimized at the BHLYP/6-311G(d,p), B3LYP/6-311G(d,p), and MP2/6-31G(d) levels, respectively, with single-point calculations for energy at the BHLYP/6-311++G(3df,2p), G3(MP2), and QCISD(T)/6-311G(d,p) levels. The enthalpies of formation for the species CF(3)CHCl(2), CF(3)CHFCl, CF(3)CCl(2), and CF(3)CFCl are evaluated at higher levels. With the information of the potential energy surface at BHLYP/6-311++G(3df,2p)//6-311G(d,p) level, we employ canonical variational transition-state theory with small-curvature tunneling correction to calculate the rate constants. The agreement between theoretical and experimental rate constants is good in the measured temperature range 276-382 K. The effect of fluorine substitution on reactivity of the C-H bond is discussed.     

Structures of water octamers (H2O)8: exploration on ab initio potential energy surfaces by the scaled hypersphere search method

The potential energy surface (PES) of water octamers has been explored by the scaled hypersphere search method. Among 164 minima on the PES (based on MP2/6-311++G(3df,2p)//B3LYP/6-311+G(d,p) calculations), the cubic structure with D2d symmetry has been confirmed to be the global minimum. In a thermodynamic simulation using these 164 structures, the cubic structure with S4 symmetry has the highest population at low temperature, though double rings can become dominant as temperature going up, in good accord with a recent Monte Carlo simulation using an empirical potential. A transition temperature from cubic to noncubic has significantly been underestimated when potential energy data of B3LYP/6-311+G(d,p) calculations are employed in the simulation. This serious discrepancy between the MP2 and the B3LYP results suggests an importance of dispersion interactions for discussions on thermodynamics of water octamers.     

On the flexibility of beta-peptides

The full conformational space was explored for an achiral and two chiral beta-peptide models: namely For-beta-Ala-NH2, For-beta-Abu-NH2, and For-beta-Aib-NH2. Stability and conformational properties of all three model systems were computed at different levels of theory: RHF/3-21G, B3LYP/6-311++G(d,p)//RHF/3-21G, B3LYP/6-311++G(d,p), MP2//B3LYP/6-311++G(d,p), CCSD//B3LYP/6-311++G(d,p), and CCSD(T)//B3LYP/6-311++G(d,p). In addition, ab initio E = E(phi, micro, psi) potential energy hypersurfaces of all three models were determined, and their topologies were analyzed to determine the inherent flexibility properties of these beta-peptide models. Fewer points were found and assigned than expected on the basis of Multidimensional Conformational Analysis (MDCA). Furthermore, it has been demonstrated, that the four-dimensional surface, E = E(phi, mu, psi), can be reduced into a three-dimensional one: E = E[phi, f(phi), psi]. This reduction of dimensionality of freedom of motion suggests that beta-peptides are less flexible than one would have thought. This agrees with experimental data published on the conformational properties of peptides composed of beta-amino acid residues.     

Conformational inversion-topomerization processes of ethylcyclohexane and 1,2-dimethylcyclohexane: A computational investigation

To deeply reveal the impact of the substituents and their special orientations in ring on conformational behaviors for substituted cyclohexanes, a comprehensive study of ethylcyclohexane, cis-, and trans-1,2-dimethylcyclohexanes has been carried out. All conformational structures for them were captured by the accurate ab intio method, that is, B3LYP/6-311++G(d,p) method was used for geometry optimizations, and MP2/6-311++G(d,p), G4, and CCSD(T)/6-311++G(d,p) methods were applied for the high-level single point energy refinements. Based on CCSD(T)/6-311++G(d,p) quantum results, the conformational populations of minima for these three substituted cyclohexanes were calculated by Boltzmann distribution over 300-2500 K. Additionally, the conformational inversion-topomerization pathways for them were thoroughly investigated. The complete characterization involved in their potential energy surfaces are clearly presented by three or two-dimensional schemes.     

Insight into global reaction mechanism of [C2, H4, O] system from ab initio calculations by the scaled hypersphere search method

A detailed computational study is performed on the singlet potential energy surface (PES) for possible isomerization and dissociation reactions of CH(3)CHO at the DFT (B3LYP/6-311++G(d,p)) and CCSD(T)/cc-pVTZ//B3LYP/6-311++G(d,p) levels. The pathways around the equilibrium structures can be discovered by the scaled hypersphere search (SHS) method, which enables us to make a global analysis of the PES for a given chemical composition. Fourteen isomers inclusive of 11 single-molecules and three "non-stabilized" oxygen-based ylides, 5 energetically favored complexes, and 79 interconversion transition states have been found on the singlet PES. Four lowest lying isomers with thermodynamic stability are also kinetically stable with respect to metastable intermediates. It was revealed that vinyl alcohols, which could be generated by the tautomerization of acetaldehyde, could undergo dissociation to form acetylene and water. In addition, recombination channels between some fragments, such as H(2)CO + (1)CH(2) and (1)CHOH + (1)CH(2), are energetically accessible via collision complex or oxygen-based ylides. Most of available unimolecular decompositions are found to be responsible for favorable hydrogen abstraction processes.     

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等价电子相同的分子的势能面进行了比较.}     

Quantum chemistry study of H+(H2O)8: a global search for its isomers by the scaled hypersphere search method, and its thermal behavior

The structures of the protonated water cluster H+(H2O)8 have been globally explored by the scaled hypersphere search method. On the Hartree-Fock potential energy surface 174 isomers were found, among which 168 were computed to be minima at the B3LYP/6-31+G** level, and their energies were further refined at the level of MP2/6-311++G(3df,2p). The global minimum on the potential energy surface computed at the B3LYP/6-31+G** level shows a cagelike structure with the "Eigen" motif, while the lowest-free-energy isomer has a five-membered-ring structure at 170 K and a chain form at 273 K. The present results are well in line with previous experimental findings. In addition, the ADMP (atom-centered density matrix propagation) simulation indicates that the extra proton in the lowest-free-energy isomer (170 K), which has a five-membered ring and the "Zundel" feature, is often in an asymmetrical hydrogen bond.