Ethylene conversion to ethylidyne on Pd(111) and Pt(111): A first-principles-based kinetic Monte Carlo study

Modified Activation-Relaxation Technique (ARTn) Method Tuned for Efficient Identification of Transition States in Surface Reactions

Exploring potential energy surfaces (PES) is essential for unraveling the underlying mechanisms of chemical reactions and material properties. While the activation-relaxation technique (ARTn) is a state-of-the-art method for identifying saddle points on PES, it often faces challenges in complex energy landscapes, especially on surfaces. In this study, we introduce iso-ARTn, an enhanced ARTn method that incorporates constraints on an orthogonal hyperplane and employs an adaptive active volume. By leveraging a neural network potential (NNP) to conduct an exhaustive saddle point search on the Pt(111) surface with 0.3 monolayers of surface oxygen coverage, iso-ARTn achieves a success rate that is 8.2% higher than the original ARTn, with 40% fewer force calls. Moreover, this method effectively finds various saddle points without compromising the success rate. Combined with kinetic Monte Carlo simulations for event table construction, iso-ARTn with NNP demonstrates the capability to reveal structures consistent with experimental observations. This work signifies a substantial advancement in the investigation of PES, enhancing both the efficiency and breadth of saddle point searches.

Bridging the size gap between experiment and theory: large-scale DFT calculations on realistic sized Pd particles for acetylene hydrogenation

Metal nanoparticles, often supported on metal oxide promoters, are a cornerstone of heterogeneous catalysis. Experimentally, size effects are well-established and are manifested through changes to catalyst selectivity, activity and durability. Density Functional Theory (DFT) calculations have provided an attractive way to study these effects and rationalise the change in nanoparticle properties. However such computational studies are typically limited to smaller nanoparticles (approximately up to 50 atoms) due to the large computational cost of DFT. How well can such simulations describe the electronic properties of the much larger nanoparticles that are often used in practice? In this study, we use the ONETEP code, which is able to achieve more favourable computational scaling for metallic nanoparticles, to bridge this size gap. We present DFT calculations on entire Pd and Pd carbide nanoparticles of more than 300 atoms (approximately 2.5 nm diameter), and find major differences in the electronic structure of such large nanoparticles, in comparison to the commonly investigated smaller clusters. These differences are also manifested in the calculated chemical properties such as adsorption energies for C2H2, C2H4 and C2H6 on the pristine Pd and PdC x nanoparticles which are significantly larger (up to twice in value) for the ∼300 atoms structures. Furthermore, the adsorption of C2H2 and C2H4 on PdC x nanoparticles becomes weaker as more C is introduced in the Pd lattice whilst the impact of C concentration is also observed in the calculated reaction energies towards the hydrogenation of C2H2, where the formation of C2H6 is hindered. Our simulations show that PdC x nanoparticles of about 5% C per atom fraction and diameter of 2.5 nm could be potential candidate catalysts of high activity in hydrogenation reactions. The paradigm presented in this study will enable DFT to be applied on similar sized metal catalyst nanoparticles as in experimental investigations, strengthening the synergy between simulation and experiment in catalysis.

Selective formation of acetate intermediate prolongs robust ethylene removal at 0 °C for 15 days

Efficient ethylene (C₂H₄) removal below room temperatures, especially near 0  °C, is of great importance to suppress that the vegetables and fruits spoil during cold-chain transportation and storage. However, no catalysts have been developed to fulfill the longer-than-2-h C₂H₄ removal at this low temperature effectively. Here we prepare gold-platinum (Au-Pt) nanoalloy catalysts that show robust C₂H₄ (of 50 ppm) removal capacity at 0 °C for 15 days (360 h). We find, by virtue of operando Fourier transformed infrared spectroscopy and online temperature-programmed desorption equipped mass spectrometry, that the Au-Pt nanoalloys favor the formation of acetate from selective C₂H₄ oxidation. And this on-site-formed acetate intermediate would partially cover the catalyst surface at 0 °C, thus exposing active sites to prolong the continuous and effective C₂H₄ removal. We also demonstrate, by heat treatment, that the performance of the used catalysts will be fully recovered for at least two times.

Atomistic simulations on the carbidisation processes in Pd nanoparticles

The formation of interstitial PdC _x nanoparticles (NPs) is investigated through DFT calculations. Insights on the mechanisms of carbidisation are obtained whilst the material's behaviour under conditions of increasing C-concentration is examined. Incorporation of C atoms in the Pd octahedral interstitial sites is occurring through the [111] facet with an activation energy barrier of 19.3-35.7 kJ mol^-1 whilst migration through the [100] facet corresponds to higher activation energy barriers of 124.5-127.4 kJ mol^-1. Furthermore, interstitial-type diffusion shows that C will preferentially migrate and reside at the octahedral interstitial sites in the subsurface region with limited mobility towards the core of the NP. For low C-concentrations, migration from the surface into the interstitial sites of the NPs is thermodynamically favored, resulting in the formation of interstitial carbide. Carbidisation reaction energies are exothermic up to 11-14% of C-concentration and slightly vary depending on the shape of the structure. The reaction mechanisms turn to endothermic for higher concentration levels showing that C will preferentially reside on the surface making the interstitial carbide formation unfavorable. As experimentally observed, our simulations confirm that there is a maximum concentration of C in Pd carbide NPs opening the way for further computational investigations on the activity of Pd carbides in directed catalysis.

Radical-Mediated C-C Coupling of Alcohols Induced by Plasmonic Hot Carriers

The C-C coupling reactions of aliphatic alcohols to aromatics and larger-mass compounds have large endothermicities and activation energies, calling for catalysts operating at high temperatures. Here, we demonstrate that plasmon-excited nanoparticles catalyze the C-C coupling of aliphatic alcohols at room temperature to produce polyaromatic hydrocarbons and graphene oxide. The conversion is quenched by radical and electron scavengers and by the surface passivation of metals, suggesting that the reaction proceeds through alkoxy, peroxyl, hydroxyalkyl, and alkyl radical intermediates created by the metal to molecule transfer of plasmonic hot carriers. Besides being the first realization of C-C coupling of aliphatic alcohols at room temperature, the result constitutes a rare example of an endothermic plasmon-induced reaction producing new bonds and a new method for photogenerating graphene derivatives. More importantly, the result demonstrates the facile generation of organic radicals directly from alcohols, which may be used as precursors for radical-based organic reactions.

Understanding the origin of structure sensitivity in hydrodechlorination of trichloroethylene on a palladium catalyst

Mechanistic understanding on the origin of structure sensitivity in the hydrodechlorination of trichloroethylene.

Controlling Hydrocarbon (De)Hydrogenation Pathways with Bifunctional PtCu Single-Atom Alloys

The conversions of surface-bound alkyl groups to alkanes and alkenes are important steps in many heterogeneously catalyzed reactions. On the one hand, while Pt is ubiquitous in industry because of its high activity toward C-H activation, many Pt-based catalysts tend to overbind reactive intermediates, which leads to deactivation by carbon deposition and coke formation. On the other hand, Cu binds intermediates more weakly than Pt, but activation barriers tend to be higher on Cu. We examine the reactivity of ethyl, the simplest alkyl group that can undergo hydrogenation and dehydrogenation via β-elimination, and show that isolated Pt atoms in Cu enable low-temperature hydrogenation of ethyl, unseen on Cu, while avoiding the decomposition pathways on pure Pt that lead to coking. Furthermore, we confirm the predictions of our theoretical model and experimentally demonstrate that the selectivity of ethyl (de)hydrogenation can be controlled by changing the surface coverage of hydrogen.

Dehydrogenation of Ethylene on Supported Palladium Nanoparticles: A Double View from Metal and Hydrocarbon Sides

Adsorption of ethylene on palladium, a key step in various catalytic reactions, may result in a variety of surface-adsorbed species and formation of palladium carbides, especially under industrially relevant pressures and temperatures. Therefore, the application of both surface and bulk sensitive techniques under reaction conditions is important for a comprehensive understanding of ethylene interaction with Pd-catalyst. In this work, we apply in situ X-ray absorption spectroscopy, X-ray diffraction and infrared spectroscopy to follow the evolution of the bulk and surface structure of an industrial catalysts consisting of 2.6 nm supported palladium nanoparticles upon exposure to ethylene under atmospheric pressure at 50 °C. Experimental results were complemented by ab initio simulations of atomic structure, X-ray absorption spectra and vibrational spectra. The adsorbed ethylene was shown to dehydrogenate to C2H3, C2H2 and C2H species, and to finally decompose to palladium carbide. Thus, this study reveals the evolution pathway of ethylene on industrial Pd-catalyst under atmospheric pressure at moderate temperatures, and provides a conceptual framework for the experimental and theoretical investigation of palladium-based systems, in which both surface and bulk structures exhibit a dynamic nature under reaction conditions.

Efficient Implementation of Cluster Expansion Models in Surface Kinetic Monte Carlo Simulations with Lateral Interactions: Subtraction Schemes, Supersites, and the Supercluster Contraction

While lateral interaction models for reactions at surfaces have steadily gained popularity and grown in terms of complexity, their use in chemical kinetics has been impeded by the low performance of current kinetic Monte Carlo (KMC) algorithms. The origins of the additional computational cost in KMC simulations with lateral interactions are traced back to the more elaborate cluster expansion Hamiltonian, the more extensive rate updating, and to the impracticality of rate-catalog-based algorithms for interacting adsorbate systems. Favoring instead site-based algorithms, we propose three ways to reduce the cost of KMC simulations: (1) representing the lattice energy by a smaller Supercluster Hamiltonian without loss of accuracy, (2) employing the subtraction schemes for updating key quantities in the simulation that undergo only small, local changes during a reaction event, and (3) applying efficient search algorithms from a set of established methods (supersite approach). The cost of the resulting algorithm is fixed with respect to the number of lattice sites for practical lattice sizes and scales with the square of the range of lateral interactions. The overall added cost of including a complex lateral interaction model amounts to less than a factor 3. Practical issues in implementation due to finite numerical accuracy are discussed in detail, and further suggestions for treating long-range lateral interactions are made. We conclude that, while KMC simulations with complex lateral interaction models are challenging, these challenges can be overcome by modifying the established variable step-size method by employing the supercluster, subtraction, and supersite algorithms (SSS-VSSM). © 2019 Wiley Periodicals, Inc.

Hydrocarbon decomposition kinetics on the Ir(111) surface

The kinetics of the thermal decomposition of hydrocarbons on the Ir(111) surface is determined using kinetic Monte Carlo (kMC) and rate equations simulations, both based on the density functional theory (DFT) calculated energy barriers of the involved reaction processes. This decomposition process is important for understanding the early stages of epitaxial graphene growth where the deposited hydrocarbon acts as a carbon feedstock for graphene formation. The methodology of the kMC simulations and the rate equation approaches is discussed and a comparison between the results obtained from both approaches is made in the case of the temperature programmed decomposition of ethylene for different initial coverages. The theoretical results are verified against experimental data from in situ X-ray photoelectron spectroscopy (XPS) experiments. Both theoretical approaches give reasonable results; however we find that, as expected, rate equations are less reliable at high coverages. We find that the agreement between experiment and theory can be improved in all cases if slight adjustments are made to the energy barriers in order to account for the intrinsic errors in DFT. Finally we extend our approach to the case where hydrocarbon species are dosed onto the substrate continuously, as in the chemical vapour deposition (CVD) graphene growth method. For ethylene and methane the thermal decomposition mechanism is determined, and it is found that in both cases the formation of C monomers is to be expected, which is limited by the presence of hydrogen atoms.

Approaching complexity of alkyl hydrogenation on Pd via density-functional modelling

Pd is widely used to catalyse hydrogenation and dehydrogenation reactions. One of them is the hydrogenation of ethylene, which includes the transformation of ethyl species to ethane. Herein, by means of density-functional calculations we address several still insufficiently understood factors affecting the latter process. In particular, we shed light on the following aspects of hydrogenation of alkyls on Pd: (i) the mechanistic details of how subsurface H accelerates the reaction on a (111) surface; (ii) the role of nanoparticle edges; and (iii) the influence of a common spectator ethylidyne, [triple bond, length as m-dash]C-CH₃. These factors are identified as significant for the height of the ethyl hydrogenation barrier on Pd. Moreover, we show that butyl hydrogenation on Pd is also governed by very similar interactions, which suggests a broader applicability of our conclusions. This study highlights the complexity of alkyl hydrogenation and analyses the factors that need to be taken into account for a more realistic description of the hydrogenation processes on metal surfaces.

Mechanistic Insights into Ethylene Transformations on Ir(111) by Density Functional Calculations and Microkinetic Modeling

Ethylidyne, ethane, and carbon monomer formations from ethylene over Ir(111) at different coverages are investigated using density functional theory methods. Two possible reaction mechanisms for ethylidyne formation are investigated. The calculations show that vinyl prefers the dehydrogenation to yield vinylidene (M2) over the hydrogenation to produce ethylidene (M1) kinetically and thermodynamically at 1/9 (1/3) ML. Ethylidyne formation could be a competitive side reaction of ethylene hydrogenation, however, the ethylidyne species does not directly participate in the ethylene hydrogenation mechanism. The mechanism for C monomer formation is also studied. Microkinetic modeling shows that the ethylene hydrogenation reactivity decreases in the sequence Ir(111)>Rh(111)>Pd(111)>Pt(111) under typical hydrogenation conditions. The catalytic activity of ethylene hydrogenation decreases with increased stability of ethylene adsorption and reaction barrier of the rate-limiting step.

Investigation on the conversion of ethylene to ethylidyne on Pt(100) and Pd(100) using density functional theory

The comprehensive formation network of ethylidyne (CH₃C) from ethylene (CH₂CH₂) is investigated on Pt(100) and Pd(100) using the density functional theory method. The structural and energetic features of all intermediate products were considered. We found that the trend of the activation barriers in each pathway on Pt(100) and Pd(100) are the same, whereas the barriers on Pt(100) are higher than that on Pd(100). The activation barriers of 1,2-H shift reactions are relatively high compared with the other reactions. We screened three possible pathways and selected the optimal route as CH₂CH₂(ethylene) → CH₂CH(vinyl) → CH₂C(vinylidene) → CH₃C(ethylidyne).

Ethylene decomposition on Ir(111): initial path to graphene formation

The complete mechanism behind the thermal decomposition of ethylene (C₂H₄) on Ir(111), which is the first step of graphene growth, is established for the first time employing a combination of experimental and theoretical methods.

Surface carbon species formation from ethylene decomposition on Pd(100): a first-principles-based kinetic Monte Carlo study

Based on the activation barriers and reaction energies from periodic density functional calculations, we conducted kinetic Monte Carlo (kMC) simulations of surface carbon species formation from ethylene decomposition on a Pd(100) surface.

Molecular understandings on the activation of light hydrocarbons over heterogeneous catalysts

Due to the depletion of petroleum and the recent shale gas revolution, the dropping of the price for light alkanes makes alkanes an attractive feedstock for the production of light alkenes and other valuable chemicals. Understanding the mechanism for the activation of C-H bonds in hydrocarbons provides fundamental insights into this process and a guideline for the optimization of catalysts used for the processing of light alkanes. In the last two decades, density functional theory (DFT) has become a powerful tool to explore elementary steps and mechanisms of many heterogeneously catalyzed processes at the atomic scale. This review describes recent progress on computational understanding of heterogeneous catalytic dehydrogenation reactions of light alkanes. We start with a short description on basic concepts and principles of DFT as well as its application in heterogeneous catalysis. The activation of C-H bonds over transition metal and alloy surfaces are then discussed in detail, followed by C-H activation over oxides, zeolites and catalysts with single atoms as active sites. The origins of coking formation are also discussed followed by a perspective on directions of future research.

Self-assembly and chemical reactivity of alkenes on platinum nanoparticles

Stable platinum nanoparticles were synthesized by the self-assembly of alkene derivatives onto the platinum surface, possibly forming platinum-vinylidene (Pt═C═CH-) or -acetylide (Pt-C≡) interfacial bonds as a result of dehydrogenation and transformation of the olefin moieties catalyzed by platinum. Transmission electron microscopic measurements showed that the nanoparticles were well-dispersed without apparent agglomeration, indicating effective passivation of the nanoparticles by the ligands, and the average core was estimated to be 1.34 ± 0.39 nm. FTIR measurements showed the emergence of a new vibrational band at 2023 cm(-1), which was ascribed to the formation of Pt-H and C≡C from the dehydrogenation of alkene ligands on platinum surfaces. Consistent behaviors were observed in photoluminescence measurements, where the emission profiles were similar to those of alkyne-functionalized Pt nanoparticles that arose from intraparticle charge delocalization between the particle-bound acetylene moieties. Selective reactivity with imine derivatives further confirmed the formation of Pt═C═CH- or Pt-C≡ interfacial linkages, as manifested in NMR and electrochemical measurements. Further structural insights were obtained by X-ray absorption near-edge spectroscopy and extended X-ray absorption fine structure analysis, where the coordinate numbers and bond lengths of the Pt-Pt and Pt-C linkages suggested that the metal-ligand interfacial bonds were in the intermediate between those of Pt-C≡ and Pt-Csp(2).

Elucidation of Intermediates and Mechanisms in Heterogeneous Catalysis Using Infrared Spectroscopy

Infrared spectroscopy has a long history as a tool for the identification of chemical compounds. More recently, various implementations of infrared spectroscopy have been successfully applied to studies of heterogeneous catalytic reactions with the objective of identifying intermediates and determining catalytic reaction mechanisms. We discuss selective applications of these techniques with a focus on several heterogeneous catalytic reactions, including hydrogenation, deNOx, water-gas shift, and reverse-water-gas shift. The utility of using isotopic substitutions and other techniques in tandem with infrared spectroscopy is discussed. We comment on the modes of implementation and the advantages and disadvantages of the various infrared techniques. We also note future trends and the role of computational calculations in such studies. The infrared techniques considered are transmission Fourier transform infrared spectroscopy, infrared reflection-absorption spectroscopy, polarization-modulation infrared reflection-absorption spectroscopy, sum-frequency generation, diffuse reflectance infrared Fourier transform spectroscopy, attenuated total reflectance, infrared emission spectroscopy, photoacoustic infrared spectroscopy, and surface-enhanced infrared absorption spectroscopy.