scholarly journals Effect of ammonia and water molecule on OH + CH3OH reaction under tropospheric condition

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Mohamad Akbar Ali ◽  
M. Balaganesh ◽  
Faisal A. Al-Odail ◽  
K. C. Lin

AbstractThe rate coefficients for OH + CH3OH and OH + CH3OH (+ X) (X = NH3, H2O) reactions were calculated using microcanonical, and canonical variational transition state theory (CVT) between 200 and 400 K based on potential energy surface constructed using CCSD(T)//M06-2X/6-311++G(3df,3pd). The results show that OH + CH3OH is dominated by the hydrogen atoms abstraction from CH3 position in both free and ammonia/water catalyzed ones. This result is in consistent with previous experimental and theoretical studies. The calculated rate coefficient for the OH + CH3OH (8.8 × 10−13 cm3 molecule−1 s−1), for OH + CH3OH (+ NH3) [1.9 × 10−21 cm3 molecule−1 s−1] and for OH + CH3OH (+ H2O) [8.1 × 10−16 cm3 molecule−1 s−1] at 300 K. The rate coefficient is at least 8 order magnitude [for OH + CH3OH(+ NH3) reaction] and 3 orders magnitude [OH + CH3OH (+ H2O)] are smaller than free OH + CH3OH reaction. Our calculations predict that the catalytic effect of single ammonia and water molecule on OH + CH3OH reaction has no effect under tropospheric conditions because the dominated ammonia and water-assisted reaction depends on ammonia and water concentration, respectively. As a result, the total effective reaction rate coefficients are smaller. The current study provides a comprehensive example of how basic and neutral catalysts effect the most important atmospheric prototype alcohol reactions.

2019 ◽  
Vol 21 (3) ◽  
pp. 1408-1416 ◽  
Author(s):  
Junxiang Zuo ◽  
Qixin Chen ◽  
Xixi Hu ◽  
Hua Guo ◽  
Daiqian Xie

A global potential energy surface for the O(3P) + C2H2reaction is developed and the quasi-classical trajectory study on the potential energy surface reproduce the rate coefficient and product branching ratio.


Gas-phase dissociation of fluorine ( 1 Ʃ + g ) molecules in an agron bath at 3000 K was studied by using the 3D Monte Carlo classical trajectory (3DMCCT) method. To assess the importance of the potential energy surface (PES) in such calculations, three surfaces, with a fixed, experimentally determined F 2 dissociation energy, were constructed. These surfaces span the existing experimental uncertainties in the shape of the F 2 potential. The first potential was the widest and softest; in the second potential the anharmonicity was minimized. The intermediate potential was constructed to ‘localize’ anharmonicity in the energy range in which the collisions are most reactive. The remaining parameters for each PES were estimated from the best available data on interatomic potentials. By using the single uniform ensemble (SUE) method (Kutz, H. D. & Burns, G. J. chem. Phys . 72, 3652-3657 (1980)), large ensembles of trajectories (LET) were generated for the PES. Two such ensembles consisted of 30000 trajectories each and the third of 26200. It was found that the computed one-way-flux equilibrium rate coefficients (Widom, B. Science 148, 1555-1560 (1965)) depend in a systematic way upon the anharmonicity of the potential, with the most anharmonic potential yielding the largest rate coefficient. Steady-state reaction-rate constants, which correspond to experimentally observable rate constants, were calculated by the SUE method. It was determined that this method yields (for a given trajectory ensemble, PES and translational temperature) a unique steady-state rate constant, independent of the initial, arbitrarily chosen, state (Tolman, R. C. The principles of statistical mechanics , p. 17. Oxford University Press (1938)) of the LET, and consequently independent of the corresponding initial value of the reaction rate coefficient. For each initial state of the LET, the development of the steady-state rate constant from the equilibrium rate coefficient was smooth, monotonic, and consistent with the detailed properties of the PES. It was found that, although the increased anharmonicity of the F 2 potential enhanced the equilibrium rate coefficients, it also enhanced the non-equilibrium effects. As a result, the steady-state rate constants were found to be insensitive to the variation of the PES. Thus, the differences among the steady-state rate constants for the three potentials were of the order of their standard errors, which was about 15% or less. On the other hand, the calculated rate constants exceeded the experimental rate constant by a factor of five to six. Because within the limitations of classical mechanics the calculations were ab initio , it was tentatively concluded that the discrepancy of five to six is due to the use of classical mechanics rather than details of the PES structure.


2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Mohamed A. Abdel-Rahman ◽  
Tarek M. El-Gogary ◽  
Nessreen Al-Hashimi ◽  
Mohamed F. Shibl ◽  
Kazunari Yoshizawa ◽  
...  

Abstract In this work, a theoretical investigation of thermochemistry and kinetics of the oxidation of bifunctional 2-Methoxyethanol (2ME) biofuel using methyl radical was introduced. Potential-energy surface for various channels for the oxidation of 2ME was studied at density function theory (M06-2X) and ab initio CBS-QB3 levels of theory. H-atom abstraction reactions, which are essential processes occurring in the initial stages of the combustion or oxidation of organic compounds, from different sites of 2ME were examined. A similar study was conducted for the isoelectronic n-butanol to highlight the consequences of replacing the ϒ CH2 group by an oxygen atom on the thermodynamic and kinetic parameters of the oxidation processes. Rate coefficients were calculated from the transition state theory. Our calculations show that energy barriers for n-butanol oxidation increase in the order of α ‹ O ‹ ϒ ‹ β ‹ ξ, which are consistent with previous data. However, for 2ME the energy barriers increase in the order α ‹ β ‹ ξ ‹ O. At elevated temperatures, a slightly high total abstraction rate is observed for the bifunctional 2ME (4 abstraction positions) over n-butanol (5 abstraction positions).


2018 ◽  
Vol 478 (3) ◽  
pp. 3952-3960 ◽  
Author(s):  
Andrey K Belyaev ◽  
Svetlana A Yakovleva

ABSTRACT Inelastic rate coefficients for 686 partial processes in low-energy Ba + H, Ba+ + H−, Ba++ H and Ba2+ + H− collisions are calculated. These data are needed for the non-local thermodynamic equilibrium (non-LTE) modelling of Ba i and Ba ii spectra, especially in cool stellar atmospheres. The calculations of the rate coefficients are performed by means of the quantum model approach, based on the asymptotic semi-empirical method for the electronic structure calculations and on multichannel formulas for the non-adiabatic nuclear dynamical calculations. The inelastic rate coefficients for all transitions between the 17 lowest covalent states and one ionic molecular state in Ba + H and Ba+ + H− collisions, as well as the inelastic rate coefficients for all transitions between the 19 lowest covalent states and one ionic molecular state in Ba+ + H and Ba2+ + H− collisions are calculated. In Ba+ + H− collisions, the highest rate coefficients correspond to the mutual neutralization processes into the   Ba(6s6p1P°), Ba(6s7s3S) and   Ba(6s7s1S) final states, with the largest value of 5.93 × 10−8 cm3 s−1 at T = 6000 K for the process Ba+ + H− →   Ba(6s7s3S) + H. The highest rate coefficient for excitation and de-excitation processes in Ba + H collisions corresponds to the   Ba(6s7s1S) →  Ba(6s7s3S) transition, with the value of 7.62 × 10−9 cm3 s−1 at T = 6000 K. In Ba2+ + H− collisions, the highest rate coefficients correspond to the neutralization processes into the Ba+( 7p2P°), Ba+( 4f 2F°), Ba+( 6d 2D) and Ba+( 7s 2S) final states. The highest neutralization rate has the value of 3.96 × 10−8 cm3 s−1 at T = 6000 K for the Ba2+ + H− → Ba+( 7p 2P°) + H process. The largest rate coefficient for excitation and de-excitation processes in Ba+ + H collisions corresponds to the Ba+(7s 2S) → Ba+( 6p 2P°) transition, with the value of 1.23 × 10−9  cm3 s−1 at T = 6000 K.


2020 ◽  
Vol 20 (12) ◽  
pp. 7429-7458 ◽  
Author(s):  
Luc Vereecken ◽  
Barbara Nozière

Abstract. A large data set of rate coefficients for H migration in peroxy radicals is presented and supplemented with literature data to derive a structure–activity relationship (SAR) for the title reaction class. The SAR supports aliphatic RO2 radicals; unsaturated bonds and β-oxo substitutions both endocyclic and exocyclic to the transition state ring; and α-oxo (aldehyde), –OH, –OOH, and –ONO2 substitutions, including migration of O-based hydrogen atoms. Also discussed are –C(=O)OH and –OR substitutions. The SAR allows predictions of rate coefficients k(T) for a temperature range of 200 to 450 K, with migrations spans ranging from 1,4 to 1,9-H shifts depending on the functionalities. The performance of the SAR reflects the uncertainty of the underlying data, reproducing the scarce experimental data on average to a factor of 2 and the wide range of theoretical data to a factor of 10 to 100, depending also on the quality of the data. The SAR evaluation discusses the performance in multi-functionalized species. For aliphatic RO2, we also present some experimental product identification that validates the expected mechanisms. The proposed SAR is a valuable tool for mechanism development and experimental design and guides future theoretical work, which should allow for rapid improvements of the SAR in the future. Relative multi-conformer transition state theory (rel-MC-TST) kinetic theory is introduced as an aid for systematic kinetic studies.


Author(s):  
Francisco J. Valdés-Parada ◽  
Didier Lasseux ◽  
Stephen Whitaker

Abstract Diffusion and reaction in porous media have been studied extensively due to the wide range of applications in which this transport phenomenon is involved. In particular, in chemical reactor engineering, reactive mass transfer is crucial to understand the performance of porous catalyst particles immersed in chemical reactors. Due to the disparity of characteristic lengths between the pores and the porous particles, this type of process is usually modeled by means of effective-medium equations, in which the solid and fluid phases are conceived as a pseudo-continuum. For conditions in which the pore-scale Thiele modulus (or Kinetic number) is much smaller than unity, it is reasonable to assume that the effective diffusivity involved in the effective-medium model is only a function of the porous medium geometry. However, a long debate has existed in the literature concerning the extensive use of this assumption for situations in which the Kinetic Number does not satisfy the above mentioned constraint. In addition, the functionality of the effective reaction rate coefficient with the Kinetic number has not been sufficiently studied. In this work we address these issues by means of the volume averaging method. Our analysis is focused on cases in which the Kinetic number can reach values up to 1. Interestingly, for this particular condition, the use of the intrinsic diffusivity tensor is justified. In addition, by means of Maclaurin series expansions, the effective reaction rate coefficient is shown to be acceptably approximated as a first-order function. These two conclusions for the effective medium coefficients constitute the major contributions from this work. In addition, the predictions from the upscaled model are validated by comparison with direct numerical simulations under steady and transient conditions.


2020 ◽  
Vol 494 (4) ◽  
pp. 5675-5681 ◽  
Author(s):  
Sanchit Chhabra ◽  
T J Dhilip Kumar

ABSTRACT Molecular ions play an important role in the astrochemistry of interstellar and circumstellar media. C3H+ has been identified in the interstellar medium recently. A new potential energy surface of the C3H+–He van der Waals complex is computed using the ab initio explicitly correlated coupled cluster with the single, double and perturbative triple excitation [CCSD(T)-F12] method and the augmented correlation consistent polarized valence triple zeta (aug-cc-pVTZ) basis set. The potential presents a well of 174.6 cm−1 in linear geometry towards the H end. Calculations of pure rotational excitation cross-sections of C3H+ by He are carried out using the exact quantum mechanical close-coupling approach. Cross-sections for transitions among the rotational levels of C3H+ are computed for energies up to 600 cm−1. The cross-sections are used to obtain the collisional rate coefficients for temperatures T ≤ 100 K. Along with laboratory experiments, the results obtained in this work may be very useful for astrophysical applications to understand hydrocarbon chemistry.


Author(s):  
Daniel C Morris ◽  
Stuart W Prescott ◽  
Jason B Harper

A series of ionic liquids based on the 1-alkyl-3-methylimidazolium cations were examined as components of the solvent mixture for a bimolecular substitution process. The effects on both the rate coefficient...


2017 ◽  
Author(s):  
Siripina Vijayakumar ◽  
Avinash Kumar ◽  
Balla Rajakuma

Abstract. Temperature dependent rate coefficients for the gas phase reactions of Cl atoms with 4-hexen-3-one and 5-hexen-2-one were measured over the temperature range of 298–363 K relative to 1-pentene, 1,3-butadiene and isoprene. Gas Chromatography (GC) was used to measure the concentrations of the organics. The derived temperature dependent Arrhenius expressions are k4-hexen-3-one+Cl (298–363 K) = (2.82 ± 1.76)×10−12exp [(1556 ± 438)/T] cm3 molecule−1 s−1 and k5-hexen-2-one+Cl (298–363 K) = (4.6 ± 2.4)×10−11exp[(646 ± 171)/T] cm3 molecule−1 s−1. The corresponding room temperature rate coefficients are (5.54 ± 0.41)×10−10 cm3 molecule−1 s−1 and (4.00 ± 0.37)×10−10 cm3 molecule−1 s−1 for the reactions of Cl atoms with 4-hexen-3-one and 5-hexen-2-one respectively. To understand the mechanism of Cl atom reactions with unsaturated ketones, computational calculations were performed for the reactions of Cl atoms with 4-hexen-3-one, 5-hexen-2-one and 3-penten-2-one over the temperature range of 275–400 K using Canonical Variational Transition state theory (CVT) with Small Curvature Tunneling (SCT) in combination with CCSD(T)/6-31+G(d, p)//MP2/6-311++G(d, p) level of theory. Atmospheric implications, reaction mechanism and feasibility of the title reactions are discussed in this manuscript.


2020 ◽  
Author(s):  
Zoi Salta ◽  
Jacopo Lupi ◽  
Vincenzo Barone ◽  
Oscar Ventura

<div> Elucidation of the oxidation mechanism of naturally emitted reduced sulfur compounds, especially dimethyl sulfide, plays a central role in understanding background acid precipitation in the natural environment. Most frequently, theoretical studies of the addition and H-elimination reactions of dimethyl sulfide with hydroxyl radicals are studied considering the presence of oxygen that further reacts with the radicals formed in the initial steps. Although the reaction of intermediate species with additional hydroxyl radicals has been considered as part of the global mechanism of oxidation, few if any attention has been dedicated to the possibility of reactions of the initial radicals with a second •OH molecule. In this work we performed a computational study using quantum-chemical methods, of the mechanism of H-abstraction from dimethyl sulfide under normal atmospheric conditions and in reaction chambers at different O2 partial pressure, including complete absence of oxygen. Additionally, important rate coefficients were computed using canonical and variational transition state theory. The rate coefficient for abstraction affords a 4.72 x 10-12 cm3 molecule1 s-1 value, very close to the most recent experimental one (4.13 x 10-12 cm3 molecule-1 s-1). According to our best results, the initial methyl thiomethyl radical was obtained at -25.2 kcal/mol (experimentally -22.4 kcal/mol), and four important paths were identified on the potential energy surface. From the interplay of thermochemical and kinetic arguments, it was possible to demonstrate that the preferred product of the reaction of dimethyl sulfide with two hydroxyl radicals, is actually dimethyl sulfoxide. </div><div> </div>


Sign in / Sign up

Export Citation Format

Share Document