The Excited-State N-H Tautomerization Rate in Free-Base Corroles.

Corrole is a tetrapyrrolic dye with a structure that resembles porphyrin, apart from a single missing carbon. The absence of this carbon results in the re-arrangement of the double bonds within the macrocycle, and the presence of three pyrrolic protons in the central cavity in its free-base form. These protons lead to the existence of two distinct tautomeric structures that exist in a dynamic equilibrium. Although the ground-state energies of the tautomers are similar, the excited states show a significant difference in energy which unbalances the equilibrium between the tautomers and results in rapid excited-state tautomerization, favouring one tautomeric species over the other. Although the excited-state tautomerization process has been known for a long time, very few studies have been performed on it, leaving many key aspects of the process poorly understood. Herein we show how ultrafast photoluminescence can be used to experimentally determine the rates of excited-state tautomerization and activation energies of three free-base corrole derivatives thus allowing us to completely describe the excited-state dynamics of the unusual excited state of free-base corrole and opening the door to the development of new materials that can exploit its unique characteristics.

High-Pressure Limit and Pressure-Dependent Rate Rules for ß-Scission Reaction Class of Hydroperoxyl Alkyl Hydroperoxyl Radicals (•P(OOH)2) in Normal-Alkyl Cyclohexanes Combustion.

Chemical kinetic studies of the ß-scission reaction class of hydroperoxyl alkyl hydroperoxyl radicals (•P(OOH)2) from normal-alkyl cyclohexanes are carried out systematically through high-level ab initio calculations. Geometry optimizations and frequency calculations for all species involved in the reactions are performed at the B3LYP/CBSB7 level of theory. Electronic single-point energy calculations are calculated at the CBS-QB3 level of theory. Rate constants for the reactions of ß-scission, in the temperature range of 500-1500 K and the pressure range of 0.01-100 atm, are calculated using transition state theory (TST) and Rice-Ramsberger-Kassel-Marcus/Master-Equation (RRKM/ME) theory taking asymmetric Eckart tunneling corrections and the one-dimensional hindered rotor approximation into consideration. The rate rules are obtained by averaging the rate constants of the representative reactions of this class. These rate rules can greatly assist in constructing more accurate low-temperature combustion mechanisms for normal-alkyl cyclohexanes.

Photobleaching Kinetics and Effect of Solvent in the Photophysical Properties of Indocyanine Green for Photodynamic Therapy.

Indocyanine green is an attractive molecule for photodynamic therapy due to its near infrared absorption, resulting in a higher tissue penetration. However, its quantum yields of the triplet and singlet state have been reported to be low and then, reactive oxygen species are unlikely to be formed. Aiming to understand the ICG role in photodynamic response, its photobleaching behavior in solution has been studied under distinct conditions of CW laser irradiation at 780 and 808 nm, oxygen saturations and solvents. Sensitizer bleaching and photoproduct formation were measured by absorption spectroscopy and analyzed using the PDT bleaching macroscopic model to extract physical parameters. ICG photobleaching occurs even at lower oxygen concentrations, indicating that the molecule presents more than one way of degradation. Photoproducts were produced even in solution of less than 4 % oxygen saturation for both solvents and excitation wavelengths. Also, the amplitude of absorption related to J-dimers was increased during irradiation, but only in 50 % PBS solution. The formation of photoproducts was enhanced in the presence of J-type dimers under low oxygen concentration, and the quantum yields of triplet and singlet states were one order of magnitude and two times higher, respectively, when compared to ICG in distilled H2 O.

Reaction of Amino Acids with Ferrate(VI): Impact of the Carboxylic Group on the Primary Amine Oxidation Kinetics and Mechanism.

Ferrate (Fe(VI)) is a novel oxidant that can be used to mitigate disinfection byproduct (DBP) precursors. However, the reaction of Fe(VI) with organic nitrogen, which is a potential precursor of potent nitrogenous DBPs, remains largely unexplored. The present work aimed to identify the kinetics and products for the reaction of Fe(VI) with primary amines, notably amino acids. A new kinetic model involving ionizable intermediates was proposed and can describe the unusual pH effect on the Fe(VI) reactivity toward primary amines and amino acids. The Fe(VI) oxidation of phenylalanine produced a mixture of nitrile, nitrite/nitrate, amide, and ammonia, while nitroalkane was an additional product in the case of glycine. The product distribution for amino acids significantly differed from that of uncarboxylated primary amines that mainly generate nitriles. A general reaction pathway for primary amines and amino acids was proposed and notably involved the formation of imines, the degradation of which was affected by the presence of a carboxylic group. In comparison, ozonation led to higher yields of nitroalkanes that could be readily converted to potent halonitroalkanes during chlor(am)ination. Based on this study, Fe(VI) can effectively mitigate primary amine-based, nitrogenous DBP precursors with little formation of toxic halonitroalkanes.

Thermodynamic Two-Site Surface Reaction Model for Predicting Munition Constituent Reduction Kinetics with Iron (Oxyhydr)oxides.

Iron (oxyhydr)oxides comprise a significant portion of the redox-active fraction of soils and are key reductants for remediation of sites contaminated with munition constituents (MCs). Previous studies of MC reduction kinetics with iron oxides have focused on the concentration of sorbed Fe(II) as a key parameter. To build a reaction kinetic model, it is necessary to predict the concentration of sorbed Fe(II) as a function of system conditions and the redox state. A thermodynamic framework is formulated that includes a generalized double-layer model that utilizes surface acidity and surface complexation reactions to predict sorbed Fe(II) concentrations that are used for fitting MC reduction kinetics. Monodentate- and bidentate Fe(II)-binding sites are used with individual oxide sorption characteristics determined through data fitting. Results with four oxides (goethite, hematite, lepidocrocite, and ferrihydrite) and four nitro compounds (NB, CN-NB, Cl-NB, and NTO) from six separate studies have shown good agreement when comparing observed and predicted surface area-normalized rate constants. While both site types are required to reproduce the experimental redox titration, only the monodentate site concentration controls the MC reaction kinetics. This model represents a significant step toward predicting the timescales of MC degradation in the subsurface.

Kinetics of Na- and K- uranyl arsenate dissolution.

We integrated aqueous chemistry analyses with geochemical modeling to determine the kinetics of the dissolution of Na and K uranyl arsenate solids (UAs(s)) at acidic pH. Improving our understanding of how UAs(s) dissolve is essential to predict transport of U and As, such as in acid mine drainage. At pH 2, Na0.48H0.52(UO2)(AsO4)(H2O)2.5(s) (NaUAs(s)) and K0.9H0.1(UO2)(AsO4)(H2O)2.5(s) (KUAs(s)) both dissolve with a rate constant of 3.2 × 10-7 mol m-2 s-1, which is faster than analogous uranyl phosphate solids. At pH 3, NaUAs(s) (6.3 × 10-8 mol m-2 s-1) and KUAs(s) (2.0 × 10-8 mol m-2 s-1) have smaller rate constants. Steady-state aqueous concentrations of U and As are similarly reached within the first several hours of reaction progress. This study provides dissolution rate constants for UAs(s), which may be integrated into reactive transport models for risk assessment and remediation of U and As contaminated waters.

Shedding light on the calculation of electrode electroactive area and heterogeneous electron transfer rate constants at graphite screen-printed electrodes.

We present in detail the most known and commonly used methods for the calculation of electrode electroactive area ([Formula: see text]) and heterogeneous electron transfer rate constants ([Formula: see text]). The correct procedure for the calculation of these parameters is often disregarded due to either lack of a minimum theoretical background or oversimplification of each method's limitations and prerequisites. The aim of this work is to provide the theoretical background as well as a detailed guide for the implementation of these measurements by impressing upon the electrochemists the parameters that need to be considered so that the obtained results are safe and useful. Using graphite screen-printed electrodes, [Formula: see text], and [Formula: see text] were calculated with different methods and techniques. Data are compared and discussed.

Chemical repair of damaged leucine and tryptophane by thiophenols at close to diffusion-controlled rates: Mechanisms and kinetics.

Thiophenols are chemical species with multiple desirable biological properties, including their primary and secondary antioxidant capacity. In this work, the repairing antioxidant activity of eight different thiophenols has been investigated for damaged leucine and tryptophane. The investigation was carried out employing quantum mechanical and transition state methods to calculate the thermodynamic and kinetic data of the reactions involved, while simulating the biological conditions at physiological pH and aqueous and lipidic medium. The analysis of the atomic charges and the spin densities at each of the points on the potential energy surface was the tool that allowed the elucidation of the reaction mechanisms through which thiophenols repair the oxidative damage caused to the amino acids leucine and tryptophan. It was found that thiophenols can repair leucine via a hydrogen atom transfer mechanism in a manner which is similar to the one used by glutathione to repair the carbon-centered radicals of guanosine. In addition, thiophenols can also restore tryptophane, a nitrogen-centered radical, via proton-coupled electron transfer and single electron transfer mechanisms. Moreover, both processes occur at close to diffusion-controlled rates.

Machine Learning-Assisted QSAR Models on Contaminant Reactivity Toward Four Oxidants: Combining Small Data Sets and Knowledge Transfer.

To develop predictive models for the reactivity of organic contaminants toward four oxidants─SO4•-, HClO, O3, and ClO2─all with small sample sizes, we proposed two approaches: combining small data sets and transferring knowledge between them. We first merged these data sets and developed a unified model using machine learning (ML), which showed better predictive performance than the individual models for HClO (RMSEtest: 2.1 to 2.04), O3 (2.06 to 1.94), ClO2 (1.77 to 1.49), and SO4•- (0.75 to 0.70) because the model "corrected" the wrongly learned effects of several atom groups. We further developed knowledge transfer models for three pairs of the data sets and observed different predictive performances: improved for O3 (RMSEtest: 2.06 to 2.01)/HClO (2.10 to 1.98), mixed for O3 (2.06 to 2.01)/ClO2 (1.77 to 1.95), and unchanged for ClO2 (1.77 to 1.77)/HClO (2.1 to 2.1). The effectiveness of the latter approach depended on whether there was consistent knowledge shared between the data sets and on the performance of the individual models. We also compared our approaches with multitask learning and image-based transfer learning and found that our approaches consistently improved the predictive performance for all data sets while the other two did not. This study demonstrated the effectiveness of combining small, similar data sets and transferring knowledge between them to improve ML model performance.

Application of Nucleotide-Based Kinetic Modeling Approaches to Predict Antibiotic Resistance Gene Degradation during UV- and Chlorine-Based Wastewater Disinfection Processes: From Bench- to Full-Scale.

This study investigated antibiotic resistance gene (ARG) degradation kinetics in wastewaters during bench- and full-scale treatment with UV light and chlorine─with the latter maintained as free available chlorine (FAC) in low-ammonia wastewater and converted into monochloramine (NH2Cl) in high-ammonia wastewater. Twenty-three 142-1509 bp segments (i.e., amplicons) of seven ARGs (blt, mecA, vanA, tet(A), ampC, blaNDM, blaKPC) and the 16S rRNA gene from antibiotic resistant bacteria (ARB) strains Bacillus subtilis, Staphylococcus aureus, Enterococcus faecium, Escherichia coli, Pseudomonas aeruginosa, and Klebsiella pneumoniae were monitored as disinfection targets by qPCR. Rate constants for ARG and 16S rRNA gene amplicon degradation by UV, FAC, and NH2Cl were measured in phosphate buffer and used to expand and validate several recently developed approaches to predict DNA segment degradation rate constants based solely on their nucleotide contents, which were then applied to model ARG degradation during bench-scale treatment in buffer and wastewater matrixes. Kinetics of extracellular and intracellular ARG degradation by UV and FAC were well predicted up to ∼1-2-log10 elimination, although with decreasing accuracy at higher levels for intracellular genes, while NH2Cl yielded minimal degradation under all conditions (agreeing with predictions). ARB inactivation kinetics varied substantially across strains, with intracellular gene degradation lagging cell inactivation in each case. ARG degradation levels observed during full-scale disinfection at two wastewater treatment facilities were consistent with bench-scale measurements and predictions, where UV provided ∼1-log10 ARG degradation, and chlorination of high-ammonia wastewater (dominated by NH2Cl) yielded minimal ARG degradation.

Modeling the Reduction Kinetics of Munition Compounds by Humic Acids.

Dissolved organic matter (DOM) comprises a sizeable portion of the redox-active constituents in the environment and is an important reductant for the abiotic transformation of nitroaromatic compounds and munition constituents (NACs/MCs). Building a predictive kinetic model for these reactions would require the energies associated with both the reduction of the NACs/MCs and the oxidation of the DOM. The heterogeneous and unknown structure of DOM, however, has prohibited reliable determination of its oxidation energies. To overcome this limitation, humic acids (HAs) were used as model DOM, and their redox moieties were modeled as a collection of quinones of different redox potentials. The reduction and oxidation energies of the NACs/MCs and hydroquinones, respectively, via hydrogen atom transfer (HAT) reactions were then calculated quantum chemically. HAT energies have been used successfully in a linear free energy relationship (LFER) to predict second-order rate constants for NAC reduction by hydroquinones. Furthermore, a linear relationship between the HAT energies and the reduction potentials of quinones was established, which allows estimation of hydroquinone reactivity (i.e., rate constants) from HA redox titration data. A training set of three HAs and two NACs/MCs was used to generate a mean HA redox profile that successfully predicted reduction kinetics in multiple HA/MC systems.

New Light on an Old Story: Breaking Kasha's Rule in Phosphorescence Mechanism of Organic Boron Compounds and Molecule Design.

In this work, the phosphorescence mechanism of (E)-3-(((4-nitrophenyl)imino)methyl)-2H-thiochroman-4-olate-BF2 compound (S-BF2) is investigated theoretically. The phosphorescence of S-BF2 has been reassigned to the second triplet state (T2) by the density matrix renormalization group (DMRG) method combined with the multi-configurational pair density functional theory (MCPDFT) to approach the limit of theoretical accuracy. The calculated radiative and non-radiative rate constants support the breakdown of Kasha's rule further. Our conclusion contradicts previous reports that phosphorescence comes from the first triplet state (T1). Based on the revised phosphorescence mechanism, we have purposefully designed some novel compounds in theory to enhance the phosphorescence efficiency from T2 by replacing substitute groups in S-BF2. Overall, both S-BF2 and newly designed high-efficiency molecules exhibit anti-Kasha T2 phosphorescence instead of the conventional T1 emission. This work provides a useful guidance for future design of high-efficiency green-emitting phosphors.

The Homogeneous Gas-Phase Formation Mechanism of PCNs from Cross-Condensation of Phenoxy Radical with 2-CPR and 3-CPR: A Theoretical Mechanistic and Kinetic Study.

Chlorophenols (CPs) and phenol are abundant in thermal and combustion procedures, such as stack gas production, industrial incinerators, metal reclamation, etc., which are key precursors for the formation of polychlorinated naphthalenes (PCNs). CPs and phenol can react with H or OH radicals to form chlorophenoxy radicals (CPRs) and phenoxy radical (PhR). The self-condensation of CPRs or cross-condensation of PhR with CPRs is the initial and most important step for PCN formation. In this work, detailed thermodynamic and kinetic calculations were carried out to investigate the PCN formation mechanisms from PhR with 2-CPR/3-CPR. Several energetically advantageous formation pathways were obtained. The rate constants of key elementary steps were calculated over 600~1200 K using the canonical variational transition-state theory (CVT) with the small curvature tunneling (SCT) contribution method. The mechanisms were compared with the experimental observations and our previous works on the PCN formation from the self-condensation of 2-CPRs/3-CPRs. This study shows that naphthalene and 1-monochlorinated naphthalene (1-MCN) are the main PCN products from the cross-condensation of PhR with 2-CPR, and naphthalene and 2-monochlorinated naphthalene (2-MCN) are the main PCN products from the cross-condensation of PhR with 3-CPR. Pathways terminated with Cl elimination are preferred over those terminated with H elimination. PCN formation from the cross-condensation of PhR with 3-CPR can occur much easier than that from the cross-condensation of PhR with 2-CPR. This study, along with the study of PCN formation from the self-condensation 2-CPRs/3-CPRs, can provide reasonable explanations for the experimental observations that the formation potential of naphthalene is larger than that of 1-MCN using 2-CP as a precursor, and an almost equal yield of 1-MCN and 2-MCN can be produced with 3-CP as a precursor.

Multiplex Label-Free Kinetic Characterization of Antibodies for Rapid Sensitive Cardiac Troponin I Detection Based on Functionalized Magnetic Nanotags.

Express and highly sensitive immunoassays for the quantitative registration of cardiac troponin I (cTnI) are in high demand for early point-of-care differential diagnosis of acute myocardial infarction. The selection of antibodies that feature rapid and tight binding with antigens is crucial for immunoassay rate and sensitivity. A method is presented for the selection of the most promising clones for advanced immunoassays via simultaneous characterization of interaction kinetics of different monoclonal antibodies (mAb) using a direct label-free method of multiplex spectral correlation interferometry. mAb-cTnI interactions were real-time registered on an epoxy-modified microarray glass sensor chip that did not require activation. The covalent immobilization of mAb microdots on its surface provided versatility, convenience, and virtually unlimited multiplexing potential. The kinetics of tracer antibody interaction with the "cTnI­capture antibody" complex was characterized. Algorithms are shown for excluding mutual competition of the tracer/capture antibodies and selecting the optimal pairs for different assay formats. Using the selected mAbs, a lateral flow assay was developed for rapid quantitative cTnI determination based on electronic detection of functionalized magnetic nanoparticles applied as labels (detection limit­0.08 ng/mL, dynamic range > 3 orders). The method can be extended to other molecular biomarkers for high-throughput screening of mAbs and rational development of immunoassays.

Quantitative structure-activity relationship for the oxidation of organic contaminants by peracetic acid using GA-MLR method.

Peracetic acid (PAA) is considered as an effective and powerful oxidant for eliminating organic contaminants in wastewater treatment. The second-order rate constant (kapp) for the reaction of PAA with organic contaminants is practically important for evaluating their removal efficiency in wastewater treatment, but only limited numbers of kapp values are available. In this study, 70 organic compounds with various structures were selected, and the kapp of PAA with each organic compound was used to develop two quantitative structure-activity relationship (QSAR) models based on three kinds of descriptors including constitutional, quantum chemical, and the PaDEL descriptors. The genetic algorithm (GA) was applied to select the molecular descriptors, then the models developed by multiple linear regression (MLR). The most important descriptors that explain the reactivity of organic compounds with PAA are the EHOMO for the model with the constitutional and quantum chemical descriptors. The maxHdsCH and minHdCH2 are two most important descriptors for the model with only PaDEL descriptors. The developed models can be used to predict kapp for a wide range of organic contaminants. The accuracy of the developed models was proved by the internal, external validation and the Y-scrambling technique. The developed QSAR models using the GA-MLR method can be used as a screening tool for predicting the elimination of organic contaminants by PAA and increasing the understanding of chemical pollutant fate.

TPPS4-Sensitized Photooxidation of Micropollutants-Singlet Molecular Oxygen Kinetic Study.

Visible light-sensitized oxidation of micropollutants (MPs) in the presence of meso-tetrakis(4-sulfonatophenyl)porphyrin photosensitizers was studied. In order to explore the role of type I (ROS generation) or type II (singlet oxygen) photooxidation, radical scavengers were used to obtain insight into the mechanism of photodegradation. It was revealed that singlet oxygen is the main ROS taking part in TPPS4- sensitized photooxidation of micropollutants. The interaction of MPs with 1O2 in deuterium oxide (D2O) was investigated by measuring the phosphorescence lifetime of 1O2. The rate constant (kq) for the total (physical and chemical) quenching of 1O2 by MPs was determined in a D2O buffer (pD 7, 9 and 10.8). The rate constants of singlet oxygen quenching and reaction with MPs were determined, and the rate constant of excited TPPS4 quenching by MPs was also estimated.

Reactivity of Chlorine Radicals (Cl• and Cl2•-) with Dissolved Organic Matter and the Formation of Chlorinated Byproducts.

Chlorine radicals, including Cl• and Cl2•-, can be produced in sunlight waters (rivers, oceans, and lakes) or water treatment processes (e.g., electrochemical and advanced oxidation processes). Dissolved organic matter (DOM) is a major reactant with, or a scavenger of, Cl• and Cl2•- in water, but limited quantitative information exists regarding the influence of DOM structure on its reactivity with Cl• and Cl2•-. This study aimed at quantifying the reaction rates and the formation of chlorinated organic byproducts produced from Cl• and Cl2•- reactions with DOM. Laser flash photolysis experiments were conducted to quantify the second-order reaction rate constants of 19 DOM isolates with Cl• (kDOM-Cl•) and Cl2•- (kDOM-Cl2•-), and compare those with the hydroxyl radical rate constants (kDOM-•OH). The values for kDOM-Cl• ((3.71 ± 0.34) × 108 to (1.52 ± 1.56) × 109 MC-1 s-1) were orders of magnitude greater than the kDOM-Cl2•- values ((4.60 ± 0.90) × 106 to (3.57 ± 0.53) × 107 MC-1 s-1). kDOM-Cl• negatively correlated with the weight-averaged molecular weight (MW) due to the diffusion-controlled reactions. DOM with high aromaticity and total antioxidant capacity tended to react faster with Cl2•-. During the same experiments, we also monitored the formation of chlorinated byproducts through the evolution of total organic chlorine (TOCl) as a function of chlorine radical oxidant exposure (CT value). Maximum TOCl occurred at a CT of 4-8 × 10-12 M·s for Cl• and 1.1-2.2 × 10-10 M·s for Cl2•-. These results signify the importance of DOM in scavenging chlorine radicals and the potential risks of producing chlorinated byproducts of unknown toxicity.

Mechanistic Insight into the Reactivities of Aqueous-Phase Singlet Oxygen with Organic Compounds.

Singlet oxygen (1O2) is a selective reactive oxygen species that plays a key role for the fate of various organic compounds in the aquatic environment under sunlight irradiation, engineered water oxidation systems, atmospheric water droplets, and biomedical systems. While the initial rate-determining charge-transfer reaction mechanisms and kinetics of 1O2 have been studied extensively, no comprehensive studies have been performed to elucidate the reaction mechanisms with organic compounds that have various functional groups. In this study, we use density functional theory calculations to determine elementary reaction mechanisms with a wide variety of organic compounds. The theoretically calculated aqueous-phase free energies of activation of single electron transfer and 1O2 addition reactions are compared to the experimentally determined rate constants in the literature to determine linear free-energy relationships. The theoretically calculated free energies of activation for the groups of phenolates and phenols show excellent correlations with the Hammett constants that accept electron densities by through-resonance. The dominant elementary reaction mechanism is discussed for each group of compounds. As a practical implication, we demonstrate the fate of environmentally relevant organic compounds induced by photochemically produced intermediate species at different pH and evaluate the impact of predicting rate constants to the half-life.

Do epoch lengths of hypnotic depth indicators affect estimated of blood-brain equilibration rate constants of propofol?

This study aimed to investigate the effect of epoch length of hypnotic depth indicators on the blood-brain equilibration rate constant (ke0) estimates of propofol. Propofol was administered by zero-order infusion (1.5, 3.0, 6, and 12 mg·kg-1·h-1) for one hour in 63 healthy volunteers. The ke0 of propofol was estimated using an effect-compartment model linking pharmacokinetics and pharmacodynamics, in which response variables were electroencephalographic approximate entropy (ApEn) or bispectral index (BIS) (n = 32 each for propofol infusion rates of 6 and 12 mg·kg-1·h-1). Epoch lengths of ApEn were 2, 10, 30, and 60 seconds (s). The correlations between plasma propofol concentrations (Cp) and BIS and ApEn 2, 10, 30, and 60 s were determined, as was the Ce associated with 50% probability of unconsciousness (Ce50,LOC). The pharmacokinetics of propofol were well described by a three-compartment model. The correlation coefficient between Cp and ApEn 2, 10, 30, and 60 s were -0.64, -0.54, -0.39, and -0.26, respectively, whereas correlation coefficient between Cp and BIS was -0.74. The blood-brain equilibration half-life based on the ke0 estimates for ApEn at 2, 10, 30, 60 s and BIS were 4.31, 3.96, 5.78. 6.54, 5.09 min, respectively, whereas the Ce50,LOC for ApEn at 2, 10, 30, 60 s and BIS were 1.55, 1.47, 1.28, 1.04, and 1.55 µg·ml-1, respectively. Since ke0, which determines the onset of drug action, varies according to the epoch length, it is necessary to consider the epoch length together when estimating ke0.

Variational transition state theory rate constants and H/D kinetic isotope effects for CH3 + CH3 OCOH reactions.

The rate constants and H/D kinetic isotope effect for hydrogen abstraction reactions involving isotopomers of methyl formate by methyl radical are computed employing methods of the variational transition state theory (VTST) with multidimensional tunneling corrections. The energy paths were built with a dual-level method using the moller plesset second-order perturbation theory (MP2) method as the low-level and complete basis set (CBS) extrapolation as the high-level energy method. Benchmark calculations with the CBSD-T approach give an enthalpy of reaction at 0 K for R1 (-4.5 kcal/mol) and R2 (-4.2 kcal/mol) which are in good agreement with the experiment, that is, -4.0 and - 4.8 kcal/mol. For the reactional paths involving the isotopomers CH3 + CH3 OCOH → CH4 + CH3 OCO and CH3 + CH3 OCOD → CH3 D + CH3 OCO, the value of kH /kD (T = 455 K) using the canonical VTST/small-curvature tunneling approximation method is 6.7 in close agreement with experimental value (6.2). © 2019 Wiley Periodicals, Inc.