Theoretical study of the oxidation mechanisms of naphthalene initiated by hydroxyl radicals: the OH-addition pathway.

The oxidation mechanisms of naphthalene by OH radicals under inert (He) conditions have been studied using density functional theory along with various exchange-correlation functionals. Comparison has been made with benchmark CBS-QB3 theoretical results. Kinetic rate constants were correspondingly estimated by means of transition state theory and statistical Rice-Ramsperger-Kassel-Marcus (RRKM) theory. Comparison with experiment confirms that, on the OH-addition reaction pathway leading to 1-naphthol, the first bimolecular reaction step has an effective negative activation energy around -1.5 kcal mol(-1), whereas this step is characterized by an activation energy around 1 kcal mol(-1) on the OH-addition reaction pathway leading to 2-naphthol. Effective rate constants have been calculated according to a steady state analysis upon a two-step model reaction mechanism. In line with experiment, the correspondingly obtained branching ratios indicate that, at temperatures lower than 410 K, the most abundant product resulting from the oxidation of naphthalene by OH radicals must be 1-naphthol. The regioselectivity of the OH(•)-addition onto naphthalene decreases with increasing temperatures and decreasing pressures. Because of slightly positive or even negative activation energies, the RRKM calculations demonstrate that the transition state approximation breaks down at ambient pressure (1 bar) for the first bimolecular reaction steps. Overwhelmingly high pressures, larger than 10(5) bar, would be required for restoring to some extent (within ∼5% accuracy) the validity of this approximation for all the reaction channels that are involved in the OH-addition pathway. Analysis of the computed structures, bond orders, and free energy profiles demonstrate that all reaction steps involved in the oxidation of naphthalene by OH radicals satisfy Leffler-Hammond's principle. Nucleus independent chemical shift indices and natural bond orbital analysis also show that the computed activation and reaction energies are largely dictated by alterations of aromaticity, and, to a lesser extent, by anomeric and hyperconjugative effects.

KiSThelP: a program to predict thermodynamic properties and rate constants from quantum chemistry results.

Kinetic and Statistical Thermodynamical Package (KiSThelP) is a cross-platform free open-source program developed to estimate molecular and reaction properties from electronic structure data. To date, three computational chemistry software formats are supported (Gaussian, GAMESS, and NWChem). Some key features are: gas-phase molecular thermodynamic properties (offering hindered rotor treatment), thermal equilibrium constants, transition state theory rate coefficients (transition state theory (TST), variational transition state theory (VTST)) including one-dimensional (1D) tunnelling effects (Wigner, and Eckart) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate constants, for elementary reactions with well-defined barriers. KiSThelP is intended as a working tool both for the general public and also for more expert users. It provides graphical front-end capabilities designed to facilitate calculations and interpreting results. KiSThelP enables to change input data and simulation parameters directly through the graphical user interface and to visually probe how it affects results. Users can access results in the form of graphs and tables. The graphical tool offers customizing of 2D plots, exporting images and data files. These features make this program also well-suited to support and enhance students learning and can serve as a very attractive courseware, taking the teaching content directly from results in molecular and kinetic modelling.

Theoretical study of H-abstraction reactions from CH3Cl and CH3Br molecules by ClO and BrO radicals.

The rate constants of the H-abstraction reactions from CH(3)Cl and CH(3)Br molecules by ClO and BrO radicals have been estimated over the temperature range of 300-2500 K using four different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using B3LYP and MP2 methods combined with the cc-pVTZ basis set. Single-point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using the cc-pVTZ and cc-pVQZ basis sets. Canonical transition-state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature. In order to choose the appropriate levels of theory with chlorine- and bromine-containing species, the reference reaction Cl ((2)P(3/2)) + CH(3)Cl → HCl + CH(2)Cl (R(ref)) was first theoretically studied because its kinetic parameters are well-established from numerous experiments, evaluation data, and theoretical studies. The kinetic parameters of the reaction R(ref) have been determined accurately using the CCSD(T)/cc-pVQZ//MP2/cc-pVTZ level of theory. This level of theory has been used for the rate constant estimation of the reactions ClO + CH(3)Cl (R(1)), ClO + CH(3)Br (R(2)), BrO + CH(3)Cl (R(3)), and BrO + CH(3)Br (R(4)). Six-parameter Arrhenius expressions have been obtained by fitting to the computed rate constants of these four reactions (including cis and trans pathways) over the temperature range of 300-2500 K.

Thermochemical data and additivity group values for ten species of o-xylene low-temperature oxidation mechanism.

o-Xylene could be a good candidate to represent the family of aromatic hydrocarbons in a surrogate fuel. This study uses computational chemistry to calculate standard enthalpies of formation at 298 K, Δ(f)H°(298 K), standard entropies at 298 K, S°(298 K), and standard heat capacities C(p)°(T) over the temperature range 300 K to 1500 K for ten target species present in the low-temperature oxidation mechanism of o-xylene: o-xylene (1), 2-methylbenzyl radical (2), 2-methylbenzylperoxy radical (3), 2-methylbenzyl hydroperoxide (4), 2-(hydroperoxymethyl)benzyl radical (5), 2-(hydroperoxymethyl)benzaldehyde (6), 1-ethyl-2-methylbenzene (7), 2,3-dimethylphenol (8), 2-hydroxybenzaldehyde (9), and 3-hydroxybenzaldehyde (10). Δ(f)H°(298 K) values are weighted averages across the values calculated using five isodesmic reactions and five composite calculation methods: CBS-QB3, G3B3, G3MP2, G3, and G4. The uncertainty in Δ(f)H°(298 K) is also evaluated. S°(298 K) and C(p)°(T) values are calculated at B3LYP/6-311G(d,p) level of theory from molecular properties and statistical thermodynamics through evaluation of translational, rotational, vibrational, and electronic partition functions. S°(298 K) and C(p)°(300 K) values are evaluated using the rigid-rotor-harmonic-oscillator model. C(p)°(T) values at T ≥ 400 K are calculated by treating separately internal rotation contributions and translational, external rotational, vibrational, and electronic contributions. The thermochemical properties of six target species are used to develop six new additivity groups taking into account the interaction between two substituents in ortho (ORT/CH2OOH/ME, ORT/ET/ME, ORT/CHO/OH, ORT/CHO/CH2OOH) or meta (MET/CHO/OH) positions, and the interaction between three substituents (ME/ME/OH123) located one beside the other (positions numbered 1, 2, 3) for two- or three-substituted benzenic species. Two other additivity groups are also developed using the thermochemical properties of benzenic species taken from the literature: the C/CB/H2/OO and the CB/CO groups. These groups extend the capacities of the group additivity method to deal with substituted benzenic species.

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2)).

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.

Theoretical study of the gas-phase reactions of iodine atoms ((2)P(3/2)) with H(2), H(2)O, HI, and OH.

The rate constants of the reactions of iodine atoms with H(2), H(2)O, HI, and OH have been estimated using 39, 21, 13, and 39 different levels of theory, respectively, and have been compared to the available literature values over the temperature range of 250-2500 K. The aim of this methodological work is to demonstrate that standard theoretical methods are adequate to obtain quantitative rate constants for the reactions involving iodine-containing species. Geometry optimizations and vibrational frequency calculations are performed using three methods (MP2, MPW1K, and BHandHLYP) combined with three basis sets (cc-pVTZ, cc-pVQZ, and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVnZ (n = T, Q, and 5), aug-cc-pVnZ (n = T, Q, and 5), 6-311G(d,p), 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Canonical transition state theory with a simple Wigner tunneling correction is used to predict the rate constants as a function of temperature. CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory provide accurate kinetic rate constants when compared to available literature data. The use of the CCSD(T)/cc-pVQZ//MP2/cc-pVTZ and CCSD(T)/6-311++G(3df,3pd) levels of theory allows one to obtain a better agreement with the literature data for all reactions with the exception of the I + H(2) reaction R(1) . This computational procedure has been also used to predict rate constants for some reactions where no available experimental data exist. The use of quantum chemistry tools could be therefore extended to other elements and next applied to develop kinetic networks involving various fission products, steam, and hydrogen in the absence of literature data. The final objective is to implement the kinetics of gaseous reactions in the ASTEC (Accident Source Term Evaluation Code) code to improve speciation of fission transport, which can be transported along the Reactor Coolant System (RCS) of a Pressurized Water Reactor (PWR) in case of a severe accident.

A CASPT2 theoretical study of the kinetics of the 2-, 3-, and 4-methylbenzylperoxy radical isomerization.

The rate constants of the 2-, 3-, and 4-methylbenzylperoxy isomerization reactions have been computed using the elaborated CASPT2 method. Geometry optimizations and vibrational frequency calculations are performed with two methods (B3LYP and MPW1K) combined with the cc-pVDZ and 6-31+G(d,p) basis sets, respectively. Single-point energy calculations are performed at the CASPT2/ANO-L-VDZP//B3LYP/cc-pVDZ level of theory as recommended by Canneaux et al. (J. Phys. Chem. A 2008, 112, 6045). Canonical transition-state theory with a simple Wigner tunneling correction is used to predict the high-pressure limit rate constants as a function of temperature. They are given by the following relations for the 2-, 3-, and 4-methylbenzylperoxy (MBP) (1,3s) isomerizations and for the 2-methylbenzylperoxy (1,6p) isomerization, respectively: k(2-MBP(1,3s))(600-2000 K) (in s(-1)) = (3.33 x 10(10))T(0.79) exp((-142.6 in kJ mol(-1))/RT); k(3-MBP(1,3s))(600-2000 K) (in s(-1)) = (0.74 x 10(10))T(0.79) exp((-130.7 in kJ mol(-1))/RT); k(4-MBP(1,3s))(600-2000 K) (in s(-1)) = (1.12 x 10(10))T(0.79) exp((-133.6 in kJ mol(-1))/RT); k(2-MBP(1,6p))(600-2000 K) (in s(-1)) = (5.10 x 10(8))T(0.85) exp((-87.1 in kJ mol(-1))/RT). These parameters can be used in the thermokinetic models involving aromatic compounds at high pressure. In the case of the 2-methylbenzylperoxy radical, the (1,6p) H-atom transfer reaction is consistently the most important channel over the studied temperature range, and the (1,3s) H-atom transfer reaction is not energetically favored.

A theoretical study of the kinetics of the benzylperoxy radical isomerization.

The rate constant of the benzylperoxy isomerization reaction has been computed using 54 different levels of theory and has been compared to the experimental value reported at 773 K. The aim of this methodology work is to demonstrate that standard theoretical methods are not adequate to obtain quantitative rate constants for the reaction under study. The use of the elaborated CASPT2 method is essential to estimate a quantitative rate constant. Geometry optimizations and vibrational frequency calculations are performed using three different methods (B3LYP, MPW1K, and MP2) and six different basis sets (6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), and cc-pVDZ). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the 6-31G(d,p) basis set, and with the CASPT2 level of theory with the ANO-L-VDZP basis set. Canonical transition-state theory with a simple Wigner tunneling correction is used to predict the high-pressure limit rate constants as a function of temperature. We recommend the use of the CASPT2/ANO-L-VDZP//B3LYP/cc-pVDZ level of theory to compute the temperature dependence of the rate constant of the four-center isomerization of the benzylperoxy radical. It is given by the following relation: k(600-2000 K) (in s (-1)) = (1.29 x 10 (10)) T (0.79) exp[(-133.1 in kJ mol (-1))/ RT]. These parameters can be used in the thermokinetic models involving aromatic compounds at high pressure. This computational procedure can be extended to predict rate constants for other similar reactions where no available experimental data exist.