DFT and CASPT2 study of two thermal reactions of nitromethane: C–N bond cleavage and nitro-to-nitrite isomerization. An example of the inverse symmetry breaking deficiency in density functional calculations of an homolytic dissociation

The reaction paths of nitromethane leading to the dissociation products or isomerization to methyl nitrite have been computationally investigated at the CAS-SCF and DFT levels of theory. Additionally, the CAS-SCF wave functions were used as reference in a second-order perturbation treatment, CASPT2, in order to obtain a good estimate for the activation energy of each reaction path. Both methods predict the isomerization as a concerted reaction. However, the behavior of the two approximations with respect to dissociation is rather different; while CASPT2 predicts a barrier height of (≈59 kcal/mol) in good accordance with the experimental activation energy (59.0 kcal/mol), B3-LYP/6-31G^* calculations overestimate the barrier for more than 30 kcal/mol. The DFT prediction of the dissociation channel exhibits inverse symmetry breaking, dissociating to the unphysical absurd CH₃^δ+ plus NO₂^δ−.     

Theoretical study on the reaction pathways of HFCO + H2O

For the reaction of methanoyl fluoride with water, both optimized structures and vibrational wavenumbers of reaction intermediates, transition structures and product complexes were calculated and characterized with theory at the MP2/6-311++G(d,p) level. Including a catalytic path and concerted and stepwise hydrolysis paths, possible reaction mechanisms were also investigated. The catalytic reaction of HFCO yielding HF and CO has the smallest activation barrier, 29.6 kcal/mol, whereas for the concerted hydrolysis 33.0 kcal/mol is required to overcome the barrier to form transoid HCOOH + HF, which is less than for the stepwise counterpart, 42.0 kcal/mol.     

Mechanisms of DABCO‐ and DMAP‐catalyzed [2 + 4] cycloaddition reactions of methylallenoate with methyleneindolonone: A DFT study

The mechanisms and stereoselectivities of the [2 + 4] cycloaddition reaction of methylallenoate R1 with methyleneindolonone R2 catalyzed by DABCO (Equation 1) and DMAP (Equation 2) organocatalysts have been examined with density functional theory (M06‐2X) calculations. Several possible reaction pathways (paths 1a, 1b, and 1c for Equation 1 and paths 2a and 2b for Equation 2) were located and compared. The results of our study reveal that for both reactions, three reaction stages have been characterized: nucleophilic addition of the catalyst ( DABCO or DMAP ) to R1 (Stage I ), addition of the other reactant R2 (Stage II ), intramolecular cycloaddition and liberation of the catalyst ( DABCO or DMAP ) afforded the final product (Stage III ). For the DABCO ‐catalyzed cycloaddition, we predict that path 1a leading to P(E) is the most energy favorable pathway among the three possible pathways. The carbon–carbon bond formation step is the rate‐determining step (ΔG ^?=23.6 kcal/mol). With DMAP catalyst, the same reaction gave P(Z) as the major product. The barrier for the rate‐determining step (addition of R1 to DMAP ) is 25.8 kcal/mol. The calculated results are in agreement with the experimental findings. Moreover, for both reactions, the analysis of global reactivity indexes has been carried out to examine the role of catalyst. The present study should provide a general mechanistic framework for the rational design of this kind of reactions.     

Identification of the transition states in the inversion of 1,4-benzodiazepines with the ab initio replica path method

The inversion of four 1,4-benzodiazepines was investigated with the ab initio "replica path method" with density functional theory at the B3LYP/6-31G* level. The reaction path constructed with this method for the inversion provides an approximate transition state (TS) geometry, which, upon further TS optimization, leads to the TS geometry characterized by a single vibrational frequency. 1,4-Benzodiazepines lacking a 5-phenyl ring have a single reaction path for the inversion with Cs symmetry at the TS. In contrast, the inversion of benzodiazepines with a 5-phenyl ring, such as the peripheral benzodiazepine receptor ligand 4'-chlorodiazepam (Ro5-4864) and its N1-desmethyl analog (Ro5-2752), can proceed through multiple reaction paths having a TS with or without Cs symmetry. Notably, the replica path method found a path connecting two asymmetric TSs of 4'-chlorodiazepam via a symmetrical TS, suggesting that these inversion paths can be readily crossed from one to another. The stabilization energies gained by 4'-chlorodiazepam and its N1-desmethyl analog from the breaking of Cs symmetry at the TS were calculated to be 0.10 and 0.07 kcal/mol, respectively. The origin of the broken symmetry of Cs was traced to the coupling of the puckering of the diazepine ring with the rotation of the chlorophenyl ring. These results show the advantages of the replica path method for locating the TSs as well as for constructing the reaction paths for the inversion of 1,4-benzodiazepines.     

二氟亚烷基卡宾与甲醛环加成反应机理的理论研究

用二阶微扰理论研究单重态二氟亚烷基卡宾与甲醛发生的环加成反应机理,采用MP2/6-31G^*方法计算了势能面上各驻点的构型参数、振动频率和能量.结果表明,单重态二氟亚烷基卡宾与甲醛的环加成反应主要有两种反应通道,通道1中,两个反应物经a,b和c三条反应途径生成三元环构型的产物P1,其中途径c是主反应途径,该途径有两步组成:(Ⅰ)二氟亚烷基卡宾与甲醛生成了1个富能中间体(INT1c),是无势垒放热反应,放出能量为219.18kJ/mol;(Ⅱ)中间体(INT1c)异构化为产物二氟亚烷基环氧乙烷,其势垒为134.71kJ/mol.通道2的反应途径由三步组成:(Ⅰ)反应物首先生成了1个富能中间体(INT1b),为无势垒的放热反应,放出的能量142.77kJ/mol;(Ⅱ)中间体(INT1b)异构化成另一中间体(INT2),其势垒为22.31kJ/mol;(Ⅲ)中间体(INT2)异构化成四元环构型产物P2,其势垒为11.98kJ/mol.     

A DFT study of SiH(4) activation by Cp(2)LnH

A theoretical study of SiH(4) activation by Cp(2)LnH complexes for the entire series of lanthanides has been carried out at the DFT-B3PW91 level of theory. The reaction paths corresponding to H/H exchange and silylation, formation of Cp(2)Ln(SiH(3)), have been computed. They both occur via a single-step sigma-bond metathesis mechanism. For the athermal H/H exchange reaction, the calculated activation barrier averages 1.8 kcal.mol(-)(1) relative to the precursor adduct Cp(2)LnH(eta(2)-SiH(4)) for all lanthanide elements. The silylation path is slightly exogenic (DeltaE approximately -6.5 kcal.mol(-1)) with an activation barrier averaging 5.2 kcal.mol(-1) relative to the precursor adduct where SiH(4) is bonded by two Si-H bonds. Both pathways are therefore thermally accessible. The H/H exchange path is calculated to be kinetically more favorable whereas the silylation reaction is thermodynamically preferred. The reactivity of this familly of lanthanide complexes with SiH(4) contrasts strongly with that obtained previously with CH(4). The considerably lower activation barrier for silylation relative to methylation is attributed to the ability of Si to become hypervalent.     

F+CH_3OH碰撞反应机机理和反应势能面

以MP2（full）/6-311 + +g(d,p)水平上详细研究了氟原子与甲醇抽氢反应的 多通道反应机理，得到了各条通道中涉及的驻点的构型和振动频率及其能量，给出 了两张完整的反应势能面，结果表明，氟原子从C原子上抽氢时有一条明显的最低 能量通道，而从氧原子上抽氢时要涉及多条分支通道和多个驻点构型，给出了各分 支通道势能面示意图，结果表明以形成五元环状过渡态通道为优势通道，计算得到 经途径1生成CH_2OH时反应放热170.62kJ/mol，经分支途径6生成CH_3O自由基时反 应放热119.4 kJ/mol，此结果与实验值一致。     

Explanation of the unusual temperature dependence of the atmospherically important OH + H(2)S --> H(2)O + HS reaction and prediction of the rate constant at combustion temperatures

Rate constants for the OH + H2S --> H2O + HS reaction, which is important for both atmospheric chemistry and combustion, are calculated by direct dynamics with the M06-2X density functional using the MG3S basis set. Energetics are compared to high-level MCG3/3//MC-QCISD/3 wave function theory and to results obtained by other density functionals. We employ canonical variational transition-state theory with multidimensional tunneling contributions and scaled generalized normal-mode frequencies evaluated in redundant curvilinear coordinates with anharmonicity included in the torsion. The transition state has a quantum mechanically distinguishable, nonsuperimposable mirror image that corresponds to a separate classical reaction path; the effect of the multiple paths is examined through use of a symmetry number and by torsional methods. Calculations with the reference-potential Pitzer-Gwinn treatment of the torsional mode agree with experiment, within experimental scatter, and predict a striking temperature dependence of the activation energy, increasing from -0.1 kcal/mol at 200 K to 0.2, 1.0, 3.4, and 9.8 kcal/mol at 300, 500, 1000, and 2400 K. The unusual temperature dependence arises from a dynamical bottleneck at an energy below reactants, following an addition complex on the reaction path with a classical binding energy of 4.4 kcal/mol. As a way to check the mechanism, kinetic isotope effects of the OH + D2S and OD + D2S reactions have been predicted.     

Reaction mechanism of platinum dimer cation with ammonia based on the relativistic density functional study

The gas‐phase reactions between Pt and NH₃ have been investigated using the relativistic density functional approach (ZORA‐PW91/TZ2P). The quartet and doublet potential energy surfaces of Pt + NH₃ have been explored. The minimum energy reaction path proceeds through the following steps: Pt(⁴Σ_u) + NH₃ → q‐1 → d‐2 → d‐3 → d‐4 → d‐Pt₂NH⁺ + H₂. In the whole reaction pathway, the step of d‐2 → d‐3 is the rate‐determining step with a energy barrier of 36.1 kcal/mol, and exoergicity of the whole reaction is 12.0 kcal/mol. When Pt₂NH⁺ reacts with NH₃ again, there are two rival reaction paths in the doublet state. One is degradation of NH and another is loss of H₂. In the case of degradation of NH, the activation energy is only 3.4 kcal/mol, and the overall reaction is exothermic by 8.9 kcal/mol. Thus, this reaction is favored both thermodynamically and kinetically. However, in the case of loss of H₂, the rate‐determining step's energy barrier is 64.3 kcal/mol and the overall reaction is endothermic by 8.5 kcal/mol, so it is difficult to take place. Predicted relative energies and barriers along the suggested reaction paths are in reasonable agreement with experimental observations. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A computational investigation on the photochemistry of the popular spin-trap agent N-tert-butyl-α-phenylnitrone (PBN) and thermal isomerization pathways of its photoproduct oxaziridine

Computational investigation on the low-lying photo-excited states of N-tert-butyl-α-phenylnitrone (PBN), a well-known spin-trap agent, has revealed its photo-product (oxaziridine) formation channel. The S₀-S₂ vertical excitation in PBN is subsequently followed by a non-radiative decay pathway through S₂/S₁ and S₀/S₁ conical intersections (CIs) with CNO-kinked structures, situated around 23 kcal/mol and 45 kcal/mol below the vertically excited S₂ state, respectively. The reverse photo-process of PBN formation involves photo-excitation of oxaziridine to its S₂ and S₃ photo-excited states. The forward photo-isomerization leads to the trans-oxaziridine with a backside CNO kink (trans-OX_B) while the reverse path studied by us, connects its front-side CNO-kinked analogue (trans-OX_F) with the PBN. Our search for the reverse thermal reaction paths from these two oxazirdines has led to their corresponding transition states, one at 35 kcal/mol and the other at 27 kcal/mol above trans-OX_F and trans-OX_B geometries, respectively. They lead to two different isomers (E and Z) of PBN which supports the reported nature of products from the trans-oxaziridine in this thermal reaction. The inversion path of the chiral nitrogen atom of this N-tert-butyl-oxaziridine (barrier 21 kcal/mol) has also been tracked. This reaction path has been compared with that of the N-methyl (barrier 30 kcal/mol) and N-acyl (barrier 10.5 kcal/mol) oxaziridine analogues.     

On the dimerization process of nitroso compounds

Summary Many organic C-nitroso compounds R-NO form stable dimers with a covalent NN bond. To gain insight into the dimerization reaction 2 R-NO (R-NO)₂ a theoretical study of the dimerization to atrans-form was performed using HNO as a model compound. Complete geometry optimizations were carried out at the HF, MP2 and QCISD levels using a 6–31G* basis. In the stationary points energies were calculated at the MP4(SDTQ) and QCISD(T) levels. For the equilibrium structure of the monomer and dimers stable RHF solutions were found, whereas for the TS UHF and UMPn calculations were applied. Extensive spin contamination was found in the UHF wavefunction, and projections up tos+4 were invoked. Relative energies were corrected for differences in ZPE. Calculations were made (a) for the least-motion path (C _2h symmetry) and (b) for a path with complete relaxation of all internal coordinates. Along the latter path a TS having virtuallyC _i symmetry was found. Along path (a) an activation energy of around 150 kcal/mol was predicted, in conformity with a symmetry forbidden reaction. On the relaxed path (b) the barrier to dimerization was estimated to be 10.7 kcal/mol at the MP4(SDTQ)//MP2 level, and 10.9 kcal/mol at the QCISD(T)//QCISD level. Unscaled ZPE corrections, calculated at the SCF level, changed these values to 12.7 and 12.9 kcal/mol, respectively. The reaction energy for the dimerization process is predicted to be – 17.2 kcal/mol at the MP4(SDTQ)//MP2 level corrected for ZPE. Calculations at the G1 level gave a corresponding value of – 16.4 kcal/mol. The equilibrium constant for the association to thetrans dimer is estimated to beK _p =259 atm, indicating that the dimer should be an observable species in the gas phase.     

Ab initio studies of isomerization and dissociation reactions of methyl peroxynitrate

Using the CCSD(T)/cc-pVDZ//B3LYP/6-311G(2d,2p) method, we calculated the detailed potential energy surfaces (PESs) for the unimolecular isomerization and decomposition of methyl peroxynitrate (CH₃O₂NO₂). The results show that there are the two most stable isomers, IS1a and IS1b, which are a pair of mirror image isomers. From IS1a and IS1b, different isomerization and unimolecular decomposition reaction channels have been studied and discussed. Among them, the predominant thermal decomposition pathways are those leading to CH₃O₂ + NO₂ and cis-CH₃ONO + O₂. The former is the lowest-energy path through the direct O–N bond rupture in IS1a or IS1b. The PES along the O–N bond in IS1a has been scanned, where the energy of IS1a reaches maximum value of 23.5 kcal/mol when the O–N bond is stretched to about 2.8 Å. This energy is 2.7 kcal/mol larger than the O–N bond dissociation energy (BDE) and 2.8 kcal/mol larger than the experimental active energy. In addition, because the energy barriers of IS1a isomerization to IS2a are 23.8 kcal/mol, close to the 20.8 kcal/mol O–N BDE in IS1a or IS1b, the isomerization reaction may compete with the direct bond rupture dissociation reaction.     

Ab Initio MO Studies on the Reaction Mechanism for Carbonyl Insertion Catalyzed by Carbonyl Cobalt Complex

IntroductionSome of ah l'nito MO studies on elementary reactions catalyzed by organometallic compounds have been published[' 41, which are of momentous significance for better understanding the catalytic action and developing the catalytic agent. At the same time, there still aremany unknown aspects that remain to be studied. Carbonyl insertion is one of the key reac.tlons, and it can go in three possible pathsL'--'j, i. e., carbonyl migration, ethyl migrationand concerted reaction paths. Th…     

HNCO+OH——NH~2+CO~2反应理论研究

用从头算UHF/6-31G基组研究了异氰酸和羟基生成氨基和二氧化碳即HNCO+OH--NH~2+CO~2的反应机理.优化得到了反应途径上的过渡态和中间体,并通过振动分析对过渡态和中间体进行了确认.在UMP4/6-31G水平上计算了它们的能量,同时对零点能进行了较正.计算结果表明:此反应是多步反应,先后通过3个过渡态(TS1,TS2,TS3),2个内旋转位垒(TSI,TSII),4个中间体(IM1,IM2,IM3,IM4),其中,IM3--TS2这一步为整个反应的决速步骤,速控步的活化能为202.388kJ/mol.与异氰酸和羟基作用的另一反应通道(即HNCO+OH--H~2O+NCO)的活化能(69.038kJ/mol)比较,可看出所研究反应通道为次要反应通道,这与实验结果是一致的。     

Theoretical study of the mechanism and rate constant of the B + CO2 reaction

The different stationary points on the potential energy surface relative to the title reaction have been reinvestigated at the B3LYP/aug-cc-pVDZ level with relative energies computed at the CCSD(T)/aug-cc-pVTZ level with B3LYP/aug-cc-pVDZ optimized geometries and by using the G3B3 composite method. Two entrance channels have been identified. The first one corresponds to boron addition at one of the oxygen atoms of the CO 2 molecule leading to trans-BOCO, which is found to be about 27 kcal/mol exothermic with a potential energy barrier of 16.4 kcal/mol (G3B3). The second channel, which has not been identified in previous theoretical works, corresponds to a direct insertion of the boron atom into a CO bond and leads to OBCO. The B + CO 2 --> OBCO step is found to be about 84 kcal/mol exothermic and needs to overcome a potential energy barrier of only 3.6 kcal/mol (G3B3). The rate constant at 300 K of the insertion step, calculated by using TST theory with G3B3 calculated activation energy value, is 5.4 10 (-14) cm (3) molecule (-1) s (-1), in very good agreement with the experimental data ((7.0 +/- 2.8) 10 (-14) cm (3) molecule (-1) s (-1), DiGiuseppe, T. G.; Davidovits, P. J. Chem. Phys. 1981, 74, 3287). The one corresponding to the addition process is found to be several orders of magnitude smaller because of a much higher potential energy barrier. The addition channel would not contribute to the title reaction even at high temperature. A modified Arrhenius equation has been fitted in the 300-1000 K temperature range, which might be useful for chemical models.     

Ab initio rate constants from hyperspherical quantum scattering: application to H+C2H6 and H+CH3OH

The dynamics and kinetics of the abstraction reactions of H atoms with ethane and methanol have been studied using a quantum mechanical procedure. Bonds being broken and formed are treated with explicit hyperspherical quantum dynamics. The ab initio potential energy surfaces for these reactions have been developed from a minimal number of grid points (average of 48 points) and are given by analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimized using the second order perturbation theory method with a correlation consistent polarized valence triple zeta basis set. Single point energies are calculated on the optimized geometries with the coupled cluster theory and the same basis set. The reaction of H with C2H6 is endothermic by 1.5 kcal/mol and has a vibrationally adiabatic barrier of 12 kcal/mol. The reaction of H with CH3OH presents two reactive channels: the methoxy and the hydroxymethyl channels. The former is endothermic by 0.24 kcal/mol and has a vibrationally adiabatic barrier of 13.29 kcal/mol, the latter reaction is exothermic by 7.87 kcal/mol and has a vibrationally adiabatic barrier of 8.56 kcal/mol. We report state-to-state and state-selected cross sections together with state-to-state rate constants for the title reactions. Thermal rate constants for these reactions exhibit large quantum tunneling effects when compared to conventional transition state theory results. For H+CH3OH, it is found that the CH2OH product is the dominant channel, and that the CH3O channel contributes just 2% at 500 K. For both reactions, rate constants are in good agreement with some measurements.     

Theoretical kinetic study of two competing paths in the H2O2 + O(P) → HO2 + OH reaction

Ab initio-TST calculations were carried out to study the kinetics of the title reaction. The H atom and the OH abstraction paths leading to the same products HO₂ and OH have been considered. The ZPE and BSSE corrected classical barrier heights were predicted to be 7.4 and 17.3 kcal/mol, respectively. Calculated thermal rate constants over the temperature range 300–5000 K showed that the H-abstraction path was the most likely to occur for temperatures below 2500 K which confirms the result found in a previous study [Y. Tarchouna, M. Bahri, N. Jaïdane, Z. Ben Lakdar, J. Mol. Struct. (Theochem), 189 (2003) 664]. The contribution of OH abstraction path to the reaction was predicted to be important for high temperatures.     

Theoretical predictions of the initial decomposition steps of dimethylnitramine

The structures and energies of the reactants, products, and transition states of the initial steps in the gas-phase decomposition of dimethylnitramine (DMNA) have been determined by quantum chemical calculations at the B3LYP density-functional theory, MP2, and G2 levels. The pathways considered are NO2 elimination, HONO elimination, and nitro-nitrite rearrangement. The NO2 elimination is predicted to be the main channel of the gas-phase decomposition of DMNA in accord with experiment. The values of the Arrhenius parameters, log A=16.6+/-0.5 and Ea=40.0+/-0.6 kcal/mol, for the N-NO2 bond-fission reaction were obtained using a canonical variational theory with B3LYP energies and frequencies. The HONO-elimination channel has the next lowest activation energy of 44.7+/-0.5 kcal/mol (log A=13.6+/-0.5) and is characterized by a five-member transition-state configuration in which a hydrogen atom from one of the methyl groups is transferred to an oxygen atom of NO2. Tunneling contributions to the rate of this reaction have been estimated. The nitro-nitrite rearrangement reaction occurs via a transition state in which both oxygen atoms of NO2 are loosely bound to the central nitrogen atom, for which Rice-Ramsperger-Kassel-Marcus theory predicts log A=14.4+/-0.6 and Ea=54.1+/-0.8 kcal/mol.     

Density functional calculations of ATP systems. 2. ATP hydrolysis at the active site of actin

The hydrolysis of adenosine 5'-triphosphate (ATP) at the active site of actin has been studied using density functional calculations. The active site is modeled by the triphosphate tail of ATP, an Mg cation, surrounding water molecules, and the nearby protein residues. Four reaction paths have been followed by constraining coordinates that allow phosphate stretching, nucleophilic attack of the catalytic water, and OH(-) formation via water deprotonation. The lowest-energy barrier (21.0 kcal/mol) is obtained for a dissociative reaction where the terminal phosphate breaks on approaching the catalytic water, followed by proton release via a proton wire mechanism. A higher barrier (39.6 kcal/mol) results for an associative reaction path where OH(-) is formed first, with a pentacoordinated phosphorus atom (P-O distances 2.1 A). Stretching the terminal bridging P-O bond results in bond rupture at 2.8 A with an energy barrier of 28.8 kcal/mol. The residues Gln137 and His161 are not important in the reactions, but insight into their roles in vivo has been obtained. The favored coordination of the end products H(2)PO(4)(-) and ADP(3-) includes a hydrogen bond and an O-Mg-O bridge between the phosphates as well as a hydrogen bond between H(2)PO(4)(-) and the Ser14 side chain. The total energy is 2.1 kcal/mol lower than in the initial reactants. Classical simulations of ATP- and ADP.P(i)-actin show few hydrolysis-induced differences in the protein structure, indicating that phosphate migration is necessary for a change in conformation.