Ab Initio Cyclic Voltammetry on Cu(111), Cu(100) and Cu(110) in Acidic, Neutral and Alkaline Solutions

Electrochemical Nitrogen Reduction to Ammonia at Ambient Condition on the (111) Facets of Transition Metal Carbonitrides

We conducted Density Functional Theory calculations to investigate a class of materials with the goal of enabling nitrogen activation and electrochemical ammonia production under ambient conditions. The source of protons at the anode could originate from either water splitting or H2, but our specific focus was on the cathode reaction, where nitrogen is reduced into ammonia. We examined the conventional associative mechanism, dissociative mechanism, and Mars-van Krevelen mechanism on the (111) facets of the NaCl-type structure found in early transition metal carbonitrides, including Ti, V, Cr, Zr, Nb, Mo, Hf, Ta, Sc, Y, and W. We explored the catalytic activity by calculating the free energy of all intermediates along the reaction pathway and constructing free energy diagrams to identify the steps that determine the reaction's feasibility. Additionally, we closely examined the potential for catalyst poisoning within the electrochemical environment, considering the bias required to drive the reaction. Furthermore, we assessed the likelihood of catalyst decomposition and the potential for catalyst regeneration among the most intriguing carbonitrides. Our findings revealed that the only carbonitride catalyst considered here exhibiting both activity and stability, capable of self-regeneration and nitrogen-to-ammonia activation, is NbCN with a low potential-determining step energy of 0.58 eV. This material can facilitate ammonia formation via a mixed associative-MvK mechanism. In contrast, other carbonitrides of this crystallographic orientation are likely to undergo decomposition, reverting to their parent metals under operational conditions.

Tailoring Cu Electrodes for Enhanced CO2 Electroreduction through Plasma Electrolysis in Non-Conventional Phosphorus-Oxoanion-Based Electrolytes

This study presents a green, ultra-fast, and facile technique for the fabrication of micro/nano-structured and porous Cu electrodes through in-liquid plasma electrolysis using phosphorous-oxoanion-based electrolytes. Besides the preferential surface faceting, the Cu electrodes exhibit unique surface structures, including octahedral nanocrystals besides nanoporous and microporous structures, depending on the employed electrolyte. The incorporation of P-atoms into the Cu surfaces is observed. The modified Cu electrodes display increased roughness, leading to higher current densities for CO2 electroreduction reaction. The selectivity of the modified Cu electrodes towards C2 products is highest for the Cu electrodes treated in Na2 HPO3 and Na3 PO4 electrolytes, whereas those treated in Na2 H2 PO2 produce the most H2 . The Cu electrode treated in Na3 PO4 produces ethylene (23 %) at -1.1 V vs. RHE, and a comparable amount of acetaldehyde (15 %) that is typically observed for Cu(110) single crystals. The enhanced selectivity is attributed to several factors, including the surface morphology, the incorporation of phosphorus into the Cu structure, and the formation of Cu(110) facets. Our results not only advance our understanding of the influence of the electrolyte's nature on the plasma electrolysis of Cu electrodes, but also underscores the potential of in-liquid plasma treatment for developing efficient Cu electrocatalysts for sustainable CO2 conversion.

Interfacial Chemistry in the Electrocatalytic Hydrogenation of CO2 over C-Supported Cu-Based Systems

Operando soft and hard X-ray spectroscopic techniques were used in combination with plane-wave density functional theory (DFT) simulations to rationalize the enhanced activities of Zn-containing Cu nanostructured electrocatalysts in the electrocatalytic CO2 hydrogenation reaction. We show that at a potential for CO2 hydrogenation, Zn is alloyed with Cu in the bulk of the nanoparticles with no metallic Zn segregated; at the interface, low reducible Cu(I)-O species are consumed. Additional spectroscopic features are observed, which are identified as various surface Cu(I) ligated species; these respond to the potential, revealing characteristic interfacial dynamics. Similar behavior was observed for the Fe-Cu system in its active state, confirming the general validity of this mechanism; however, the performance of this system deteriorates after successive applied cathodic potentials, as the hydrogen evolution reaction then becomes the main reaction pathway. In contrast to an active system, Cu(I)-O is now consumed at cathodic potentials and not reversibly reformed when the voltage is allowed to equilibrate at the open-circuit voltage; rather, only the oxidation to Cu(II) is observed. We show that the Cu-Zn system represents the optimal active ensembles with stabilized Cu(I)-O; DFT simulations rationalize this observation by indicating that Cu-Zn-O neighboring atoms are able to activate CO2, whereas Cu-Cu sites provide the supply of H atoms for the hydrogenation reaction. Our results demonstrate an electronic effect exerted by the heterometal, which depends on its intimate distribution within the Cu phase and confirms the general validity of these mechanistic insights for future electrocatalyst design strategies.

Surface Water as an Initial Proton Source for the Electrochemical CO Reduction Reaction on Copper Surfaces

We have employed in situ electrochemical shell-isolated nanoparticle-enhanced Raman spectroscopy (SHINERS) and density functional theory (DFT) calculations to study the CO reduction reaction (CORR) on Cu single-crystal surfaces under various conditions. Coadsorbed and structure-/potential-dependent surface species, including *CO, Cu-O_ad , and Cu-OH_ad , were identified using electrochemical spectroscopy and isotope labeling. The relative abundance of *OH follows a "volcano" trend with applied potentials in aqueous solutions, which is yet absent in absolute alcoholic solutions. Combined with DFT calculations, we propose that the surface H₂ O can serve as a strong proton donor for the first protonation step in both the C₁ and C₂ pathways of CORR at various applied potentials in alkaline electrolytes, leaving adsorbed *OH on the surface. This work provides fresh insights into the initial protonation steps and identity of key interfacial intermediates formed during CORR on Cu surfaces.

Structure of the (Bi)carbonate Adlayer on Cu(100) Electrodes

(Bi)carbonate adsorption on Cu(100) in 0.1 M KHCO₃ has been studied by in situ scanning tunneling microscopy. Coexistence of different ordered adlayer phases with ( 2 ${\sqrt{2}}$ ×6 2 ${\sqrt{2}}$ )R45° and (4×4) unit cells was observed in the double layer potential regime. The adlayer is rather dynamic and undergoes a reversible order-disorder phase transition at 0 V vs. the reversible hydrogen electrode. Density functional calculations indicate that the adlayer consists of coadsorbed carbonate and water molecules and is strongly stabilized by liquid water in the adjacent electrolyte.

Screening Surface Structure-Electrochemical Activity Relationships of Copper Electrodes under CO2 Electroreduction Conditions

Understanding how crystallographic orientation influences the electrocatalytic performance of metal catalysts can potentially advance the design of catalysts with improved efficiency. Although single crystal electrodes are typically used for such studies, the one-at-a-time preparation procedure limits the range of secondary crystallographic orientations that can be profiled. This work employs scanning electrochemical cell microscopy (SECCM) together with co-located electron backscatter diffraction (EBSD) as a screening technique to investigate how surface crystallographic orientations on polycrystalline copper (Cu) correlate to activity under CO₂ electroreduction conditions. SECCM measures spatially resolved voltammetry on polycrystalline copper covering low overpotentials of CO₂ conversion to intermediates, thereby screening the different activity from low-index facets where H₂ evolution is dominant to high-index facets where more reaction intermediates are expected. This approach allows the acquisition of 2500 voltammograms on approximately 60 different Cu surface facets identified with EBSD. The results show that the order of activity is (111) < (100) < (110) among the Cu primary orientations. The collection of data over a wide range of secondary orientations leads to the construction of an "electrochemical-crystallographic stereographic triangle" that provides a broad comprehension of the trends among Cu secondary surface facets rarely studied in the literature, [particularly (941) and (741)], and clearly shows that the electroreduction activity scales with the step and kink density of these surfaces. This work also reveals that the electrochemical stripping of the passive layer that is naturally formed on Cu in air is strongly grain-dependent, and the relative ease of stripping on low-index facets follows the order of (100) > (111) > (110). This allows a procedure to be implemented, whereby the oxide is removed (to an electrochemically undetectable level) prior to the kinetic analyses of electroreduction activity. SECCM screening allows for the most active surfaces to be ranked and prompts in-depth follow-up studies.

Addressing Dynamics at Catalytic Heterogeneous Interfaces with DFT-MD: Anomalous Temperature Distributions from Commonly Used Thermostats

Density functional theory-based molecular dynamics (DFT-MD) has been widely used for studying the chemistry of heterogeneous interfacial systems under operational conditions. We report frequently overlooked errors in thermostated or constant-temperature DFT-MD simulations applied to study (electro)catalytic chemistry. Our results demonstrate that commonly used thermostats such as Nosé-Hoover, Berendsen, and simple velocity-rescaling methods fail to provide a reliable temperature description for systems considered. Instead, nonconstant temperatures and large temperature gradients within the different parts of the system are observed. The errors are not a "feature" of any particular code but are present in several ab initio molecular dynamics implementations. This uneven temperature distribution, due to inadequate thermostatting, is well-known in the classical MD community, where it is ascribed to the failure in kinetic energy equipartition among different degrees of freedom in heterogeneous systems (Harvey et al. J. Comput. Chem. 1998, 726-740) and termed the flying ice cube effect. We provide tantamount evidence that interfacial systems are susceptible to substantial flying ice cube effects and demonstrate that the traditional Nosé-Hoover and Berendsen thermostats should be applied with care when simulating, for example, catalytic properties or structures of solvated interfaces and supported clusters. We conclude that the flying ice cube effect in these systems can be conveniently avoided using Langevin dynamics.

Structural Evolution of Pt, Au and Cu Anodes by Electrolysis up to Contact Glow Discharge Electrolysis in Alkaline Electrolytes*

Applying a voltage to metal electrodes in contact with aqueous electrolytes results in the electrolysis of water at voltages above the decomposition voltage and plasma formation in the electrolyte at much higher voltages referred to as contact glow discharge electrolysis (CGDE). While several studies explore parameters that lead to changes in the I–U characteristics in this voltage range, little is known about the evolution of the structural properties of the electrodes. Here we study this aspect on materials essential to electrocatalysis, namely Pt, Au, and Cu. The stationary I–U characteristics are almost identical for all electrodes. Detailed structural characterization by optical microscopy, scanning electron microscopy, and electrochemical approaches reveal that Pt is stable during electrolysis and CGDE, while Au and Cu exhibit a voltage‐dependent oxide formation. More importantly, oxides are reduced when the Au and Cu electrodes are kept in the electrolysis solution after electrolysis. We suspect that H₂O₂ (formed during electrolysis) is responsible for the oxide reduction. The reduced oxides (which are also accessible via electrochemical reduction) form a porous film, representing a possible new class of materials in energy storage and conversion studies.

Electronic and Optical Properties of Polythiophene Molecules and Derivatives

The electronic and optical properties of polythiophene (PT) for polymer light-emitting diodes (PLEDs) were calculated using density functional theory (DFT) and time-dependent DFT. We calculated the electronic and optical properties of thiophene and PT polymers with degrees of polymerization (DP) from 2 to 30 monomers (T1–T30) and their derivatives. The associated highest occupied molecular orbital (HOMO) energy, lowest unoccupied molecular orbital (LUMO) energy, band gaps, electron orbitals, and molecular structures were determined. As the DP increased, the LUMO energy gradually decreased, and the HOMO energy gradually increased. The band gap of PT approached 2 eV as the DP of the PT polymer increased from 1 to 30. The calculations and exchange–correlation functional were verified against values in the literature and experimental data from cyclic voltammetry (redox potential) and ultraviolet-visible, photoluminescence, and ultraviolet photoelectron spectra. The color of PT PLEDs can be adjusted by controlling the DP of the polymer and the substituents.

Electrochemical biomass upgrading: degradation of glucose to lactic acid on a copper(ii) electrode

Biomass upgrading - the conversion of biomass waste into value-added products - provides a possible solution to reduce global dependency on nonrenewable resources. This study investigates the possibility of green biomass upgrading for lactic acid production by electrochemically-driven degradation of glucose. Herein we report an electrooxidized copper(ii) electrode which exhibits a turnover frequency of 5.04 s^-1 for glucose conversion. Chronoamperometry experiments under varied potentials, alkalinity, and electrode preparation achieved a maximum lactic acid yield of 23.3 ± 1.2% and selectivity of 31.1 ± 1.9% (1.46 V vs. RHE, 1.0 M NaOH) for a room temperature and open-to-atmosphere reaction. Comparison between reaction conditions revealed lactic acid yield depends on alkalinity and applied potential, while pre-oxidation of the copper had a negligible effect on yield. Post-reaction cyclic voltammetry studies indicated no loss in reactivity for copper(ii) electrodes after a 30 hour reaction. Finally, a mechanism dependent on solvated Cu²⁺ species is proposed as evidenced by similar product distributions in electrocatalytic and thermocatalytic systems.

Development of a single crystal sample holder for interfacing ultrahigh vacuum and electrochemical experimentation

Electrocatalyst surfaces prepared under ultrahigh vacuum (UHV) conditions can create model surfaces to better connect theoretical calculations with experimental studies. The development of a single crystal sample holder and inert electrochemical cells prepared with modularity and chemical stability in mind would allow for expensive single crystals to be reused indefinitely in both UHV and electrochemical settings. This sample holder shows reproducible surface preparations for single crystal samples and consistent electrochemical experiments without the introduction of impurities into the surface. The presented setup has been used as a critical piece for the characterization of Cu(111) surfaces under CO₂ electrochemical reduction reaction conditions as a test case.

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.

Selectivity Control of Cu Nanocrystals in a Gas-Fed Flow Cell through CO2 Pulsed Electroreduction

In this study, we have taken advantage of a pulsed CO₂ electroreduction reaction (CO₂RR) approach to tune the product distribution at industrially relevant current densities in a gas-fed flow cell. We compared the CO₂RR selectivity of Cu catalysts subjected to either potentiostatic conditions (fixed applied potential of −0.7 V_RHE) or pulsed electrolysis conditions (1 s pulses at oxidative potentials ranging from E_an = 0.6 to 1.5 V_RHE, followed by 1 s pulses at −0.7 V_RHE) and identified the main parameters responsible for the enhanced product selectivity observed in the latter case. Herein, two distinct regimes were observed: (i) for E_an = 0.9 V_RHE we obtained 10% enhanced C₂ product selectivity (FE_C2H4 = 43.6% and FE_C2H5OH = 19.8%) in comparison to the potentiostatic CO₂RR at −0.7 V_RHE (FE_C2H4 = 40.9% and FE_C2H5OH = 11%), (ii) while for E_an = 1.2 V_RHE, high CH₄ selectivity (FE_CH4 = 48.3% vs 0.1% at constant −0.7 V_RHE) was observed. Operando spectroscopy (XAS, SERS) and ex situ microscopy (SEM and TEM) measurements revealed that these differences in catalyst selectivity can be ascribed to structural modifications and local pH effects. The morphological reconstruction of the catalyst observed after pulsed electrolysis with E_an = 0.9 V_RHE, including the presence of highly defective interfaces and grain boundaries, was found to play a key role in the enhancement of the C₂ product formation. In turn, pulsed electrolysis with E_an = 1.2 V_RHE caused the consumption of OH^– species near the catalyst surface, leading to an OH-poor environment favorable for CH₄ production.

The Potential of Zero Charge and the Electrochemical Interface Structure of Cu(111) in Alkaline Solutions

Copper (Cu) is a unique electrocatalyst, which is able to efficiently oxidize CO at very low overpotentials and reduce CO₂ to valuable fuels with reasonable Faradaic efficiencies. Yet, knowledge of its electrochemical properties at the solid/liquid interface is still scarce. Here, we present the first two-stranded correlation of the potential of zero free charge (pzfc) of Cu(111) in alkaline electrolyte at different pH values through application of nanosecond laser pulses and the corresponding interfacial structure changes by in situ electrochemical scanning tunneling microscopy imaging. The pzfc of Cu(111) at pH 13 is identified at -0.73 V_SHE in the apparent double layer region, prior to the onset of hydroxide adsorption. It shifts by (88 ± 4) mV to more positive potentials per decreasing pH unit. At the pzfc, Cu(111) shows structural dynamics at both pH 13 and pH 11, which can be understood as the onset of surface restructuring. At higher potentials, full reconstruction and electric field dependent OH adsorption occurs, which causes a remarkable decrease in the atomic density of the first Cu layer. The expansion of the Cu-Cu distance to 0.3 nm generates a hexagonal Moiré pattern, on which the adsorbed OH forms a commensurate (1 × 2) adlayer structure with a steady state coverage of 0.5 monolayers at pH 13. Our experimental findings shed light on the true charge distribution and its interrelation with the atomic structure of the electrochemical interface of Cu.

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.

Why do RuO2 electrodes catalyze electrochemical CO2 reduction to methanol rather than methane or perhaps neither of those?

The electrochemical CO₂ reduction reaction (CO₂RR) on RuO₂ and RuO₂-based electrodes has been shown experimentally to produce high yields of methanol, formic acid and/or hydrogen while methane formation is not detected. This CO₂RR selectivity on RuO₂ is in stark contrast to copper metal electrodes that produce methane and hydrogen in the highest yields whereas methanol is only formed in trace amounts. Density functional theory calculations on RuO₂(110) where only adsorption free energies of intermediate species are considered, i.e. solvent effects and energy barriers are not included, predict however, that the overpotential and the potential limiting step for both methanol and methane are the same. In this work, we use both ab initio molecular dynamics simulations at room temperature and total energy calculations to improve the model system and methodology by including both explicit solvation effects and calculations of proton–electron transfer energy barriers to elucidate the reaction mechanism towards several CO₂RR products: methanol, methane, formic acid, CO and methanediol, as well as for the competing H₂ evolution. We observe a significant difference in energy barriers towards methane and methanol, where a substantially larger energy barrier is calculated towards methane formation than towards methanol formation, explaining why methanol has been detected experimentally but not methane. Furthermore, the calculations show why RuO₂ also catalyzes the CO₂RR towards formic acid and not CO(g) and methanediol, in agreement with experimental results. However, our calculations predict RuO₂ to be much more selective towards H₂ formation than for the CO₂RR at any applied potential. Only when a large overpotential of around −1 V is applied, can both formic acid and methanol be evolved, but low faradaic efficiency is predicted because of the more facile H₂ formation.

Energy barriers are calculated for the electrochemical CO₂ reduction reaction on the RuO₂(110) surface towards methanol, methane, formic acid, methanediol, CO and the competing H₂ formation and compared with experimental literature.

Fingerprint Voltammograms of Copper Single Crystals under Alkaline Conditions: A Fundamental Mechanistic Analysis

A critical step toward the systematic development of electrocatalysts is the determination of the microscopic structure and processes at the electrified solid/electrolyte interface. The major challenges toward this end for experiment and computations are achieving sufficient cleanliness and modeling the complexity of electrochemical systems, respectively. In this sense, benchmarks of well-defined model systems are sparse. This work presents a rigorous joint experimental-theoretical study on the single-crystal (SC) Cu/aqueous interface. Within typical computational uncertainties, we find quantitative agreement between simulated and experimentally measured voltammograms, which allows us to unequivocally identify the *OH adsorption feature in the fingerprint region of Cu(110), Cu(100), and Cu(111) SCs under alkaline conditions. We find the inclusion of hydrogen evolution reaction kinetics in the theoretical model to be crucial for an accurate steady-state description that gives rise to a negligible H* coverage. A purely thermodynamic description of the H* coverage through a Pourbaix analysis would incorrectly lead to a H* adsorption peak. The presented results establish a fundamental benchmark for all electrochemical applications of Cu.

Adsorption processes on a Pd monolayer-modified Pt(111) electrode

Specific adsorption of anions is an important aspect in surface electrochemistry for its influence on reaction kinetics in either a promoted or inhibited fashion. Perchloric acid is typically considered as an ideal electrolyte for investigating electrocatalytic reactions due to the lack of specific adsorption of the perchlorate anion on several metal electrodes. In this work, cyclic voltammetry and computational methods are combined to investigate the interfacial processes on a Pd monolayer deposited on Pt(111) single crystal electrode in perchloric acid solution. The “hydrogen region” of this Pd_MLPt(111) surface exhibits two voltammetric peaks: the first “hydrogen peak” at 0.246 V_RHE actually involves the replacement of hydrogen by hydroxyl, and the second “hydrogen peak” H_II at 0.306 V_RHE appears to be the replacement of adsorbed hydroxyl by specific perchlorate adsorption. The two peaks merge into a single peak when a more strongly adsorbed anion, such as sulfate, is involved. Our density functional theory calculations qualitatively support the peak assignment and show that anions generally bind more strongly to the Pd_MLPt(111) surface than to Pt(111).

Specific adsorption of anions is an important aspect in surface electrochemistry for its influence on reaction kinetics in either a promoted or inhibited fashion.