Exploring the stereodynamics and microscopic mechanism of the O(3P) + CH4, CD4→ OH + CH3, OD + CD3 combustion reactions

Potential energy surfaces for singlet and triplet states of the LiH2+ system and quasi-classical trajectory cross sections for H + LiH+ and H+ + LiH

A new set of six accurate ab initio potential energy surfaces (PESs) is presented for the first three singlet and triplet states of LiH2+ (1,21A', 11A'', 1,23A', and 13A'' states, where four of them are investigated for the first time), which have allowed new detailed studies gaining a global view on this interesting system. These states are relevant for the study of the most important reactions of lithium chemistry in the early universe. More than 45 000 energy points were calculated using the multi-reference configuration interaction level of theory using explicitly correlated methods (ic-MRCI-F12), and the results obtained for each individual electronic state were fitted to an analytical function. Using quasiclassical trajectories and considering the initial diatomic fragment in the ground rovibrational state, we have determined the integral cross sections for the H + LiH+(X2Σ+, C2Π) and H+ + LiH(X1Σ+, B1Π) reactions. In these calculations all available reaction channels were considered: the chemically most important H or H+ transfer/abstraction as well as atom exchange and collision induced dissociation for up to 1.0 eV of collision energy.

Mechanism and Kinetics of Methane Combustion. Part II: Potential Energy Surface for Hydrogen-Abstraction Reaction of CH4 + O(3P)

Methane combustion plays an important role in various fields such as combustion chemistry and atmospheric chemistry of the stratosphere. Highly accurate study of its initial reaction remains a key challenge. Here, through extensive studies with a state-of-the-art ab initio and neural network method, we present a potential energy surface of the O(³P) + CH₄ → OH + CH₃ reaction on the ground state 1³A and the first excited state 2³A. In this work, the energies of 10 167 points covering all important regions are obtained with state-averaged complete active space self-consistent field calculations and then fitted using the Levenberg-Marquardt algorithm with a root-mean-square error of 0.391 and 0.442 kcal/mol for the 1³A and 2³A states, respectively. This study explores the characteristics of the radical van der Waals (VdW) complex and reveals a detailed mechanism of the methane combustion initial reaction. Within the scope of this mechanism, this surface gives a fairly accurate description of the regions around the saddle point, conical intersection, and vdW wells in the entrance for efficient computational simulations. As a theoretical study on a prototypical polyatomic reaction, it is hopeful that this work will modify our understanding of the primary process in hydrocarbon combustion.

Surface Structure of Alkyl/Fluoroalkylimidazolium Ionic-Liquid Mixtures

The gas-liquid interface of ionic liquids (ILs) is critically important in many applications, for example, in supported IL phase (SILP) catalysis. Methods to investigate the interfacial structure in these systems will allow their performance to be improved in a rational way. In this study, reactive-atom scattering (RAS), surface tension measurements, and molecular dynamics (MD) simulations were used to study the vacuum interface of mixtures of partially fluorinated and normal alkyl ILs. The underlying aim was to understand whether fluorinated IL ions could be used as additives to modify the surface structure of one of the most widely used families of alkyl ILs. The series of ILs 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([C_nmim][Tf₂N]) with n = 4-12 were mixed with a fixed-length, semiperfluorinated analogue (1H,1H,2H,2H-perfluorooctyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([C₈mimF₁₃][Tf₂N]), forming [C_nmim]_(1-x)[C₈mimF₁₃]_x[Tf₂N] mixtures, where x is the bulk mole fraction of the fluorinated component. The RAS-LIF method combined O-atom projectiles with laser-induced fluorescence (LIF) detection of the product OH as a measure of surface exposure of the alkyl chains. For [C₈mim]_(1-x)[C₈mimF₁₃]_x[Tf₂N] mixtures, RAS-LIF OH yields are below those expected from stoichiometry. There are quantitatively consistent negative deviations from linearity of the surface tension. Both results imply that the lower-surface-tension fluoroalkyl material dominates the surface. A similar deficit is found for alkyl chain lengths n = 4, 6, 8, and 12 and for all (nonzero) x investigated by RAS-LIF. Accessible-surface-area (ASA) analyses of the MD simulations for [C_nmim]_(1-x)[C₈mimF₁₃]_x[Tf₂N] mixtures qualitatively reproduce the same primary effect of fluoro-chain predominance of the surface over most of the range of n. However, there are significant quantitative discrepancies between MD ASA predictions and experiment relating to the strength of any n-dependence of the relative alkyl coverage at fixed x, and on the x-dependence at fixed n. These discrepancies are discussed in the context of detailed examinations of the surface structures predicted in the MD simulations. Potential explanations, beyond experimental artifacts, include inadequacies in the classical force fields used in the MD simulations or the inability of simple ASA algorithms to capture dynamical factors that influence RAS-LIF yields.

Thermal Rate Constants for the O(3P) + CH4 → OH + CH3 Reaction: The Effects of Quantum Tunneling and Potential Energy Barrier Shape

The rate constants and kinetic isotope effects for the O(³P) + CH₄ reaction have been investigated with the quantum instanton method in full dimensionality. The calculated rate constants are in good agreement with the experimental values above 400 K, below which the measured values are scattered. Compared to other theoretical approaches, the quantum instanton method predicts the largest quantum tunneling effect, so it gives the largest rate constants at low temperatures. The calculated kinetic isotope effects are always much larger than 1 and increase with decreasing temperature, due to the zero-point energy and quantum tunneling. Our calculations on different potential energy surfaces demonstrate that the potential energy barrier shape dominates the magnitude of quantum tunneling and has a great effect on the kinetic isotope effect.