Thermochemical analysis and kinetics aspects for a chemical model for camphene ozonolysis

Abstract

In this work, a chemical model for the camphene ozonolysis, leading to carbonyl final products, is proposed and discussed on the basis of the thermochemical properties and kinetic data obtained at density functional theory levels of calculation. The mechanism is initiated by the electrophilic attack of ozone to the double bond in camphene leading to a 1,2,3-trioxolane intermediate, which decomposes to peroxy radicals and carbonyl compounds in a total of 10 elementary reactions. The thermodynamic properties (enthalpy and entropies differences) are calculated at 298 K. For the thermochemical evaluation, theoretical calculations are performed with the B3LYP, MPW1PW91, and mPW1K density functionals and the basis sets 6-31G(d), 6-31G(2d,2p), 6-31+G(d,p), and 6-31+G(2d,2p). Eventually, single point calculations adopting the 6-311++G(2d,2p) basis set are performed in order to improve the electronic energies. The enthalpy profiles suggest highly exothermic reactions for the individual steps, with a global enthalpy difference of -179.18 kcal mol(-1), determined at the B3LYP∕6-31+G(2d,2p) level. The Gibbs free energy differences for each step, at 298 K, calculated at the B3LYP∕6-311++G(2d,2p)∕∕B3LYP∕6-31+G(2d,2p) level, are used to estimate the composition of a final product mixture under equilibrium conditions as 58% of camphenilone and 42% of 6,6-dimethyl-ɛ-caprolactone-2,5-methylene. For the reaction kinetics, the bimolecular O(3) + camphene step is assumed to be rate determining in the global mechanism. A saddle point for the ozone addition to the double bond is located and rate constants are determined on the basis of the transition state theory. This saddle point is well represented by a loosely bound structure and corrections for the basis set superposition error (BSSE) are calculated, either by considering the effect over the geometry optimization procedure (here referred as CP1 procedure), or the effect of the BSSE over the electronic energy of a previously optimized geometry, included a posteriori (here referred as CP2). The rate constants, calculated at 298 K from the data obtained at the mPW1K∕6-31+G(d,p), CP1∕B3LYP∕∕6-31+G(2d,2p), and CP2∕B3LYP∕∕6-31+G(2d,2p) levels (3.62 × 10(-18), 1.12 × 10(-18), and 1.39 × 10(-18) cm(3) molecule(-1) s(-1)), are found in good agreement with the available experimental data at the same temperature, 0.9 × 10(-18) cm(3) molecule(-1) s(-1) [R. Atkinson, S. M. Aschmann, and J. Arey, Atmos. Environ. 24, 2647 (1990)]. The importance of the BSSE corrections for the final rate constants must be pointed out. Furthermore, this work will contribute to a better understanding of the chemistry of monoterpenes in the atmosphere, as well as the implications for the phenomena of pollution.

A Computational Study of the Ozonolysis of Phenanthrene

“…Furthermore, the bond length of O-O bond cleavage in TS 2 A is 2.109 Å, using B3LYP/6-31G(2df,p) method. In general, the bond lengths are in good agreement with the ozonolysis of monoterpenes [43][44][45].…”

OH- and O3-initiated atmospheric degradation of camphene: temperature dependent rate coefficients, product yields and mechanisms

“…There have also been computational studies on the kinetics of the ozonolysis of camphene. 18,19 In addition, Peeters et al 20 have estimated the rate coefficient for the OH-camphene reaction using a Structure-Activity Relationship (SAR).…”

“…Reissell et al 23 have measured the acetone yield from the reaction of OH with camphene at 298 K and 740 Torr in a Teon chamber using GC/FID and GC/MS. Finally, Oliveira and Bauerfeldt 18 in a computational study present a chemical model for the ozonolysis of camphene leading to carbonyl nal products.…”

Ab initio investigation of 2,2′‐bis(4‐trifluoromethylphenyl)‐5,5′‐bithiazole for the design of efficient organic field‐effect transistors

“…The molecular geometry in the ground state has been optimized within the framework of density functional theory (DFT) using exchange functional given by Becke and correlation functional collectively given by Lee et al, expressed as B3LYP in combination with 6‐31G(d,p) basis sets . The dipole moment has also been calculated for the ground‐state geometry.…”