Development of a Detailed Kinetic Model for the Oxidation of n-Butane in the Liquid Phase.

The chemistry underlying liquid-phase oxidation of organic compounds, the main cause of their aging, is characterized by a free-radical chain reaction mechanism. The rigorous simulation of these phenomena requires the use of detailed kinetic models that contain thousands of species and reactions. The development of such models for the liquid phase remains a challenge as solvent-dependent thermokinetic parameters have to be provided for all the species and reactions of the model. Therefore, accurate and high-throughput methods to generate these data are required. In this work, we propose new methods to generate these data, and we apply them for the development of a detailed chemical kinetic model for n-butane autoxidation, which is then validated against literature data. Our approach for model development is based on the work of Jalan et al. [J. Phys. Chem. B 2013, 117, 2955-2970] who used Gibbs free energies of solvation [ΔsolvG(T)] to correct the data of the gas-phase kinetic model. In our approach, an equation of state (EoS) is used to compute ΔsolvG as a function of temperature for all the chemical species in the mechanism. Currently, ΔsolvG(T) of free radicals cannot be computed with an EoS and it was calculated for their parent molecule (H-atom added on the radical site). Theoretical calculations with the implicit solvent model were performed to quantify the impact of this assumption and showed that it is acceptable for radicals in n-butane and probably in all n-alkanes. New rate rules were proposed for the most important reactions of the model, based on theoretical calculations and the literature data. The developed detailed kinetic model for n-butane autoxidation is the first proposed model in the literature and was validated against the experimental data from the literature. Simulations showed that the main autoxidation products, sec-butyl hydroperoxides and 2-butanol, are produced from H-abstractions from n-butane by sec-C4H9OO radicals and the C4H9OO + C4H9OO reaction, respectively. The uncertainty of the product ratio ("butanone + 2-butanol"/"2-butoxy + 2-butoxy") of the latter reaction remains high in the literature, and our simulations suggest a 1:1 ratio in n-butane solvent.

Theoretical study of the pyrolysis of ß-1,4-xylan: a detailed investigation on unimolecular concerted reactions.

A theoretical study of the thermal decomposition of ß-1,4-xylan, a model polymer of hemicelluloses, is proposed for the first time. A mechanism based on unimolecular concerted reactions is elaborated in a comprehensive way. Elementary reactions, such as dehydrations, retro-aldol, retro Diels-Alder, retro-ene, glycosidic bond fissions, isomerizations, etc., are applied to ß-1,4-xylan, as well as to the fragments formed. At each stage of the construction of the mechanism, the fragments previously retained are decomposed and the low energy paths are selected to define new fragments. Energy barriers are computed at the CBS-QB3 level of theory and rate coefficients of important reactions are calculated. It is shown that the main reaction pathways can be modelled by reactions involving two specific fragments, which react in closed sequences, similarly to chain-propagating reactions. The proposed reaction scheme allows to predict important species observed during the pyrolysis of xylan, such as aldehydes or CO. In addition, we show that dehydrations require high activation energy and cannot compete with the other reactions. Therefore, it seems difficult to explain, by means of unimolecular homogeneous gas phase reactions, the significant formation of specific species such as furfural as reported by several authors.

Kinetic modeling of the thermal destruction of lewisite.

Organoarsenic compounds have been widely used as pesticides and chemical agents. Lewisite (C2H2AsCl3), a blister agent, is a model of such compounds. A comprehensive detailed kinetic mechanism of combustion has been developed based on theoretical investigations. A benchmark allowed to select an appropriate methodology able to deal with such a heavy atom as As with precision and reasonable computational times. The density functional theory (DFT) method ωB97X-D was found to give the best results on target data. Core pseudo potentials were used for arsenic with the cc-pVTZ-PP basis set, whereas Def2-TZVP basis set was used for other atoms. The mechanism of the decomposition of lewisite includes all reactions involved in thermal decomposition and combustion mechanisms, including molecular and radical intermediates, and the decomposition reactions of small species containing arsenic. Simulation shows that lewisite decomposition starts around 700 K and is very little sensitive to the presence of oxygen since the radical reactions involve mainly very reactive Cl-atoms as chain carriers. The main reaction paths have been derived. As experimentally observed, AsCl3 is the main arsenic product produced almost in one-to-one yield, whereas acetylene is an important hydrocarbon product in pyrolysis. In combustion, several arsenic oxides, eventually chlorinated, are produced, which toxicity need to be assessed.

Combustion and Pyrolysis Kinetics of Chloropicrin.

Chloropicrin (CCl3NO2) is widely used in agriculture as a pesticide, weed-killer, fungicide or nematicide. It has also been used as a chemical agent during World War I. The precise understanding of its combustion chemistry for destruction processes or in the event of accidental fire of stored reserves is a major safety issue. A detailed chemical kinetic model for the combustion and pyrolysis of chloropicrin is proposed for the first time. A large number of thermo-kinetic parameters were calculated using quantum chemistry and reaction rate theory. The model was validated against experimental pyrolysis data available in the literature. It was shown that the degradation of chloropicrin is ruled by the breaking of the C-N bond followed by the oxidation of the trichloromethyl radical by NO2 through the formation of the adduct CCl3ONO, which can decompose to NO, chlorine atom, and phosgene. Phosgene is much more stable than chloropicrin and its decomposition starts at much higher temperatures. Combustion and pyrolysis simulations were also compared and demonstrated that the addition of oxygen has very little effect on the reactivity or product distribution due to the absence of hydrogen atoms in chloropicrin.

Intramolecular effects on the kinetics of unimolecular reactions of ß-HOROO˙ and HOQ˙OOH radicals.

A theoretical study describing the influence of intramolecular effects on the energy barriers and rate constants of unimolecular reactions involving ß-HOROO˙ and HOQ˙OOH radicals is proposed. The reactions considered are HO2˙ elimination, the Waddington mechanism, H-shift, cyclic ether formation and ß-scission. All the calculations are performed at the CBS-QB3 level of theory along with canonical transition state theory and statistical thermodynamics, including a specific treatment of hindered rotors. Several structural parameters are investigated, such as the location of the hydroxyl function in the cyclic transition states or the substitution of H atoms by alkyl groups on carbon atoms involved in the reaction coordinate. It is shown that these molecular systems involve numerous transition states, especially for reactions such as 1,5 or 1,6 H-shift, and that, a priori simplification is not possible. It is also shown that the position of the -OH group in the transition state can largely modify both the barrier heights and the rate constants. However, opposite trends can be observed depending on the competition between energetic and entropic effects. Similar observations are made when H atoms are replaced by methyl or alkyl groups. These results can largely be explained by intramolecular effects such as hydrogen bonds, stabilization effects (from -OH or -CH3 groups), steric influences and by the coupling between them. The last point renders the classic establishment of the structure-reactivity relationship challenging.

Quantum Chemical Study of the Thermochemical Properties of Organophosphorous Compounds.

Organophosphorous compounds are involved in many toxic compounds such as fungicides, pesticides, or chemical warfare nerve agents. The understanding of the decomposition chemistry of these compounds in the environment is largely limited by the scarcity of thermochemical data. Because of the high toxicity of many of these molecules, experimental determination of their thermochemical properties is very difficult. In this work, standard gas-phase thermodynamic data, i.e., enthalpies of formation (ΔfH298°), standard entropies (S298°), and heat capacities (Cp°(T)), were determined using quantum chemical calculations and more specifically the CBS-QB3 composite method, which was found to be the best compromise between precision and calculation time among high accuracy composite methods. A large number of molecules was theoretically investigated, involving trivalent and pentavalent phosphorus atoms, and C, H, O, N, S, and F atoms. These data were used to propose 83 original groups, used in the semiempirical group contribution method proposed by Benson. Thanks to these latter group values, thermochemical properties of several nerve agents, common pesticides and herbicides have been evaluated. Bond dissociations energies (BDE), useful for the analysis the thermal stability of the compounds, were also determined in several molecules of interest.

Theoretical kinetic study of the reactions of cycloalkylperoxy radicals.

Reactions of alkyl radicals with oxygen are key reactions in the low-temperature oxidation of hydrocarbons, but they have not been extensively studied yet in the case of cycloalkanes. Isomerizations of cycloalkylperoxy radicals and formation of cyclic ethers are especially important. In the present work, a theoretical study of the gas-phase reactions of cyclopentylperoxy and cyclohexylperoxy radicals has been carried out by means of quantum chemical calculations at the CBS-QB3 level. Computations on cyclopentylperoxy decomposition pathways are reported here for the first time. Thermochemical data have been obtained by means of isodesmic reactions, and the contribution of hindered rotors has been explicitly taken into account. Transition state theory has been used to calculate rate constants for all the elementary reactions. Three-parameter Arrhenius expressions have been derived in the temperature range 300-1000 K. Tunneling effects have been accounted for in the case of H-atom transfers. Our results compare well with experimental data and previous calculations available in the literature. In particular, the predicted rate constants for processes involving cyclohexylperoxy radicals, which have been introduced in a reaction mechanism scheme proposed before, exhibit excellent agreement with experiments at low and intermediate temperatures.

Theoretical kinetic study of thermal unimolecular decomposition of cyclic alkyl radicals.

Whereas many studies have been reported on the reactions of aliphatic hydrocarbons, the chemistry of cyclic hydrocarbons has not been explored extensively. In the present work, a theoretical study of the gas-phase unimolecular decomposition of cyclic alkyl radicals was performed by means of quantum chemical calculations at the CBS-QB3 level of theory. Energy barriers and high-pressure-limit rate constants were calculated systematically. Thermochemical data were obtained from isodesmic reactions, and the contribution of hindered rotors was taken into account. Classical transition state theory was used to calculate rate constants. The effect of tunneling was taken into account in the case of CH bond breaking. Three-parameter Arrhenius expressions were derived in the temperature range of 500-2000 K at atmospheric pressure, and the CC and CH bond breaking reactions were studied for cyclic alkyl radicals with a ring size ranging from three to seven carbon atoms, with and without a lateral alkyl chain. For the ring-opening reactions, the results clearly show an increase of the activation energy as the pi bond is being formed in the ring (endo ring opening) in contrast to the cases in which the pi bond is formed on the side chain (exo ring opening). These results are supported by analyses of the electronic charge density that were performed with Atoms in Molecules (AIM) theory. For all cycloalkyl radicals considered, CH bond breaking exhibits larger activation energies than CC bond breaking, except for cyclopentyl for which the ring-opening and H-loss reactions are competitive over the range of temperatures studied. The theoretical results compare rather well with the experimental data available in the literature. Evans-Polanyi correlations for CC and CH beta-scissions in alkyl and cycloalkyl free radicals were derived. The results highlight two different types of behavior depending on the strain energy in the reactant.

Detailed kinetic study of the ring opening of cycloalkanes by CBS-QB3 calculations.

This work reports a theoretical study of the gas-phase unimolecular decomposition of cyclobutane, cyclopentane and cyclohexane by means of quantum chemical calculations. A biradical mechanism has been envisaged for each cycloalkane, and the main routes for the decomposition of the biradicals formed have been investigated at the CBS-QB3 level of theory. Thermochemical data(DeltaHf(o), S(o), Cp(o)) for all the involved species have been obtained by means of isodesmic reactions. The contribution of hindered rotors has also been included. Activation barriers of each reaction have been analyzed to assess the energetically most favorable pathways for the decomposition of biradicals. Rate constants have been derived for all elementary reactions using transition-state theory at 1 atm and temperatures ranging from 600 to 2000 K. Global rate constant for the decomposition of the cyclic alkanes in molecular products have been calculated. Comparison between calculated and experimental results allowed us to validate the theoretical approach. An important result is that the rotational barriers between the conformers, which are usually neglected, are of importance in decomposition rate of the largest biradicals. Ring strain energies (RSE) in transition states for ring opening have been estimated and show that the main part of RSE contained in the cyclic reactants is removed upon the activation process.