Significant quantum effects in hydrogen activation

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.

Investigating Spillover Energy as a Descriptor for Single-Atom Alloy Catalyst Design

The identification of thermodynamic descriptors of catalytic performance is essential for the rational design of heterogeneous catalysts. Here, we investigate how spillover energy, a descriptor quantifying whether intermediates are more stable at the dopant or host metal sites, can be used to design single-atom alloys (SAAs) for formic acid dehydrogenation. Using theoretical calculations, we identify NiCu as a SAA with favorable spillover energy and demonstrate that formate intermediates produced after the initial O-H activation are more stable at Ni sites where rate-determining C-H activation occurs. Surface science experiments demonstrated that NiCu(111) SAAs are more reactive than Cu(111) while they still follow the formate reaction pathway. However, reactor studies of silica-supported NiCu SAA nanoparticles showed only a modest improvement over Cu resulting from surface coverage effects. Overall, this study demonstrates the potential of engineering SAAs using spillover energy as a design parameter and highlights the importance of adsorbate-adsorbate interactions under steady-state operation.

Sustained Hydrogen Spillover on Pt/Cu(111) Single-Atom Alloy: Dynamic Insights into Gas-Induced Chemical Processes

Hydrogen spillover, involving the surface migration of dissociated hydrogen atoms from active metal sites to the relatively inert catalyst support, plays a crucial role in hydrogen-involved catalytic processes. However, a comprehensive understanding of how H atoms are driven to spill over from active sites onto the catalyst support is still lacking. Here, we examine the atomic-scale perspective of the H spillover process on a Pt/Cu(111) single atom alloy surface using machine-learning accelerated molecular dynamics calculations based on density functional theory. Our results show that when an impinging H2 dissociates at an active Pt site, the Pt atom undergoes deactivation due to the dissociated hydrogen atoms that attach to it. Interestingly, collisions between H2 and sticking H atoms facilitate H spillover onto the host Cu, leading to the reactivation of the Pt atom and the realization of a continuous H spillover process. This work underscores the importance of the interaction between gas molecules and adsorbates as a driving force in elucidating chemical processes under a gaseous atmosphere, which has so far been underappreciated in thermodynamic studies.

Geometric and Electronic Effects in the Binding Affinity of Imidazole-Based N-Heterocyclic Carbenes to Cu(100)- and Ag(100)-Based Pd and Pt Single-Atom Alloy Surfaces

We have conducted nonlocal periodic density functional theory (DFT) calculations of N-heterocyclic carbenes (NHCs) adsorbed to Pd/Cu(100), Pt/Cu(100), Pd/Ag(100), and Pt/Ag(100) single atom alloys (SAAs) utilizing the nonlocal optPBE-vdW functional. NHCs with electron donating groups (EDGs) are predicted to bind more strongly to the SAA surface compared to NHCs functionalized with electron withdrawing groups (EWGs). Our calculations show that NHCs typically bind to SAA geometries containing a small space between the heteroatom sites for the SAAs considered. Generally, this pattern is predicted to persist for a single NHCs or for a pair of NHCs bound to the SAA surfaces. Approximate linear relationships between NMR-based parameters and NHC-SAA binding energies are uncovered. We predict that the binding of NHCs to SAA surfaces is composition-dependent and heteroatom geometry dependent.

One Decade of Computational Studies on Single-Atom Alloys: Is In Silico Design within Reach?

ConspectusSingle-Atom alloys (SAAs) are an emerging class of materials consisting of a coinage metal (Cu, Ag, and Au) doped, at the single-atom limit, with another metal. As catalysts, coinage metals are rarely very active on their own, but when they are, they exhibit high selectivity. On the other hand, transition metals are usually very active but not as selective. Incorporating transition metals (guest elements) into coinage metals (host material) is therefore appealing for combining the activity and selectivity of each constituent in a balanced way. Additionally, first-principles calculations have shown that single atoms embedded in the surface of a coinage metal can exhibit emergent properties. Here, we describe how computational studies based on density functional theory (DFT) and kinetic Monte Carlo (KMC) simulations, often undertaken in close collaboration with experimental research groups, have shaped, over the past decade, the way we understand SAA catalysis.This Account reviews our contributions in elucidating the stability of SAAs, their electronic structure, and the way adsorbates interact and react on SAA catalytic surfaces. By studying in detail the processes that affect the stability of the SAA phase, we have shown that out of several bimetallic combinations of coinage metals with prominent Pt-group metals only PtCu and PdCu are stable surface alloys under vacuum. However, more surface alloy structures are possible in the presence of adsorbates because the latter can stabilize, via strong binding, dopants in the surface of the material. More interestingly, a large number of these surface alloys are resistant to the aggregation of dopant atoms into clusters, thereby favoring the SAA structure. These major results from DFT calculations serve as a guide for experimentalists to explore new SAA catalysts. Further analysis has shown that SAAs have a unique electronic structure with a very sharp d-band feature close to the Fermi level, analogous to the electronic structure of molecular entities. This is one of the reasons that SAAs are particularly sought after: although they are metallic nanoparticles, they have properties akin to those of homogeneous catalysts. In this context, we have contributed extensive screening studies, focusing on molecular fragments of catalytic relevance on a range of SAAs, which have driven the identification of new catalysts. We have also explored the rich chemistry of two-adsorbate systems via kinetic modeling, demonstrating how a spectator species with greater affinity for the dopant can modulate the reactivity of the catalyst via the so-called (punctured) molecular cork effect.Since the first experimental characterization of SAAs about a decade ago, theoretical models have been able to support and explain various experimental observations. These models have served as benchmarks for assessing the predictive capability of the underlying theoretical methods. In turn, the predictions that have been delivered have guided and continue to guide the experimental research efforts in the field. These advancements show that the in silico design of new SAA catalysts is now within reach.

Elucidation of the atomic-scale processes of dissociative adsorption and spillover of hydrogen on the single atom alloy catalyst Pd/Cu(111)

Hydrogen spillover is a crucial process in the selective hydrogenation reactions on Pd/Cu single atom alloy catalysts. In this study, we report the atomic-scale perspective of these processes on the...

Dynamics of Initial Hydrogen Spillover from a Single Atom Platinum Active Site to the Cu(111) Host Surface: The Impact of Substrate Electron-Hole Pairs

The initial impulsive diffusion of hot hydrogen atoms resulted from the dissociative chemisorption of H₂ at atomically dispersed Pt atoms embedded in Cu(111) is investigated using ab initio molecular dynamics. Upon dissociation, one of the two hydrogen atoms tends to roam away from the dissociation site while the other remains trapped. It is shown that the fraction of diffusion and the average diffusion length increase with the incident energy and H₂ vibrational excitation, due apparently to the increased initial kinetic energy of the hot atoms. Most importantly, the strong interaction with surface electron-hole pairs, modeled using an electronic friction model, is shown to play an important role in rapid energy dissipation and significant retardation of the impulsive diffusion.

The role of oxygenated species in the catalytic self-coupling of MeOH on O pre-covered Au(111)

The oxidation of alcohols plays a central role in the valorisation of biomass, in particular when performed with a non-toxic oxidant such as O₂. Aerobic oxidation of methanol on gold has attracted attention lately and the main steps of its mechanism have been described experimentally. However, the exact role of O and OH on each elementary step and the effect of the interactions between adsorbates are still not completely understood. Here we investigate the mechanism of methanol oxidation to HCOOCH₃ and CO₂. We use Density Functional Theory (DFT) to assess the energetics of the underlying pathways, and subsequently build lattice kinetic Monte Carlo (kMC) models of increasing complexity, to elucidate the role of different oxygenates. Detailed comparisons of our simulation results with experimental temperature programmed desorption (TPD) spectra enable us to validate the mechanism and identify rate determining steps. Crucially, taking into account dispersion (van der Waals forces) and adsorbate-adsorbate lateral interactions are both important for reproducing the experimental data.

Anharmonic calculations of vibrational spectra for molecular adsorbates: A divide-and-conquer semiclassical molecular dynamics approach

The vibrational spectroscopy of adsorbates is becoming an important investigation tool for catalysis and material science. This paper presents a semiclassical molecular dynamics method able to reproduce the vibrational energy levels of systems composed by molecules adsorbed on solid surfaces. Specifically, we extend our divide-and-conquer semiclassical method for power spectra calculations to gas-surface systems and interface it with plane-wave electronic structure codes. The Born-Oppenheimer classical dynamics underlying the semiclassical calculation is full dimensional, and our method includes not only the motion of the adsorbate but also those of the surface and the bulk. The vibrational spectroscopic peaks related to the adsorbate are accounted together with the most coupled phonon modes to obtain spectra amenable to physical interpretations. We apply the method to the adsorption of CO, NO, and H₂O on the anatase-TiO₂ (101) surface. We compare our semiclassical results with the single-point harmonic estimates and the classical power spectra obtained from the same trajectory employed in the semiclassical calculation. We find that CO and NO anharmonic effects of fundamental vibrations are similarly reproduced by the classical and semiclassical dynamics and that H₂O adsorption is fully and properly described in its overtone and combination band relevant components only by the semiclassical approach.

How coverage influences thermodynamic and kinetic isotope effects for H2/D2 dissociative adsorption on transition metals

Hydrogen isotope effects are influenced by adsorbate coverage: at high coverages, isotope effects are lower than at low coverages. This helps to rationalize observed isotope effects, allowing more precise elucidation of reaction mechanisms.

First principles rates for surface chemistry employing exact transition state theory: application to recombinative desorption of hydrogen from Cu(111)

We present first principles calculations of the reactive flux for thermal recombinative desorption of hydrogen from Cu(111).

Understanding the Role of Rutile TiO2 Surface Orientation on Molecular Hydrogen Activation

Titanium oxide (TiO₂) has been widely used in many fields, such as photocatalysis, photovoltaics, catalysis, and sensors, where its interaction with molecular H₂ with TiO₂ surface plays an important role. However, the activation of hydrogen over rutile TiO₂ surfaces has not been systematically studied regarding the surface termination dependence. In this work, we use density functional theory (PBE+U) to identify the pathways for two processes: the heterolytic dissociation of H₂ as a hydride–proton pair, and the subsequent H transfer from Ti to near O accompanied by reduction of the Ti sites. Four stoichiometric surface orientations were considered: (001), (100), (110), and (101). The lowest activation barriers are found for hydrogen dissociation on (001) and (110), with energies of 0.56 eV and 0.50 eV, respectively. The highest activation barriers are found on (100) and (101), with energies of 1.08 eV and 0.79 eV, respectively. For hydrogen transfer from Ti to near O, the activation barriers are higher (from 1.40 to 1.86 eV). Our results indicate that the dissociation step is kinetically more favorable than the H transfer process, although the latter is thermodynamically more favorable. We discuss the implications in the stability of the hydride–proton pair, and provide structures, electronic structure, vibrational analysis, and temperature effects to characterize the reactivity of the four TiO₂ orientations.

Novel catalytically active Pd/Ru bimetallic nanoparticles synthesized by Bacillus benzeovorans

Bacillus benzeovorans assisted and supported growth of ruthenium (bio-Ru) and palladium/ruthenium (bio-Pd@Ru) core@shell nanoparticles (NPs) as bio-derived catalysts. Characterization of the bio-NPs using various electron microscopy techniques and high-angle annular dark field (HAADF) analysis confirmed two NP populations (1–2 nm and 5–8 nm), with core@shells in the latter. The Pd/Ru NP lattice fringes, 0.231 nm, corresponded to the (110) plane of RuO_2. While surface characterization using X-ray photoelectron spectroscopy (XPS) showed the presence of Pd(0), Pd(II), Ru(III) and Ru(VI), X-ray absorption (XAS) studies of the bulk material confirmed the Pd speciation (Pd(0) and Pd(II)- corresponding to PdO), and identified Ru as Ru(III) and Ru(IV). The absence of Ru–Ru or Ru–Pd peaks indicated Ru only exists in oxide forms (RuO₂ and RuOH), which are surface-localized. X ray diffraction (XRD) patterns did not identify Pd-Ru alloying. Preliminary catalytic studies explored the conversion of 5-hydroxymethyl furfural (5-HMF) to the fuel precursor 2,5-dimethyl furan (2,5-DMF). Both high-loading (9.7 wt.% Pd, 6 wt.% Ru) and low-loading (2.4 wt.% Pd, 2 wt.% Ru) bio-derived catalysts demonstrated high conversion efficiencies (~95%) and selectivity of ~63% (~20% better than bio-Ru NPs) and 58%, respectively. These materials show promising future scope as efficient low-cost biofuel catalysts.

Single-Atom Alloys as a Reductionist Approach to the Rational Design of Heterogeneous Catalysts

Heterogeneous catalysts are workhorses in the industrial production of most commodity and specialty chemicals, and have widespread energy and environmental applications, with the annual market value of the catalysts themselves reaching almost $20 billion in 2018. These catalysts are complex, comprising multicomponent materials and multiple structures, making their rational design challenging, if not impossible. Furthermore, typical active metals like Pt, Pd, and Rh are expensive and can be susceptible to poisoning by CO, coking, and they are not always 100% selective. Efforts to use these elements sparingly and improve their selectivity has led to recent identification of single-atom heterogeneous catalysts in which individual transition metal atoms anchored on oxide or carbon-based supports are excellent catalysts for reactions like the CO oxidation, water-gas shift, alcohol dehydrogenation, and steam reforming. In this Account, we describe a new class of single-atom heterogeneous catalysts, namely, Single-Atom Alloys (SAAs) that comprise catalytically active elements like Pt, Pd, and Ni alloyed in more inert host metals at the single-atom limit. These materials evolved by complementary surface science and scanning probe studies using single crystals, and catalytic evaluation of the corresponding alloy nanoparticles with compositions informed by the surface science findings. The well-defined nature of the active sites in SAAs makes accurate modeling with theory relatively easy, enabling the rational design of SAA catalysts via a complementary three-prong approach, encompassing surface science model catalysts, theory, and real catalyst synthesis and testing under industrially relevant conditions. SAAs constitute one of just a few examples of when heterogeneous catalyst design has been guided by an understanding of fundamental surface processes. The Account starts by describing scanning tunneling microscopy studies of highly dilute alloys formed by doping small amounts of a catalytically active element into a more inert host metal. We first discuss hydrogenation reactions in which dissociation of H₂ is often rate limiting. Results indicate how the SAA geometry allows the transition state and the binding site of the reaction intermediates to be decoupled, which enables both facile dissociation of reactants and weak binding of intermediates, two key factors for efficient and selective catalysis. These results were exploited to design the first PtCu SAA hydrogenation catalysts which showed high selectivity, stability and resistance to poisoning in industrially relevant hydrogenation reactions, such as the selective conversion of butadiene to butenes. Model studies also revealed spillover of hydrogen atoms from the Pt site where dissociation of H₂ occurs to Cu sites where selective hydrogenation is facilitated in a bifunctional manner. We then discuss selective dehydrogenations on SAAs demonstrating that they enable efficient C-H activation, while being resistant to coking that plagues typical Pt catalysts. SAA PtCu nanoparticle catalysts showed excellent stability in butane dehydrogenation for days-on-stream at 400 °C. Another advantage of SAA catalysts is that on many alloy combinations CO, a common catalyst poison, binds more weakly to the alloy than the pure metal. We conclude by discussing recent theory results that predict the energetics of many key reaction steps on a wide range of SAAs and the exciting possibilities this reductionist approach to heterogeneous catalysis offers for the rational design of new catalysts.

Lonely Atoms with Special Gifts: Breaking Linear Scaling Relationships in Heterogeneous Catalysis with Single-Atom Alloys

We discuss a simple yet effective strategy for escaping traditional linear scaling relations in heterogeneous catalysis with highly dilute bimetallic alloys known as single-atom alloys (SAAs). These systems, in which a reactive metal is atomically dispersed in a less reactive host, were first demonstrated with the techniques of surface science to be active and selective for hydrogenation reactions. Informed by these early results, PdCu and PtCu SAA nanoparticle hydrogenation catalysts were shown to work under industrially relevant conditions. To efficiently survey the many potential metal combinations and reactions, simulation is crucial for making predictions about reactivity and guiding experimental focus on the most promising candidate materials. This recent work reveals that the high surface chemical heterogeneity of SAAs can result in significant deviations from Brønsted-Evans-Polanyi scaling relationships for many key reaction steps. These recent insights into SAAs and their ability to break linear scaling relations motivate discovery of novel alloy catalysts.

Decisive role of nuclear quantum effects on surface mediated water dissociation at finite temperature

Water molecules adsorbed on inorganic substrates play an important role in several technological applications. In the presence of light atoms in adsorbates, nuclear quantum effects (NQEs) influence the structural stability and the dynamical properties of these systems. In this work, we explore the impact of NQEs on the dissociation of water wires on stepped Pt(221) surfaces. By performing ab initio molecular dynamics simulations with van der Waals corrected density functional theory, we note that several competing minima for both intact and dissociated structures are accessible at finite temperatures, making it important to assess whether harmonic estimates of the quantum free energy are sufficient to determine the relative stability of the different states. We thus perform ab initio path integral molecular dynamics (PIMD) in order to calculate these contributions taking into account the conformational entropy and anharmonicities at finite temperatures. We propose that when adsorption is weak and NQEs on the substrate are negligible, PIMD simulations can be performed through a simple partition of the system, resulting in considerable computational savings. We then calculate the full contribution of NQEs to the free energies, including also anharmonic terms. We find that they result in an increase of up to 20% of the quantum contribution to the dissociation free energy compared with the harmonic estimates. We also find that the dissociation process has a negligible contribution from tunneling but is dominated by zero point energies, which can enhance the rate of dissociation by three orders of magnitude. Finally we highlight how both temperature and NQEs indirectly impact dipoles and the redistribution of electron density, causing work function changes of up to 0.4 eV with respect to static estimates. This quantitative determination of the change in the work function provides a possible approach to determine experimentally the most stable configurations of water oligomers on the stepped surfaces.

Density functional theory study of thermodynamic and kinetic isotope effects of H2/D2 dissociative adsorption on transition metals

Thermodynamic/kinetic isotope effects for H₂/D₂ dissociative adsorption calculated on metal surfaces offer a means to identify active sites.

Controlling Hydrogen Activation, Spillover, and Desorption with Pd-Au Single-Atom Alloys

Key descriptors in hydrogenation catalysis are the nature of the active sites for H2 activation and the adsorption strength of H atoms to the surface. Using atomically resolved model systems of dilute Pd-Au surface alloys and density functional theory calculations, we determine key aspects of H2 activation, diffusion, and desorption. Pd monomers in a Au(111) surface catalyze the dissociative adsorption of H2 at temperatures as low as 85 K, a process previously expected to require contiguous Pd sites. H atoms preside at the Pd sites and desorb at temperatures significantly lower than those from pure Pd (175 versus 310 K). This facile H2 activation and weak adsorption of H atom intermediates are key requirements for active and selective hydrogenations. We also demonstrate weak adsorption of CO, a common catalyst poison, which is sufficient to force H atoms to spill over from Pd to Au sites, as evidenced by low-temperature H2 desorption.

Quantum and classical dynamics of reactive scattering of H2 from metal surfaces

State-of-the-art theoretical models allow nowadays an accurate description of H₂/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.

Toward a Database of Chemically Accurate Barrier Heights for Reactions of Molecules with Metal Surfaces

Being able to calculate reaction barrier heights to within chemical accuracy (errors < 1 kcal/mol) is crucial to the accurate modeling of chemical reactions. Although accurate databases exist that can help theorists with benchmarking new electronic structure theories on gas-phase chemical reactions, no such databases exist for reactions of molecules with metal surfaces. Nonetheless, most chemicals are made in heterogeneously catalyzed processes, of which many take place over metal particles. Presently, barrier heights for molecule-metal surface reactions have been determined with chemical accuracy for only two systems, that is, H2 + Cu(111) and H2 + Cu(100). This has been done with semiempirically determined density functionals, which were fitted through comparisons of dynamics results with molecular beam sticking probabilities. The prospects of extending the database with chemically accurate data for other molecule-metal reactions, either with the use of semiempirical density functional theory or with first-principles theory, are discussed.

Kinetic modelling of heterogeneous catalytic systems

The importance of heterogeneous catalysis in modern life is evidenced by the fact that numerous products and technologies routinely used nowadays involve catalysts in their synthesis or function. The discovery of catalytic materials is, however, a non-trivial procedure, requiring tedious trial-and-error experimentation. First-principles-based kinetic modelling methods have recently emerged as a promising way to understand catalytic function and aid in materials discovery. In particular, kinetic Monte Carlo (KMC) simulation is increasingly becoming more popular, as it can integrate several sources of complexity encountered in catalytic systems, and has already been used to successfully unravel the underlying physics of several systems of interest. After a short discussion of the different scales involved in catalysis, we summarize the theory behind KMC simulation, and present the latest KMC computational implementations in the field. Early achievements that transformed the way we think about catalysts are subsequently reviewed in connection to latest studies of realistic systems, in an attempt to highlight how the field has evolved over the last few decades. Present challenges and future directions and opportunities in computational catalysis are finally discussed.