Nicolas
Dubouis
abc and
Alexis
Grimaud
*abc
aChimie du Solide et de l’Energie, Collège de France, UMR 8260, 75231 Paris Cedex 05, France
bRéseau sur le Stockage Electrochimique de l’Energie (RS2E), CNRS FR3459, 33 rue Saint Leu, 80039 Amiens Cedex, France
cSorbonne Université, Paris, France. E-mail: alexis.grimaud@college-de-france.fr
First published on 10th September 2019
The production of sustainable hydrogen with water electrolyzers is envisaged as one of the most promising ways to match the continuously growing demand for renewable electricity storage. While so far regarded as fast when compared to the oxygen evolution reaction (OER), the hydrogen evolution reaction (HER) regained interest in the last few years owing to its poor kinetics in alkaline electrolytes. Indeed, this slow kinetics not only may hinder the foreseen development of the anionic exchange membrane water electrolyzer (AEMWE), but also raises fundamental questions regarding the parameters governing the reaction. In this perspective, we first briefly review the fundamentals of the HER, emphasizing how studies performed on model electrodes allowed for achieving a good understanding of its mechanism under acidic conditions. Then, we discuss how the use of physical descriptors capturing the sole properties of the catalyst is not sufficient to describe the HER kinetics under alkaline conditions, thus forcing the catalysis community to adopt a more complex picture taking into account the electrolyte structure at the electrochemical interface. This work also outlines new techniques, such as spectroscopies, molecular simulations, or chemical approaches that could be employed to tackle these new fundamental challenges, and potentially guide the future design of practical and cheap catalysts while also being useful to a wider community dealing with electrochemical energy storage devices using aqueous electrolytes.
Since its discovery in 1789,3 and because it simply consists of two consecutive proton–electron transfers with no side reactions, the HER has probably been the most studied reaction in electrocatalysis.4 Until recently, numerous studies have been devoted to understanding from the material point of view how the physical properties of metallic surfaces govern the HER kinetics, as largely discussed in recent reviews.5,6 This large body of research led to a good understanding of the chemical nature of the active sites and its correlation with the HER activity.5,6 Furthermore, realizing that the HER kinetics is drastically affected by a change of pH from acidic to alkaline, the field has recently seen a regain of interest for fundamental studies on the HER, with the overall goal to master the HER activity by a fine tuning of the active site–electrolyte interactions. In light of these recent studies, questions regarding the role of the solvent structure and dynamics at the electrode–electrolyte interface, the effective nature of the reactants, and the use of spectator-additives in the HER kinetics arose and are now at the edge of the research focus. Moreover, not only such an understanding of the structure of the electrolyte in the vicinity of the active site became of prime importance for the HER, but a complete understanding of such fine effects will surely benefit unlocking the massive development of various electrochemical processes such as those involving CO2/N2/O2 reduction7–9 and their selectivity as well as developing aqueous electrolytes for Li-ion batteries.10–13
In order to reflect this past development in the field of the HER, we will first recall how the hydrogen binding energy (HBE) has emerged as a physical descriptor to rationalize the HER activity on the surface of a wide variety of catalysts in acid. We will then emphasize how recent observations highlight that some discrepancies exist regarding the universality of the HBE as a descriptor for the HER, before discussing new models that recently emerged to explain the pH/electrolyte dependence of the HER as well as new physical descriptors capturing the catalyst–electrolyte interactions. Finally, we will discuss how new experimental and computational techniques combined with strategies relying on a better understanding of interfacial interactions could unravel the factors controlling the HER kinetics.
(Volmer) H3O+ + e− + * = H* + H2O | (1.1) |
(Heyrovsky) H* + H3O+ + e− = H2 + H2O | (1.2) |
(Tafel) 2H* = H2 | (1.3) |
Owing to its outstanding electrocatalytic activity, most of the studies therefore focus on the HER mechanism on the surface of Pt in acidic media. A typical cyclic voltammogram (CV) recorded for a polycrystalline platinum (Pt(pc)) electrode in acidic media is shown in Fig. 1.
Fig. 1 Cyclic voltammogram recorded on a Pt(pc) electrode at 50 mV s−1 in a 0.5 M H2SO4 solution degassed with argon. |
This typical cyclic voltammogram exhibits 4 regions. At high potentials (>0.65 V vs. RHE), the Pt surface is oxidized/reduced with the subsequent formation of Pt–OH and/or Pt-oxide that can be concomitant with the specific adsorption of anions on the surface of the platinum. The second region at potentials comprised between 0.35 V and 0.65 V vs. RHE is referred to as the “double layer region” as no faradaic process occurs. At potentials between 0.05 V and 0.35 V vs. RHE, the third region is usually denoted as the hydrogen underpotential deposition (HUPD) region and is characterized by the discharge of protons on the Pt surface following the reaction: H+ + e− + Pt* = Pt − HUPD. Note that the precise nature of this phenomenon was recently reviewed, and may involve cation adsorption as well as OH desorption, which will be discussed later on.15–17 Finally, at more negative potentials, H2 is evolved in the HER region.
Three main features can be seen in the HUPD region: one broad envelope, a couple of redox peaks at ∼0.20 V vs. RHE followed by a third peak at ∼0.10 V vs. RHE. While these phenomena were already described by comparing CVs recorded on single-crystal surfaces in the 60's and tentatively assigned to the hydrogen UPD on different facets of platinum,18 this assignment was confirmed in the 80's owing to the development by Clavilier of new methods for preparing high-quality surface single crystals.19 Hence, the large envelope is related to the HUPD on Pt(111) facets, while the peaks at ∼0.10 V vs. RHE and ∼0.20 V vs. RHE are related to HUPD on the Pt(110) and the Pt(100) sites, respectively. It should be noted that the underpotential deposition of hydrogen is not unique to platinum and is known to occur on other metallic surfaces.20 Due to the ultra-fast kinetics of the HER on platinum electrodes, the faradaic current is often limited by the mass-transport of H2 generated during the HER, even when a rotating-disk electrode (RDE) apparatus is employed.21,22 This limitation hinders a proper fitting with microkinetic models14 and has thus prevented the catalysis community from achieving a fine understanding of the HER mechanisms under acidic conditions. Despite this limitation, by measuring the HER and hydrogen oxidation reaction (HOR) currents in the micropolarization zone (small potential window across the RHE potential) at low temperature, Marković et al. were able to discriminate different activation energies for the HER on different Pt-single crystal facets, with the activation energy ranging from 9.5 kJ mol−1 for Pt(110) to 18 kJ mol−1 for Pt(111) (Fig. 2a).21
Fig. 2 (a) Arrhenius plot of exchange current densities for Pt(110) (circle), Pt(100) (square) and Pt(111) obtained by cyclic voltammetry in a 50 mM H2SO4 solution (reproduced from ref. 21 with permission from the American Chemical Society, Copyright (1997)). (b) Snapshot of a MD-DFT simulation showing the presence of two types of adsorbed hydrogen on a Pt(111) electrode. Reproduced from ref. 23 with permission from Elsevier, Copyright (2007). |
From the Tafel slope measurements, which are themselves subject to practical limitations (ohmic drop correction, gas saturation, impurities, fitting parameters, etc.), the recombination step (Tafel step with a slope of ∼30 mV per decade) was determined to be the rds for Pt(110) while the Heyrovsky (Tafel slope of ∼40 mV per decade) step is believed to be the rds for Pt(100). For Pt(111), the value obtained (74 mV per decade) does not match any of the expected values from microkinetics analysis. For polycrystalline platinum, Tafel slope values of around 30 mV per decade are usually measured, which match well with the expected values for the Tafel step being the rds.4
While the growth of high quality single crystals has been tremendous for the understanding of the HER on Pt, the development of computational methods has also brought a deeper understanding for the HER/HOR. For instance, combining molecular dynamics (MD) with density-functional theory (DFT), the HOR was predicted to follow a Heyrovsky–Volmer mechanism on Pt(111) (Fig. 2b).23 Furthermore, these calculations also confirmed the presence of two types of hydrogen atoms adsorbed on the Pt(111) surface. The first one is directly adsorbed from the H2 molecule on the bridge sites of the Pt slab and is found not to be reactive, thus its tentative assignment to the HUPD. The second one, adsorbed on top of the Pt(111) slab corresponds to a proton in solution being reductively adsorbed on the femtosecond time scale on the electrode surface. While theoretical insights help the understanding of the mechanism in a very short timescale, they must be validated with experimental results to assess their validity. Toward that goal, femto spectro-electrochemical techniques were recently used to correlate the activation energies measured by Marković et al.21 with the rate of electron transfer along the Pt–HUPD bond, showing that this transfer occurs at the fs timescale.24 Doing so, the interfacial solvent reorganization which is taking place at a greater timescale prior to any electron transfer was also shown to be critical for the HUPD kinetics.24
Fig. 3 Dependence of the exchange current density of the HER on the (a) measured (reproduced from ref. 26 with permission from Elsevier, Copyright (1972)) and (b) calculated (reproduced from ref. 23 with permission from Electrochemical Society, Inc., Copyright (2005)) metal hydrogen binding energy and (c) difference of work-function between a hydrogenated and a wetted surface (reproduced from ref. 28 with permission from PCCP Owner Societies, Copyright (2017)). |
Nevertheless, as mentioned by Trasatti, one of the challenges pertaining to this approach is the reliability in measuring the physical properties of pure metals omitting the presence of the electrolyte and adsorbed species, alike in the work of Trasatti, or even of a passivation layer (e.g. Mo or W are covered by an oxide layer while platinum is covered by the HUPD). In light of these limitations and thanks to recent advances in ultra-high vacuum (UHV) techniques, the work-function of realistic surfaces was recently reinvestigated by Zeradjanin et al. as a physical descriptor for the HER including the role of adsorbed water.28 Considering that the free energy of activation ΔG* for proton adsorption is a combination of (1) the work to transfer a proton from the Outer-Helmholtz plane (OHP) to the Inner-Helmholtz Plane (IHP), (2) the proton desolvation and the energy resulting from the difference of potential between the metal and the IHP and (3) the energy to transfer electrons from the metal (at its Fermi level) to the IHP, it was thus qualitatively predicted that the proton adsorption energy becomes more endergonic (weaker E(M − H)) when the electrode (with an interfacial water layer) work function is increased.28 As a direct consequence, increasing the work function should increase ΔG* and thus decrease the exponential term of the electron transfer rate constant, usually described within the framework of the transition-state theory (TST) as: where κel is the electronic transmission coefficient (0 ≤ κel ≤ 1) which accounts for the tunneling-probability of the electron increasing with the electrode–reactant coupling, Γn is the nuclear tunneling factor (Γn ≥ 1) which corrects the rate expression taking into account molecules/ions that react without entirely surmounting the classical electrochemical free energy barrier, νn is the nuclear frequency factor which represents the frequency at which the reacting species approach the transition state thanks to the solvent and bond coordinate motions and ΔG* is the activation free energy.28 However, as deduced from this expression, increasing the work function not only affects the activation free energy but also results in an increase of the pre-exponential constant (greater reactant–electrode coupling resulting in larger values of both κel and Γn), which is often omitted but counterbalances the increase in activation energy. Thus, a more precise volcano plot is obtained when plotting the exchange current density versus the difference of work functions between hydrogenated and non-hydrogenated metals in the presence of an interfacial water layer (Fig. 3c). Finally, it is worth mentioning that this discussion on the work function stands only to assess the ease with which hydrogen binds to a given material. Nevertheless, it does not capture the effect of other parameters such as modification of the electrolyte and/or interfacial interactions which can control the HER kinetics for the best catalysts.
Interestingly, despite some approximations and limitations, the HBE theory remains largely used by the catalysis community since it can easily act as a guide for the discovery of new catalysts. For instance, it was used to explain the great performances of new cheap catalysts, such as MoS2 for the HER,29,30 and helped to clarify the nature of active sites which led to the design of new chalcogenide compounds enriched in active sites.31
Hence, the HER under acidic conditions is now relatively well understood. Nevertheless, when increasing the pH of the electrolyte from acidic to alkaline conditions, the HER and HOR performances of numerous catalysts were reported to dramatically decrease. For instance, platinum is known to undergo a decrease from 2 to 3 orders of magnitude of the HER exchange current density from pH = 0 to pH = 13.32,33 This pH effect thus renders the use of the HBE theory to find a simple activity descriptor based on a sole material's properties difficult. The second part of this perspective will therefore focus on recent studies discussing the preponderant role of other adsorbed species, such as cations or hydroxide anions, as well as the reorganization of the electrolyte on the surface of the electrode in the HER kinetics.
Fig. 4 (a) CVs recorded at different pH values on a polycrystalline platinum electrode (reproduced from ref. 35 with permission from Springer Nature, Copyright (2015)) from which the hydrogen binding energy is extracted, (b) overpotential for the HER at −1 mA cm−2 on a polycrystalline platinum electrode in different buffered electrolytes, plotted versus the HBE obtained from electrochemical measurements (reproduced from ref. 35 with permission from Springer Nature, Copyright (2015)), (c) calculated free energy diagram for the HER at pH = 0 and no overpotential for Pt, MoS2 and Au surfaces (data are from ref. 27) and (d) cyclic-voltammogram at pH = 2 (H2SO4 18 mM + K2SO4 82 mM degassed with Ar) on Pt, amorphous MoS2 and Au electrodes with a scan rate of 50 mV s−1 with oxidation currents in yellow and reduction current in blue, dashed lines indicating the j = 0 axis. |
While being attractive, some concerns have thus been recently raised about the HBE theory.37 Indeed, as the H+/H2 couple is found to be reversible on the Pt surface for instance,21 all the HER elementary steps are reversible and the HOR intermediates are certainly identical.32 Therefore, the HBE should control the HER and the HOR in a similar fashion. However, exploring other materials with relatively good HER performances and a small HBE value such as MoS2 (ref. 29 and 38) (Fig. 4c), it is observed that MoS2 poorly catalyzes the HOR (Fig. 4d). Hence, the absence of reversibility on a low-HBE material suggests that the energy and the nature of the intermediates for the HOR and HER reactions are not identical, which may originate from changes of some material properties (defects, electronic properties, etc.) or of the double-layer structure/composition as a function of the applied potential, thus preventing the use of HBE as a sole descriptor of the HER kinetics.
Furthermore, as previously mentioned, HBE measurement by either UHV or electrochemical measurements (HUPD) involves hydrogen atoms that are likely to be different from those involved as reactive intermediates in the HER. One obvious manifestation of this difference is that while the HUPD deposition potential on Pt(111) follows a Nernstian behavior with pH (60 mV per unit) and so the HBE determined from electrochemical measurements is not expected to change with pH, the HER activity drastically decreases by 2–3 orders of magnitude changing from pH = 0 to pH = 13 for the same surface.33 Finally, the effect of buffer solution should also be investigated when performing measurements at intermediate pH values, as it can also impact the proton transfer kinetics.39
At this point, we have shown how the identification of reaction intermediates on single-crystal model electrodes has triggered the quest for finding universal descriptors for the HER activity. While the HBE was initially seen as a good candidate, this theory cannot predict how the electrode–electrolyte interfacial structure impacts the performances of a given surface. Indeed, puzzled by the severe loss of activity of Pt when used under neutral or alkaline conditions, the electrocatalysis community recently investigated questions previously overlooked: what is the rate determining step for the HER under alkaline conditions? Which is the proton donor? What is the role of spectator species? How is the catalyst surface modified?
(Volmer) H2O + e− + * = H* + OH− | (1.4) |
(Heyrovsky) H2O + e− + H* = H2 + OH− | (1.5) |
Moreover, it was observed that the HER kinetics over Pt(111) under alkaline conditions is drastically enhanced by the presence of oxophilic groups, such as Pt-islands (defects) and Ni(OH)2 on the catalyst surface (Fig. 5a),33 which could be explained as originating from an easier H2O dissociation. Following this observation, it was confirmed that the HER activity for Pt(111) decorated by different 3d transition metal hydroxides M(OH)2 follows the Brønsted–Evans–Polanyi principle and the oxophilic groups' affinity to OHads should be neither too strong (surface poisoning) nor too weak (no binding) to promote the HER kinetics.40
Fig. 5 (a) HER activity in alkaline solutions measured on a Pt(111) single crystal electrode, customized with different oxophilic groups and (b) “bi-functional” mechanism responsible for the enhanced activity when groups that promote water dissociation are present on the Pt(111) surface (reproduced from ref. 33 with permission from AAAS, Copyright (2011)). |
Thus, in alkaline electrolytes, the cleavage of water O–H bond and the transport of OH− from the catalyst surface to the bulk of the electrolyte are likely to be part of the limiting step for this so-called bi-functional mechanism. A similar conclusion was recently made for copper electrodes for which a larger HER activity was measured with more oxophilic surface created by applying a mechanical or electrochemical surface treatment.41 Overall, switching from H3O+ to H2O as a proton donor could simply hamper the Volmer step,42–44 explaining the difference in activity between acidic and alkaline conditions. Indeed, as stated earlier, the Volmer step in alkaline electrolytes includes a water O–H bond dissociation step, while it is not true for acidic electrolytes. If this dissociation is in a fast equilibrium, the HER kinetics should only be determined by the proton–electron transfer rate. Nevertheless, as the activation energy for water auto-ionization (around 75 kJ mol−1 at 300 K5,45) is found to be larger when compared to the one for the HER in alkaline electrolytes (around 30–45 kJ mol−1 for Pt catalysts5,32), this could suggest that the water auto-ionization step is coupled with the electron transfer (as recently supported by DFT calculations46), and thus is partly responsible for the slow HER kinetics measured in alkaline electrolytes.
Alike for the HBE, recent studies however interrogate the importance of the surface oxophilicity and of the bi-functional mechanism in the HER kinetics.36,47,48 Indeed, comparing the HER activity of Pt/C, PtNi/C and PtNi/C treated in acid to remove Ni from the surface of the nanoparticles revealed that while the presence of nickel increases the HER kinetics in alkaline media when compared to Pt/C,47 removing it from the surface does not hinder the HER and HOR kinetics for PtNi/C. Thus, the role of the 3d metal was rather assigned to a modification of the HBE. Additionally, it was demonstrated that in unbuffered solutions, the two reduction events previously attributed to a different kinetics for the H3O+ and the H2O reduction42 may simply arise from the modification of the surface pH with the applied potential.49 Hence, assuming that the surface pH is solely determined by the applied potential on Pt electrodes (Nernst equation), the peculiar electrochemical response observed at large current densities is well explained by the modification of the surface pH. However, this peculiar shape has been observed for a wide variety of surfaces on which the H+/H2 couple is not reversible.38,42 Also, it may not explain the role of spectator species observed in the study of alkaline HER kinetics.
Indeed, even more revealing than the surface decoration of Pt with M(OH)2, the addition of Li+ cations in the electrolyte was found to be beneficial for the HER in the presence of surface oxophilic groups (Fig. 5a).33 This effect was initially rationalized by the anchoring of Li+ to the M(OH)2 moieties promoting the water-dissociation through non-covalent interactions (Fig. 5b).33 While the role of these “interfacial” interactions was previously discussed for the oxygen reduction reaction (ORR) and the HOR,50 this first observation for the HER highlighted the role of “spectator” species from the electrolyte and prompted novel fundamental studies.5 In particular, unlike the previous model which was based on non-specific adsorption of cations,50 recent DFT and electrochemical measurements on Pt single crystal electrodes revealed that cation specific adsorption can occur.16,17,51 Hence, both computational and experimental studies were carried out on Pt(111), Pt(100) and Pt(110) to tentatively explain the unusual non-Nernstian shift of 50 mV (vs. SHE) for the UPD peak observed for the two latter facets16 and revealed that considering the sole adsorption of hydrogen in the HUPD region cannot reproduce the experimental cyclic voltammograms recorded for Pt(110) and Pt(100).15,16 Rather, OH adsorption should also be taken into account (or in fact the co-adsorption with H2O, as discussed later), resulting in the equation
xOHδ+ad + (x + 1)H+ + [(1 + δ)x + 1]e− = H* + xH2O, | (1.6) |
Later on, the effect of the specific adsorption of cations onto the Pt(100) and Pt(110) was investigated to explain the remaining charge on OHad. Indeed, K+ adsorption was found to shift the OH−/H2O adsorption (and so desorption) to a higher potential by destabilizing OHad on the surface of the electrode.16 This destabilization was ascribed to an increase of the Pt–O bond lengths (Fig. 6a) as well as to a partial charge retention by the adsorbed cation, which is compensated by the presence of a negative charge on the adsorbed OHad and H2Oad (eqn (1.6) and Fig. 6b). Both of these effects eventually lead to a more positive potential for the adsorption of OH from the bulk of the electrolyte as the concentration of K+ is increased, explaining the non-Nernstian shift of the “hydrogen” peaks with the KOH concentration. This hypothesis was further generalized to more complex Pt surfaces such as Pt(533) that contains Pt(111) terraces and (100) steps.17 Finally, the role of the cation has been investigated on Pt(110) steps for which a non-Nernstian positive shift of the “UPD peaks” was only observed for pH > 3 following the ionic radius trend in the alkali series: Li+ < Na+ < K+ < Cs+ (Fig. 6c).17 The absence of the cation effect at low pH was explained as the consequence of the potential of zero charge (pzc) being more negative than the potentials of the UPD region for Pt(110), so no cation is expected to be adsorbed in the steps.17 In contrast, at higher pH the influence of the cation was explained by a size effect: the larger the cation, the more charge it retains (as the charge retained is proportional to the change in surface-normal dipole moment) and the more OHad is destabilized.17,51
Fig. 6 (a) Pt–O bond lengths and (b) Bader charge for adsorbed OH and H2O molecules on Pt(100) calculated using DFT at different adsorbed potassium coverages (reproduced from ref. 16 with permission from the American Chemical Society, Copyright (2016)), and (c) experimental position of the OHad/H* peak on Pt(533)ad at different pH values with different alkali cation contents (reproduced from ref. 17 with permission from John Wiley and Sons, Copyright (2017)). |
Unfortunately, for most of these fundamental studies, the impact of the H*/OHad exchange on the HER kinetics is not discussed. The few studies focusing on the role of alkali concentration and pH in the HER performances of Pt or PtNi52,53 could show that even if modest, the nature of the cation impacts the HER kinetics as follows: Li+ > Na+ > K+. This trend corresponds to a smaller destabilization of the OHad and so virtually a greater interfacial oxophilicity (and easier water dissociation, as shown by the group of Marković40,41) of the surface when the cation is specifically adsorbed, as stated by Koper's group.17,51 If not specifically adsorbed, the cation effect could be rationalized by the formation of an OHad–(H2O)–AM+ adduct, as discussed above, as the cation would promote the removal of OHad from this adduct and its transport out of the double-layer. This effect is expected to be enhanced in the presence of more acidic cations which would also explain the greater HER activity measured in the presence of smaller alkali cations. Overall, further studies will be necessary to better understand this interfacial structure–reactivity dependence.
While most of the research studies carried out recently have been focusing on the role of the cationic species at the solid/liquid interface, the seminal work conducted by Trasatti more than 50 years ago suggested that the water orientation at the interface could be responsible for the discrepancies in activity over different metallic surfaces.26 Pitfalls exist in experimentally assessing such sensitive effects and the accuracy of theoretical calculations suffers as well from the high computational cost of including explicit solvent molecules.54 However, while it was shown by UHV and DFT calculation that interfacial water does not drastically affect the HUPD energetics,15,16 it is expected to largely influence the alkali cations and hydroxide adsorption.16 Therefore, investigating the water structure (orientation, chemical environment) and dynamics (lifetime in the IHP, structure reorganization along charge transfer) at the electrode–electrolyte interfaces will be key to unravel the mechanisms by which interfacial interactions impact the HER activity. These observations on the role of the spectator species lead us to speculate that the use of bare physical descriptors such as the HBE cannot fully capture the HER kinetics which requires inclusion of dynamic effects occurring at the electrochemical interface. We believe that, within the experimental framework developed by Zeradjanin et al.,28 it would be critical to reinvestigate, using different electrolytes, physical descriptors such as work function difference between non-hydrogenated and hydrogenated wetted surfaces.
Fig. 7 (a) Effect of the pH and pzfc on the double-layer rigidity during the HER, (b) plot of HER exchange current density as a function of the potential of zero charge for several metallic electrodes (pzc values are taken from ref. 61 and exchange current densities under acidic conditions from ref. 26). |
While the pzfc model should properly work for surfaces having similar properties (i.e. similar active sites), a more precise picture of the influence of the pzc on the proton-transfer kinetics for a given surface must also consider eventual modifications of the pre-exponential factor of the kinetic constant which can reflect a modification of the proton/electron tunneling or the H+ concentration at the interface, among numerous parameters.28,60 Other models considering the double-layer ordering were thus developed to explain this pH-dependent activity. Hence, it was suggested that the activation energy for proton adsorption is mainly driven by a loss of entropy while protons are transferred from the bulk of the electrolyte to the OHP, a loss which increases with the pH.62 This entropy loss is compensated by the enthalpy term with a delay, so that the variation in the Gibbs free energy remains null, satisfying the equilibrium conditions.62 Such a proposition relies on the assumption that the liquid water structure between the OHP and the electrode is not affected by pH variations,62 which should be discussed in light of the study discussed above on the influence of interfacial water structure on the HER activity (pzfc theory).57
Furthermore, some computational efforts attempting to explicitly model the electrode–electrolyte interface are worth highlighting. For instance, a seminal study on the water dissociation on a Pd surface by Filhol et al. proposed that at least three distinct water layers at the electrochemical interface are required for the H2O reductive (resp. oxidative) dissociation: a layer with proton donors (resp. hydroxide donors), a layer to accept the generated hydroxide (resp. protons) by H-bonding and a layer to diffuse the as-generated species in solution.63 Moreover, thanks to recent advances in DFT calculations, a reduced charge of +0.7e for protons located in the OHP at the Pt(111)–H2O interface was found, which was explained by a charge transfer between the electrode and the interfacial solvent molecules.64 Interestingly, a bulk behavior (protons exhibiting a unit charge) was found for protons located in the third layer of water. These differences in charge highlight the importance of understanding the electrode–electrolyte interface at the molecular scale to precisely assess the charge transfer energetics. Finally, we should recall that all these models share the same assumption: water self-ionization is in a fast equilibrium and H3O+ should be considered as the reactant at every pH. However, since the H2O self-ionization activation energy exceeds by a factor of ∼2 the one measured for the HER on Pt,46 the validity of this hypothesis has to be demonstrated. Overall, both water-dissociation (bifunctional mechanism) and interfacial electric field (pzfc theory) can simultaneously play a role in the HER kinetics in alkaline media.
Naively speaking, considering all the parameters reviewed above, an ideal catalyst should possess the following properties. First, the active site must be able to bind H neither too strongly nor too weakly (HBE theory). Regarding its pzc, it should be high enough to enhance the pre-exponential factor of the rate constant (more protons at the interface and more coupled to the electrode), but low enough to minimize the double-layer rigidity and facilitate charge transport through the double layer. With these properties presumably competing with each other, finding an optimal catalyst with a single active site might prove difficult. Rather combining with the help of synthesis routes such as (de)alloying, surface decoration, etc. different active sites with different functionalities is probably a more promising strategy to follow, bearing in mind that some species contained in the electrolyte are taking part in the reaction mechanism. Taking into account the surface properties, non-reactive species such as Lewis acids can help tune the energies of the intermediates by for instance promoting the H2O self-ionization and reducing the overall HER activation energy. Nevertheless, these effects have mostly been overlooked till recently, the reason being the relative lack of experimental and theoretical studies that could draw a molecular understanding of the catalyst–electrolyte structure. Thus, in the last part of this perspective, in situ spectroscopic techniques and computational methods recently developed and that could prove useful to probe the impact of water interfacial structure on its reactivity will be highlighted.
Thus, the determination of the water structure (density, orientation, etc.) in the first molecular layers at the surface of the electrode, and especially in the IHP, requires the use of other in situ analytical techniques. Among them, surface X-ray scattering was employed to determine the potential-dependent density and orientation of water molecules in NaF aqueous solution at the Ag(111)–electrolyte interface. Doing so, it was demonstrated that under polarization, the first layer of water at the interface exhibits a larger density than in the bulk of the electrolyte, which implied a disruption of the hydrogen bonding network.68 Additionally, this observation was correlated with the presence of a large electric field (∼107 V cm−1) which attracts and orientates water molecules at the interface, switching from hydrogen pointing toward the electrode (H-down) at potentials more negative than the pzc to oxygen pointing toward the surface (H-up) at potentials above the pzc. This potential-dependent orientation of water molecules at the interface was also confirmed by laser-jump measurements on Pt single crystals.55 Furthermore, combining X-ray absorption spectroscopy (XAS) measurements at the oxygen K-edge in total fluorescence yield (TFY, bulk sensitive) and total electron yield (TEY, surface sensitive) modes with ab initio simulations, useful information about the water structure at a gold interface upon polarization could be gained.56 Hence, at potentials more negative than the pzc, the H-down orientation was confirmed and the amount of dangling bonds was found to increase, leading to a disruption of the water H-bond network (Fig. 8a and b).56 Finally, a very recent study employing in situ Raman spectroscopy in the presence of shiners confirmed the decrease of hydrogen-bond donor water molecules at decreasing potential for a gold electrode (Fig. 8b),69 as well as the densification of the interface and the favored “one H-down” orientation toward the surface upon negative polarization, while 2H-down water molecules could be spotted only at very negative potentials.
Fig. 8 (a) Experimental O K-edge total electron yield XAS spectra of water, collected at a gold electrode under polarization and (b) population of H-bonded water molecules deduced from these spectra at different polarization (adapted from ref. 56 with permission from AAAS, Copyright (2014)). (c) Number of hydrogen-bond donors of interfacial water molecules at different potentials on a gold electrode calculated from in situ Raman measurement (reproduced from ref. 69 with permission from Springer Nature, Copyright (2019)). |
While being very useful, these techniques often require the use of model non-reactive electrodes, as in the examples cited above. Nevertheless, it is worth mentioning that similar techniques were also employed on more practical systems attempting to probe the presence/absence of chemical transformations of the catalyst surface. Hence, using in situ XAS experiments on a platinum electrode in sulfuric acid, doubts were raised about the nature of the phenomena occurring in the oxidation region for the Pt electrochemistry.70 Unlike the commonly accepted explanation that the anodic current in this region is, at least partly, related to the formation of Pt–O or Pt–OH groups on the surface, results obtained by in situ XAS experiments suggest that the change observed in the oxygen K-edge XAS spectra upon positive polarization is rather due to the specific adsorption of sulfates.70 This conclusion was further confirmed by in situ APXPS measurements which revealed the specific adsorption of anions onto platinum which is only oxidized at potentials higher than +1 V vs. SHE.71 In addition to these measurements focusing on Pt under anodic polarization, a recent operando study on Pt in alkaline electrolytes72 confirmed the absence of PtO or PtO2 species during the HER. Rather, the formation of Pt–H bonds that are not stripped in the UPD region was observed. Unfortunately, in the presence of faradaic electrochemical processes, information about the solvent structure is often difficult to obtain and emphasis was paid to the modifications occurring on the surface of the electrodes.
Moreover, the dynamics observed for the defect reorganization in the adlayer is much slower than in the bulk electrolyte (∼1 ns vs. ∼5 ps), which is very likely to impact electrocatalytic reactions. Indeed, it can result in a larger concentration of hydronium ions in the first adlayer and can also hinder the diffusion of adsorbed species on the Pt surface, thus explaining that the Tafel step is the rds for the HER on Pt in acidic media.73 Furthermore, conducting similar simulations, the same group also revealed that the transfer of a water molecule from the bulk of the electrolyte to the adlayer requires a reorganization of the adlayer. Hence, this reorganization results in slow dynamics for the exchange of water molecules from the bulk to the interface, with a lifetime of ∼40 ns for water molecules in the adlayer.74 While these solvent reorganization phenomena happen at timescales (ps/ns) faster than the HER in alkaline media, it was recently observed that even if the “pure electron transfer” timescale for the HER (under acidic conditions) is less than 3 ps, it requires a pre-organized solvent–electrolyte interfacial structure to occur,24 which advocates that faster solvent reorganization kinetics would increase the HER rate. Altogether, these results highlight the importance of considering entropic effects when transition-state energies are calculated for electrochemical processes.
Fig. 9 (a) Snapshots from MD simulation showing water structure at Pt(100) and Pt(111) interfaces (reproduced from ref. 73 with permission from the National Academy of Sciences, Copyright (2013)) and (b) hydrogen-bond lifetime at the Pt(111) surface (blue) compared to bulk water (red) obtained by ab initio MD simulations (reproduced from ref. 75 with permission from the National Academy of Sciences, Copyright (2017)). |
Computational chemistry also proves useful to investigate ion-separation dynamics at the interface and more specifically the kinetics of the Grotthuss OH−/H3O+ shuttling mechanism. Using ab initio calculations for water–ions, while regular MD was used for cation–anion (Na+/I−) interactions, and considering the Madelung potential as the water collective contribution to the electric potential, Kattirtzi et al. showed that these two families of ions do not respond similarly at the electrode/electrolyte interface.75 Indeed, the activation energy for ion recombination at the interface is larger in the case of water–ions than for classical Na+/I− ions and the charge separation rate for classical ions is slower by more than one order of magnitude at the interface compared to the bulk, while water–ion separation remains equally fast. Interestingly, the activation energy for both systems is found to be about 1.5 kBT greater at the interface than in the bulk of the electrolyte. Nevertheless, in the case of water–ions, the flux of reactants is increased thanks to a larger lifetime for hydrogen bonds (Fig. 9b), which results in a neutral effect on the water–ion recombination kinetics.75 We believe that such a fine understanding of the water structure at the electrode–electrolyte interface would help to rationalize experimental observations from a molecular point of view. For instance, following the evolution of the hydrogen bond structure and dynamics upon polarization may provide a better understanding of the pzfc theory previously discussed. Moreover, these calculations so far limited to Pt and other model electrodes could be extended to other materials to provide new descriptors (H-bond lifetime at the catalyst surface, Grotthuss kinetics at the interface, etc.) for the HER.
Among the few studies carried out on ILs, most of them were focusing on the impact of moisture contamination on the cathodic stability of these ILs,78–80 with the aim to generate H2 through H2O reduction with low-cost catalysts81,82 or to construct a hydrogen-based reference electrode for ILs.83 Interestingly, following this strategy it was found that electrochemical interfaces are enriched in water upon polarization, following the hydrophilicity of the ions in solution (more water molecules at the positive electrode as IL anions are more hydrophilic).79 Furthermore, alike for aqueous electrolytes, the redox couple H+/H2 was found to be reversible on a Pt electrode in TFSI-based ILs.83 However, under these conditions, the Volmer step was identified as the rds which differs from the step usually found to be rate determining in acidic aqueous electrolytes.83 Similarly, using H2O as a proton donor, the activity of non-precious low carbon steel electrodes was found to outweigh that of Pt in BF4-based ILs.82 Altogether, these observations confirm the strong dependence of the proton donor structure on the HER mechanism and activity and certainly open further experimental explorations.
In contrast with ILs, a deeper understanding on the HER mechanism was developed in organic electrolytes. Interestingly, when using a strong acid such as HClO4 in acetonitrile, the proton reduction occurs at less negative potentials than water reduction on Au and Pt as well as remains independent of the water content,84,85 unlike for the HOR on Pt which is impacted by the water concentration.85,86 Rather, when water is used as the source of protons in an organic electrolyte, its reduction can be observed at more negative potentials (<−0.8 V vs. SHE on Pt(pc)) and strongly depends on the nature of the cation employed: the stronger the cation–H2O interactions (Li+ > Na+ for instance), the less negative the onset potential for water reduction is.84,85,87 It is no surprise that in the presence of large hydrophobic cations such as organic-ammonium, the water reduction can be nearly suppressed.84,87 Furthermore, in these organic–H2O electrolytes, the reduction of water follows on the different facets of Pt the same trend as in aqueous alkaline solutions, suggesting that the HER kinetics is controlled by the cleavage of the O–H in the Volmer step.87 Moreover, measuring the 1H NMR spectra of H2O in these electrolytes helped reveal the role of alkali cations in the proton acidification following the strength of the cation–water interaction, Li+ being the strongest.84 As previously discussed, in situ techniques such as sum frequency generation88 (SFG), Raman spectroscopy,89,90 FTIR spectroscopy,85,86 or surface X-ray diffraction91 (XRD) are currently employed to appreciate how these constraints on the water structure in the bulk of the electrolyte translate to the electrochemical interface. While all these techniques converge to the fact that polarizing the electrode negatively results in an enrichment of water at the interface, some discrepancies about the role of the cation,85,88 or water orientation88,91 remain.
Finally, working in organic electrolytes has enabled researchers not only to tune the nature and the amount of cation–water and water–water interactions, but also helped unravel the impact of the nature of the proton donor on the proton reduction (Fig. 10, left). Indeed, comparing two different proton donors possessing similar pKa values but different geometries, Jackson et al. demonstrated that the use of a bulkier proton donor slows down the HER in acetonitrile, while in an aqueous environment no such steric effect was observed for the HER, highlighting the critical role of water in the overall proton transfer,92 as later on rationalized by theoretical calculations.93 Similarly, adding crown-ethers known for their ability to chelate alkali cations to organic electrolytes has been shown to allow fine tuning of water–cation interactions and thus the water reactivity in organic electrolytes (Fig. 10, right).84 Interestingly, other electrocatalytic reactions such as the OER are already taking benefit of this approach to isolate reaction-intermediates that could not be observed in aqueous electrolytes due to their short lifetime.94,95 We thus believe that characterizing water-involving electrochemical reactions out-of-water could bring a wider comprehension on how water reacts. For instance, the role of the Grotthuss mechanism in the HER could be tested by comparing the reactivity of “single” water molecules isolated by working at low H2O content with molecules undergoing water–water interactions simply by increasing the H2O concentration in organic electrolytes (Fig. 10, bottom). Similarly, the influence of the proton donor structure (hydronium vs. water, or buffered electrolyte)33,42,49 could be studied (Fig. 10, top).
Fig. 10 Illustration of how water structure could be constrained by using an inert matrix (light green), and how the study of the water reactivity in all these cases could bring a fundamental understanding of several parameters in the HER kinetics. The top image is reproduced from ref. 42 with permission from Springer Nature, Copyright (2013). The image on the right is reproduced from ref. 84 with permission from the American Chemical Society, Copyright (2018). The bottom image is reproduced from ref. 63 with permission from John Wiley and Sons, Copyright (2006). The image on the left is reproduced from ref. 92 with permission from the American Chemical Society, Copyright (2016). |
In addition to these operando and theoretical studies, simple chemical strategies can be designed to gain deeper understandings at the molecular level.96 For that, decoupling all the different interactions at the interface (water–catalyst, water–water, ion–water, ion–catalyst) is critical. While challenging in aqueous solutions, the difficulty of decoupling these interactions can be alleviated by studying the HER under “non-usual” conditions such as in aprotic solvents or confined structures. This new comprehension could help understand the dynamic role of the water structure in the HER and pave the way towards the design of more efficient and stable HER catalysts. Moreover, not only this understanding is critical for the design of better water electrolyzers but, in light of recent studies on aqueous superconcentrated electrolytes for batteries10–12,97,98 or the role of water and proton impurities on the growth of solid–electrolyte interphases (SEI),13,99 this knowledge will certainly benefit the field of energy storage and conversion devices as a whole. Finally, achieving an in depth understanding of water structure and interfacial interactions is critical for other electrochemical reactions such as the OER,100 ORR,9 and CO2 (ref. 7) and N2 (ref. 8) reduction. In conclusion, the field is only at its infancy regarding the control at a molecular level of the water structure at the electrochemical interfaces, and no doubt that major achievements/understandings will be accomplished toward that goal in the near future.
This journal is © The Royal Society of Chemistry 2019 |