Theoretical studies on cycloaddition reactions between the 2-aza-1,3-butadiene cation and olefins

Density functional (B3LYP) calculations, using the 6-31G basis set, have been employed to study the title reactions. For the model reaction (H(2)C=C-NH(+)=CH(2) + H(2)C=CH(2)), a complex has been formed with 6.2 kcal/mol of stabilization energy and the transition state is 4.0 kcal/mol above this complex, but 2.1 kcal/mol below the reactants. However, the substituent effects are quite remarkable. When ethene is substituted by electron-withdrawing group CN, the reaction could also yield six-membered-ring products, but the energy barriers are all more than 7 kcal/mol, which shows that CN group unfavors the reaction. The other substituents, such as CH(3)O and CH(3) groups, have also been considered in the present work, and the results show that they are favorable for the formation of six-membered-ring adducts. The calculated results have been rationalized with frontier orbital interaction and topological analysis.     

Mechanisms of the reactions of W AND W+ with H2O: computational studies

The mechanism of the reactions of W and W(+) with the water molecule have been studied for several lower-lying electronic states of tungsten centers at the CCSD(T)/6-311G(d,p)+SDD and B3LYP/6-31G(d,p)+SDD levels of theory. It is shown that these reactions are essentially multistate processes, during which lower-lying electronic states of the systems cross several times. They start with the formation of initial prereaction M(H(2)O) complexes with M-H(2)O bonding energies of 9.6 and 48.2 kcal/mol for M = W and W(+), followed by insertion of the metal center into an O-H bond with 20.0 and 53.3 kcal/mol barriers for neutral and cationic systems, respectively. The overall process of M + H(2)O --> t-HM(OH) is calculated to be highly exothermic, 48.4 and 48.8 kcal/mol for M = W and W(+). From the HM(OH) intermediate the reaction may proceed via several different channels, among which the stepwise HM(OH) --> HMO + H --> (H)(2)MO and concerted HM(OH) --> (H)(2)MO pathways are more favorable and can compete (energetically) with each other. For the neutral system (M = W), the concerted process is the most favorable, whereas for the charged system (M = W(+)), the stepwise pathway is slightly more favorable. From the energetically most favorable intermediate (H)(2)MO the reactions proceed via H(2)-molecule formation with a 53.1 kcal/mol activation barrier for the neutral system. For the cationic system, H-H formation and dissociation is an almost barrierless process. The overall reaction of W and W(+) with the water molecule leading to H(2) + MO formation is found to be exothermic by 48.2 and 39.8 kcal/mol, respectively. In the gas phase with the collision-less conditions the reactions W((7)S) + H(2)O --> H(2) + WO((3)Sigma(+)), and W(+)((6)D) + H(2)O --> H(2) + WO(+)((4)Sigma(+)) are expected to proceed via a 10.4 and 5.1 kcal/mol overall energy barrier corresponding to the first O-H dissociation at the TS1. On the basis of these PESs, we predict kinetic rate constants for the reactions of W and W(+) with H(2)O.     

Quantum chemical simulation of the reaction of methane with gold(I) acetylacetonate and aquaacetylacetonate complexes

The interaction between methane and gold(I) acetylacetonate via electrophilic substitution (reaction (I)) and oxidative addition (reaction (II)) is simulated. In both cases, the formation of the products is thermodynamically favorable: the decrease in energy is 31 kcal/mol for reaction (I) and 26 kcal/mol for reaction (II). The product of reaction (II) is additionally stabilized by Au-H interaction. Both reactions have a low activation barrier and proceed via the formation of structurally different methane complexes reducing the energy of the system by 9.3 kcal/mol for reaction (I) and by 10.9 kcal/mol for reaction (II). The complex [Au(H₂O)(acac)] is also capable of forming methane complexes. These complexes result from a thermally neutral reaction and turn into products after overcoming a low energy barrier. The structure of the complex activating methane in the gold-rutin system is deduced from the data obtained.     

Chemo- and periselectivity in the addition of [OsO2(CH2)2] to ethylene: a theoretical study

Quantum chemical calculations by using density functional theory at the B3LYP level have been carried out to elucidate the reaction course for the addition of ethylene to [OsO2(CH2)2] (1). The calculations predict that the kinetically most favorable reaction proceeds with an activation barrier of 8.1 kcal mol(-1) via [3+2] addition across the O=Os=CH2 moiety. This reaction is -42.4 kcal mol(-1) exothermic. Alternatively, the [3+2] addition to the H2C=Os=CH2 fragment of 1 leads to the most stable addition product 4 (-72.7 kcal mol(-1)), yet this process has a higher activation barrier (13.0 kcal mol(-1)). The [3+2] addition to the O=Os=O fragment yielding 2 is kinetically (27.5 kcal mol(-1)) and thermodynamically (-7.0 kcal mol(-1)) the least favorable [3+2] reaction. The formal [2+2] addition to the Os=O and Os=CH2 double bonds proceeds by initial rearrangement of 1 to the metallaoxirane 1 a. The rearrangement 1-->1 a and the following [2+2] additions have significantly higher activation barriers (>30 kcal mol(-1)) than the [3+2] reactions. Another isomer of 1 is the dioxoosmacyclopropane 1 b, which is 56.2 kcal mol(-1) lower in energy than 1. The activation barrier for the 1-->1 b isomerization is 15.7 kcal mol(-1). The calculations predict that there are no energetically favorable addition reactions of ethylene with 1 b. The isomeric form 1 c containing a peroxo group is too high in energy to be relevant for the reaction course. The accuracy of the B3LYP results is corroborated by high level post-HF CCSD(T) calculations for a subset of species.     

Reliable predictions of the thermochemistry of boron-nitrogen hydrogen storage compounds: BxNxHy, x = 2, 3

Thermochemical data calculated using ab initio molecular orbital theory are reported for 16 BxNxHy compounds with x = 2, 3 and y > or = 2x. Accurate gas-phase heats of formation were obtained using coupled cluster with single and double excitations and perturbative triples (CCSD(T)) valence electron calculations extrapolated to the complete basis set (CBS) limit with additional corrections including core/valence, scalar relativistic, and spin-orbit corrections to predict the atomization energies and scaled harmonic frequencies to correct for zero point and thermal energies and estimate entropies. Computationally cheaper calculations were also performed using the G3MP2 and G3B3 variants of the Gaussian 03 method, as well as density functional theory (DFT) using the B3LYP functional. The G3MP2 heats of formation are too positive by up to approximately 6 kcal/mol as compared with CCSD(T)/CBS values. The more expensive G3B3 method predicts heats of formation that are too negative as compared with the CCSD(T)/CBS values by up to 3-4 kcal/mol. DFT using the B3LYP functional and 6-311+G** basis set predict isodesmic reaction energies to within a few kcal/mol compared with the CCSD(T)/CBS method so isodesmic reactions involving BN compounds and the analogous hydrocarbons can be used to estimate heats of formation. Heats of formation of c-B3N3H12 and c-B3N3H6 are -95.5 and -115.5 kcal/mol at 298 K, respectively, using our best calculated CCSD(T)/CBS approach. The experimental value for c-B3N3H6 appears to be approximately 7 kcal/mol too negative. Enthalpies, entropies, and free energies are calculated for many dehydrocoupling and dehydrogenation reactions that convert BNH6 to alicyclic and cyclic oligomers and H2(g). Generally, the reactions are highly exothermic and exergonic as well because of the release of 1 or more equivalents of H2(g). For c-B3N3H12 and c-B3N3H6, available experimental data for sublimation and vaporization lead to estimates of their condensed phase 298 K heats of formation: DeltaHf degrees [c-B3N3H12(s)] = -124 kcal/mol and DeltaHf degrees [c-B3N3H6(l)] = -123 kcal/mol. The reaction thermochemistries for the dehydrocoupling of BNH6(s) to c-B3N3H12(s) and the dehydrogenation of c-B3N3H12(s) to c-B3N3H6(l) are much less exothermic compared with the gas-phase reactions due to intermolecular forces which decrease in the order BNH6 > cyclo-B3N3H12 > cyclo-B3N3H6. The condensed phase reaction free energies are less negative compared with the gas-phase reactions but are still too favorable for BNH6 to be regenerated from either c-B3N3H12 or c-B3N3H6 by just an overpressure of H2.     

Computational study of reaction pathways for the formation of indium nitride from trimethylindium with HN3: comparison of the reaction with NH3 and that on TiO2 rutile (110) surface

The reactions of trimethylindium (TMIn) with HN3 and NH3 are relevant to the chemical vapor deposition of indium nitride thin film. The mechanisms and energetics of these reactions in the gas phase have been investigated by density functional theory and ab initio calculations using the CCSD(T)/Lanl2dz//B3LYP/Lanl2dz and CCSD(T)/Lanl2dz//MP2/Lanl2dz methods. The results of both methods are in good agreement for the optimized geometries and relative energies. These results suggest that the reaction with HN3 forms a new stable product, dimethylindiumnitride, CH3-In=N-CH3 via another stable In(CH3)2N3 (dimethylindium azide, DMInA) intermediate. DMInA may undergo unimolecular decomposition to form CH3InNCH3 by two main possible pathways: (1) a stepwise decomposition process through N2 elimination followed by CH3 migration from In to the remaining N atom and (2) a concerted process involving the concurrent CH3 migration and N2 elimination directly giving N2+CH3InNCH3. The reaction of TMIn with NH3 forms a most stable product DMInNH2 following the initial association and CH4-elimination reaction. The required energy barrier for the elimination of the second CH4 molecule from DMInNH2 is 74.2 kcal/mol. Using these reactions, we predict the heats of formation at 0 K for all the products and finally for InN which is 123+/-1 kcal/mol predicted by the two methods. The gas-phase reaction of HN3 with TMIn is compared with that occurring on rutile TiO2 (110). The most noticeable difference is the high endothermicity of the gas-phase reaction for InN production (53 kcal/mol) and the contrasting large exothermicity (195 kcal/mol) released by the low-barrier Langmuir-Hinshelwood type processes following the adsorption of TMIn and HN3 on the surface producing a horizontally adsorbed InN(a), Ti-NIn-O(a), and other products, CH4(g)+N2(g)+2CH3O(a) [J. Phys. Chem. B 2006, 110, 2263].     

Unimolecular dissociation of the propargyl radical intermediate of the CH+C2H2 and C+C2H3 reactions

This paper examines the unimolecular dissociation of propargyl (HCCCH2) radicals over a range of internal energies to probe the CH+HCCH and C+C2H3 bimolecular reactions from the radical intermediate to products. The propargyl radical was produced by 157 nm photolysis of propargyl chloride in crossed laser-molecular beam scattering experiments. The H-loss and H2 elimination channels of the nascent propargyl radicals were observed. Detection of stable propargyl radicals gave an experimental determination of 71.5 (+5-10) kcal/mol as the lowest barrier to dissociation of the radical. This barrier is significantly lower than predictions for the lowest barrier to the radical's dissociation and also lower than calculated overall reaction enthalpies. Products from both H2+HCCC and H+C3H2 channels were detected at energies lower than what has been theoretically predicted. An HCl elimination channel and a minor C-H fission channel were also observed in the photolysis of propargyl chloride.     

Thermolysis of alkyl sulfoxides and derivatives: a comparison of experiment and theory.

Gas-phase activation data were obtained for model sulfoxide elimination reactions. The activation enthalpy for methyl 3-phenylpropyl sulfoxide is 32.9 +/- 0.9 kcal/mol. Elimination by methyl vinyl sulfoxide to form acetylene has an enthalpic barrier of 41.6 +/- 0.8 kcal/mol and that of 3-phenylpropyl methanesulfinate to form hydrocinnamaldehyde is 34.6 +/- 0.6 kcal/mol. Calculations at the MP2/6-311+G(3df,2p)//MP2/6-31G(d,p) level for simplified models of these reactions provide barriers of 32.3, 40.3, and 32.7 kcal/mol, respectively. A series of other compounds are examined computationally, and it is shown that the substituent effects on the sulfoxide elimination reaction are much more straightforward to interpret if DeltaH data are available in addition to the usually determined DeltaH++. The activation enthalpy of the reverse addition reaction is also subject to structural variation and can usually be rationalized on the basis of nucleophilicity of the sulfur or polarity matching between the sulfenic acid and olefin derivative.     

Reaction HFCO + 2H2O — Calculated Mechanisms

Three possible reaction mechanisms of methanoyl fluoride with 2H₂O include a concerted and a stepwise hydrolysis of HFCO into HCOOH + HF, and a pure catalytic decomposition of HFCO into HF + CO. Among these, the two H₂O molecules acting as catalyst to decompose HFCO has the lowest calculated barrier, 25.1 kcal/mol with respect to the reactant‐adduct complex, whereas the barriers for the concerted and stepwise hydrolytic reactions in which one H₂O acts as a reactant and the other H₂O as catalyst are similar, 30.8 kcal/mol for concerted and 29.9 kcal/mol for stepwise. The formation of transoid HCOOH in the hydrolysis of HFCO is more favorable than cisoid HCOOH.     

Theoretical analysis of concerted and stepwise mechanisms of the hetero-Diels–Alder reaction of butadiene with formaldehyde and thioformaldehyde

The concerted and stepwise mechanisms of the hetero-Diels–Alder reaction of butadiene with formaldehyde and thioformaldehyde were studied by a CASSCF molecular orbital method. The energy barrier of the concerted reaction of butadiene with formaldehyde is about 21 kcal/mol higher than that of butadiene with thioformaldehyde at the CAS-MP2 calculation level. For the stepwise reaction paths, the energy barrier for the first step process of the reaction of butadiene with formaldehyde is about 17 kcal/mol above that of butadiene with thioformaldehyde. The concerted pathways for both systems are more favorable by 9–12 kcal/mol than the stepwise pathways. The electronic mechanisms for the concerted reactions of both reaction systems are also discussed by a CiLC analysis.     

DFT study on the mechanism of water-assisted dihydrogen elimination in group 6 octahedral metal hydride complexes

The reaction of water with octahedral bis-, tris- and tetrakis-(phosphine)tungsten, (phosphine)molybdenum and (phosphine)chromium complexes has been studied using B3LYP/def2-TZVPP level of DFT to elucidate dissociative, associative and hydride migratory insertion mechanisms for hydrogen elimination. In the dissociative mechanism, phosphine dissociation requires 19.3-28.5 kcal mol(-1) of energy. The phosphine-water ligand exchange is endergonic due to poor coordination ability of water to group 6 metals (binding energy 8.8-15.5 kcal mol(-1)). The ligand exchange leads to intermolecular M-HH(2)O dihydrogen interaction and facilitates dihydrogen elimination (G(act) = 6.8-15.5 kcal mol(-1)). In the associative mechanism, a water molecule in the first solvation shell interacts with the M-H bond through a dihydrogen bond (interaction energy 2.7-4.0 kcal mol(-1)) and leads to the elimination of H(2) by forming a hydroxide complex. Compared to the dissociative mechanism, G(act) of associative mechanisms are ～22 kcal mol(-1) higher. In the hydride migratory insertion mechanism, the hydride ligand shifts to the CO ligand (G(act) = 25.4-30.4 kcal mol(-1)) to afford a formyl complex and subsequently the H-H bond coupling occurs between formyl and water ligand (G(act) = 2.8-4.4 kcal mol(-1)). In many cases, the migratory insertion mechanism can simultaneously operate with the dissociative mechanism as a minor pathway, whereas owing to high G(act) value, the associative mechanism can be described as inactive in the reaction. The general argument that dihydrogen elimination is preceded by the formation of a dihydrogen intermediate is not applicable for the systems studied herein as the most favoured dissociative mechanism does not pass through such an intermediate. On the other hand, irrespective of the mechanisms, dihydrogen elimination invariably occurs with the formation of a dihydrogen bonded transition state. Our results also suggest that group 6 octahedral metal hydride complexes are attractive targets for the design of water splitting reactions.     

Mechanisms of the reactions of W and W+ with NOx (x=1, 2): a computational study

The mechanisms of the reactions of W and W+ with NOx (x=1, 2) were studied at the CCSD(T)/[SDD+6-311G(d)]//B3LYP/[SDD+6-31G(d)] level of theory. It was shown that the insertion pathway of the reaction W(7S)+NO2(2A1) is a multistate process, which involves several lower lying electronic states of numerous intermediates and transition states, and leads to oxidation, WO(3Sigma)+NO(2Pi), and/or nitration, WN(4Sigma)+O2(3Sigmag-), of the W-center. Oxidation products WO(3Sigma)+NO(2Pi) lie 87.6 kcal/mol below the reactants, while the nitration channel is only 31.0 kcal/mol exothermic. Furthermore, it was shown that nitration of W with NO2 is kinetically less favorable than its oxidation. The addition-dissociation pathway of the reaction W(7S)+NO2(2A1) proceeds via the octet (ground) state potential energy surface of the reaction, requires 3.3 kcal/mol barrier, and leads exclusively to oxidation products. Calculations show that oxidation of the W+ cation by NO2 is a barrierless process in the gas phase, proceeds exclusively via the insertion pathway, and is exothermic by 82.9 kcal/mol. The nitration of W+ by NO2 is only 14.1 kcal/mol exothermic and could be accessible only under high-temperature conditions. Reactions of M=W/W+ with NO are also barrierless processes in the gas phase and lead to the N-O insertion product NMO, which are 105.4 and 77.4 kcal/mol lower than the reactants for W and W+, respectively.     

The reaction of tricarbon with acetylene: an ab initio/RRKM study of the potential energy surface and product branching ratios

Ab initio calculations of the potential energy surface for the C3(1Sigmag+)+C2H2(1Sigmag+) reaction have been performed at the RCCSD(T)/cc-pVQZ//B3LYP/6-311G(d,p) + ZPE[B3LYP/6-311G(d,p)] level with extrapolation to the complete basis set limit for key intermediates and products. These calculations have been followed by statistical calculations of reaction rate constants and product branching ratios. The results show the reaction to begin with the formation of the 3-(didehydrovinylidene)cyclopropene intermediate i1 or five-member ring isomer i7 with the entrance barriers of 7.6 and 13.8 kcal/mol, respectively. i1 rearranges to the other C5H2 isomers, including ethynylpropadienylidene i2, singlet pentadiynylidene i3, pentatetraenylidene i4, ethynylcyclopropenylidene i5, and four- and five-member ring structures i6, i7, and i8 by ring-closure and ring-opening processes and hydrogen migrations. i2, i3, and i4 lose a hydrogen atom to produce the most stable linear isomer of C5H with the overall reaction endothermicity of approximately 24 kcal/mol. H elimination from i5 leads to the formation of the cyclic C5H isomer, HC2C3, +H, 27 kcal/ mol above C3+C2H2. 1,1-H2 loss from i4 results in the linear pentacarbon C5+H2 products endothermic by 4 kcal/mol. The H elimination pathways occur without exit barriers, whereas the H2 loss from i4 proceeds via a tight transition state 26.4 kcal/mol above the reactants. The characteristic energy threshold for the reaction under single collision conditions is predicted be in the range of approximately 24 kcal/mol. Product branching ratios obtained by solving kinetic equations with individual rate constants calculated using RRKM and VTST theories for collision energies between 25 and 35 kcal/mol show that l-C5H+H are the dominant reaction products, whereas HC2C3+H and l-C5+H2 are minor products with branching ratios not exceeding 2.5% and 0.7%, respectively. The ethynylcyclopropenylidene isomer i5 is calculated to be the most stable C5H2 species, more favorable than triplet pentadiynylidene i3t by approximately 2 kcal/mol.     

Theoretical study on reactions of nitroethylene with the Si(100)-2 x 1 surface

Possible reaction pathways of nitroethylene with the Si(100)-2 x 1 surface have been investigated by unrestricted density functional theory. The facile occurrence of the studied reactions was demonstrated by the low activation energies of the rate-determining steps (1.07-5.23 kcal/mol). It was found that the [4 + 2] cycloaddition reaction of nitroethylene is most kinetically favorable. The isomerization reactions of the addition products were also investigated. The [3 + 2] cycloaddition product may further undergo a rearrangement by overcoming a 12.37 kcal/mol activation energy barrier into an isomer, with an oxygen atom of the nitryl group inserted between two silicon atoms of the Si(100) surface.     

Computational study on the mechanisms and energetics of trimethylindium reactions with H2O and H2S

The reactions of trimethylindium (TMIn) with H2O and H2S are relevant to the chemical vapor deposition of indium oxide and indium sulfide thin films. The mechanisms and energetics of these reactions in the gas phase have been investigated by density functional theory and ab initio calculations using the CCSD(T)/[6-31G(d,p)+Lanl2dz]//B3LYP/[6-31G(d,p)+Lanl2dz] and CCSD(T)/[6-31G(d,p)+Lanl2dz] //MP2/[6-31G(d,p)+Lanl2dz] methods. The results of both methods are in good agreement for the optimized geometries and relative energies. When TMIn reacts with H2O and H2S, initial molecular complexes [(CH3)3In:OH2 (R1)] and [(CH3)3In:SH2 (R2)] are formed with 12.6 and 3.9 kcal/mol binding energies. Elimination of a CH4 molecule from each complex occurs with a similar energy barrier at TS1 (19.9 kcal/mol) and at TS3 (22.1 kcal/mol), respectively, giving stable intermediates (CH3)2InOH and (CH3)2InSH. The elimination of the second CH4 molecule from these intermediate products, however, has to overcome very high and much different barriers of 66.1 and 53.2 kcal/mol, respectively. In the case of DMIn with H2O and H2S reactions, formation of both InO and InS is exothermic by 3.1 and 30.8 kcal/mol respectively. On the basis of the predicted heats of formation of R1 and R2 at 0 K and -20.1 and 43.6 kcal/mol, the heats of formation of (CH3)2InOH, (CH3)2InSH, CH3InO, CH3InS, InO, and InS are estimated to be -20.6, 31.8, and 29.0 and 48.4, 35.5, and 58.5 kcal/mol, respectively. The values for InO and InS are in good agreement with available experimental data. A similar study on the reactions of (CH3)2In with H2O and H2S has been carried out; in these reactions CH3InOH and CH3InSH were found to be the key intermediate products.     

Radical-molecule reactions HCO/HOC + C2H2: mechanistic study

A detailed computational study is performed on the unknown radical-molecule reactions between HCO/HOC and acetylene (C2H2) at the CCSD(T)/6-311G(2d,p)//B3LYP/6-311G(d,p)+ZPVE, Gaussian-3//B3LYP/6-31G(d), and Gaussian-3//MP2(full)/6-31G(d) levels. For the HCO + C2H2 reaction, the most favorable pathway is direct C-addition forming the intermediate HC=CHCH=O followed by a 1,3-H-shift leading to H2C=CHC=O, which finally dissociates to the product C2H3 + CO. The overall reaction barrier is 13.8, 10.5, and 11.3 kcal/mol, respectively, at the three levels. The quasi-direct H-donation process to produce C2H3 + CO with barriers of 14.0, 14.1, and 14.1 kcal/mol is less competitive. Thus only at higher temperatures could the HCO + C2H2 reaction play a role. In contrast, the HOC + C2H2 reaction can barrierlessly generate C2H3 + CO via the quasi-direct H-donation mechanism proceeding via a prereactive complex with OH...C2 hydrogen bonding. This is suggestive of the potential importance of the HOC + C2H2 reaction in both combustion and interstellar processes. However, the direct C-addition channel is much less competitive. For both reactions, the possible formation of the intriguing interstellar molecules propadiene and propynal is also discussed. The present theoretical study represents the first attempt to probe the reaction mechanism between HOC and pi-systems. Future laboratory investigations on both reactions (particularly HOC + C2H2) are recommended.     

Pathway from N-Alkylglycine to Alkylisonitrile Catalyzed by Iron(II) and 2-Oxoglutarate-Dependent Oxygenases

N-alkylisonitrile, a precursor to isonitrile-containing lipopeptides, is biosynthesized by decarboxylation-assisted -N≡C group (isonitrile) formation by using N-alkylglycine as the substrate. This reaction is catalyzed by iron(II) and 2-oxoglutarate (Fe/2OG) dependent enzymes. Distinct from typical oxygenation or halogenation reactions catalyzed by this class of enzymes, installation of the isonitrile group represents a novel reaction type for Fe/2OG enzymes that involves a four-electron oxidative process. Reported here is a plausible mechanism of three Fe/2OG enzymes, Sav607, ScoE and SfaA, which catalyze isonitrile formation. The X-ray structures of iron-loaded ScoE in complex with its substrate and the intermediate, along with biochemical and biophysical data reveal that -N≡C bond formation involves two cycles of Fe/2OG enzyme catalysis. The reaction starts with an Fe^IV-oxo-catalyzed hydroxylation. It is likely followed by decarboxylation-assisted desaturation to complete isonitrile installation.     

Theoretical analysis of concerted and stepwise mechanisms of Diels-Alder reactions of butadiene with silaethylene and disilene

The concerted and the stepwise mechanisms of the Diels-Alder reactions of butadiene with silaethylene and disilene were studied by ab initio MO methods. For the reaction of butadiene and silaethylene, an asymmetric concerted process that is almost stepwise and two stepwise processes were located. For the first step of the stepwise process, the C-Si bond formation is more favorable than the C-C bond formation. The activation energy barrier of the concerted transition state is only 0.89 kcal/mol lower than that of the first-step transition state of the C-Si bond formation for the stepwise process by the CASPT2 calculation level. For the reaction of butadiene and disilene, the activation energy barrier of the concerted-type transition state constrained with Cs symmetry is about 9 kcal/mol higher than that of the stepwise transition state by the CASSCF method. The energy barrier of the first step of the stepwise reaction disappears at the CASPT2/6-311++G(d,p) calculation level including the nondynamical correlation energy, although the reaction of the butadiene with disilene occurs through the stepwise-like process.     

Effects of ion-pairing and hydration on the SNAr reaction of the F- with p-chlorobenzonitrile in aprotic solvents

Theoretical ab initio calculations including liquid phase optimizations were used to investigate the S(N)Ar reaction of the fluoride ion with p-chlorobenzonitrile in dimethyl sulfoxide solution. The effect of the counter ion and hydration of the fluoride ion with one water molecule was analyzed. The calculations indicate that the gas-phase S(N)Ar reaction is more favorable than the corresponding S(N)2 reactions involving fluoride ion and 2-chlorobutane. However, the substantially higher solvent effect on the S(N)Ar reaction makes the nucleophilic substitution on the aromatic ring less favorable than the aliphatic reaction in the liquid phase. For the anhydrous tetrabutylammonium fluoride system, the theoretical free energy barrier of 22.1 kcal mol(-1) is close to the experimental one of 24.4 kcal mol(-1). The smaller tetramethylammonium cation strongly associates with the fluoride ion and increases the barrier by 5 kcal mol(-1). Similarly, just one water molecule hydrating the fluoride ion has the same effect. An analysis of the reaction involving the ion pair and the free anion in different polarity media predicts an unexpected behavior and indicates there is an ideal solvent polarity for each counter ion.     

Theoretical studies of proton-transfer reactions of 2-hydroxypyridine--(H2O)n (n = 0-2) in the ground and excited states

The potential energy profiles for proton-transfer reactions of 2-hydroxypyridine and its complexes with water were determined by MP2, CASSCF and MR-CI calculations with the 6-31G** basis set. The tautomerization reaction between 2-hydroxypyridine (2HP) and 2-pyridone (2PY) does not take place at room temperature because of a barrier of approximately 35 kcal/mol for the ground-state pathway. The water-catalyzed enol-keto tautomerization reactions in the ground state proceed easily through the concerted proton transfer, especially for the two-water complex. The S1 tautomerization between the 2HP and 2PY monomers has a barrier of 18.4 kcal/mol, which is reduced to 5.6 kcal/mol for the one-water complex and 6.4 kcal/mol for the two-water complex. The results reported here predict that the photoinduced tautomerization reaction between the enol and keto forms involves a cyclic transition state having one or two water molecules as a bridge.