Search

To elucidate interfacial dynamics during electrocatalytic reactions, it is crucial to understand the adsorption behavior of organic molecules on catalytic electrodes within the electric double layer (EDL). However, the EDL structure in aqueous environments remains intricate when it comes to the electrochemical amination of acetone, using methylamine as a nitrogen source. Specifically, the interactions of acetone and methylamine with the copper electrode in water remain unclear, posing challenges in the prediction and optimization of reaction outcomes. In this study, initial investigations employed impedance spectroscopy at the potential of zero charge to explore the surface preconfiguration. Here, the capacitance of the EDL was utilized as a primary descriptor to analyze the adsorption tendencies of both acetone and methylamine. Acetone shows an increase in the EDL capacitance, while methylamine shows a decrease. Experiments are interpreted using combined grand canonical density functional theory and ab initio molecular dynamics to delve into the microscopic configurations, focusing on their capacitance and polarizability. Methylamine and acetone have larger molecular polarizability than water. Acetone shows a partial hydrophobic character due to the methyl groups, forming a distinct adlayer at the interface and increasing the polarizability of the liquid interface component. In contrast, methylamine interacts more strongly with water due to its ability to both donate and accept hydrogen bonds, leading to a more significant disruption of the hydrogen bond network. This disruption of the hydrogen network decreases the local polarizability of the interface and decreases the effective capacitance. Our findings underscore the pivotal role of EDL capacitance and polarizability in determining the local reaction environment, shedding light on the fundamental processes important for electro-catalysis.

Cu dissolution has been identified as the dominant process that causes cathode degradation and losses even under cathodic conditions involving methylamine. Despite extensive experimental research, our fundamental and theoretical understanding of the atomic-scale mechanism for Cu dissolution under electrochemical conditions, eventually coupled with surface restructuring processes, is limited. Here, driven by the observation that the working Cu electrode is corroded using mixtures of acetone and methylamine even under reductive potential conditions (-0.75 V vs. RHE), we employed Grand Canonical density functional theory to understand this dynamic process under potential from a microscopic perspective. We show that amine ligands in solution directly chemisorb on the electrode, coordinate with the metal center, and drive the rearrangement of the copper surface by extracting Cu as adatoms in low coordination positions, where other amine ligands can coordinate and stabilize a surface copper-ligand complex, finally forming a detached Cu-amine cationic complex in solution, even under negative potential conditions. Calculations predict that dissolution would occur for a potential of -1.1 V vs. RHE or above. Our work provides a fundamental understanding of Cu dissolution facilitated by surface restructuring in amine solutions under electroreduction conditions, which is required for the rational design of durable Cu-based cathodes for electrochemical amination or other amine involving reduction processes.

In the context of carbon reutilization, the direct electroreduction of captured CO2 (c-CO2RR) appears as an appealing approach since it avoids the energetically costly separation of CO2 from the capture agent. In this process, CO2 is directly reduced from its captured form. Here, we investigate the influence of the capture agent and proton source on that reaction from a combination of theory and experiment. Specifically, we consider methoxide-captured CO2, NH3-captured CO2, and bicarbonate on silver electrocatalysts. We show that the proton source plays a key role in the interplay of the chemistries for the electroreduction of protons, free CO2, and captured CO2. Our density functional theory calculations, including the influence of the potential, demonstrate that a proton source with smaller pKa improves the reactivity for c-CO2RR, but also increases the selectivity toward the hydrogen evolution reaction (HER) on silver surfaces. Since c-CO2RR requires an additional chemical protonation step, the influence of the proton source is stronger than that of the HER. However, c-CO2RR cannot compete with the HER on Ag, Experimentally, the dominant product observed is H2 with low amounts of CO being produced. Through a rotating cylinder electrode cell of well-defined mass-transport properties, we conclude that although methanol solvent exhibits a lower HER activity, HER remains dominant over c-CO2RR. Our work suggests that methoxide is a potential alternative capture agent to NH3 for direct reduction of captured CO2, though challenges in catalyst design, particularly in reducing the onset potential of c-CO2RR to surpass the HER, remain to be addressed.

Heterogeneous nanocatalysts play a crucial role in both the chemical and energy industries. Despite substantial advancements in theoretical, computational and experimental studies, identifying their active sites remains a major challenge. Here we utilize atomic electron tomography to determine the three-dimensional atomic structure of PtNi and Mo-doped PtNi nanocatalysts for the electrochemical oxygen reduction reaction. We then employ the experimental atomic structures as input to first-principles-trained machine learning to identify the active sites of the nanocatalysts. Through the analysis of the structure–activity relationships, we formulate an equation termed the local environment descriptor, which balances the strain and ligand effects to provide physical and chemical insights into active sites in the oxygen reduction reaction. The ability to determine the three-dimensional atomic structure and chemical composition of realistic nanoparticles, combined with machine learning, could transform our fundamental understanding of the active sites of catalysts and guide the rational design of optimal nanocatalysts. (Figure presented.)

The dynamic restructuring of Cu has been observed under electrochemical conditions, and it has been hypothesized to underlie the unique reactivity of Cu toward CO2 electroreduction. Roughening is one of the key surface phenomena for Cu activation, whereby numerous atomic vacancies and adatoms form. However, the atomic structure of such surface motifs in the presence of relevant adsorbates has remained elusive. Here, we explore the chemical space of Cu surface restructuring under coverage of CO and H in realistic electroreduction conditions, by combining grand canonical DFT and global optimization techniques, from which we construct a potential-dependent grand canonical ensemble representation. The regime of intermediate and mixed CO and H coverage─where structures exhibit some elevated surface Cu─is thermodynamically unfavorable yet kinetically inevitable. Therefore, we develop a quasi-kinetic Monte Carlo simulation to track the system's evolution during a simulated cathodic scan. We reveal the evolution path of the system across coverage space and identify the accessible metastable structures formed along the way. Chemical bonding analysis is performed on the metastable structures with elevated Cu*CO species to understand their formation mechanism. By molecular dynamics simulations and free energy calculations, the surface chemistry of the Cu*CO species is explored, and we identify plausible mechanisms via which the Cu*CO species may diffuse or dimerize. This work provides rich atomistic insights into the phenomenon of surface roughening and the structure of involved species. It also features generalizable methods to explore the chemical space of restructuring surfaces with mixed adsorbates and their nonequilibrium evolution.

The decomposition of dimethyl methyl phosphonate (DMMP), a simulant for the nerve agent sarin, was investigated on Cu4/TiO2(110) and K/Cu4/TiO2(110) surfaces using a combination of near-ambient-pressure X-ray photoelectron spectroscopy (NAP-XPS) and density functional theory calculations (DFT). Mass-selected Cu4 clusters and potassium (K) atoms were deposited onto TiO2(110) as a metal catalyst and alkali promoter to improve the reactivity and recyclability of the TiO2 surface after exposure to DMMP. Surface reaction products resulting from decomposition of DMMP were probed by NAP-XPS measurements of phosphorus (P) 2p and carbon 1s core-level spectra. The Cu4/TiO2(110) surface is found to be very active for DMMP decomposition with highly reduced P-species observed even at room temperature (RT). The codeposition of K atoms and Cu4 clusters further improves the reactivity with no intact DMMP detectable. Temperature-dependent measurements show that the presence of K atoms promotes the removal of residual P-species at temperatures > 600 K. Detailed DFT calculations were performed to determine the surface structures and energetically accessible pathways for DMMP decomposition on Cu4/TiO2(110) and K/Cu4/TiO2(110) surfaces. The calculations show that DMMP and P-containing reaction products preferentially bind to the TiO2 surface, while the molecular fragments, i.e., methoxy and methyl, bind to both the Cu4 clusters and TiO2. The Cu4 clusters make the P-O, O-C, and P-C bond cleavages of DMMP markedly more exothermic. The Cu4 clusters are highly fluxional with atomic structures that depend on the configuration of fragments bound to them. Finally, the manifold of P 2p chemical shifts calculated for a large number of energetically favorable configurations of decomposition products is in good agreement with the observed XPS spectra and provides an alternative way of interpreting incompletely resolved core-level spectra using an ensemble of observed structures.

Hydrazine assisted water electrolysis provides an attractive pathway for low-voltage hydrogen production while at the same time mitigating the hazardous hydrazine environmental pollutants. Herein we report the design and synthesis of Ru decorated Ag nanoparticles (NPs) where the Ag-Ru interfaces act as highly effective bifunctional electrocatalysts for the hydrazine oxidation reaction (HzOR) and the hydrogen evolution reaction (HER). The electrocatalysts with Ag-Ru interfaces demonstrate improved HzOR performance with lower overpotential, enhanced mass activity (MA) and highly selective oxidation of hydrazine into N2. Density functional theory (DFT) computations reveal the Ag-Ru interfaces feature higher barrier for N-N bond cleavage and easier N2 desorption, contributing to the electrocatalytic activity and selectivity. At the same time, improved HER performance is also observed due to the more favourable hydrogen desorption. Together, by employing the Ru decorated Ag NPs as electrocatalysts for both HzOR and HER, the hydrazine assisted water electrolyser delivers record-high performance with a current density of 100 mA cm−2 at an ultralow cell voltage of 16 mV and a high current density of 983 ± 30 mA cm−2 at a cell voltage of 0.45 V without any IR compensation.

Addressing sustainable energy storage remains crucial for transitioning to renewable sources. While Li-ion batteries have made significant contributions, enhancing their capacity through alternative materials remains a key challenge. Micro-crystalline silicon is a promising anode material due to its tenfold higher theoretical capacity compared to conventional graphite. However, its substantial volumetric expansion during cycling impedes practical application due to mechanical failure and rapid capacity fading. We propose a novel approach to mitigate this issue by incorporating trace amounts of aluminum into the micro-crystalline silicon electrode using ball milling. We employ density functional theory (DFT) to establish a theoretical framework elucidating how grain boundary sliding, a key mechanism involved in preventing mechanical failure, is facilitated by the presence of trace aluminum at grain boundaries. This, in turn, reduces stress accumulation within the material, reducing the likelihood of failure. To validate our theoretical predictions, we conducted capacity retention experiments on undoped and Al-doped micro-crystalline silicon samples. The results demonstrate significantly reduced capacity fading in the doped sample, corroborating the theoretical framework and showcasing the potential of aluminum doping for improved Li-ion battery performance.

CO2 electroreduction (CO2 R) operating in acidic media circumvents the problems of carbonate formation and CO2 crossover in neutral/alkaline electrolyzers. Alkali cations have been universally recognized as indispensable components for acidic CO2 R, while they cause the inevitable issue of salt precipitation. It is therefore desirable to realize alkali-cation-free CO2 R in pure acid. However, without alkali cations, stabilizing *CO2 intermediates by catalyst itself at the acidic interface poses as a challenge. Herein, we first demonstrate that a carbon nanotube-supported molecularly dispersed cobalt phthalocyanine (CoPc@CNT) catalyst provides the Co single-atom active site with energetically localized d states to strengthen the adsorbate-surface interactions, which stabilizes *CO2 intermediates at the acidic interface (pH=1). As a result, we realize CO2 conversion to CO in pure acid with a faradaic efficiency of 60 % at pH=2 in flow cell. Furthermore, CO2 is successfully converted in cation exchanged membrane-based electrode assembly with a faradaic efficiency of 73 %. For CoPc@CNT, acidic conditions also promote the intrinsic activity of CO2 R compared to alkaline conditions, since the potential-limiting step, *CO2 to *COOH, is pH-dependent. This work provides a new understanding for the stabilization of reaction intermediates and facilitates the designs of catalysts and devices for acidic CO2 R.

The sulfur reduction reaction (SRR) plays a central role in high-capacity lithium sulfur (Li-S) batteries. The SRR involves an intricate, 16-electron conversion process featuring multiple lithium polysulfide intermediates and reaction branches1-3. Establishing the complex reaction network is essential for rational tailoring of the SRR for improved Li-S batteries, but represents a daunting challenge4-6. Herein we systematically investigate the electrocatalytic SRR to decipher its network using the nitrogen, sulfur, dual-doped holey graphene framework as a model electrode to understand the role of electrocatalysts in acceleration of conversion kinetics. Combining cyclic voltammetry, in situ Raman spectroscopy and density functional theory calculations, we identify and directly profile the key intermediates (S8, Li2S8, Li2S6, Li2S4 and Li2S) at varying potentials and elucidate their conversion pathways. Li2S4 and Li2S6 were predominantly observed, in which Li2S4 represents the key electrochemical intermediate dictating the overall SRR kinetics. Li2S6, generated (consumed) through a comproportionation (disproportionation) reaction, does not directly participate in electrochemical reactions but significantly contributes to the polysulfide shuttling process. We found that the nitrogen, sulfur dual-doped holey graphene framework catalyst could help accelerate polysulfide conversion kinetics, leading to faster depletion of soluble lithium polysulfides at higher potential and hence mitigating the polysulfide shuttling effect and boosting output potential. These results highlight the electrocatalytic approach as a promising strategy for tackling the fundamental challenges regarding Li-S batteries.

In acidic conditions, the electroreduction of CO or CO2 (noted CO(2)RR) on metal surfaces is conventionally hindered by intense competition with the hydrogen evolution reaction (HER). In this study, we present first-principles calculations of a mechanism wherein the formation of H-induced Cu adatoms on Cu(111) serves as a pivotal trigger for CORR in acidic environments. Through an analysis of the grand canonical surface state population, we elucidate that these newly formed adatoms create an array of active sites essential for both CO adsorption and subsequent reduction. Our ensemble-based kinetic models unveil the role of adatoms, enhancing the HER while simultaneously initiating CORR. Notably, the cumulative activity of the HER and CORR is contingent upon the combination of various surface states, with their individual contributions varying based on the electrode potential and pH. The interplay between surface state dynamics and electrochemical activity sheds new light on the potential-dependent nature of the active site and reaction kinetics governing CORR on Cu(111) in acidic media.

Selective hydrogenation of α,β-unsaturated aldehydes to produce unsaturated alcohols remains a challenge in catalysis. Here, we explore, on the basis of first-principles simulations, single-atom alloy (SAA) catalysts on copper as a class of catalytic materials to enhance the selectivity for C═O bond hydrogenation in unsaturated aldehydes by controlling the binding strength of the C═C and C═O bonds. We show that on SAA of early transition metals such as Ti, Zr, and Hf, the C═O binding mode of acrolein is favored but the strong binding renders subsequent hydrogenation and desorption impossible. On SAA of late-transition metals, on the other hand, the C═C binding mode is favored and C═C bond hydrogenation follows, resulting in the production of undesired saturated aldehydes. Mid-transition metals (Cr and Mn) in Cu(111) appear as the optimal systems, since they favor acrolein adsorption via the C═O bond but with a moderate binding strength, compatible with catalysis. Additionally, acrolein migration from the C═O to the C═C binding mode, which would open the low energy path for C═C bond hydrogenation, is prevented by a large barrier for this process. SAA of Cr in Cu appears as an optimal candidate, and kinetic simulations show that the selectivity for propenol formation is controlled by preventing the acrolein migration from the more stable C═O to the less stable C═C binding mode and subsequent H-migration and by the formation of the O-H bond from the monohydrogenated intermediate. Dilute alloy catalysts therefore enable tuning the binding strength of intermediates and transition states, opening control of catalytic activity and selectivity.

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.

Noble metals supported on reducible oxides, like CoOx and TiOx, exhibit superior activity in many chemical reactions, but the origin of the increased activity is not well understood. To answer this question we studied thin films of CoOx supported on an Au(111) single crystal surface as a model for the CO oxidation reaction. We show that three reaction regimes exist in response to chemical and topographic restructuring of the CoOx catalyst as a function of reactant gas phase CO/O2 stoichiometry and temperature. Under oxygen-lean conditions and moderate temperatures (≤150 °C), partially oxidized films (CoOx1) forms containing Co3+ species that are effective in a wide temperature range. Resonant photoemission spectroscopy (ResPES) revealed the unique role of Co3+ sites in catalyzing the CO oxidation. Density function theory (DFT) calculations provided deeper insights into the pathway and free energy barriers for the reactions on these oxide phases. These findings in this work highlight the versatility of catalysts and their evolution to form different active phases, both topological and chemically, in response to reaction conditions exposing a new paradigm in the catalyst structure during operation.

Despite the wide applications, the ab initio modeling of the ceria based catalyst is challenging. The partial occupation in the 4f orbitals creates a fundamental challenge for commonly used density functional theory (DFT) methods, including semilocal functionals with Hubbard U correction to force localization and hybrid functionals. In this work, we benchmark the random phase approximation (RPA) for ceria surface properties, including surface energy and hydrogenation energy, compared to the results utilizing the DFT + U approach or hybrid functionals. We show that, for the latter approaches, different surface properties require opposite directions of parameter tuning. This forms a dilemma for the parameter based DFT methods, as the improvement of a certain property by tuning parameters will inevitably lead to the worsening of other properties. Our results suggest that the parameter-free many-body perturbation theory methods exemplified by RPA are a promising strategy to escape the dilemma and provide highly accurate descriptions, which will allow us to better understand the catalytic reactions in ceria related systems.

Hydrazine-assisted water electrolysis offers a feasible path for low-voltage green hydrogen production. Herein, the design and synthesis of ultrathin RhRu0.5 -alloy wavy nanowires as bifunctional electrocatalysts for both the anodic hydrazine oxidation reaction (HzOR) and the cathodic hydrogen evolution reaction (HER) is reported. It is shown that the RhRu0.5 -alloy wavy nanowires can achieve complete electrooxidation of hydrazine with a low overpotential and high mass activity, as well as improved performance for the HER. The resulting RhRu0.5 bifunctional electrocatalysts enable, high performance hydrazine-assisted water electrolysis delivering a current density of 100 mA cm-2 at an ultralow cell voltage of 54 mV and a high current density of 853 mA cm-2 at a cell voltage of 0.6 V. The RhRu0.5  electrocatalysts further demonstrate a stable operation at a high current density of 100 mA cm-2 for 80 hours of testing period with little irreversible degradation. The overall performance greatly exceeds that of the previously reported hydrazine-assisted water electrolyzers, offering a pathway for efficiently converting hazardous hydrazine into molecular hydrogen.

As a showcase for experimentally accessible single-site catalysts, well-defined and characterized isolated Ir sites anchored on a covalent triazine framework (CTF) were investigated with computational methods.

We report the size-dependent activity and stability of supported Pt1,4,7,8 for electrocatalytic hydrogen evolution reaction, and show that clusters outperform polycrystalline Pt in activity, with size-dependent stability. To understand the size effects, we use DFT calculations to study the structural fluxionality under varying potentials. We show that the clusters can reshape under H coverage and populate an ensemble of states with diverse stoichiometry, structure, and thus reactivity. Both experiment and theory suggest that electrocatalytic species are hydridic states of the clusters (≈2 H/Pt). An ensemble-based kinetic model reproduces the experimental activity trend and reveals the role of metastable states. The stability trend is rationalized by chemical bonding analysis. Our joint study demonstrates the potential- and adsorbate-coverage-dependent fluxionality of subnano clusters of different sizes and offers a systematic modeling strategy to tackle the complexities.

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)4 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.

A combination of density functional theory (DFT) and experiments with atomically size-selected Ptn clusters deposited on indium-tin oxide (ITO) electrodes was used to examine the effects of applied potential and Ptn size on the electrocatalytic activity of Ptn (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 Ptn size such that Pt7/ITO and Pt8/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 (Hupd) results in Ptn/ITO (n = 4, 7, and 8) adsorbing ∼2H atoms/Pt atom at the HER threshold potential, equal to ca. double the Hupd 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 Pt1/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 PtnHx/ITO structures to correctly predict the activity vs Ptn size and applied potential. For the small clusters, spillover of Hads from the clusters to the ITO support is significant, resulting in a competing channel for loss of Hads, particularly at slow potential scan rates.

We report the size-dependent activity and stability of FTO-supported Pt1,4,7,8 for electrocatalytic hydrogen evolution reaction, and show that clusters are significantly more active than polycrystalline Pt, but also have stability under HER conditions that is size-dependent. To understand the mechanistic origin of the size effects, we carried out detailed DFT-based theory, accounting for the structural fluxionality under varying potentials. We show that, even under coverage of the simplest adsorbate H, the clusters can undergo drastic changes and populate a grand canonical ensemble of hydride states with diverse distributions of stoichiometry, structure, and thus reactivity. Both experiment and theory find that electrocatalysis mainly results from heavily hydrogenated states of the clusters (~2 H/Pt), and an ensemble-based kinetic model reproduces the experimental activity trend and provides detailed insights into the contributions from metastable configurations. The size-dependent stability trend is rationalized by chemical bonding analysis. The comparison of experiment and theory demonstrates the potential- and adsorbate-coverage-dependent fluxionality of subnano clusters of different sizes and offers a systematic modeling strategy to tackle the complexities.

A combination of density functional theory (DFT) and experiments with atomically size-selected Ptn clusters deposited on indium-tin oxide (ITO) electrodes was used to examine the effects of applied potential and Ptn size on the electrocatalytic activity of Ptn (n = 1, 4, 7, 8) for the hydrogen evolution reaction (HER). Activity is found to be negligible for isolated Pt atoms on ITO, increasing rapidly with Ptn size, such that Pt7/ITO and Pt8/ITO have roughly double the activity per Pt atom compared atoms in the surface layer of polycrystalline Pt. Both DFT and experiment find that hydrogen under-potential deposition (Hupd) results in Ptn/ITO (n = 4, 7, 8) adsorbing ~2 H atoms/Pt atom at the HER threshold potential, equal to ca. double the Hupd observed for bulk Pt or Pt 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 Pt1/ITO, where H adsorption at the HER threshold potential is energetically unfavorable. Theory combines global optimization with grand canonical approaches for the influence of potential, uncovering that several metastable structures contribute to HER, changing with the applied potential. It is hence critical to include reactions of the ensemble of energetically accessible PtnHx/ITO structures to correctly predict the activity vs. Ptn size and applied potential. For the small clusters, spillover of Hads from the clusters to the ITO support is significant, resulting in a competing channel for loss of Hads, particularly at slow potential scan rates.

The oxygen reduction reaction (ORR) plays a key role in renewable energy transformation processes. Unfortunately, it is inherently sluggish, which greatly limits its industrial application. Sub-nano-cluster-decorated electrode interfaces are promising candidate ORR electrocatalysts. However, understanding the nature of the active sites on these catalysts under electrocatalytic conditions presents a formidable challenge for both experiment and theory, due to their dynamic fluxional character. Here, we combine global optimization with the electronic Grand Canonical DFT to elucidate the structure and dynamics of subnano Ptnclusters deposited on electrified graphite. We show that, under electrochemical conditions, these clusters exist as statistical ensembles of multiple states, whose fluxionality is greatly affected by the applied potential, electrolyte, and adsorbate coverage. The results reveal the presence of potential-dependent active sites and, hence, reaction energetics.

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.

Copper (Cu) remains the most important metal catalyst for the carbon dioxide reduction reaction (CO2 RR) into C2 products. Due to limited evidence from in situ experiments, mechanistic studies are often performed in the framework of density functional theory (DFT), using functionals at the generalized gradient approximation (GGA) level, which have fundamental difficulties to correctly describe CO adsorption and surface stability. We employ the adiabatic connection fluctuation dissipation theorem within the random phase approximation (RPA), in combination with the linearized Poisson-Boltzmann equation to describe solvation effects, to investigate the mechanism of CO2 RR on the Cu(100) facet. Qualitatively different from the DFT-GGA results, RPA results propose the formation of *OCCHO as the potential determining step towards C2 products. The results suggest that it is important to use more accurate methods like RPA when modeling reactions involving multiple CO-related species like CO2 RR.

Dynamic restructuring of bimetallic catalysts plays a crucial role in their catalytic activity and selectivity. In particular, catalyst pretreatment with species such as carbon monoxide and oxygen has been shown to be an effective strategy for tuning the surface composition and morphology. Mechanistic and kinetic understanding of such restructuring is fundamental to the chemistry and engineering of surface active sites but has remained challenging due to the large structural, chemical, and temporal degrees of freedom. Here, we combine time-resolved temperature-programmed infrared reflection absorption spectroscopy, ab initio thermodynamics, and machine-learning molecular dynamics to uncover previously unidentified timescale and kinetic parameters of in situ restructuring in Pd/Au(111), a highly relevant model system for dilute Pd-in-Au nanoparticle catalysts. The key innovation lies in utilizing CO not only as a chemically sensitive probe of surface Pd but also as an agent that induces restructuring of the surface. Upon annealing in vacuum, as-deposited Pd islands became encapsulated by Au and partially dissolved into the subsurface, leaving behind isolated Pd monomers on the surface. Subsequent exposure to 0.1 mbar CO enabled Pd monomers to repopulate the surface up to 373 K, above which complete Pd dissolution occurred by 473 K, with apparent activation energies of 0.14 and 0.48 eV, respectively. These restructuring processes occurred over the span of ∼1000 s at a given temperature. Such a minute-timescale dynamics not only elucidates the fluxional nature of alloy catalysts but also presents an opportunity to fine-tune the surface under moderate temperature and pressure conditions.

Electrocatalysis plays a key role in sustainable energy conversion and storage. It is critical to model the grand canonical treatment of electrons, which accounts for the electrochemical potential explicitly, at the atomic scale and understand these reactions at electrified interfaces. However, such a grand canonical treatment for electrocatalytic modeling is in practice restricted to a treatment of electronic structure with density functional theory, and more accurate methods are in many cases desirable. Here, we develop an original workflow combining the grand canonical treatment of electrons with many-body perturbation theory in its random phase approximation (RPA). Using the potential dependent adsorption of carbon monoxide on the copper (100) facet, we show that the grand canonical RPA energetics provide the correct on-top Cu geometry for CO at reducing potential. The match with experimental results is significantly improved compared to the functionals at the generalized gradient approximation level, which is the most commonly used approximation for electrochemical applications. We expect this development to pave the way to further electrochemical applications using RPA.

The development of new catalyst materials for energy-efficient chemical synthesis is critical as over 80% of industrial processes rely on catalysts, with many of the most energy-intensive processes specifically using heterogeneous catalysis. Catalytic performance is a complex interplay of phenomena involving temperature, pressure, gas composition, surface composition, and structure over multiple length and time scales. In response to this complexity, the integrated approach to heterogeneous dilute alloy catalysis reviewed here brings together materials synthesis, mechanistic surface chemistry, reaction kinetics, in situ and operando characterization, and theoretical calculations in a coordinated effort to develop design principles to predict and improve catalytic selectivity. Dilute alloy catalysts─in which isolated atoms or small ensembles of the minority metal on the host metal lead to enhanced reactivity while retaining selectivity─are particularly promising as selective catalysts. Several dilute alloy materials using Au, Ag, and Cu as the majority host element, including more recently introduced support-free nanoporous metals and oxide-supported nanoparticle "raspberry colloid templated (RCT)" materials, are reviewed for selective oxidation and hydrogenation reactions. Progress in understanding how such dilute alloy catalysts can be used to enhance selectivity of key synthetic reactions is reviewed, including quantitative scaling from model studies to catalytic conditions. The dynamic evolution of catalyst structure and composition studied in surface science and catalytic conditions and their relationship to catalytic function are also discussed, followed by advanced characterization and theoretical modeling that have been developed to determine the distribution of minority metal atoms at or near the surface. The integrated approach demonstrates the success of bridging the divide between fundamental knowledge and design of catalytic processes in complex catalytic systems, which can accelerate the development of new and efficient catalytic processes.

Graphene has received tremendous interest in both chemical and physical fields. Among different edges of the graphene system, the zigzag edge terminated graphene nanoribbons (ZGNRs) show unique magnetic properties in the antiferromagnetic (AFM) state. However, to date, the understanding of ZGNR chemical properties is mainly based on the partial radical concept, and in previous studies, the energy differences between the ferromagnetic (FM) and AFM states are smaller than experimental evidence. Here, we report that the strongly constrained and appropriately normed functional gives a significantly larger energy difference, which matches the experimental observation. Furthermore, utilizing the energetics in the large difference case, we propose a conceptual supplement to the previous partial radical concept: the overall stabilization of the AFM state compared to the nonmagnetic (NM) state consists of two parts that affect the adsorption energy conversely. The NM-FM energy differences will strengthen the adsorption, being in line with the previous partial radical concept. The FM-AFM energy differences will instead weaken the adsorption. We perform calculations of H, OH, and LiS radical adsorption energies on ZGNRs to show that this weakening effect is numerically non-negligible: at least a ∼0.2 eV difference in the adsorption energies is found. We expect that this refinement of the partial radical concept can provide a more comprehensive understanding of the chemical properties of ZGNRs. The differences in adsorption energies for the H, OH, and LiS radicals found here lead to significant changes in the predicted reactivity of the ZGNR models.

With the emergence of methods for computing rate constants for elementary reaction steps of catalytic reactions, benchmarking their accuracy becomes important. The unimolecular dehydrogenation of adsorbed formate on metal surfaces serves as a prototype for comparing experiment and theory. Previously measured pre-exponential factors for CO2 formation from formate on metal surfaces, including Cu(110), are substantially higher than expected from the often used value of k B T/h, or ∼6 × 1012 s-1, suggesting that the entropy of the transition state is higher than that of the adsorbed formate. Herein, the rate constant parameters for formate decomposition on Au(110) and Cu(110) are addressed quantitatively by both experiment and theory and compared. A pre-exponential factor of 2.3 × 1014 s-1 was obtained experimentally on Au(110). DFT calculations revealed the most stable configuration of formate on both surfaces to be bidentate and the transition states to be less rigidly bound to the surface compared to the reactant state, resulting in a higher entropy of activation and a pre-exponential factor exceeding k B T/h. Though reasonable agreement is obtained between experiment and theory for the pre-exponential factors, the activation energies determined experimentally remain consistently higher than those computed by DFT using the GGA-PBE functional. This difference was largely erased when the metaGGA-SCAN functional was applied. This study provides insight into the underlying factors that result in the relatively high pre-exponential factors for unimolecular decomposition on metal surfaces generally, highlights the importance of mobility for the transition state, and offers vital information related to the direct use of DFT to predict rate constants for elementary reaction steps on metal surfaces.