Kinetics of hydrogen abstraction reactions of butene isomers by OH radical

The rate coefficients of H-abstraction reactions of butene isomers by the OH radical were determined by both canonical variational transition-state theory and transition-state theory, with potential energy surfaces calculated at the CCSD(T)/6-311++G(d,p)//BH&HLYP/6-311G(d,p) level and CCSD(T)/6-311++G(d,p)//BH&HLYP/cc-pVTZ level and quantum mechanical tunneling effect corrected by either the small-curvature tunneling method or the Eckart method. While 1-butene contains allylic, vinylic, and alkyl hydrogens that can be abstracted to form different butene radicals, results reveal that s-allylic H-abstraction channels have low and broad energy barriers, and they are the most dominant channels which can occur via direct and indirect H-abstraction channels. For the indirect H-abstraction s-allylic channel, the reaction can proceed via forming two van der Waals prereactive complexes with energies that are 2.7-2.8 kcal mol(-1) lower than that of the entrance channel at 0 K. Assuming that neither mixing nor crossover occurs between different reaction pathways, the overall rate coefficient was calculated by summing the rate coefficients of the s-allyic, methyl, and vinyl H-abstraction paths and found to agree well with the experimentally measured OH disappearance rate. Furthermore, the rate coefficients of p-allylic H abstraction of cis-2-butene, trans-2-butene, and isobutene by the OH radical were also determined at 300-1500 K, with results analyzed and compared with available experimental data.     

Thermal decomposition of ethanol. 4. Ab initio chemical kinetics for reactions of H atoms with CH3CH2O and CH3CHOH radicals

The potential energy surfaces of H-atom reactions with CH(3)CH(2)O and CH(3)CHOH, two major radicals in the decomposition and oxidation of ethanol, have been studied at the CCSD(T)/6-311+G(3df,2p) level of theory with geometric optimization carried out at the BH&HLYP/6-311+G(3df,2p) level. The direct hydrogen abstraction channels and the indirect association/decomposition channels from the chemically activated ethanol molecule have been considered for both reactions. The rate constants for both reactions have been calculated at 100-3000 K and 10(-4) Torr to 10(3) atm Ar pressure by microcanonical VTST/RRKM theory with master equation solution for all accessible product channels. The results show that the major product channel of the CH(3)CH(2)O + H reaction is CH(3) + CH(2)OH under atmospheric pressure conditions. Only at high pressure and low temperature, the rate constant for CH(3)CH(2)OH formation by collisonal deactivation becomes dominant. For CH(3)CHOH + H, there are three major product channels; at high temperatures, CH(3)+CH(2)OH production predominates at low pressures (P < 100 Torr), while the formation of CH(3)CH(2)OH by collisional deactivation becomes competitive at high pressures and low temperatures (T < 500 K). At high temperatures, the direct hydrogen abstraction reaction producing CH(2)CHOH + H(2) becomes dominant. Rate constants for all accessible product channels in both systems have been predicted and tabulated for modeling applications. The predicted value for CH(3)CHOH + H at 295 K and 1 Torr pressure agrees closely with available experimental data. For practical modeling applications, the rate constants for the thermal unimolecular decomposition of ethanol giving key accessible products have been predicted; those for the two major product channels taking place by dehydration and C-C breaking agree closely with available literature data.     

A single transition state serves two mechanisms. The branching ratio for CH2O*- + CH3Cl on improved potential energy surfaces

The reaction of formaldehyde radical anion with methyl chloride, CH2O*- + CH3Cl, is an example in which a single transition state leads to two products: substitution at carbon (Sub(C), CH3CH2O* + Cl-) and electron transfer (ET, CH2O + CH3* + Cl-). The branching ratio for this reaction has been studied by ab initio molecular dynamics (AIMD). The energies of transition states and intermediates were computed at a variety of levels of theory and compared to accurate energetics calculated by the G3 and CBS-QB3 methods. A bond additivity correction has been constructed to improve the Hartree-Fock potential energy surface (BAC-UHF). A satisfactory balance between good energetics and affordable AIMD calculations can be achieved with BH&HLYP/6-31G(d) and BAC-UHF/6-31G(d) calculations. Approximately 200 ab initio classical trajectories were calculated for each level of theory with initial conditions sampled from a thermal distribution at 298 K at the transition state. Three types of trajectories were distinguished: trajectories that go directly to ET product, trajectories that go to Sub(C) product, and trajectories that initially go into the Sub(C) valley and then dissociate to ET products. The BH&HLYP/6-31G(d) calculations overestimate the number of nonreactive and direct ET trajectories because the transition state is too early. For the BH&HLYP and BAC-UHF methods, about one-third of the trajectories that initially go into the Sub(C) valley dissociate to ET products, compared to just over half with UHF/6-31G(d) in the earlier study. This difference can be attributed to a better value for the calculated energy release from the initial transition state and to an improved Sub(C) --> ET barrier height with the BH&HLYP and BAC-UHF methods.     

Theoretical study on the kinetics of OH radical reactions with CH3OOH and CH3CH2OOH

The mechanisms and kinetics studies of the OH radical with alkyl hydroperoxides CH(3)OOH and CH(3)CH(2)OOH reactions have been carried out theoretically. The geometries and frequencies of all the stationary points are calculated at the UBHandHLYP/6-311G(d,p) level, and the energy profiles are further refined by interpolated single-point energies method at the MC-QCISD level of theory. For two reactions, five H-abstraction channels are found and five products (CH(3)OO, CH(2)OOH, CH(3)CH(2)OO, CH(2)CH(2)OOH, and CH(3)CHOOH) are produced during the above processes. The rate constants for the CH(3)OOH/CH(3)CH(2)OOH + OH reactions are corrected by canonical variational transition state theory within 250-1500 K, and the small-curvature tunneling is included. The total rate constants are evaluated from the sum of the individual rate constants and the branching ratios are in good agreement with the experimental data. The Arrhenius expressions for the reactions are obtained.     

Computational study of static first hyperpolarizability of donor-acceptor substituted (E)-benzaldehyde phenylhydrazone

Various hybrid functionals (B3LYP, B97-2, PBE0, BMK, BH&HLYP, CAM-B3LYP, and LC-ωPBE) implemented in density functional theory were applied to give estimate of static first hyperpolarizabilty (β(0)) of (E)-benzaldehyde phenylhydrazone designated as (E)-BPH. Against those of MP2 computations as a function of the underlying density functional, good agreement was obtained with the BH&HLYP and CAM-B3LYP functionals. The LC-ωPBE functional and the B3LYP, PBE0, B97-2, and BMK functionals underestimated and overestimated β(0), respectively. The basis set effect on the calculated β(0) was also investigated. It turned out that the 6-311+G(2d,p) basis set provided excellent converged value of β(0). On the basis of the calculated results, we investigated the substituent effect on β(0) of donor-acceptor (D-A) substituted (E)-BPH systematically by using the BH&HLYP and CAM-B3LYP computations with the 6-311+G(2d,p) basis set. We proposed a Zwitterion structure to explain the calculated trend in the substituent effect and the enhanced hyperpolarizability of type II compounds (A-(E)-BPH-D) than type I compounds (D-(E)-BPH-A). Natural bonding orbital analysis carried out at BH&HLYP/6-311+G(2d,p)//B3LYP/6-31G(2df,p) level of theory substantiated the claim.     

Thermal decomposition of methyl butanoate: ab initio study of a biodiesel fuel surrogate

In this paper, we report a detailed analysis of the breakdown kinetic mechanism for methyl butanoate (MB) using theoretical approaches. Electronic structures and structure-related molecular properties of reactants, intermediates, products, and transition states were explored at the BH&HLYP/cc-pVTZ level of theory. Rate constants for the unimolecular and bimolecular reactions in the temperature range of 300-2500 K were calculated using Rice-Ramsperger-Kassel-Marcus and transition state theories, respectively. Thirteen pathways were identified leading to the formation of small compounds such as CH(3), C(2)H(3), CO, CO(2), and H(2)CO. For the initial formation of MB radicals, H, CH(3), and OH were considered as reactive radicals participating in hydrogen abstraction reactions. Kinetic simulation results for a high temperature pyrolysis environment show that MB radicals are mainly produced through hydrogen abstraction reactions by H atoms. In addition, the C(O)OCH(3) = CO + CH(3)O reaction is found to be the main source of CO formation. The newly computed kinetic sub-model for MB breakdown is recommended as a core component to study the combustion of oxygenated species.     

Atmospheric chemistry of ethyl propionate

Ethyl propionate is a model for fatty acid ethyl esters used as first-generation biodiesel. The atmospheric chemistry of ethyl propionate was investigated at 980 mbar total pressure. Relative rate measurements in 980 mbar N(2) at 293 ± 0.5 K were used to determine rate constants of k(C(2)H(5)C(O)OC(2)H(5) + Cl) = (3.11 ± 0.35) × 10(-11), k(CH(3)CHClC(O)OC(2)H(5) + Cl) = (7.43 ± 0.83) × 10(-12), and k(C(2)H(5)C(O)OC(2)H(5) + OH) = (2.14 ± 0.21) × 10(-12) cm(3) molecule(-1) s(-1). At 273-313 K, a negative Arrhenius activation energy of -3 kJ mol(-1) is observed.. The chlorine atom-initiated oxidation of ethyl propionate in 980 mbar N(2) gave the following products (stoichiometric yields): ClCH(2)CH(2)C(O)OC(2)H(5) (0.204 ± 0.031), CH(3)CHClC(O)OC(2)H(5) (0.251 ± 0.040), and C(2)H(5)C(O)OCHClCH(3) (0.481 ± 0.088). The chlorine atom-initiated oxidation of ethyl propionate in 980 mbar of N(2)/O(2) (with and without NO(x)) gave the following products: ethyl pyruvate (CH(3)C(O)C(O)OC(2)H(5)), propionic acid (C(2)H(5)C(O)OH), formaldehyde (HCHO), and, in the presence of NO(x), PAN (CH(3)C(O)OONO(2)). The lack of acetaldehyde as a product suggests that the CH(3)CH(O)C(O)OC(2)H(5) radical favors isomerization over decomposition. From the observed product yields, we conclude that H-abstraction by chlorine atoms from ethyl propionate occurs 20.4 ± 3.1%, 25.1 ± 4.0%, and 48.1 ± 8.8% from the CH(3)-, -CH(2)-, and -OCH(2)- groups, respectively. The rate constant and branching ratios for the reaction between ethyl propionate and the OH radical were investigated theoretically using quantum mechanical calculations and transition state theory. The stationary points along the reaction path were optimized using the CCSD(T)-F12/VDZ-F12//BH&HLYP/aug-cc-pVTZ level of theory; this model showed that OH radicals abstract hydrogen atoms primarily from the -OCH(2)- group (80%).     

Effect of substituents on the stabilities of multiply-substituted carbon-centered radicals

The bond dissociation energies (BDEs) and radical stabilization energies (RSEs) which result from 166 reactions that lead to carbon-centered radicals of the type ˙CH(2)X, ˙CHXY and ˙CXYZ, where X, Y and Z are any of the fourteen substituents H, F, Cl, NH(2), OH, SH, CH[double bond, length as m-dash]CH(2), C[triple bond, length as m-dash]CH, BH(2), CHO, COOH, CN, CH(3), and CF(3), were calculated using spin-restricted and -unrestricted variants of the double-hybrid B2-PLYP method with the 6-311+G(3df,2p) basis set. The interactions of substituents X, Y, and Z in both the radicals (˙CXYZ) and in the precursor closed-shell molecules (CHXYZ), as well as the extent of additivity of such interactions, were investigated by calculating radical interaction energies (RIEs), molecule interaction energies (MIEs), and deviations from additivity of RSEs (DARSEs) for a set of 152 reactions that lead to di- (˙CHXY) and tri- (˙CXYZ) substituted carbon-centered radicals. The pairwise quantities describing the effects of pairs of substituents in trisubstituted systems, namely pairwise MIEs (PMIEs), pairwise RIEs (PRIEs) and deviations from pairwise additivity of RSEs (DPARSEs), were also calculated for the set of 61 reactions that lead to trisubstituted radicals (˙CXYZ). Both ROB2-PLYP and UB2-PLYP were found to perform quite well in predicting the quantities related to the stabilities of carbon-centered radicals when compared with available experimental data and with the results obtained from the high-level composite method G3X(MP2)-RAD. Particular selections of substituents or combinations of substituents from the current test set were found to lead to specially stable radicals, increasing the RSEs to a maximum of +68.2 kJ mol(-1) for monosubstituted radicals ˙CH(2)X (X = CH[double bond, length as m-dash]CH(2)), +131.7 kJ mol(-1) for disubstituted radicals ˙CHXY (X = NH(2), Y = CHO), and +177.1 kJ mol(-1) for trisubstituted radicals ˙CXYZ (X = NH2, Y = Z = CHO).     

Dual-level direct dynamics studies for the reactions of OH radical with bromine-substituted ethanes

The dynamic properties of the multichannel hydrogen abstraction reactions of CH(3)CH(2)Br + OH --> products and CH(3)CHBr(2) + OH --> products are studied by dual-level direct dynamics method. For each reaction, three reaction channels, one for alpha-hydrogen abstraction and two for beta-hydrogen abstractions, have been identified. The minimum energy paths (MEPs) of both the reactions are calculated at the Becke's half-and-half (BH&H)-Lee-Yang-Parr (LYP)/6-311G(d, p) level and the energy profiles along the MEPs are further refined with interpolated single-point energies (ISPE) method at the G2M(RCC5)//BH&H-LYP level. There are complexes with energies less than those of the reactants or products located at the entrance or exit channels, which indicates that the reactions may proceed via an indirect mechanism. By canonical variational transition-state theory (CVT) the rate constants are calculated incorporating the small-curvature tunneling (SCT) correction in the temperature range of 220-2000 K. The agreement of the rate constants with available experimental values for two reactions is good in the measured temperature range. The calculated results show that alpha-hydrogen abstraction channel is the major reaction pathway in the lower temperature for two reactions, while the contribution of beta-hydrogen abstraction will increase with the increase in temperature.     

Computations of 36 tautomer/isomer equilibria of different lactams

Gas-phase energies of 36 tautomer/isomer pairs of 18 six-membered N-heterocyclic compounds were computed quantum chemically. Among the considered B3LYP, BH&HLYP, BH&HLYP(G), and PW6B95 DFT functionals, the latter two provide accurate tautomer/isomer pair energies with root-mean-square deviations (rmsd) relative to experiments of 0.2 and 0.3 kcal/mol, respectively. Since only few (namely five) experimental data are available, 15 tautomer/isomer pair energies were computed with the very precise QCISD(T)(quadruple-ζ) method serving as reference. Relative to this reference the PW6B95 DFT functional is slightly superior to the BH&HLYP(G) functional, yielding an rmsd of 0.7 and 0.8 kcal/mol, respectively. In contrast to BH&HLYP(G), the PW6B95 DFT functional yields also accurate tautomer/isomer pair energies if zwitterionic structures are involved. The tautomer/isomer pair states possess different amounts of aromaticity. This is characterized by nucleus-independent chemical shift (NICS) values. The tautomer/isomer pair reference energies, from which the energies computed with PW6B95 are subtracted, correlate linearly with the corresponding differences in the NICS values. This correlation is used to construct a correction term for the pair energies computed with PW6B95, yielding tautomer/isomer pair energies with rmsd of 0.3 kcal/mol with respect to the more CPU time demanding QCISD(T)(quadruple-ζ) method.     

Assessment of experimental bond dissociation energies using composite ab initio methods and evaluation of the performances of density functional methods in the calculation of bond dissociation energies

Composite ab initio CBS-Q and G3 methods were used to calculate the bond dissociation energies (BDEs) of over 200 compounds listed in CRC Handbook of Chemistry and Physics (2002 ed.). It was found that these two methods agree with each other excellently in the calculation of BDEs, and they can predict BDEs within 10 kJ/mol of the experimental values. Using these two methods, it was found that among the examined compounds 161 experimental BDEs are valid because the standard deviation between the experimental and theoretical values for them is only 8.6 kJ/mol. Nevertheless, 40 BDEs listed in the Handbook may be highly inaccurate as the experimental and theoretical values for them differ by over 20 kJ/mol. Furthermore, 11 BDEs listed in the Handbook may be seriously flawed as the experimental and theoretical values for them differ by over 40 kJ/mol. Using the 161 cautiously validated experimental BDEs, we then assessed the performances of the standard density functional (DFT) methods including B3LYP, B3P86, B3PW91, and BH&HLYP in the calculation of BDEs. It was found that the BH&HLYP method performed poorly for the BDE calculations. B3LYP, B3P86, and B3PW91, however, performed reasonably well for the calculation of BDEs with standard deviations of about 12.1-18.0 kJ/mol. Nonetheless, all the DFT methods underestimated the BDEs by 4-17 kJ/mol in average. Sometimes, the underestimation by the DFT methods could be as high as 40-60 kJ/mol. Therefore, the DFT methods were more reliable for relative BDE calculations than for absolute BDE calculations. Finally, it was observed that the basis set effects on the BDEs calculated by the DFT methods were usually small except for the heteroatom-hydrogen BDEs.     

Dissociative photoionization mechanism of methanol isotopologues (CH3OH, CD3OH, CH3OD and CD3OD) by iPEPICO: energetics, statistical and non-statistical kinetics and isotope effects

The dissociative photoionization of energy selected methanol isotopologue (CH(3)OH, CD(3)OH, CH(3)OD and CD(3)OD) cations was investigated using imaging Photoelectron Photoion Coincidence (iPEPICO) spectroscopy. The first dissociation is an H/D-atom loss from the carbon, also confirmed by partial deuteration. Somewhat above 12 eV, a parallel H(2)-loss channel weakly asserts itself. At photon energies above 15 eV, in a consecutive hydrogen molecule loss to the first H-atom loss, the formation of CHO(+)/CDO(+) dominates as opposed to COH(+)/COD(+) formation. We see little evidence for H-atom scrambling in these processes. In the photon energy range corresponding to the B[combining tilde] and C[combining tilde] ion states, a hydroxyl radical loss appears yielding CH(3)(+)/CD(3)(+). Based on the branching ratios, statistical considerations and ab initio calculations, this process is confirmed to take place on the first electronically excited ?(2)A' ion state. Uncharacteristically, internal conversion is outcompeted by unimolecular dissociation due to the apparently weak Renner-Teller-like coupling between the X[combining tilde] and the ? ion states. The experimental 0 K appearance energies of the ions CH(2)OH(+), CD(2)OH(+), CH(2)OD(+) and CD(2)OD(+) are measured to be 11.646 ± 0.003 eV, 11.739 ± 0.003 eV, 11.642 ± 0.003 eV and 11.737 ± 0.003 eV, respectively. The E(0)(CH(2)OH(+)) = 11.6454 ± 0.0017 eV was obtained based on the independently measured isotopologue results and calculated zero point effects. The 0 K heat of formation of CH(2)OH(+), protonated formaldehyde, was determined to be 717.7 ± 0.7 kJ mol(-1). This yields a 0 K heat of formation of CH(2)OH of -11.1 ± 0.9 kJ mol(-1) and an experimental 298 K proton affinity of formaldehyde of 711.6 ± 0.8 kJ mol(-1). The reverse barrier to homonuclear H(2)-loss from CH(3)OH(+) is determined to be 36 kJ mol(-1), whereas for heteronuclear H(2)-loss from CH(2)OH(+) it is found to be 210 kJ mol(-1).     

Unimolecular decomposition of ethyl hydroperoxide: ab initio/Rice-Ramsperger-Kassel-Marcus theoretical prediction of rate constants

Alkyl hydroperoxides are found to be important intermediates in the combustion and oxidation processes of hydrocarbons. However, studies of ethyl hydroperoxide (CH(3)CH(2)OOH) are limited. In this work, kinetics and mechanisms for unimolecular decomposition of CH(3)CH(2)OOH have been investigated. The potential energy surface of decomposition reactions have first been predicted at the CCSD(T)/6-311+G(3df,2p)//B3LYP/6-311G(d,p) level. The results show that the formation of CH(3)CH(2)O + OH via O-O direct bond dissociation is dominant, the branching ratio of which is over 99% in the whole temperature range from 300 to 1000 K, and its rate constant can be expressed as k1 = 9.26 × 10(52)T(-11.91)exp(-26879/T) s(-1) at 1 atm. The rate constants of the reaction CH(3)CH(2)OOH → CH(3)CH(2)O + OH at different temperatures and pressures have been calculated, which can help us to comprehend the reactions of CH(3)CH(2)OOH at experimental conditions.     

A theoretical investigation of the gas-phase oxidation reaction of the saturated tert-butyl radical.

The radical-radical reaction mechanisms and dynamics of ground-state atomic oxygen [O(3P)] with the saturated tert-butyl radical (t-C4H9) are investigated using the density functional method and the complete basis set model. Two distinctive reaction pathways are predicted to be in competition: addition and abstraction. The barrierless addition of O(3P) to t-C4H9 leads to the formation of an energy-rich intermediate (OC4H9) on the lowest doublet potential energy surface, which undergoes subsequent direct elimination or isomerization-elimination leading to various products: C3H6O + CH3, iso-C4H8O + H, C3H7O + CH2, and iso-C4H8 + OH. The respective microscopic reaction processes examined with the aid of statistical calculations, predict that the major addition pathway is the formation of acetone (C3H6O) + CH3 through a low-barrier, single-step cleavage. For the direct, barrierless H-atom abstraction mechanism producing iso-C4H8 (isobutene) + OH, which was recently reported in gas-phase crossed-beam investigations, the reaction is described in terms of both an abstraction process (major) and a short-lived addition dynamic complex (minor).     

Negative-ion photoelectron spectroscopy, gas-phase acidity, and thermochemistry of the peroxyl radicals CH(3)OO and CH(3)CH(2)OO

Methyl, methyl-d(3), and ethyl hydroperoxide anions (CH(3)OO(-), CD(3)OO(-), and CH(3)CH(2)OO(-)) have been prepared by deprotonation of their respective hydroperoxides in a stream of helium buffer gas. Photodetachment with 364 nm (3.408 eV) radiation was used to measure the adiabatic electron affinities: EA[CH(3)OO, X(2)A' '] = 1.161 +/- 0.005 eV, EA[CD(3)OO, X(2)A' '] = 1.154 +/- 0.004 eV, and EA[CH(3)CH(2)OO, X(2)A' '] = 1.186 +/- 0.004 eV. The photoelectron spectra yield values for the term energies: Delta E(X(2)A' '-A (2)A')[CH(3)OO] = 0.914 +/- 0.005 eV, Delta E(X(2)A' '-A (2)A')[CD(3)OO] = 0.913 +/- 0.004 eV, and Delta E(X(2)A' '-A (2)A')[CH(3)CH(2)OO] = 0.938 +/- 0.004 eV. A localized RO-O stretching mode was observed near 1100 cm(-1) for the ground state of all three radicals, and low-frequency R-O-O bending modes are also reported. Proton-transfer kinetics of the hydroperoxides have been measured in a tandem flowing afterglow-selected ion flow tube (FA-SIFT) to determine the gas-phase acidity of the parent hydroperoxides: Delta(acid)G(298)(CH(3)OOH) = 367.6 +/- 0.7 kcal mol(-1), Delta(acid)G(298)(CD(3)OOH) = 367.9 +/- 0.9 kcal mol(-1), and Delta(acid)G(298)(CH(3)CH(2)OOH) = 363.9 +/- 2.0 kcal mol(-1). From these acidities we have derived the enthalpies of deprotonation: Delta(acid)H(298)(CH(3)OOH) = 374.6 +/- 1.0 kcal mol(-1), Delta(acid)H(298)(CD(3)OOH) = 374.9 +/- 1.1 kcal mol(-1), and Delta(acid)H(298)(CH(3)CH(2)OOH) = 371.0 +/- 2.2 kcal mol(-1). Use of the negative-ion acidity/EA cycle provides the ROO-H bond enthalpies: DH(298)(CH(3)OO-H) = 87.8 +/- 1.0 kcal mol(-1), DH(298)(CD(3)OO-H) = 87.9 +/- 1.1 kcal mol(-1), and DH(298)(CH(3)CH(2)OO-H) = 84.8 +/- 2.2 kcal mol(-1). We review the thermochemistry of the peroxyl radicals, CH(3)OO and CH(3)CH(2)OO. Using experimental bond enthalpies, DH(298)(ROO-H), and CBS/APNO ab initio electronic structure calculations for the energies of the corresponding hydroperoxides, we derive the heats of formation of the peroxyl radicals. The "electron affinity/acidity/CBS" cycle yields Delta(f)H(298)[CH(3)OO] = 4.8 +/- 1.2 kcal mol(-1) and Delta(f)H(298)[CH(3)CH(2)OO] = -6.8 +/- 2.3 kcal mol(-1).     

Bond dissociation energies and radical stabilization energies associated with model peptide-backbone radicals

Bond dissociation energies (BDEs) and radical stabilization energies (RSEs) have been calculated for a series of models that represent a glycine-containing peptide-backbone. High-level methods that have been used include W1, CBS-QB3, U-CBS-QB3, and G3X(MP2)-RAD. Simpler methods used include MP2, B3-LYP, BMK, and MPWB1K in association with the 6-311+G(3df,2p) basis set. We find that the high-level methods produce BDEs and RSEs that are in good agreement with one another. Of the simpler methods, RBMK and RMPWB1K achieve good accuracy for BDEs and RSEs for all the species that were examined. For monosubstituted carbon-centered radicals, we find that the stabilizing effect (as measured by RSEs) of carbonyl substituents (CX=O) ranges from 24.7 to 36.9 kJ mol(-1), with the largest stabilization occurring for the CH=O group. Amino groups (NHY) also stabilize a monosubstituted alpha-carbon radical, with the calculated RSEs ranging from 44.5 to 49.5 kJ mol(-1), the largest stabilization occurring for the NH2 group. In combination, NHY and CX=O substituents on a disubstituted carbon-centered radical produce a large stabilizing effect ranging from 82.0 to 125.8 kJ mol(-1). This translates to a captodative (synergistic) stabilization of 12.8 to 39.4 kJ mol(-1). For monosubstituted nitrogen-centered radicals, we find that the stabilizing effect of methyl and related (CH2Z) substituents ranges from 25.9 to 31.7 kJ mol(-1), the largest stabilization occurring for the CH3 group. Carbonyl substituents (CX=O) destabilize a nitrogen-centered radical relative to the corresponding closed-shell molecule, with the calculated RSEs ranging from -30.8 to -22.3 kJ mol(-1), the largest destabilization occurring for the CH=O group. In combination, CH2Z and CX=O substituents at a nitrogen radical center produce a destabilizing effect ranging from -19.0 to -0.2 kJ mol(-1). This translates to an additional destabilization associated with disubstitution of -18.6 to -7.8 kJ mol(-1).     

The Stability of α‐Hydroperoxyalkyl Radicals

High‐level ab initio and Born–Oppenheimer molecular dynamic calculations have been carried out on a series of hydroperoxyalkyl (α‐QOOH) radicals with the aim of investigating the stability and unimolecular decomposition mechanism into QO+OH of these species. Dissociation was shown to take place through rotation of the C?O(OH) bond rather than through elongation of the CO?OH bond. Through the C?O(OH) rotation, the unpaired electron of the radical overlaps with the electron density on the O?OH bond, and from this overlap the C=O π bond forms and the O?OH bond breaks spontaneously. The CH₂OOH, CH(CH₃)OOH, CH(OH)OOH, and α‐hydroperoxycycloheptadienyl radical were found to decompose spontaneously, but the CH(CHO)OOH has a decomposition energy barrier of 5.95 kcal mol^?1 owing to its steric and electronic features. The systems studied in this work provide the first insights into how structural and electronic effects govern the stabilizing influence on elusive α‐QOOH radicals.     

A computational study of the thermochemistry of bromine- and iodine-containing methanes and methyl radicals

Bromo- and iodomethanes and the corresponding halogenated methyl radicals have been investigated by ab initio methods. Geometries and vibrational frequencies were derived with quadratic configuration interaction methods at the QCISD/6-311G(d,p) level of theory, and energies via QCISD(T)/6-311+G(3df,2p). Core electrons were represented with relativistic effective potentials. Anharmonicity of the out-of-plane bending modes in the methyl radicals was taken into account by numerical integration of the Schr?dinger equation with potentials derived from relaxed scans of these modes. The results are in good accord with experimental data where available. Thermochemistry derived via isodesmic reactions referenced to CH3, CH4, and monohalomethanes yields excellent accord with new experiments on dihalomethanes and provides recommendations for the more poorly characterized tri- and tetrahalomethanes and halomethyl radicals. For the methanes CH2Br2, CHBr3, CBr4, CH2I2, CHI3, CI4, CH2BrI, CHBr2I, and CHBrI2 we compute DeltafH degrees (298) values of 4.3, 51.6, 110.6, 108.1, 208.5, 321.3, 56.8, 104.8, and 157.1 kJ mol(-1), respectively. For the methyl radicals CH2Br, CHBr2, CBr3, CH2I, CHI2, CI3, CHBrI, CBr2I, and CBrI2 we compute DeltafH degrees (298) values of 166.6, 191.7, 224.0, 217.2, 290.4, 369.1, 241.6, 320.8, and 272.3 kJ mol(-1), respectively. Recommended confidence limits are +/-3 kJ mol(-1) per Br or I atom. Trends in these values and the corresponding C-H bond strengths are discussed and compared with prior experiments, empirical estimation schemes, and ab initio calculations.     

Computational study of the halogen atom-benzene complexes

The structures of halogen atom-benzene complexes were investigated by modern DFT and ab initio computational methods. The spectroscopic properties of the complexes are also predicted and are in good agreement with experiment where such data have been reported. The fluorine atom-benzene complex is predicted to be a sigma complex due to the strength of a C-F bond. The chlorine atom-benzene complex is predicted to have an eta(1) pi complex structure, which is only slightly more favorable (1.1 kcal/mol with the BH&HLYP/6-311++G method including the ZPE correction) than a sigma complex but is significantly more stable (4.4 kcal/mol with the BH&HLYP/6-311++G method including the ZPE correction) than the eta(6) pi complex. The bromine and iodine benzene complexes are also predicted to prefer an eta(1) pi complex structure.     

Ab initio direct dynamics studies on the reactions of chlorine atom with CH3-nFnCH2OH (n = 1-3)

The hydrogen abstraction reactions of Cl atom with a series of fluorinated alcohols, i.e., CH(3-n)F(n)CH(2)OH + Cl (n = 1-3) (R1-R3) have been studied systematically by ab initio direct dynamics method and the canonical variational transition state theory (CVT). The potential energy surface information is calculated at the MP2/6-311G(d,p) level. Energies along the minimum energy paths are improved by a series of single-point calculations at the higher modified GAUSSIAN-2 (G2M) level of theory. Theoretical analysis shows that three kinds of hydrogen atoms can be abstracted from the reactants CH(2)FCH(2)OH and CHF(2)CH(2)OH, and for CF(3)CH(2)OH, two possible pathways are found. The rate constants for each reaction channel are evaluated by CVT with the small-curvature tunneling correction (SCT) over a wide range of temperature from 200 to 2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values for the reactions CHF(2)CH(2)OH + Cl and CF(3)CH(2)OH + Cl. However, for the reaction CH(2)FCH(2)OH + Cl, there is negative temperature dependence below 500 K, which is different from the experimental fitted. It is shown that in the low temperature ranges, the three reactions all proceed predominantly via H-abstraction from the methylene positions, and with the increase of the temperature the H-abstraction channels from the fluorinated-methyl positions should be taken into account, while the H-abstraction channels from the hydroxyl groups are negligible over the whole temperature ranges. Also, the reactivity decreases substantially with fluorine substitution at the methyl position of alcohol.