Pt8 cluster on alumina under a pressure of hydrogen: Support-dependent reconstruction from first-principles global optimization

Nature of Zirconia on a Copper Inverse Catalyst Under CO2 Hydrogenation Conditions

The growing concern over the escalating levels of anthropogenic CO2 emissions necessitates effective strategies for its conversion to valuable chemicals and fuels. In this research, we embark on a comprehensive investigation of the nature of zirconia on a copper inverse catalyst under the conditions of CO2 hydrogenation to methanol. We employ density functional theory calculations in combination with the Grand Canonical Basin Hopping method, enabling an exploration of the free energy surface including a variable amount of adsorbates within the relevant reaction conditions. Our focus centers on a model three-atom Zr cluster on a Cu(111) surface decorated with various OH, O, and formate ligands, noted Zr3Ox (OH)y (HCOO)z/Cu(111), revealing major changes in the active site induced by various reaction parameters such as the gas pressure, temperature, conversion levels, and CO2/H2 feed ratios. Through our analysis, we have unveiled insights into the dynamic behavior of the catalyst. Specifically, under reaction conditions, we observe a large number of composition and structures with similar free energy for the catalyst, with respect to changing the type, number, and binding sites of adsorbates, suggesting that the active site should be regarded as a statistical ensemble of diverse structures that interconvert.

Amorphous nickel hydroxide shell tailors local chemical environment on platinum surface for alkaline hydrogen evolution reaction

In analogy to natural enzymes, an elaborated design of catalytic systems with a specifically tailored local chemical environment could substantially improve reaction kinetics, effectively combat catalyst poisoning effect and boost catalyst lifetime under unfavourable reaction conditions. Here we report a unique design of 'Ni(OH)2-clothed Pt-tetrapods' with an amorphous Ni(OH)2 shell as a water dissociation catalyst and a proton conductive encapsulation layer to isolate the Pt core from bulk alkaline electrolyte while ensuring efficient proton supply to the active Pt sites. This design creates a favourable local chemical environment to result in acidic-like hydrogen evolution reaction kinetics with a lowest Tafel slope of 27 mV per decade and a record-high specific activity and mass activity in alkaline electrolyte. The proton conductive Ni(OH)2 shell can also effectively reject impurity ions and retard the Oswald ripening, endowing a high tolerance to solution impurities and exceptional long-term durability that is difficult to achieve in the naked Pt catalysts. The markedly improved hydrogen evolution reaction activity and durability in an alkaline medium promise an attractive catalyst material for alkaline water electrolysers and renewable chemical fuel generation.

Elucidation of the Active Site for the Oxygen Evolution Reaction on a Single Pt Atom Supported on Indium Tin Oxide

Single-atom catalysts (SACs) have attracted attention for their high catalytic activity and selectivity, but the nature of their active sites under realistic reaction conditions, involving various ligands, is not well-understood. In this study, we use density functional theory calculations and grand canonical basin hopping to theoretically investigate the active site for the oxygen evolution reaction (OER) on a single Pt atom supported on indium tin oxide, including the influence of the electrochemical potential. We show that the ligands on the Pt atom change from Pt-OH in the absence of electrochemical potential to PtO(OH)₄ in electrochemical conditions. This change of the chemical state of Pt is associated with a decrease of 0.3 V for the OER overpotential. This highlights the importance of accurately identifying the nature of the active site under reaction conditions and the impact of adsorbates on the electrocatalytic activity. This theoretical investigation enhances our understanding of SACs for the OER.

Electrocatalytic Hydrogen Evolution at Full Atomic Utilization over ITO-Supported Sub-nano-Ptn Clusters: High, Size-Dependent Activity Controlled by Fluxional Pt Hydride Species

A combination of density functional theory (DFT) and experiments with atomically size-selected Pt_n clusters deposited on indium-tin oxide (ITO) electrodes was used to examine the effects of applied potential and Pt_n size on the electrocatalytic activity of Pt_n (n = 1, 4, 7, and 8) for the hydrogen evolution reaction (HER). Activity is found to be negligible for isolated Pt atoms on ITO, increasing rapidly with Pt_n size such that Pt₇/ITO and Pt₈/ITO have roughly double the activity per Pt atom compared to atoms in the surface layer of polycrystalline Pt. Both the DFT and experiment find that hydrogen under-potential deposition (H_upd) results in Pt_n/ITO (n = 4, 7, and 8) adsorbing ∼2H atoms/Pt atom at the HER threshold potential, equal to ca. double the H_upd observed for Pt bulk or nanoparticles. The cluster catalysts under electrocatalytic conditions are hence best described as a Pt hydride compound, significantly departing from a metallic Pt cluster. The exception is Pt₁/ITO, where H adsorption at the HER threshold potential is energetically unfavorable. The theory combines global optimization with grand canonical approaches for the influence of potential, uncovering the fact that several metastable structures contribute to the HER, changing with the applied potential. It is hence critical to include reactions of the ensemble of energetically accessible Pt_nH_x/ITO structures to correctly predict the activity vs Pt_n size and applied potential. For the small clusters, spillover of H_ads from the clusters to the ITO support is significant, resulting in a competing channel for loss of H_ads, particularly at slow potential scan rates.

Atomic-Scale Mechanism of Platinum Catalyst Restructuring under a Pressure of Reactant Gas

Heterogeneous catalysis is key for chemical transformations. Understanding how catalysts' active sites dynamically evolve at the atomic scale under reaction conditions is a prerequisite for accurately determining catalytic mechanisms and predictably developing catalysts. We combine in situ time-dependent scanning tunneling microscopy observations and machine-learning-accelerated first-principles atomistic simulations to uncover the mechanism of restructuring of Pt catalysts under a pressure of carbon monoxide (CO). We show that a high CO coverage at a Pt step edge triggers the formation of atomic protrusions of low-coordination Pt atoms, which then detach from the step edge to create sub-nano-islands on the terraces, where under-coordinated sites are stabilized by the CO adsorbates. The fast and accurate machine-learning potential is key to enabling the exploration of tens of thousands of configurations for the CO-covered restructuring catalyst. These studies open an avenue to achieve an atomic-scale understanding of the structural dynamics of more complex metal nanoparticle catalysts under reaction conditions.

Global optimization of extraframework ensembles in zeolites: structural analysis of extraframework aluminum species in MOR and MFI zeolites

Metal-modified zeolites are versatile catalytic materials with a wide range of industrial applications. Their catalytic behaviour is determined by the nature of externally introduced cationic species, i.e., its geometry, chemical composition, and location within the zeolite pores. Superior catalyst designs can be unlocked by understanding the confinement effect and spatial limitations of the zeolite framework and its influence on the geometry and location of such cationic active sites. In this study, we employ the genetic algorithm (GA) global optimization method to investigate extraframework aluminum species and their structural variations in different zeolite matrices. We focus on extraframework aluminum (EFAl) as a model system because it greatly influences the product selectivity and catalytic stability in several zeolite catalyzed processes. Specifically, the GA was used to investigate the configurational possibilities of EFAl within the mordenite (MOR) and ZSM-5 frameworks. The xTB semi-empirical method within the GA was employed for an automated sampling of the EFAl-zeolite space. Furthermore, geometry refinement at the density functional theory (DFT) level of theory allowed us to improve the most stable configurations obtained from the GA and elaborate on the limitations of the xTB method. A subsequent ab initio thermodynamics analysis (aiTA) was chosen to predict the most favourable EFAl structure(s) under the catalytically relevant operando conditions.

Hydrogen-Induced Restructuring of a Cu(100) Electrode in Electroreduction Conditions

The rearrangement of Cu surfaces under electrochemical conditions is known to play a key role in the surface activation for major electrocatalytic reactions. Despite the extensive experimental insights into such rearrangements, from surface-sensitive spectroscopy and microscopy, the spatial and temporal resolution of these methods is insufficient to provide an atomistic picture of the electrochemical interface. Theoretical characterization has also been challenged by the diversity of restructuring configurations, surface stoichiometry, adsorbate configurations, and the effect of the electrode potential. Here, atomistic insight into the restructuring of the electrochemical interface is gained from first principles. Cu(100) restructuring under varying applied potentials and adsorbate coverages is studied by grand canonical density functional theory and global optimization techniques, as well as ab initio molecular dynamics and mechanistic calculations. We show that electroreduction conditions cause the formation of a shifted-row reconstruction on Cu(100), induced by hydrogen adsorption. The reconstruction is initiated at 1/6 ML H coverage, when the Cu-H bonding sufficiently weakens the Cu-Cu bonds between the top- and sublayer, and further stabilized at 1/3 ML when H adsorbates fill all the created 3-fold hollow sites. The simulated scanning tunneling microscopy (STM) images of the calculated reconstructed interfaces agree with experimental in situ STM. However, compared to the thermodynamic prediction, the onsets of reconstruction events in the experiment occur at more negative applied voltages. This is attributed to kinetic effects in restructuring, which we describe via different statistical models, to produce the potential- and pH-dependent surface stability diagram. This manuscript provides rich atomistic insight into surface restructuring in electroreduction conditions, which is required for the understanding and design of Cu-based materials for electrocatalytic processes. It also offers the methodology to study the problem of in situ electrode reconstruction.

Theoretical Perspective on Operando Spectroscopy of Fluxional Nanocatalysts

Improvements in operando spectroscopy have enabled the catalysis community to investigate the dynamic nature of catalysts under operating conditions with increasing detail. Still, the highly dynamic nature of some catalysts, such as fluxional supported subnano clusters, presents a formidable challenge even for the most state-of-the-art techniques. The reason is that such fluxional catalytic interfaces contain a variety of thermally accessible states. Operando spectroscopies used in catalysis generally fall into two categories: ensemble-based techniques, which provide spectra containing the signals of the entire ensemble of states of the catalyst and are not necessarily dominated by the most active species, and localized techniques, which provide atomistic-level information about the dynamics of active sites in a very small area, which might not include the most active species. Combining many different kinds of techniques can provide detailed insight; however, we propose that effective utilization of specific computational techniques and approaches within the fluxionality paradigm can fill the gap and enable atomistic characterization of the most relevant catalytic sites.

Structural dynamics of Ru clusters during nitrogen dissociation in ammonia synthesis

The dynamic evolution of catalyst structures greatly influences the reactivity, especially sub-nanometer clusters, exhibiting complex configurational fluctuation. In the present work, we study the structural dynamics of a Ru₁₉ cluster during the dissociation of N₂ and calculate the reaction free energies using ab initio molecular dynamics (AIMD). Our AIMD calculation predicts a peak-shaped reaction entropy curve due to the adsorption-induced phase transition of the Ru₁₉ cluster. The low melting points of sub-nanometer clusters make it possible to activate N₂ at low temperatures. This work demonstrates that the dynamic changes of cluster structures have a non-negligible effect on reaction free energy and offer an opportunity for achieving ammonia synthesis under mild conditions.

Size-dependent phase transitions boost catalytic activity of sub-nanometer gold clusters

The characterization and identification of the dynamics of cluster catalysis are crucial to unraveling the origin of catalytic activity. However, the dynamical catalytic effects during the reaction process remain unclear. Herein, we investigate the dynamic coupling effect of elementary reactions with the structural fluctuations of sub-nanometer Au clusters with different sizes using ab initio molecular dynamics and the free energy calculation method. It was found that the adsorption-induced solid-to-liquid phase transitions of the cluster catalysts give rise to abnormal entropy increase, facilitating the proceeding of reaction, and this phase transition catalysis exists in a range of clusters with different sizes. Moreover, clusters with different sizes show different transition temperatures, resulting in a non-trivial size effect. These results unveil the dynamic effect of catalysts and help understand cluster catalysis to design better catalysts rationally.

Co and Ni Incorporated γ-Al2O3 (110) Surface: A Density Functional Theory Study

Investigation into the state and mechanisms of the active metal substitution into the γ-Al2O3 support is the basis for design of many catalysts. Periodic density functional theory (DFT) +U calculations were used to investigate the surface properties of transition metals Co3+ and Ni3+ cations substitute for the Al3+ cations of γ-Al2O3 (110) surface. It was found that the substitution energy of one Al3+ substituted by Co3+ and Ni3+ are −61 and −57 kJ/mol, respectively. The Co and Ni preferentially substitute the tetrahedral Al sites instead of the octahedral Al sites. Using thermodynamics, the Al atoms in the top layer of γ-Al2O3 (110) can be 100% substituted by Co and Ni. Ni is easier to substitute the Al atom than Co. There is no obvious structural distortion that occurs after Co and Ni substituted all the top layer Al atoms. While the band gaps of the substituted surface become narrower, resulting in the increase of surface Lewis acidity. In addition, the oxygen vacancy formation energies of the Co and Ni substituted surface are 268 and 53 kJ/mol, respectively. The results provide interface structure and physical chemistry properties of metal-doped catalysts.

Real-time dynamics and structures of supported subnanometer catalysts via multiscale simulations

Understanding the performance of subnanometer catalysts and how catalyst treatment and exposure to spectroscopic probe molecules change the structure requires accurate structure determination under working conditions. Experiments lack simultaneous temporal and spatial resolution and could alter the structure, and similar challenges hinder first-principles calculations from answering these questions. Here, we introduce a multiscale modeling framework to follow the evolution of subnanometer clusters at experimentally relevant time scales. We demonstrate its feasibility on Pd adsorbed on CeO2(111) at various catalyst loadings, temperatures, and exposures to CO. We show that sintering occurs in seconds even at room temperature and is mainly driven by free energy reduction. It leads to a kinetically (far from equilibrium) frozen ensemble of quasi-two-dimensional structures that CO chemisorption and infrared experiments probe. CO adsorption makes structures flatter and smaller. High temperatures drive very rapid sintering toward larger, stable/metastable equilibrium structures, where CO induces secondary structure changes only.

Achieving industrial ammonia synthesis rates at near-ambient conditions through modified scaling relations on a confined dual site

The production of ammonia through the Haber-Bosch process is regarded as one of the most important inventions of the 20th century. Despite significant efforts in optimizing the process, it still consumes 1 to 2% of the worldwide annual energy for the high working temperatures and pressures. The design of a catalyst with a high activity at milder conditions represents another challenge for this reaction. Herein, we combine density functional theory and microkinetic modeling to illustrate a strategy to facilitate low-temperature and -pressure ammonia synthesis through modified energy-scaling relationships using a confined dual site. Our results suggest that an ammonia synthesis rate two to three orders of magnitude higher than the commercial Ru catalyst can be achieved under the same reaction conditions with the introduction of confinement. Such strategies will open pathways for the development of catalysts for the Haber-Bosch process that can operate at milder conditions and present more economically viable alternatives to current industrial solutions.

Reagent-Triggered Isomerization of Fluxional Cluster Catalyst via Dynamic Coupling

Metallic cluster catalysts have many thermodynamically accessible isomers with diverse active sites and low reaction barriers, and lately a strong hypothesis emerged that the many catalyst states collectively drive the catalysis. However, it remained a hypothesis that catalyst isomerization is actually kinetically feasible under the current reaction conditions. Using high-temperature dynamics simulations and sampling, a range of orientations, and vibrational energy distributions, we probe how thermal effects and molecular events affect cluster catalyst dynamics. We show that even such a delicate affair as the dissociation or scattering of a methane molecule on the heavy and thus slow Pt₁₃ cluster triggers substantial isomerization of the catalyst, far beyond thermal at 700 K. A kinetic coupling between the methane activity and cluster catalyst dynamics is observed. In return, the thermal dynamics of the cluster affects the methane reaction and scattering probabilities. Hence, molecular events at the surfaces of fluxional cluster catalysts should facilitate the population of an ensemble of catalyst states under the current reaction conditions, with implications for available active sites, reaction mechanisms, and apparent rates.