Quantum Tunnelling Driven H2 Formation on Graphene

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

Nonadiabatic Tunneling in Chemical Reactions

Quantum tunneling can have a dramatic effect on chemical reaction rates. In nonadiabatic reactions such as electron transfers or spin crossovers, nuclear tunneling effects can be even stronger than for adiabatic proton transfers. Ring-polymer instanton theory enables molecular simulations of tunneling in full dimensionality and has been shown to be far more reliable than commonly used separable approximations. First-principles instanton calculations predict significant nonadiabatic tunneling of heavy atoms even at room temperature and give excellent agreement with experimental measurements for the intersystem crossing of two nitrenes in cryogenic matrix isolation, the spin-forbidden relaxation of photoexcited thiophosgene in the gas phase, and singlet oxygen deactivation in water at ambient conditions. Finally, an outlook of further theoretical developments is discussed.

Robust Gaussian Process Regression Method for Efficient Tunneling Pathway Optimization: Application to Surface Processes

Simulation of surface processes is a key part of computational chemistry that offers atomic-scale insights into mechanisms of heterogeneous catalysis, diffusion dynamics, and quantum tunneling phenomena. The most common theoretical approaches involve optimization of reaction pathways, including semiclassical tunneling pathways (called instantons). The computational effort can be demanding, especially for instanton optimizations with an ab initio electronic structure. Recently, machine learning has been applied to accelerate reaction-pathway optimization, showing great potential for a wide range of applications. However, previous methods still suffer from numerical and efficiency issues and were not designed for condensed-phase reactions. We propose an improved framework based on Gaussian process regression for general transformed coordinates, which has improved efficiency and numerical stability, and we propose a descriptor that combines internal and Cartesian coordinates suitable for modeling surface processes. We demonstrate with 11 instanton optimizations in three representative systems that the improved approach makes ab initio instanton optimization significantly cheaper, such that it becomes not much more expensive than a classical transition-state theory rate calculation.

A personal perspective of the present status and future challenges facing thermal reaction rate theory

Reaction rate theory has been at the center of physical chemistry for well over one hundred years. The evolution of the theory is not only of historical interest. Reliable and accurate computation of reaction rates remains a challenge to this very day, especially in view of the development of quantum chemistry methods, which predict the relevant force fields. It is still not possible to compute the numerically exact rate on the fly when the system has more than at most a few dozen anharmonic degrees of freedom, so one must consider various approximate methods, not only from the practical point of view of constructing numerical algorithms but also on conceptual and formal levels. In this Perspective, I present some of the recent analytical results concerning leading order terms in an ℏ2m series expansion of the exact rate and their implications on various approximate theories. A second aspect has to do with the crossover temperature between tunneling and thermal activation. Using a uniform semiclassical transmission probability rather than the "primitive" semiclassical theory leads to the conclusion that there is no divergence problem associated with a "crossover temperature." If one defines a semiclassical crossover temperature as the point at which the tunneling energy of the instanton equals the barrier height, then it is a factor of two higher than its previous estimate based on the "primitive" semiclassical approximation. In the low temperature tunneling regime, the uniform semiclassical theory as well as the "primitive" semiclassical theory were based on the classical Euclidean action of a periodic orbit on the inverted potential. The uniform semiclassical theory wrongly predicts that the "half-point," which is the energy at which the transmission probability equals 1/2, for any barrier potential, is always the barrier energy. We describe here how augmenting the Euclidean action with constant terms of order ℏ2 can significantly improve the accuracy of the semiclassical theory and correct this deficiency. This also leads to a deep connection with and improvement of vibrational perturbation theory. The uniform semiclassical theory also enables an extension of the quantum version of Kramers' turnover theory to temperatures below the "crossover temperature." The implications of these recent advances on various approximate methods used to date are discussed at length, leading to the conclusion that reaction rate theory will continue to challenge us both on conceptual and practical levels for years to come.

Engineering nanoscale H supply chain to accelerate methanol synthesis on ZnZrOx

Metal promotion is the most widely adopted strategy for enhancing the hydrogenation functionality of an oxide catalyst. Typically, metal nanoparticles or dopants are located directly on the catalyst surface to create interfacial synergy with active sites on the oxide, but the enhancement effect may be compromised by insufficient hydrogen delivery to these sites. Here, we introduce a strategy to promote a ZnZrO_x methanol synthesis catalyst by incorporating hydrogen activation and delivery functions through optimized integration of ZnZrO_x and Pd supported on carbon nanotube (Pd/CNT). The CNT in the Pd/CNT + ZnZrO_x system delivers hydrogen activated on Pd to a broad area on the ZnZrO_x surface, with an enhancement factor of 10 compared to the conventional Pd-promoted ZnZrO_x catalyst, which only transfers hydrogen to Pd-adjacent sites. In CO₂ hydrogenation to methanol, Pd/CNT + ZnZrO_x exhibits drastically boosted activity-the highest among reported ZnZrO_x-based catalysts-and excellent stability over 600 h on stream test, showing potential for practical implementation.

First-Principles Microkinetic Modeling Unravelling the Performance of Edge-Decorated Nanocarbons for Hydrogen Production from Methane

The doping of graphitic and nanocarbon structures with nonmetal atoms allows for the tuning of surface electronic properties and the generation of new active sites, which can then be exploited for several catalytic applications. In this work, we investigate the direct conversion of methane into H₂ and C₂H_x over Klein-type zigzag graphene edges doped with nitrogen, boron, phosphorus and silicon. We combine Density Functional Theory (DFT) and microkinetic modeling to systematically investigate the reaction network and determine the most efficient edge decoration. Among the four edge-decorated nanocarbons (EDNCs) investigated, N-EDNC presented an outstanding performance for H₂ production at temperatures over 900 K, followed by P-EDNC, Si-EDNC and B-EDNC. The DFT and microkinetic analysis of the enhanced desorption rate of atomic hydrogen reveal the presence of an Eley-Rideal mechanism, in which P-EDNC showed higher activity for H₂ production in this scenario. Coke deposition resistance in the temperature range between 900 and 1500 K was evaluated, and we compared the selectivity toward H₂ and C₂H₄ production. The N-EDNC and P-EDNC active sites showed strong resistance to carbon poisoning, whereas Si-EDNC showed higher propensity to regenerate its active sites at temperatures over 1100 K. This work shows that decorated EDNCs are promising metal-free catalysts for methane conversion into H₂ and short-length alkenes.

Heavy-Atom Quantum Tunnelling in Spin Crossovers of Nitrenes

We simulate two recent matrix-isolation experiments at cryogenic temperatures, in which a nitrene undergoes spin crossover from its triplet state to a singlet state via quantum tunnelling. We detail the failure of the commonly applied weak-coupling method (based on a linear approximation of the potentials) in describing these deep-tunnelling reactions. The more rigorous approach of semiclassical golden-rule instanton theory in conjunction with double-hybrid density-functional theory and multireference perturbation theory does, however, provide rate constants and kinetic isotope effects in good agreement with experiment. In addition, these calculations locate the optimal tunnelling pathways, which provide a molecular picture of the reaction mechanism. The reactions involve substantial heavy-atom quantum tunnelling of carbon, nitrogen and oxygen atoms, which unexpectedly even continues to play a role at room temperature.

Nonadiabatic instanton rate theory beyond the golden-rule limit

Fermi's golden rule describes the leading-order behaviour of the reaction rate as a function of the diabatic coupling. Its asymptotic (ℏ →0) limit is the semiclassical golden-rule instanton rate theory, which rigorously approximates nuclear quantum effects, lends itself to efficient numerical computation and gives physical insight into reaction mechanisms. However the golden rule by itself becomes insufficient as the strength of the diabatic coupling increases, so higher-order terms must be additionally considered. In this work we give a first-principles derivation of the next-order term beyond the golden rule, represented as a sum of three components. Two of them lead to new instanton pathways that extend the golden-rule case and, among other factors, account for the effects of recrossing on the full rate. The remaining component derives from the equilibrium partition function and accounts for changes in potential energy around the reactant and product wells due to diabatic coupling. The new semiclassical theory demands little computational effort beyond a golden-rule instanton calculation. It makes it possible to rigorously assess the accuracy of the golden-rule approximation and sets the stage for future work on general semiclassical nonadiabatic rate theories.