Conformational analysis of 4,4‐dimethyl‐4‐silathiane and its S‐oxides

4,4‐Dimethyl‐4‐silathiane and its S‐oxides [n = 0 ( 1 ), 1 ( 2 ), 2 ( 3 )] were studied experimentally by variable temperature dynamic NMR spectroscopy down to 103 K and the frozen ring inversion was revealed for all three compounds. The barriers for the degenerate ring inversion in 1 and 3 were measured to be 4.8 and 5.0 kcal/mol at the coalescence temperatures of 111 and 116 K, respectively, and practically coincide with the calculated barriers of 4.60 kcal/mol in 1 and 4.46 kcal/mol in 3 . The frozen equilibrium mixture 2‐ax/2‐eq contains 37% of the 2‐ax and 63% of the 2‐eq conformer. The ring inversion barrier proved to be ca. 4.8 kcal/mol. Calculations at the B3LYP/6‐311+G(d,p) level of theory showed the 2‐ax conformer to be 0.90 kcal/mol more stable than the 2‐eq conformer in the gas phase whereas in solution the relative stability of the conformers calculated using the PCM model at the same level of theory is inverted to become 0.19 (in CHCl₃) or 0.36 kcal/mol (in DMSO) in favor of the 2‐eq conformer. The chair–chair interconversion mechanism of sulfoxide 2 includes two intermediate energetically equivalent 1,4‐twist forms and the 2,5‐boat transition state: 2‐ax (chair) ? 2 (1,4‐twist) ? [ 2 (2,5‐boat)]^≠ ? 2 (1,4‐twist) ? 2‐eq (chair). The calculated ring inversion barriers are 5.1 ( 2‐ax → 2‐eq ) and 4.2 kcal/mol ( 2‐eq → 2‐ax ) in the gas phase, and 4.03 and 4.22 kcal/mol, respectively, in chloroform. Copyright © 2011 John Wiley & Sons, Ltd.     

Conformational analysis of 4,4‐dimethyl‐1‐(trifluoromethylsulfonyl)‐1,4‐azasilinane and 2,2,6,6‐tetramethyl‐4‐(trifluoromethylsulfonyl)‐1,4,2,6‐oxazadisilinane

4,4‐Dimethyl‐1‐(trifluoromethylsulfonyl)‐1,4‐azasilinane 1 and 2,2,6,6‐tetramethyl‐4‐(trifluoromethylsulfonyl)‐1,4,2,6‐oxazadisilinane 2 were studied by variable temperature dynamic ¹H, ¹³C, ¹⁹F NMR spectroscopy and theoretical calculations at the DFT (density functional theory) and MP2 (Møller‐Plesset 2) levels of theory. Both kinetic (barriers to ring inversion) and thermodynamic data (frozen conformational equilibria) could be obtained for the two compounds. The computations revealed two minima on the potential energy surface for molecules 1 and 2 corresponding to the rotamers with the CF₃SO₂ group directed ‘inward’ and ‘outward’ the ring, the latter being 0.2–0.4 kcal/mol (for 1 ) and 1.1 kcal/mol (for 2 ) more stable than the former. The vibrational calculations at the DFT and MP2 levels of theory give the values of the free energy difference ΔG^o for the ‘inward’ ‘outward’ equilibrium consistent with those determined from the experimentally measured ratio of the rotamers. The structure of crystalline compound 2 was ascertained by X‐ray diffraction analysis. Copyright © 2011 John Wiley & Sons, Ltd.     

Proton transfer reactions of hydrazine‐boranes

Hydrazine‐borane and hydrazine‐diborane contain, respectively, 15.4 and 16.9 wt% of hydrogen and are potential materials for hydrogen storage. In this work we present the gas‐phase complexation energies, acidities, and basicities of hydrazine‐borane and hydrazine‐bisborane calculated at MP2/6‐311 + G(d,p) level. We also report the release of dihydrogen from both protonated complexes (ΔG_{hydrazine‐borane} = ?20.9 kcal/mol and ΔG_{hydrazine‐bisborane} = ?27.2 kcal/mol) which is much more exergonic than from analogues amine‐boranes. The addition of the first BH₃ to the hydrazine releases 17.1 kcal/mol, and the second addition releases 15.8 kcal/mol. The attachment of BH₃ also increases the N―H acidity of hydrazine by 46.3 kcal/mol. It was found that the B―H deprotonation leads to intramolecular rearrangement. The basicity values for hydrazine‐borane and ‐bisborane are 180 and 172.8 kcal/mol, respectively. For both complexes the protonation centres are located at the boron moiety. The protonated structure of hydrazine‐bisborane is cyclic and can be described as H₂ captured between a negatively charged B―H hydrogen and positive boron (B―H??H2??B). Atoms in molecules analysis are used to investigate bond paths in concerning structures. Copyright © 2015 John Wiley & Sons, Ltd.     

Chemical reactivity of lithium‐doped fullerenes

The addition of free radicals and the 1,3 dipolar cycloaddition onto pristine and lithium‐doped C₆₀ were studied by means of the Perdew–Burke–Ernzerhof (PBE) and M06‐2X density functionals. In all cases, lithium increased the reactivity even though for the 1,3 dipolar cycloaddition onto C₆₀ the change observed with respect to bare C₆₀ was minimal. Both functionals employed gave similar encapsulation energies for Li@C₆₀ namely, 33.1 and 38.2 kcal/mol at the PBE/6‐31G* and M06‐2X/6‐31G*, respectively. However, the increased reactivity because of lithium doping determined at the PBE level is smaller as compared with that computed with the M06‐2X functional, whereas that determined at the second‐order Møller–Plesset (MP2) level is the largest one. For example, using the M06‐2X functional the binding energy of fluorine to Li@C₆₀ is 28.5 kcal/mol larger than that determined for C₆₀, whereas at the PBE/6‐31G* level it is predicted to be increased by 24.7 kcal/mol. The results clearly suggest that Li@C₆₀ is a much better free radical scavenger than C₆₀. Finally, the complex hindered rotations of lithium inside C₆₀ are expected to be strongly inhibited because lithium doping increases the well depth between the cage center and the equilibrium position near the addition site of the lithium atom. Copyright © 2011 John Wiley & Sons, Ltd.     

Conformational analysis of 3,3‐dimethyl‐3‐silathiane, 2,3,3‐trimethyl‐3‐silathiane and 2‐trimethylsilyl‐3,3‐dimethyl‐3‐silathiane—preferred conformers,barriers to ring inversion and substituent effects

The first conformational analysis of 3‐silathiane and its C‐substituted derivatives, namely, 3,3‐dimethyl‐3‐silathiane 1 , 2,3,3‐trimethyl‐3‐silathiane 2 , and 2‐trimethylsilyl‐3,3‐dimethyl‐3‐silathiane 3 was performed by using dynamic NMR spectroscopy and B3LYP/6‐311G(d,p) quantum chemical calculations. From coalescence temperatures, ring inversion barriers ΔG^≠ for 1 and 2 were estimated to be 6.3 and 6.8 kcal/mol, respectively. These values are considerably lower than that of thiacyclohexane (9.4 kcal/mol) but slightly higher than the one of 1,1‐dimethylsilacyclohexane (5.5 kcal/mol). The conformational free energy for the methyl group in 2 (?ΔG° = 0.35 kcal/mol) derived from low‐temperature ¹³C NMR data is fairly consistent with the calculated value. For compound 2 , theoretical calculations give ΔE value close to zero for the equilibrium between the 2 ‐Me_ax and 2 ‐Me_eq conformers. The calculated equatorial preference of the trimethylsilyl group in 3 is much more pronounced (?ΔG° = 1.8 kcal/mol) and the predominance of the 3 ‐SiMe_{3 eq} conformer at room temperature was confirmed by the simulated ¹H NMR and 2D NOESY spectra. The effect of the 2‐substituent on the structural parameters of 2 and 3 is discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

Theoretical calculations for neighboring group participation in gas‐phase elimination kinetics of 2‐hydroxyphenethyl chloride and 2‐methoxyphenethyl chloride

Theoretical calculation of the kinetics and mechanisms of gas‐phase elimination of 2‐hydroxyphenethyl chloride and 2‐methoxyphenethyl chloride has been carried out at the MP2/6‐31G(d,p), B3LYP/6‐31G(d,p), B3LYP/6‐31 + G(d,p), B3PW91/6‐31G(d,p) and CCSD(T) levels of the theory. The two substrates undergo parallel elimination reactions. The first process of elimination appears to proceed through a three‐membered cyclic transition state by the anchimeric assistance of the aromatic ring to produce the corresponding styrene product and HCl. The second process of elimination occurs through a five‐membered cyclic transition state by participation of the oxygen of o‐OH or the o‐OCH₃ to yield in both cases benzohydrofuran. The B3PW91/6‐31G(d,p) method was found to be in good agreement with the experimental kinetic and thermodynamic parameters for both substrates in the two reaction channels. However, some differences in the performance of the different methods are observed. NBO analysis of the pyrolysis of both phenethyl chlorides implies a C? Cl bond polarization, in the sense of C^δ+…Cl^δ?, which is a rate‐determining step for both parallel reactions. Synchronicity parameters imply polar transition states of these elimination reactions. Copyright © 2009 John Wiley & Sons, Ltd.     

P‐Heterocyclic silylenes: a survey of stability with density functional theory

The effects of phosphorous atom on the stability, multiplicity, and reactivity of six‐member cyclic silylenes are investigated at B3LYP/AUG‐cc‐pVTZ//B3LYP/6‐31+G* and MP2/6‐311++G**//B3LYP/6‐31+G* coupled with appropriate isodesmic reactions. From a thermodynamic point of view, 1H‐2‐silaphosphinine‐2‐ylidene ( 1_a ) and 1H‐4‐silaphosphinine‐4‐ylidene ( 2_a ) are relatively the most stable with singlet–triplet energy gaps (ΔE_S–T) of 37.0 and 28.1 kcal/mol, respectively. The calculated energy barrier for the 1,2‐H shift of 1_a to the corresponding 2‐silapyridine ( 1 ) is 26.5 kcal/mol, which is lower than the 28.8 kcal/mol required for the 1,4‐H shift of 2_a to the corresponding 4‐silapyridine ( 2 ). In contrast to the previous reports, isodesmic reactions indicate that π‐donor/σ‐donor phosphorous destabilizes the singlet while stabilizes the triplet state. Both 1_a and 2_a silylenes appear invulnerable to the head‐to‐head as well as the head‐to‐tail dimerization, inviting experimental explorations. Copyright © 2011 John Wiley & Sons, Ltd.     

Joint theoretical and experimental study of the gas‐phase elimination kinetics of tert‐butyl ester of carbamic,N, N‐dimethylcarbamic,N‐hydroxycarbamic acids and 1‐(tert‐butoxycarbonyl)‐imidazole

The gas‐phase elimination kinetics of the title compounds were carried out in a static reaction system and seasoned with allyl bromide. The working temperature and pressure ranges were 200–280 °C and 22–201.5 Torr, respectively. The reactions are homogeneous, unimolecular, and follow a first‐order rate law. These substrates produce isobutene and corresponding carbamic acid in the rate‐determining step. The unstable carbamic acid intermediate rapidly decarboxylates through a four‐membered cyclic transition state (TS) to give the corresponding organic nitrogen compound. The temperature dependence of the rate coefficients is expressed by the following Arrhenius equations: for tert‐butyl carbamate logk₁ (s^?1) = (13.02 ± 0.46) – (161.6 ± 4.7) kJ/mol(2.303 RT)^?1, for tert‐butyl N‐hydroxycarbamate logk₁ (s^?1) = (12.52 ± 0.11) – (147.8 ± 1.1) kJ/mol(2.303 RT)^?1, and for 1‐(tert‐butoxycarbonyl)‐imidazole logk₁ (s^?1) = (11.63 ± 0.21)–(134.9 ± 2.0) kJ/mol(2.303 RT)^?1. Theoretical studies of these elimination were performed at Møller–Plesset MP2/6‐31G and DFT B3LYP/6‐31G(d), B3LYP/6‐31G(d,p) levels of theory. The calculated bond orders, NBO charges, and synchronicity (Sy) indicate that these reactions are concerted, slightly asynchronous, and proceed through a six‐membered cyclic TS type. Results for estimated kinetic and thermodynamic parameters are discussed in terms of the proposed reaction mechanism and TS structure. Copyright © 2007 John Wiley & Sons, Ltd.     

The pKa theoretical estimation of C―H,N―H,O―H and S―H acids in dimethylsulfoxide solution

A scheme for the pK_a estimation of organic acids in dimethylsulfoxide (DMSO) solution based on quantum chemical calculations is proposed. The procedure of pK_a calculation requires several steps. The first is the calculation of the gas phase acidity of the compound. The G3MP2B3, G4MP2 as well as CBS‐QB3 composite methods made it possible to estimate values of gas phase acidities of an extensive set of structures with a high confidence level (standard deviations equal to 1.15, 1.13 and 1.29 kcal mol^?1, respectively; the test set included 91 compounds). The second step is the computation of the solvation correction with the integral equation formalism version of polarizable continuum model (IEF‐PCM)–B3LYP/6‐311+G(d,p) approximation. Within the bounds of our approach, the medium properties were covered only by the PCM model, i.e. the proposed procedure neglects specific interactions between DMSO and the solute. It was determined that the approach to pK_a estimation mentioned above is the most balanced in terms of accuracy, resource intensity and computation time cost. In the third step, the error of the pK_a calculation was decreased by correlation allowances. Correlation allowances were determined for each acid class (62 С―Н, 55 N―Н, 24 O―Н and 5 S―Н acids) in the range of 50 units in terms of logarithmic scale using the test set including 146 compounds. Seven O―H acids showing the ability to form cyclic dimers were separated into a discrete group. The proposed methodology was applied to the estimation of pK_a for trans‐ and cis‐dimethyl‐4,5‐dihydro‐3H‐pyrazol‐3,5‐dicarboxylates as well as for 5‐fluorouracil subject to competitive dissociation, the latter by N1―H or N3―H bonds. Copyright © 2014 John Wiley & Sons, Ltd.     

Insights into unusual stability of 5‐membered‐ring endocyclic benzyl carbocations in aqueous solution

The acyclic o‐oxygen benzyl carbocation 1 , the 6‐membered‐ring endocyclic o‐oxygen benzyl carbocation 2 , and the 5‐membered‐ring endocyclic o‐oxygen benzyl carbocation 3 were used as model compounds to get insights into the general phenomenon for the unusual stability of the 5‐membered‐ring endocyclic benzyl carbocations in aqueous solution. The hydride‐ion affinities of 1 , 2 , and 3 in gas phase, acetonitrile, and DMSO were calculated and compared by the density functional theory method, and 3 isodesmic reactions were designed to confirm their thermodynamic stability. What we found is that the 5‐membered‐ring endocyclic o‐oxygen stabilizes the benzyl carbocation 3 less than the acyclic o‐oxygen stabilizes the benzyl carbocation 1 in gas phase because of ring strain and through‐bond induction. However, the high solvation energies of the 5‐membered‐ring endocyclic o‐oxygen benzyl carbocation 3 not only offset the destabilizing effects of ring strain and through‐bond induction but also make it even more stable than the acyclic o‐oxygen benzyl carbocation 1 in polar solvents like acetonitrile, DMSO, and water.     

Theoretical calculations of the kinetics and mechanisms of the gas phase elimination of primary,secondary, and tertiary 2‐hydroxyalkylbenzenes

Theoretical study of the elimination kinetics of 2‐phenylethanol, 1‐phenyl‐2‐propanol, and 2‐methyl‐1‐phenyl‐2‐propanol in the gas‐phase has been carried out at the MP2/6‐31G(d,p), B3LYP/6‐31G(d,p), B3LYP/6‐31++G(d,p), MPW1PW91/6‐31G(d,p), MPW1PW91/6‐31++G(d,p), PBEPBE/6‐31G(d,p), and PBEPBE/6‐31++G(d,p) levels of theory. The three substrates undergo two parallel elimination reactions. The first elimination appears to proceed through a six‐membered cyclic transition state to give toluene and the corresponding aldehyde or ketone. The second parallel elimination takes place through a four‐membered cyclic transition state producing water and the corresponding unsaturated aromatic hydrocarbon. Results from MP2/6‐31G(d,p) and MPW1PW91/6‐31++G(d,p) methods were found to be in good agreement with the experimental kinetic and thermodynamic parameters in the formation of toluene and the corresponding carbonyl compound. However, the results for PBEPBE/6‐31G(d,p) were in better agreement with the experimental data for the second parallel reaction yielding water and the corresponding unsaturated aromatic hydrocarbon. The charge distribution differences in the TS related to the substitution by methyl groups in the substrates can account for the observed reaction rate coefficients. The synchronicity parameters imply semi‐polar transition states for these elimination reactions. Copyright © 2009 John Wiley & Sons, Ltd.     

Structure and NMR properties of 6‐substituted‐5,6‐dihydrobenzo[c]phenanthridine alkaloids

We report a preparation of new 6‐substituted‐5,6‐dihydrobenzo[c]phenanthridines by the reaction of azoles with quaternary benzo[c]phenanthridine alkaloids sanguinarine and chelerythrine. The prepared compounds have been characterized by NMR spectroscopy, mass spectrometry, and single‐crystal X‐ray diffraction. Conformational behaviors of carbazole derivatives in solution have been investigated by low‐temperature NMR experiments. Barriers to rotation around newly formed C6–N bonds were determined to be 12–13 kcal/mol. Quantum chemical calculations have been used to reproduce the experimental observations. Large structural effects on several ¹H NMR resonances were observed experimentally, analyzed by Density Functional Theory (DFT) calculations at B3LYP/6‐311+G(d,p)/PCM level, and interpreted by ring‐current effects of the benzo[c]phenanthridine and carbazole units. Copyright © 2013 John Wiley & Sons, Ltd.     

CO2 production in the bromate‐1,4‐cyclohexanedione oscillatory reaction

NMR and GC/MS spectroscopy of the organic extracts of the oscillatory bromate‐1,4‐cyclohexanedione reaction illustrate the presence of ring‐opening products 5‐(dibromomethylene)‐2(5H)‐furanone, (E)‐5,5,5‐tribromo‐4‐oxo‐2‐pentenoic acid, and dibromoacetic acid, particularly at elevated temperatures. The loss of a carbon atom from the six‐membered ring after ring opening led to gas formation and such a process became more vigorous at >60 °C, with the direct observation of bubbles in a stirred batch reactor. Gravimetric experiments confirm that the amount of carbon dioxide gas produced increases rapidly with reaction temperature. Parallel experiments suggest that the ring‐opening process involves the oxidation of brominated benzoquinones by bromate. Copyright © 2010 John Wiley & Sons, Ltd.     

Homoaromaticity in aza‐ and phosphasemibullvalenes. A computational study

The enthalpies of activation ΔH^? for the Cope rearrangement in several aza‐ and phosphasemibullvalenes have been investigated by MP4/cc‐pVDZ//MP2/cc‐pVDZ and CCSD(T)/cc‐pVDZ//MP2/cc‐pVDZ calculations. One tetraazasemibullvalene and several phosphasemibullvalenes were found to have vanishing ΔH^? values, which together with calculated large negative Nucleus Independent Chemical Shift (NICS) values and geometrical data show that these molecules have delocalized and bishomoaromatic minima. Furthermore, three azasemibullvalenes were found to have small ΔH^? values (≤2 kcal/mol) combined with large negative NICS, suggesting that they could also have bishomoaromatic minima. Copyright © 2010 John Wiley & Sons, Ltd.     

4‐Alkyl‐2,2,6,6‐tetramethyl‐1,4,2,6‐oxaazadisilinanes: synthesis,structure, and conformational analysis

4‐Alkyl‐2,2,6,6‐tetramethyl‐1,4,2,6‐oxaazadisilinanes RN[CH₂Si(Me)₂]₂O [R = Me ( 1 ), i‐Pr ( 2 )] were synthesized by two methods which provided good yields up to 84%. Low temperature NMR study of compounds ( 1 ) and ( 2 ) revealed a frozen ring inversion with the energy barriers of 8.5 and 7.7 kcal/mol at 163 and 143 K, respectively, which is substantially lower than that for their carbon analog, N‐methylmorpholine. DFT calculations performed on the example of molecule ( 1 ) showed that N? Me_ax conformer to exist in the sofa conformation with the coplanar fragment C? Si? O? Si? C, and its N? Me_eq conformer in a flattened chair conformation. Copyright © 2009 John Wiley & Sons, Ltd.     

The keto-enol equilibrium in substituted acetaldehydes: focal-point analysis and ab initio limit

High-level ab initio electronic structure calculations up to the CCSD(T) theory level, including extrapolations to the complete basis set (CBS) limit, resulted in high precision energetics of the tautomeric equilibrium in 2-substituted acetaldehydes (XH₂C-CHO). The CCSD(T)/CBS relative energies of the tautomers were estimated using CCSD(T)/aug-cc-pVTZ, MP3/aug-cc-pVQZ, and MP2/aug-cc-pV5Z calculations with MP2/aug-cc-pVTZ geometries. The relative enol (XHC?=?CHOH) stabilities (ΔE _{e,CCSD(T)/CBS}) were found to be 5.98?±?0.17, ?1.67?±?0.82, 7.64?±?0.21, 8.39?±?0.31, 2.82?±?0.52, 10.27?±?0.39, 9.12?±?0.18, 5.47?±?0.53, 7.50?±?0.43, 10.12?±?0.51, 8.49?±?0.33, and 6.19?±?0.18?kcal?mol^?1 for X?=?BeH, BH₂, CH₃, Cl, CN, F, H, NC, NH₂, OCH₃, OH, and SH, respectively. Inconsistencies between the results of complex/composite energy computations methods Gn/CBS (G2, G3, CBS-4M, and CBS-QB3) and high-level ab initio methods (CCSD(T)/CBS and MP2/CBS) were found. DFT/aug-cc-pVTZ results with B3LYP, PBE0 (PBE1PBE), TPSS, and BMK density functionals were close to the CCSD(T)/CBS levels (MAD?=?1.04?kcal?mol^?1).     

Theoretical prediction of size‐expansion effect on the C8‐site activity in the modified guanine‐cytosine analogs

Homo/hetero ring‐expanded DNA analogs have been shown to be rationally modified DNA motifs with improved physical or biological properties. In this work, using density functional theory, the stability of these artificial DNA base pairs was examined with regard to three aspects associated with DNA damage, namely deprotonation, H‐abstraction, and H‐radical addition. The effect of size expansion on C8 activity was investigated because C8‐oxidative guanine (G) is one of the most important products of DNA damage. Computational results indicate that the insertion of an aromatic spacer ring in G considerably decreases the electron density over the C8 site, leading to easier deprotonation or H‐abstraction from the C8 site and more difficult H ^. ‐radical attack on the C8 site. However, the opposite phenomenon is observed if the spacer ring is antiaromatic, because of the increased electron density over the C8 site. Moreover, these effects are more prominent the larger the aromaticity or antiaromaticity of the spacer ring. Further analyses, using natural bond orbitals (NBOs) and the nucleus‐independent chemical shift (NICS) index of aromaticity, indicate that the changes of the electron distribution over the C8 site arise because the aromatic spacer ring, involved in the conjugation structure, increases the electron delocalization from the electron‐rich imidazole ring to the diatropic six‐membered rings, while the antiaromatic spacer ring acts as an electron‐donating group, not only inhibiting the above electron delocalization, but also slightly increasing the electron density over the C8 site. The improved stability of these size‐expanded base pairs in different DNA‐damaged environments may encourage their use in practical applications. Copyright © 2009 John Wiley & Sons, Ltd.     

DFT study on abstraction reaction mechanism of oh radical with 2‐methoxyphenol

Reaction mechanism of 2‐methoxyphenol (2MP) (guaiacol) with OH radical has been performed using density functional theory methods BH&HLYP and MPW1K method with 6‐311++G(d,p) basis set. Single‐point energy calculations were done using CCSD(T)/6‐311++G(d,p). The theoretical results reveal that the hydrogen abstraction from methoxy group is found to be the dominant reaction channel with an energy barrier of 9.31 kcal/mol. Also, time‐dependent density functional theory calculations have been performed using BH&HLYP/6‐311++G(d,p) level of theory, and the results reveal that the reactions occur in ground state than the excited state. The results of reaction force profile indicate that structural rearrangements are most influential with high percentage than the relaxation process. The calculated theoretical rate constants (12.19 × 10^?11 cm³ molecule^?1 s^?1) are in good agreement with the experimental rate constant. The atmospheric lifetime of 2‐methoxyphenol with respect to OH radicals is 2.27 hours, which implies that OH radical plays an important role in the degradation of 2MP. The Wiberg bond index of the abstraction reaction reveals that the bond order is concerted, partially synchronic. The reactant‐like transition state satisfies Hammond postulate, which eventually results in an exothermic reaction, and the product‐like transition state reveals in endothermic nature.     

Thermochemistry,bond energies and internal rotor barriers of methyl sulfinic acid,methyl sulfinic acid ester and their radicals

Sulfur–Oxygen containing hydrocarbons are formed in oxidation of sulfides and thiols in the atmosphere, on aerosols and in combustion processes. Understanding their thermochemical properties is important to evaluate their formation and transformation paths. Structures, thermochemical properties, bond energies, and internal rotor potentials of methyl sulfinic acid CH₃S(?O)OH, its methyl ester CH₃S(?O)OCH₃ and radicals corresponding to loss of a hydrogen atom have been studied. Gas phase standard enthalpies of formation and bond energies were calculated using B3LYP/6‐311G (2d, p) individual and CBS‐QB3 composite methods employing work reactions to further improve accuracy of the ${\Delta} _{{\bf f}} H_{{\bf 298}}^{{\bf o}} $ . Molecular structures, vibration frequencies, and internal rotor potentials were calculated. Enthalpies of the parent molecules CH₃S(?O)OH and CH₃S(?O)OCH₃ are evaluated as ?77.4 and ?72.7 kcal mol^?1 at the CBS? QB3 level; Enthalpies of radicals C^?H₂? S(?O)? OH, CH₃? S^?(?O)₂, C^?H₂? S(?O)? OCH₃ and CH₃? S(?O)? OC^?H₂ (CBS‐QB3) are ?25.7, ?52.3, ?22.8, and ?26.8 kcal mol^?1, respectively. The CH₃C(?O)O—H bond dissociation energy is of 77.1 kcal mol^?1. Two of the intermediate radicals are unstable and rapidly dissociate. The CH₃S(?O)? O. radical obtained from the parent CH₃? S(?O)? OH dissociates into methyl radical (${\bf CH}_{{\bf 3}}^{{\bf .}} $ ) plus SO₂ with endothermicity (ΔH_rxn) of only 16.2 kcal mol^?1. The CH₃? S(?O)? OC^?H₂ radical dissociates into CH₃? S^?=O and CH₂=O with little or no barrier and an exothermicity of ?19.9 kcal mol^?1. DFT and the Complete Basis Set‐QB3 enthalpy values are in close agreement; this accord is attributed to use of isodesmic work reactions for the analysis and suggests this combination of B3LYP/work reaction approach is acceptable for larger molecules. Copyright © 2010 John Wiley & Sons, Ltd.     

Computational study on structures,thermochemical properties,and bond energies of disulfide oxygen (S–S–O)‐bridged CH3SSOH and CH3SS(=O)H and radicals

Cleavage of disulfide bonds is a common method used in linking peptides to proteins in biochemical reactions. The structures, internal rotor potentials, bond energies, and thermochemical properties (Δ_fH°, S°, and Cp(T)) of the S–S bridge molecules CH₃SSOH and CH₃SS(=O)H and the radicals CH₃SS?=O and C?H₂SSOH that correspond to H‐atom loss are determined by computational chemistry. Structure and thermochemical parameters (S° and Cp(T)) are determined using density functional Becke, three‐parameter, Lee–Yang–Parr (B3LYP)/6‐31++G (d, p), B3LYP/6‐311++G (3df, 2p). The enthalpies of formation for stable species are calculated using the total energies at B3LYP/6‐31++G (d, p), B3LYP/6‐311++G (3df, 2p), and the higher level composite CBS–QB3 levels with work reactions that are close to isodesmic in most cases. The enthalpies of formation for CH₃SSOH, CH₃SS(=O)H are ?38.3 and ?16.6 kcal mol^?1, respectively, where the difference is in enthalpy RSO–H versus RS(=O)–H bonding. The C–H bond energy of CH₃SSOH is 99.2 kcal mol^?1, and the O–H bond energy is weaker at 76.9 kcal mol^?1. Cleavage of the weak O–H bond in CH₃SSOH results in an electron rearrangement upon loss of the CH₃SSO–H hydrogen atom; the radical rearranges to form the more stable CH₃SS· = O radical structure. Cleavage of the C–H bond in CH₃SS(=O)H results in an unstable [CH₂SS(=O)H]* intermediate, which decomposes exothermically to lower energy CH₂ = S + HSO. The CH₃SS(=O)–H bond energy is quite weak at 54.8 kcal mol^?1 with the H–C bond estimated at between 91 and 98 kcal mol^?1. Disulfide bond energies for CH₃S–SOH and CH3S–S(=O)H are low: 67.1 and 39.2 kcal mol^?1. Copyright © 2011 John Wiley & Sons, Ltd.