Zooming into the Inner Helmholtz Plane of Pt(111)-Aqueous Solution Interfaces: Chemisorbed Water and Partially Charged Ions

Double-layer capacitance peaks: Origins, ion dependence, and temperature effects

Differential capacitance (Cdl) is arguably the most important lumped parameter of electrical double layers (EDLs). Two peaks in the Cdl profile have been commonly attributed to the crowding of counterions within the EDL. More recent studies have suggested that the two peaks are primarily caused by orientational polarization of interfacial water molecules. Herein, this recent perspective is extended by considering orientation-dependent adsorption free energy of water and tested at Au(111)-aqueous solution interfaces. Our comparative analysis of the ion dependency of the Cdl profile corroborates the view that the capacitance peaks are caused mainly by the saturation of the orientational polarization of interfacial water molecules. In addition, the temperature dependency of the Cdl profile is consistently interpreted as a consequence of the temperature effects on the orientational polarization of interfacial water.

Revealing the molecular interplay of coverage, wettability, and capacitive response at the Pt(111)-water solution interface under bias

While electrified interfaces are crucial for electrocatalysis and corrosion, their molecular morphology remains largely unknown. Through highly realistic ab initio molecular dynamics simulations of the Pt(111)-water solution interface in reducing conditions, we reveal a deep interconnection among electrode coverage, wettability, capacitive response, and catalytic activity. We identify computationally the experimentally hypothesised states for adsorbed hydrogen on Pt, HUPD and HOPD, revealing their role in governing interfacial water reorientation and hydrogen evolution. The transition between these two H states with increasing potential, induces a shift from a hydrophobic to a hydrophilic interface and correlates with a change in the primary electrode screening mechanism. This results in a slope change in differential capacitance, marking the onset of the experimentally observed peak around the potential of zero charge. Our work produces crucial insights for advancing electrocatalytic energy conversion, developing deep understanding of electrified interfaces.

Molecular and Electronic Structures at Electrochemical Interfaces from In Situ Resonant X-Ray Diffraction

An original approach to characterize electrochemical interfaces at the atomic level, a challenging topic toward the understanding of electrochemical reactivity, is reported. We employed in situ surface resonant X-ray diffraction experiments combined with their simulation using first-principle density functional theory calculations and were thus able to determine the molecular and electronic structures of the partially ionic layer facing the electrode surface, as well as the charge distribution in the surface metal layers. Pt(111) in an acidic medium at an applied potential excluding specific adsorption was studied. The presence of a positively charged counter layer composed of 1.60 water and 0.15 hydronium molecules per platinum surface unit cell at 2.8 Å from the oppositely charged Pt(111) surface was found. Our results give a unique insight into the water-metal interaction at the electrochemical interfaces.

Proton-Controlled Electron Injection in MoS2 During Hydrogen Evolution Revealed by Time-Resolved Spectroelectrochemistry

Monolayer MoS2 is an effective electrocatalyst for the hydrogen evolution reaction (HER). Despite significant efforts to optimize the active sites, its catalytic performance still falls short of theoretical predictions. One key factor that has often been overlooked is the electron injection from the conductive substrate into the MoS2. The charge transfer behavior at the substrate-MoS2 interface is nonclassical, exhibiting a liquid-gated electron injection behavior, the underlying mechanism of which remain under debate. To investigate this, we employ nanosecond time-resolved spectroelectrochemistry to probe the electron injection dynamics into monolayer MoS2 under operando HER conditions. Simultaneously, transient current measurements provide insights into the electron density at the substrate. By combining the electron density obtained from the MoS2 through spectroelectrochemical analysis with the electron density at the conductive substrate derived from transient current measurements, we explore the electron injection dynamics and characterize the current density potential (J-E) behavior at the substrate-MoS2 interface. Our findings show that the electron injection barrier and capability correlate strongly with proton concentration in the electrolyte. This relationship likely reflects the electron concentration-dependent conductivity of MoS2, where higher proton concentrations lead to fewer stray electrons before injection begins.

Structural Basis of Ultralow Capacitances at Metal-Nonaqueous Solution Interfaces

Metal-nonaqueous solution interfaces, a key to many electrochemical technologies, including lithium metal batteries, are much less understood than their aqueous counterparts. Herein, on several metal-nonaqueous solution interfaces, we observe capacitances that are 2 orders of magnitude lower than the usual double-layer capacitance. Combining electrochemical impedance spectroscopy, atomic force microscopy, and physical modeling, we ascribe the ultralow capacitance to an interfacial layer of 10-100 nm above the metal surface. This nanometric layer has a Young's modulus around 2 MPa, which is much softer than typical solid-electrolyte interphase films. In addition, its AC ionic conductivity is 4-to-5 orders of magnitude lower than that of the bulk electrolyte. The temperature dependencies of the AC ionic conductivity and thickness suggest that the soft layer is formed from metal-mediated, dipole-dipole interactions of the nonaqueous solvent molecules. The observed soft layer opens new avenues of modulating battery performance via rational design of ion transport, (de)solvation, and charge transfer in this interfacial region.

Tracking the surface structure and the influence of cations and anions on the double-layer region of a Au(111) electrode

We examined the electric double-layer (EDL) of a Au(111) electrode in a dilute perchloric acid solution using a combination of capacitance measurements and in situ scanning tunnelling microscopy under electrochemical conditions (ECSTM). The "camel-shaped" capacitance curve of the EDL is studied with different cations and anions, including their impact on the potential of zero charge (PZC). We show that the ECSTM images of thermally reconstructed and of the potential-induced surface reconstruction of Au(111) in perchloric acid electrolyte resemble previous work in sulphuric acid, displaying a herringbone pattern for a thermally reconstructed surface. Once the reconstruction is lifted, the Au(111) forms islands with an average of 1 atomic step height. When the potential is lowered below that of the PZC, the potential-induced surface reconstruction results in a more disoriented pattern than the thermally reconstructed surface. ECSTM images at different potentials are correlated with the voltammogram to understand the time and potential dependence of the surface. This correlation has led to the development of a potential window technique that can be used to reveal the surface structure of Au(111) based on observing the changes in PZC in the voltammogram. This method provides an indirect approach to understanding the surface structure without always relying on ECSTM. From the voltammogram, we also observed that anions (SO42-, CH3SO3-, ClO4-, F-) interact more strongly with the Au(111) surface than the alkali cations. The cation capacitance peak shape does not depend strongly on the identity of the alkali metal cation (Li+, Na+, K+). However, the anion capacitance peak depends strongly on the anion identity. It suggests that some level of specific adsorption cannot be excluded, even for anions that are traditionally not considered to adsorb specifically (perchlorate, fluoride).

Hierarchical Modeling of the Local Reaction Environment in Electrocatalysis

ConspectusElectrocatalytic reactions, such as oxygen reduction/evolution reactions and CO2 reduction reaction that are pivotal for the energy transition, are multistep processes that occur in a nanoscale electric double layer (EDL) at a solid-liquid interface. Conventional analyses based on the Sabatier principle, using binding energies or effective electronic structure properties such as the d-band center as descriptors, are able to grasp overall trends in catalytic activity in specific groups of catalysts. However, thermodynamic approaches often fail to account for electrolyte effects that arise in the EDL, including pH, cation, and anion effects. These effects exert strong impacts on electrocatalytic reactions. There is growing consensus that the local reaction environment (LRE) prevailing in the EDL is the key to deciphering these complex and hitherto perplexing electrolyte effects. Increasing attention is thus paid to designing electrolyte properties, positioning the LRE at center stage. To this end, unraveling the LRE is becoming essential for designing electrocatalysts with specifically tailored properties, which could enable much needed breakthroughs in electrochemical energy science.Theory and modeling are getting more and more important and powerful in addressing this multifaceted problem that involves physical phenomena at different scales and interacting in a multidimensional parametric space. Theoretical models developed for this purpose should treat intrinsic multistep kinetics of electrocatalytic reactions, EDL effects from subnm scale to the scale of 10 nm, and mass transport phenomena bridging scales from <0.1 to 100 μm. Given the diverse physical phenomena and scales involved, it is evident that the challenge at hand surpasses the capabilities of any single theoretical or computational approach.In this Account, we present a hierarchical theoretical framework to address the above challenge. It seamlessly integrates several modules: (i) microkinetic modeling that accounts for various reaction pathways; (ii) an LRE model that describes the interfacial region extending from the nanometric EDL continuously to the solution bulk; (iii) first-principles calculations that provide parameters, e.g., adsorption energies, activation barriers and EDL parameters. The microkinetic model considers all elementary steps without designating an a priori rate-determining step. The kinetics of these elementary steps are expressed in terms of local concentrations, potential and electric field that are codetermined by EDL charging and mass transport in the LRE model. Vital insights on electrode kinetic phenomena, i.e., potential-dependent Tafel slopes, cation effects, and pH effects, obtained from this hierarchical framework are then reviewed. Finally, an outlook on further improvement of the model framework is presented, in view of recent developments in first-principles based simulation of electrocatalysis, observations of dynamic reconstruction of catalysts, and machine-learning assisted computational simulations.

Parameter-Fitting-Free Continuum Modeling of Electric Double Layer in Aqueous Electrolyte

Electric double layers (EDLs) play fundamental roles in various electrochemical processes. Despite the extensive history of EDL modeling, there remain challenges in the accurate prediction of its structure without expensive computation. Herein, we propose a predictive multiscale continuum model of EDL that eliminates the need for parameter fitting. This model computes the distribution of the electrostatic potential, electron density, and species' concentrations by taking the extremum of the total grand potential of the system. The grand potential includes the microscopic interactions that are newly introduced in this work: polarization of solvation shells, electrostatic interaction in parallel plane toward the electrode, and ion-size-dependent entropy. The parameters that identify the electrode and electrolyte materials are obtained from independent experiments in the literature. The model reproduces the trends in the experimental differential capacitance with multiple electrode and nonadsorbing electrolyte materials (Ag(110) in NaF, Ag(110) in NaClO4, and Hg in NaF), which verifies the accuracy and predictiveness of the model and rationalizes the observed values to be due to changes in electron stability. However, our calculation on Pt(111) in KClO4 suggests the need for the incorporation of electrode/ion-specific interactions. Sensitivity analyses confirmed that effective ion radius, ion valence, the electrode's Wigner-Seitz radius, and the bulk modulus of the electrode are significant material properties that control the EDL structure. Overall, the model framework and findings provide insights into EDL structures and predictive capability at low computational cost.

Contrasting Views of the Electric Double Layer in Electrochemical CO2 Reduction: Continuum Models vs Molecular Dynamics

In the field of electrochemical CO2 reduction, both continuum models and molecular dynamics (MD) models have been used to understand the electric double layer (EDL). MD often focuses on the region within a few nm of the electrode, while continuum models can span up to the device level (cm). Still, both methods model the EDL, and for a cohesive picture of the CO2 electrolysis system, the two methods should agree in the regions where they overlap length scales. To this end, we make a direct comparison between state-of-the-art continuum models and classical MD simulations under the conditions of CO2 reduction on a Ag electrode. For continuum modeling, this includes the Poisson-Nernst-Planck formulation with steric (finite ion size) effects, and in MD the electrode is modeled with the constant potential method. The comparison yields numerous differences between the two modeling methods. MD shows cations forming two adsorbed layers, including a fully hydrated outer layer and a partial hydration layer closer to the electrode surface. The strength of the inner adsorbed layer increases with cation size (Li+ < Na+ < K+ < Cs+) and with more negative applied potentials. Continuum models that include steric effects predict CO2 to be mostly excluded within 1 nm of the cathode due to tightly packed cations, yet we find little evidence to support these predictions from the MD results. In fact, MD shows that the concentration of CO2 increases within a few Å of the cathode surface due to interactions with the Ag electrode, a factor not included in continuum models. The EDL capacitance is computed from the MD results, showing values in the range of 7-9 μF cm-2, irrespective of the electrolyte concentration, cation identity, or applied potential. The direct comparison between the two modeling methods is meant to show the areas of agreement and disagreement between the two views of the EDL, so as to improve and better align these models.

The pH-Induced Increase of the Rate Constant for HER at Au(111) in Acid Revealed by Combining Experiments and Kinetic Simulation

Origins of pH effects on the kinetics of electrocatalytic reactions involving the transfer of both protons and electrons, including the hydrogen evolution reaction (HER) considered in this study, are heatedly debated. By taking the HER at Au(111) in acid solutions of different pHs and ionic concentrations as the model systems, herein, we report how to derive the intrinsic kinetic parameters of such reactions and their pH dependence through the measurement of j-E curves and the corresponding kinetic simulation based on the Frumkin-Butler-Volmer theory and the modified Poisson-Nernst-Planck equation. Our study reveals the following: (i) the same set of kinetic parameters, such as the standard activation Gibbs free energy, charge transfer coefficient, and Gibbs adsorption energy for Had at Au(111), can simulate well all the j-E curves measured in solutions with different pH and temperatures; (ii) on the reversible hydrogen electrode scale, the intrinsic rate constant increases with the increase of pH, which is in contrast with the decrease of the HER current with the increase of pH; and (iii) the ratio of the rate constants for HER at Au(111) in x M HClO4 + (0.1 - x) M NaClO4 (pH ≤ 3) deduced before properly correcting the electric double layer (EDL) effects to the ones estimated with EDL correction is in the range of ca. 10 to 40, and even in a solution of x M HClO4 + (1 - x) M NaClO4 (pH ≤ 2) there is a difference of ca. 5× in the rate constants without and with EDL correction. The importance of proper correction of the EDL effects as well as several other important factors on unveiling the intrinsic pH-dependent reaction kinetics are discussed to help converge our analysis of pH effects in electrocatalysis.

Electrostatic interactions between charge regulated spherical macroions

We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann type theory with charge regulation boundary conditions. For a range of conditions we find symmetry breaking transitions from symmetric to asymmetric charge distribution exhibiting annealed charge patchiness, which results in like-charge attraction even in a univalent electrolyte-thus fundamentally modifying the nature of electrostatic interactions in charge-stabilized colloidal suspensions.

Key Strategies for Enhancing H2 Production in Transition Metal Oxide Based Photocatalysts

Transition metal oxides (TMOs) were one of the first photocatalysts used to produce hydrogen from water using solar energy. Despite the emergence of many other genres of photocatalysts over the years, TMO photocatalysts remain dominant due to their easy synthesis and unique physicochemical properties. Various strategies have been developed to enhance the photocatalytic activity of TMOs, but the solar-to-hydrogen (STH) conversion efficiency of TMO photocatalysts is still very low (<2 %), which is far below the targeted STH of 10 % for commercial viability. This article provides a comprehensive analysis of several widely used strategies, including oxygen defects control, doping, establishing interfacial junctions, and phase-facet-morphology engineering, that have been adopted to improve TMO photocatalysts. By critically evaluating these strategies and providing a roadmap for future research directions, this article serves as a valuable resource for researchers, students, and professionals seeking to develop efficient energy materials for green energy solutions.

Understanding the Effects of Electrode Material, Single Crystal Facet, and Electrolyte Ion on the Helmholtz Capacitance of Metal/Aqueous Solution Interfaces

The comprehensive interpretation of the measured differential Helmholtz capacitance curve is vital for advancing our understanding of the interfacial structure. While several possible physical effects contributing to the Helmholtz capacitance have been proposed theoretically, combining those factors to explain the experimentally observed potential-dependent capacitance profile remains a significant challenge. In this study, we employ ab initio molecular dynamics simulations to model various metal/solution interfaces. Our investigation primarily emphasizes the substantial effect of water chemisorption on the potential-dependent behavior of the Helmholtz capacitance. Additionally, we identify other critical factors that profoundly impact the Helmholtz capacitance: (1) Ions with low hydration energy hinder the availability of surface sites for water adsorption, resulting in a diminished enhancement of capacitance from water chemisorption. (2) Using large-sized ions leads to an expansion of the Helmholtz layer, causing a decrease in the Helmholtz capacitance. (3) Metal surfaces with higher affinity for water attract water adsorption at lower potentials, resulting in a lower peak potential for the differential Helmholtz capacitance curve.

Correlated surface-charging behaviors of two electrodes in an electrochemical cell

It is revealed herein that surface-charging behaviors of the two electrodes constituting an electrochemical cell cannot be described independently by their respective electric double-layer (EDL) properties. Instead, they are correlated in such a way that the surface-charging behavior of each electrode is determined by the EDL and the reaction kinetics at both electrodes. Two fundamental equations describing the correlated surface-charging behaviors are derived, and approximate analytical solutions are obtained at low and high current densities, respectively, to facilitate transparent understanding. Important implications of the presented conceptual analysis for theoretical and computational electrochemistry are discussed. A strategy of modulating the activity of one electrode by tuning EDL parameters of the other in a two-electrode electrochemical cell is demonstrated.

Comparing Ab Initio Molecular Dynamics and a Semiclassical Grand Canonical Scheme for the Electric Double Layer of the Pt(111)/Water Interface

The theoretical modeling of metal/water interfaces centers on an appropriate configuration of the electric double layer (EDL) under grand canonical conditions. In principle, ab initio molecular dynamics (AIMD) simulations would be the appropriate choice for treating the competing water-water and water-metal interactions and explicitly considering the atomic and electronic degrees of freedom. However, this approach only allows simulations of relatively small canonical ensembles over a limited period (shorter than 100 ps). On the other hand, computationally efficient semiclassical approaches can treat the EDL model based on a grand canonical scheme by averaging the microscopic details. Thus, an improved description of the EDL can be obtained by combining AIMD simulations and semiclassical methods based on a grand canonical scheme. By taking the Pt(111)/water interface as an example, we compare these approaches in terms of the electric field, water configuration, and double-layer capacitance. Furthermore, we discuss how the combined merits of the approaches can contribute to advances in EDL theory.