Electroreduction of Captured CO2 on Silver Catalysts: Influence of the Capture Agent and Proton Source

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.

H and CO Co-Induced Roughening of Cu Surface in CO2 Electroreduction Conditions

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.

Stable Pentagonal Layered Palladium Diselenide Enables Rapid Electrosynthesis of Hydrogen Peroxide

Electrosynthesis of hydrogen peroxide (H2O2) via the two-electron oxygen reduction reaction (2e- ORR) is promising for various practical applications, such as wastewater treatment. However, few electrocatalysts are active and selective for 2e- ORR yet are also resistant to catalyst leaching under realistic operating conditions. Here, a joint experimental and computational study reveals active and stable 2e- ORR catalysis in neutral media over layered PdSe2 with a unique pentagonal puckered ring structure type. Computations predict active and selective 2e- ORR on the basal plane and edge of PdSe2, but with distinct kinetic behaviors. Electrochemical measurements of hydrothermally synthesized PdSe2 nanoplates show a higher 2e- ORR activity than other Pd-Se compounds (Pd4Se and Pd17Se15). PdSe2 on a gas diffusion electrode can rapidly accumulate H2O2 in buffered neutral solution under a high current density. The electrochemical stability of PdSe2 is further confirmed by long device operational stability, elemental analysis of the catalyst and electrolyte, and synchrotron X-ray absorption spectroscopy. This work establishes a new efficient and stable 2e- ORR catalyst at practical current densities and opens catalyst designs utilizing the unique layered pentagonal structure motif.

The Impact of Electron Donating and Withdrawing Groups on Electrochemical Hydrogenolysis and Hydrogenation of Carbonyl Compounds

The hydrogenolysis or hydrodeoxygenation of a carbonyl group, where the C═O group is converted to a CH2 group, is of significant interest in a variety of fields. A challenge in electrochemically achieving hydrogenolysis of a carbonyl group with high selectivity is that electrochemical hydrogenation of a carbonyl group, which converts the C═O group to an alcohol group (CH-OH), is demonstrated not to be the initial step of hydrogenolysis. Instead, hydrogenation and hydrogenolysis occur in parallel, and they are competing reactions. This means that although both hydrogenolysis and hydrogenation require adding H atoms to the carbonyl group, they involve different intermediates formed on the electrode surface. Thus, revealing the difference in intermediates, transition states, and kinetic barriers for hydrogenolysis and hydrogenation pathways is the key to understanding and controlling hydrogenolysis/hydrogenation selectivity of carbonyl compounds. In this study, we aimed to identify features of reactant molecules that can affect their hydrogenolysis/hydrogenation selectivity on a Zn electrode that was previously shown to promote hydrogenolysis over hydrogenation. In particular, we examined the electrochemical reduction of para-substituted benzaldehyde compounds with substituent groups having different electron donating/withdrawing abilities. Our results show a strikingly systematic impact of the substituent group where a stronger electron-donating group promotes hydrogenolysis and a stronger electron-withdrawing group promotes hydrogenation. These experimental results are presented with computational results explaining the substituent effects on the thermodynamics and kinetics of electrochemical hydrogenolysis and hydrogenation pathways, which also provide critically needed information and insights into the transition states involved with these pathways.

Composite interfaces of g-C3N4 fragments loaded on a Cu substrate for CO2 reduction

Designing an electrocatalyst with high efficiency and product selectivity is always crucial for an electrocatalytic CO2 reduction reaction (CO2RR). Inspired by the great progress of two-dimensional (2D) nanomaterials growing on Cu surfaces and their promising CO2RR catalytic efficiencies at their interfaces, the unique performance of Cu-based 2D materials as high-efficiency and low-cost CO2RR electrocatalysts has attracted extensive attention. Herein, based on density functional theory (DFT) calculations, we proposed a composite structure of graphitic carbon nitride (g-C3N4) fragments loaded on a Cu surface to explore the CO2RR catalytic property of the interface between g-C3N4 and the Cu surface. Three composite interfaces of C3N4/Cu(111), C3N4/Cu(110) and C3N4/Cu(100) have been studied by considering the reaction sites of vertex nitrogen atoms, edge nitrogen atoms and the nearby Cu atoms. It was found that the C3N4/Cu interfaces where nitrogen atoms contact the Cu substrate present competitive CO2RR activity. Among them, C3N4/Cu(111)-N3 exhibited a better activity for CH3OH production, with a low overpotential of 0.38 V. For HCOOH and CH4 production, C3N4/Cu(111)-Cu and C3N4/Cu(100)-N1 have overpotentials of 0.26 V and 0.44 V. The electronic analysis indicates the electron transfer from the Cu substrate to the g-C3N4 fragment and mainly accumulates on the nitrogen atoms of the interface. Such charge accumulation can activate the adsorbed CO bond of CO2 and lead to lower energetic barriers of CO2RR. DFT calculations indicate that the boundary nitrogen sites reduced the energy barrier of *CHO, which is crucial for CO2RR, compared with that of the pristine Cu surface. Our study explores a new Cu-based electrocatalyst and indicates that the C3N4/Cu interface can enhance the activities and selectivity of CO2RR and open a new strategy to design high-efficiency electrocatalysts for CO2RR.

Understanding the Adiabatic Evolution of Surface States in Tetradymite Topological Insulators under Electrochemical Conditions

Nontrivial surface states in topological materials have emerged as exciting targets for surface chemistry research. In particular, topological insulators have been used as electrodes in electrocatalytic reactions. Herein, we investigate the robustness of the topological surface states and band topology under electrochemical conditions, specifically in the presence of an electric double layer. First-principles band structure calculations are performed on the electrified (111) surfaces of Bi2Te3, Bi2Se3, and Sb2Te3 using an implicit electrolyte model. Our results demonstrate the adiabatic evolution of the surface states upon surface charging. Under oxidizing potentials, the surface states are shifted upward in energy, preserving the Dirac point on the surface and the band inversion in the bulk. Conversely, under reduced potentials, hybridization is observed between the surface and bulk states, suggesting a likely breakdown of topological protection. The position of the Fermi level, which dictates the working states in catalytic reactions, should ideally be confined within the bulk bandgap. This requirement defines a potential window for the effective application of topological electrocatalysis.

A Readily Achieved Potentiostatic Method in Density Functional Theory Calculation for Improved Prediction of the Performance for Electrocatalytic Nitrogen Reduction Reaction

Accurate prediction of the catalytic performance of nitrogen reduction reaction catalysts based on density functional theory (DFT) calculation is of great significance for developing catalytic materials for nitrogen fixation. However, the applied electrode potential induced the fixation of Fermi level and solvation effect are commonly ignored in the current computational hydrogen electrode method, which leads to the large deviation between the calculation predicted limit potential and the experimentally measured limit potential. In this work, the simple external iteration method is proposed to simulate the Fermi level of the catalysts that are fixed by the applied electrode potential, along with the hybrid solvent model to describe the strong interaction, such as hydrogen bond, between the solvent molecules and the intermediates. This method allowed the theoretical and experimental limit potentials to be in good agreement, indicating the significant effect of the electrode potential and solvation in the DFT calculation. These results will guide the calculation-based prediction of other reaction systems in the field of electrocatalysis.

Halide Adsorption Enhances Electrochemical Hydrogenolysis of 5-Hydroxymethylfurfural by Suppressing Hydrogenation

Reductive upgrading of 5-hydroxymethylfurfural (HMF), a biomass-derived platform molecule, to 2,5-dimethylfuran (DMF), a biofuel with an energy density 40% greater than that of ethanol, involves hydrogenolysis of both the aldehyde (C═O) and the alcohol (C-OH) groups of HMF. It is known that when hydrogenation of the aldehyde occurs to form 2,5-bis(hydroxymethyl)furan (BHMF), BHMF cannot be further reduced to DMF. Thus, aldehyde hydrogenation must be suppressed to increase the selectivity for DMF production. Previously, it was shown that on a Cu electrode hydrogenolysis occurs mainly through proton-coupled electron transfer (PCET), where a proton from the solution and an electron from the electrode are transferred to the organic species. In contrast, hydrogenation occurs not only through PCET but also through hydrogen atom transfer (HAT), where a surface-adsorbed hydrogen atom (H*) is transferred to the organic species. This study shows that halide adsorption on Cu can effectively suppress HAT by decreasing the steady-state H* coverage on Cu during HMF reduction. As HAT enables only aldehyde hydrogenation, a striking suppression of BHMF is observed, thereby enhancing DMF production. We discuss how the identity and concentration of the halide, along with the reduction conditions (i.e., potential and pH), affect halide adsorption on Cu and identify when optimal halide coverages are achieved to maximize DMF selectivity. Our experimental results are presented alongside computational results that elucidate how halide adsorption affects the adsorption energy of hydrogen and the steady-state H* coverage on Cu, which provide an atomic-level understanding of all experimentally observed effects.

Adaptive Electric Fields Embedded Electrochemical Barrier Calculations

Electrochemical interfaces are grand canonical ensembles of varying electrons. Simulating them by standard first-principles methods is a challenging task, since the number of electrons (or charge) is fixed in the calculation. Under the constant charge framework, we developed a constant potential simulation method realized by adding an adaptive electric field to a charge neutral system. Electric field is the controlling variable. In addition, we defined an internal reversible hydrogen electrode potential (ϕ_IRHE), which can ensure the model independence of our method. To validate our method, the reaction energies of some electrochemical reactions are calculated, the results are comparable with the computational hydrogen electrode model and experiments. At last, the evolution of transition state structures and charge transfer coefficients of some electrochemical reactions on Ag(111) surface were discussed by our method.

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.

Modeling Electrochemical Processes with Grand Canonical Treatment of Many-Body Perturbation Theory

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.

Emerging Electrochemical Processes to Decarbonize the Chemical Industry

Electrification is a potential approach to decarbonizing the chemical industry. Electrochemical processes, when they are powered by renewable electricity, have lower carbon footprints in comparison to conventional thermochemical routes. In this Perspective, we discuss the potential electrochemical routes for chemical production and provide our views on how electrochemical processes can be matured in academic research laboratories for future industrial applications. We first analyze the CO₂ emission in the manufacturing industry and conduct a survey of state of the art electrosynthesis methods in the three most emission-intensive areas: petrochemical production, nitrogen compound production, and metal smelting. Then, we identify the technical bottlenecks in electrifying chemical productions from both chemistry and engineering perspectives and propose potential strategies to tackle these issues. Finally, we provide our views on how electrochemical manufacturing can reduce carbon emissions in the chemical industry with the hope to inspire more research efforts in electrifying chemical manufacturing.

Theory-guided experimental design in battery materials research

A reliable energy storage ecosystem is imperative for a renewable energy future, and continued research is needed to develop promising rechargeable battery chemistries. To this end, better theoretical and experimental understanding of electrochemical mechanisms and structure-property relationships will allow us to accelerate the development of safer batteries with higher energy densities and longer lifetimes. This Review discusses the interplay between theory and experiment in battery materials research, enabling us to not only uncover hitherto unknown mechanisms but also rationally design more promising electrode and electrolyte materials. We examine specific case studies of theory-guided experimental design in lithium-ion, lithium-metal, sodium-metal, and all-solid-state batteries. We also offer insights into how this framework can be extended to multivalent batteries. To close the loop, we outline recent efforts in coupling machine learning with high-throughput computations and experiments. Last, recommendations for effective collaboration between theorists and experimentalists are provided.

Modeling Operando Electrochemical CO2 Reduction

Since the seminal works on the application of density functional theory and the computational hydrogen electrode to electrochemical CO₂ reduction (eCO₂R) and hydrogen evolution (HER), the modeling of both reactions has quickly evolved for the last two decades. Formulation of thermodynamic and kinetic linear scaling relationships for key intermediates on crystalline materials have led to the definition of activity volcano plots, overpotential diagrams, and full exploitation of these theoretical outcomes at laboratory scale. However, recent studies hint at the role of morphological changes and short-lived intermediates in ruling the catalytic performance under operating conditions, further raising the bar for the modeling of electrocatalytic systems. Here, we highlight some novel methodological approaches employed to address eCO₂R and HER reactions. Moving from the atomic scale to the bulk electrolyte, we first show how ab initio and machine learning methodologies can partially reproduce surface reconstruction under operation, thus identifying active sites and reaction mechanisms if coupled with microkinetic modeling. Later, we introduce the potential of density functional theory and machine learning to interpret data from Operando spectroelectrochemical techniques, such as Raman spectroscopy and extended X-ray absorption fine structure characterization. Next, we review the role of electrolyte and mass transport effects. Finally, we suggest further challenges for computational modeling in the near future as well as our perspective on the directions to follow.

Challenges of modeling nanostructured materials for photocatalytic water splitting

Understanding the water splitting mechanism in photocatalysis is a rewarding goal as it will allow producing clean fuel for a sustainable life in the future. However, identifying the photocatalytic mechanisms by modeling photoactive nanoparticles requires sophisticated computational techniques based on multiscale modeling. In this review, we will survey the strengths and drawbacks of currently available theoretical methods at different length and accuracy scales. Understanding the surface-active site through Density Functional Theory (DFT) using new, more accurate exchange-correlation functionals plays a key role for surface engineering. Larger scale dynamics of the catalyst/electrolyte interface can be treated with Molecular Dynamics albeit there is a need for more generalizations of force fields. Monte Carlo and Continuum Modeling techniques are so far not the prominent path for modeling water splitting but interest is growing due to the lower computational cost and the feasibility to compare the modeling outcome directly to experimental data. The future challenges in modeling complex nano-photocatalysts involve combining different methods in a hierarchical way so that resources are spent wisely at each length scale, as well as accounting for excited states chemistry that is important for photocatalysis, a path that will bring devices closer to the theoretical limit of photocatalytic efficiency.

The Sensitivity of Metal Oxide Electrocatalysis to Bulk Hydrogen Intercalation: Hydrogen Evolution on Tungsten Oxide

Metal oxides are attracting increased attention as electrocatalysts owing to their affordability, tunability, and reactivity. However, these materials can undergo significant chemical changes under reaction conditions, presenting challenges for characterization and optimization. Herein, we combine experimental and computational methods to demonstrate that bulk hydrogen intercalation governs the activity of tungsten trioxide (WO₃) toward the hydrogen evolution reaction (HER). In contrast to the focus on surface processes in heterogeneous catalysis, we demonstrate that bulk oxide modification is responsible for experimental HER activity. Density functional theory (DFT) calculations reveal that intercalation enables the HER by altering the acid-base character of surface sites and preventing site blocking by hydration. First-principles microkinetic modeling supports that the experimental HER rates can only be explained by intercalated H_xWO₃, whereas nonintercalated WO₃ does not catalyze the HER. Overall, this work underscores the critical influence of hydrogen intercalation on aqueous cathodic electrocatalysis at metal oxides.

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. A vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pK_a values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential-dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

Thermodynamic cyclic voltammograms: peak positions and shapes

Based on a mean-field description of thermodynamic cyclic voltammograms (CVs), we analyze here in full generality, how CV peak positions and shapes are related to the underlying interface energetics, in particular when also including electrostatic double layer (DL) effects. We show in particular, how non-Nernstian behaviour is related to capacitive DL charging, and how this relates to common adsorbate-centered interpretations such as a changed adsorption energetics due to dipole-field interactions and the electrosorption valency - the number of exchanged electrons upon electrosorption per adsorbate. Using Ag(111) in halide-containing solutions as test case, we demonstrate that DL effects can introduce peak shifts that are already explained by rationalizing the interaction of isolated adsorbates with the interfacial fields, while alterations of the peak shape are mainly driven by the coverage-dependence of the adsorbate dipoles. In addition, we analyze in detail how changing the experimental conditions such as the ion concentrations in the solvent but also of the background electrolyte can affect the CV peaks via their impact on the potential drop in the DL and the DL capacitance, respectively. These results suggest new routes to analyze experimental CVs and use of those for a detailed assessment of the accuracy of atomistic models of electrified interfaces e.g. with and without explicitly treated interfacial solvent and/or approximate implicit solvent models.

Thermodynamic Cyclic Voltammograms Based on Ab Initio Calculations: Ag(111) in Halide-Containing Solutions

Cyclic voltammograms (CVs) are a central experimental tool for assessing the structure and activity of electrochemical interfaces. Based on a mean-field ansatz for the interface energetics under applied potential conditions, we here derive an ab initio thermodynamics approach to efficiently simulate thermodynamic CVs. All unknown parameters are determined from density functional theory (DFT) calculations coupled to an implicit solvent model. For the showcased CVs of Ag(111) electrodes in halide-anion-containing solutions, these simulations demonstrate the relevance of double-layer contributions to explain experimentally observed differences in peak shapes over the halide series. Only the appropriate account of interfacial charging allows us to capture the differences in equilibrium coverage and total electronic surface charge that cause the varying peak shapes. As a case in point, this analysis demonstrates that prominent features in CVs do not only derive from changes in adsorbate structure or coverage but can also be related to variations of the electrosorption valency. Such double-layer effects are proportional to adsorbate-induced changes in the work function and/or interfacial capacitance. They are thus especially pronounced for electronegative halides and other adsorbates that affect these interface properties. In addition, the analysis allows us to draw conclusions on how the possible inaccuracy of implicit solvation models can indirectly affect the accuracy of other predicted quantities such as CVs.

Challenges for density functional theory: calculation of CO adsorption on electrocatalytically relevant metals

We assess four DFT functionals, RTPSS, RPBE, SCAN and B97M-rV, for surface interactions. We find that B97M-rV predicts the correct site preference for CO binding on Ag and Au while RTPSS performs well for surface relaxations and binding of CO on Pt and Cu.

Defect-Enhanced CO2 Reduction Catalytic Performance in O-Terminated MXenes

Electrochemical carbon dioxide reduction reaction (CO₂ RR) represents a promising way to generate fuels and chemical feedstock sustainably. Recently, studies have shown that two-dimensional metal carbides and nitrides (MXenes) can be promising CO₂ RR electrocatalysts due to the alternating -C and -H coordination with intermediates that decouples scaling relations seen on transition metal catalysts. However, further by tuning the electronic and surface structure of MXenes it should still be possible to reach higher turnover number and selectivities. To this end, defect engineering of MXenes for electrochemical CO₂ RR has not been investigated to date. In this work, first-principles modelling simulations are employed to systematically investigate CO₂ RR on M₂ XO₂ -type MXenes with transition metal and carbon/nitrogen vacancies. We found that the -C-coordinated intermediates take the form of fragments (e. g., *COOH, *CHO) whereas the -H-coordinated intermediates form a complete molecule (e. g., *HCOOH, *H₂ CO). Interestingly, the fragment-type intermediates become more strongly bound when transition-metal vacancies are present on most MXenes, while the molecule-type intermediates are largely unaffected, allowing the CO₂ RR overpotential to be tuned. The most promising defective MXene is Hf₂ NO₂ containing Hf vacancies, with a low overpotential of 0.45 V. More importantly, through electronic structure analysis it could be observed that the Fermi level of the MXene changes significantly in the presence of vacancies, indicating that the Fermi level shift can be used as an ideal descriptor to rapidly predict the catalytic performance of defective MXenes. Such an evaluation strategy is applicable to other catalysts beyond MXenes, which could enhance high throughput screening efforts for accelerated catalyst discovery.

Aromatic Ester-Functionalized Ionic Liquid for Highly Efficient CO2 Electrochemical Reduction to Oxalic Acid

Electrochemical reduction of CO₂ into valuable chemicals is a significant route to utilize CO₂ resources. Among various electroreduction products, oxalic acid (H₂ C₂ O₄ ) is an important chemical for pharmaceuticals, rare earth extraction, and metal processing. Here, an aprotic aromatic ester-functionalized ionic liquid (IL), 4-(methoxycarbonyl) phenol tetraethylammonium ([TEA][4-MF-PhO]), was designed as an electrolyte for CO₂ electroreduction into oxalic acid. It exhibited a large oxalic acid partial current density of 9.03 mA cm^-2 with a faradaic efficiency (FE) of 86 % at -2.6 V (vs. Ag/Ag⁺ ), and the oxalic acid formation rate was as high as 168.4 μmol cm^-2  h^-1 , which is the highest reported value to date. Moreover, the results of density functional theory calculations demonstrated that CO₂ was efficiently activated to a -COOH intermediate by bis-active sites of the aromatic ester anion via the formation of a [4-MF-PhO-COOH]^- adduct, which finally dimerized into oxalic acid.

Demystifying the Atomistic Origin of the Electric Field Effect on Methane Oxidation

Understanding the role of an electric field on the surface of a catalyst is crucial in tuning and promoting the catalytic activity of metals. Herein, we evaluate the oxidation of methane over a Pt surface with varying oxygen coverage using density functional theory. The latter is controlled by the electrode polarization, giving rise to the non-Faradaic modification of catalytic activity phenomenon. At -1 V, the Pt(111) surface is present, while at 1 V, α-PtO₂ on Pt(111) takes over. Our results demonstrate that the alteration of the platinum oxide surface under the influence of an electrochemical potential promotes the catalytic activity of the metal oxides by lowering the activation energy barrier of the reaction.

Oxygen reduction reaction on Pt-skin Pt3V(111) fuel cell cathode: a density functional theory study

Pt-non-precious transition metals (Pt-NPTMs) alloy electrocatalysts have gained considerable attention to develop cheaper and efficient electrocatalysts for oxygen reduction reaction (ORR) in proton exchange membrane fuel cells (PEMFCs). In this report, density functional theory (DFT) has been applied to study the catalytic activity of Pt-skin Pt3V(111) electrocatalyst for ORR in PEMFCs. The results revealed that the ORR intermediates (O, OH and OOH) have lower binding energies on Pt-skin Pt3V(111) compared to pure Pt(111) surface. The ORR on Pt-skin Pt3V(111) surface proceed via OOH dissociation with an activation energy of 0.33 eV. The formation of OH is found to be the rate determining step with an activation energy of 0.64 eV, which is even lower than in pure Pt(111) surface (0.72 eV). This indicates a better performance of Pt-skin Pt3V(111) for ORR compared to pure Pt(111) surface. Moreover, the DFT results revealed that the negative formation energy of the Pt3V alloy and the positive dissolution potential shift of the surface Pt atoms revealed the better stability of Pt-skin Pt3V(111) surface over pristine Pt(111) surface. Due to the improved activity and better stability, the new Pt3V alloy electrocatalyst is very promising for the development of low-cost and efficient PEMFCs.

Revisiting the Active Sites at the MoS2/H2O Interface via Grand-Canonical DFT: The Role of Water Dissociation

MoS₂ is a promising low-cost catalyst for the hydrogen evolution reaction (HER). However, the nature of the active sites remains a subject of debate. By taking the electrochemcal potential explicitly into account using grand-canonical density functional theory (DFT) in combination with the linearized Poisson-Boltzmann equation, we herein revisit the active sites of 2H-MoS₂. In addition to the well-known catalytically active edge sites, also specific point defects on the otherwise inert basal plane provide highly active sites for HER. Given that HER takes place in water, we also assess the reactivity of these active sites with respect to H₂O. The thermodynamics of proton reduction as a function of the electrochemical potential reveals that four edge sites and three basal plane defects feature thermodynamic overpotentials below 0.2 V. In contrast to current proposals, many of these active sites involve adsorbed OH. The results demonstrate that even though H₂O and OH block "active" sites, HER can also occur on these "blocked" sites, reducing protons on surface OH/H₂O entities. As a consequence, our results revise the active sites, highlighting the so far overlooked need to take the liquid component (H₂O) of the functional interface into account when considering the stability and activity of the various active sites.

Ten Facets, One Force Field: The GAL19 Force Field for Water-Noble Metal Interfaces

Understanding the structure of the water/metal interfaces plays an important role in many areas ranging from surface chemistry to environmental processes. The size, required phase-space sampling, and the slow diffusion of molecules at the water/metal interfaces motivate the development of accurate force fields. We develop and parametrize GAL19, a novel force field, to describe the interaction of water with two facets (111 and 100) of five metals (Pt, Pd, Au, Ag, Cu). To increase transferability compared to its predecessor GAL17, the water-metal interaction is described as a sum of pairwise terms. The interaction energy has three contributions: (i) physisorption is described via a Tang and Toennies potential, (ii) chemisorption and surface corrugation rely on an attractive Gaussian term, and (iii) the angular dependence is explicitly included as a truncated Fourier series. Thirteen parameters are used for each metal surface and were fitted on 250 water adsorption energies computed at the PBE+dDsC level. The performance of GAL19 was evaluated on a set of more than 600 DFT adsorption energies for each surface, leading to an average root-mean-square deviation of only 1 kcal/mol, correctly reproducing the adsorption trends: strong on Pt and Pd but weaker on Ag, Au, and Cu. This force field was then used to simulate the water/metal interface for all ten surfaces for 1 ns. Structural analyses reveal similar tendencies for all surfaces: a first, dense water layer that is mostly adsorbed on the metal top sites and a second layer up to around 6 Å, which is less structured. On Pt and Pd, the first layer is strongly organized with water lying flat on the surface. The pairwise additive functional form allows one to simulate the water adsorption on alloys, which is demonstrated at the example of Ag/Cu and Au/Pt alloys. The water/Ag-Cu interface is predicted to be disordered with water mostly adsorbed on Cu which should exacerbate the Ag reactivity. On the contrary, incorporating Pt into Au materials leads to a structuring of the water interface. Our promising results make GAL19 an ideal candidate to get representative sampling of complex metal/water interfaces as a first step toward accurate estimation of free energies of reactions in solution at the metal interface.

Parameter-free coordination numbers for solutions and interfaces

Coordination numbers are among the central quantities to describe the local environment of atoms and are thus used in various applications such as structure analysis, fingerprints, and parameters. Yet, there is no consensus regarding a practical algorithm, and many proposed methods are designed for specific systems. In this work, we propose a scale-free and parameter-free algorithm for nearest neighbor identification. This algorithm extends the powerful Solid-Angle based Nearest-Neighbor (SANN) framework to explicitly include local anisotropy. As such, our Anisotropically corrected SANN (ASANN) algorithm provides with a fast, robust, and adaptive method for computing coordination numbers. The ASANN algorithm is applied to flat and corrugated metallic surfaces to demonstrate that the expected coordination numbers are retrieved without the need for any system-specific adjustments. The same applies to the description of the coordination numbers of metal atoms in AuCu nanoparticles, and we show that ASANN based coordination numbers are well adapted for automatically counting neighbors and the establishment of cluster expansions. Analysis of classical molecular dynamics simulations of an electrified graphite electrode reveals a strong link between the coordination number of Cs⁺ ions and their position within the double layer, a relation that is absent for Na⁺, which keeps its first solvation shell even close to the electrode.

Mononuclear Fe in N-doped carbon: computational elucidation of active sites for electrochemical oxygen reduction and oxygen evolution reactions

First principles simulations show that in Fe and N co-doped carbon, Fe coordination controls the activity for oxygen reduction and oxygen evolution reactions, and that including the electrostatic potential has a major influence at high potential.

Implicit self-consistent electrolyte model in plane-wave density-functional theory

The ab initio computational treatment of electrochemical systems requires an appropriate treatment of the solid/liquid interfaces. A fully quantum mechanical treatment of the interface is computationally demanding due to the large number of degrees of freedom involved. In this work, we develop a computationally efficient model where the electrode part of the interface is described at the density-functional theory (DFT) level, and the electrolyte part is represented through an implicit solvation model based on the Poisson-Boltzmann equation. We describe the implementation of the linearized Poisson-Boltzmann equation into the Vienna Ab initio Simulation Package, a widely used DFT code, followed by validation and benchmarking of the method. To demonstrate the utility of the implicit electrolyte model, we apply it to study the surface energy of Cu crystal facets in an aqueous electrolyte as a function of applied electric potential. We show that the applied potential enables the control of the shape of nanocrystals from an octahedral to a truncated octahedral morphology with increasing potential.

Toward a quantitative theoretical method for infrared and Raman spectroscopic studies on single-crystal electrode/liquid interfaces

In situ electrochemical infrared spectroscopy and Raman spectroscopy are powerful tools for probing potential-dependent adstructures at solid/liquid electrochemical interfaces. However, it is very difficult to quantitatively interpret the observed spectral features including potential-dependent vibrational frequency and spectral intensity, even from model systems such as single-crystal electrode/liquid interfaces. The clear understanding of electrochemical vibrational spectra has remained as a fundamental issue for four decades. Here, we have developed a method to combine computational vibrational spectroscopy tools with interfacial electrochemical models to accurately calculate the infrared and Raman spectra. We found that the solvation model and high precision level in the self-consistent-field convergence are critical elements to realize quantitative spectral predictions. This method's predictive power is verified by analysis of a classic spectroelectrochemical system, saturated CO molecules electro-adsorbed on a Pt(111) electrode. We expect that this method will pave the way to precisely reveal the physicochemical mechanism in some electrochemical processes such as electrocatalytic reactions.

Potential-Dependent CO2 Electroreduction Pathways on Cu(111) Based on an Improved Electrode/Aqueous Interface Model: Determination of the Origin of the Overpotentials

Potential-dependent CO₂ electroreduction pathways on Cu(111) are systematically studied with the aim of applying an improved electrode/aqueous interface model in this paper. The results indicate that our present defined CH₂O and CHOH pathways may be able to parallelly take place at low overpotentials. Notably, the applied potentials will not alter the optimal CO₂ reduction mechanisms. However, the presence of high overpotentials makes CO₂ electroreduction more favorable, thus explaining why high overpotentials at experiments are required during CO₂ electroreduction on Cu. Based on the potential-dependent energetics, the results suggest that COOH and CHO intermediates may be unstable at low overpotentials, in which COOH can easily change back to CO₂ and CHO can easily change back to CO, thus preventing CO₂ electroreduction. However, the high overpotentials will facilitate the formation and further electroreduction of CO and CHO. Thus, we can speculate that CO formation and then further electroreduction into CHO are the possible potential-limiting steps during CO₂ electroreduction, which are regarded as the origin of experimentally observed high overpotentials. The present comprehensive understanding on CO₂ electroreduction pathways can provide theoretical guidelines for efficiently designing Cu-based alloy electrocatalysts operated under the conditions of low overpotentials.

Catalytic Effect on CO2 Electroreduction by Hydroxyl-Terminated Two-Dimensional MXenes

Electrocatalysis represents a promising method to generate renewable fuels and chemical feedstock from the carbon dioxide reduction reaction (CO₂RR). However, traditional electrocatalysts based on transition metals are not efficient enough because of the high overpotential and slow turnover. MXenes, a family of two-dimensional metal carbides and nitrides, have been predicted to be effective in catalyzing CO₂RR, but a systematic investigation into their catalytic performance is lacking, especially on hydroxyl (-OH)-terminated MXenes relevant in aqueous reaction conditions. In this work, we utilized first-principles simulations to systematically screen and explore the properties of MXenes in catalyzing CO₂RR to CH₄ from both aspects of thermodynamics and kinetics. Sc₂C(OH)₂ was found to be the most promising catalyst with the least negative limiting potential of -0.53 V vs RHE. This was achieved through an alternative reaction pathway, where the adsorbed species are stabilized by capturing H atoms from the MXene's OH termination group. New scaling relations, based on the shared H interaction between intermediates and MXenes, were established. Bader charge analyses reveal that catalysts with less electron migration in the *(H)COOH → *CO elementary step exhibit better CO₂RR performance. This study provides new insights regarding the effect of surface functionalization on the catalytic performance of MXenes to guide future materials design.

Practical Considerations for Continuum Models Applied to Surface Electrochemistry

Modelling the electrolyte at the electrochemical interface remains a major challenge in ab initio simulations of charge transfer processes at surfaces. Recently, the development of hybrid polarizable continuum models/ab initio models have allowed for the treatment of solvation and electrolyte charge in a computationally efficient way. However, challenges remain in its application. Recent literature has reported that large cell heights are required to reach convergence, which presents a serious computational cost. Furthermore, calculations of reaction energetics require costly iterations to tune the surface charge to the desired potential. In this work, we present a simple capacitor model of the interface that illuminates how to circumvent both of these challenges. We derive a correction to the energy for finite cell heights to obtain the large cell energies at no additional computational expense. We furthermore demonstrate that the reaction energetics determined at constant charge are easily mapped to those at constant potential, which eliminates the need to apply iterative schemes to tune the system to a constant potential. These developments together represent more than an order of magnitude reduction of the computational overhead required for the application of polarizable continuum models to surface electrochemistry.