Variational Transition State Theory Research Articles

Comprehensive Multipath Variational Kinetics Study on Hydrogen Abstraction Reactions from Three Typical Dimethylcyclohexane Isomers by Hydroxyl Radicals: from the Electronic Structure to Model Applications.

Cycloalkanes serve as an important class of chemical components in both fossil and alternative transportation fuels and have attracted considerable attention from the combustion community. Hydrogen abstractions from cycloalkanes by hydroxyl radicals initiate the fuel decomposition process and trigger off the subsequent chain reactions and thus play an important role in both combustion and atmospheric chemistry. The target of this study is to fill the vacancy in kinetics data toward the H-abstraction reactions by hydroxyl radical from three typical dimethylcyclohexane isomers through first-principles and direct dynamics. The rate constants involving 18 elementary reactions in total were accurately determined by the multipath canonical variational transition state theory with the multidimensional small-curvature correction for tunneling (MP-CVT/SCT), over a broad temperature range of 200-2000 K. The significant roles of multistructural torsional anharmonicity and recrossing effects were stressed per abstraction site, while the quantum tunneling effect was found to be slight at temperatures of interest in combustion. The discrepancies observed among different reaction systems at a similar abstraction site highlight the fuel molecular effects on site-specific rate constants. The comparison results of total rate constants given by different dynamics approaches prove the importance of considering the torsional anharmonicity, recrossing, and tunneling effects, and the robust feature of the simplified MS-CVT/SCT. The calculated total constants for dimethylcyclohexane isomers by OH are consistent with those measured for methylcyclohexane and 1,4-dimethylcyclohexane at low temperatures. The branching ratio analysis confirms the predominant role of the tertiary abstraction at low-to-intermediate temperatures and its growing competition with distinct secondary abstractions as temperature increases. The calculated rate constants were eventually fitted into the analytical expressions and incorporated into the kinetic models to learn about the influences on modeling performance.

Dual-level strategy for quantitative kinetics for the reaction between ethylene and hydroxyl radical.

The atmospheric reactions are mainly initiated by hydroxyl radical (OH). Here, we choose the C2H4 + OH reaction as a model reaction for other reactions of OH with alkenes. We use the GMM(P).L//CCSD(T)-F12a/cc-pVTZ-F12 theoretical method as the benchmark results close to the approximation of CCSDTQ(P)/CBS accuracy to investigate the C2H4 + OH reaction. The rate constants for the C2H4 + OH reaction at high-pressure limit were calculated by using the dual-level strategy. It integrates the transition state theory rate constant calculated by GMM(P).L//CCSD(T)-F12a/cc-pVTZ-F12 with the canonical variational transition state theory containing small-curvature tunneling (CVT/SCT) calculated by using the M11-L functional method with the MG3S basis set. The rate constants of C2H4 + OH at different pressures were obtained by using both the system-specific quantum Rice-Ramsperger-Kassel (SS-QRRK) theory and master equation method. The calculated results uncover that both the calculated rate constants at different pressures and temperatures are quantitatively consistent with the values obtained by the experimental measurements in the C2H4 + OH reaction. We find that the post-CCSD(T) contributions to the barrier height for the C2H4 + OH reaction are significant with the calculated value of -0.38 kcal/mol. We also find that the rate determining step is only dominated by the tight transition state under atmospheric conditions, whereas previous investigations indicated that the rate constants were controlled by both the loose and tight transition states in the C2H4 + OH reaction. The present findings unravel that it is an important factor for the effect of torsional anharmonicity on quantitative kinetics.

Thermal rate constants and kinetic isotope effects of the H + H2O2 reactions: barrier height and reaction energy from single- and multireference methods.

One of the more significant sub-mechanisms of H2/O2 combustion involves the reaction of hydrogen peroxide with hydrogen atoms (H + H2O2), resulting in the production of OH + H2O (R1) and H2 + HO2 (R2) paths. Previous experimental and ab initio calculations reveal some variations in the barrier height for (R1). To improve the energetics of both (R1) and (R2), single reference and multireference ab initio methods are employed, and the rate constants and H/D kinetic isotope effects (KIEs) are calculated as a function of temperature. For (R1), the best results for the barrier height and reaction energies computed with the CASPT2(15,11)/aug-cc-pV6Z are 5.2 and - 70.3kcal.mol-1, respectively. CCSD(T)/aug-cc-pV5Z + CV (core-valence) calculations for (R2) give 9.7 and - 15.6kcal.mol-1 to those parameters. The CVT/SCT rate constants of both paths agree well with the fitted rate constants from uncertainty-weighted statistical analysis of the 14-mechanism of H2/O2. The kinetic isotopic effect (kH/kD) for the reaction D + H2O2 → DH + HO2 was found to be 0.47, which is in excellent agreement with the experimental value of 0.43. The structures of reactants, transition state, and products of (R1) and (R2) are calculated with the aug-cc-pVTZ basis set and M062X DFT, CCSD(T), and CASSCF methods. The barrier heights and reaction energies of (R1) and (R2) are computed using the M06-2X, CCSD(T), MRCI, and CASPT2 methods and various basis sets. The rate constants are calculated with the variational transition state theory including multidimensional tunneling corrections (VTST-MT), with potential energy surfaces built by the M06-2X/aug-cc-pVTZ approach.

Unraveling the Catalytic Mechanism of β-Cyclodextrin in the Vitamin D Formation.

Previous experimental studies have shown that the isomerization reaction of previtamin D3 (PreD3) to vitamin D3 (VitD3) is accelerated 40-fold when it takes place within a β-cyclodextrin dimer, in comparison to the reaction occurring in conventional isotropic solutions. In this study, we employ quantum mechanics-based molecular dynamics (MD) simulations and statistical multistructural variational transition state theory to unveil the origin of this acceleration. We find that the conformational landscape in the PreD3 isomerization is highly dependent on whether the system is encapsulated. In isotropic media, the triene moiety of the PreD3 exhibits a rich torsional flexibility. However, when encapsulated, such a flexibility is limited to a more confined conformational space. In both scenarios, our calculated rate constants are in close agreement with experimental results and allow us to identify the PreD3 flexibility restriction as the primary catalytic factor. These findings enhance our understanding of VitD3 isomerization and underscore the significance of MD and environmental factors in biochemical modeling.

Comparative Study on the H-Abstraction Reactions of Isopropyl Acetate and Propyl Acetate by HO2 and OH Radicals.

Isopropyl acetate (IPA) and propyl acetate (PA) are recognized as promising biofuels suitable for applications as fuel additives and biodiesel models. The H-abstraction reactions with radicals stand out as the fundamental initiating reactions in the combustion kinetic models for IPA and PA. In the present work, the kinetic calculations of IPA and PA plus HO2 and OH radicals were investigated at M06-2X/cc-pVTZ//G4, M08-HX/maug-cc-pVTZ, and CCSD(T)/jul-cc-pVTZ levels. The thermodynamic calculations were obtained based on the G4 and CBS-APNO methods. Rate coefficients were calculated using both transition state theory and canonical variational transition state theory with tunneling correction at the temperature range of 250-2000 K. The total rate constants for the IPA + OH system were fitted as follows: k = 0.4674 × T3.927 exp(2128/T) (cm3 mol-1 s-1), and for the PA + OH system, the total rate constants were determined using the following equation: k = 0.0161 × T4.373 exp(2220/T) (cm3 mol-1 s-1). The rate coefficients of IPA + OH reactions determined based on the M08-HX/maug-cc-pVTZ level effectively replicate the experimental data, while H-abstraction rate coefficients of PA + OH by the CCSD(T)/jul-cc-pVTZ method accurately reproduce the experimental data. Refining the H-abstraction rate coefficients in the kinetic mechanism of PA, as proposed by Dayma et al. [Proc. Combust. Inst. 37 (2019) 429-436], has been achieved through incorporating the present calculated data, leading to the development of a revised mechanism. The validation of the updated mechanism against jet-stirred reactor data is presented, showcasing its effective performance in predicting JSR data.

Theoretical Study on Ethylamine Dissociation Reactions Using VRC-VTST and SS-QRRK Methods.

Barrierless bond dissociation reactions play an important role in fuel combustion. In this work, the pressure-dependent dissociation rate constants of ethylamine (EA) are accurately determined using variable-reaction-coordinate variational transition-state theory combined with the system-specific quantum Rice-Ramsperger-Kassel method. Before the kinetics calculations, the performances of four density functional theory methods in describing the bond dissociation of EA are evaluated against the benchmark method, FIC-MRCISD(T)+Q/cc-pVTZ, and the MN15-L/cc-pVTZ method is the best choice. By comparison of the Gibbs free energies and the rate constants for the bond dissociation reactions of EA, ethanol, and propane, the influence of functional groups on the reaction kinetics is discussed. The kinetics calculations show that the dissociation rate constants of EA are sensitive to pressure at low pressures and high temperatures, and the dominant channel is the reaction that yields C2H5 and NH2 radicals. A literature combustion model of EA is updated with our calculations, and the satisfactory agreement between the model predictions and reported ignition delay times of EA suggests the reliability of our calculations.

First-principles investigation of hydrogen-related reactions on (100)–(2 × 1)[formula omitted]H diamond surfaces

Hydrogen radical attacks and subsequent hydrogen migrations are considered to play an important role in the atomic-scale mechanisms of diamond chemical vapour deposition growth. We perform a comprehensive analysis of the reactions involving H-radical and vacancies on H-passivated diamond surfaces exposed to hydrogen radical-rich atmosphere. By means of first principles calculations—density functional theory and climbing image nudged elastic band method—transition states related to these mechanisms are identified and characterised. In addition, accurate reaction rates are computed using variational transition state theory. Together, these methods provide—for a broad range of temperatures and hydrogen radical concentrations—a picture of the relative likelihood of the migration or radical attack processes, along with a statistical description of the hydrogen coverage fraction of the (100) H-passivated surface, refining earlier results via a more thorough analysis of the processes at stake. Additionally, the migration of H-vacancy is shown to be anisotropic, and occurring preferentially across the dimer rows of the reconstructed surface. The approach used in this work can be generalised to other crystallographic orientations of diamond surfaces or other semiconductors.

Rapid Atmospheric Reactions between Criegee Intermediates and Hypochlorous Acid.

Hypochlorous acid (HOCl) is a paramount compound in the atmosphere due to its significant contribution to both tropospheric oxidation capacity and ozone depletion. The main removal routes for HOCl are photolysis and the reaction with OH during the daytime, while these processes are unimportant during the nighttime. Here, we report the rapid reactions of Criegee intermediates (CH2OO and anti/syn-CH3CHOO) with HOCl by using high-level quantum chemical methods as the benchmark; their accuracy is close to coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis limit by extrapolation [CCSDT(Q)/CBS]. Their rate constants have been calculated by using a dual-level strategy; this combines conventional transition state theory calculated at the benchmark level with variational transition state theory with small-curvature tunneling by a validated density functional method. We find that the introduction of the methyl group into Criegee intermediates not only affects their reactivities but also exerts a remarkable influence on anharmonicity. The calculated results uncover that anharmonicity increases the rate constants of CH2OO + HOCl by a factor of 18-5, while it is of minor importance in the anti/syn-CH3CHOO + HOCl reaction at 190-350 K. The present findings reveal that the loose transition state for anti-CH3CHOO and HOCl is a rate-determining step at 190-350 K. We also find that the reaction of Criegee intermediates with HOCl contributes significantly to the sink of HOCl during the nighttime in the atmosphere.

A theoretical and experimental study of 2-ethylfuran + OH reaction

2-Ethylfuran (EF2), having similar thermophysical properties as gasoline, can be produced from biomass and has received increasing attention in recent years. The reaction of EF2 with OH radical plays a fundamental role in combustion chemistry of furan fuels. This work provides a comprehensive investigation on the reaction kinetics of the title reaction. The performance of several density functional theory methods is evaluated against the DLPNOCCSD(T)/CBS(T-Q) benchmark method. The results show that M08-HX/jun-cc-pVTZ is the best one for both abstraction and addition reactions with the lowest average unsigned deviations of 0.61 and 0.63 kcal mol−1, respectively. Multi-structural variational transition state theory combined with small-curvature tunneling approximation (MS-CVT/SCT) is employed to calculate the high-pressure limit (HPL) rate constants over a wide temperature range of 500−2000 K. The system-specific quantum Rice-Ramsperger-Kassel (SS-QRRK) method is used to obtain the pressure-dependent rate constants for addition reactions at several selected pressures from 0.01 to 100 atm. We have measured the rate constants of EF2 + OH reaction behind reflected shock wave at temperatures of 985−1342 K and pressure around 1 bar. The measured rate constants exhibit a positive temperature dependence and may be described by a three-parameter Arrhenius expression: k=1.28×10−10×T6.89exp(5318T) (cm3 mol−1 s–1). The calculated results are in good agreement with measurements and the total rate constants exhibit a non-monotonic temperature dependence due to the competition between the abstraction and addition channels. At the HPL, the addition reaction is dominant, and the branching ratio is larger than abstraction reaction at 500–1800 K and is 0.92 at 500 K. Finally, to evaluate the performance of our calculations on the ignition delay time of EF2, three available kinetic models (Somers et al., Xu et al. and Li et al.) of EF2 are updated and the results show that the predictions of Li's model are improved.

Theoretical study on iso‐pentanol oxidation chemistry: Fuel radical isomerization and decomposition kinetics and mechanism development

AbstractThis study undertakes a detailed theoretical investigation into the iso‐pentanol radical isomerization and decomposition kinetics and the mechanism development of the iso‐pentanol oxidation. The CCSD(T)/CBS//M08‐HX/6‐311+G(2df,2p) method was adopted to calculate the reaction potential energy surface. The reaction rate coefficients were calculated by variational transition state theory (VTST) with multistructural torsional (MS‐T) partition function and small curvature tunneling (SCT) correction. Moreover, the pressure‐dependent rate coefficients were determined using the system‐specific quantum Rice‐Ramsperger‐Kassel theory (SS‐QRRK). The variational and tunneling effects were discussed, and the dominant reaction channels were identified. It reveals that the isomerization reactions play a significant role at low temperatures, while the decomposition reactions dominate the high‐temperature regime. Notably, the quantitative rate expressions for iso‐pentanol radical decomposition reactions were also obtained. Furthermore, a new kinetic model incorporating the calculated rate coefficients was constructed, exhibiting satisfactory prediction performance on ignition delay times and improved predictive accuracy of species mole fractions. This work provides accurate rate data of isomerization and decomposition kinetics and contributes to a more comprehensive understanding of the iso‐pentanol oxidation mechanism.

New insights into the mechanism and kinetics of the addition reaction of unsaturated Criegee intermediates to CF3COOH and tropospheric implications.

In this work, we have investigated the mechanism, thermochemistry and kinetics of the reaction of syn-cis-CH2RzCRyCO+O- (where Rz, Ry = H, CH3-) unsaturated Criegee intermediates (CIs) with CF3COOH using quantum chemical methods. The rate coefficients for the barrierless reactions were calculated using variable reaction coordinate variational transition state theory (VRC-VTST). For the syn-cis-CH2RzCRyCO+O- conformation in which conjugated CC and CO double bonds are aligned with each other, we propose a new pathway for the unidirectional addition of an OC-OH molecule (CF3COOH) to the CC double bond of syn-cis-CH2RzCRyCO+O-. The rate coefficient for the 1,4-CC addition reaction at 298 K is ∼10-10 to 10-11 cm3 s-1, resulting in the formation of CF3C(O)OCH2CRzRyCOOH trifluoroacetate alkyl allyl hydroperoxide (TFAAAH) as a new transitory adduct. It can act as a precursor for the formation of secondary organic aerosols (SOAs). This novel TFAAAH hydroperoxide was identified through a detailed quantum chemical study of the 1,4-addition mechanism and will provide new insights into the significance of the 1,4-addition reaction of unsaturated Cls with trace tropospheric gases on -CRzCH2 vinyl carbon atoms.

Toward an Accurate Black-Box Tool for the Kinetics of Gas-Phase Reactions Involving Barrier-less Elementary Steps.

An enhanced computational protocol has been devised for the accurate characterization of gas-phase barrier-less reactions in the framework of the reaction-path (RP) and variable reaction coordinate variational transition-state theory. In particular, the synergistic combination of density functional theory and Monte Carlo sampling to optimize reactive fluxes led to a reliable yet effective computational workflow. A black-box strategy has been developed for selecting the most suited density functional with reference to a high-level one-dimensional reference potential. At the same time, different descriptions of hindered rotations are automatically selected, depending on the corresponding harmonic frequencies along the RP. The performance of the new tool is investigated by means of two prototypical reactions involving different degrees of static and dynamic correlation, namely, H2S + Cl and CH3 + CH3. The remarkable agreement of the computed kinetic parameters with the available experimental data confirms the accuracy and robustness of the proposed approach. Together with their intrinsic interest, these results also pave the way toward systematic investigations of gas-phase reactions involving barrier-less elementary steps by a reliable, user-friendly tool, which can be confidently used also by nonspecialists.

Theoretical Kinetics Predictions for Reactions on the NH2O Potential Energy Surface.

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH2O potential energy surface, including the NH2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (∼0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H2 channel contributes 10% of the total NH2 + O flux at combustion temperatures, although this channel is not included in modern NH3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH3 oxidation and NOx formation. The present results for NH2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH3 and thermal DeNOx.

Kinetics of the Simplest Criegee Intermediate Reaction with Water Vapor: Revisit and Isotope Effect.

The kinetics of the simplest Criegee intermediate (CH2OO) reaction with water vapor was revisited. By improving the signal-to-noise ratio and the precision of water concentration, we found that the kinetics of CH2OO involves not only two water molecules but also one and three water molecules. Our experimental results suggest that the decay of CH2OO can be described as d[CH2OO]/dt = -kobs[CH2OO]; kobs = k0 + k1[water] + k2[water]2 + k3[water]3; k1 = (4.22 ± 0.48) × 10-16 cm3 s-1, k2 = (10.66 ± 0.83) × 10-33 cm6 s-1, k3 = (1.48 ± 0.17) × 10-50 cm9 s-1 at 298 K and 300 Torr with the respective Arrhenius activation energies of Ea1 = 1.8 ± 1.1 kcal mol-1, Ea2 = -11.1 ± 2.1 kcal mol-1, Ea3 = -17.4 ± 3.9 kcal mol-1. The contribution of the k3[water]3 term becomes less significant at higher temperatures around 345 K, but it is not ignorable at 298 K and lower temperatures. By quantifying the concentrations of H2O and D2O with a Coriolis-type direct mass flow sensor, the kinetic isotope effect (KIE) was investigated at 298 K and 300 Torr and KIE(k1) = k1(H2O)/k1(D2O) = 1.30 ± 0.32; similarly, KIE(k2) = 2.25 ± 0.44 and KIE(k3) = 0.99 ± 0.13. These mild KIE values are consistent with theoretical calculations based on the variational transition state theory, confirming that the title reaction has a broad and low barrier, and the reaction coordinate involves not only the motion of a hydrogen atom but also that of an oxygen atom. Comparing the results recorded under 300 Torr (N2 buffer gas) with those under 600 Torr, a weak pressure effect of k3 was found. From quantum chemistry calculations, we found that the CH2OO + 3H2O reaction is dominated by the reaction pathways involving a ring structure consisting of two water molecules, which facilitate the hydrogen atom transfer, while the third water molecule is hydrogen-bonded outside the ring. Furthermore, analysis based on dipole capture rates showed that the CH2OO(H2O) + (H2O)2 and CH2OO(H2O)2 + H2O pathways will dominate in the three water reaction.

Mechanistic study on the effect of ammonia co-firing with pulverized coal on NO formation and reduction

In this study, the heterogeneous formation and reduction mechanisms of NO by ammonia (NH3) were determined to evaluate the performance of NH3 co-firing with coal in NO emissions by density functional theory and canonical variational transition-state theory. NO formation is initialized by the dissociative adsorption of O2 on the char, and the subsequent dehydrogenation reactions after NH3 adsorption are prioritized. NO reduction via char synergism with NH3 are also elucidated by the NNH dissociation or ONNH conversion. Thermodynamic and kinetic results indicate that introducing NH3 has the advantage of promoting heterogeneous NO reduction, which is characteristic of the synergistic effect afforded by the properties of char and reducing agent NH3. The reduction effect of NH3 on NO is more significant than its oxidation effect for NO formation at elevated temperatures, which provides theoretical evidence for the potential advantages of NH3 co-firing on NO emissions in coal-fired boilers.

A density functional theory study of the molecular structure, reactivity, and spectroscopic properties of 2-(2-mercaptophenyl)-1-azaazulene tautomers and rotamers

Five stable tautomer and rotamers of the 2-(2-Mercaptophenyl)-1-azaazulene (thiol, thione, R1, R2, and R3) molecules were studied using density functional theory (DFT). The geometries of the studied tautomer and rotamers were fully optimized at the B3LYP/6-31G(d,p) level. Thermodynamic calculations were performed at M06-2X/6-311G++(2d,2p) and ωB97XD/6-311G++(2d,2p) in the gas phase and ethanol solution conditions modeled by the solvation model based on density (SMD). The kinetic constant of tautomer and rotamers conversion was calculated in the temperature range 270–320 K using variational transition state theory (VTST) accompanied by one-dimensional wigner tunneling correction. Energy refinement at CCSD(T)/6–311++G(2d,2p) in the gas phase has been calculated. All the studied DFT methods qualitatively give similar tautomer stability orders in the gas phase. The ethanol solvent causes some reordering of the relative stability of 2-(2-Mercaptophenyl)-1-azaazulene conformers. The transition states for the 2-(2-Mercaptophenyl)-1-azaazulene tautomerization and rotamerization processes were also determined. The reactivity, electric dipole moment, and spectroscopic properties of the studied tautomer and rotamers were computed. The hyper-Rayleigh scattering (βHRS), and depolarization ratio (DR) exhibited promising optical properties when nonlinear optical properties were calculated.

Polyrate 2023: A computer program for the calculation of chemical reaction rates for polyatomics. New version announcement

Polyrate 2023: A computer program for the calculation of chemical reaction rates for polyatomics. New version announcement

TUMME 2023: Tsinghua University Minnesota Master Equation program. New version announcement

This is an updated version of TUMME (Tsinghua University Minnesota Master Equation program), which is a computer program for setting up and solving master equations for chemical kinetics of unimolecular and bimolecular reactions. The master equation is discretized in terms of energy bins, and rate constants are calculated from chemically significant eigenmodes. TUMME has interfaces to Gaussian 16, Polyrate 2023, and MSTor 2023 output files that allow the master equation code to obtain the microcanonical flux coefficients needed for the kernel of the master equation as calculated by conventional transition state theory (TST), variational transition state theory (VTST) with various tunneling methods, or multi-structural or multi-path VTST (MS-VTST or MP-VTST) with various tunneling methods. The tunneling methods supported include zero-curvature tunneling (ZCT), small-curvature tunneling (SCT), large-curvature tunneling (LCT), and microcanonically optimized multidimensional tunneling (μOMT). For mechanisms involving only unimolecular isomerization (no bimolecular pairs), TUMME can set up and solve a conservative master equation for both rate constants and time-dependent energy-bin populations. For mechanisms involving bimolecular pairs, TUMME 2023 can set up and solve two kinds of master equation: (i) a nonconservative master equation for calculating rate constants of bimolecular reactions and (ii) a conservative master equation that includes bimolecular association in the transition matrix and that can be used for calculating the time evolution of the concentration of a pseudo-first-order bimolecular reactant. TUMME is written in double precision with Python 3; quadruple and octuple precision are also available for some subtasks in C++. The Python code can run in serial or parallel (MP or MPI), and C++ code can run on a single processor or on multiple processors with OpenMP. The program includes a manual and a tutorial. New version program summaryProgram Title: TUMME 2023CPC Library link to program files:https://doi.org/10.17632/whcnvm2mc9.2Developer's repository link:https://doi.org/10.5281/zenodo.7943283Licensing Provisions: Apache-2.0 for the program and CC-BY-4.0 for the manualProgramming languages: Python 3 and C++Does the new version supersede the previous version?: YesJournal reference of previous version: Comput. Phys. Commun. 270 (2022) 108140 [1]External libraries: Numpy, Scipy, Numba, mpi4py (optional), modified mpack (optional), qd (optional), omp (optional)Interfaces to other programs: Gaussian 16 [2], Polyrate 2023 [3], and/or MSTor 2023 [4]Nature of problem: Characterize and calculate rate constants for a temperature-dependent and pressure-dependent complex reaction systemSolution method: Solve the energy-grained master equation based on the chemically significant eigenmodes to get phenomenological rate constants and time evolutions of the populations [1,5–7].Reason for new version: The new version increases the functionality and accuracy, and it offers an easier way to run the program. In addition, some bugs have been corrected.Summary of revisions:1.Added an option (available in high precision) to treat pseudo-first-order bimolecular reactions to obtain the time evolution of the reactants and intermediates. [8]2.Added an inverse-Laplace-transform option in the REACTION block.3.Abandoned rounding the sizes of energy bins to an integer wavenumber (ESOT).4.Changed the way to run TUMME by providing an executable bash script.5.Slightly adjusted the strategy for eigenmode mergers to deal with some unexpected cases.6.Adjusted the output for the partition functions.7.Added options for large-curvature tunneling (LCT) and microcanonically optimized multidimensional tunneling (μOMT).8.Added the keyword pyrfile to support a Polyrate 2023 [3] file that ends with a suffix .fu100.9.Fixed some bugs.10.Updated the manual and added a tutorial to the program distribution, so the distribution now has both a manual and a tutorial.

Kinetics of the Reactions of Chlorinated Very Short-Lived Substances (VSLSs) with Chlorine Atoms.

Chlorinated very short-lived substances (VSLSs), which are not controlled by the Montreal Protocol, are of current concern with regard to recovery of stratospheric ozone. Further study is needed on the temperature dependences of chlorinated VSLSs relevant to atmospheric conditions. Here, the kinetics of chlorinated VSLSs, such as chloroform (CHCl3), dichloromethane (CH2Cl2), dichloroethane (CH2ClCH2Cl), and trichloroethene (C2HCl3) reacting with chlorine atoms, were investigated between 180 and 400 K, expanding the range of temperatures relative to previous studies. RRKM/Master Equation and Canonical Variational Transition State Theory were utilized to calculate the rate coefficients using the MultiWell suite of programs. CCSD(T), QCISD(T), and M062X with aug-cc-pV(T+d)Z levels of theory were used to calculate the kinetic parameters. Arrhenius equations obtained from fits to the calculated rate coefficients are k1 = (2.66 ± 0.7) × 10-12 exp [(-927 ± 131)/T] cm3 molecule-1 s-1, k2 = (8.99 ± 0.3) × 10-12 exp [(-957 ± 19)/T] cm3 molecule-1 s-1, k3 = (1.51 ± 0.16) × 10-11 exp [(-714 ± 54)/T] cm3 molecule-1 s-1, and k4 = (9.17 ± 1.8) × 10-12 exp [(612 ± 101)/T] cm3 molecule-1 s-1 for the reactions of CHCl3, CH2Cl2, CH2ClCH2Cl, and C2HCl3 with Cl atoms, respectively. The rate coefficients for the reactions of chlorinated VSLSs with Cl atoms from this study are compared with the most recent recommended values from the NASA/JPL and IUPAC evaluations and with literature values. The reactivity trends of the reactions are discussed.

Cl atoms-initiated degradation of 1-Chlorobutane and 2-Chlorobutane: Kinetics, product analysis and atmospheric implications

The relative rate method was employed to investigate the kinetics of the Cl-initiated reactions of 1-chlorobutane (1-CB) and 2-chlorobutane (2-CB) over 263–363 K, and the measured rate coefficients at room temperature are (1.04 ± 0.24) × 10−10 and (5.84 ± 0.27) × 10−11 cm3 molecule−1 s−1, respectively. The Arrhenius equations for the title reactions were derived to be k1-CB + Cl (T = 263–363 K) = (2.77 ± 0.72) × 10−11 exp [(422 ± 79)/T] and k2-CB + Cl (T = 263–363 K) = (1.40 ± 0.32) × 10−11 exp [(415 ± 70)/T] cm3 molecule−1 s−1, respectively. The products were analysed qualitatively using gas chromatography-mass spectrometry (GC-MS), and the reaction mechanism was proposed for the reactions. The rate coefficients for the title reactions were calculated computationally over the temperature range of 200–400 K using canonical variational transition state theory with appropriate tunnelling corrections at CCSD(T)/6–311++G(2d,2p)//BHandHLYP/6–311++G(2d,2p) level of theory to complement our experimentally measured kinetic parameters. The experimental and theoretical data obtained were used to evaluate the impact of the studied molecules in the troposphere.