Clusters of hydrated methane sulfonic acid CH3SO3H.(H2O)n (n = 1-5): a theoretical study

Ab initio and density functional methods have been used to examine the structures and energetics of the hydrated clusters of methane sulfonic acid (MSA), CH3SO3H.(H2O)n (n = 1-5). For small clusters with one or two water molecules, the most stable clusters have strong cyclic hydrogen bonds between the proton of OH group in MSA and the water molecules. With three or more water molecules, the proton transfer from MSA to water becomes possible, forming ion-pair structures between CH3SO3- and H3O+ moieties. For MSA.(H2O)3, the energy difference between the most stable ion pair and neutral structures are less than 1 kJ/mol, thus coexistence of neutral and ion-pair isomers are expected. For larger clusters with four and five water molecules, the ion-pair isomers are more stable (>10 kJ/mol) than the neutral ones; thus, proton transfer takes place. The ion-pair clusters can have direct hydrogen bond between CH3SO3- and H3O+ or indirect one through water molecule. For MSA.(H2O)5, the energy difference between ion pairs with direct and indirect hydrogen bonds are less than 1 kJ/mol; namely, the charge separation and acid ionization is energetically possible. The calculated IR spectra of stable isomers of MSA.(H2O)n clusters clearly demonstrate the significant red shift of OH stretching of MSA and hydrogen-bonded OH stretching of water molecules as the size of cluster increases.     

Theoretical investigation of gas phase ethanol–(water)n (n = 1–5) clusters and comparison with gas phase pure water clusters (water)n (n = 2–6)

Various properties (such as optimal structures, structural parameters, hydrogen bonds, natural bond orbital charge distributions, binding energies, electron densities at hydrogen bond critical points, cooperative effects, and so on) of gas phase ethanol–(water)_n (n = 1–5) clusters with the change in the number of water molecules have been systematically explored at the MP2/aug‐cc‐pVTZ//MP2/6‐311++G(d,p) computational level. The study of optimal structures shows that the most stable ethanol‐water heterodimer is the one where exists one primary hydrogen bond (O? H…O) and one secondary hydrogen bond (C? H …O) simultaneously. The cyclic geometric pattern formed by the primary hydrogen bonds, where all the molecules are proton acceptor and proton donor simultaneously, is the most stable configuration for ethanol–(water)_n (n = 2–4) clusters, and a transition from two‐dimensional cyclic to three‐dimensional structures occurs at n = 5. At the same time, the cluster stability seems to correlate with the number of primary hydrogen bonds, because the secondary hydrogen bond was extremely weaker than the primary hydrogen bond. Furthermore, the comparison of cooperative effects between ethanol–water clusters and gas phase pure water clusters has been analyzed from two aspects. First of all, for the cyclic structure, the cooperative effect in the former is slightly stronger than that of the latter with the increasing of water molecules. Second, for the ethanol–(water)₅ and (water)₆ structure, the cooperative effect in the former is also correspondingly stronger than that of the latter except for the ethanol–(water)₅ book structure. © 2012 Wiley Periodicals, Inc.     

Exploration of the (ethanol)4–water heteropentamers potential energy surface by simulated annealing and Ab initio molecular dynamics

A DFT electronic structure study of the (ethanol)₄–water heteropentamers at the B3LYP/6‐31+G(d) model chemistry was carried out. To get determine possible configurations, the potential energy surface (PES) was explored with two methods: simulated annealing and ab initio molecular dynamics. The results suggest that the PES is very flat. A total of 81 stable structures were determined. These structures were classified into 16 different geometric patterns according to geometric criteria like the number of hydrogen bonds and their spatial arrangement: cyclic, bicyclic, or lineal patterns. Thermodynamic stability was used for defining the order of such classification. Hydrogen bonds are mutually disturbed due to the existence of cooperative effects. Cooperativity affects the nature of the hydrogen bonds and the overall stability of the ethanol–water system given that the strongest interactions are markedly covalent and the most stable geometric pattern corresponds to the pentagonal arrangement. These observations were supported by the analysis of the loss of atomic charge of the hydrogen atoms involved in hydrogen bonds. These hydrogen bonds were classified as primary and secondary hydrogen bonds: O? H ··· O and C? H ··· O, respectively. For comparative purposes, some (ethanol)₅, (methanol)₅, and (methanol)₄–water clusters were characterized in this study. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Dissociation chemistry of hydrogen halides in water

To understand the mechanism of aqueous acid dissociation, which plays a fundamental role in aqueous chemistry, the ionic dissociation of HX acids (X=F, Cl, Br, and I) in water clusters up to hexamer is examined using density functional theory and M?ller-Plesset second-order perturbation methods (MP2). Further accurate analysis based on the coupled clusters theory with singles and doubles excitations agrees with the MP2 results. The equilibrium structures, binding energies, electronic properties, stretching frequencies, and rotational constants of HX(H(2)O)(n) and X(-)(H(3)O)(+)(H(2)O)(n-1) are calculated. The dissociated structures of HF and HCl can be formed for n>/=4, while those of HBr and HI can be formed for n>/=3. Among these, the dissociated structures of HX (X=Cl, Br, and I) are more stable than the undissociated ones for n>/=4, while such cases for HF would require much more than six water molecules, in agreement with previous reports. The IR spectra of stable clusters including anharmonic frequencies are predicted to facilitate IR experimental studies. Undissociated systems have X-H stretching modes which are highly redshifted by hydration. Dissociated hydrogen halides show three characteristic OH stretching modes of hydronium moiety, which are redshifted from the OH stretching modes of water molecules.     

Ab initio investigation of the lower-energy candidate structures for (H<Subscript>2</Subscript>O)<Subscript>10</Subscript><Superscript>+</Superscript> water cluster

Low-lying structures of water cationic clusters and the compounds with the OH radical have become a hot topic in recent years. We here investigate the cluster \( {\left({\mathrm{H}}_2\mathrm{O}\right)}_{10}^{+} \) and calculate its ideal structures by the quantum chemical calculation together with the particle swarm optimization method. We analyzed the properties of the obtained lower-energy isomers of \( {\left({\mathrm{H}}_2\mathrm{O}\right)}_{10}^{+} \). Their energies are further re-optimized and demonstrated at three different methods with two basis sets. Based on our numerical calculations, a new cage-like structure of \( {\left({\mathrm{H}}_2\mathrm{O}\right)}_{10}^{+} \) with the lowest energy is obtained at MP2/aug-cc-pVDZ level. Our results showed the comparison of energy order at different conditions and demonstrated the influence of temperature on the relative Gibbs energy and IR spectra. Moreover, we also contained the molecule orbitals to discuss the stability of these representative isomers.     

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala₂₅). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala₂₅. DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅. This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala₂₅. Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.     

In search of CS2(H2O)(n=1-4) clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS(2)(H2O)n, where n=1-4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS(2)(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS(2)(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS(2) in the atmosphere, the concentrations of these clusters are predicted to be on the order of 10(5) CS(2)(H2O) clusters.cm(-3) and 10(2) CS(2)(H2O)(2) clusters.cm(-3) at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS(2)(H(2)O)(n) structures predicted to be present in a low-temperature molecular beam experiment.     

High-level ab initio calculations for the four low-lying families of minima of (H2O)20. I. Estimates of MP2/CBS binding energies and comparison with empirical potentials

We report estimates of complete basis set (CBS) limits at the second-order M?ller-Plesset perturbation level of theory (MP2) for the binding energies of the lowest-lying isomers within each of the four major families of minima of (H(2)O)(20). These were obtained by performing MP2 calculations with the family of correlation-consistent basis sets up to quadruple zeta quality, augmented with additional diffuse functions (aug-cc-pVnZ, n=D, T, Q). The MP2/CPS estimates are -200.1 (dodecahedron, 30 hydrogen bonds), -212.6 (fused cubes, 36 hydrogen bonds), -215.0 (face-sharing pentagonal prisms, 35 hydrogen bonds), and -217.9 kcal/mol (edge-sharing pentagonal prisms, 34 hydrogen bonds). The energetic ordering of the various (H(2)O)(20) isomers does not follow monotonically the number of hydrogen bonds as in the case of smaller clusters such as the different isomers of the water hexamer. The dodecahedron lies ca. 18 kcal/mol higher in energy than the most stable edge-sharing pentagonal prism isomer. The TIP4P, ASP-W4, TTM2-R, AMOEBA, and TTM2-F empirical potentials also predict the energetic stabilization of the edge-sharing pentagonal prisms with respect to the dodecahedron, albeit they universally underestimate the cluster binding energies with respect to the MP2/CBS result. Among them, the TTM2-F potential was found to predict the absolute cluster binding energies to within <1% from the corresponding MP2/CBS values, whereas the error for the rest of the potentials considered in this study ranges from 3% to 5%.     

Theoretical calculation of structures and proton transfer in hydrated ammonia-hydrogen chloride clusters

Ab initio molecular orbital calculations have been performed to investigate the structures and quantum effects of the proton motion in NH(3):HCl:(H(2)O)(n) (n = 0-3) clusters using a MP2/aug-cc-pVDZ level of theory. Three new stable structures and one transition-state structure are investigated for these clusters. The detailed analyses of the intermolecular interactions suggest that three-body interactions play an important role to determine the relative stability in each size of cluster. The quantum effects of the proton motion result in frequency shifts for proton-stretching modes. Our one-dimensional and two-dimensional models fairly closely reproduce the experimental proton-stretching vibrational frequency of the NH(3):HCl cluster. The most stable isomer for n = 1 has a proton-transfer structure, which is weakened by the quantum effects of the proton motion.     

Geometries and stabilities of various configurations of benzene dimer: details of novel V-shaped structure revealed

Benzene dimer configurations namely T-shaped, parallel-displaced, sandwich, and V-shaped that were proposed by experimental studies are investigated using second- and fourth-order Møller–Plesset perturbation theory. The MP2 method with aug-cc-pVDZ and aug-cc-pVTZ basis sets unequivocally shows that the parallel-displaced configuration is considerably more stable than T-shaped structure. On the other hand, the MP4(SDTQ)/aug-cc-pVDZ level predicts that the T-shaped and parallel-displaced configurations are nearly isoenergetic, which is parallel to the previous results of estimated CCSD(T)/CBS level reported recently. The lowest energy T-shaped configuration is stabilized by 0.17 kcal/mol over the parallel-displaced configuration at the MP4(SDTQ)/aug-cc-pVDZ level. Although the structures of all the four different types of configurations are found to be stable at both MP2 and full MP4 methods, the V-shaped configuration is the least stable among them. The calculated interaction energy of ?2.3 kcal/mol for the lowest energy T-shaped structure at the MP4(SDTQ)/aug-cc-pVDZ level is in good agreement with the experimental value of ?2.4 ± 0.4 kcal/mol. We conclude that the MP4(SDTQ) with a reasonably good basis set can be used for systems involving π–π interactions to obtain qualitative and quantitative results.     

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.     

Ring-Stacking Water Clusters: Morphology and Stabilities

The structures and interaction energies of water clusters with ring stacking motifs are studied by using ab initio calculations. The structures of the water clusters are constructed by stacking either single rings or multi-rings of tetramer, pentamer, and hexamer. We found that, in the single-ring-stacking motif, the most stable isomers exhibit an alternative clockwise-anticlockwise stacking pattern. We also show that four-layer single-ring-stacking isomers are not energetically favorable in comparison with those of two-layer multi-ring-stacking isomers. The relative stability of the isomers is also analyzed in terms of H-bond strength and elastic distortions of the water molecules.     

Electron propagator theory study of N-/O-methylglycine conformers

The low-lying conformers of N-/O-methylglycine are studied by ab initio calculations at the B3LYP, MP3, and MP4(SDQ) levels of theory with the aug-cc-pVDZ basis set. The conformers having the intramolecular hydrogen bonds N-H...O=C or O-H...N are more stable than the others. Vertical ionization energies for the valence molecular orbitals of each conformer predicted with the electron propagator theory in the partial third-order quasiparticle approximation are in good agreement with the experimental data available in the literatures. The relative energies of the conformers and comparison between the simulated and the experimental photoelectron spectra demonstrate that there are at least three and two conformers of N- and O-methylglycine, respectively, in the gas-phase experiments. The intramolecular hydrogen bonding O-H...N effects on the molecular electronic structures are discussed for the glycine methyl derivatives, on the basis of the ab initio electronic structure calculations, natural orbital bond, and atoms-in-molecules analyses. The intramolecular hydrogen bonding O-H...N interactions hardly affect the electronic structures of the O-NH2-CH2-C(=O)-O-CH3 and alpha-methylated NH2-CH2-C(CH3)OOH conformers, while the similar intramolecular interactions lead to the significantly lower-energy levels of the highest occupied molecular orbitals for the N-(CH3-NH-CH2-COOH) and beta-methylated (NH2-CH2-CH2-COOH) conformers.     

Characterization of solvated electrons in hydrogen cyanide clusters: (HCN)n- (n=3, 4)

Theoretical studies of the solvated electrons (HCN)n- (n=3, 4) reveal a variety of electron trapping possibilities in the (HCN)n (n=3, 4) clusters. Two isomers for (HCN)3- and four isomers for (HCN)4- are obtained at the MP2aug-cc-pVDZ+dBF (diffusive bond functions) level of theory. In view of vertical electron detachment energies (VDEs) at the CCSD(T) level, the excess electron always "prefers" locating in the center of the system, i.e., the isomer with higher coordination number shows larger VDE value. However, the most stable isomers of the solvated electron state (HCN)3- and (HCN)4- are found to be the linear Cinfinitynu and Dinfinityh structures, respectively, but not the fullyl symmetric structures which have the largest VDE values.     

Theoretical studies on photoelectron and IR spectral properties of Br2.-(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Becke's half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.     

Furan-formic acid dimers: an ab initio and matrix isolation study

The dimers formed by formic acid (FA) and furan are investigated by ab initio methods and matrix isolation spectroscopy. Nine complexes with binding energies between -3.91 and -0.82 kcal/mol (MP2/6-311++G(d,p) + ZPE + BSSE) are identified. Another five weaker bound complexes are localized at lower level of theory only. The binding in the furan-FA dimers can be described in terms of OH...O, C=O...H, HO...H, CH...O, OH...pi, and CH...pi interactions. Therefore, the furan-FA complexes are classified in two types: (1) the dimers where the OH hydrogen of formic acid interacts with the furan molecule and (2) the dimers where the main interactions of FA with the furan molecule are via the less acidic CH hydrogen. Duning's and Pople's triple and double basis sets were used to study the dependence of the geometries and energies of the complexes from the basis set. BSSE (basis set superposition error) counterpoise corrections (CP) were included during the geometry optimizations of all dimers at the MP2/6-31G(d,p) level of theory. Matrix isolation spectroscopy allowed us to record the IR spectrum of aggregates between FA and furan. By comparison of the experimental IR spectrum with calculated IR spectra of a variety of complexes, it was possible to identify the most stable furan-FA dimer as the major product of the aggregation.     

High-level ab initio calculations for the four low-lying families of minima of (H2O)20. II. Spectroscopic signatures of the dodecahedron, fused cubes, face-sharing pentagonal prisms, and edge-sharing pentagonal prisms hydrogen bonding networks

We report the first harmonic vibrational spectra for each of the lowest lying isomers within the four major families of minima of (H2O)20, namely, the dodecahedron, fused cubes, face-sharing pentagonal prisms, and edge-sharing pentagonal prisms. These were obtained at the second-order Moller-Plesset perturbation level of theory (MP2) with the augmented correlation consistent basis set of double zeta quality (aug-cc-pVDZ) at the corresponding minimum energy geometries. The computed infrared (IR) spectra are the first ones obtained from first principles for these clusters. They were found to contain spectral features, which can be directly mapped onto the distinctive spectroscopic signatures of their constituent tetramer, pentamer, and octamer fragments. The dodecahedron spectra show the richest structure in the OH stretching region and are associated with the most redshifted OH vibrations with respect to the monomer. The lowest lying edge-sharing pentagonal prism isomer displays intense IR active vibrations that are redshifted by approximately 600 cm(-1) with respect to the water monomer. Furthermore the most redshifted, IR-active OH stretching vibrations for all four networks correspond to hydrogen bonded OH groups, which exhibit the following two common characteristics: (i) they belong to fragments which have a "free" OH stretch and (ii) they act as donors to a neighboring water molecule along a "dimerlike" (strong) hydrogen bond. The zero-point energy corrected MP2/CBS (complete basis set) limit binding energies D(0) for the four isomers are -163.1 kcal/mol (edge-sharing pentagonal prism), -160.1 kcal/mol (face-sharing pentagonal prism), -157.5 kcal/mol (fused cubes), and -148.1 kcal/mol (dodecahedron).     

甲醇团簇的多光子电离质谱及其从头算

用355 nm激光对脉冲分子束超声膨胀冷却的甲醇分子进行多光子电离, 飞行时间质谱仪观测到除甲醇碎片离子外的质子化甲醇团簇(CH3OH)nH+(n=1-16), 且离子的种类及相对强度与激光相对于脉冲分子束的延时无关, 取决于团簇离子内在结构的稳定性. 结合从头算密度泛函理论, 在B3LYP/6-31G(d)基组水平上优化得到了(CH3OH)n和(CH3OH)nH+(n=1-4)的稳定构型. 振动频谱分析显示, 团簇中最强的红外振动模主要来自氢键H伸缩振动的贡献. 团簇电离后发生于团簇内的质子转移反应也可能与激光电离引起的与氢键有关的振动模激发密切相关.     

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.     

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等价电子相同的体系进行了比较。