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

Non-Nernstian Effects in Theoretical Electrocatalysis

Electrocatalysis is one of the principal pathways for the transition to sustainable chemistry, promising greater energy efficiency and reduced emissions. As the field has grown, our theoretical understanding has matured. The influence of the applied potential on reactivity has developed from the first-order predictions based on the Nernst equation to the implicit inclusion of second-order effects including the interaction of reacting species with the interfacial electric field. In this review, we explore these non-Nernstian field effects in electrocatalysis, aiming to both understand and exploit them through theory and computation. We summarize the critical distinction between Nernstian and non-Nernstian effects and outline strategies to address the latter in theoretical studies. Subsequently, we examine the specific energetic contributions of the latter on capacitive and faradaic processes separately. We also underscore the importance of considering non-Nernstian effects in catalyst screening and mechanistic analysis. Finally, we provide suggestions on how to experimentally unravel these effects, offering insights into practical approaches for advancing the field.

Best practices of modeling complex materials in electrocatalysis, exemplified by oxygen evolution reaction on pentlandites

Pentlandites are natural ores with structural properties comparable to that of [FeNi] hydrogenases. While this class of transition-metal sulfide materials - (Fe,Ni)9S8 - with a variable Fe : Ni ratio has been proven to be an active electrode material for the hydrogen evolution reaction, it is also discussed as electrocatalyst for the alkaline oxygen evolution reaction (OER), corresponding to the bottleneck of anion exchange membrane electrolyzers for green hydrogen production. Despite the experimental evidence for the use of (Fe,Ni)9S8 as an OER catalyst, a detailed investigation of the elementary reaction steps, including consideration of adsorbate coverages and limiting steps under anodic polarizing conditions, is still missing. We address this gap in the present manuscript by gaining atomistic insights into the OER on an Fe4.5Ni4.5S8(111) surface through density functional theory calculations combined with a descriptor-based analysis. We use this system to introduce best practices for modeling this rather complex material by pointing out hidden pitfalls that can arise when using the popular computational hydrogen electrode approach to describe electrocatalytic processes at the electrified solid/liquid interface for energy conversion and storage.

On the Thermodynamic Equivalence of Grand Canonical, Infinite-Size, and Capacitor-Based Models in First-Principle Electrochemistry

First principles-based computational and theoretical methods are constantly evolving trying to overcome the many obstacles towards a comprehensive understanding of electrochemical processes on an atomistic level. One of the major challenges has been the determination of reaction energetics under a constant potential. Here, a theoretical framework was proposed applying standard electronic structure methods and extrapolating to the infinite-cell size limit where reactions do not alter the potential. Today, electronically grand canonical modifications to electronic structure methods, holding the potential constant by varying the number of electrons in a finite simulation cell, become increasingly popular. In this perspective, we show that these two schemes are thermodynamically equivalent. Further, we link these methods to capacitive models of the interface, in the limit that the capacitance of the charging components (whether continuum or atomistic) are equal and invariant along the reaction pathway. We benchmark the three approaches with an example of alkali cation adsorption on Pt(111) showing that all three approaches converge in the cases of Li, Na and K. For Cs, however, strong deviation from the ideal conditions leads to a spread in the respective results. We discuss the latter by highlighting the cases of broken equivalence and assumptions among the approaches.

Ab Initio-Based Modeling of Thermodynamic Cyclic Voltammograms: A Benchmark Study on Ag(100) in Bromide Solutions

Experimental cyclic voltammograms (CVs) measured in the slow scan rate limit can be entirely described in terms of the thermodynamic equilibrium quantities of the electrified solid-liquid interface. They correspondingly serve as an important benchmark for the quality of first-principles calculations of interfacial thermodynamics. Here, we investigate the partially drastic approximations made presently in computationally efficient calculations for the well-defined showcase of an Ag(100) model electrode in Br-containing electrolytes, where the nontrivial part of the CV stems from the electrosorption of Br ions. We specifically study the entanglement of common approximations in the treatment of solvation and field effects, as well as in the way macroscopic averages of the two key quantities, namely, the potential-dependent adsorbate coverage and electrosorption valency, are derived from the first-principles energetics. We demonstrate that the combination of energetics obtained within an implicit solvation model and a perturbative second order account of capacitive double layer effects with a constant-potential grand-canonical Monte Carlo sampling of the adsorbate layer provides an accurate description of the experimental CV. However, our analysis also shows that error cancellation at lower levels of theory may equally lead to good descriptions even though key underlying physics such as the disorder-order transition of the Br adlayer at increasing coverages is inadequately treated.

Fundamental Drivers of Electrochemical Barriers

We find that ion creation and destruction dominate the behavior of electrochemical reaction barriers, through grand-canonical electronic structure calculations of proton deposition on transition metal surfaces. We show that barriers respond to potential in a nonlinear manner and trace this to the continuous degree of electron transfer as an ion is created or destroyed. This explains both Marcus-like curvature and Hammond-like shifts. Across materials, we find the barrier energy to be driven primarily by the charge presented on the surface, which, in turn, is dictated by the native work function, a fundamentally different driving force than in nonelectrochemical systems.

Controlled Electrochemical Barrier Calculations without Potential Control

The knowledge of electrochemical activation energies under applied potential conditions is a prerequisite for understanding catalytic activity at electrochemical interfaces. Here, we present a new set of methods that can compute electrochemical barriers with accuracy comparable to that of constant potential grand canonical approaches, without the explicit need for a potentiostat. Instead, we Legendre transform a set of constant charge, canonical reaction paths. Additional straightforward approximations offer the possibility to compute electrochemical barriers at a fraction of computational cost and complexity, and the analytical inclusion of geometric response highlights the importance of incorporating electronic as well as the geometric degrees of freedom when evaluating electrochemical barriers.

Electrochemistry from the atomic scale, in the electronically grand-canonical ensemble

The ability to simulate electrochemical reactions from first-principles has advanced significantly in recent years. Here, we discuss the atomistic interpretation of electrochemistry at three scales: from the electronic structure to elementary processes to constant-potential reactions. At each scale, we highlight the importance of the grand-canonical nature of the process and show that the grand-canonical energy is the natural thermodynamic state variable, which has the additional benefit of simplifying calculations. We show that atomic forces are the derivative of the grand-potential energy when the potential is fixed. We further examine the meaning of potential at the atomic scale and its link to the chemical potential and discuss the link between charge transfer and potential in several situations.

A transferable prediction model of molecular adsorption on metals based on adsorbate and substrate properties

Surface adsorption is one of the fundamental processes in numerous fields, including catalysis, the environment, energy and medicine. The development of an adsorption model which provides an effective prediction of binding energy in minutes has been a long term goal in surface and interface science. The solution has been elusive as identifying the intrinsic determinants of the adsorption energy for various compositions, structures and environments is non-trivial. We introduce a new and flexible model for predicting adsorption energies to metal substrates. The model is based on easily computed, intrinsic properties of the substrate and adsorbate, which are the same for all the considered systems. It is parameterised using machine learning based on first-principles calculations of probe molecules (e.g., H₂O, CO₂, O₂, N₂) adsorbed to a range of pure metal substrates. The model predicts the computed dissociative adsorption energy to metal surfaces with a correlation coefficient of 0.93 and a mean absolute error of 0.77 eV for the large database of molecular adsorption energies provided by Catalysis-Hub.org which have a range of 15 eV. As the model is based on pre-computed quantities it provides near-instantaneous estimates of adsorption energies and it is sufficiently accurate to eliminate around 90% of candidates in screening study of new adsorbates. The model, therefore, significantly enhances current efforts to identify new molecular coatings in many applied research fields.

Understanding the Electric Double-Layer Structure, Capacitance, and Charging Dynamics

Significant progress has been made in recent years in theoretical modeling of the electric double layer (EDL), a key concept in electrochemistry important for energy storage, electrocatalysis, and multitudes of other technological applications. However, major challenges remain in understanding the microscopic details of the electrochemical interface and charging mechanisms under realistic conditions. This review delves into theoretical methods to describe the equilibrium and dynamic responses of the EDL structure and capacitance for electrochemical systems commonly deployed for capacitive energy storage. Special emphasis is given to recent advances that intend to capture the nonclassical EDL behavior such as oscillatory ion distributions, polarization of nonmetallic electrodes, charge transfer, and various forms of phase transitions in the micropores of electrodes interfacing with an organic electrolyte or ionic liquid. This comprehensive analysis highlights theoretical insights into predictable relationships between materials characteristics and electrochemical performance and offers a perspective on opportunities for further development toward rational design and optimization of electrochemical systems.

Ab Initio Simulations of Water/Metal Interfaces

Structures and processes at water/metal interfaces play an important technological role in electrochemical energy conversion and storage, photoconversion, sensors, and corrosion, just to name a few. However, they are also of fundamental significance as a model system for the study of solid-liquid interfaces, which requires combining concepts from the chemistry and physics of crystalline materials and liquids. Particularly interesting is the fact that the water-water and water-metal interactions are of similar strength so that the structures at water/metal interfaces result from a competition between these comparable interactions. Because water is a polar molecule and water and metal surfaces are both polarizable, explicit consideration of the electronic degrees of freedom at water/metal interfaces is mandatory. In principle, ab initio molecular dynamics simulations are thus the method of choice to model water/metal interfaces, but they are computationally still rather demanding. Here, ab initio simulations of water/metal interfaces will be reviewed, starting from static systems such as the adsorption of single water molecules, water clusters, and icelike layers, followed by the properties of liquid water layers at metal surfaces. Technical issues such as the appropriate first-principles description of the water-water and water-metal interactions will be discussed, and electrochemical aspects will be addressed. Finally, more approximate but numerically less demanding approaches to treat water at metal surfaces from first-principles will be briefly discussed.

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

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.