scholarly journals Direct Kinetic Measurements and Master Equation Modelling of the Unimolecular Decomposition of Resonantly-Stabilized CH2CHCHC(O)OCH3 Radical and an Upper Limit Determination for CH2CHCHC(O)OCH3 + O2 Reaction

2020 ◽  
Vol 234 (7-9) ◽  
pp. 1251-1268 ◽  
Author(s):  
Satya Prakash Joshi ◽  
Prasenjit Seal ◽  
Timo Theodor Pekkanen ◽  
Raimo Sakari Timonen ◽  
Arrke J. Eskola

AbstractMethyl-Crotonate (MC, (E)-methylbut-2-enoate, CH3CHCHC(O)OCH3) is a potential component of surrogate fuels that aim to emulate the combustion of fatty acid methyl ester (FAME) biodiesels with significant unsaturated FAME content. MC has three allylic hydrogens that can be readily abstracted under autoignition and combustion conditions to form a resonantly-stabilized CH2CHCHC(O)OCH3 radical. In this study we have utilized photoionization mass spectrometry to investigate the O2 addition kinetics and thermal unimolecular decomposition of CH2CHCHC(O)OCH3 radical. First we determined an upper limit for the bimolecular rate coefficient of CH2CHCHC(O)OCH3 + O2 reaction at 600 K (k ≤ 7.5 × 10−17 cm3 molecule−1 s−1). Such a small rate coefficient suggest this reaction is unlikely to be important under combustion conditions and subsequent efforts were directed towards measuring thermal unimolecular decomposition kinetics of CH2CHCHC(O)OCH3 radical. These measurements were performed between 750 and 869 K temperatures at low pressures (<9 Torr) using both helium and nitrogen bath gases. The potential energy surface of the unimolecular decomposition reaction was probed at density functional (MN15/cc-pVTZ) level of theory and the electronic energies of the stationary points obtained were then refined using the DLPNO-CCSD(T) method with the cc-pVTZ and cc-pVQZ basis sets. Master equation simulations were subsequently carried out using MESMER code along the kinetically important reaction pathway. The master equation model was first optimized by fitting the zero-point energy corrected reaction barriers and the collisional energy transfer parameters $\Delta{E_{{\text{down}},\;{\text{ref}}}}$ and n to the measured rate coefficients data and then utilize the constrained model to extrapolate the decomposition kinetics to higher pressures and temperatures. Both the experimental results and the MESMER simulations show that the current experiments for the thermal unimolecular decomposition of CH2CHCHC(O)OCH3 radical are in the fall-off region. The experiments did not provide definite evidence about the primary decomposition products.

2020 ◽  
Vol 234 (7-9) ◽  
pp. 1233-1250 ◽  
Author(s):  
Arrke J. Eskola ◽  
Mark A. Blitz ◽  
Michael J. Pilling ◽  
Paul W. Seakins ◽  
Robin J. Shannon

AbstractThe rate coefficient for the unimolecular decomposition of CH3OCH2, k1, has been measured in time-resolved experiments by monitoring the HCHO product. CH3OCH2 was rapidly and cleanly generated by 248 nm excimer photolysis of oxalyl chloride, (ClCO)2, in an excess of CH3OCH3, and an excimer pumped dye laser tuned to 353.16 nm was used to probe HCHO via laser induced fluorescence. k1(T,p) was measured over the ranges: 573–673 K and 0.1–4.3 × 1018 molecule cm−3 with a helium bath gas. In addition, some experiments were carried out with nitrogen as the bath gas. Ab initio calculations on CH3OCH2 decomposition were carried out and a transition-state for decomposition to CH3 and H2CO was identified. This information was used in a master equation rate calculation, using the MESMER code, where the zero-point-energy corrected barrier to reaction, ΔE0,1, and the energy transfer parameters, ⟨ΔEdown⟩ × Tn, were the adjusted parameters to best fit the experimental data, with helium as the buffer gas. The data were combined with earlier measurements by Loucks and Laidler (Can J. Chem.1967, 45, 2767), with dimethyl ether as the third body, reinterpreted using current literature for the rate coefficient for recombination of CH3OCH2. This analysis returned ΔE0,1 = (112.3 ± 0.6) kJ mol−1, and leads to $k_{1}^{\infty}(T)=2.9\times{10^{12}}$ (T/300)2.5 exp(−106.8 kJ mol−1/RT). Using this model, limited experiments with nitrogen as the bath gas allowed N2 energy transfer parameters to be identified and then further MESMER simulations were carried out, where N2 was the buffer gas, to generate k1(T,p) over a wide range of conditions: 300–1000 K and N2 = 1012–1025 molecule cm−3. The resulting k1(T,p) has been parameterized using a Troe-expression, so that they can be readily be incorporated into combustion models. In addition, k1(T,p) has been parametrized using PLOG for the buffer gases, He, CH3OCH3 and N2.


2013 ◽  
Vol 91 (4) ◽  
pp. 241-247 ◽  
Author(s):  
Yan Zhao ◽  
Haitao Sun ◽  
Renjun Wang ◽  
Fei Gao

The reaction mechanisms of vinyl acetate with OH radicals in the atmosphere have been studied using the density functional theory method. The geometry parameters and frequencies of all of the stationary points are calculated at the MPWB1K level with the 6-31G(d,p) basis sets. The single-point energy calculations are carried out at the MPWB1K/6-311+G(3df,2pd) level. The detailed profiles of the potential energy surfaces for the reactions are constructed. Two OH addition and three H abstraction reaction pathways are considered for the reaction of vinyl acetate with OH radicals. The theoretical study shows that the most energetically favorable isomer is that of OH addition to the terminal carbon positions (C1 atom). The α-ester rearrangement, which is characteristic of ester oxidation processes, is confirmed to be thermodynamically and kinetically favorable. The main products of the OH-initiated atmospheric oxidation of vinyl acetate are formaldehyde, formic acetic anhydride, and acetic acid.


2014 ◽  
Vol 92 (7) ◽  
pp. 598-604 ◽  
Author(s):  
Yan Zhao ◽  
Xiaomin Sun ◽  
Wenxing Wang ◽  
Laixiang Xu

The chlorine-initiated reaction mechanism of ethyl formate in the atmosphere was investigated using the density functional theory method. The geometry parameters and frequencies of all of the stationary points were calculated at the B3LYP/ 6-31G(d,p) level. The single-point energy calculations were carried out at different levels, including MP2/6-31G(d), MP2/6-311++G(d,p), and CCSD(T)/6-31G(d). A detailed oxidation mechanism is provided and discussed. Present results show that α-ester rearrangement reaction and the O2 direct abstraction from IM6 (HC(O)OCH(O)CH3) are the more favorable pathway and are competitive. The 1,4-H shift isomerization of IM6 proved to be feasible under general atmospheric conditions. The decomposition of IM18 (CH3CH2OC(O)O) is favorable both thermodynamically and kinetically. Canonical variational transition theory with small-curvature tunneling correction was employed to predict the rate constants. The overall rate constant of ethyl formate at 298 K is 8.63 × 10−12 cm3 molecule−1 s−1. The Arrhenius equations of rate constants at the temperature range of 200–380 K were fitted.


2013 ◽  
Vol 101 (9) ◽  
pp. 561-570
Author(s):  
R. Polly ◽  
B. Schimmelpfennig ◽  
M. Flörsheimer ◽  
Th. Rabung ◽  
T. Kupcik ◽  
...  

Summary Sorption plays a major role in the safety assessment of nuclear waste disposal. In the present theoretical study we focused on understanding the interaction of trivalent lanthanides and actinides (La3+, Eu3+ and Cm3+) with the corundum (110) surface. Optimization of the structures were carried out using density functional theory with different basis sets. Additionally, Møller-Plesset perturbation theory of second order was used for single point energy calculations. We studied the structure of different inner-sphere complexes depending on the surface deprotonation and the number of water molecules in the first coordination shell. The most likely structure of the inner-sphere complex (tri- or tetradentate) was predicted. For the calculations we used a cluster model for the surface. By deprotonating the cluster a chemical environment at elevated pH values was mimicked. Our calculations predict the highest stability for a tetradentate inner-sphere surface complexes with five water molecules remaining in the first coordination sphere of the metal ions. The formation of the inner-sphere complexes is favored when a coordination takes place with at most one deprotonated surface aluminol group located beneath the inner-sphere complex. The mutual interaction between sorbing metal ions at the surface is studied as well. The minimal possible distance between two inner-sphere sorbed metal ions at the surface was determined to be 530 pm.


2004 ◽  
Vol 69 (1) ◽  
pp. 13-33 ◽  
Author(s):  
Dimitrios A. Pantazis ◽  
Athanassios C. Tsipis ◽  
Constantinos A. Tsipis

Density functional calculations at the B3LYP level of theory, using the 6-31G(d) and 6-311+G(3df) basis sets, provide a satisfactory description of the geometric and energetic reaction profile of the Fe + CO2 → FeO + CO reaction. The reaction is predicted to be endothermic by 23.24 kcal/mol at the B3LYP/6-311+G(3df)//B3LYP/6-31G(d) level of theory and to proceed by formation of either a Fe(η2-OCO) or a Fe(η3-OCO) intermediate. The Fe(η2-OCO) intermediate in the 5A' ground state is weakly bound with respect to Fe(5D) and CO2 dissociation products by 0.78 (2.88) kcal/mol at the B3LYP/6-31G(d) (B3LYP/6-311+ G(3df)//B3LYP/6-31G(d)) levels of theory. In contrast, the Fe(η3-OCO) intermediate in the 5A1 ground state is unbound with respect to Fe(5D) and CO2 dissociation products by 8.27 (11.15) kcal/mol at the same levels of theory. However, both intermediates are strongly bound relative to the separated Fe+(6D) and [CO2]- anion; the computed bond dissociation energies for the Fe(η2-OCO) and Fe(η3-OCO) intermediates are 207.33 and 198.28 kcal/mol in terms of ∆E0 at the B3LYP/6-31G(d), respectively. In the Fe(η2-OCO) and Fe(η3-OCO) intermediates, an intramolecular insertion reaction of the Fe atom to O-C bond takes place yielding the isomeric OFe(η1-CO) and OFe(η1-OC) products, respectively, with a relatively low activation barrier of 25.24 (21.69) and 26.36 (23.38) kcal/mol at the B3LYP/6-31G(d) (B3LYP/6-311+G(3df)//B3LYP/6-31G(d)) levels of theory, respectively. The calculated structures, relative stability and bonding properties of all stationary points are discussed with respect to computed electronic and spectroscopic properties, such as charge density distribution and harmonic vibrational frequencies.


2013 ◽  
Vol 12 (01) ◽  
pp. 1250101 ◽  
Author(s):  
YA-NA WENG ◽  
XIAO-JUAN YAN ◽  
SHU-JIN LI

The mechanism of the reaction CF3CHFO2 + NO was investigated using ab initio and density functional theory (DFT). The optimized geometries for all stationary points on the reaction energy surface were calculated using MP2 and B3LYP methods with the aug-cc-pVDZ basis set. Single-point energy calculations were performed using the coupled cluster method with single, double and perturbative triple configurations, CCSD(T). The most important energy minima on the potential energy surface (PES) were found corresponding to two conformers of the peroxynitrite association adducts, cis- CF3CHFOONO and trans- CF3CHFOONO , and the nitrate, CF3CHFONO2 . The radical pairs ( CF3CHFO + NO2 ) and the nitrate are formed through the breaking of the peroxy bond of trans- CF3CHFOONO and the rearrangement of cis- CF3CHFOONO , respectively. The nitrate can be decomposed to carbonylated species ( CF3CHO or CF3CFO ), nitryl fluoride (NO2F), nitrous acid (HONO), and radical pairs ( CF3CHFO + NO2 ), which are of potential atmospheric importance.


Author(s):  
O.A. Nurkenov ◽  
L.K. Abulyaissova ◽  
G.Sh. Zhaksybayeva

Structural and spectroscopic properties of quinolizidine alkaloids lupinine and epilupinine stereoisomers were studied theoretically. The influence of the calculation method and structural change in the molecule on the results of geometry and other properties of compounds was considered. The equilibrium geometry, harmonic vibrational frequencies and infrared intensities were obtained by means of density functional theory (DFT/B3LYP) calculations with the splitvalence medium-sized 6-31G(d) basis and Dunning’s correlation consistent basis set cc-pVDZ. From the optimized structure of the (+)-lupinine and (+)-epilupinine molecules, geometric parameters were compared with the literature X-ray experimental data. Structural and vibrational parameters for the (-)-lupinine and (-)-epilupinine molecules are predicted by theoretical calculations at B3LYP/6-31G(d) and B3LYP/cc-pVDZ levels of theory. Some physical characteristics for the title compounds, such as total electronic energy, zero-point energy, rotational constants and dipole moments were also defined by DFT methods. The thermodynamic functions of the title compounds were performed at the same theory levels. Stationary points are identified by the solution of the oscillatory problem.


2015 ◽  
Vol 8 (3) ◽  
pp. 2197-2221
Author(s):  
Theraviyum Chithambarathanu ◽  
M. Darathi ◽  
J. DaisyMagdaline ◽  
S. Gunasekaran

The molecular vibrations of Trichloro isocyanuric acid (C3Cl3N3O3) and Trithio cyanuric acid (C3H3N3S3) have been investigated in polycrystalline sample at room temperature by Fourier Transform Infrared (FT-IR) and FT-Raman spectroscopies in the region 4000-450 cm-1 and 4000-50 cm-1 respectively, which provide a wealth of structural information about the molecules. The spectra are interpreted with the aid of normal co-ordinate analysis following full structure optimization and force field calculations based on density functional theory   (DFT) using standard B3LYP / 6-311++ G (d, p) basis set for investigating the structural and spectroscopic properties. The vibrational frequencies are calculated and the scaled values are compared with experimental FT-IR and FT-Raman spectra. The scaled theoretical wave numbers shows very good agreement with experimental ones. The complete vibrational assignments are performed on the basis of potential energy distribution (PED) of vibrational modes, calculated with scaled quantum (SQM) method. Stability of the molecule arising from hyper conjugative interactions, charge delocalization has been analyzed using natural bond orbital (NBO) analysis. The results show that change in electron density (ED) in σ* and π* anti-bonding orbitals and second order delocalization   energy (E2) confirm the occurrence of Intra molecular Charge Transfer (ICT) within the molecule. The thermodynamic properties like heat capacity, entropy, enthalpy and zero point energy have been calculated for the molecule. The frontier molecular orbitals have been visualized and the HOMO-LUMO energy gap has been calculated. The Molecular Electrostatic Potential (MEP) analysis reveals the sites for electrophilic attack and nucleophilic reactions in the molecule.


2015 ◽  
Vol 8 (2) ◽  
pp. 2122-2134
Author(s):  
Sarvendra Kumar ◽  
Rajesh Kumar ◽  
Jayant Teotia ◽  
M. K. Yadav

In the present work, UV- Visible spectra of 2-Chloro-3,4-Dimethoxybenzaldehyde (2,3,4-CDMB) compound  have been carried out experimentally and theoretically. The ultraviolet absorption spectrum of title compound in three solvents (Acetone, Diethyl Ether, CCl4) of different polarity were examined in the range of 200–500 nm. The structure of the molecule was optimized and the structural characteristics were determined by HF and DFT (B3LYP) methods with 6-31+G(d,p) and 6-311++G(d,p) as basis sets. The excitation energy, wavelength corresponds to absorption maxima () and oscillator strength (f) are calculated by Time-Dependent Density Functional Theory (TD-DFT) using B3LYP/6-31+G(d,p) and B3LYP/6-311++G(d,p) as basis sets. The electric dipole moment (μ), polarizability (α) and the first hyperpolarizability (β ) have been computed to evaluate the non-linear optical (NLO) response of the investigated compound by HF and DFT (B3LYP) with already mentioned basis sets. Thermodynamic functions of the title compound at different temperatures were also calculated.


2018 ◽  
Author(s):  
Oscar A. Douglas-Gallardo ◽  
David A. Sáez ◽  
Stefan Vogt-Geisse ◽  
Esteban Vöhringer-Martinez

<div><div><div><p>Carboxylation reactions represent a very special class of chemical reactions that is characterized by the presence of a carbon dioxide (CO2) molecule as reactive species within its global chemical equation. These reactions work as fundamental gear to accomplish the CO2 fixation and thus to build up more complex molecules through different technological and biochemical processes. In this context, a correct description of the CO2 electronic structure turns out to be crucial to study the chemical and electronic properties associated with this kind of reactions. Here, a sys- tematic study of CO2 electronic structure and its contribution to different carboxylation reaction electronic energies has been carried out by means of several high-level ab-initio post-Hartree Fock (post-HF) and Density Functional Theory (DFT) calculations for a set of biochemistry and inorganic systems. We have found that for a correct description of the CO2 electronic correlation energy it is necessary to include post-CCSD(T) contributions (beyond the gold standard). These high-order excitations are required to properly describe the interactions of the four π-electrons as- sociated with the two degenerated π-molecular orbitals of the CO2 molecule. Likewise, our results show that in some reactions it is possible to obtain accurate reaction electronic energy values with computationally less demanding methods when the error in the electronic correlation energy com- pensates between reactants and products. Furthermore, the provided post-HF reference values allowed to validate different DFT exchange-correlation functionals combined with different basis sets for chemical reactions that are relevant in biochemical CO2 fixing enzymes.</p></div></div></div>


Sign in / Sign up

Export Citation Format

Share Document