Theoretical Investigation of the Adsorbate and Potential-Induced Stability of Cu Facets During Electrochemical CO2 and CO Reduction.

The activity and product selectivity of electrocatalysts for reactions like the carbon dioxide reduction reaction (CO2RR) are intimately dependent on the catalyst's structure and composition. While engineering catalytic surfaces can improve performance, discovering the key sets of rational design principles remains challenging due to limitations in modeling catalyst stability under operating conditions. Herein, we perform first-principles density functional calculations adopting implicit solvation methods with potential control to study the influence of adsorbates and applied potential on the stability of different facets of model Cu electrocatalysts. Using coverage dependencies extracted from microkinetic models, we describe an approach for calculating potential and adsorbate-dependent contributions to surface energies under reaction conditions, where Wulff constructions are used to understand the morphological evolution of Cu electrocatalysts under CO2RR conditions. We identify that CO*, a key reaction intermediate, exhibits higher kinetically and thermodynamically accessible coverages on (100) relative to (111) facets, which can translate into an increased relative stabilization of the (100) facet during CO2RR. Our results support the known tendency for increased (111) faceting of Cu nanoparticles under more reducing conditions and that the relative increase in (100) faceting observed under CO2RR conditions is likely attributed to differences in CO* coverage between these facets.

Evaluating the stability and activity of dilute Cu-based alloys for electrochemical CO2 reduction.

Cu-based catalysts currently offer the most promising route to actively and selectively produce value-added chemicals via electrochemical reduction of CO2 (eCO2R); yet further improvements are required for their wide-scale deployment in carbon mitigation efforts. Here, we systematically investigate a family of dilute Cu-based alloys to explore their viability as active and selective catalysts for eCO2R through a combined theoretical-experimental approach. Using a quantum-classical modeling approach that accounts for dynamic solvation effects, we assess the stability and activity of model single-atom catalysts under eCO2R conditions. Our calculations identify that the presence of eCO2R intermediates, such as CO*, H*, and OH*, may dynamically influence the local catalyst surface composition. Additionally, we identify through binding energy descriptors of the CO*, CHO*, and OCCO* dimer intermediates that certain elements, such as group 13 elements (B, Al, and Ga), enhance the selectivity of C2+ species relative to pure Cu by facilitating CO dimerization. The theoretical work is corroborated by preliminary testing of eCO2R activity and selectivity of candidate dilute Cu-based alloy catalyst films prepared by electron beam evaporation in a zero-gap gas diffusion electrode-based reactor. Of all studied alloys, dilute CuAl was found to be the most active and selective toward C2+ products like ethylene, consistent with the theoretical predictions. We attribute the improved performance of dilute CuAl alloys to more favorable dimerization reaction energetics of bound CO species relative to that on pure Cu. In a broader context, the results presented here demonstrate the power of our simulation framework in terms of rational catalyst design.

Chemical Modifications of Ag Catalyst Surfaces with Imidazolium Ionomers Modulate H2 Evolution Rates during Electrochemical CO2 Reduction.

Bridging polymer design with catalyst surface science is a promising direction for tuning and optimizing electrochemical reactors that could impact long-term goals in energy and sustainability. Particularly, the interaction between inorganic catalyst surfaces and organic-based ionomers provides an avenue to both steer reaction selectivity and promote activity. Here, we studied the role of imidazolium-based ionomers for electrocatalytic CO2 reduction to CO (CO2R) on Ag surfaces and found that they produce no effect on CO2R activity yet strongly promote the competing hydrogen evolution reaction (HER). By examining the dependence of HER and CO2R rates on concentrations of CO2 and HCO3-, we developed a kinetic model that attributes HER promotion to intrinsic promotion of HCO3- reduction by imidazolium ionomers. We also show that varying the ionomer structure by changing substituents on the imidazolium ring modulates the HER promotion. This ionomer-structure dependence was analyzed via Taft steric parameters and density functional theory calculations, which suggest that steric bulk from functionalities on the imidazolium ring reduces access of the ionomer to both HCO3- and the Ag surface, thus limiting the promotional effect. Our results help develop design rules for ionomer-catalyst interactions in CO2R and motivate further work into precisely uncovering the interplay between primary and secondary coordination in determining electrocatalytic behavior.

Electrocatalytic Hydrogenation of Biomass-Derived Organics: A Review.

Sustainable energy generation calls for a shift away from centralized, high-temperature, energy-intensive processes to decentralized, low-temperature conversions that can be powered by electricity produced from renewable sources. Electrocatalytic conversion of biomass-derived feedstocks would allow carbon recycling of distributed, energy-poor resources in the absence of sinks and sources of high-grade heat. Selective, efficient electrocatalysts that operate at low temperatures are needed for electrocatalytic hydrogenation (ECH) to upgrade the feedstocks. For effective generation of energy-dense chemicals and fuels, two design criteria must be met: (i) a high H:C ratio via ECH to allow for high-quality fuels and blends and (ii) a lower O:C ratio in the target molecules via electrochemical decarboxylation/deoxygenation to improve the stability of fuels and chemicals. The goal of this review is to determine whether the following questions have been sufficiently answered in the open literature, and if not, what additional information is required:(1)What organic functionalities are accessible for electrocatalytic hydrogenation under a set of reaction conditions? How do substitutions and functionalities impact the activity and selectivity of ECH?(2)What material properties cause an electrocatalyst to be active for ECH? Can general trends in ECH be formulated based on the type of electrocatalyst?(3)What are the impacts of reaction conditions (electrolyte concentration, pH, operating potential) and reactor types?

Toward Engineering of Solution Microenvironments for the CO2 Reduction Reaction: Unraveling pH and Voltage Effects from a Combined Density-Functional-Continuum Theory.

Engineering the electrolyte microenvironment represents an attractive route to tuning the selectivity of electrocatalytic reactions beyond catalyst composition and morphology. However, harnessing the full potential of this approach requires understanding the interplay between voltage, electrolyte composition, and adsorbate binding within the electrical double layer, which is absent from the usual theoretical approaches. In this work, we apply a recently developed density functional theory (DFT)-continuum approach based on the effective screening medium method and reference interaction site model (ESM-RISM) to explore electrolyte effects with an enhanced description of the electrochemical interface. Applying this method to the binding of CO adsorbates in potassium-containing electrolytes on copper, a problem of direct relevance to CO2 electroreduction to value-added products, we show that the interdependence of voltage and pH leads to an unexpected change in adsorption site preference on Cu(001) terraces. Our findings highlight the often-overlooked importance of the electrical double-layer structure for predicting catalyst operation.

The Nature of Hydrogen Adsorption on Platinum in the Aqueous Phase.

The thermodynamic state of H2 adsorbed on Pt in the aqueous phase was determined by kinetic analysis of H2 reacting with D2 O to HDO, HD, and D2 , and by DFT-based ab initio molecular dynamics simulations of H2 adsorption on Pt(111), Pt(110), and Pt nanoparticles. Dissociative adsorption of H2 on Pt is significantly weakened in the aqueous phase compared to adsorption at gas-solid interfaces. Water destabilizes the adsorbed H atoms, decreasing the heat of adsorption by 19-22 kJ m o l H 2 - 1 while inducing an additional entropy loss of 50-70 J m o l H 2 - 1 K-1 . Upon dissociative adsorption of H2 , the average distance of water from the Pt surface increases and the liquid adopts a structure that is more ordered than before close to the Pt surface, which limits the translation mobility of the adsorbed H atoms. The presence of hydrated hydronium ions next to the Pt surface further lowers the H-Pt bond strength.

Applied Potentials in Variable-Charge Reactive Force Fields for Electrochemical Systems.

An atomic description of water dynamics and electrochemical properties at electrode-electrolyte interfaces is presented using molecular dynamics with the third generation of the charge-optimized many-body (COMB3) potential framework. Externally applied potentials in electrochemical applications were simulated by offsetting electronegativity on electrode atoms. This approach is incorporated into the variable charge scheme within COMB3 and is used to investigate electrochemical systems consisting of two Cu electrodes and a water electrolyte with varying concentrations of hydroxyls (OH-) and protons (H+). The interactions between the electronegativity offset method and the charge equilibration method in a variable charge scheme are analyzed. In addition, a charge equilibration method for electrochemical applications is proposed, where the externally applied potentials are treated by the electronegativity offset on the electrodes thus enforcing charge neutrality on the electrolyte. This method is able to qualitatively capture the relevant electrochemistry and predict consistently correct voltages with precalibration.

Effect of Surface Chemistry on Water Interaction with Cu(111).

The interfacial dynamics of water in contact with bare, oxidized, and hydroxylated copper surfaces are examined using classical molecular dynamics (MD) simulations. A third-generation charge-optimized many-body (COMB3) potential is used in the MD simulations to investigate the adsorption of water molecules on Cu(111), and the results are compared to the findings of density functional theory (DFT) calculations. The adsorption energies and structures predicted by COMB3 are generally consistent with those determined with DFT. The COMB3 potential is then used to investigate the wetting behavior of water nanodroplets on Cu(111) at 20, 130, and 300 K. At room temperature, the simulations predict that the spreading rate of the base radius, R0, of a water droplet with a diameter of about 1.5 nm exhibits a spreading rate of R0 ≈ t(0.16) and a final base radius of 3.5 nm. At 20 and 130 K, water droplets are predicted to retain their structure after adsorption on Cu(111) and to undergo minimal spreading in agreement with scanning tunneling microscopy data. When the same water droplet encounters a reconstructed, oxidized Cu(111) surface, the classical MD simulations predict wetting with a spreading rate of R ≈ t(0.14) and a final base radius of 3.0 nm. Similarly, our MD simulations predict a spreading rate of R ≈ t(0.14) and a final base radius of 2.5 nm when water encounters OH-covered Cu(111). These results indicate that oxidation and hydroxylation cause a reduction in the degree of spreading and final base radius that is directly associated with a decreased spreading rate for water nanodroplets on copper.

Effects of Sublattice Symmetry and Frustration on Ionic Transport in Garnet Solid Electrolytes.

We use rigorous group-theoretic techniques and molecular dynamics to investigate the connection between structural symmetry and ionic conductivity in the garnet family of solid Li-ion electrolytes. We identify new ordered phases and order-disorder phase transitions that are relevant for conductivity optimization. Ionic transport in this materials family is controlled by the frustration of the Li sublattice caused by incommensurability with the host structure at noninteger Li concentrations, while ordered phases explain regions of sharply lower conductivity. Disorder is therefore predicted to be optimal for ionic transport in this and other conductor families with strong Li interaction.

Poisoning effect of adsorbed CO during CO2 electroreduction on late transition metals.

Copper cathodes, at sufficiently negative potentials, are selective for hydrocarbon production during the electrochemical reduction of carbon dioxide. Other metals, such as Pt, Fe, Ni and Co, produce low to zero hydrocarbons. We employ density functional theory to examine the coverage of reaction intermediates under CO2 electroreduction conditions. A detailed thermodynamic analysis suggests that a high coverage of adsorbed CO at relevant reduction potentials blocks the metal surface sites for H adsorption, preventing C-H bond formation. The potential-dependent energetics of H adsorption and CO formation are highly sensitive to the surface coverage of the adsorbed species. The formation of surface carbon as a competing adsorption intermediate is also explored at relevant reduction potentials. CO2 electroreduction to hydrocarbons over metals active for the thermal reduction process (Fe, Ni, Co, Pt) would require a H supply for C-H bond formation that is competitive with CO* and C* at the surface.

Density functional theory study of carbon dioxide electrochemical reduction on the Fe(100) surface.

Carbon dioxide electroreduction offers the possibility of producing hydrocarbon fuels using energy from renewable sources. Herein, we use density functional theory to analyze the feasibility of CO2 electroreduction on a Fe(100) surface. Experimentally, iron is nonselective for hydrocarbon formation. A simplistic analysis of low-coverage reaction intermediate energies for the paths to produce CH4 and CH3OH from CO2 suggests Fe(100) could be more active than Cu(111), currently the only metallic catalyst to show selectivity towards hydrocarbon formation. We consider a series of impediments to CO2 electroreduction on Fe(100) including O*/OH* (* denotes surface bound species) blockage of active surface sites; competitive adsorption effects of H*, CO* and C*; and iron carbide formation. Our results indicate that under CO2 electroreduction conditions, Fe(100) is predicted to be covered in C* or CO* species, blocking any C-H bond formation. Further, bulk Fe is predicted to be unstable relative to FeCx formation at potentials relevant to CO2 electroreduction conditions.

Effects of strain, d-band filling, and oxidation state on the surface electronic structure and reactivity of 3d perovskite surfaces.

Trends in the dissociative oxygen adsorption energy and oxygen vacancy formation energy on cubic LaBO(3) and SrBO(3) perovskite (001) surfaces (where B = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, and Cu) and their dependence on strain, d-band filling, and oxidation state were examined using density functional theory in the generalized gradient approximation. The effects of strain were found to be small compared to the effects of d-band filling and oxidations state. Electronic structure descriptors such as the d-band center of the B-atom were identified for trends in the dissociative oxygen adsorption energy and for the oxygen vacancy formation energy. A chemical correlation between these two reaction energies was also identified showing the trends in these reaction energies are not independent of each other.

Effects of strain, d-band filling, and oxidation state on the bulk electronic structure of cubic 3d perovskites.

The properties of the d-band structure of the transition metal atom in cubic LaBO(3) and SrBO(3) perovskites (where B = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, and Cu) and their dependence on strain, d-band filling, and oxidation state were investigated using density functional theory calculations and atom-projected density of states. The strain dependence of the d-band width is shown to depend systematically on the size of the B atom. We show that the transition metal d-band width and center are linearly correlated with each other in agreement with a rectangular band model. A simple matrix element formalism based on the solid state table can readily predict the strain dependence of the d-band width.