Temperature and Isotope Dependent Kinetics of Nickel-Catalyzed Oxidation of Methane by Ozone

Designing main-group catalysts for low-temperature methane combustion by ozone

The catalytic combustion of methane at a low temperature is becoming increasingly key to controlling unburned CH₄ emissions from natural gas vehicles and power plants, although the low activity of benchmark platinum-group-metal catalysts hinders its broad application. Based on automated reaction route mapping, we explore main-group elements catalysts containing Si and Al for low-temperature CH₄ combustion with ozone. Computational screening of the active site predicts that strong Brønsted acid sites are promising for methane combustion. We experimentally demonstrate that catalysts containing strong Bronsted acid sites exhibit improved CH₄ conversion at 250 °C, correlating with the theoretical predictions. The main-group catalyst (proton-type beta zeolite) delivered a reaction rate that is 442 times higher than that of a benchmark catalyst (5 wt% Pd-loaded Al₂O₃) at 190 °C and exhibits higher tolerance to steam and SO₂. Our strategy demonstrates the rational design of earth-abundant catalysts based on automated reaction route mapping.

Kinetics of O3 with Ca+ and Its Higher Oxides CaOn+ (n = 1-3) and Updates to a Model of Meteoric Calcium in the Mesosphere and Lower Thermosphere

The room-temperature rate constants and product branching fractions of CaO_n⁺ (n = 0-3) + O₃ are measured using a selected ion flow tube apparatus. Ca⁺ + O₃ produces CaO⁺ + O₂ with k = 9 ± 4 × 10^-10 cm³ s^-1, within uncertainty equal to the Langevin capture rate constant. This value is significantly larger than several literature values. Most likely, those values were underestimated due to the reformation of Ca⁺ from the sequential chemistry of higher calcium oxide cations with O₃, as explored here. A rate constant of 8 ± 3 × 10^-10 cm³ s^-1 is recommended. Both CaO⁺ and CaO₂⁺ react near the capture rate constant with ozone. The CaO⁺ reaction yields both CaO₂⁺ + O₂ (0.80 ± 0.15 branching) and Ca⁺ + 2O₂. Similarly, the CaO₂⁺ reaction yields both CaO₃⁺ + O₂ (0.85 ± 0.15 branching) and CaO⁺ + 2O₂. CaO₃⁺ + O₃ yield CaO₂⁺ + 2O₂ at 2 ± 1 × 10^-11 cm³ s^-1, about 2% of the capture rate constant. The results are supported using density functional calculations and statistical modeling. In general, CaO_n⁺ + O₃ yield CaO_n+1⁺ + O₂, the expected oxidation. Some fraction of CaO_n+1⁺ is produced with sufficient internal energy to further dissociate to CaO_n-1⁺ + O₂, yielding the same products as the oxidation of O₃ by CaO_n⁺. Mesospheric Ca and Ca⁺ concentrations are modeled as functions of day, latitude, and altitude using the Whole Atmosphere Community Climate Model (WACCM); incorporating the updated rate constants improves agreement with concentrations derived from lidar measurements.

Gas-Phase Reactivity of Ozone with Lanthanide Ions (Sm+, Nd+) and Their Higher Oxides

The kinetics of SmO_n⁺ (n = 0-2) and NdO_n⁺ (n = 0-2) with O₃ are measured using a selected-ion flow tube. Reaction of Nd⁺ to yield NdO⁺ + O₂ occurs rapidly, with a rate constant near the capture-controlled limit of ∼8 × 10^-10 cm³ s^-1. NdO⁺ reacts at ∼40% of the capture limit to yield NdO₂⁺ with little temperature dependence from 200 to 400 K. NdO₂⁺ likely reacts very slowly (k ∼ 10^-13 cm³ s^-1) to yield NdO⁺ + 2O₂, does not react to yield NdO₃⁺, and associates slowly (k ∼ 10^-12 cm³ s^-1) to yield NdO₂⁺(O₃)_1-3. Reaction of Sm⁺ also yields SmO⁺ at near the capture limit at all temperatures, but a significant fraction (∼50%) of the SmO⁺ is produced in excited states that are long-lived compared to the millisecond time scale of the experiment. These states are evidently resistant to both radiative and collisional relaxation. The excited-state production is likely due to a spin-conservation constraint on the reaction, despite the large spin-orbit coupling typical for lanthanide-containing species. Ground-state SmO⁺ reacts inefficiently (k = 2 × 10^-11 (T/300)^-2.5 cm³ s^-1) to yield SmO₂⁺ + O₂, while the excited-state SmO^+* reacts at the capture limit, with branching to yield Sm⁺ + 2O₂ (ΔH_r,0K = 148.7 ± 0.4 kJ mol^-1 for ground-state SmO⁺) approximately 60% of the time, the remainder forming SmO₂⁺, which further reacts with O₃ to yield SmO⁺ at about 1% of the collisional value.

Cyclotrimerization of Acetylene under Thermal Conditions: Gas-Phase Kinetics of V+ and Fe+ + C2H2

The kinetics of successive reactions of acetylene (C₂H₂) initiated on either vanadium or iron atomic cations have been investigated under thermal conditions using the variable-ion source and temperature-adjustable selected-ion flow tube apparatus. Consistent with the literature results, the reaction of Fe⁺ + C₂H₂ primarily yields Fe⁺(m/z = (C₂H₂)₃); however, analysis via quantum chemical calculations and statistical modeling shows that the mechanism does not form benzene upon the third acetylene addition. The kinetics are more consistent with successive addition of three acetylene molecules, yielding Fe⁺(C₂H₂)₃, followed by an addition of a fourth acetylene molecule, initiating cyclotrimerization, yielding either Fe⁺(C₂H₂) + neutral benzene or Fe⁺(Bz) + acetylene, where Bz is a benzene ligand. In contrast, the reaction of V⁺ + C₂H₂ yields products via successive associations V⁺(m/z = (C₂H₂)_n) either with or without a bimolecular step involving loss of one H₂ and V⁺C₂(m/z = (C₂H₂)_m), where n and m extend at least up to 11 under conditions of 0.32 Torr at 300 K. Stabilized V⁺(Bz) is not a significant intermediate in the association mechanism. We propose a plausible mechanism for the generation of neutral benzene in this reaction and compare with the Fe⁺ results. The reaction steps that produce benzene result in turnover of the single-atom catalyst, and the large hydrocarbons produced that remain associated to the catalyst are proposed to be polycyclic aromatic hydrocarbons.

Bonding, Thermodynamics, and Dissociation Dynamics of NiO+ and NiS+ Determined by Photofragment Imaging and Theory

We use photofragment ion imaging and ab initio calculations to determine the bond strength and photodissociation dynamics of the nickel oxide (NiO⁺) and nickel sulfide (NiS⁺) cations. NiO⁺ photodissociates broadly from 20350 to 32000 cm^-1, forming ground state products Ni⁺(²D) + O(³P) below ∼29000 cm^-1. Above this energy, Ni⁺(⁴F) + O(³P) products become accessible and dominate over the ground state channel. In certain images, product spin-orbit levels are resolved, and spin-orbit propensities are determined. Image anisotropy and the results of MRCI calculations suggest NiO⁺ photodissociates via a 3 ⁴Σ^- ← X ⁴Σ^- transition above the Ni⁺(⁴F) threshold and via 3 ⁴Σ^-, 2 ⁴Σ^-, and/or 2 ⁴Π and 3 ⁴Π excited states below the ⁴F threshold. The photodissociation spectrum of NiS⁺ from 19900 to 23200 cm^-1 is highly structured, with ∼12 distinct vibronic peaks, each containing underlying substructure. Above 21600 cm^-1, the Ni⁺(²D_5/2) + S(³P) and Ni⁺(²D_3/2) + S(³P) product spin-orbit channels compete, with a branching ratio of ∼2:1. At lower energy, Ni⁺(²D_5/2) is formed exclusively, and S(³P₂) and S(³P₁) spin-orbit channels are resolved. MRCI calculations predict the ground state of NiS⁺ to be one of two nearly degenerate states, the 1 ⁴Σ^- and 1 ⁴Δ. Based on images and spectra, the ground state of NiS⁺ is assigned as ⁴Δ_7/2, with the 1 ⁴Σ_3/2^- and 1 ⁴Σ_1/2^- states 81 ± 30 and 166 ± 50 cm^-1 higher in energy, respectively. The majority of the photodissociation spectrum is assigned to transitions from the 1 ⁴Δ state to two overlapping, predissociative excited ⁴Δ states. Our D₀ measurements for NiO⁺ (D₀ = 244.6 ± 2.4 kJ/mol) and NiS⁺ (D₀ = 240.3 ± 1.4 kJ/mol) are more precise and closer to each other than previously reported values. Finally, using a recent measurement of D₀(NiS), we derive a more precise value for IE (NiS): 8.80 ± 0.02 eV (849 ± 1.7 kJ/mol).

Large Anharmonic Effects on Tunneling and Kinetics: Reaction of Propane with Muonium

Calculations of kinetic isotope effects (KIEs) provide challenging tests of quantal mass effects on reaction rates, and muonium KIEs are the most challenging. Here, we show that it can be very important to include reaction-coordinate-dependent vibrational anharmonicity along the whole reaction path to calculate tunneling probabilities and KIEs. For the reaction of propane with Mu, this decreases both the height and width of the vibrationally adiabatic potential barrier, with both effects increasing the rate constants. Our results agree well with the experimental observations.

A study of the reactions of Ni+ and NiO+ ions relevant to planetary upper atmospheres

The reactions between Ni+(2D) and O3, O2, N2, CO2 and H2O were studied at 294 K using the pulsed laser ablation at 532 nm of a nickel metal target in a fast flow tube, with mass spectrometric detection of Ni+ and NiO+. The rate coefficient for the reaction of Ni+ with O3 is k(294 K) = (9.7 ± 2.1) × 10-10 cm3 molecule-1 s-1; the reaction proceeds at the ion-permanent dipole enhanced Langevin capture rate with a predicted T-0.16 dependence. Electronic structure theory calculations were combined with Rice-Ramsperger-Kassel-Markus theory to extrapolate the measured recombination rate coefficients to the temperature and pressure conditions of planetary upper atmospheres. The following low-pressure limiting rate coefficients were obtained for T = 120-400 K and He bath gas (in cm6 molecule-2 s-1, uncertainty ±σ at 180 K): log10(k, Ni+ + N2) = -27.5009 + 1.0667log10(T) - 0.74741(log10(T))2, σ = 29%; log10(k, Ni+ + O2) = -27.8098 + 1.3065log10(T) - 0.81136(log10(T))2, σ = 32%; log10(k, Ni+ + CO2) = -29.805 + 4.2282log10(T) - 1.4303(log10(T))2, σ = 28%; log10(k, Ni+ + H2O) = -24.318 + 0.20448log10(T) - 0.66676(log10(T))2, σ = 28%). Other rate coefficients measured (at 294 K, in cm3 molecule-1 s-1) were: k(NiO+ + O) = (1.7 ± 1.2) × 10-10; k(NiO+ + CO) = (7.4 ± 1.3) × 10-11; k(NiO+ + O3) = (2.7 ± 1.0) × 10-10 with (29 ± 21)% forming Ni+ as opposed to NiO2+; k(NiO2+ + O3) = (2.9 ± 1.4) × 10-10, with (16 ± 9)% forming NiO+ as opposed to ONiO2+; and k(Ni+·N2 + O) = (7 ± 4) × 10-12. The chemistry of Ni+ and NiO+ in the upper atmospheres of Earth and Mars is then discussed.

Thermal activation of methane by MgO+: temperature dependent kinetics, reactive molecular dynamics simulations and statistical modeling

The kinetics methane activation (MgO⁺ + CH₄) was studied experimentally and computationally by running and analyzing reactive atomistic simulations.

On the Role of Hydrogen Atom Transfer (HAT) in Thermal Activation of Methane by MnO+: Entropy vs. Energy

The temperature dependent kinetics and product branching fractions of first-row transition metal oxide cation MnO+ with CH4 and CD4 at temperatures between 200 and 600 K are measured using a selected-ion flow tube apparatus. Likely reaction mechanisms are determined by comparison of temperature dependent kinetics to statistical modeling along calculated reaction coordinates. The data is well-modeled with the reaction proceeding over a rate limiting four-centered transition state leading to an insertion intermediate, similar to reactions of NiO+ and FeO+, and showing characteristics of proton-coupled electron transfer (PCET). However, a more direct pathway traversing a transition state of hydrogen atom transfer (HAT) character to a hydroxyl intermediate is found to possibly be competitive, especially with increasing temperature. While uncertainties in calculated energetics limit quantitative assessment of the role of HAT at thermal energies, it is clear that this mechanism becomes increasingly prevalent in higher energy regimes.