Theoretical Investigation of the Potential-Dependent CO Adsorption on Copper Electrodes

Despite the importance of CO adsorption in many electrocatalytic reaction mechanisms, there has been little investigation of the dependence of the free energy of CO adsorption on the applied potential. Herein, we report on the potential-dependent adsorption of CO on Cu electrodes using a grand-canonical density functional theory approach. We demonstrate that, within the working potential range of electrocatalytic CO2 reduction on Cu(111) and Cu(100), the CO adsorption strength can change by over 0.1 eV. Our analyses explain the potential dependence through an interfacial capacitance loss upon CO adsorption as well as orbital relaxation induced by the electrode potential. Via sensitivity analysis with respect to two electrolyte model parameters (solvent dielectric constant and Debye screening length), we find that the surface excess charge density is a useful descriptor of the CO adsorption free energy.

Interfacial Distribution of Nafion Ionomer Thin Films on Nitrogen-Modified Carbon Surfaces

Chemically modified carbon supports for the cathode catalyst layers of polymer electrolyte fuel cells (PEFCs) show considerable promise for boosting the oxygen reduction reaction. This study evaluated the ionomer distribution of Nafion ionomer thin films on nitrogen (N)-modified carbon surfaces along their depth direction. Neutron reflectivity (NR) measurements performed using the double-contrast technique with H2O and D2O revealed that the introduction of N functional groups to carbon thin films promoted ionomer adsorption onto the surface under wet conditions (22 °C, 85% relative humidity). Molecular dynamics (MD) simulations conducted to verify the origin of the robust contact between the ionomer and N-modified carbon surface revealed an ionomer adsorption mechanism on the N-modified carbon surfaces, which involved Coulomb interactions between the positively charged carbon surface and the ionomer side chains with negatively charged sulfonic acid groups. The positive surface charge, which was determined using the contents of the N functional groups estimated by X-ray photoelectron spectroscopy, was found to be sufficient as an impetus for ionomer adsorption. This strategy involving NR measurements and MD simulations can provide insights into the solid-ionomer interfacial structures in a cathode catalyst layer and can therefore be extensively employed in studies on PEFCs.

Unraveling Molecular Fingerprints of Catalytic Sulfur Poisoning at the Nanometer Scale with Near-Field Infrared Spectroscopy

Fundamental understanding of catalytic deactivation phenomena such as sulfur poisoning occurring on metal/metal-oxide interfaces is essential for the development of high-performance heterogeneous catalysts with extended lifetimes. Unambiguous identification of catalytic poisoning species requires experimental methods simultaneously delivering accurate information regarding adsorption sites and adsorption geometries of adsorbates with nanometer-scale spatial resolution, as well as their detailed chemical structure and surface functional groups. However, to date, it has not been possible to study catalytic sulfur poisoning of metal/metal-oxide interfaces at the nanometer scale without sacrificing chemical definition. Here, we demonstrate that near-field nano-infrared spectroscopy can effectively identify the chemical nature, adsorption sites, and adsorption geometries of sulfur-based catalytic poisons on a Pd(nanodisk)/Al2O3 (thin-film) planar model catalyst surface at the nanometer scale. The current results reveal striking variations in the nature of sulfate species from one nanoparticle to another, vast alterations of sulfur poisoning on a single Pd nanoparticle as well as at the assortment of sulfate species at the active metal-metal-oxide support interfacial sites. These findings provide critical molecular-level insights crucial for the development of long-lifetime precious metal catalysts resistant toward deactivation by sulfur.

First-principles hydration free energies of oxygenated species at water-platinum interfaces

The hydration free energy of atoms and molecules adsorbed at liquid-solid interfaces strongly influences the stability and reactivity of solid surfaces. However, its evaluation is challenging in both experiments and theories. In this work, a machine learning aided molecular dynamics method is proposed and applied to oxygen atoms and hydroxyl groups adsorbed on Pt(111) and Pt(100) surfaces in water. The proposed method adopts thermodynamic integration with respect to a coupling parameter specifying a path from well-defined non-interacting species to the fully interacting ones. The atomistic interactions are described by a machine-learned inter-atomic potential trained on first-principles data. The free energy calculated by the machine-learned potential is further corrected by using thermodynamic perturbation theory to provide the first-principles free energy. The calculated hydration free energies indicate that only the hydroxyl group adsorbed on the Pt(111) surface attains a hydration stabilization. The observed trend is attributed to differences in the adsorption site and surface morphology.

Sulfate, Bisulfate, and Hydrogen Co-adsorption on Pt(111) and Au(111) in an Electrochemical Environment

The co-adsorption of sulfate, bisulfate and hydrogen on Pt(111) and Au(111) electrodes was studied based on periodic density functional calculations with the aqueous electrolyte represented by both explicit and implicit solvent models. The influence of the electrochemical control parameters such as the electrode potential and pH was taken into account in a grand-canonical approach. Thus, phase diagrams of the stable coadsorption phases as a function of the electrochemical potential and Pourbaix diagrams have been derived which well reproduce experimental findings. We demonstrate that it is necessary to include explicit water molecules in order to determine the stable adsorbate phases as the (bi)sulfate adsorbates rows become significantly stabilized by bridging water molecules.

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.

Resolving the chemical identity of H2SO4 derived anions on Pt(111) electrodes: they're sulfate

Understanding how electrolyte composition controls electrocatalytic reactions requires molecular-level insight into electrode/electrolyte interaction. Perhaps the most basic aspect of this interaction, the speciation of the interfacial ion, is often controversial for even relatively simple systems. For example, for Pt(111) in 0.5 M H₂SO₄ it has long been debated whether the adsorbed anion is SO₄^2-, HSO₄^- or an H₃O⁺SO₄^2- ion pair. Here we apply interface-specific vibrational sum frequency (VSF) spectroscopy and theory to this problem and perform an isotope exchange study: we collect VSF spectra of Pt(111) in H₂SO₄(H₂O) and D₂SO₄(D₂O) as a function of bias and show that at all potentials they are identical. This is the most direct spectroscopic evidence to date that SO₄^2- is the dominant adsorbate, despite the fact that at 0.5 M H₂SO₄ bulk solution is dominated by HSO₄^-. This approach is based on the unique selection rule of the VSF spectroscopy and thus offers a new way of accessing general electrode/electrolyte interaction in electrocatalysis.

Interfacial Study of Nickel-Modified Pt(111) Surfaces in Phosphate-Containing Solutions: Effect on the Hydrogen Evolution Reaction

The surface modification of electrodes attracts great interest in electrocatalysis. It has often been observed that deposition of foreign adatoms on the surface of an electrode can originate a significant enhancement in the catalytic activity. For example, it has been reported that nickel deposits on Pt surfaces improve the rate of the hydrogen evolution reaction (HER, Nature Energy 2017, 2, 17031). During the deposition process of such metal adlayers, the pH and the nature of the ions in the electrolyte play an important role. Phosphate species are typically used to prepare buffer solutions in a wide range of pH. Therefore, electrolytes containing phosphate species are used in a large number of applications. However, the effect of phosphate on platinum surface modification with nickel deposits has not been studied yet. In this work, new data about the interaction of phosphate with nickel adatoms deposited on Pt(111) at pH 5 is investigated using cyclic voltammetry and infrared spectroscopy. The results show that, when nickel is in solution, the phosphate ions are adsorbed at lower potentials than in the absence of nickel. In addition, Laser-Induced Temperature Jump Technique demonstrates that nickel facilitates the adsorption of phosphate because of a shift of the potential of zero charge (pzc) toward negative potentials. This increases the magnitude of the positive electric field on the electrode surface, at a given potential E>pzc, facilitating the adsorption of anions. CO displacement technique has been also employed to obtain additional information about co-adsorbed phosphate on nickel adlayers. Finally, the HER has been studied at pH 5 in the presence of nickel, with and without phosphate in the bulk solution.

Two-dimensional π-conjugated metal-organic nanosheets as single-atom catalysts for the hydrogen evolution reaction

The design and synthesis of electrocatalysts with controllable active sites are highly desirable for producing hydrogen. Inspired by the experimental realization of two-dimensional metal-organic framework (2D-MOF) nanosheets, a computational screening was performed to find out the optimal single-atom catalysts for the hydrogen evolution reaction (HER) from 36 2D-MOFs composed of transition metal and benzene-like molecular ligands. Through first-principles calculations, six configurations were identified as potential catalysts for the HER with performance comparable to those of Pt-based catalysts due to the near-zero Gibbs free energy and favourable activation barriers originating from their appropriate d-orbital bonding mechanism with atomic hydrogen.

Maximal Predictability Approach for Identifying the Right Descriptors for Electrocatalytic Reactions

Density functional theory (DFT) calculations are being routinely used to identify new material candidates that approach activity near fundamental limits imposed by thermodynamics or scaling relations. DFT calculations are associated with inherent uncertainty, which limits the ability to delineate materials (distinguishability) that possess high activity. Development of error-estimation capabilities in DFT has enabled uncertainty propagation through activity-prediction models. In this work, we demonstrate an approach to propagating uncertainty through thermodynamic activity models leading to a probability distribution of the computed activity and thereby its expectation value. A new metric, prediction efficiency, is defined, which provides a quantitative measure of the ability to distinguish activity of materials and can be used to identify the optimal descriptor(s) ΔG_opt. We demonstrate the framework for four important electrochemical reactions: hydrogen evolution, chlorine evolution, oxygen reduction and oxygen evolution. Future studies could utilize expected activity and prediction efficiency to significantly improve the prediction accuracy of highly active material candidates.

Speciation of Adsorbed Phosphate at Gold Electrodes: A Combined Surface-Enhanced Infrared Absorption Spectroscopy and DFT Study

Despite the significance of phosphate buffer solutions in (bio)electrochemistry, detailed adsorption properties of phosphate anions at metal surfaces remain poorly understood. Herein, phosphate adsorption at quasi-Au(111) surfaces prepared by a chemical deposition technique has been systematically investigated over a wide range of pH by surface-enhanced infrared absorption spectroscopy in the ATR configuration (ATR-SEIRAS). Two different pH-dependent states of adsorbed phosphate are spectroscopically detected. Together with DFT calculations, the present study reveals that pKa for adsorbed phosphate species at the interface is much lower than that for phosphate species in the bulk solution; the dominant phosphate anion, H2PO4(-) at 2 < pH < 7 or HPO4(2-) at 7 < pH < 12, undergoes deprotonation upon adsorption and transforms into the adsorbed HPO4 or PO4, respectively. This study leads to a conclusion different than earlier spectroscopic studies have reached, highlighting the capability of the ATR-SEIRAS technique at electrified metal-solution interfaces.

Kinetically induced irreversibility in electro-oxidation and reduction of Pt surface

A mean field kinetic model was developed for electrochemical oxidations and reductions of Pt(111) on the basis of density functional theory calculations, and the reaction mechanisms were analyzed. The model reasonably describes asymmetric shapes of cyclic voltammograms and small Tafel slopes of relevant redox reactions observed in experiments without assuming any unphysical forms of rate equations. Simulations using the model indicate that the oxidation of Pt(111) proceeds via an electrochemical oxidation from Pt to PtOH and a disproportionation reaction from PtOH to PtO and Pt, while its reduction proceeds via two electrochemical reductions from PtO to PtOH and from PtOH to Pt.

Non-covalent interactions at electrochemical interfaces: one model fits all?

The shift with increasing concentration of alkali-metal cations of the potentials of both the spike and the hump observed in the cyclic voltammograms of Pt(111) electrodes in sulfuric acid solutions is shown to obey the simple model recently developed by us to explain the effect of non-covalent interactions at the electrical double layer. The results suggest that the model, originally developed to describe the effect of alkali-metal cations on the cyclic voltammogram of cyanide-modified Pt(111) electrodes, is of general applicability and can explain quantitatively the effect of cations on the properties of the electrical double layer.

Effect of thermal motion on catalytic activity of nanoparticles in polar solvent

In this study, we propose that electrode potential fluctuations due to the thermal motion of the solvent may serve to enhance the catalytic activity of nanostructures. The proposed model uses a simple, Marcus-type treatment of the statistical behavior of the solvent and the Butler-Volmer law for the instantaneous catalytic rate as a function of the electrode potential. The rapid development of probing techniques with high spatial and temporal resolution will help to further confirm and characterize the dynamical properties of nanostructures.

First-principles prediction of the equilibrium shape of nanoparticles under realistic electrochemical conditions

A first-principles model of the electrochemical double layer is applied to study surface energies and surface coverage under realistic electrochemical conditions and to determine the equilibrium shape of metal nanoparticles as a function of applied potential. The potential bias is directly controlled by adding electronic charge to the system, while total energy calculations and thermodynamic relations are used to predict electrodeposition curves and changes in surface energies and coverage. This approach is applied to Pt surfaces subject to hydrogen underpotential deposition. The shape of Pt nanoparticles under a cathodic scan is shown to undergo an octahedric-to-cubic transition, which is more pronounced in alkaline media due to the interaction energy of the pH-dependent surface charge with the surface dipole.