Potential energy surface for activation of methane by Pt(+): a combined guided ion beam and DFT study

A guided-ion beam tandem mass spectrometer is used to study the reactions of Pt(+) with methane, PtCH(2)(+) with H(2) and D(2), and collision-induced dissociation of PtCH(4)(+) and PtCH(2)(+) with Xe. These studies experimentally probe the potential energy surface for the activation of methane by Pt(+). For the reaction of Pt(+) with methane, dehydrogenation to form PtCH(2)(+) + H(2) is exothermic, efficient, and the only process observed at low energies. PtH(+), formed in a simple C-H bond cleavage, dominates the product spectrum at high energies. The observation of a PtH(2)(+) product provides evidence that methane activation proceeds via a (H)(2)PtCH(2)(+) intermediate. Modeling of the endothermic reaction cross sections yields the 0 K bond dissociation energies in eV (kJ/mol) of D(0)(Pt(+)-H) = 2.81 +/- 0.05 (271 +/- 5), D(0)(Pt(+)-2H) = 6.00 +/- 0.12 (579 +/- 12), D(0)(Pt(+)-C) = 5.43 +/- 0.05 (524 +/- 5), D(0)(Pt(+)-CH) = 5.56 +/- 0.10 (536 +/- 10), and D(0)(Pt(+)-CH(3)) = 2.67 +/- 0.08 (258 +/- 8). D(0)(Pt(+)-CH(2)) = 4.80 +/- 0.03 eV (463 +/- 3 kJ/mol) is determined by measuring the forward and reverse reaction rates for Pt(+) + CH(4) right harpoon over left harpoon PtCH(2)(+) + H(2) at thermal energy. We find extensive hydrogen scrambling in the reaction of PtCH(2)(+) with D(2). Collision-induced dissociation (CID) of PtCH(4)(+), identified as the H-Pt(+)-CH(3) intermediate, with Xe reveals a bond energy of 1.77 +/- 0.08 eV (171 +/- 8 kJ/mol) relative to Pt(+) + CH(4). The experimental thermochemistry is favorably compared with density functional theory calculations (B3LYP using several basis sets), which also establish the electronic structures of these species and provide insight into the reaction mechanism. Results for the reaction of Pt(+) with methane are compared with those for the analogous palladium system and the differences in reactivity and mechanism are discussed.     

Ab initio studies on mechanisms of cycloaddition reaction between dichloro-vinylidene and ethylene

Mechanisms of cycloaddition reaction between singlet dichloro-vinylidene (R1) and ethylene (R2) have been investigated with MP2 and B3LYP /6-31G* methods, including geometry optimization, vibrational analysis, and energy for the involved stationary points on the potential energy surface. CCSD(T)/6-31G* single point calculations are also applied to the geometries from both methods. CCSD(T) relative energies for the stationary points predicted by MP2 and B3LYP are in a very good agreement. Three reaction pathways are located. Along the first one, one intermediate (INT1) is firstly generated, which then rearranges into a three-membered ring compound (P1) via a small barrier of 5.4 kJ/mol; the other two paths share the other intermediate (INT2), which isomerizes into two four-membered ring compounds (P2 and P3) via a chlorine and a hydrogen transfer, respectively. The barriers of the latter two paths are significantly higher by approximately 25 kJ/mol than that of the former (27.2 and 28.8 vs 5.4 kJ/mol), the major reaction is therefore the formation of P1.     

Experimental and theoretical studies of potassium cation interactions with the acidic amino acids and their amide derivatives

The binding of K(+) to aspartic acid (Asp), glutamic acid (Glu), asparagine (Asn), and glutamine (Gln) is examined in detail by studying the collision-induced dissociation (CID) of the four potassium cation-bound amino acid complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Formed by electrospray ionization, these complexes have energy-dependent CID cross sections that are analyzed to provide 0 K bond energies after accounting for unimolecular decay rates, internal energy of reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each K(+)(L) complex are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at B3LYP, B3P86, and MP2(full) levels using a 6-311+G(2d,2p) basis set. Theoretical bond dissociation energies are in good agreement with the experimental values. This coordinated examination of both experimental work and quantum chemical calculations allows for a comprehensive understanding of the molecular interactions of K(+) with the Asx and Glx amino acids. K(+) binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 9+/-1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the K(+) binding affinity for the longer-chain amino acids (Glx) is enhanced by 5+/-1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced. Further, a detailed comparison between experimental and theoretical results reveals interesting differences in the binding of K(+) and Na(+) to these amino acids.     

Quantum mechanical investigation of the inner-sphere reorganization energy of cyclooctatetraene/cyclooctatetraene radical anion. Part I

The inner-sphere reorganization energy of the electron self-exchange of the couple cyclooctatetraene/cyclooctatetraene radical anion has been investigated by quantum mechanical calculations. The more stable Jahn Teller distorted B2g conformation of the radical anion has been used in this study. Two different theories have been applied in this first part. The harmonic approximation in the classical Marcus scheme has been modified by using projected force constants, which are obtained from the complete force constant matrix and the geometry changes of the molecule during the ET (introduced by Mikkelsen). A different approach (introduced by Nelsen) combines the different energies of the neutral and radical anion with and without relaxation corresponding to the vertical ionization potential and the vertical electron affinity. The electronic energies of the neutral molecule and the radical anion differ dramatically applying three different levels of quantum mechanical calculations (UAM1, UB3LYP, PMP2 with three different basis sets with and without diffuse functions). Nevertheless the Nelsen method gives almost consistent results for the inner-sphere reorganization energies: 120.1 kJ/mol for semiempirical UAM1 method, 159.3 kJ/mol, 156.4 kJ/mol and 158.3 kJ/mol for density functional UB3LYP/6-31G*, UB3LYP/6-31++G* and UB3LYP/AUG-cc-pVDZ calculations and 192.5 kJ/mol for ab-initio PMP2/6-31G* investigations, respectively. These values are in agreement with earlier experimental work supposing the total reorganization energy to be larger than 38 kcal/mol assuming an electron self-exchange rate of 10(4) M(-1) s(-1). The simple harmonic approximation of Marcus relation has not yet been applied for a molecule like cyclooctatetraene with large torsional geometry changes. Using the projected force constants after scaling, considerably different results for the inner-sphere reorganization energy have been calculated: 738.1 kJ/mol for the UB3LYP/6-31G*, 743.3 kJ/mol for UB3LYP/6-31++G* and 759.1 kJ/mol for UB3LYP/AUG-cc-pVDZ level of theory. Comparison with our concentration dependent EPR experiments are controversial to the earlier experimental results, but the latter supports the assumption that the electron self-exchange occurs in a time scale so that the molecules cannot complete their vibrational motions. Therefore the projected Marcus relation is not valid for cyclooctatetraene/cyclooctatetraene radical anion including a large torsional change during the electron transfer.     

Mechanistic study of the deamination reaction of guanine: a computational study

The mechanism for the deamination of guanine with H(2)O, OH(-), H(2)O/OH(-) and for GuaH(+) with H(2)O has been investigated using ab initio calculations. Optimized geometries of the reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), B3LYP/6-31G(d), and B3LYP/6-31+G(d) levels of theory. Energies were also determined at G3MP2, G3MP2B3, G4MP2, and CBS-QB3 levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Thermodynamic properties (ΔE, ΔH, and ΔG), activation energies, enthalpies, and Gibbs free energies of activation were also calculated for each reaction investigated. All pathways yield an initial tetrahedral intermediate and an intermediate in the last step that dissociates to products via a 1,3-proton shift. At the G3MP2 level of theory, deamination with OH(-) was found to have an activation energy barrier of 155 kJ mol(-1) compared to 187 kJ mol(-1) for the reaction with H(2)O and 243 kJ mol(-1) for GuaH(+) with H(2)O. The lowest overall activation energy, 144 kJ mol(-1) at the G3MP2 level, was obtained for the deamination of guanine with H(2)O/OH(-). Due to a lack of experimental results for guanine deamination, a comparison is made with those of cytosine, whose deamination reaction parallels that of guanine.     

尿酸质子转移反应的理论研究I——分子内质子转移研究

使用量子化学中的Hartree-Fock方法和密度泛函理论中的B3LYP方法，分别在3－21G^*和6－31G（d)水平上，计算了尿酸分子从三羰基异构体向三羟基异构体的转化。结果表明，转化过程经历了单羟基和双羟基异构体2种中间物和3种过渡态时的分子内质子转移（IPT），转移中的H原邻近的N，O和C原子形成了具有四元环结构的过渡态。随着IPT的进行，N－H键逐渐被削弱和断裂，O－H键则逐渐生成。3个反应的活化能分别为190.3kJ/mol,181.4kJ/mol和249.9kJ/mol(B3LYP/6-31G(d))。较高的活化能表明在室温下，无催化剂的IPT难以进行。     

F原子与瞬态自由基CH_2SH反应的理论研究

用量子化学从头算和密度泛函理论（DFT）对F原子与自由基CH_2SH在势能面上 的反应进行了研究。在B3LYP/6-311G水平上计算出了各物种的优化构型、振动频率 和零点振动能(ZPVE);各物种的总能量由B3LYP/6-311 + G(2df, pd)//B3LYP/6- 311G计算,另外对反应物和产物还计算了其G3能量。结果表明：首先F通过与C或S 结合的两种途径与CH_2SH相配位,再通过H（4）原子转移形成甲基,然后甲基再旋 转,甲基中H（4）原子最终与F结合,反应产物为HF和CH_2S。反应为放热反应,分 别为ΔH_r = -370.7 kJ/mol (DFT)和-396.94 kJ/mol (G3)。此外依据计算出的反 应热,可得自由基·CH_2SH的生成热Δ_fH°_(298.15) = 146.44 kJ/mol (DFT), 而Δ_fH°_0 = 167.36 kJ/mol (G3)。它们与以前的实验和理论值是一致的。     

Theoretical study on the nitration of methane by acyl nitrate catalyzed by H-ZSM5 zeolite

A theoretical study on the nitration of methane by acyl nitrate catalyzed by HZSM-5 zeolite is reported. The zeolite was represented by a "double ring" 20T cluster. The calculations were performed at the DFT/X3LYP/6-31G** and MP2/6-31G** levels. The first step of the mechanism involves the protonation of the acyl nitrate by the zeolite and the formation of a nitronium-like ion. The reaction proceeds through a concerted step with the attack of the methane molecule by the nitronium-like ion and the simultaneous transfer of a proton from the methane molecule to the zeolite, thus reconstructing the acidic site. The activation energies for the first and second steps of this reaction are, respectively, 14.09 and 10.14 kcal/mol at X3LYP/6-31G** level and 16.68 and 13.85 kcal/mol at the MP2/6-31G**.     

Structures of cationized proline analogues: evidence for the zwitterionic form

The structures of lithiated and sodiated alpha-methyl-proline (alpha-Me-Pro) and structural isomers, both with and without a water molecule, are investigated using blackbody infrared radiative dissociation (BIRD) and density functional theory. From the BIRD kinetic data measured as a function of temperature, combined with master equation modeling of these data, threshold dissociation energies for the loss of a water molecule from these clusters are obtained. These energies are 77.5 +/- 0.5 and 53 +/- 1 kJ/mol for lithiated and sodiated alpha-Me-Pro, respectively. For the nonzwitterionic isomer, proline methyl ester, these values are 3.0-4.5 kJ/mol higher. These results provide compelling experimental evidence that alpha-Me-Pro is zwitterionic in these clusters. Theory at the temperature corrected B3LYP/6-311++G**//B3LYP/6-31++G** level indicates that the salt-bridge or zwitterionic forms of lithiated and sodiated alpha-Me-Pro are between 17 and 23 kJ/mol lower in energy than the nonzwitterionic or charge-solvated forms and that attachment of a single water molecule does not significantly change the structure or the relative energies of these clusters. The proton affinity of proline is 8 kJ/mol higher than that of alpha-Me-Pro, indicating that lithiated and sodiated singly hydrated proline should also be zwitterionic.     

Accurate measurement of the relative bond energies of CO and H2O ligands in Fe+ mono- and bis-ligated complexes

The systems Fe(H(2)O)(n) (+)/CO[bond]H(2)O and Fe(CO)(n) (+)/CO[bond]H(2)O (n = 1 and 2) were investigated in a triple cell Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. Using mixtures of CO with a very small amount of water, the ligand exchange equilibrium was reached, allowing experimental determination of the relevant equilibrium constants and free energies of reaction. Quantum chemical calculations at the B3LYP level of theory on the reactant and product species allowed us to determine the entropic terms and to derive the relative bond energies of CO and H(2)O in the mono- and bis-ligated complexes. For n = 1, H(2)O is more strongly bound to Fe(+) than CO by 4.1 +/- 1.6 kJ x mol(-1) at 298 K. For n = 2, at the same temperature, H(2)O is more strongly bound than CO to (H(2)O)Fe(+) by 7.6 +/- 1.6 kJ x mol(-1), and to (CO)Fe(+) by more than 20.1 kJ x mol(-1).     

Structure and reactivity of V(2)O(5): bulk solid, nanosized clusters, species supported on silica and alumina, cluster cations and anions

Vanadyl bond dissociation energies are calculated by density functional theory (DFT). While the hybrid (B3LYP) functional results are close to the available reference data, gradient corrected functionals (BP86, PBE) yield large errors (about 50 to 100 kJ mol(-1)), but reproduce trends correctly. PBE calculations on a V(20)O(62)H(24) cluster model for the (001) surface of V(2)O(5) crystals virtually reproduce periodic slab calculations. The low bond dissociation energy (formation of oxygen surface defect) of 113 kJ mol(-1)(B3LYP) is due to substantial structure relaxations leading to formation of V-O-V bonds between the V(2)O(5) layers of the crystal. This relaxation cannot occur in polyhedral (V(2)O(5))(n) clusters and also not for V(2)O(5) species supported on silica or alumina (represented by cage-type models) for which bond dissociation energies of 250-300 kJ mol(-1) are calculated. The OV(OCH(3))(3) molecule and its dimer are also considered. Radical cations V(2)O(5)(+) and V(4)O(10)(+) have very low bond dissociation energies (22 and 14 kJ mol(-1), respectively), while the corresponding radical anions have higher dissociation energies (about 330 kJ mol(-1)) than the neutral clusters. The bond dissociation energies of the closed shell V(3)O(7)(+) cation (165 kJ mol(-1)) and the closed shell V(3)O(8)(-) anion (283 kJ mol(-1)) are closest to the values of the neutral clusters. This makes them suitable for gas phase studies which aim at comparisons with V(2)O(5) species on supporting oxides.     

Bioinorganic chemistry modeled with the TPSSh density functional

In this work, the TPSSh density functional has been benchmarked against a test set of experimental structures and bond energies for 80 transition-metal-containing diatomics. It is found that the TPSSh functional gives structures of the same quality as other commonly used hybrid and nonhybrid functionals such as B3LYP and BP86. TPSSh gives a slope of 0.99 upon linear fitting to experimental bond energies, whereas B3LYP and BP86, representing 20% and 0% exact exchange, respectively, give linear fits with slopes of 0.91 and 1.07. Thus, TPSSh eliminates the large systematic component of the error in other functionals, reducing rms errors from 46-57 to 34 kJ/mol. The nonhybrid version of the functional, TPSS, gives a slope of 1.08, similar to BP86, implying that using 10% exact exchange is the main reason for the success of TPSSh. Typical bioinorganic reactions were then investigated, including spin inversion and electron affinity in iron-sulfur clusters, and breaking or formation of bonds in iron proteins and cobalamins. The results show that differences in reaction energies due to exact exchange can be much larger than the usually cited approximately 20 kJ/mol, sometimes exceeding 100 kJ/mol. The TPSSh functional provides energies approximately halfway between nonhybrids BP86 and TPSS, and 20% exact exchange hybrid B3LYP: Thus, a linear correlation between the amount of exact exchange and the numeric value of the reaction energy is observed in all these cases. For these reasons, TPSSh stands out as a most promising density functional for use and further development within the field of bioinorganic chemistry.     

Ab initio study of bonding between nucleophilic species and the nitroso group

The bonding between anionic nucleophiles and the nitroso group has been studied in the common nitrosating agents nitroso chloride (ONCl), nitroso bromide (ONBr), nitroso thiocyanate (ONSCN), and dinitrogen trioxide (N2O3) in aqueous solution. A variety of theoretical methods were employed, including ab initio, density functional theory (DFT), and composite theoretical techniques, with solvent effects described using the polarizable continuum model (PCM). Experimental nitroso bond heterolytic dissociation free energies were accurately reproduced with a number of composite theoretical methods, the most successful being CBS-Q and G2MP2, with average errors of 3.1 and 3.4 kJ mol(-1), respectively. Using the MP2 and B3LYP methods, calculations were made with correlation consistent basis sets up to quadruple-zeta, extrapolated to the complete basis set (CBS) limit. The MP2/CBS calculations were accurate to around 10 kJ mol(-1), while the B3LYP/CBS calculations routinely overpredicted experimental bond free energies by ca. 40 kJ mol(-1). It is therefore highly recommended that B3LYP energies are not used for nitroso compounds, although other results demonstrate that the B3LYP method provides a good account of nitroso compound geometries, frequencies, and entropies. Single-point CBS energy calculations using MP2/aug-cc-pVQZ geometries and frequencies showed that the MP4(SDTQ) and QCISD(T) methods provide a slight improvement over MP2 at the CBS limit, although the inclusion of triple excitations is necessary to achieve this improvement in accuracy. Enthalpy-entropy compensation was also discovered, with an average isoequilibrium temperature of 825 K. This relatively large isoequilibrium temperature indicates that enthalpic effects dominate over entropic ones.     

A Theoretical Study on the Equilibrium Structures and Relative Stabilities of H_2SO

1INTRODUCTIONTheconformationandvariousphysicalprope-rtiesofdichalcogencompoundsR-X-Y-R?(X,Y=O,S,Se,andTe)havebeenatopicofseveralexperi-mentalandtheoreticalinvestigations[1].CompoundscontainingS=SandS=Obondshavelongbeenpro-posedasintermediatesinorganicsynthesisand,onoccasion,asstableentities[2～18].Despitethecom-monlyusedrepresentationofS–Obondinsulfoxi-desandothersulfinylderivativesasS=O,sulfoxidesareinmanywaysbestdescribedasylideswithhighlypolarizedS–Oσ-bondbecauseofelectro-sta…     

单分子水对CH_2SH+NO_2反应机理影响的理论研究

在B3LYP/6-311++G(2df,p)水平下对单分子水参与下的CH_2SH+NO_2反应的微观机理进行了研究.为了获得更准确的能量信息,采用HL复合方法和CCSD(T)/aug-ccpvtz方法进行单点能校正.结果表明,加入单分子水后的CH_2SH+NO_2反应体系,共经过10条不同的反应路径,得到6种反应产物.与裸反应(CH_2SH+NO_2)相比,水分子在反应中起到了明显的正催化作用.不仅使生成产物trans-HONO的能垒(-52.84kJ·mol~(-1))降低了176.94kJ·mol~(-1),而且不需经过复杂的重排和异构化过程便可得到产物cis-HONO.在生成产物cis-HONO通道(Path3和Path4)中,活化能垒分别为143.65和126.70kJ·mol~(-1),而其裸反应的活化能垒却高达238.34kJ·mol~(-1).生成HNO_2的通道中(Path5和Path6)活化能垒分别为295.23和-42.19kJ·mol~(-1).其中Path6的无势垒过程使HNO_2也成为该反应的主要产物.另外,单分子水还可通过氢迁移的方式直接参与CH_2SH+NO_2的反应,活化能垒(TS7-TS10)分别为-10.62,151.03,186.22和155.10kJ·mol~(-1).除直接抽氢通道中的(Path8-Path10)外,其余反应通道均为放热反应,在热力学上是可行的.     

Molecular structure of the trans and cis isomers of metal-free phthalocyanine studied by gas-phase electron diffraction and high-level quantum chemical calculations: NH tautomerization and calculated vibrational frequencies

The molecular structure of the trans isomer of metal-free phthalocyanine (H2Pc) is determined using the gas electron diffraction (GED) method and high-level quantum chemical calculations. B3LYP calculations employing the basis sets 6-31G**, 6-311++G**, and cc-pVTZ give two tautomeric isomers for the inner H atoms, a trans isomer having D2h symmetry and a cis isomer having C2v symmetry. The trans isomer is calculated to be 41.6 (B3LYP/6-311++G**, zero-point corrected) and 37.3 kJ/mol (B3LYP/cc-pVTZ, not zero-point corrected) more stable than the cis isomer. However, Hartree-Fock (HF) calculations using different basis sets predict that cis is preferred and that trans does not exist as a stable form of the molecule. The equilibrium composition in the gas phase at 471 degrees C (the temperature of the GED experiment) calculated at the B3LYP/6-311++G** level is 99.8% trans and 0.2% cis. This is in very good agreement with the GED data, which indicate that the mole fraction of the cis isomer is close to zero. The transition states for two mechanisms of the NH tautomerization have been characterized. A concerted mechanism where the two H atoms move simultaneously yields a transition state of D2h symmetry and an energy barrier of 95.8 kJ/mol. A two-step mechanism where a trans isomer is converted to a cis isomer, which is converted into another trans isomer, proceeds via two transition states of C(s) symmetry and an energy barrier of 64.2 kJ/mol according to the B3LYP/6-311++G** calculation. The molecular geometry determined from GED is in very good agreement with the geometry obtained from the quantum chemical calculations. Vibrational frequencies, IR, and Raman intensities have been calculated using B3LYP/6-311++G**. These calculations indicate that the molecule is rather flexible with six vibrational frequencies in the range of 20-84 cm(-1) for the trans isomer. The cis isomer might be detected by infrared matrix spectroscopy since the N-H stretching frequencies are very different for the two isomers.     

Conformers and intramolecular hydrogen bonding of the oxalic acid monomer and its anions

Density functional theory, B3LYP/6‐31G** and B3LYP/6‐311+G(2d,p), and ab initio MP2/6‐31G** calculations have been carried out to investigate the conformers, transition states, and energy barriers of the conformational processes of oxalic acid and its anions. QCISD/6‐31G** geometrical optimization is also performed in the stable forms. Its calculated energy differences between the two most stable conformers are very near to the related observed value at 7.0 kJ/mol. It is found that the structures and relative energies of oxalic acid conformers predicted by these methods show similar results, and that the conformer L1 (C_2h) with the double‐interfunctional‐groups hydrogen bonds is the most stable conformer. The magnitude of hydrogen bond energies depends on the energy differences of various optimized structures. The hydrogen bond energies will be about 32 kJ/mol for interfunctional groups, 17 kJ/mol for weak interfunctional groups, 24 kJ/mol for intra‐COOH in (COOH)₂, and 60 kJ/mol for interfunctional groups in (COOH)COO⁻¹ ion if calculated using the B3LYP/6‐311+G(2d,p) method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 541–551, 2000     

Theoretical study on keto–enol tautomerism and isomerization in pyruvic acid

Isomerization and tautomerism of 12 isomers of pyruvic acid including 4 keto and 8 enol forms were studied at the MP2 and B3LYP levels of theory using 6‐311++G(2df,p) basis set, separately. Activation energy (E_a), imaginary frequency (υ), and Gibbs free energy (ΔG^#) of the considered isomerization and tautomerism reactions were calculated. Interconversion of the enol forms proceeds through two paths: (i) proton transfer and (ii) internal rotation. Activation energies for the proton transfer paths were in the range of 125–145 kJ/mol and for the internal rotation paths were in the range of 5–45 kJ/mol. Keto–enol tautomerism of pyruvic acid proceeds only through proton transfer route and their activation energies were in the range of 200–300 kJ/mol. Effect of microhydration on the transition state structures and activation energies was also investigated. It was found that the presence of a water molecule catalyzes the isomerization and tautomerism reactions of pyruvic acid so that the activation energies decrease. © 2013 Wiley Periodicals, Inc.     

苯并氧化呋咱稳定性和异构化的DFT和ab initio研究

运用B3LYP/6-31G(d)密度泛函理论（DFT）方法对苯并氧化呋咱、邻二亚硝基苯及其间的异构化反应进行了计算研究。结果表明,苯并氧化呋咱的分子总能量比邻二亚硝基苯的低;由苯并氧化呋咱异构为邻二亚硝基苯的正向反应活化能（Ea+=51.0kJ/mol）,与文献实测值（58.6kJ/mol）较接近,而其逆向反应活化能（Ea-=4.6kJ/mol）很小,从而揭示了苯并氧化呋咱比邻二亚硝基苯更稳定·此外,进行了HF/3-21G、HF/6-31G(d)和MP2/6-31G(d)//6-31G(d)水平下相应的计算,发现B3LYP-DFT的结果较abinitio为优。谐振动频率的B3LYP/6-31G(d)计算还支持了邻二亚硝基苯为苯并氧化呋咱“自-自”互变重排反应的中间体。