Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Toward a quantitative theoretical method for infrared and Raman spectroscopic studies on single-crystal electrode/liquid interfaces

Yuan Fang a, Jin-Chao Dong a, Song-Yuan Ding *a, Jun Cheng a, Juan Miguel Feliu b, Jian-Feng Li a and Zhong-Qun Tian a
aState Key Laboratory of Physical Chemistry of Solid Surfaces (PCOSS), Collaborative Innovation Center of Chemistry for Energy Materials (iChEM), Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China. E-mail:
bInstituto Universitario de Electroquímica, Universidad de Alicante, Carretera San Vicente del Raspeig s/n, E-03690 San Vicente del Raspeig, Alicante, Spain

Received 28th October 2019 , Accepted 9th December 2019

First published on 10th December 2019


In situ electrochemical infrared spectroscopy and Raman spectroscopy are powerful tools for probing potential-dependent adstructures at solid/liquid electrochemical interfaces. However, it is very difficult to quantitatively interpret the observed spectral features including potential-dependent vibrational frequency and spectral intensity, even from model systems such as single-crystal electrode/liquid interfaces. The clear understanding of electrochemical vibrational spectra has remained as a fundamental issue for four decades. Here, we have developed a method to combine computational vibrational spectroscopy tools with interfacial electrochemical models to accurately calculate the infrared and Raman spectra. We found that the solvation model and high precision level in the self-consistent-field convergence are critical elements to realize quantitative spectral predictions. This method's predictive power is verified by analysis of a classic spectroelectrochemical system, saturated CO molecules electro-adsorbed on a Pt(111) electrode. We expect that this method will pave the way to precisely reveal the physicochemical mechanism in some electrochemical processes such as electrocatalytic reactions.


The determination of adstructures at electrochemical (EC) solid/liquid interfaces is a fundamental issue in fuel cells, metal/alloy plating and corrosion, etc.1–3 Vibrational spectroscopies can be utilized to provide fingerprint information about adstructures with high spectral resolution and have been developed for the characterization of EC interfacial adstructures by infrared (IR) spectroscopy since the mid-1960s4 and by Raman spectroscopy since the mid-1970s.5 Nevertheless, the vibrational frequencies and the intensities of EC-IR and EC-Raman spectra strongly depend on the applied potential, electrode materials, coverage of adsorbates and coadsorbed species, and thus are too complicated to be clearly interpreted in most of cases.6,7 For instance, with the Stark tuning slope (STS), the slope of the vibrational frequency as a function of the applied potential, it is difficult to precisely quantify potential-dependent behaviours of electroadsorption. Even for single-crystal electrodes with structurally well-defined surfaces, as perfect model systems for studying electroadsorption and electrocatalytic reactions,8–14 researchers still faced challenges in unambiguously assigning the vibrational modes in EC-IR and/or EC-Raman spectra. The above issue exposes the need for developing comprehensive and highly precise computational tools that facilitate the interpretation of EC-vibrational spectra.

First-principles computational methods based on cluster models with metal clusters to mimic the electrodes have been employed for calculating potential-dependent vibrational spectra of electroadsorbates.15–18 However, the cluster models in EC simulations are usually too simplified to exactly consider important effects on electrified single-crystal electrodes, such as the intermolecular interactions between electroadsorbates and periodic lattice structures. Therefore, slab models with periodic structures have been widely employed to study electroadsorption configurations and electrocatalytic mechanisms, especially for single-crystal or nanocrystal electrodes.19–29

At present, many studies have calculated the vibrational frequencies of adsorbates on slab models without the solvation model, which leads to the calculated STSs being much smaller than the experimental values.30–33 Moreover, studies have rarely calculated EC-IR or EC-Raman intensities,34 although crucial information, such as the orientation and coverage of adsorbates, adsorbate–adsorbate interactions, and charge-transfer interactions between adsorbates and substrates, could be extracted from the experimental analysis of the spectral intensities.35,36 Some efforts to understand the intensities were made based on the analysis of the symmetry of groups and the surface selection rule. These approaches can determine whether the vibrational modes of adstructures are either IR- or Raman-active, or both, but fail to predict potential-dependent relative and absolute intensities. In fact, accurately calculating EC-IR and EC-Raman intensities based on slab models by simultaneously considering the electrified surfaces and the surface solvation effect is a longstanding and difficult task.

In this article, we report a new method for quantitatively predicting not only the vibrational frequencies but also the intensities of EC-IR and EC-Raman spectra from single-crystal electrodes. By coupling the surface charged method to mimic electrified surfaces (Scheme 1a),22,37,38 the implicit solvation model (Scheme 1b),28,29 and the finite difference method, the calculation of vibrational frequencies and intensities is successfully achieved (Scheme 1c; see also the ESI for details).

image file: c9sc05429d-s1.tif
Scheme 1 The computational method for EC-IR and EC-Raman spectra: (a) surface charged method, (b) implicit solvation model and (c) finite difference method. μ, F, X, A, υk and GZ represent the surface dipole moment, atomic force, external atomic coordinates, dynamical matrix at the Γ point, vibrational frequency of the kth normal-mode coordinate Qk, and external electric field along the Z axis, respectively. Δx and ΔGZ represent the step size of X and GZ for the finite difference. Image (a) adapted from ref. 38 with permission, copyright (2011) Royal Society of Chemistry; image (b) adapted from ref. 29 with permission, copyright (2014) American Institute of Physics.

Results and discussion

Simulation of EC-IR spectra of the Pt(111)(2 × 2)-3CO adstructure

To validate the present computational method, we investigated the classic EC system of carbon monoxide (CO) adsorbed onto Pt(111) electrodes.8,39–46 Villegas and Weaver observed the Pt(111)(2 × 2)-3CO adstructure (Fig. 1e) with three CO molecules adsorbed onto a p(2 × 2) Pt(111) in a unit cell by performing EC-STM at potentials below 0.44 V vs. the standard hydrogen electrode (SHE) in a CO-saturated 0.1 M HClO4 aqueous solution.44 Two possible structures were proposed for the adstructure. One was Pt(111)(2 × 2)-3CO α1 with one CO molecule adsorbed at an atop site (COL) and two CO molecules at hollow sites (COM) in the unit cell (Fig. 1a). The other one was Pt(111)(2 × 2)-3CO α2 with a CO molecule adsorbed at a bridge site (COB) and two CO molecules at near-top sites (Fig. 1b). Additionally, in experimental EC-IR spectra, the authors found a strong band from 2066 to 2074 cm−1 and a relatively weak band from 1780 to 1798 cm−1 (Fig. 1h). The two bands were assigned to the COL and COM stretching modes, respectively. Finally, they concluded that the most likely adstructure was the α1 rather than the α2 adstructure.47,48
image file: c9sc05429d-f1.tif
Fig. 1 Pt(111)(2 × 2)-3CO adstructures at 0.04–0.44 V vs. SHE in a CO-saturated 0.1 M HClO4 aqueous solution: (a and b) the two possible adstructures (2 × 2)-3CO α1 and α2, respectively, proposed by Villegas et al.; (c and d) the respective optimized structures; (e) the measured EC-STM pattern;44 (f and g) the corresponding calculated patterns of the α1 and α2 structures; (h) the measured EC-IR spectra;47 and (i and j) the corresponding calculated spectra of the α1 and α2 adstructures. Images (a), (b) and (e) adapted from ref. 44 with permission, copyright (1994) American Institute of Physics; image (h) adapted from ref. 47 with permission, copyright (2000) American Chemical Society.

In our calculated EC-IR spectra of the α1 structure (Fig. 1i), the bands centred at 2065–2074 cm−1 and 1789–1798 cm−1 are assigned to the stretching mode of one COL molecule and the in-phase stretching mode of two COM molecules, respectively. The localized or combinational CO stretching modes are shown in Fig. S1. The calculated IR relative intensity I(COL)/I(COM) is ca. 1.2, which is close to the measured value of 1.5. Hence, the calculated vibrational frequencies and intensities of the α1 structure suitably agree with the values in the experimental EC-IR spectra. In addition, the calculated STM image in Fig. 1f contains one line of the brightest spots, one line of the second-brightest spots and one line of the darker spots along the A1 direction, which is consistent with the experimental image (Fig. 1e). The α2 adstructure was unstable in the potential range from 0.04 to 0.44 V vs. SHE. Moreover, the frequency of the COB stretching mode (ca. 1827 cm−1), the IR relative intensity I(COL)/I(COB) (4.1) and the STM pattern of the α2 adstructure (Fig. 1g and j) disagreed with the experimental observations. Therefore, this α2 structure was discarded.

Importance of the solvation model for quantitatively predicting STSs of potential-dependent vibrational frequencies

Consideration of the solvation effect is a fundamental requirement for the reliable prediction of STSs. The STSs calculated from the COL and COM bands involving the implicit solvation model are 21 cm−1 V−1 and 24 cm−1 V−1, respectively. However, the calculated STSs of the two bands without the solvation model equal only ca. 1 cm−1 V−1 (Fig. S3a), and compared with the measured slopes (22 cm−1 V−1 and 43 cm−1 V−1), these values are remarkably underestimated. We found that solvation could significantly increase the effective electrostatic field across the interfaces, further polarize the adsorbates, and finally result in much larger STSs than those in a vacuum. This role of solvation effect on STSs could also be elucidated by the monotonous increase in the absolute values of STSs as the increase in the relative permittivity εre of solvents as shown in Fig. S4b. Furthermore, compared with the measured STS of the relatively weak band (43 cm−1 V−1), the calculated STS of the COM band (24 cm−1 V−1) is underestimated. This finding appears analogous to the “CO/Pt(111) puzzle”,49,50 which refers to the underestimated energy of the 2π* orbital from COM by DFT methods, and to the overestimation of the interaction between COM molecules and Pt substrates in a vacuum environment. Moreover, we found that the overestimated interaction was unaffected by the charged surface and the implicit solvation, which might in turn cause the underestimation of the STS of the COM band in the EC environment.

Simulation of EC-Raman spectra of the Pt(111)(2 × 2)-3CO adstructure

EC-Raman spectroscopy is suitable for the characterization of vibrational modes occurring at low spectral frequencies. These vibrations directly reflect chemical complexation between adsorbed molecules and transition metal substrates. Therefore, we performed shell-isolated nanoparticle-enhanced Raman spectroscopy (SHINERS) on CO/Pt(111) EC interfaces in a CO-saturated 0.1 M HClO4 aqueous solution (Fig. 2a, experimental details in the Experimental section). The calculated vibrational frequencies and relative intensities of the Raman spectra of the α1 adstructure in Fig. 2b consistently agreed with the ones found by SHINERS. The Raman bands centred at ca. 2071, 1793, 473 and 394 cm−1 are assigned to the COL, COM, PtCL and PtCM stretching modes (Fig. S1), respectively.
image file: c9sc05429d-f2.tif
Fig. 2 (a) Experimental EC-Raman spectra of a Pt(111) single-crystal electrode measured by SHINERS at 0.14–0.44 V vs. SHE in a CO-saturated 0.1 M HClO4 aqueous solution. (b) Calculated EC-Raman spectra of the (2 × 2)-3CO α1 adstructure. Potential-dependent NBO occupancies of the (c) 1π and 2π* orbitals of the COL and COM bonds and the (d) σ and σ* orbitals of the PtCL bond. The slopes of the NBO occupancies are listed in Table S1. The detailed NBO wavefunctions are listed in Fig. S6 and S7.

It is interesting to note that the STSs of the PtCL and PtCM bands are negative (−8 and −6 cm−1 V−1, respectively), while the COL and COM counterparts are positive, indicating that the CO bond strengthens and the PtC bond weakens as the potential increases. To understand the molecular mechanism behind this trend, we analysed the electronic structure of the α1 adstructure. The projected electronic densities of states (pDOSs) in Fig. S5a and b show that the hybridized 5σ orbital of the CO molecule donates electrons to the Pt(111) surface to strengthen the PtC bonds, while the hybridized CO 2π* orbital accepts the back-donated electrons from the Pt d orbitals to weaken the CO bonds. Moreover, the natural bond orbitals (NBOs) in the periodic implementation were computed to quantitatively analyse the chemical-bonding response to the applied potential.51Fig. 2c and d show that the NBO occupancies of the CO 2π* and PtCL σ orbitals shift negatively as the potential increases, which induces the strengthening of the COL and COM bonds and the weakening of the PtCL bonds. In addition, the variations in the NBO occupancies of the 5σ, 5σ*, and 1π orbitals of CO and the PtCL σ* orbital are insensitive to the applied potential (Table S1).

Importance of high-level precision in the SCF convergence for quantitatively predicting vibrational intensities

To precisely predict the vibrational intensities, the precision level of the convergence in the self-consistent-field (SCF) calculation (εSCF) should be as high as 1 × 10−9 eV. For example, according to symmetry analysis, modes (ii) and (iii), assigned to the in-phase and out-of-phase combinational stretching modes of the two COM molecules, respectively, should originally be IR-active and IR-inactive. However, the two modes become pronouncedly IR-inactive and IR-active, respectively, if εSCF is 1 × 10−6 eV (Fig. 3a) which is the value by default in the calculation of vibrational frequencies. Similarly, several misleading relative intensities of Raman bands were predicted with larger εSCF values (Fig. 3b). This problem might be the reason for the scarce publications reporting the IR or Raman intensities at EC interfaces. In fact, all the spectra in Fig. 1 and 2 were calculated while including the implicit solvation model and an εSCF of 1 × 10−9 eV. With the as-developed method in hand, we can start predicting the EC-IR and EC-Raman spectra of some unknown adstructures at EC interfaces.
image file: c9sc05429d-f3.tif
Fig. 3 (a) Calculated EC-IR and (b) EC-Raman spectra of the (2 × 2)-3CO α1 adstructure at 0.45 V vs. SHE in a CO-saturated 0.1 M HClO4 aqueous solution under different SCF convergence criteria from 10−6 to 10−10 eV. The blue lines are the measured EC-IR in (a) and EC-Raman spectra in (b) at 0.44 V vs. SHE. The symbols (i), (ii), (iii), (iv), (v) and (vi) correspond to the COL stretching mode, symmetric COM stretching mode, antisymmetric COM stretching mode, PtCL stretching mode, CO rotation mode and PtCM stretching mode, respectively.

Prediction of the EC-IR and EC-Raman spectra of metastable CO/Pt(111) adstructures

Other than the well-known (2 × 2)-3CO α1 adstructure, Jung et al. proposed two unitarily atop-coordinated adstructures, (2 × 2)-3CO β and (1 × 1)-CO, at 0.36–0.45 V vs. SHE (Fig. 4a and c). However, neither the featured COL stretching bands (∼2094 cm−1 and ∼2110 cm−1) nor the PtCL stretching band (482 cm−1) of the two adstructures in Fig. 4g–j can be observed in the measured IR spectra in Fig. 1h and the measured Raman spectra in Fig. 2a (∼2074 cm−1 and ∼474 cm−1).
image file: c9sc05429d-f4.tif
Fig. 4 Two metastable adstructures at 0.36–0.45 V vs. SHE in a CO-saturated 0.1 M HClO4 aqueous solution: (a and c) the measured EC-STM images of the (2 × 2)-3CO β and (1 × 1)-CO adstructures, respectively;45,46 (b and d) the corresponding calculated EC-STM images; (e and f) the corresponding optimized adstructures; (g and h) the corresponding calculated EC-IR spectra; (i and j) the corresponding calculated EC-Raman spectra. Image (a) adapted from ref. 45 with permission, copyright (2006) Royal Society of Chemistry; image (c) adapted from ref. 46 with permission, copyright (2007) American Chemical Society.

Accordingly, the absolute values of the binding energies (|EB|) of the (2 × 2)-3CO β and the (1 × 1)-CO adstructures are approximately 0.39 and 0.59 eV smaller, respectively, than those of the (2 × 2)-3CO α1 adstructure (Fig. S8), which indicates that the (2 × 2)-3CO β and (1 × 1)-CO adstructures are metastable and might occupy domains of the Pt(111) electrode down to the nanometer scale. Thus, the (2 × 2)-3CO α1 adstructure dominantly contribute to the EC-IR and EC-Raman intensities by wide-field spectroscopic measurement, while the two metastable adstructures might provide only a negligible contribution to the wide-field spectra and the vibrational bands assigned to the metastable adstructures cannot be observed. It would be desirable to characterize the metastable adstructures in the minute domains by employing tip-enhanced Raman spectroscopy (TERS) at solid/liquid EC interfaces with ultrahigh sensitivity and spatial resolution or by developing a method for nanoscale IR spectroscopy in liquid.


A theoretical method combining the computational tools of vibrational spectroscopy and interfacial-EC models (surface charge method and implicit solvation model) was developed for quantitatively predicting EC-IR and EC-Raman spectra, not only the vibrational frequencies but also the spectral intensities. The validity of the combined method was well demonstrated by the study of the well-known Pt(111)(2 × 2)-3CO α1 adstructure. We demonstrated that the precision level of the convergence in the SCF calculation should be as high as 1 × 10−9 eV rather than the commonly used 1 × 10−6 eV for the quantitative prediction of the relative intensities of EC-IR and EC-Raman spectra. The as-developed method can be straightforwardly applied to reveal the electro-adstructures of general molecules on general single-crystal electrodes, such as bimetallic electrodes, semiconductor electrodes, and some emerging 2D material electrodes. In addition, the implicit solvation model may be further combined with explicit water molecules in the inner Helmholtz plane for more complicated reactant-water or oxyanion-water adstructures coadsorbed on electrified surfaces. This work may provide a new opportunity to deeply reveal reaction mechanisms mediated by coadsorbed water molecules and to study the potential-induced phase transition of a structured-water layer in future.

Computational section

The metal surfaces were modelled as p(2 × 2) and p(√19 × √19) Pt(111) slabs with a lattice constant of 3.97 Å and a thickness of 7 layers in a periodic box of 108.75 Å (the middle three layers were frozen in its bulk position and detailed geometries are given in Appendix I of the ESI). The first-principles computations were performed using the Vienna ab initio simulation package (VASP) with projector augmented wave (PAW) pseudopotentials.52,53 The alternative revision of the Perdew–Burke–Ernzerhof functional RPBE was employed to exactly calculate the chemisorption energies of CO molecules on the Pt(111) surfaces.54 We considered a 500 eV cut-off energy for reciprocal space mesh size and a 5 × 5 × 1 and 2 × 2 × 1 Γ-centred K-point mesh for p(2 × 2) and p(√19 × √19) unit cells, respectively. The SCF convergence criterion was set to 1 × 10−9 eV. The force of optimized geometries should be less than 0.01 eV Å−1 for converged vibrational frequencies of the systems (Table S2).

The implicit solvation model implemented in VASPsol based on the linearized Poisson–Boltzmann equation (also named the linearized polarizable continuum model, LinearPCM) was employed to mimic the surface solvation effect (see the ESI for details).28 The relative permittivity εre and Debye screening length λ were, respectively, set to 78 and 9.5 Å for ionic strength I = 0.1 M to describe the 0.1 M HClO4 aqueous solution in experiments. The calculated potential of zero charge (PZC) in this article under the adstructure √19 × √19 (a 0.68 monolayer) for the coverage (CO per surface Pt atoms), is 1.22 V vs. SHE and is close to the measured PZC, 1.1 ± 0.04 V.55,56 The precise prediction of PZC suggests that the implicit solvation model is reasonable for mimicking the CO/Pt(111) interface at a high CO coverage. In addition, more accurate implicit solvation models such as nonlinear electrochemical soft-sphere model should be employed for more accurate description of solid/liquid interfaces with low ionic strengths,57,58 because the LinearPCM model might fail to simulate the ionic response in low ionic strengths and the capacitance of solvent near PZC.

Experimental section

We employed an Xplora Raman instrument to perform in situ spectroelectrochemical experiments. The excitation laser wavelength of the Xplora instrument was 638 nm (the laser power about 6 mW) and a 50× magnification long working distance (8 mm) objective was used. Before the Raman experiment, the Raman frequencies were calibrated using a Si wafer and the experimental Raman spectra of the experiment were collected during 60 s for a single-spectrum curve. A three-electrode spectroelectrochemical cell was employed for the electrochemical Raman experiment. A platinum wire was employed as a counter electrode and a saturated calomel electrode (SCE) was used as a reference electrode. A Pt(111) single-crystal electrode was used as a working electrode assembled with Au@SiO2 nanoparticles to enhance the Raman signal.11,59 The CO electro-adsorption experiment was performed in a CO-saturated 0.1 M HClO4 aqueous solution.

Conflicts of interest

There are no conflicts to declare.


The authors acknowledge funding support from the National Natural Science Foundation of China (21403179, 21727807, 21872115, 21533006, and 21621091) and MOST of China (2015CB932300). The first-principles calculations in this work were performed on the supercomputing system in the big data system and computing platform of the State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University. Y. F. and S. Y. D. sincerely thank Prof. Philippe Sautet and Dr Stephan N. Steinmann for discussion on the surface charge method, Prof. Georg Kresse for discussion on the simulation of vibrational spectra, Prof. Jinlong Yang and Prof. Xin Xu for discussion on the “CO puzzle” issue, and Prof. Jun Li for discussion on the NBO analysis in the periodic implementation.

Notes and references

  1. V. R. Stamenkovic, D. Strmcnik, P. P. Lopes and N. M. Markovic, Nat. Mater., 2017, 16, 57–69 CrossRef CAS PubMed.
  2. H. Zhang, H. Ning, J. Busbee, Z. Shen, C. Kiggins, Y. Hua, J. Eaves, J. Davis III, T. Shi, Y. T. Shao, J. M. Zuo, X. Hong, Y. Chan, S. Wang, P. Wang, P. Sun, S. Xu, J. Liu and P. V. Braun, Sci. Adv., 2017, 3, e1602427 CrossRef PubMed.
  3. H. T. Dinh, J. Kuever, M. Mussmann, A. W. Hassel, M. Stratmann and F. Widdel, Nature, 2004, 427, 829–832 CrossRef CAS PubMed.
  4. H. B. Mark and B. S. Pons, Anal. Chem., 1966, 38, 119–121 CrossRef CAS.
  5. M. Fleischmann, P. J. Hendra and A. J. McQuillan, Chem. Phys. Lett., 1974, 26, 163–166 CrossRef CAS.
  6. M. J. Weaver and S. A. Wasileski, Langmuir, 2001, 17, 3039–3043 CrossRef CAS.
  7. S. A. Wasileski and M. J. Weaver, Faraday Discuss., 2002, 121, 285–300 RSC.
  8. S. G. Sun, J. Clavilier and A. Bewick, J. Electroanal. Chem. Interfacial Electrochem., 1988, 240, 147–159 CrossRef CAS.
  9. H. Ogasawara and M. Ito, Chem. Phys. Lett., 1994, 221, 213–218 CrossRef CAS.
  10. K. I. Ataka, T. Yotsuyanagi and M. Osawa, J. Phys. Chem., 1996, 100, 10664–10672 CrossRef CAS.
  11. J. F. Li, S. Y. Ding, Z. L. Yang, M. L. Bai, J. R. Anema, X. Wang, A. Wang, D. Y. Wu, B. Ren, S. M. Hou, T. Wandlowski and Z. Q. Tian, J. Am. Chem. Soc., 2011, 133, 15922–15925 CrossRef CAS PubMed.
  12. N. R. Honesty and A. A. Gewirth, J. Raman Spectrosc., 2012, 43, 46–50 CrossRef CAS.
  13. Y. F. Huang, P. J. Kooyman and M. T. Koper, Nat. Commun., 2016, 7, 12440 CrossRef CAS PubMed.
  14. J. C. Dong, X. G. Zhang, V. Briega Martos, X. Jin, J. Yang, S. Chen, Z. L. Yang, D. Y. Wu, J. M. Feliu, C. T. Williams, Z. Q. Tian and J. F. Li, Nat. Energy, 2018, 4, 60–67 CrossRef.
  15. S. A. Wasileski, M. T. M. Koper and M. J. Weaver, J. Am. Chem. Soc., 2002, 124, 2796–2805 CrossRef CAS PubMed.
  16. D. Y. Wu, X. M. Liu, S. Duan, X. Xu, B. Ren, S. H. Lin and Z. Q. Tian, J. Phys. Chem. C, 2008, 112, 4195–4204 CrossRef CAS.
  17. S. Y. Ding, B. J. Liu, Q. N. Jiang, D. Y. Wu, B. Ren, X. Xu and Z. Q. Tian, Chem. Commun., 2012, 48, 4962–4964 RSC.
  18. K. Iida, M. Noda and K. Nobusada, J. Chem. Phys., 2014, 141, 124124 CrossRef PubMed.
  19. E. Santos, A. Lundin, K. Pötting, P. Quaino and W. Schmickler, Phys. Rev. B: Condens. Matter Mater. Phys., 2009, 79, 235436 CrossRef.
  20. A. Y. Lozovoi, A. Alavi, J. Kohanoff and R. M. Lynden-Bell, J. Chem. Phys., 2001, 115, 1661–1669 CrossRef CAS.
  21. J. K. Nørskov, J. Rossmeisl, A. Logadottir, L. Lindqvist, J. R. Kitchin, T. Bligaard and H. Jonsson, J. Phys. Chem. B, 2004, 108, 17886–17892 CrossRef.
  22. J. S. Filhol and M. Neurock, Angew. Chem., Int. Ed., 2006, 45, 402–406 CrossRef CAS PubMed.
  23. S. Venkatachalam and T. Jacob, Z. Phys. Chem., 2007, 221, 1393–1406 CrossRef CAS.
  24. S. Schnur and A. Groß, Catal. Today, 2011, 165, 129–137 CrossRef CAS.
  25. J. Le, M. Iannuzzi, A. Cuesta and J. Cheng, Phys. Rev. Lett., 2017, 119, 016801 CrossRef PubMed.
  26. J. W. Halley, A. Mazzolo, Y. Zhou and D. Price, J. Electroanal. Chem., 1998, 450, 273–280 CrossRef CAS.
  27. H. F. Wang and Z. P. Liu, J. Phys. Chem. C, 2009, 113, 17502–17508 CrossRef CAS.
  28. K. Mathew, R. Sundararaman, K. Letchworth-Weaver, T. Arias and R. G. Hennig, J. Chem. Phys., 2014, 140, 084106 CrossRef PubMed.
  29. K. Mathew and R. G. Hennig, arXiv preprint arXiv:1601.03346, 2016.
  30. P. Deshlahra, E. E. Wolf and W. F. Schneider, J. Phys. Chem. A, 2009, 113, 4125–4133 CrossRef CAS PubMed.
  31. A. Y. Lozovoi and A. Alavi, J. Electroanal. Chem., 2007, 607, 140–146 CrossRef CAS.
  32. S. Duan, P. P. Fang, F. R. Fan, I. Broadwell, F. Z. Yang, D. Y. Wu, B. Ren, C. Amatore, Y. Luo, X. Xu and Z. Q. Tian, Phys. Chem. Chem. Phys., 2011, 13, 5441–5449 RSC.
  33. J. F. Li, Y. J. Zhang, A. V. Rudnev, J. R. Anema, S. B. Li, W. J. Hong, P. Rajapandiyan, J. Lipkowski, T. Wandlowski and Z. Q. Tian, J. Am. Chem. Soc., 2015, 137, 2400–2408 CrossRef CAS PubMed.
  34. R. Jinnouchi, T. Hatanaka, Y. Morimoto and M. Osawa, Phys. Chem. Chem. Phys., 2012, 14, 3208–3218 RSC.
  35. L. W. H. Leung, A. Wieckowski and M. J. Weaver, J. Phys. Chem., 1988, 92, 6985–6990 CrossRef CAS.
  36. J. R. Lombardi and R. L. Birke, J. Chem. Phys., 2007, 126, 244709 CrossRef PubMed.
  37. S. N. Steinmann, C. Michel, R. Schwiedernoch and P. Sautet, Phys. Chem. Chem. Phys., 2015, 17, 13949–13963 RSC.
  38. M. Mamatkulov and J. S. Filhol, Phys. Chem. Chem. Phys., 2011, 13, 7675–7684 RSC.
  39. M. Neurock, M. Janik and A. Wieckowski, Faraday Discuss., 2008, 140, 363–378 RSC.
  40. B. Peng, H. F. Wang, Z. P. Liu and W. B. Cai, J. Phys. Chem. C, 2010, 114, 3102–3107 CrossRef CAS.
  41. M. Liu, Y. Pang, B. Zhang, P. De Luna, O. Voznyy, J. Xu, X. Zheng, C. T. Dinh, F. Fan, C. Cao, F. P. de Arquer, T. S. Safaei, A. Mepham, A. Klinkova, E. Kumacheva, T. Filleter, D. Sinton, S. O. Kelley and E. H. Sargent, Nature, 2016, 537, 382–386 CrossRef CAS PubMed.
  42. S. Gao, Y. Lin, X. Jiao, Y. Sun, Q. Luo, W. Zhang, D. Li, J. Yang and Y. Xie, Nature, 2016, 529, 68–71 CrossRef CAS PubMed.
  43. S. C. Chang and M. J. Weaver, Surf. Sci., 1990, 238, 142–162 CrossRef CAS.
  44. I. Villegas and M. J. Weaver, J. Chem. Phys., 1994, 101, 1648–1660 CrossRef CAS.
  45. C. Jung, B. Ku, J. Kim and C. K. Rhee, Chem. Commun., 2006, 45, 2191–2193 RSC.
  46. C. Jung, J. Kim and C. K. Rhee, Langmuir, 2007, 23, 9495–9500 CrossRef CAS PubMed.
  47. A. Rodes, R. Gómez, J. M. Feliu and M. J. Weaver, Langmuir, 2000, 16, 811–816 CrossRef CAS.
  48. A. López-Cudero, A. Cuesta and C. Gutiérrez, J. Electroanal. Chem., 2005, 579, 1–12 CrossRef.
  49. P. J. Feibelman, B. Hammer, J. K. Nørskov, F. Wagner, M. Scheffler, R. Stumpf, R. Watwe and J. Dumesic, J. Phys. Chem. B, 2001, 105, 4018–4025 CrossRef CAS.
  50. L. Schimka, J. Harl, A. Stroppa, A. Grüneis, M. Marsman, F. Mittendorfer and G. Kresse, Nat. Mater., 2010, 9, 741–744 CrossRef CAS PubMed.
  51. B. D. Dunnington and J. R. Schmidt, J. Chem. Theory Comput., 2012, 8, 1902–1911 CrossRef CAS PubMed.
  52. G. Kresse and J. Furthmuller, Phys. Rev. B: Condens. Matter Mater. Phys., 1996, 54, 11169–11186 CrossRef CAS PubMed.
  53. G. Kresse and D. Joubert, Phys. Rev. B: Condens. Matter Mater. Phys., 1999, 59, 1758–1775 CrossRef CAS.
  54. B. Hammer, L. B. Hansen and J. K. Nørskov, Phys. Rev. B: Condens. Matter Mater. Phys., 1999, 59, 7413–7421 CrossRef.
  55. A. Cuesta, Surf. Sci., 2004, 572, 11–22 CrossRef CAS.
  56. M. J. Weaver, Langmuir, 1998, 14, 3932–3936 CrossRef CAS.
  57. R. Sundararaman and K. Schwarz, J. Chem. Phys., 2017, 146, 084111 CrossRef PubMed.
  58. R. Sundararaman, K. Letchworth-Weaver and K. A. Schwarz, J. Chem. Phys., 2018, 148, 144105 CrossRef PubMed.
  59. J. F. Li, Y. F. Huang, Y. Ding, Z. L. Yang, S. B. Li, X. S. Zhou, F. R. Fan, W. Zhang, Z. Y. Zhou, D. Y. Wu, B. Ren, Z. L. Wang and Z. Q. Tian, Nature, 2010, 464, 392–395 CrossRef CAS PubMed.


Electronic supplementary information (ESI) available. See DOI: 10.1039/c9sc05429d

This journal is © The Royal Society of Chemistry 2020