Quantum reactive scattering of H + hydrocarbon reactions

Electroassisted Propane Dehydrogenation at Low Temperatures: Far beyond the Equilibrium Limitation

Propylene production by propane dehydrogenation (PDH) generally requires high temperatures due to thermodynamic equilibrium limitations. This study developed a novel type of catalytic system for low-temperature PDH by combining a surface protonics methodology with intermetallic active sites. By application of an electric current, the intermetallic Pt-In/TiO2 catalyst gave a propylene yield of 10.2% with high selectivity, even at 250 °C, where the thermodynamic equilibrium yield was only 0.15%. Electroassisted proton collisions with propane allowed an unusual reaction pathway for low-temperature PDH. Alloying of Pt with In drastically enhanced the activity and selectivity due to the increased electron density of Pt.

Enhanced activity of catalysts on substrates with surface protonic current in an electrical field - a review

It has over the last few years been reported that the application of a DC electric field and resulting current over a bed of certain catalyst-support systems enhances catalytic activity for several reactions involving hydrogen-containing reactants, and the effect has been attributed to surface protonic conductivity on the porous ceramic support (typically ZrO2, CeO2, SrZrO3). Models for the nature of the interaction between the protonic current, the catalyst particle (typically Ru, Ni, Co, Fe), and adsorbed reactants such as NH3 and CH4 have developed as experimental evidence has emerged. Here, we summarize the electrical enhancement and how it enhances yield and lowers reaction temperatures of industrially important chemical processes. We also review the nature of the relevant catalysts, support materials, as well as essentials and recent progress in surface protonics. It is easily suspected that the effect is merely an increase in local vs. nominal set temperature due to the ohmic heating of the electrical field and current. We address this and add data from recent studies of ours that indicate that the heating effect is minor, and that the novel catalytic effect of a surface protonic current must have additional causes.

Calculations of quantum tunnelling rates for muonium reactions with methane, ethane and propane

Thermal rate constants for Mu + CH₄, Mu + C₂H₆ and Mu + C₃H₈ and their equivalent reactions with H were evaluated with ab initio instanton rate theory. The potential-energy surfaces are fitted using Gaussian process regression to high-level electronic-structure calculations evaluated around the tunnelling pathway. This method was able to successfully reproduce various experimental measurements for the rate constant of these reactions. However, it was not able to reproduce the faster-than-expected rate of Mu + C₃H₈ at 300 K reported by Fleming et al. [Phys. Chem. Chem. Phys., 2015, 17, 19901 and Phys. Chem. Chem. Phys., 2020, 22, 6326]. Analysis of our results indicates that the kinetic isotope effect at this temperature is not significantly influenced by quantum tunnelling. We consider many possible factors for the discrepancy between theory and experiment but conclude that in each case, the instanton approximation is unlikely to be the cause of the error. This is in part based on the good agreement we find between the instanton predictions and new multiconfigurational time-dependent Hartree (MCTDH) calculations for Mu + CH₄ using the same potential-energy surface. Further experiments will therefore be needed to resolve this issue.

Recent progress in use and observation of surface hydrogen migration over metal oxides

Hydrogen migration over a metal oxide surface is an extremely important factor governing the activity and selectivity of various heterogeneous catalytic reactions. Passive migration of hydrogen governed by a concentration gradient is called hydrogen spillover, which has been investigated broadly for a long time. Recently, well-fabricated samples and state-of-the-art measurement techniques such as operando spectroscopy and electrochemical analysis have been developed, yielding findings that have elucidated the migration mechanism and novel utilisation of hydrogen spillover. Furthermore, great attention has been devoted to surface protonics, which is hydrogen migration activated by an electric field, as applicable for novel low-temperature catalysis. This article presents an overview of catalysis related to hydrogen hopping, sophisticated analysis techniques for hydrogen migration, and low-temperature catalysis using surface protonics.

Rate constants and kinetic isotope effects for H-atom abstraction reactions by muonium in the Mu + propane and Mu + n-butane reactions from 300 K to 435 K: challenges for theory

This paper reports measurements of the temperature dependence of the rate constants for H-atom abstraction reactions from propane and n-butane by the light isotopic H-atom muonium (Mu), k_Mu(T), over temperatures in the range 300 K to 435 K.

Spiers Memorial Lecture : Introductory lecture: quantum dynamics of chemical reactions

This Spiers Memorial Lecture discusses quantum effects that can be calculated and observed in the chemical reactions of small molecules.

Surface Protonics Promotes Catalysis

Catalytic steam reforming of methane for hydrogen production proceeds even at 473 K over 1 wt% Pd/CeO₂ catalyst in an electric field, thanks to the surface protonics. Kinetic analyses demonstrated the synergetic effect between catalytic reaction and electric field, revealing strengthened water pressure dependence of the reaction rate when applying an electric field, with one-third the apparent activation energy at the lower reaction temperature range. Operando–IR measurements revealed that proton conduction via adsorbed water on the catalyst surface occurred during electric field application. Methane was activated by proton collision at the Pd–CeO₂ interface, based on the inverse kinetic isotope effect. Proton conduction on the catalyst surface plays an important role in methane activation at low temperature. This report is the first describing promotion of the catalytic reaction by surface protonics.

Testing the Palma-Clary Reduced Dimensionality Model Using Classical Mechanics on the CH4 + H → CH3 + H2 Reaction

Application of exact quantum scattering methods in theoretical reaction dynamics of bimolecular reactions is limited by the complexity of the equations of nuclear motion to be solved. Simplification is often achieved by reducing the number of degrees of freedom to be explicitly handled by freezing the less important spectator modes. The reaction cross sections obtained in reduced-dimensionality (RD) quantum scattering methods can be used in the calculation of rate coefficients, but their physical meaning is limited. The accurate test of the performance of a reduced-dimensionality method would be a comparison of the RD cross sections with those obtained in accurate full-dimensional (FD) calculations, which is not feasible because of the lack of complete full-dimensional results. However, classical mechanics allows one to perform reaction dynamics calculations using both the RD and the FD model. In this paper, an RD versus FD comparison is made for the 8-dimensional Palma-Clary model on the example of four isotopologs of the CH4 + H → CH3 + H2 reaction, which has 12 internal dimensions. In the Palma-Clary model, the only restriction is that the methyl group is confined to maintain C3v symmetry. Both RD and FD opacity and excitation functions as well as differential cross sections were calculated using the quasiclassical trajectory method. The initial reactant separation has been handled according to our one-period averaging method [ Nagy et al. J. Chem. Phys. 2016, 144, 014104 ]. The RD and FD excitation functions were found to be close to each other for some isotopologs, but in general, the RD reactivity parameters are lower than the FD reactivity parameters beyond statistical error, and for one of the isotopologs, the deviation is significant. This indicates that the goodness of RD cross sections cannot be taken for granted.

Accounting for conformational flexibility and torsional anharmonicity in the H + CH3CH2OH hydrogen abstraction reactions: a multi-path variational transition state theory study

This work reports a detailed theoretical study of the hydrogen abstraction reactions from ethanol by atomic hydrogen. The calculated thermal rate constants take into account torsional anharmonicity and conformational flexibility, in addition to the variational and tunneling effects. Specifically, the kinetics calculations were performed by using multi-path canonical variational transition state theory with least-action path tunneling corrections, to which we have added the two-dimensional non-separable method to take into account torsional anharmonicity. The multi-path thermal rate constant is expressed as a sum over conformational reaction channels. Each of these channels includes all the transition states that can be reached by internal rotations. The results show that, in the interval of temperatures between 250 and 2500 K, the account for multiple paths leads to higher thermal rate constants with respect to the single path approach, mainly at low and at high temperatures. In addition, torsional anharmonicity enhances the slope of the Arrhenius plot in this range of temperatures. Finally, we show that the incorporation of tunneling into the hydrogen abstraction reactions substantially changes the contribution of each of the transition states to the conformational reaction channel.

Rate constants for the slow Mu + propane abstraction reaction at 300 K by diamagnetic RF resonance

The rate constant for the slow Mu + propane abstraction reaction has been determined by diamagnetic RF resonance. The curves show simulations of the μSR resonance signal. This study provides an important new test of reaction rate theory for the alkanes.

Quantum effects in the abstraction reaction by H atoms of primary and secondary hydrogens in n-C4H10: a test of a new potential energy surface construction method

Recently, von Horsten et al. suggested an efficient method to construct two dimensional potential energy surfaces (PESs) for use in quantum scattering simulations, in which they utilised the minimum energy path (MEP) and assumed harmonic behaviour of the PES near the MEP. In the same paper, the authors applied this method to various H-abstraction reactions from C1-C3 alkane molecules. In this work we demonstrate an alternative PES construction method, and apply it to the more challenging H-abstraction from n-C(4)H(10) reactions. The geometry optimizations and frequency calculations are done at the MP2/cc-pVTZ level of theory, while the energies are calculated with the CCSD(T) method with the same basis set. The calculations give adiabatic energy barrier heights of 45.9 kJ mol(-1) and 34.4 kJ mol(-1) for the primary and secondary hydrogens in n-C(4)H(10). When compared to purely classical transition state theory, quantum scattering calculations show that quantum tunnelling and zero-point effects have large contributions at low temperatures, typically below 500 K. The branching ratio study suggests that the abstraction of secondary hydrogen in n-C(4)H(10) dominates the overall reaction rate at low temperatures. The rate constants for the two abstraction channels become more comparable as the temperature increases.

Ab initio methods for reactive potential surfaces

Case studies of ten reactions using a variety of standard electronic structure methods are presented. These case studies are used to illustrate the usefulness and shortcomings of these standard methods for various classes of reactions. Limited comparisons with experiment are made. The reactions studied include four radical-radical combinations, H + CH(3)--> CH(4), CH(3) + CH(3)--> C(2)H(6), H + HCO --> H(2)CO and CH(3) + HCO --> CH(3)CHO, three abstraction reactions, H + HO(2)--> H(2) + O(2), H + HCO --> H(2) + CO and CH(3) + HCO --> CH(4) + CO, a radical-molecule addition, H + HCCH --> C(2)H(3), and two molecular decompositions, H(2)CO --> H(2) + CO and CH(3)CHO --> CH(4) + CO. The electronic structure methods used are DFT, MP2, CCSD(T), QCISD(T), CASSCF, CASPT2, and CAS+1+2+QC.

Temperature, composition, and hydrogen isotope effect in the hydrogenation of CO on amorphous ice surface at 10–20K

An experiment on the addition reaction of a D atom (deuteration) to CO on a cold ice surface is performed by deuterium atom exposure of three types of samples (pure solid CO, CO-capped H2O ice, and CO-H2O mixed ice) at 10-20 K. The variation of IR absorption spectra for the samples was measured by a Fourier transform infrared spectrometer during exposure to deuterium atoms. Reactions on pure solid CO were observed only at 10 K, while reactions on CO-capped H2O ice and CO-H2O mixed ice were observed to proceed even at 20 K. This indicates that the coexistence of H2O at the surface raises the reactive temperature. In addition, the experiment on H atom exposure was also carried out at 15 K to compare the reaction rate constant between the H and D atoms. The ratio of reaction rate constant kD/kH obtained is about 0.08 at 15 K. The authors provide information on the potential energy for the H+CO reaction at the surface by using the ratio kD/kH and by a model calculation of the potential tunneling with the asymmetric Eckart potential.

Torsional anharmonicity in transition state theory calculations

We present a new application for the Torsional Path Integral Monte Carlo (TPIMC) method in which the TPI partition functions are introduced into the calculation of Transition State Theory (TST) rate constants. In this way, an explicit treatment of torsional anharmonicity is included in the TST calculations and the magnitude of these effects can be assessed. The new method is tested on the C(2)H(6) + H hydrogen abstraction reaction and concerted hydrogen transfer in the carbonic acid dimer, for which we have developed torsional potential energy surfaces. For the C(2)H(6) + H reaction the rate constants are halved at room temperature on including a treatment of torsional anharmonicity, while the effects are found to be much smaller for the hydrogen transfer reaction in the carbonic acid dimer.

Reduced dimensionality quantum dynamics of Cl + CH4? HCl + CH3 on an ab initio potential

We study the reaction Cl + CH(4)--> HCl + CH(3) using a 2-D potential energy surface obtained by fitting a double Morse analytical function to high level (CCSD(T)/cc-pVTZ//MP2/cc-pVTZ)ab initio data. Dynamics simulations are performed in hyperspherical coordinates with the close-coupled equations being solved using R-matrix propagation. Quantum contributions from spectator modes are included via a harmonic zero-point correction to the ab initio data prior to fitting the potential. This is the first time this method has been applied to a heavy-light-heavy reaction and the first time it has been used to study differential cross sections. We find thermal rate constants and state-to-state differential cross sections which are in good agreement with experimental data. We discuss the applicability of our method to the study of kinetic isotope effects (KIEs), which we derive for the CH(4)/CD(4) substitution. The calculated KIE compares favourably with experiment. Finally, we discuss the sensitivity of the results of dynamics simulations on the accuracy of the fitted potential.

Understanding the kinetics of spin-forbidden chemical reactions

Many chemical reactions involve a change in spin-state and are formally forbidden. This article summarises a number of previously published applications showing that a form of Transition State Theory (TST) can account for the kinetics of these reactions. New calculations for the emblematic spin-forbidden reaction HC + N(2) are also reported. The observed reactivity is determined by two factors. The first is the critical energy required for reaction to occur, which in spin-forbidden reactions is often defined by the relative energy of the Minimum Energy Crossing Point (MECP) between potential energy surfaces corresponding to the different spin states. The second factor is the probability of hopping from one surface to the other in the vicinity of the crossing region, which is largely defined by the spin-orbit coupling matrix element between the two electronic wavefunctions. The spin-forbidden transition state theory takes both factors into account and gives good results. The shortcomings of the theory, which are largely analogous to those of standard TST, are discussed. Finally, it is shown that in cases where the surface-hopping probability is low, the kinetics of spin-forbidden reactions will be characterised by unusually unfavourable entropies of activation. As a consequence, reactions involving a spin-state change can be expected to compete poorly with spin-allowed reactions at high temperatures (or energies).

Theories of reactive scattering

This paper is an overview of the theory of reactive scattering, with emphasis on fully quantum mechanical theories that have been developed to describe simple chemical reactions, especially atom-diatom reactions. We also describe related quasiclassical trajectory applications, and in all of this review the emphasis is on methods and applications concerned with state-resolved reaction dynamics. The review first provides an overview of the development of the theory, including a discussion of computational methods based on coupled channel calculations, variational methods, and wave packet methods. Choices of coordinates, including the use of hyperspherical coordinates are discussed, as are basis set and discrete variational representations. The review also summarizes a number of applications that have been performed, especially the two most comprehensively studied systems, H+H2 and F+H2, along with brief discussions of a large number of other systems, including other hydrogen atom transfer reactions, insertion reactions, electronically nonadiabatic reactions, and reactions involving four or more atoms. For each reaction we describe the method used and important new physical insight extracted from the results.