To model the behavior of the H adsorption region, we carried out DFT calculations on three-layer slab models using DMol3 (Biovia), with geometry optimization with semicore pseudopotentials and energy calculations with all-electron scalar relativistic accuracy.7 This approach has been found to provide realistic adsorption energies. The present work tests its ability to provide realistic trends in the Fermi energy EFas a function of H coverage alone, since it has been shown by Ishikawa and coworkers that the (111) surface with underpotentially deposited hydrogen HUPD is highly hydrophobic.7 This approach can also accurately model the vacuum experiment, in which the work function change is measured vs. H coverage.8
We also plot the integrated H adsorption charges taken directly from CV, which is the most straightforward way to compare with the calculations. We carried out calculations for a Pt(111)-6x6 3-layer model with sequentially added H atoms, up to unity coverage (36 atoms), with the further addition of 12 H atoms, corresponding to overpotentially deposited H (HOPD). This region is difficult to evaluate with CV, although it is observable spectroscopically.9 The coverage can be converted straightforwardly to a charge, but EF must be empirically transformed to match the electrochemical potential. Interestingly, the calculated curve has a delayed onset, similar to that observed in vacuum.8 Calculations were carried out for a Pt/Pt3Co/Pt3Co model, in which the second and third layers both have the bulk composition (Fig. 1). The curve is shifted downward due to the effect of Co. Finally, calculations were carried out for a realistic Pt/Co/Pt3Co(111) model. Since these are extremely time-consuming, only the bare surface and the unity coverage cases have successfully completed thus far, but these show a further downward shift.
The present results show for the first time that the electrochemical behavior of a Pt-Co alloy (111) single crystal can be understood from first principles. Calculations are also in progress to model the OH adsorption region, as well as the catalytic reactions.
This work was supported by funds for the “Superlative, Stable, and Scalable Performance Fuel Cell” (SPer-FC) project from the New Energy and Industrial Technology Development Organization (NEDO) of Japan.
References
1. T. Toda, H. Igarashi, H. Uchida, and M. Watanabe, J. Electrochem. Soc.146, 3750 (1999).
2. G. Shi, H. Yano, D. A. Tryk, A. Iiyama and H. Uchida, ACS Catal., 7, 267 (2017).
3. G. Shi, H. Yano, D. A. Tryk, S, Noharaa, and H. Uchida, Phys. Chem. Chem. Phys., 21, 2861 (2019).
4. S. Kobayashi, M. Wakisaka, D. A. Tryk, A. Iiyama, and H. Uchida,J. Phys. Chem. C,121, 11234 (2017).
5. S. Kobayashi, M. Aoki, M. Wakisaka, T. Kawamoto, R. Shirasaka, K. Suda, D. A. Tryk, J. Inukai, T. Kondo, and H. Uchida, ACS Omega, 3, 154 (2018).
6. V. R. Stamenkovic, B. Fowler, B. S. Mun, G. F. Wang, P. N. Ross, C. A. Lucas, and N. M. Markovic, Science, 315, 493 (2007).
7. Y. Ishikawa, J. J. Mateo, D. A. Tryk, and C. R. Cabrera, J. Electroanal. Chem., 607, 37 (2007).
8. K. Christmann, G. Ertl, and T. Pignet, Surf. Sci., 54, 365 (1976).
9. K. Kunimatsu, H. Uchida, M. Osawa, and M. Watanabe, J. Electroanal. Chem.,587, 299 (2006).
Fig. 1. (Upper panel) CV for Pt(111) (black) and Pt-skin/Pt73Co27(111) (blue) electrodes in N2-saturated 0.1 M HClO4 at 50 mV s-1. (Lower panel) Integrated charges (solid lines) based on the CV, with calculated points (open symbols) for Pt/Pt3Co/Pt3Co(111) (red circles) and Pt/Co/Pt3Co (blue squares) models.
