M.
Favaro
abc,
C.
Valero-Vidal
ad,
J.
Eichhorn
c,
F. M.
Toma
c,
P. N.
Ross
e,
J.
Yano
bf,
Z.
Liu
gh and
E. J.
Crumlin
*ad
aAdvanced Light Source, Lawrence Berkeley National Laboratory, One Cyclotron Rd, Berkeley, CA 94720, USA. E-mail: ejcrumlin@lbl.gov
bJoint Center for Artificial Photosynthesis, Lawrence Berkeley National Laboratory, One Cyclotron Rd, Berkeley, CA 94720, USA
cChemical Sciences Division, Lawrence Berkeley National Laboratory, One Cyclotron Rd, Berkeley, CA 94720, USA
dJoint Center for Energy Storage Research, Lawrence Berkeley National Laboratory, One Cyclotron Rd, Berkeley, CA 94720, USA
eMaterials Sciences Division, Lawrence Berkeley National Laboratory, One Cyclotron Rd, Berkeley, CA 94720, USA
fMolecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, One Cyclotron Rd, Berkeley, CA 94720, USA
gState Key Laboratory of Functional Materials for Informatics, Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences, Shanghai 200050, People's Republic of China
hDivision of Condensed Matter Physics and Photon Science, School of Physical Science and Technology, ShanghaiTech University, Shanghai 200031, China
First published on 7th March 2017
Understanding the interplay between surface chemistry, electronic structure, and reaction mechanism of the catalyst at the electrified solid/liquid interface will enable the design of more efficient materials systems for sustainable energy production. The substantial progress in operando characterization, particularly using synchrotron based X-ray spectroscopies, provides the unprecedented opportunity to uncover surface chemical and structural transformations under various (electro)chemical reaction environments. In this work, we study a polycrystalline platinum surface under oxygen evolution conditions in an alkaline electrolyte by means of ambient pressure X-ray photoelectron spectroscopy performed at the electrified solid/liquid interface. We elucidate previously inaccessible aspects of the surface chemistry and structure as a function of the applied potential, allowing us to propose a reaction mechanism for oxygen evolution on a platinum electrode in alkaline solutions.
Previous studies have been conducted to indentify the role of Pt oxides in the oxygen evolution reaction pathway (also known as oxygen discharge).8,24–26 Krasilshchikov24,26,28 demonstrated that at high overpotentials the corresponding current density (Tafel slope 120 mV dec−1) depends exponentially on the thickness of the Pt surface oxide. Other studies report much higher Tafel slopes (up to 160 mV dec−1)13,24,26 in both acidic and alkaline environments. Such an increase is often attributed to an increase of the surface oxide layer thickness and Pt dissolution.13,15–26,29 In alkaline media, the reversible binding of a hydroxide ion (OH−) coupled to a one electron oxidation is thought to precede the removal of one proton and one electron to form surface oxide species.30 This hypothesis has been mostly supported so far by surface structure and chemical characterizations far from operando conditions,6,31,32 thereby limiting an unambiguous investigation of the role of Pt oxides in the oxygen evolution catalytic cycle. Moving towards realistic operating conditions, Arrigo et al. recently studied Pt nanoparticles deposited on an acidic-polymer membrane using ambient pressure X-ray photoelectron spectroscopy (APXPS) in a humidified (solid/gas interface) environment.27 The authors found that the first stage of Pt oxidation involves the formation of a surface Pt oxide27,33 where atomic oxygen (O) is weakly bound to Pt (Pt–O).27,34
However, the detailed chemistry, structure, and dynamics of the surface oxide(s) on Pt at an electrified solid/liquid interface as a function of the applied potential and during oxygen evolution as a function of time still remains unclear. A definitive chemical speciation of such species and their role in the reaction pathway has not been completely achieved yet, in particular on Pt surfaces in alkaline conditions.
To achieve a molecular-level comprehension of the OER mechanism, we studied a polycrystalline Pt electrode surface under anodic conditions in an alkaline electrolyte, by means of operando APXPS performed with “tender” X-rays (hν = 4 keV) directly at the electrified solid/liquid interface.35,36 We report the speciation of porous Pt surface oxides as a function of applied potential, thereby elucidating the differences between thermodynamic predictions (Pourbaix diagrams)37 and the actual surface composition under polarization. Moreover, we study the Pt surface chemistry as a function of time during oxygen evolution, thereby observing the evolution of the different Pt oxidation states (as well as their spatial localization in the near-surface region) under non-equilibrium conditions (initial transient) and upon reaching steady state conditions. Based on these findings, we propose an OER mechanism in terms of a catalytic (Tolman) cycle38 for the various oxidized Pt species. We show that the rate-determining step is constituted by activated Pt sites selectively binding OH−. In addition, by coupling the experimental results with numerical simulations of the photoelectron intensity, we reveal that the active sites during OER are located in the sub-surface region, whereas Pt oxides such as Pt(II)O and Pt(IV)O2 are present in the high-surface outer layer and are not actively contributing to the oxygen evolution.
In order to study the evolution of the current density as a function of applied potential and time (by holding the applied potential at a constant value), we performed both cyclic voltammetry (CV) and chronoamperometry (CA) measurements (Fig. 1b), respectively. We explored different potentials along the anodic branch of the CV from open circuit potential (OCP) up to +900 mV (OER conditions, see Fig. 1b). In addition, the reversibility of the Pt surface chemistry and structure was assessed by returning back to OCP after OER conditions, defined as “OCP after OER”. The overpotential (η) has been computed using the standard Nernstian potential for water oxidation for a pH equal to 13.9 (+205 mV vs. Ag/AgCl/Cl(sat.)−).
To approximate the current density of the thin electrolyte layer under operando conditions formed by the “dip & pull” procedure, we conducted a separate electrochemical experiment that preserved the same “dip & pull” geometries, we masked the bottom of the electrode that was immersed in the bulk electrolyte and compared to an unmasked electrode (see schematization reported in ESI Fig. 1†). The thickness of electrolyte layer on the Pt WE was equal to 21.3 nm. This value was calculated via the inverse Beer–Lambert relation using the attenuation of the Pt 4f core level signal passing from hydrated conditions (saturated environment before “dip & pull” procedure) to the final solid/liquid interface configuration.31,36 We have determined that a bulk (unmasked) current density of about 1.0 mA cm−2 at +900 mV (overpotential η = +695 mV, OER conditions) corresponds to ca. 0.3 mA cm−2 in the nanometric-thick electrolyte layer on the WE surface (ESI Fig. 1b†). This is due to the predominant contribution of the current rising from the part of the WE immersed in the bulk electrolyte. This observation was reaffirmed by electrochemical impedance spectroscopy (EIS) measurements which extracted the polarization resistance (Rp) for the masked and unmasked configurations (see ESI Fig. 2 and 3, Tables 1 and 2 and Discussion 1 for further details†). The ratio between the two configurations polarization resistances (unmasked to masked) is equal to 3.37, thereby confirming that the total current density associated with the thin electrolyte film was about a third of the overall measured current.
Upon establishing current flow within the nanometric electrolyte layer, we can use the XPS data in Fig. 2 to confirm the applied potential directly at the probing location.35–39 O 1s and K 2p core levels as a function of the applied potential are reported in Fig. 2a–c, respectively. In particular, Fig. 2b reports a magnification of the O 1s multipeak fitting procedure for the low BE tail. Notably, at different applied potentials, the area of the OHsol− (hydroxide ions in solution) and K 2p (solely depicting the K+ ions in solution) spectral feature remains constant throughout the experiment. Indeed, the liquid phase water (LPW)/OHsol− and LPW/K 2p concentration ratio is equal to 51.4 ± 2.1 as the applied potential is changed from OCP to OER conditions (Fig. 2d). Our measurement is therefore in agreement with the expected value for a 1.0 M electrolyte concentration, which is 55.51 water to one KOH.
The binding energy (BE) trends of the O 1s LPW, OHsol− and K 2p3/2 core levels present in the liquid layer as a function of the applied potential is reported in Fig. 2d. The BEs for elements in the electrolyte (BEelectrolyte) obey the relation:
BEelectrolyte = BE0 − eE | (1) |
E is the applied potential and BE0 is the BE measured in absence of an applied potential. According to eqn (1), the aforementioned core levels shift toward lower BEs with respect to the Fermi level of the sample (grounded to the analyzer, see schematization reported in Fig. 1a) for applied anodic potentials, thereby confirming that potential control is achieved in the measuring spot (1:1 shift).
From the O 1s low BE tail in Fig. 2b, it is possible to observe that the area of the Pt–O component (centered at 530.4 eV) increases as well as the adsorbed OH− (OHabs− at 531.3 eV).
Finally, the area of the chemisorbed water component (H2Ochem centered at 532.4 eV) also remains constant as the potential changes from OCP to +500 mV, with an increase under OER at +900 mV. Interestingly, upon removing the potential (“OCP after OER”), the peak area does not appreciably change. Such observation is consistent with an increase of the electrode surface area under catalytic conditions, which will be detailed in a following section.
As proposed by Krasilshchikov28 and Damjanovic,24,26 it is well-accepted that under alkaline conditions the OER on Pt starts with the nucleophilic attack of OH− ion on Pt followed by a fast electron transfer to the metastable configuration where OH− is then chemisorbed on an “activated” platinum site (Ptδ).24,26 The formal positive charge of the Ptδ (set between Pt metal and Pt(II)) in the Ptδ–OHads active site, as suggested by Wagner et al.,42 could explain the spectral position of its corresponding photoelectron peak, centered between the metal and Pt(II) spectral components.
Interestingly, while the integrated area of this component remains unchanged with applied potentials below the OER threshold, it drastically increases under OER conditions. This finding is supported by a similar trend observed in the O 1s spectrum for the OHabs− (Fig. 2a and b). At OCP after OER, the intensity of this component is comparable to the initial OCP state. This reversibility suggests that the formation of this component is driven by high overpotentials (i.e. under OER conditions). Since oxygen evolution in alkaline environments is proceeded by the deprotonation of adsorbed OH−, we conclude that the observed spectral feature, Ptδ–OHads, is representing the active species in the OER catalytic cycle.
Together with Ptδ–OHads, the amount of Pt(II) species increases with the applied potentials, where we observe a sharp increase once the potential is raised to +900 mV (OER conditions). Under this applied potential, the formation of Pt(IV)O2 oxide becomes prominent concurrently with the previously discussed Pt species (Fig. 3).
The evolution of the outer surface morphology as a function of the applied potential was studied by ex situ AFM (Fig. 3c). The surface undergoes a clear reconstruction upon OER conditions transitioning from a flat metallic surface at OCP and low applied potentials (as confirmed by the low surface roughness (SR) values) to a high surface area oxidic system after catalytic conditions, characterized by a much higher SR then the pristine surface. Therefore, we believe that the so-called “gelly structure” reported in literature for Pt surfaces under OER conditions in acidic electrolytes,27,43,44 also occurs when the reaction is driven in alkaline environments. Using the previously described simulations, we can investigate the chemistry of the mesoporous (“gelly”) structure, which is composed of a Pt(IV)O2 outer layer and an inner layer (on top of the Ptδ–OHads active phase) constituted by a mixture of Pt(II) phases consisting of Pt(II)(OH)2 and Pt(II)O.
On the other hand, the experimental results clearly shows that up to +500 mV Ptδ–OHads and Pt(II) species are present on the surface (Fig. 3). Notably, the Pt(IV) spectral component is eventually observed only under the OER potential (+900 mV), therefore at a much higher potentials than the thermodynamically predicted in the Pourbaix diagram at the same pH. Such a deviation can be due to kinetic effects that are not taken into account in a merely thermodynamic investigation of the stability of Pt phases. In addition, at OCP after OER, the surface chemistry exhibits a partial chemical hysteresis. While the Ptδ–OHads feature decreases (as a result of the potential removal, as previously discussed), it is still possible to observe both Pt(II) and Pt(IV) spectral fingerprints. On the contrary, the Pourbaix diagram predicts only Pt(II) (present as Pt(II)(OH)2) as the stable phase at this pH value.
It is important to consider the possibility for Pt(IV) dissolution in the form of a Pt(IV) aquocomplex as indicated by the Pourbaix diagram. Based on our data, we do not observe any Pt dissolution during the potentiostatic measurements. However, one must consider that reversing (or removing) potentials could induce Pt dissolution on the order of tens of nM as reported by Lopes et al.15 The reported dissolution rates are considerably below the detection limit of the APXPS (about 100 mM). Therefore, although we would not be able to detect this process, a Pt dissolution cannot be excluded.50,51
This observation allows us to propose that the applied overpotential (above the thermodynamic oxidation potential of water) has two main effects. First, it positively polarizes the electrode surface, which lowers the binding energy of hydroxyl ions and increase their adsorption rate. Second, the OH− may electrostatically be driven through the oxidic porous passivation layer.52 Therefore, from the experimental observations reported here we provide motivation for reevaluating the kinetic limitations (overpotentials) of OER on Pt to take into account both adsorption/electron transfer rate and mass transport effects through the modified electrode under catalytic conditions.
Moreover, we have elucidated the differences between thermodynamic predictions of the Pt phase stability and the actual surface chemistry under polarization. The high chemical sensitivity of XPS probing the solid/liquid electrified interface has allowed us to re-interpret the Pourbaix diagram of Pt. Thus motivating the need for these types of operando experiments and the development of new thermodynamic, kinetic, and computational models for electrified interfaces under realistic operating conditions.53 Expanding this experimental strategy for the study of different materials and catalytic reactions under operando conditions,31,32,54 will open a path for new discoveries that will deepen our knowledge towards the accelerated development of novel tailored materials.
Prior its introduction into the experimental chamber, the electrolyte (KOH 1.0 M aqueous solution) was outgassed for at least 30 min at low pressure (around 10 Torr) in a dedicated off-line chamber. Then, once the manipulator (Fig. 1a) and the outgassed electrolyte were placed into the APXPS experimental chamber, the pressure was carefully lowered down to the water vapor pressure (between 16 and 20 Torr, at r.t.).
To create a solid/liquid interface (WE/electrolyte), all three electrodes were immersed into the electrolyte. Then they were slowly extracted from the electrolyte solution by rising the manipulator at a constant vertical rate. Following this procedure, a thin layer of liquid electrolyte film is formed on the Pt electrode (typical values ranges between 10 and 30 nm); the three electrodes were then positioned at the focal point of the analyzer, giving the possibility of performing XPS investigation of the solid/liquid interface as a function of the applied potential. During the measurements, the bottom parts of the electrodes were kept into the bulk electrolyte, in order to ensure electrical continuity between the thin electrolyte layer on the WE surface and the bulk electrolyte itself.
The working (WE) and counter electrodes (CE, Pt polycrystalline foils, 99.99%, thickness 0.5 mm, Aldrich) were polished to a mirror finish with silicon carbide papers of decreasing grain size (Struers, grit: 2400 and 4000). The samples were then cleaned by two cycles of ultrasonic treatment in a mixture of MilliQ water/ethanol (Aldrich, 1:1) for 10 min. A third ultrasonic cleaning was then conducted in pure MilliQ water for 15 min, followed by a thoroughly rinsing and drying in N2 stream. As reference electrode (RE), a miniaturized leakless Ag/AgCl/Cl(sat.)− RE (ET072-1, eDAQ) was used (standard electrode potential with respect to the normal hydrogen electrode, NHE). All the potentials reported in this work are referred to this RE unless differently stated.
An electrochemical cleaning procedure was conducted by holding the Pt WE at −1200 mV (within the HER region) for 30 minutes, in order to obtain a pure metallic surface and a homogeneous surface, increasing in this way the experimental reproducibility. The cleaning procedure was stopped after reaching the current open circuit potential (OCP) of the cell from the cathodic side. The data acquisition started only when the OCP reached a stable value over the observation time (about 30 minutes). All the voltammetric measurements were conducted at a scan rate of 20 mV s−1.
Electrochemical impedance spectroscopy (EIS) measurements were carried out with the aim of electrochemically evaluate the interfacial properties of the solid/liquid interface. EIS measurements were conducted by measuring the response in current after applying to the WE a sinusoidal wave of potential of 10 mV peak-to-peak voltage amplitude, at OCP (−200 mV) and OER conditions (+900 mV). The frequency of the sinusoidal potential wave was ranging from 105 Hz to 5 mHz, scanning 10 frequencies per decade.
For the data acquisition, the pass energy of the Scienta analyzer (R4000 HiPP-2) was set to 200 eV, using a step of 100 meV and a dwell time of 300 ms. Under these conditions, the total resolution (X-rays and analyzer) was equal to about 250 meV at room temperature (r.t.) and at 4 keV. The measurements were taken using a photon energy of 4 keV at r.t. and in normal emission (NE), at a pressure in the experimental chamber matching the water vapor tension at r.t., between 16 and 18 Torr. The detection stage in the analyzer (multichannel plate, MCP) was under high vacuum conditions (∼10−7 Torr) throughout the all experiment, thanks to the differentially pumped electrostatic lens system based on four independent pumping stages. To limit the evaporation from the electrochemical cell we place a larger outgassed solution buffer in the analysis chamber to facilitate the fully humidified environment at room temperature. The calibration of the binding energy (BE) scale was carried out using the Au 4f photoelectron peak as reference (4f7/2 BE = 84.0 eV) as well as the Fermi edge, from a clean gold polycrystalline surface measured under all the investigated experimental conditions.
All the fit reported in this work have been carried out using a Doniach–Šunjić shape for the Au and Pt 4f photoelectron peaks, whereas a symmetrical Voigt function (G/L ratio ranging from 85/15 to 75/25) was used to fit the O 1s photoelectron peaks (after Shirley background subtraction).36 During the fitting procedure, the Shirley background was optimized as well together with the photoelectron intensity of the spectral components, increasing in this manner the precision and reliability of the fitting procedure.36,55,56 The χ2 minimization was ensured by the use of a nonlinear least squares routine, with increased stability over simplex minimization.36,55,56
Data processing of the raw images was performed using customized Igor procedures.
From observations already reported in literature about the oxidic surface structure of Pt under OER conditions,27,44 the simulations were built adopting a conformal layer-by-layer structure “buried” underneath an electrolyte (KOH 1.0 M) film characterized by an average thickness of 21.3 ± 0.9 nm. Three crucial aspects needs to be highlighted:
- The simulations were performed using net and well-defined layer/layer interfaces. This is an approximation of the real configuration, where most likely gradients and mixed regions may exist at the layer/layer interfaces. However, the spatial extension of such regions are typically of about 1–2 nm, a dimension that is typically difficult to access and discriminate with tender X-rays. Therefore, the approximation used here is reasonable at the light of the features offered by the APXPS technique;
- The simulations were conducted by using the well-known growth modes: layer-by-layer or Frank-Van der Merwe, layer-island or Stranski–Krastanov and island formation or Volmer–Weber.61,62 Moreover, for a chosen model, we have investigated different spatial localization of the oxidized Pt chemistries (see summary reported in ESI Fig. 4†). The model that fitted better with the experimental simulations was found by using a chi-square (χ2) statistical test. The quantities used for its evaluation were the Ptδ–OHads/Pt0 (R1), Pt (II)/Pt0 (R2) and Pt(IV)O2/Pt0 (R3) ratios, determined experimentally through the multipeak fitting procedure of Pt 4f photoelectron peak (EXP) and from the numerical simulations (SIM). Therefore:
For the final evaluation of model/configuration, the space of parameters was explored using a variational principle: for each model, the minimum value of the χ2 was achieved by properly adjusting the layer thickness of the different oxidized Pt species. This procedure allows to determine both the structure and the corresponding configuration (layer thicknesses) that best match with the experimental results;
- The relative surface area (RSA) for the top most layer in the calculations have been derived directly from the experimental ex situ AFM surface characterizations reported in Fig. 3.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7ta00409e |
This journal is © The Royal Society of Chemistry 2017 |