Heck‐Type Reactions of Imine Derivatives: A DFT Study

The mechanism of a recently discovered intramolecular Heck‐type coupling of oximes with aryl halides (Angew. Chem. Int. Ed. 2007 , 46, 6325) was systematically studied by using density functional methods enhanced with a polarized continuum solvation model. The overall catalytic cycle of the reaction was found to consist of four steps: oxidative addition, migratory insertion, β‐H elimination, and catalyst regeneration, whereas an alternative base‐promoted C? H activation pathway was determined to be less favorable. Migratory insertion was found to be the rate determining step in the catalytic cycle. The apparent activation barrier of migratory insertion of the (E)‐oxime was +20.5 kcal mol^?1, whereas the barrier of (Z)‐oxime was as high as +32.7 kcal mol^?1. However, (Z)‐oxime could isomerize to form the more active (E)‐oxime with the assistance of K₂CO₃, so that both the (E)‐ and (Z)‐oxime substrates could be transformed to the desired product. Our calculations also indicated that the Z product was predominant in the equilibrium of the isomerization of the imine double bond, which constituted the reason for the good Z‐selectivity observed for the reaction. Furthermore, we examined the difference between the intermolecular Heck‐type reactions of imines and of olefins. It was found that in the intermolecular Heck‐type coupling of imines, the apparent activation barrier of migratory insertion was as high as +35 kcal mol^?1, which should be the main obstacle of the reaction. The analysis also revealed the main problem for the intermolecular Heck‐type reactions of imines, which was that the breaking of a C?N π bond was much more difficult than the breaking of a C?C π bond. After systematic examination of a series of substituted imines, (Z)‐N‐amino imine and N‐acetyl imine were found to have relatively low barriers of migratory insertion, so that they might be possible substrates for intermolecular Heck‐type coupling.     

Equilibrium Between Acetaldehyde, Acetyl Halide and Corresponding 1-Haloethyl Acetates in Gas Chromatograph Injector. Calculation of Enthalpies of Reaction

Gas chromatographic analysis of 1-chloroethyl acetate and 1-bromoethyl acetate revealed that the equilibrium between acetaldehyde, acetyl halide and the corresponding 1-haloethyl acetates exists in the injector of the chromatograph. Analyses were performed under strictly isothermal conditions of both injector and column at different temperatures. The results allowed calculation of the enthalpy of the reaction: CH₃CHO + CH₃COX → CH₃COOCHXCH₃ Enthalpies calculated for the reaction are as follows: ?17.3 [kcal mol^?1] (X=Cl) and ?18.5 [kcal mol^?1] (X=Br).     

Theoretical Investigations on the Mechanism of Hetero‐Diels–Alder Reactions of Brassard’s Diene and 1, 3‐Butadiene Catalyzed by a Tridentate Schiff Base Titanium(IV) Complex

The mechanism of the hetero‐Diels–Alder reactions of Brassard’s diene and 1,3‐butadiene catalyzed by a titanium(IV) complex of a tridentate Schiff base was investigated by DFT and ONIOM methods. The calculations indicate that the mechanism of the reaction is closely related to the nucleophilicity–electrophilicity between diene and carbonyl substrates. A stepwise pathway is adopted for Brassard’s diene, and the step corresponding to the formation of the C? C bond is predicted to be the rate‐determining step with a free‐energy barrier of 8.4 kcal mol^?1. For 1,3‐butadiene, the reaction takes place along a one‐step, two‐stage pathway with a free‐energy barrier of 14.9 kcal mol^?1. For Brassard’s diene as substrate, the OCH₃ and OSi(CH₃)₃ substituents may play a key role in the formation of the transition state and zwitterionic intermediate by participating in charge transfer from Brassard’s diene to formaldehyde. The combination of the phenyl groups at the amino alcohol moiety and the ortho‐tert‐butyl group of the salicylaldehyde moiety in the chiral tridentate Schiff base ligand plays an important role in the control of the stereoselectivity, which is in agreement with experimental observations.     

σ‐Insertive Mechanism versus Concerted Non‐insertive Mechanism in the Intramolecular Hydroamination of Aminoalkenes Catalyzed by Phenoxyamine Magnesium Complexes: A Synthetic and Computational Study

The phenoxyamine magnesium complexes [{ONN}MgCH₂Ph] ( 4 a : {ONN}=2,4‐tBu₂‐6‐(CH₂NMeCH₂CH₂NMe₂)C₆H₂O^?; 4 b : {ONN}=4‐tBu‐2‐(CH₂NMeCH₂CH₂NMe₂)‐6‐(SiPh₃)C₆H₂O^?) have been prepared and investigated with respect to their catalytic activity in the intramolecular hydroamination of aminoalkenes. The sterically more shielded triphenylsilyl‐substituted complex 4 b exhibits better thermal stability and higher catalytic activity. Kinetic investigations using complex 4 b in the cyclisation of 1‐allylcyclohexyl)methylamine ( 5 b ), respectively, 2,2‐dimethylpent‐4‐en‐1‐amine ( 5 c ), reveal a first‐order rate dependence on substrate and catalyst concentration. A significant primary kinetic isotope effect of 3.9±0.2 in the cyclisation of 5 b suggests significant N?H bond disruption in the rate‐determining transition state. The stoichiometric reaction of 4 b with 5 c revealed that at least two substrate molecules are required per magnesium centre to facilitate cyclisation. The reaction mechanism was further scrutinized computationally by examination of two rivalling mechanistic pathways. One scenario involves a coordinated amine molecule assisting in a concerted non‐insertive N?C ring closure with concurrent amino proton transfer from the amine onto the olefin, effectively combining the insertion and protonolysis step to a single step. The alternative mechanistic scenario involves a reversible olefin insertion step followed by rate‐determining protonolysis. DFT reveals that a proton‐assisted concerted N?C/C?H bond‐forming pathway is energetically prohibitive in comparison to the kinetically less demanding σ‐insertive pathway (ΔΔG^≠=5.6 kcal mol^?1). Thus, the σ‐insertive pathway is likely traversed exclusively. The DFT predicted total barrier of 23.1 kcal mol^?1 (relative to the {ONN}Mg pyrrolide catalyst resting state) for magnesium?alkyl bond aminolysis matches the experimentally determined Eyring parameter (ΔG^≠=24.1(±0.6) kcal mol^?1 (298 K)) gratifyingly well.     

Fragmentation of oxime and silyl oxime ether odd‐electron positive ions by the McLafferty rearrangement: new insights on structural factors that promote α,β fragmentation

The McLafferty rearrangement is an extensively studied fragmentation reaction for the odd‐electron positive ions from a diverse range of functional groups and molecules. Here, we present experimental and theoretical results of 12 model compounds that were synthesized and investigated by GC‐TOF MS and density functional theory calculations. These compounds consisted of three main groups: carbonyls, oximes and silyl oxime ethers. In all electron ionization mass spectra, the fragment ions that could be attributed to the occurrence of a McLafferty rearrangement were observed. For t‐butyldimethylsilyl oxime ethers with oxygen in a β‐position, the McLafferty rearrangement was accompanied by loss of the t‐butyl radical. The various mass spectra showed that the McLafferty rearrangement is relatively enhanced compared with other primary fragmentation reactions by the following factors: oxime versus carbonyl, oxygen versus methylene at the β‐position and ketone versus aldehyde. Calculations predict that the stepwise mechanism is favored over the concerted mechanism for all but one compound. For carbonyl compounds, C–C bond breaking was the rate‐determining step. However, for both the oximes and t‐butyldimethylsilyl oxime ethers with oxygen at the β‐position, the hydrogen transfer step was rate limiting, whereas with a CH₂ group at the β‐position, the C–C bond breaking was again rate determining. n‐Propoxy‐acetaldehyde, bearing an oxygen atom at the β‐position, is the only case that was predicted to proceed through a concerted mechanism. The synthesized oximes exist as both the (E)‐ and (Z)‐isomers, and these were separable by GC. In the mass spectra of the two isomers, fragment ions that were generated by the McLafferty rearrangement were observed. Finally, fragment ions corresponding to the McLafferty reverse charge rearrangement were observed for all compounds at varying relative ion intensities compared with the conventional McLafferty rearrangement. Copyright © 2012 John Wiley & Sons, Ltd.     

Evidence for a Rapid Degenerate Hetero‐Cope‐Type Rearrangement in [Cp*W(S)2S‐CH2‐CHCH2]

Treatment of the salt [PPh₄]⁺[Cp*W(S)₃]^? ( 6 ) with allyl bromide gave the neutral complex [Cp*W(S)₂S‐CH₂‐CH?CH₂] ( 7 ). The product 7 was characterized by an X‐ray crystal structure analysis. Complex 7 features dynamic NMR spectra that indicate a rapid allyl automerization process. From the analysis of the temperature‐dependent NMR spectra a Gibbs activation energy of ΔG^≠ (278 K)≈13.7±0.1 kcal mol^?1 was obtained [ΔH^≠≈10.4±0.1 kcal mol^?1; ΔS^≠≈?11.4 cal mol^?1 K^?1]. The DFT calculation identified an energetically unfavorable four‐membered transition state of the “forbidden” reaction and a favorable six‐membered transition state of the “Cope‐type” allyl rearrangement process at this transition‐metal complex core.     

Oxidative Dehydrogenation of Propane over a VO2‐Exchanged MCM‐22 Zeolite: A DFT Study

The adsorption and the mechanism of the oxidative dehydrogenation (ODH) of propane over VO₂‐exchanged MCM‐22 are investigated by DFT calculations using the M06‐L functional, which takes into account dispersion contributions to the energy. The adsorption energies of propane are in good agreement with those from computationally much more demanding MP2 calculations and with experimental results. In contrast, B3LYP binding energies are too small. The reaction begins with the movement of a methylene hydrogen atom to the oxygen atom of the VO₂ group, which leads to an isopropyl radical bound to a HO? V? O intermediate. This step is rate determining with the apparent activation energy of 30.9 kcal mol^?1, a value within the range of experimental results for ODH over other silica supports. In the propene formation step, the hydroxyl group is the more reactive group requiring an apparent activation energy of 27.7 kcal mol^?1 compared to that of the oxy group of 40.8 kcal mol^?1. To take the effect of the extended framework into account, single‐point calculations on 120T structures at the same level of theory are performed. The apparent activation energy is reduced to 28.5 kcal mol^?1 by a stabilizing effect caused by the framework. Reoxidation of the catalyst is found to be important for the product release at the end of the reaction.     

The fragmentations of [M-H]- anions derived from underivatised peptides. The side-chain loss of H2S from Cys. A joint experimental and theoretical study

Loss of H₂S is the characteristic Cys side‐chain fragmentation of the [M? H]^? anions of Cys‐containing peptides. A combination of experiment and theory suggests that this reaction is initiated from the Cys enolate anion as follows: RNH‐^?C(CH₂SH)CONHR′ Ø [RNHC(?CH₂)CONHR′ (HS^?)] Ø [RNHC(?CH₂)CO‐HNR′‐H]^?+H₂S. This process is facile. Calculations at the HF/6‐31G(d)//AM1 level of theory indicate that the initial anion needs only ≥20.1 kcal mol^?1 of excess energy to effect loss of H₂S. Loss of CH₂S is a minor process, RNHCH(CH₂SH)CON^?‐R′ Ø RNHCH(CH₂S^?)CONHR′ Ø RNH ^?CHCONHR+CH₂S, requiring an excess energy of ≥50.2 kcal mol^?1. When Cys occupies the C‐terminal end of a peptide, the major fragmentation from the [M–H]^? species involves loss of (H₂S+CO₂). A deuterium‐labelling study suggests that this could either be a charge‐remote reaction (a process which occurs remote from and uninfluenced by the charged centre in the molecule), or an anionic reaction initiated from the C‐terminal CO₂^? group. These processes have barriers requiring the starting material to have an excess energy of ≥79.6 (charge‐remote) or ≥67.1 (anion‐directed) kcal mol^?1, respectively, at the HF/6‐31G(d)//AM1 level of theory. The corresponding losses of CH₂O and H₂O from the [M? H]^? anions of Ser‐containing peptides require ≥35.6 and ≥44.4 kcal mol^?1 of excess energy (calculated at the AM1 level of theory), explaining why loss of CH₂O is the characteristic side‐chain loss of Ser in the negative ion mode. Copyright © 2003 John Wiley & Sons, Ltd.     

C-H bond activation of methane in aqueous solution: a hybrid quantum mechanical/effective fragment potential study

In this study, we investigated the C? H bond activation of methane catalyzed by the complex [PtCl₄]^2?, using the hybrid quantum mechanical/effective fragment potential (EFP) approach. We analyzed the structures, energetic properties, and reaction mechanism involved in the elementary steps that compose the catalytic cycle of the Shilov reaction. Our B3LYP/SBKJC/cc‐pVDZ/EFP results show that the methane activation may proceed through two pathways: (i) electrophilic addition or (ii) direct oxidative addition of the C? H bond of the alkane. The electrophilic addition pathway proceeds in two steps with formation of a σ‐methane complex, with a Gibbs free energy barrier of 24.6 kcal mol^?1, followed by the cleavage of the C? H bond, with an energy barrier of 4.3 kcal mol^?1. The activation Gibbs free energy, calculated for the methane uptake step was 24.6 kcal mol^?1, which is in good agreement with experimental value of 23.1 kcal mol^?1 obtained for a related system. The results shows that the activation of the C? H bond promoted by the [PtCl₄]^2? catalyst in aqueous solution occurs through a direct oxidative addition of the C? H bond, in a single step, with an activation free energy of 25.2 kcal mol^?1, as the electrophilic addition pathway leads to the formation of a σ‐methane intermediate that rapidly undergoes decomposition. The inclusion of long‐range solvent effects with polarizable continuum model does not change the activation energies computed at the B3LYP/SBKJC/cc‐pVDZ/EFP level of theory significantly, indicating that the large EFP water cluster used, obtained from Monte Carlo simulations and analysis of the center‐of‐mass radial pair distribution function, captures the most important solvent effects. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Thermal decomposition of formaldehyde diperoxide in aqueous solution

Thermal decomposition of formaldehyde diperoxide (1,2,4,5-tetraoxane) in aqueous solution with an initial concentration of 6.22 × 10^?3 M was studied in the temperatures range from 403 to 439 K. The reaction was found to follow first-order kinetic law, and formaldehyde was the major decomposition product. The activation parameters of the initial step of the reaction (ΔH ^≠ = 15.25 ± 0.5 kcal mol^?1, ΔS ^≠ = ?47.78 ± 0.4 cal mol^?1K^?1, E _a = 16.09 ± 0.5 kcal mol^?1) support a mechanism involving homolytic rupture of one peroxide bond in the 1,2,4,5-tetraoxane molecule with participation of the solvent and formation of a diradical intermediate.     

Aktivierungsenergie und Mechanismus der intramolekularen Wasserstoff-Wanderung bei Radikal-Kationen von Keto/Enol-Tautomeren in der Gasphase

The collisional activation mass spectra prove that non-decomposing ionized methyl acetate [CH₃COOCH₃]⁺? and its enolic isomer [CH₂?C(OH)OCH₃]⁺? exist as stable species in potential wells. It is shown, however, that prior to CH₃O? loss the decomposing [CH₂?C(OH)OCH₃]⁺? ion isomerizes via a rate determining symmetry forbidden [1.3] hydrogen rearrangement to ionized methyl acetate. The alternative mode of two consecutive formally symmetry allowed [1.2] hydrogen migrations can be certainly excluded for this isomerization. The activation energy of such hydrogen rearrangements is of the order of 41–83 kcal · mol^?1 depending on the electronic nature of the cations (“open” or “closed” shell systems).     

Reduction Process of Tetraplatin in the Presence of Deoxyguanosine Monophosphate (dGMP): A Computational DFT Study

The reduction mechanism of [Pt^IV(dach)Cl₄] (dach=diaminocyclohexyl) in the presence of dGMP was studied. The first step is substitution of a chloro ligand by dGMP, followed by nucleophilic attack of a phosphate or sugar oxygen atom to the C8‐position of guanine. Subsequent reduction forms the [Pt^II(dach)Cl₂] complex. The whole process is completed by a hydrolysis. Two different pathways for the substitution reaction were examined: a direct associative and a Basolo–Pearson autocatalytic mechanism. All the explored structures were optimized at the B3LYP‐D3/6‐31G(d) level and by using the COSMO solvation model with Klamt's radii. Single‐point energetics was determined at the B3LYP‐GD3BJ/6‐311++G(2df,2pd)/PCM/scaled‐UAKS level. Activation barriers were used for an estimation of the rate constants and these were compared with experimental values. It was found that the rate‐determining step is the nucleophilic attack with a slightly faster performance in the 3′‐dGMP branch than in the case of 5′‐dGMP with activation barriers of 21.1 and 20.4 kcal mol^?1 (experimental: 23.8 and 23.2 kcal mol^?1). The reduction reaction is connected with an electron flow from guanine. The product of the reduction reaction is a chelate structure, which dissociates within the last reaction step, that is, a hydrolysis reaction. The whole redox process (substitution, reduction, and hydrolysis) is exergonic by 34 and 28 kcal mol^?1 for 5′‐dGMP and 3′‐dGMP, respectively.     

Heterolysis of H2 Across a Classical Lewis Pair, 2,6‐Lutidine⋅BCl3: Synthesis,Characterization, and Mechanism

We report that 2,6‐lutidine?trichloroborane (Lut?BCl₃) reacts with H₂ in toluene, bromobenzene, dichloromethane, and Lut solvents producing the neutral hydride, Lut?BHCl₂. The mechanism was modeled with density functional theory, and energies of stationary states were calculated at the G3(MP2)B3 level of theory. Lut?BCl₃ was calculated to react with H₂ and form the ion pair, [LutH⁺][HBCl₃^?], with a barrier of ΔH^≠=24.7 kcal mol^?1 (ΔG^≠=29.8 kcal mol^?1). Metathesis with a second molecule of Lut?BCl₃ produced Lut?BHCl₂ and [LutH⁺][BCl₄^?]. The overall reaction is exothermic by 6.0 kcal mol^?1 (Δ_rG°=?1.1). Alternate pathways were explored involving the borenium cation (LutBCl₂⁺) and the four‐membered boracycle [(CH₂{NC₅H₃Me})BCl₂]. Barriers for addition of H₂ across the Lut/LutBCl₂⁺ pair and the boracycle B?C bond are substantially higher (ΔG^≠=42.1 and 49.4 kcal mol^?1, respectively), such that these pathways are excluded. The barrier for addition of H₂ to the boracycle B?N bond is comparable (ΔH^≠=28.5 and ΔG^≠=32 kcal mol^?1). Conversion of the intermediate 2‐(BHCl₂CH₂)‐6‐Me(C₅H₃NH) to Lut?BHCl₂ may occur by intermolecular steps involving proton/hydride transfers to Lut/BCl₃. Intramolecular protodeboronation, which could form Lut?BHCl₂ directly, is prohibited by a high barrier (ΔH^≠=52, ΔG^≠=51 kcal mol^?1).     

Ethylene Epoxidation with Nitrous Oxide over Fe–BTC Metal–Organic Frameworks: A DFT Study

The epoxidation of ethylene with N₂O over the metal‐organic framework Fe–BTC (BTC=1,3,5‐benzentricarboxylate) is investigated by means of density functional calculations. Two reaction paths for the production of ethylene oxide or acetaldehyde are systematically considered in order to assess the efficiency of Fe–BTC for the selective formation of ethylene oxide. The reaction starts with the decomposition of N₂O to form an active surface oxygen atom on the Fe site of Fe–BTC, which subsequently reacts with an ethylene molecule to form an ethyleneoxy intermediate. This intermediate can then be selectively transformed either by 1,2‐hydride shift into the undesired product acetaldehyde or into the desired product ethylene oxide by way of ring closure of the intermediate. The production of ethylene oxide requires an activation energy of 5.1 kcal mol^?1, which is only about one‐third of the activation energy of acetaldehyde formation (14.3 kcal mol^?1). The predicted reaction rate constants for the formation of ethylene oxide in the relevant temperature range are approximately 2–4 orders of magnitude higher than those for acetaldehyde. Altogether, the results suggest that Fe–BTC is a good candidate catalyst for the epoxidation of ethylene by molecular N₂O.     

Gas‐Phase Energetics of Reductive Elimination from a Palladium(II) N‐Heterocyclic Carbene Complex

Energy‐resolved collision‐induced dissociation experiments using tandem mass spectrometry are reported for an phenylpalladium N‐heterocyclic carbene (NHC) complex. Reductive elimination of an NHC ligand as a phenylimidazolium ion involves a barrier of 30.9(14) kcal mol^?1, whereas competitive ligand dissociation requires 47.1(17) kcal mol^?1. The resulting three‐coordinate palladium complex readily undergoes reductive C? C coupling to give the phenylimidazolium π complex, for which the binding energy was determined to be 38.9(10) kcal mol^?1. Density functional calculations at the M06‐L//BP86/TZP level of theory are in very good agreement with experiment. In combination with RRKM modeling, these results suggest that the rate‐determining step for the direct reductive elimination process switches from the C? C coupling step to the fragmentation of the resulting σ complex at low activation energy.     

The Unexpected Mechanism Underlying the High‐Valent Mono‐Oxo‐Rhenium(V) Hydride Catalyzed Hydrosilylation of CN Functionalities: Insights from a DFT Study

In this study, we theoretically investigated the mechanism underlying the high‐valent mono‐oxo‐rhenium(V) hydride Re(O)HCl₂(PPh₃)₂ ( 1 ) catalyzed hydrosilylation of C?N functionalities. Our results suggest that an ionic S_N2‐Si outer‐sphere pathway involving the heterolytic cleavage of the Si?H bond competes with the hydride pathway involving the C?N bond inserted into the Re?H bond for the rhenium hydride ( 1 ) catalyzed hydrosilylation of the less steric C?N functionalities (phenylmethanimine, PhCH=NH, and N‐phenylbenzylideneimine, PhCH=NPh). The rate‐determining free‐energy barriers for the ionic outer‐sphere pathway are calculated to be ～28.1 and 27.6 kcal mol^?1, respectively. These values are slightly more favorable than those obtained for the hydride pathway (by ～1–3 kcal mol^?1), whereas for the large steric C?N functionality of N,1,1‐tri(phenyl)methanimine (PhCPh=NPh), the ionic outer‐sphere pathway (33.1 kcal mol^?1) is more favorable than the hydride pathway by as much as 11.5 kcal mol^?1. Along the ionic outer‐sphere pathway, neither the multiply bonded oxo ligand nor the inherent hydride moiety participate in the activation of the Si?H bond.     

Theoretical Elucidation of the Mechanism of Cleavage of the Aromatic C?C Bond in Quinoxaline by a Tungsten‐Based Complex [W(PMe3)4(η2‐CH2PMe2)H]

The aromatic C? C bond cleavage by a tungsten complex reported recently by Sattler and Parkin 15 offers fresh opportunities for the functionalization of organic molecules. The mechanism of such a process has not yet been determined, which appeals to computational assistance to understand how the unstrained C? C bond is activated at the molecular level. 16 , 17 In this work, by performing density functional theory calculations, we studied various possible mechanisms of cleavage of the aromatic C? C bond in quinoxaline (QoxH) by the W‐based complex [W(PMe₃)₄(η²‐CH₂PMe₂)H]. The calculated results show that the mechanism proposed by Sattler and Parkin involves an overall barrier of as high as 42.0 kcal mol^?1 and thus does not seem to be consistent with the experimental observation. Alternatively, an improved mechanism has been presented in detail, which involves the removal and recoordination of a second PMe₃ ligand on the tungsten center. In our new mechanism, it is proposed that the C? C cleavage occurs prior to the second C? H bond addition, in contrast to Sattler and Parkin’s mechanism in which the C? C bond is broken after the second C? H bond addition. We find that the rate‐determining step of the reaction is the ring‐opening process of the tungsten complex with an activation barrier of 28.5 kcal mol^?1 after the first PMe₃ ligand dissociation from the metal center. The mono‐hydrido species is located as the global minimum on the potential‐energy surface, which is in agreement with the experimental observation for this species. The present theoretical results provide new insight into the mechanism of the remarkable C? C bond cleavage.     

Probing Reaction Mechanism of [1,5]‐Migration in Pyrrolium and Pyrrole Derivatives: Activation of a Stronger Bond in Electropositive Groups Becomes Easier

The [1,5]‐migration reaction has attracted considerable attention from experimentalists and theoreticians for decades. Although it has been extensively investigated in various systems, studies on pyrrolium derivatives are underdeveloped. Herein, a theoretical study on the reaction mechanism of [1,5]‐migration in both pyrrolium and pyrrole derivatives is presented. The results reveal lower activation barriers in [1,5]‐migration of electropositive groups (AuPMe₃ and SnH₃) in pyrrolium derivatives, although the bond dissociation energies of the Au?N bond (98.8 kcal mol^?1) and Sn?N bond (81.7 kcal mol^?1) are larger than that of the N?F bond (57.6 kcal mol^?1). The unexpectedly lower activation barriers (4.5 and 4.9 kcal mol^?1 for AuPMe₃ and SnH₃, respectively) for [1,5]‐migration of electropositive groups, in comparison with the [1,5]‐fluorine shift, can be attributed to aromaticity stabilizing the transition states, as revealed by significantly negative nucleus‐independent chemical shift (NICS) values. Further studies indicate that charge distribution and frontier molecular orbitals also play some roles in [1,5]‐migration of pyrrolium derivatives.     

Mechanism for the gas‐phase hydrogen fluoride‐mediated decomposition of peroxyacetyl nitrate (PAN) studied by DFT method

Density functional theory has been used to study the mechanism of the decomposition of peroxyacetyl nitrate (CH₃C(O)OONO₂) in hydrogen fluoride clusters containing one to three hydrogen fluoride molecules at the B3LYP/6‐311++G(d,p) and B3LYP/6‐311+G(3df,3pd) levels. The calculations clarify some of the uncertainties in the mechanism of PAN decomposition in the gas phase. The energy barrier decreases from 30.5 kcal mol^?1 (single hydrogen fluoride) to essentially 18.5 kcal mol^?1 when catalyzed by three hydrogen fluoride molecules. As the size of the hydrogen fluoride cluster is increased, PAN shows increasing ionization along the O? N bond, consistent with the proposed predissociation in which the electrophilicity of the nitrogen atom is enhanced. This reaction is found to proceed through an attack of a fluorine to the PAN nitrogen in concert with a proton transfer to a PAN oxygen. On the basis of our calculations, an alternative reaction mechanism for the decomposition of PAN is proposed. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Computational study on the mechanism for the gas‐phase reaction of dimethyl disulfide with OH

The mechanisms for the reaction of CH₃SSCH₃ with OH radical are investigated at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311++G(d,p) level of theory. Five channels have been obtained and six transition state structures have been located for the title reaction. The initial association between CH₃SSCH₃ and OH, which forms two low‐energy adducts named as CH₃S(OH)SCH₃ (IM1 and IM2), is confirmed to be a barrierless process, The S? S bond rupture and H? S bond formation of IM1 lead to the products P1(CH₃SH + CH₃SO) with a barrier height of 40.00 kJ mol^?1. The reaction energy of Path 1 is ?74.04 kJ mol^?1. P1 is the most abundant in view of both thermodynamics and dynamics. In addition, IMs can lead to the products P2 (CH₃S + CH₃SOH), P3 (H₂O + CH₂S + CH₃S), P4 (CH₃ + CH₃SSOH), and P5 (CH₄ + CH₃SSO) by addition‐elimination or hydrogen abstraction mechanism. All products are thermodynamically favorable except for P4 (CH₃ + CH₃SSOH). The reaction energies of Path 2, Path 3, Path 4, and Path 5 are ?28.42, ?46.90, 28.03, and ?89.47 kJ mol^?1, respectively. Path 5 is the least favorable channel despite its largest exothermicity (?89.47 kJ mol^?1) because this process must undergo two barriers of TS5 (109.0 kJ mol^?1) and TS6 (25.49 kJ mol^?1). Hopefully, the results presented in this study may provide helpful information on deep insight into the reaction mechanism. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011