The13C and D kinetic isotope effects in the reaction of CH4 with Cl

State-to-state integral cross sections for the Cl + CHD3 (vCH = 0,1) → HCl + CD3 reactions

To precisely study the effect of CH stretch excitation, we calculate the first seven-dimensional state-to-state integral cross sections (ICSs) for the Cl + CHD3(vCH = 0,1) → HCl + CD3 reactions using the time-dependent wave-packet approach on a highly accurate fundamental invariant neural network potential energy surface. Compared to experiments that focus solely on CD3(v2 = 0), our theoretical calculations give a more global picture. It is found that the vibrational enhancement factors of the total ICS are less than 1 only in the very low collision energy region (<0.55 kcal/mol); moreover, the CH stretch-excited reaction also produces vibrationally cold HCl. These two findings indicate that Polanyi's rules apply to this late-barrier polyatomic reaction, not only in terms of the energy requirement in the breaking bond but also in the energy disposal in the newly formed bond. In addition, although the experimental and theoretical product HCl(v1 = 1) branching fractions for CD3(v2 = 0) are in quantitative agreement, the new quantum vibrational enhancement factor for CD3(v2 = 0) is ∼2.7 times greater than the experimental results.

Accurate Quantum Dynamics Calculations for the Cl + CH4/CHD3/CD4 Reaction Rates

Full-dimensional quantum dynamics simulations of the reaction of Cl with methane and its isotopomers are reported. Thermal rate constants are computed for the Cl + CH4 → HCl + CH3, Cl + CHD3 → HCl + CD3, and Cl + CD4 → DCl + CD3 reactions. Temperatures between 200 and 500 K are considered. In this temperature range, excellent agreement with the experiment is obtained. A detailed analysis of the kinetic isotope effect reveals the crucial importance of the CH3/CD3 umbrella motion. Comparison with approximate ring-polymer molecular dynamics simulations shows significant differences depending on the isotope studied.

Isotope Effect and Heavy-Light-Heavy Reactivity Oscillation in the Cl + CHD3/CHT3 Reaction

Previous experiments and theories have shown the existence of heavy-light-heavy (HLH) reactivity oscillation in the Cl + CH4 reaction and anticipated that similar oscillations should exist in many HLH reactions involving polyatomic reagents. However, the total reaction probabilities for the Cl + CHD3 → HCl + CD3 reaction exhibit only a step-like feature, and the total reaction probabilities for Cl + CHT3 → HCl + CT3 do not show any structure at all. Here, we report seven-dimensional state-to-state quantum dynamics studies for this reaction on the FI-NN PES, and we demonstrate that HLH reactivity oscillations also exist in these two reactions, manifesting as peaks in the reaction probabilities for low product rotational states. These oscillations, however, are obscured in the total reaction probability because of the higher excitation of j ≥ 2 product rotational states. Furthermore, the isotope replacement of nonreactive hydrogen with deuterium and tritium significantly enhances reactivity at collision energies above 0.112 eV, indicating an inverse secondary isotope effect on the probabilities, which is proved to be also caused by HLH mass combination. We also demonstrate that the highly rotational excitation of CHD3 substantially enhances reactivity and the HLH oscillations, similar to HLH triatomic reactions. These observations are completely different from those in the H + CHD3 reaction, which is also a late-barrier reaction. Therefore, the HLH mass combination is very important, which affects not only the reactivity oscillation but also the amplitude and product rotational state distribution and makes the initial rotation excitation play a pivotal role in the reaction.

Hydrogen/Deuterium Transfer from Anisole to Methoxy Radicals: A Theoretical Study of a Deuterium-Labeled Drug Model

Recently, deuterium-labeled drugs, such as deutetrabenazine, have attracted considerable attention. Consequently, understanding the reaction mechanisms of deuterium-labeled drugs is crucial, both fundamentally and for real applications. To understand the mechanisms of H- and D-transfer reactions, in this study, we used deuterated anisole as a deutetrabenazine model and computationally considered the nuclear quantum effects of protons, deuterons, and electrons. We demonstrated that geometrical differences exist in the partially and fully deuterated methoxy groups and hydrogen-bonded structures of intermediates and transition states due to the H/D isotope effect. The observed geometrical features and electronic structures are ascribable to the different nuclear quantum effects of protons and deuterons. Primary and secondary kinetic isotope effects (KIEs) were calculated for H- and D-transfer reactions from deuterated and undeuterated anisole, with the calculated primary KIEs in good agreement with the corresponding experimental data. These results reveal that the nuclear quantum effects of protons and deuterons need to be considered when analyzing the reaction mechanisms of H- and D-transfer reactions and that a theoretical approach that directly includes nuclear quantum effects is a powerful tool for the analysis of H/D isotope effects in H- and D-transfer reactions.

Improved Constraints on Global Methane Emissions and Sinks Using δ 13C-CH4

We study the drivers behind the global atmospheric methane (CH4) increase observed after 2006. Candidate emission and sink scenarios are constructed based on proposed hypotheses in the literature. These scenarios are simulated in the TM5 tracer transport model for 1984-2016 to produce three-dimensional fields of CH4 and δ 13C-CH4, which are compared with observations to test the competing hypotheses in the literature in one common model framework. We find that the fossil fuel (FF) CH4 emission trend from the Emissions Database for Global Atmospheric Research 4.3.2 inventory does not agree with observed δ 13C-CH4. Increased FF CH4 emissions are unlikely to be the dominant driver for the post-2006 global CH4 increase despite the possibility for a small FF emission increase. We also find that a significant decrease in the abundance of hydroxyl radicals (OH) cannot explain the post-2006 global CH4 increase since it does not track the observed decrease in global mean δ 13C-CH4. Different CH4 sinks have different fractionation factors for δ 13C-CH4, thus we can investigate the uncertainty introduced by the reaction of CH4 with tropospheric chlorine (Cl), a CH4 sink whose abundance, spatial distribution, and temporal changes remain uncertain. Our results show that including or excluding tropospheric Cl as a 13 Tg/year CH4 sink in our model changes the magnitude of estimated fossil emissions by ∼20%. We also found that by using different wetland emissions based on a static versus a dynamic wetland area map, the partitioning between FF and microbial sources differs by 20 Tg/year, ∼12% of estimated fossil emissions.

An accurate potential energy surface and ring polymer molecular dynamics study of the Cl + CH4→ HCl + CH3reaction

Thermal rate coefficients for the Cl + CH₄/CD₄reactions were studied on a new full-dimensional accurate potential energy surface with the spin–orbit corrections considered in the entrance channel.

Full-dimensional analytical potential energy surface describing the gas-phase Cl + C2H6 reaction and kinetics study of rate constants and kinetic isotope effects

Within the Born–Oppenheimer approximation a full-dimensional analytical potential energy surface, PES-2017, was developed for the gas-phase hydrogen abstraction reaction between the chlorine atom and ethane, which is a nine body system.

Light-Driven C-H Oxygenation of Methane into Methanol and Formic Acid by Molecular Oxygen Using a Perfluorinated Solvent

The chlorine dioxide radical (ClO₂^. ) was found to act as an efficient oxidizing agent in the aerobic oxygenation of methane to methanol and formic acid under photoirradiation. Photochemical oxygenation of methane occurred in a two-phase system comprising perfluorohexane and water under ambient conditions (298 K, 1 atm). The yields of methanol and formic acid were 14 and 85 %, respectively, with a methane conversion of 99 % without formation of the further oxygenated products such as CO₂ and CO. Ethane was also photochemically converted into ethanol (19 %) and acetic acid (80 %). The methane oxygenation is initiated by the photochemical Cl-O bond cleavage of ClO₂^. to generate Cl^. and O₂ . The produced Cl^. reacts with CH₄ to form a methyl radical (CH₃^. ). Finally, the oxygenated products such as methanol and formic acid were given by the radical chain reaction. A fluorous solvent plays an important role of inhibiting the deactivation of reactive radical species such as Cl^. and CH₃^. .

Glacial/interglacial wetland, biomass burning, and geologic methane emissions constrained by dual stable isotopic CH4 ice core records

Atmospheric methane (CH4) records reconstructed from polar ice cores represent an integrated view on processes predominantly taking place in the terrestrial biogeosphere. Here, we present dual stable isotopic methane records [δ13CH4 and δD(CH4)] from four Antarctic ice cores, which provide improved constraints on past changes in natural methane sources. Our isotope data show that tropical wetlands and seasonally inundated floodplains are most likely the controlling sources of atmospheric methane variations for the current and two older interglacials and their preceding glacial maxima. The changes in these sources are steered by variations in temperature, precipitation, and the water table as modulated by insolation, (local) sea level, and monsoon intensity. Based on our δD(CH4) constraint, it seems that geologic emissions of methane may play a steady but only minor role in atmospheric CH4 changes and that the glacial budget is not dominated by these sources. Superimposed on the glacial/interglacial variations is a marked difference in both isotope records, with systematically higher values during the last 25,000 y compared with older time periods. This shift cannot be explained by climatic changes. Rather, our isotopic methane budget points to a marked increase in fire activity, possibly caused by biome changes and accumulation of fuel related to the late Pleistocene megafauna extinction, which took place in the course of the last glacial.

Temperature dependence of the Cl atom reaction with deuterated methanes

Kinetic isotope effect (KIE) and reaction rate coefficients, k1-k4, for the gas-phase reaction of Cl atoms with (12)CH3D (k1), (12)CH2D2 (k2), (12)CHD3 (k3), and (12)CD4 (k4) over the temperature range 223-343 K in 630 Torr of synthetic air are reported. Rate coefficients were measured using a relative rate technique with (12)CH4 as the primary reference compound. Fourier transform infrared spectroscopy was used to monitor the methane isotopologue loss. The obtained KIE values were (12)CH3D: KIE1(T) = (1.227 ± 0.004) exp((43 ± 5)/T); (12)CH2D2: KIE2(T) = (1.14 ± 0.20) exp((191 ± 60)/T); (12)CHD3: KIE3(T) = (1.73 ± 0.34) exp((229 ± 60)/T); and (12)CD4: KIE4(T) = (1.01 ± 0.3) exp((724 ± 19)/T), where KIEx(T) = kCl+(12)CH4(T)/kx(T). The quoted uncertainties are at the 2σ (95% confidence) level and represent the precision of our data. The following Arrhenius expressions and 295 K rate coefficient values (in units of cm(3) molecule(-1) s(-1)) were derived from the above KIE using a rate coefficient of 7.3 × 10(-12) exp(-1280/T) cm(3) molecule(-1) s(-1) for the reaction of Cl with (12)CH4: k1(T) = (5.95 ± 0.70) × 10(-12) exp(-(1323 ± 50)/T), k1(295 K) = (6.7 ± 0.8) × 10(-14); k2(T) = (6.4 ± 1.3) × 10(-12) exp(-(1471 ± 60)/T), k2(295 K) = (4.4 ± 0.9) × 10(-14); k3(T) = (4.2 ± 1.0) × 10(-12) exp(-(1509 ± 60)/T), k3(295 K) = (2.53 ± 0.6) × 10(-14); and k4(T) = (7.13 ± 2.3) × 10(-12) exp(-(2000 ± 120)/T), k4(295 K) = (0.81 ± 0.26) × 10(-14). The reported uncertainties in the pre-exponential factors are 2σ and include estimated systematic errors in our measurements and the uncertainty in the reference reaction rate coefficient. The results from this study are compared with previously reported room-temperature rate coefficients for each of the deuterated methanes as well as the available temperature dependent data for the Cl atom reactions with CH3D and CD4. A two-dimensional atmospheric chemistry model was used to examine the implications of the present results to the atmospheric lifetime and vertical variation in the loss of the deuterated methane isotopologues. The relative contributions of the reactions of OH, Cl, and O((1)D) to the loss of the isotopologues in the stratosphere were also examined. The results of the calculations are described and discussed.

Quantum rate coefficients and kinetic isotope effect for the reaction Cl + CH4 → HCl + CH3 from ring polymer molecular dynamics

Thermal rate coefficients and kinetic isotope effect have been calculated for prototypical heavy-light-heavy polyatomic bimolecular reactions Cl + CH4/CD4 → HCl/DCl + CH3/CD3, using a recently proposed quantum dynamics approach: ring polymer molecular dynamics (RPMD). Agreement with experimental rate coefficients, which are quite scattered, is satisfactory. However, differences up to 50% have been found between the RPMD results and those obtained from the harmonic variational transition-state theory on one of the two full-dimensional potential energy surfaces used in the calculations. Possible reasons for such discrepancy are discussed. The present work is an important step in a series of benchmark studies aimed at assessing accuracy for RPMD for chemical reaction rates, which demonstrates that this novel method is a quite reliable alternative to previously developed techniques based on transition-state theory.

Potential energy surface, kinetics, and dynamics study of the Cl+CH4→HCl+CH3 reaction

A modified and recalibrated potential energy surface for the gas-phase Cl+CH4-->HCl+CH3 reaction is reported and tested. It is completely symmetric with respect to the permutation of the four methane hydrogen atoms and is calibrated with respect to updated experimental and theoretical stationary point properties and experimental forward thermal rate constants. From the kinetics point of view, the forward and reverse thermal rate constants and the activation energies were calculated using the variational transition-state theory with semiclassical transmission coefficients over a wide temperature range of 150-2500 K. The theoretical results reproduce the available experimental data, with a small curvature of the Arrhenius plot which indicates the role of tunneling in this hydrogen abstraction reaction. A dynamics study was also performed on this PES using quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories. First, we found a noticeable internal energy in the coproduct methyl radical, both in the ground-state [CH4 (v=0)] and vibrationally excited [CH4 (v=1)] reactions. This CH3 internal energy was directly precluded in some experiments or oversimplified in previous theoretical studies using pseudotriatomic models. Second, our QCT calculations give HCl rotational distributions slightly hotter than those in experiment, but correctly describing the experimental trend of decreasing the HCl product rotation excitation in going from HCl (v'=0) to HCl (v'=1) for the CH4 (v=1) reaction. Third, the state specific scattering distributions present a reasonable agreement with experiment, although they tend to make the reaction more forward and backward scattered than found experimentally probably because of the hotter rotational distribution and the deficiencies of the QCT methods.

Study of the Carbon-13 and Deuterium Kinetic Isotope Effects in the Cl and OH Reactions of CH4 and CH3Cl

Relative rate experiments have been carried out for three isotopologues of chloromethane and their reactions with Cl atoms and OH radicals. The OH and Cl reaction rates of CH2DCl and CHD2Cl were measured by long-path FTIR spectroscopy relative to CH3Cl at 298+/-2 K and 1013+/-10 hPa in purified air. The FTIR spectra were fitted using a nonlinear least squares spectral fitting method including measured high-resolution infrared spectra as references. The relative reaction rates defined by alpha=klight/kheavy were determined to be kOH+CH3Cl/kOH+CH2DCl=1.41+/-0.05, kOH+CH3Cl/kOH+CHD2Cl=2.03+/-0.05, kCl+CH3Cl/kCl+CH2DCl=1.42+/-0.04, and kCl+CH3Cl/kCl+CHD2Cl=2.27+/-0.04. The carbon-13 and deuterium kinetic isotope effects in the OH and Cl reactions of CH3Cl were investigated further using variational transition state theory, and the results were compared to similar calculations performed for the CH4+OH/Cl reaction systems. The calculations show that the order of magnitude difference for the carbon-13 kinetic isotope effect in the OH reaction of CH3Cl compared to CH4 reported by Gola et al. (Atmos. Chem. Phys. 2005, 5, 2395) can be explained by the lower barrier to internal rotation of the OH radical in the transition state of the CH4+OH reaction than in the CH3Cl+OH reaction. The deuterium kinetic isotope effects can be explained in terms of combined variational effects and tunneling.