Reinvestigation of the Deceptively Simple Reaction of Toluene with OH and the Fate of the Benzyl Radical: The "Hidden" Routes to Cresols and Benzaldehyde.

In a previous work, we have investigated the initial steps of the reaction of toluene with the hydroxyl radical using several quantum chemical approaches including density functional and composite post-Hartree-Fock models. Comparison of H-abstraction from the methyl group and additions at different positions of the phenyl ring showed that the former reaction channel is favored at room temperature. This conclusion appears at first sight incompatible with the experimental observation of a lower abundance of the product obtained from abstraction (benzaldehyde) with respect to those originating from addition (cresols). Further reactions of the intermediate radicals with oxygen, water, and additional OH radicals are explored in this paper through theoretical calculations on more than 120 species on the corresponding potential energy surface. The study of the addition reactions, to obtain the cresols through hydroxy methylcyclodienyl intermediate radicals, showed that only in the case of o-cresol the reaction proceeds by addition of O2 to the ring, internal H-transfer, and hydroperoxyl abstraction and not through direct H-abstraction. For both p- and m-cresol, instead, the reaction occurs through a higher-energy direct H-abstraction, thus explaining in part the observed larger concentration of the ortho isomer in the final products. It was also found that the benzyl radical, formed by H-abstraction from the methyl group, is able to react further if additional OH is present. Two reaction paths leading to o-cresol, two leading to p-cresol, and one leading to m-cresol were determined. Moreover, in this situation, the benzyl radical is predicted to produce benzyl alcohol, as was found in some experiments. The commonly accepted route to benzaldehyde was found to be not the energetically favored one. Instead, a route leading to the benzoyl radical (and ultimately to benzoic acid) with the participation of one water molecule was clearly more favorable, both thermodynamically and kinetically.

Computational Evidence Suggests That 1-Chloroethanol May Be an Intermediate in the Thermal Decomposition of 2-Chloroethanol into Acetaldehyde and HCl.

The dehalogenation of 2-chloroethanol (2ClEtOH) in the gas phase with and without the participation of catalytic water molecules has been investigated using methods rooted into the density functional theory. The well-known HCl elimination leading to vinyl alcohol (VA) was compared to the alternative elimination route toward oxirane and shown to be kinetically and thermodynamically more favorable. However, the isomerization of VA to acetaldehyde in the gas phase, in the absence of water, was shown to be kinetically and thermodynamically less favorable than the recombination of VA and HCl to form the isomeric 1-chloroethanol (1ClEtOH) species. At the ωB97X-D/cc-pVTZ level of calculation, this species is more stable than 2ClEtOH by about 6 kcal mol-1 at 298 K, and the reaction barrier for VA to 1ClEtOH is 23 kcal mol-1 versus 55 kcal mol-1 for the direct transformation of VA to acetaldehyde. In a successive step, 1ClEtOH can decompose directly to acetaldehyde and HCl with a lower barrier (29 kcal mol-1) than that of VA to the same products (55 kcal mol-1). The calculations were repeated using a single ancillary water molecule (W) in the complexes 2ClEtOH_W and 1ClEtOH_W. The latter adduct is now more stable than 2ClEtOH_W by about 8 kcal mol-1 at 298 K, implying that the water molecule increased the already higher stability of 1ClEtOH in the gas phase. However, this catalytic water molecule lowers dramatically the barrier for the interconversion of VA to acetaldehyde (from 55 to 7 kcal mol-1). This barrier is now smaller than the one for the conversion to 1ClEtOH (which also decreases, but not so much, from 23 to 13 kcal mol-1). Thus, it is concluded that while 1ClEtOH may be a plausible intermediate in the gas phase dehalogenation of 2ClEtOH, it is unlikely that it plays a major role in water complexes (or, by inference, aqueous solution). It is also shown that neither in the gas phase nor in the cluster with one water molecule, the oxirane path is more favorable than the VA alcohol path. Additionally, a direct conversion of 2ClEtOH to 1ClEtOH through a transition state which resembles a VA molecule in a complex with a chlorine atom and a hydrogen atom on both sides of this planar species was found. This reaction path has also lower activation energy than the conversion to oxirane but not as low as the conversion to VA.

Theoretical study on the mechanism of the reaction of CF(3)S with NO(2).

The singlet potential energy surface for the CF3S + NO2 reaction has been theoretically investigated using the B3LYP/6-311+G(3df) level of theory. The geometries, vibrational frequencies, and zero-point energies of all stationary points involved in the title reaction have been examined. More accurate energies of stationary points were obtained using ab initio G3//B3LYP and CBS-QB3 composite methods. The results show that the initial addition of CF3S with NO2 leads to CF3SNO2 or CF3SONO intermediates, which are formed without an electronic barrier. CF3SNO2 can easily isomerizes to CF3SONO, while CF3SONO readily isomerizes to CF3S(O)NO or dissociates to CF3SO + NO, which are the major products of the title reaction. Reaction channels leading to the formation the CF3O + SNO and CF2S + FNO2 products are highly improbable processes due to high energy barriers involved. We have also computed heats of formation for CF3SNO2, CF3SONO, and CF3S(O)NO intermediates. It was found that the most stable is the cis-perpendicular form of CF3SONO isomer with DeltaH(f,0)0 = -243.6 kcal mol-1.

Effect of halogenation on the mechanism of the atmospheric reactions between methylperoxy radicals and NO. A computational study.

The mechanism of the reactions between the halogenated methylperoxy radicals, CHX(2)O(2) (X = F, Cl), and NO is investigated by using ab initio and density functional quantum mechanical methods. Comparison is made with the mechanism of the CH(3)O(2) + NO reaction. The most important energy minima in the potential energy surface are found to be the two conformers of the halogenated methyl peroxynitrite association adducts, CHX(2)OONOcp and CHX(2)OONOtp, and the halogenated methyl nitrates, CHX(2)ONO(2). The latter are suggested to be formed through the one-step isomerization of the peroxynitrite adduct and may lead upon decomposition to carbonylated species, CX(2)O + HONO and CHXO + XNO(2). The ambiguous issue of the unimolecular peroxynitrite to nitrate isomerization is reconsidered, and the possibility of a triplet transition state involvement in the ROONOtp <--> RONO(2) rearrangement is examined. The overall calculations and the detailed correlation with the methyl system show the significant effect of the halogenation on the lowering of the entrance potential energy well which corresponds to the formation of the peroxynitrites. The increased attractive character of the potential energy surface found upon halogenation combined with the increased exothermicity of the CHX(2)O(2) + NO --> CHX(2)O + NO(2) reaction are suggested to be the important factors contributing to the enhanced reactivity of the halogenated reactions relative to CH(3)O(2) + NO. The calculated heat of formation values indicate the large stabilization of the fluorinated derivatives.

Ab initio characterization of (CH3IO3) isomers and the CH3O2 + IO reaction pathways.

The geometries, harmonic vibrational frequencies, relative energetics, and enthalpies of formation of (CH(3)IO(3)) isomers and the reaction CH(3)O(2) + IO have been investigated using quantum mechanical methods. Optimization has been performed at the MP2 level of theory, using all electron and effective core potential, ECP, computational techniques. The relative energetics has been studied by single-point calculations at the CCSD(T) level. Methyl iodate, CH(3)OIO(2), is found to be the lowest-energy isomer showing particular stabilization. The two nascent association minima, CH(3)OOOI and CH(3)OOIO, show similar stabilities, and they are considerably higher located than CH(3)OIO(2). Interisomerization barriers have been determined, along with the transition states involved in various pathways of the reaction CH(3)O(2) + IO.

Quantum mechanical investigation of the atmospheric reaction CH3O2 + NO.

The important stationary points on the potential energy surface of the reaction CH(3)O(2) + NO have been investigated using ab initio and density functional theory techniques. The optimizations were carried out at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) levels of theory while the energetics have been refined using the G2MP2, G3//B3LYP, and CCSD(T) methodologies. The calculations allow the proper characterization of the transition state barriers that determine the fate of the nascent association conformeric minima of methyl peroxynitrite. The main products, CH(3)O + NO(2), are formed through either rearrangement of the trans-conformer to methyl nitrate and its subsequent dissociation or via the breaking of the peroxy bond of the cis-conformer to CH(3)O + NO(2) radical pair. The important consequences of the proposed mechanism are (a) the allowance on energetic grounds for nitrate formation parallel to radical propagation under favorable external conditions and (b) the confirmation of the conformational preference of the homolytic cleavage of the peroxy bond, discussed in previous literature.

Hydrolysis by phospholipase D of phospholipids in solution state or adsorbed on a silica matrix.

Phospholipases D (PLD) catalyse hydrolysis and transphosphatidylation reactions in phospholipids. In the present study, the hydrolytic activity for cabbage PLD was investigated with five different substrates (dipalmitoylphosphatidylethanolamine (DPPE), dipalmitoylphosphatidylcholine (DPPC), didecanoylphosphatidylcholine (DDPC), 1-alkyl-2-acetyl-sn-glycero-3-phosphocholine and lyso-phosphatidylcholine (lyso-PC)) in solution or adsorbed on a silica matrix. In the specific buffer solutions, where the substrates were proved to form large multilamellar polydisperse aggregates, PLD showed preference for DPPC > DPPE > DDPC > 1-alkyl-2-acetyl-sn-glycero-3-phosphocholine > lyso-PC. When the substrates were adsorbed on the silica matrix, PLD hydrolysed 1-alkyl-2-acetyl-sn-glycero-3-phosphocholine and lyso-PC, DDPC, but not DPPC or DPPE. Theoretical studies of the simplest possible adducts between the phospholipids and the silica matrix were performed. Examination of local geometries of DPPC showed a significant blocking of the P-O-X bond-prone to hydrolysis, which could possibly block the access of PLD. Immobilization of phospholipids could be applied for improving the yield of reactions catalysed by PLD as well as for performing a targeted production of short-chain length phosphatidic acid analogs.

Synthesis, structural characterization, and computational studies of novel diiodine adducts with the heterocyclic thioamides N-methylbenzothiazole-2-thione and benzimidazole-2-thione: implications with the mechanism of action of antithyroid drugs.

Reaction of N-methylbenzothiazole-2-thione (C8H7NS2 or NMBZT) with diiodine produced the charge-transfer (ct) complex [(NMBZT).I2] (1). NMBZT reacts with diiodine in the presence of FeCl3 in a molar ratio of 3:6:1 and forms the ionic complex [[(NMBZT)2I+].[FeCl4]-] (2) together with [[(NMBZT)2I+].[I7]-] (2a) iodonium salt. The reaction of benzimidazole-2-thione (C7H6N2S or MBZIM) with diiodine on the other hand results in the formation of the ct [[(MBZIM)2I]+[I3]-].[(MBZIM).I2] (3) compound. The compounds have been characterized by elemental analyses, DTA-TG, FT-Raman, FT-IR, UV-vis, and 1H NMR spectroscopies, and X-ray crystal structure determinations. Compound 1, C8H7I2NS2, is orthorhombic with a space group Pna2(1) and a = 12.5147(13) angstroms, b = 22.536(3) angstroms, c = 4.2994(5) angstroms, and Z = 4. Compound 2, C16H14Cl4FeIN2S4, is monoclinic, space group C2/c, a = 35.781(2) angstroms, b = 7.4761(5) angstroms, c = 18.4677(12) angstroms, beta = 107.219(1) degrees, and Z = 8. Compound 3, C21H18I6N6S3, monoclinic, space group P2(1)/n, a = 14.0652(11) angstroms, b = 22.536(3) angstroms, c = 4.2994(5) angstroms, beta = 99.635(7) degrees, and Z = 4, consists of two component moieties cocrystallized, one neutral which contains the benzimidazole-2-thione (MBIZM) ligand bonded with an iodine atom through sulfur, forming a compound with a "spoke" structure [(MBZIM)I2] 3a, while the other is the ionic complex [[(MBZII)2I+].[I3]-] (3b). The X-ray crystal structure of 1 shows a bond between the thione-sulfur atom and one of the iodine atoms in an essentially planar arrangement. In the cation of 2, an iodine is coordinated by two thione-sulfur atoms in a linear arrangement but the molecule is not planar. For the first time in the solid state a spoke-ionic mixed complex has been characterized in 3. One component of the structure is a molecular diiodine adduct, i.e., [(MBZIM)I2] (3a), with a linear coordination geometry in a decidedly planar arrangement, and the other component is an ionic adduct [[(MBZIM)2I]+.[I3]-] (3b) with the cation having an arrangement similar to that found for 1. Theoretical calculations using density functional (DFT) and ab initio Hartree-Fock theory have been carried out for 1 and 3a,b. The results are consistent with the experimental data. Conclusions on the behavior of a thioamide, when used as an antithyroid drug, have also been made.

Computational studies of (HIO3) isomers and the HO2 + IO reaction pathways.

The geometries, harmonic vibrational frequencies, relative energetics, and enthalpies of formation of (HIO3) isomers have been examined using quantum mechanical methods. At all levels of theory employed, MP2, B3LYP, and CCSD(T), the lowest energy structure is found to be the HOIO2 form, which shows particular stability. The two isomers HOOOI and HOOIO are closely located at the CCSD(T) level of theory. The higher energy structure is HIO3. The interisomerization transition states have been determined, along with the transition states involved in the various pathways of the reaction HO2 + IO.