
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39251693
71703
10.1038/s41598-024-71703-w
Article
Hydrogen adsorption on fcc metal surfaces towards the rational design of electrode materials
Lousada Cláudio M. cmlp@kth.se

Kotasthane Atharva M.
https://ror.org/026vcq606 grid.5037.1 0000 0001 2158 1746 Department of Materials Science and Engineering, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
9 9 2024
9 9 2024
2024
14 2097210 4 2024
30 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The successful large-scale implementation of hydrogen as an energy vector requires high performance electrodes and catalysts made of abundant materials. Rational materials design strategies are the most efficient means of reaching this goal. Here we present a study on the adsorption of H-atoms onto fcc transition metal surfaces and propose descriptors for the rational design of electrodes and catalysts by means of correlations between fundamental properties of the materials and among other properties, their experimentally measured performance as hydrogen evolution electrodes (HEE). A large set of quantum mechanical modelling data at the DFT level was produced, covering the adsorption of H-atoms onto the most stable surfaces (100), (110) and (111) of: Ag, Au, Co, Cu, Ir, Ni, Pd, Pt and Rh. For each material and surface, a coverage dependent set of minimum energy structures was produced and chemical potentials for adsorption of H-atoms were obtained. Averaging procedures are here proposed to approach modelling to the experiments. Several correlations between the computed data and experimentally measured quantities are done to validate our methodology: surface plane dependent adsorption energies, chemical potentials and experimentally determined surface energies and work functions. We search for descriptors of catalytic activity by testing correlations between the DFT data obtained from our averaging procedures and experimental data on HEE performance. Our methodology allows us to obtain linear correlations between the adsorption energy of H-atoms and the exchange current density (i0) in a HEE, avoiding the volcano-like plots. We show that the chemical potential has limitations as a descriptor of i0 because it reaches an early plateau in terms of i0. Simple quantities obtained from database data such as the first stage electronegativity (χ) as devised by Mulliken has a strong linear correlation i0. With a quantity we denominate modified second-stage electronegativity (χ2m) we can reproduce the typical volcano plot in a correlation with i0. A theoretical and conceptual framework is presented. It shows that both χ and χ2m, that depend on the first ionization potential, second ionization potential and electron affinity of the elements can be used as descriptors in rational design of electrodes or of catalysts for hydrogen systems.

Subject terms

Electrocatalysis
Electronic structure
Heterogeneous catalysis
Erasmus+Royal Institute of TechnologyOpen access funding provided by Royal Institute of Technology.

issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Dissociative hydrogen adsorption or H-atom adsorption at metal surfaces is a phenomenon that determines the performance of innumerous systems of high technological and scientific importance. The most obvious connection between this atomic scale process with practical applications is within hydrogen production and usage. This encompasses an array of electrocatalytic systems such as electrodes in diverse electrolysis reactors1 or in fuel cells2, but also photocatalytic systems where there has been a considerable effort to develop H2(g) production from the sunlight assisted dissociation of H2O at the surface of semiconductors doped with metallic inclusions that function as co-catalysts3,4. H-atom adsorption is also a key parameter in corrosion and hydrogen absorption and embrittlement by metals5,6. Additionally, diverse important catalytic processes involve H-atom transfer reactions between molecular species at the surface of high-performance catalysts7.

Rational materials design is a strategy for the development of new functional, high-performance materials that is gaining momentum because of the link it can provide between properties and function8. This topic has attracted considerable attention because of the developments in modern experimental and computational techniques that allow more detailed studies of materials at the atomic scale9. To correlate atomic scale properties with processes to generate descriptors of performance for usage in rational design is of utmost importance for the successful development of better materials10. Improving the materials used in high performance applications is an important task that requires access to high-quality data that can be rationalized and correlated to produce descriptors for the rational design of those materials11. Especially relevant to this task, is to first develop correlations between properties and function for simpler materials—single component with simpler geometries—before embarking on the study of complex multicomponent-multi-structure systems12,13.

Despite the relevance of H-atom adsorption at metal surfaces, there is no account that systematically investigates the effects of the substrate in terms of both surface planes and increasing coverage of H-atoms on the adsorption energies (ΔEads) and on the driving force for adsorption. Understanding these effects and the origin for the driving force for adsorption and for the formation of chemical bonds is essential to understand and predict performance of materials in real systems. The success of the d band model for this task is unquestionable14,15. However, this model requires a rather complex analysis to be put in practice16. Because of its inherent complexity, its physical significance may not be directly understood by researchers more oriented towards engineering that focus on system design and upscaling. For these tasks, the existence of simpler descriptors of performance that can be used in rational materials design and in the first stages of materials selection, that can link simple properties with performance, and that can be understood in terms of simple quantities that are intrinsic to the materials is important12,17,18.

The available bulk literature data on H-atom adsorption onto different surfaces of common transition metals has been obtained with different methods and by different authors often under different conditions as detailed in the Supplementary Information section to this paper. This has an impact especially in the computational data whose quality relies on the cancelation of errors of differential quantities such as ΔEads. Computational data obtained with different computational methodologies has different error magnitudes and deviates from the real values in different directions—some methods systematically underestimate while others overestimate—which results in large data dispersions that do not allow accurate analysis, rationalization, and derivation of performance descriptors. Supplementary Information Figure S1 shows how ΔEads data for adsorption of H-atoms on Cu as a function of surface coverage defined as fraction of monolayer (ML), varies depending on the data source.

In general, for transition metals with fcc structure, the low Miller index surfaces (100), (110), and (111) are the most stable and most frequently occurring in polycrystals19. For practical applications these are thus the most relevant surfaces in terms of the performance of the systems where these materials are used for their surface properties. For hydrogen systems that depend on surface reactions, the adsorption energy of the hydrogen atom (ΔEadsH) is among the most significant parameters to consider. This quantity depends on both the geometrical and electronic properties of the adsorbent surfaces. In general, the traditional decoupling of geometric and electronic effects in adsorption can be seen as a necessary oversimplification, because the electronic structure of the surfaces depends on their geometries and vice-versa. Properties that depend on the surface plane such as work function (W) among others, are both determined by the surface structure and by bulk properties such as the ionization potential (IP), electron affinity (EA) and consequently by the electronegativity (χ). We have previously shown that for transition metal and lanthanide oxides in solution, very good correlations between experimentally measured catalytic activity and computational data could be obtained against the IP, EA and χ of the cations. These correlations could predict catalytic properties such as energy barriers and strength of adsorption of products and intermediates in processes where free radicals were generated at the surfaces12. Based on this principle, we have further derived a descriptor of O-atom adsorption onto clean and doped Al surfaces with different dopants. The descriptor is a modified Mulliken electronegativity that involves the second ionization potential IP2 of the dopant atom and their first EA20. The very good predictive power of this function which is based on textbook data, allowed us to predict trends in adsorption energy of O-atoms using quantities solely intrinsic to the dopant atoms at the surface. The systematic analysis of the electronic structure of the surface-dopant-O-atom complex revealed how the bands and newly formed electronic states sequentially follow a trend that depends on IP and EA and explains why these quantities could be used for devising predictive descriptors of adsorption20. The success of this descriptor of adsorption is because it includes quantities dependent on filled or valence states, the IP2, and on empty states, the EA. This makes the descriptor able to account for donation and back donation upon formation of chemical bonds at the surface. Because for adsorption of the same adsorbate onto different surfaces, the reference electron energy of the adsorbate is the same, the descriptor will grasp the relative differences between filled and empty states for different surfaces or different materials.

In this work we tested the usage of descriptors for predicting trends in adsorption of H-atoms on a series of transition metals of technological relevance that have fcc structure. We show that the first-stage electronegativity and the modified second-stage electronegativity of the metals are good descriptors of trends in adsorption at (111), (110) and (100) surfaces of fcc metals, and that these descriptors can be used in tasks of rational design or materials selection. The first-stage electronegativity of the metals correlates linearly with experimentally measured exchange current density, avoiding the volcano-like plot. These results also highlight the importance of averaging over different adsorbate coverages and over different surface planes for the modelling data to be closer to the experimental conditions in engineered catalytic reactor systems.

Methods

Literature data

Methodologies for literature data collection on hydrogen adsorption are detailed in the Supplementary Information to this paper.

Computational details

Density functional theory (DFT) calculations were performed using the Vienna ab initio simulation package (VASP 5.4.1)21 with the Perdew–Burke–Ernzerhof (refs 22,23.) (PBE) exchange–correlation functional with ultrasoft pseudopotentials of the projector augmented wave24,25 (PAW) type. The PBE functional has shown good accuracy for describing the adsorption of H-atoms onto metal surfaces26–28, and also for surface and bulk properties of metals29,30. We tested the inclusion of descriptors of van der Waals interactions which are important in molecular adsorption, using the D3 correction31 as PBE-D3. However, the adsorption energies of H-atoms at the surfaces of these transition metals computed at the PBE level are orders of magnitude larger than the D3 correction to those energies which results in negligible relative differences— < 0.002 eV on average—in the adsorption energies obtained with and without the D3 correction. Because of this and to facilitate direct comparison with the literature, the data here presented has been computed at the PBE level. The energies herein reported are electronic energies at 0 K which allow accurate direct comparisons between preferred binding sites with similar chemical environments32–34. The (100), (110) and (111) fcc surfaces were modelled with p(2 × 2) surfaces using the supercell models shown in Fig. 1.Fig. 1 Supercell models with p (2 × 2) symmetry used to study adsorption of H-atoms onto the fcc transition metals: Ag, Au, Co, Cu, Ir, Ni, Pd, Pt, Rh.

The supercells with periodic boundary conditions have a thickness of five atomic layers for (100) and (110) along the axis perpendicular to the surface planes and four atomic layers for (111), all with vacuum layers of 15 Å. These models have been employed in our previous works for the study of adsorption of H-atoms and are sufficient to reduce spurious interactions between the periodically repeating images27,28,35. Adsorption of H-atoms was investigated at one of the sides of the slabs. During optimization of the geometries, for (100) and (110), the three bottom layers were kept fixed while the remaining top layers as well as the H-atoms were allowed to relax. For (111) the two bottom layers were kept fixed while the remaining top layers and the H-atoms were allowed to relax. For geometry optimizations a plane wave cutoff of 500 eV and a Γ centered k-point mesh of (6 × 6 × 1) in the Monkhorst–Pack sampling scheme were used with a Gaussian smearing of 0.1 eV width. The self-consistent field (SCF) electronic energies were considered as converged when the change was smaller than 1 × 10−5 eV between cycles and the force acting on each of the atoms smaller than 0.002 eV Å−1. The study of the adsorption of H-atoms started with a single H-atom per supercell (0.25 ML). At this stage we searched for adsorption site preference. Subsequently the coverage was increased sequentially until one monolayer (1 ML). The definition of ML here used is that recommended by IUPAC36.

The adsorption energies of the H-atoms were determined as1 ΔEadsH=E(nH -M\_slab)-nEH+EM_slab

where E(nH-M_slab) is the energy of n H-atoms bound to a specific surface of a metal M; EH is the energy of an H-atom in vacuum and EM_slab is the energy of the supercell containing the bare surface of M. A more negative value for ΔEadsH implies stronger adsorption. The values of ΔEadsH are normalized per H-atom considering the following correspondence between ML and the number of adsorbed H-atoms per supercell: 0.25 ML = 1H; 0.5 ML = 2H: 0.75 ML = 3H and 1 ML = 4H. Two different average quantities have been devised to express the ΔEadsH. The average over the coverage for each surface plane for each metal is defined as2 ΔE¯adsH100,ΔE¯adsH110,ΔE¯adsH111=ΔEadsH0.25+ΔEadsH0.5+ΔEadsH0.75+ΔEadsH14100,110,111

where ΔĒadsH(100), ΔĒadsH(110), ΔĒadsH(111) are the average adsorption energies for the H-atom for each of the surfaces for a given metal; and ΔEadsH(0.25), ΔEadsH(0.5), ΔEadsH(0.75), and ΔEadsH(1) are the corresponding adsorption energies of the H-atom for each surface coverage—0.25 to 1 ML—expressed per H-atom. Because the differences in the surface energies for the three low Miller index surfaces of fcc metals are in general not large37, in a polycrystal all these surfaces can be present, even if in different extents, but will thus be responsible for the overall catalytic activity of the material. Hence, the average adsorption energy of the H-atoms for each metal here studied (ΔĒadsH) was obtained as3 ΔE¯adsH=ΔE¯adsH100+ΔE¯adsH110+ΔE¯adsH1113

The electronic energy dependent chemical potential µi is defined as a function of ML fraction up to 1 ML as4 μi=∂ΔEads∂MLi100,110,(111)

where ΔEads is the adsorption energy of i H atoms. This expression does not include the zero-point vibrational energies because in previous works we found that this quantity does not contribute to changes in the trends of adsorption nor absorption for similar geometries28,38. The expression does not include also corrections for pressure nor temperature. The values of µ were thus determined from plots of ΔEadsHi vs ML for each surface leading to the quantities µi(100), µi(110) and µi(111) which are the coverage dependent chemical potentials for H-atom adsorption at the corresponding surfaces of each metal. Because each 0.25 increase in ML corresponds to an H-atom, the µi expressed per ML can be easily expressed per number of H-atoms by computing µi/4. The average chemical potential µavg for adsorption of H at the surface of each metal here considered is determined as5 μavg=μi(100)+μi110+μi(111)3

Results and discussion

Analysis of literature data

It is thus expectable that, if for a given material the electronic properties of each surface, for example (100), (110) and (111), change in a consistent way and can be linked to certain trends in adsorption, the intrinsic properties of the materials such as IP, EA, and χ can also be related to the type of surface by a function. Correlations between surface plane dependent work function W, for the (100), (110) and (111) surfaces and bulk IP and EA that are independent of the surface are shown in the Supplementary Information Figure SI2 for a series of transition metals with fcc structure. The correlations between W and IP, EA and χ are not linear as expected because the electronic states of the surfaces that influence W depend on complex quantum mechanical effects that do not correlate linearly with IP, EA and χ. However, the fact that the data of Figure SI2 varies for the different surface planes in a tendentially linear fashion suggests that it is possible to correlate certain aspects of the performance of the surfaces with bulk properties of the materials if we obtain an adequate parameter that “adapts” the bulk properties χ, EA and IP, to the specific surface plane, taking into account a specific surface process, that itself it is surface plane dependent. Such simplification that does not depend on descriptors that need complex computational data would allow a fast track for evaluation of performance and for the design of materials and systems that are aimed at specific surface performance.

Coverage dependent adsorption energies

Adsorbate coverage is an important parameter that determines the adsorption energy39,40. We have previously demonstrated that even for defective surfaces, adsorbate coverage has an important role both in the formation and in the magnitude of dipoles and in the strength of adsorption17,41,42. Long-range adsorbate–adsorbate effects have been determined to amount to 15 Å in Cu for coverages as low as 1/16 ML and 1/32 ML. For S-atoms, these interactions can affect the magnitude of the adsorption energy by 0.4 eV and the magnitude and sign of local dipole momenta—from + 1.217 D for adsorption at neighboring sites, to − 2.390 D for adsorption at a distance of 18 Å17. For applications in high performance catalytic reactors however, higher coverages closer to 1 ML have to be considered. We started by studying the coverage dependence of ΔEadsH from 0.25 to 1 ML. The resulting data is shown in Supplementary Information Figure SI3 and the obtained chemical potentials µi are summarized in Table 1. Table 1 Electronic energy dependent chemical potentials µi expressed as coverage of H-atoms at the corresponding surfaces from 0.25 to 1 ML. µi is defined in Eq. (4).

	µi (ev)/ML	b (eV)	
(hkl)	100	110	111	100	110	111	
Ag	 − 8.166	 − 8.856	 − 4.230	0.210	0.588	 − 0.604	
Au	 − 9.011	 − 7.534	 − 6.392	0.442	 − 0.123	 − 0.339	
Co	 − 11.199	 − 10.613	 − 4.615	0.114	 − 0.581	 − 1.726	
Cu	 − 9.325	 − 10.004	 − 4.749	0.009	0.553	 − 1.381	
Ir	 − 11.016	 − 10.667	 − 10.464	 − 0.125	 − 0.076	 − 0.068	
Ni	 − 12.618	 − 11.337	 − 7.837	1.215	0.506	 − 0.309	
Pd	 − 11.384	 − 8.571	 − 8.722	0.968	 − 0.089	 − 0.142	
Pt	 − 10.559	 − 11.131	 − 10.282	 − 0.163	 − 0.093	 − 0.139	
Rh	 − 11.142	 − 11.511	 − 9.388	0.554	0.827	 − 0.121	

The dependence of ΔEads with ML is linear for all metals and surface planes except for Cu(111) and Co(111) most likely due to spin effects as these surfaces are known to undergo adsorbate induced transitions that can lead to non-zero spin momenta and to localized dipoles upon adsorption17,43–45. The extent of surface reconstruction for the same surface planes at the same coverage of H-atoms is different for different metals as shown in Figure SI4 for Ag and Co. This highlights the importance of studying both different coverages and also allowing the models to relax during geometry optimization to find the real minima. The chemical potentials expressed as a function of surface energy reveal how the process of adsorption resembles that of reconstructing the bulk material: chemical bonds between a surface and an adsorbate which are more similar to the bonds in the bulk metal lead to stronger correlations with the surface energy. The correlations between µi and the surface energies of the fcc metals are shown in Fig. 2.Fig. 2 µi (eV/ML) as a function of the surface energies (eV) for the series of fcc metals here studied. The surface energy of Co(110) is non-existent in the literature and has not been included.

There is a strong dependency of µi with the surface energy for the sets of surfaces (110) and (111), while for the set (100) the correlation is weak. The 110 surface is the more open of the three which leads to more degrees of freedom for bonding with H-atoms in terms of their disposition ie. top, bridge hollow etc., as highlighted in Figure SI4. The effect of the changes in binding site preference reflects in the correlation between surface coverage and µi. However, the 110 surface is not the surface for which the binding site preference changes the most with coverage. For the 100 set of surfaces, the changes in adsorption site preference as the coverage of H-atoms changes are the most extreme from the complete set of metals here studied. These effects can be attributed to adsorbate–adsorbate effects that change both magnitude and direction of local dipole momenta at different coverages17. These results show the importance of allowing adsorbate geometry relaxation in computations of adsorption dependent such quantities. The consequences of the more extensive changes in binding site preference as a function of coverage (100), are different correlations between the adsorption energy and the descriptors of adsorption that we are devising in this work as it will be shown in the next section.

Descriptors of adsorption energy

Simple descriptors of ΔEads that predict relative differences in this quantity in catalytic reactors could boost the design of catalysts46–48. Predictors of relative ΔEadsH can thus help design better electrodes for electrolysis and fuel cells among other systems where ΔEadsH is of relevance46,49. In this context, we tested how well our averaging methodology can grasp the non-linearities observed in experiments by correlating experimentally measured Ws of the metals and our average adsorption energies of the H-atom ΔĒadsH (eV) obtained according to Eq. (3). The resulting data is shown in Fig. 3. The data shows that our averaging strategy grasps experimental conditions better than other modelling approaches that focus on single coverages as the linear correlations show.Fig. 3 Left panel: experimentally measured work functions W for polycrystals of the metals50 and their correlation with the average adsorption energies of the H-atom ΔĒadsH (eV) obtained according to Eq. (3), over coverage and surface planes. Average (avg) W obtained by averaging the W of the (100), (110) and (111) surfaces retrieved from literature51.

Additionally, the correlation between ΔĒadsH and the two sets of experimental W shows that averaging adsorption energies over coverage and over the surfaces with lower surface energies is a good approach to correlate intrinsic physical properties with adsorption performance. This is because in the polycrystals the three surfaces (100), (110) and (111) can occur because the differences in their energies are small37. Averaging over data obtained for each surface is thus a better approximation for polycrystals. The linear correlations between ΔEadsH at each coverage for each type of surface and the W of each surface are weak—0.11 < R2 < 0.48—for most of the surfaces and coverages as the collection of data in Supplementary Information Figure SI5 shows. The exceptions are the higher coverages 0.75 and 1 ML at the (111) surfaces and the 1 ML coverage at the (110) surfaces. This shows that W cannot account for the shifts in adsorbate binding site preference that occur when the coverage increases, as we have already discussed above. The set of (100) surfaces are those where these effects are more significant because of the larger number of possible adsorbate binding sites. This demonstrates that geometric and electronic effects go hand in hand and methods for prediction of adsorption must account for these. Averaging over both coverage and surface planes, and finding the energy minima for each geometry at each coverage accounts substantially for both effects as Fig. 3 demonstrates.

Descriptors of catalytic activity

The exchange current density, i0, in electrochemical reactions is analogous to the rate constant in chemical reactions. For a homologous series of reactions, Brønsted, Evans and Polanyi (BEP) have shown that there is a linear correlation between the transition state energies—that in turn affect the corresponding reaction rate constants—and the adsorption energies in what is now known as the BEP relation52,53. Several authors have explored this correlation for the successful prediction of catalytic performance of metal surfaces54, including ourselves, in our case for surfaces of oxides in solution12. Here we explore the correlation between i0 at the hydrogen evolution electrode (HEE) and the ΔEadsH with an approach fundamentally similar to what was tested previously, but with important differences55. The difficulties that the authors of reference55 encountered are due to the usage of data obtained with different methods for the determination of ΔEadsH which in turn have different errors. Additionally, the study considered only ¼ ML and 1 ML. In the present study we consider the overall average of the ΔEadsH, defined as ΔĒadsH that in turn is the average of the ΔEadsH obtained for each surface plane also averaging over the different coverages as described in Eq. (3). This makes the quantity ΔĒadsH grasp both surface coverage effects on the ΔEadsH and the effects of the surface index in this quantity.

In Fig. 4, the more reactive metals are at the left of the plot while the less reactive and unreactive are at the right as previously observed for a volcano plot of the activities of these metals55. However, the close to linear dependency of i0 with ΔĒadsH can be attributed to the fact that as observed for the Ws data of Fig. 3, the averaging of the adsorption energies in terms of surface planes and in terms of surface coverages grasps the experimental conditions to a better extent that if only selected surface planes or coverages are used. This is because of the highly dynamical nature of the reactions of evolution of hydrogen whose rate of formation depends on the ΔEadsH but the rate of desorption from the surface depends on the surface coverage as well on the relation between ΔEadsH and surface coverage, that we express here as the surface coverage dependent chemical potential, µi. We have previously succeed in using a similar averaging approach for the prediction of the formation of H2(g) at the surface of Cu(110), where averaging over different adsorbate conformation energies and coverages was necessary for accounting for the contributions of adsorbate dispositions that can occur in the real systems because these have similar energies35,41,42,56.Fig. 4 Dependency of the exchange current density at the hydrogen evolution electrode log (i0/A cm−2) with the adsorption energy of the H-atom ΔĒadsH (eV) obtained according to Eq. (3). The linear fit has a slope m =  − 2.6, an intercept at origin b =  − 10.59 and a goodness of fit of 0.60.

The quadratic correlation between µavg and i0 in Fig. 5 shows that a volcano-like trend of quadratic nature may apply to these two quantities. Fitting the data with linear functions leads to very poor correlations for elements whose µavg are more negative than − 9 eV/ML and a quadratic function leads to better goodness of fit. The data clearly shows two groups of materials in terms of correlation between µavg and i0. It shows also that the driving force for adsorption reaches a plateau in terms of its effect in the electrode performance.Fig. 5 Average chemical potential for H-atom adsorption (µavg) and correlation with the exchange current density at the hydrogen evolution electrode (i0). µavg were determined from Eq. (5). Experimental values of i0 retrieved from55 and references therein. Less active catalysts (circles), more reactive catalysts (diamonds).

We previously used primary quantities for the prediction of oxygen adsorption and dissociation on pure and doped Al surfaces20 and in the description of catalytic activity of oxides in solution12, correlating computational with experimental data. The first of those quantities that we will test here is the first-stage electronegativity (χ) as defined by Mulliken. The second quantity is what we called “modified second-stage electronegativity” where instead of the quantities used for the second-stage electronegativity in the definition by Mulliken, we employed the second ionization potential IP2 of the metal and the one electron affinity EA. The details are given elsewhere12,20. Both expressions are6 χ=IP+EA2

7 χ2m=IP2+EA4

where χ2m is the modified second-stage electronegativity, IP is the first ionization potential of the metal, IP2 is the second ionization potential of the metal—surface independent quantities–and EA is the electron affinity of the metal—a surface independent quantity as well. IP, IP2 and EA are quantities—for the elements—retrieved from common databases. The details of the choice of IP2 instead of IP (first ionization energies) for devising χ2m are given in our original work. In general, IP2 together with EA accounts for the different extents of donation-backdonation that occurs when the H-atom binds to the surface57–61, and the relaxation of the remaining orbitals in the d-band of the metal as described in the Newns-Anderson model. For a given adsorbate, depending on the metal, the balance between donation-backdonation of electron density will change. Our descriptor (χ2m) grasps the correct phenomenology of these processes for dopants on Al surfaces and their effect on oxygen adsorption20. For the set of fcc metals here study, the plots of the exchange current density (i0) in a hydrogen evolution electrode versus χ and χ2m are shown in Fig. 6.Fig. 6 First-stage electronegativity χ (eV) and modified second-stage electronegativity χ2m (eV) determined according to Eqs. (6) and (7) and their correlations with the exchange current density i0. χ and χ2m are obtained from the well-known quantities IP and EA. Experimental values of i0 retrieved from55 and references therein.

The data of Fig. 6 shows an excellent correlation between the first-stage electronegativity of the metals (χ) and i0. This shows that χ can in fact be used as a linear descriptor in the search for catalysts for HEE. Further support to this observation is also given by χ2m in Fig. 6. The “sideways volcano plot” resultant from correlating χ2m with i0 provides considerable insight into the phenomenological description of bonding: donation–backdonation, between the H-atom and the surfaces as it will be discussed in the next section.

Theoretical and conceptual framework of the descriptors of catalytic activity

In the Newns–Andersson model the energetic positioning of the d-band center with respect to the orbitals of the adsorbate explains certain trends in catalytic activity. The direct correlation between the model and a surface process is with regard to adsorption energies, that in turn relate with frequency of events or with rates by the BEP relation leading to an indirect correlation between the d-band derived quantities and the frequency or rates of events. Bonding between an adsorbate to a series of surfaces can thus be rationalized in terms of how high or low the d-band center and/or edge for the different materials lay with respect to one another. However, this model does not account directly for the role of electrons of the metal that lay in other orbitals such as p and s and nor accounts for effects of Pauli exclusion for example. The relationship between the energies of the d-band and the p and s bands of the metal is considered only indirectly because certain shifts in the d-band can lead to shifts in the p and s energies, but this is not a universal rule, it is in fact far from it16,62,63. The details of this model and how the chemical bonding between different surfaces and adsorbates come into play have been extensively described in reference texts62,64,65.

In our model here presented we use quantities that despite their simplicity contain the necessary information for the description of donation and backdonation from and to adsorbates. Many of the data here computed is correlated with polycrystalline and “poly-Miller-index” surfaces. The averaging over coverage and surface plane that we have employed is necessary to obtain such correlations which in turn show that there are primary electronic properties of the materials which are surface plane independent and still determinant for the description of the activity of the metal as a HEE. However, such primary quantities require an algebraic workout to obtain a descriptor that contains the correct assessment of electron donation and backdonation to and from the adsorbate, and that weighs these two processes of donation and backdonation. The typical quantities devised for example in conceptual DFT66, are related with where the electron orbitals are located energetically. As we previously employed for predicting reactivity of transition metal and actinide oxides12, this allows a description of hardness, softness and similar properties that depend on the work necessary to attach and detach electrons. However, the driving forces for the loss or gain of electrons are equally well described by the IP, IP2 and EA provided the correct algebraic expression is used.

The plots correlating the primary and secondary quantities used to determine the first-stage and second-stage electronegativities, χ and χ2m respectively, and the average adsorption energies of H-atoms, ΔĒadsH are given Fig. 7.Fig. 7 First ionization potential IP, second ionization potential IP2 and electron affinity EA and their correlation with the average adsorption energies ΔĒadsH of the H-atom to the corresponding metals obtained according to Eq. 3. The values of IP, IP2 and EA used here are retrieved from common databases for the elements and are surface independent properties.

The first ionization potential IP accounts for the energy of the highest occupied electronic states of the metal, with no distinction between orbitals, d or p. IP accounts for the energy—or work necessary if we wish to express so—for the electron transfer from the metal to occur. The electron affinity EA accounts for the ability of the metal to accommodate electrons. This quantity can be positive or negative but for metals it is mostly positive, meaning that energy is released when the metal atom accepts an incoming electron. The second ionization potential, IP2 accounts for the energy necessary to remove a second electron from the metal. This quantity includes thus the relaxation of the metal orbitals after the first electron has been removed which in turn account for electronic relaxation effects in the donation–backdonation to and from the adsorbate. IP and IP2 are largely related but include different physical descriptors of the electronic structure and this is demonstrated by the different degrees of correlation and signs of the slopes of their plots versus ΔĒadsH in Fig. 7. The first-stage electronegativity χ that contains IP and EA will thus account for the balance between the energy gained and lost in the donation-backdonation to and from the H-atom to the surfaces. The second-stage electronegativity χ2m will account for relaxation effects of the orbitals of the metal in the bonding with the adsorbate. This is because it accounts for the ease of removal of a second electron from the metal, which in turn is a quantity that contains information about the ease of relaxation of the highest energy electronic bands of the metal. Interestingly, the quantities IP and EA, despite leading to poor correlations with the adsorption energy of the H-atoms as shown in Fig. 7, are included both in the linear correlation between χ and i0 and the volcano like correlation between χ2m and i0 shown in Fig. 6 that result in strong correlations. This shows that χ and χ2m contain the essential information of electronic structure features of the metals that determine the dynamical factors that result in different i0 in the HEE. χ and χ2m have thus the potential to be used as descriptors in the search for HEE materials and other electrode or catalytic materials where adsorption of H-atoms is a key process. The simplicity of these functions means that they can be easily implemented in machine learning algorithms for catalyst screening and design.

Conclusions

The search for a large set of literature data on H-atom adsorption on fcc metals shows that data computed by different authors with different methods have uncertainties that do not allow the creation of statistically significant sets of data for usage in analysis of properties with good accuracy. Thus we present here a large set of DFT data on H-atom adsorption on the surfaces of fcc metals that is both coverage and surface plane dependent. Despite our set of data being obtained for discrete surface planes and coverages, our averaging procedure leads to data that can be correlated with experimental properties obtained for polycrystals. This shows that our averaging procedure can approximate the DFT data to the diversity of coverages and surfaces that experiments deal with. The correlations between surface energies, work functions and adsorption related properties if obtained for single surface planes and coverages without averaging are poor. This is because changing coverage implies changes in adsorbate biding site preferences for certain planes also accompanied by different extents of surface reconstruction. Our averaging methodology grasps these effects and allows the correlation between H-atom adsorption and experimental data obtained for polycrystals as shown by the linear correlations that we obtained between DFT data and exchange current density (i0) in a hydrogen evolution electrode.

We also propose descriptors of electrode activity using the first-stage (χ) and second-stage electronegativity (χ2m). Because these quantities depend on electron energies of the metals, they include the essential information that determines donation and backdonation to and from adsorbates and even effects of electronic relaxation. The result is that can χ be linearly correlated with i0 and χ2m can be correlated with i0 via a function that resembles the typical volcano plot. Because of these strong correlations, especially between i0 and χ with is strongly linear, both χ and χ2m can be used as robust descriptors of catalytic activity in search and in the design of new materials for usage as hydrogen evolution electrodes. The direct dependence of χ and χ2m on electron energy levels simplifies considerably the number of degrees of freedom to consider in the design of catalysts.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71703-w.

Acknowledgements

Cláudio M. Lousada acknowledges financial support from the Erasmus+ program. The computations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) and NAISS at the High-Performance Computing Center North (HPC2N), Umeå and at the PDC Center for High Performance Computing at the KTH—Royal Institute of Technology, Stockholm, partially funded by the Swedish Research Council through grant agreement no. 2018-05973 and partially funded by the Swedish Research Council through grant agreement no. 2022-06725.

Author contributions

CML conceptualized the study, performed the computations, treated the data, and wrote the first draft. AK performed the literature search, treated the data, and wrote the first draft. All authors have revised and approved the text.

Funding

Open access funding provided by Royal Institute of Technology.

Data availability

The data generated is available from the corresponding author upon reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Wang S Lu A Zhong CJ Hydrogen production from water electrolysis: Role of catalysts Nano Converg. 2021 8 1 4 10.1186/s40580-021-00254-x 33575919
Wang, S., Lu, A. & Zhong, C. J. Hydrogen production from water electrolysis: Role of catalysts. Nano Converg. 8(1), 4. 10.1186/s40580-021-00254-x (2021).33575919 10.1186/s40580-021-00254-x
2. Zhong CJ Nanostructured catalysts in fuel cells Nanotechnology 2010 21 6 062001 10.1088/0957-4484/21/6/062001 20065536
Zhong, C. J. et al. Nanostructured catalysts in fuel cells. Nanotechnology 21(6), 062001. 10.1088/0957-4484/21/6/062001 (2010).20065536 10.1088/0957-4484/21/6/062001
3. Cao H Tang X Tang H Yuan Y Wu J Photoinduced intermolecular hydrogen atom transfer reactions in organic synthesis Chem. Catal. 2021 1 3 523 598 10.1016/j.checat.2021.04.008
Cao, H., Tang, X., Tang, H., Yuan, Y. & Wu, J. Photoinduced intermolecular hydrogen atom transfer reactions in organic synthesis. Chem. Catal. 1(3), 523–598. 10.1016/j.checat.2021.04.008 (2021).10.1016/j.checat.2021.04.008
4. Engge Y Maulana F Nurhuda M Dissociation of water into hydrogen and oxygen through a combination of electrolysis and photocatalyst IOP Conf. Ser. Earth Environ. Sci. 2021 743 1 012054 10.1088/1755-1315/743/1/012054
Engge, Y., Maulana, F. & Nurhuda, M. Dissociation of water into hydrogen and oxygen through a combination of electrolysis and photocatalyst. IOP Conf. Ser. Earth Environ. Sci. 743(1), 012054. 10.1088/1755-1315/743/1/012054 (2021).10.1088/1755-1315/743/1/012054
5. Okonkwo PC A focused review of the hydrogen storage tank embrittlement mechanism process Int. J. Hydrog. Energy 2023 10.1016/j.ijhydene.2022.12.252
Okonkwo, P. C. et al. A focused review of the hydrogen storage tank embrittlement mechanism process. Int. J. Hydrog. Energy10.1016/j.ijhydene.2022.12.252 (2023).10.1016/j.ijhydene.2022.12.252
6. Dong L Application of atomic simulation for studying hydrogen embrittlement phenomena and mechanism in iron-based alloys Int. J. Hydrog. Energy. 2022 47 46 20288 20309 10.1016/j.ijhydene.2022.04.119
Dong, L. et al. Application of atomic simulation for studying hydrogen embrittlement phenomena and mechanism in iron-based alloys. Int. J. Hydrog. Energy. 47(46), 20288–20309. 10.1016/j.ijhydene.2022.04.119 (2022).10.1016/j.ijhydene.2022.04.119
7. Jia J Solar water splitting by photovoltaic-electrolysis with a solar-to-hydrogen efficiency over 30 Nat. Commun. 2016 7 1 13237 10.1038/ncomms13237 27796309
Jia, J. et al. Solar water splitting by photovoltaic-electrolysis with a solar-to-hydrogen efficiency over 30. Nat. Commun. 7(1), 13237. 10.1038/ncomms13237 (2016).27796309 10.1038/ncomms13237
8. Gao D Arán-Ais RM Jeon HS Roldan Cuenya B Rational catalyst and electrolyte design for CO2 electroreduction towards multicarbon products Nature Catal. 2019 2 3 198 210 10.1038/s41929-019-0235-5
Gao, D., Arán-Ais, R. M., Jeon, H. S. & Roldan Cuenya, B. Rational catalyst and electrolyte design for CO2 electroreduction towards multicarbon products. Nature Catal. 2(3), 198–210. 10.1038/s41929-019-0235-5 (2019).10.1038/s41929-019-0235-5
9. Belviso F Viewpoint: Atomic-scale design protocols toward energy, electronic, catalysis, sensing applications Inorg. Chem. 2019 58 22 14939 14980 10.1021/acs.inorgchem.9b01785 31668070
Belviso, F. et al. Viewpoint: Atomic-scale design protocols toward energy, electronic, catalysis, sensing applications. Inorg. Chem. 58(22), 14939–14980. 10.1021/acs.inorgchem.9b01785 (2019).31668070 10.1021/acs.inorgchem.9b01785
10. Cheng N Zhang L Doyle-Davis K Sun X Single-atom catalysts: From design to application Electrochem. Energy Rev. 2019 2 4 539 573 10.1007/s41918-019-00050-6
Cheng, N., Zhang, L., Doyle-Davis, K. & Sun, X. Single-atom catalysts: From design to application. Electrochem. Energy Rev. 2(4), 539–573. 10.1007/s41918-019-00050-6 (2019).10.1007/s41918-019-00050-6
11. Wang W Liu J Yu M Zhao C Wang L Acceleration of material R&D process through rational design Comput. Mater. Sci.. 2019 160 397 402 10.1016/j.commatsci.2019.01.009
Wang, W., Liu, J., Yu, M., Zhao, C. & Wang, L. Acceleration of material R&D process through rational design. Comput. Mater. Sci.. 160, 397–402. 10.1016/j.commatsci.2019.01.009 (2019).10.1016/j.commatsci.2019.01.009
12. Lousada CM Brinck T Jonsson M Application of reactivity descriptors to the catalytic decomposition of hydrogen peroxide at oxide surfaces Comput. Theor. Chem. 2015 1070 108 116 10.1016/j.comptc.2015.08.001
Lousada, C. M., Brinck, T. & Jonsson, M. Application of reactivity descriptors to the catalytic decomposition of hydrogen peroxide at oxide surfaces. Comput. Theor. Chem. 1070, 108–116. 10.1016/j.comptc.2015.08.001 (2015).10.1016/j.comptc.2015.08.001
13. Suen NT Hung SF Quan Q Zhang N Xu YJ Chen HM Electrocatalysis for the oxygen evolution reaction: Recent development and future perspectives Chem. Soc. Rev. 2017 46 2 337 365 10.1039/c6cs00328a 28083578
Suen, N. T. et al. Electrocatalysis for the oxygen evolution reaction: Recent development and future perspectives. Chem. Soc. Rev. 46(2), 337–365. 10.1039/c6cs00328a (2017).28083578 10.1039/c6cs00328a
14. Hammer B Norskov JK Why gold is the noblest of all the metals Nature 1995 376 6537 238 240 10.1038/376238a0
Hammer, B. & Norskov, J. K. Why gold is the noblest of all the metals. Nature 376(6537), 238–240. 10.1038/376238a0 (1995).10.1038/376238a0
15. Pettersson LGM Nilsson A A molecular perspective on the d-band model: Synergy between experiment and theory Topics Catal. 2014 57 1 2 13 10.1007/s11244-013-0157-4
Pettersson, L. G. M. & Nilsson, A. A molecular perspective on the d-band model: Synergy between experiment and theory. Topics Catal. 57(1), 2–13. 10.1007/s11244-013-0157-4 (2014).10.1007/s11244-013-0157-4
16. Xin H Linic S Communications: Exceptions to the d-band model of chemisorption on metal surfaces: The dominant role of repulsion between adsorbate states and metal d-states J. Chem. Phys. 2010 132 22 10.1063/1.3437609
Xin, H. & Linic, S. Communications: Exceptions to the d-band model of chemisorption on metal surfaces: The dominant role of repulsion between adsorbate states and metal d-states. J. Chem. Phys. 132, 22. 10.1063/1.3437609 (2010).10.1063/1.3437609
17. Lousada CM Johansson AJ Korzhavyi A Adsorption of hydrogen sulfide, hydrosulfide and sulfide at Cu(110)—polarizability and cooperativity effects. First stages of formation of a sulfide layer ChemPhysChem 2018 19 17 2159 2168 10.1002/cphc.201800246 29797487
Lousada, C. M., Johansson, A. J. & Korzhavyi, A. Adsorption of hydrogen sulfide, hydrosulfide and sulfide at Cu(110)—polarizability and cooperativity effects. First stages of formation of a sulfide layer. ChemPhysChem 19(17), 2159–2168. 10.1002/cphc.201800246 (2018).29797487 10.1002/cphc.201800246
18. Montemore MM Medlin JW Predicting and comparing C–M and O–M bond strengths for adsorption on transition metal surfaces J. Phys. Chem. C 2014 118 5 2666 2672 10.1021/jp5001418
Montemore, M. M. & Medlin, J. W. Predicting and comparing C–M and O–M bond strengths for adsorption on transition metal surfaces. J. Phys. Chem. C 118(5), 2666–2672. 10.1021/jp5001418 (2014).10.1021/jp5001418
19. Tran R Surface energies of elemental crystals Sci. Data 2016 3 1 160080 10.1038/sdata.2016.80 27622853
Tran, R. et al. Surface energies of elemental crystals. Sci. Data 3(1), 160080. 10.1038/sdata.2016.80 (2016).27622853 10.1038/sdata.2016.80
20. Lousada CM Korzhavyi PA Oxygen adsorption onto pure and doped Al surfaces—The role of surface dopants Phys. Chem. Chem. Phys. 2015 17 3 1667 1679 10.1039/C4CP04277H 25464123
Lousada, C. M. & Korzhavyi, P. A. Oxygen adsorption onto pure and doped Al surfaces—The role of surface dopants. Phys. Chem. Chem. Phys. 17(3), 1667–1679. 10.1039/C4CP04277H (2015).25464123 10.1039/C4CP04277H
21. Kresse G Furthmüller J Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set Phys. Rev. B 1996 54 16 11169 11186 10.1103/PhysRevB.54.11169
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54(16), 11169–11186 (1996).10.1103/PhysRevB.54.11169
22. Perdew JP Burke K Ernzerhof M Generalized gradient approximation made simple Phys. Rev. Lett. 1996 77 18 3865 3868 10.1103/PhysRevLett.77.3865 10062328
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77(18), 3865–3868 (1996).10062328 10.1103/PhysRevLett.77.3865
23. Perdew JP Burke K Ernzerhof M Generalized gradient approximation made simple Phys. Rev. Lett. 1997 78 7 1396 1396 10.1103/PhysRevLett.78.1396
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 78(7), 1396–1396 (1997).10.1103/PhysRevLett.78.1396
24. Blöchl PE Projector augmented-wave method Phys. Rev. Lett. 1994 50 24 17953 17979
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. Lett. 50(24), 17953–17979 (1994).
25. Kresse G Joubert D From ultrasoft pseudopotentials to the projector augmented-wave method Phys. Rev. B 1999 59 3 1758 1775 10.1103/PhysRevB.59.1758
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59(3), 1758–1775 (1999).10.1103/PhysRevB.59.1758
26. Hanh TTT Takimoto Y Sugino O First-principles thermodynamic description of hydrogen electroadsorption on the Pt(111) surface Surf. Sci. 2014 625 104 111 10.1016/j.susc.2014.03.006
Hanh, T. T. T., Takimoto, Y. & Sugino, O. First-principles thermodynamic description of hydrogen electroadsorption on the Pt(111) surface. Surf. Sci. 625, 104–111. 10.1016/j.susc.2014.03.006 (2014).10.1016/j.susc.2014.03.006
27. Zhao M Anderson AB Theory of hydrogen deposition and evolution on Cu(111) electrodes J. Electrochem. Soc. 2017 164 9 H691 H695 10.1149/2.1621709jes
Zhao, M. & Anderson, A. B. Theory of hydrogen deposition and evolution on Cu(111) electrodes. J. Electrochem. Soc. 164(9), H691–H695. 10.1149/2.1621709jes (2017).10.1149/2.1621709jes
28. Lousada CM Korzhavyi PA Hydrogen sorption capacity of crystal lattice defects and low Miller index surfaces of copper J. Mater. Sci.. 2020 55 15 6623 6636 10.1007/s10853-020-04459-z
Lousada, C. M. & Korzhavyi, P. A. Hydrogen sorption capacity of crystal lattice defects and low Miller index surfaces of copper. J. Mater. Sci.. 55(15), 6623–6636. 10.1007/s10853-020-04459-z (2020).10.1007/s10853-020-04459-z
29. Kurth S Perdew JP Blaha P Molecular and solid-state tests of density functional approximations: LSD, GGAs, meta-GGAs Int. J. Quantum Chem. 1999 75 4–5 889 909 10.1002/(SICI)1097-461X(1999)75:4/5<889::AID-QUA54>3.0.CO;2-8
Kurth, S., Perdew, J. P. & Blaha, P. Molecular and solid-state tests of density functional approximations: LSD, GGAs, meta-GGAs. Int. J. Quantum Chem. 75(4–5), 889–909 (1999).10.1002/(SICI)1097-461X(1999)75:4/5<889::AID-QUA54>3.0.CO;2-8
30. Haas P Tran F Blaha P Calculation of the lattice constant of solids with semilocal functionals Phys. Rev. B 2009 79 8 085104 10.1103/PhysRevB.79.085104
Haas, P., Tran, F. & Blaha, P. Calculation of the lattice constant of solids with semilocal functionals. Phys. Rev. B 79(8), 085104. 10.1103/PhysRevB.79.085104 (2009).10.1103/PhysRevB.79.085104
31. Grimme S Antony J Ehrlich S Krieg H A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H–Pu J. Chem. Phys. 2010 132 15 154104 10.1063/1.3382344 20423165
Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H–Pu. J. Chem. Phys. 132(15), 154104. 10.1063/1.3382344 (2010).20423165 10.1063/1.3382344
32. You Y-W Dissolving, trapping and detrapping mechanisms of hydrogen in bcc and fcc transition metals AIP Adv. 2013 3 1 012118 10.1063/1.4789547
You, Y.-W. et al. Dissolving, trapping and detrapping mechanisms of hydrogen in bcc and fcc transition metals. AIP Adv. 3(1), 012118. 10.1063/1.4789547 (2013).10.1063/1.4789547
33. Lousada CM Sophonrat N Yang W Mechanisms of formation of H, HO, water and of water desorption in the early stages of cellulose pyrolysis J. Phys. Chem. C 2018 122 23 12168 12176 10.1021/acs.jpcc.8b02173
Lousada, C. M., Sophonrat, N. & Yang, W. Mechanisms of formation of H, HO, water and of water desorption in the early stages of cellulose pyrolysis. J. Phys. Chem. C 122(23), 12168–12176. 10.1021/acs.jpcc.8b02173 (2018).10.1021/acs.jpcc.8b02173
34. Heinola K Ahlgren T Nordlund K Keinonen J Hydrogen interaction with point defects in tungsten Phys. Rev. B 2010 82 9 094102 10.1103/PhysRevB.82.094102
Heinola, K., Ahlgren, T., Nordlund, K. & Keinonen, J. Hydrogen interaction with point defects in tungsten. Phys. Rev. B 82(9), 094102. 10.1103/PhysRevB.82.094102 (2010).10.1103/PhysRevB.82.094102
35. Lousada CM Johansson AJ Korzhavyi PA Thermodynamics of H2O splitting and H2 formation at the Cu(110)–water interface J. Phys. Chem. C 2015 119 25 14102 14113 10.1021/acs.jpcc.5b01154
Lousada, C. M., Johansson, A. J. & Korzhavyi, P. A. Thermodynamics of H2O splitting and H2 formation at the Cu(110)–water interface. J. Phys. Chem. C 119(25), 14102–14113. 10.1021/acs.jpcc.5b01154 (2015).10.1021/acs.jpcc.5b01154
36. McNaught AD Wilkinson A Kani B Compendium of chemical terminology The Gold Book 1997 2 WileyBlackwell
McNaught, A. D. & Wilkinson, A. Compendium of chemical terminology. In The Gold Book 2nd edn (ed. Kani, B.) (WileyBlackwell, 1997).
37. Zhang J-M Ma F Xu K-W Calculation of the surface energy of FCC metals with modified embedded-atom method Appl. Surf. Sci. 2004 229 1 34 42 10.1016/j.apsusc.2003.09.050
Zhang, J.-M., Ma, F. & Xu, K.-W. Calculation of the surface energy of FCC metals with modified embedded-atom method. Appl. Surf. Sci. 229(1), 34–42. 10.1016/j.apsusc.2003.09.050 (2004).10.1016/j.apsusc.2003.09.050
38. Lousada CM Korzhavyi PA Hydrogen at symmetric tilt grain boundaries in aluminum: Segregation energies and structural features Sci. Rep. 2022 12 1 19872 10.1038/s41598-022-23535-9 36400815
Lousada, C. M. & Korzhavyi, P. A. Hydrogen at symmetric tilt grain boundaries in aluminum: Segregation energies and structural features. Sci. Rep. 12(1), 19872. 10.1038/s41598-022-23535-9 (2022).36400815 10.1038/s41598-022-23535-9
39. Lausche AC On the effect of coverage-dependent adsorbate–adsorbate interactions for CO methanation on transition metal surfaces J. Catal. 2013 307 275 282 10.1016/j.jcat.2013.08.002
Lausche, A. C. et al. On the effect of coverage-dependent adsorbate–adsorbate interactions for CO methanation on transition metal surfaces. J. Catal. 307, 275–282. 10.1016/j.jcat.2013.08.002 (2013).10.1016/j.jcat.2013.08.002
40. Kitchin JR Correlations in coverage-dependent atomic adsorption energies on Pd(111) Phys. Rev. B 2009 79 20 205412 10.1103/PhysRevB.79.205412
Kitchin, J. R. Correlations in coverage-dependent atomic adsorption energies on Pd(111). Phys. Rev. B 79(20), 205412. 10.1103/PhysRevB.79.205412 (2009).10.1103/PhysRevB.79.205412
41. Lousada CM Johansson AJ Korzhavyi PA Molecular and dissociative adsorption of water at a defective Cu(110) surface Surf. Sci. 2017 658 1 8 10.1016/j.susc.2016.11.012
Lousada, C. M., Johansson, A. J. & Korzhavyi, P. A. Molecular and dissociative adsorption of water at a defective Cu(110) surface. Surf. Sci. 658, 1–8. 10.1016/j.susc.2016.11.012 (2017).10.1016/j.susc.2016.11.012
42. Lousada CM Johansson AJ Korzhavyi PA Molecular and dissociative adsorption of water and hydrogen sulfide at perfect and defective Cu(110) surfaces Phys. Chem. Chem. Phys. 2017 19 11 8111 8120 10.1039/c6cp07732c 28267157
Lousada, C. M., Johansson, A. J. & Korzhavyi, P. A. Molecular and dissociative adsorption of water and hydrogen sulfide at perfect and defective Cu(110) surfaces. Phys. Chem. Chem. Phys. 19(11), 8111–8120. 10.1039/c6cp07732c (2017).28267157 10.1039/c6cp07732c
43. Eone, J. R. A magnetic tight-binding model: Surface properties of transition metals and cobalt nanoparticules. 2022, hal-01949869v3. https://hal.science/hal-01949869v3 (2022).
44. Per H Mats P Long-ranged adsorbate-adsorbate interactions mediated by a surface-state band J. Phys.Condens. Matter. 2000 12 1 L13 10.1088/0953-8984/12/1/103
Per, H. & Mats, P. Long-ranged adsorbate-adsorbate interactions mediated by a surface-state band. J. Phys.Condens. Matter. 12(1), L13. 10.1088/0953-8984/12/1/103 (2000).10.1088/0953-8984/12/1/103
45. Stepanyuk VS Quantum interference and long-range adsorbate-adsorbate interactions Phys. Rev. B 2003 68 20 205410 10.1103/PhysRevB.68.205410
Stepanyuk, V. S. et al. Quantum interference and long-range adsorbate-adsorbate interactions. Phys. Rev. B 68(20), 205410. 10.1103/PhysRevB.68.205410 (2003).10.1103/PhysRevB.68.205410
46. Quaino P Juarez F Santos E Schmickler W Volcano plots in hydrogen electrocatalysis—Uses and abuses Beilstein J. Nanotechnol. 2014 5 846 854 10.3762/bjnano.5.96 24991521
Quaino, P., Juarez, F., Santos, E. & Schmickler, W. Volcano plots in hydrogen electrocatalysis—Uses and abuses. Beilstein J. Nanotechnol. 5, 846–854. 10.3762/bjnano.5.96 (2014).24991521 10.3762/bjnano.5.96
47. Ooka H Huang J Exner KS The Sabatier principle in electrocatalysis: Basics, limitations, extensions Front. Energy Res. 2021 10.3389/fenrg.2021.654460
Ooka, H., Huang, J. & Exner, K. S. The Sabatier principle in electrocatalysis: Basics, limitations, extensions. Front. Energy Res.10.3389/fenrg.2021.654460 (2021).10.3389/fenrg.2021.654460
48. Laplaza R Das S Wodrich MD Corminboeuf C Constructing and interpreting volcano plots and activity maps to navigate homogeneous catalyst landscapes Nat. Protoc. 2022 17 11 2550 2569 10.1038/s41596-022-00726-2 35978038
Laplaza, R., Das, S., Wodrich, M. D. & Corminboeuf, C. Constructing and interpreting volcano plots and activity maps to navigate homogeneous catalyst landscapes. Nat. Protoc. 17(11), 2550–2569. 10.1038/s41596-022-00726-2 (2022).35978038 10.1038/s41596-022-00726-2
49. Laursen AB Electrochemical hydrogen evolution: Sabatier’s principle and the volcano plot J. Chem. Educ. 2012 89 12 1595 1599 10.1021/ed200818t
Laursen, A. B. et al. Electrochemical hydrogen evolution: Sabatier’s principle and the volcano plot. J. Chem. Educ. 89(12), 1595–1599. 10.1021/ed200818t (2012).10.1021/ed200818t
50. Trasatti S Work function, electronegativity, electrochemical behaviour of metals: III. Electrolytic hydrogen evolution in acid solutions J. Electroanal. Chem. Interfac. Electrochem. 1972 39 1 163 184 10.1016/S0022-0728(72)80485-6
Trasatti, S. Work function, electronegativity, electrochemical behaviour of metals: III. Electrolytic hydrogen evolution in acid solutions. J. Electroanal. Chem. Interfac. Electrochem. 39(1), 163–184. 10.1016/S0022-0728(72)80485-6 (1972).10.1016/S0022-0728(72)80485-6
51. Derry GN Kern ME Worth EH Recommended values of clean metal surface work functions J. Vacuum Sci. Technol. A 2015 33 6 060801 10.1116/1.4934685
Derry, G. N., Kern, M. E. & Worth, E. H. Recommended values of clean metal surface work functions. J. Vacuum Sci. Technol. A 33(6), 060801. 10.1116/1.4934685 (2015).10.1116/1.4934685
52. Bronsted JN Acid and basic catalysis Chem. Rev. 1928 5 3 231 338 10.1021/cr60019a001
Bronsted, J. N. Acid and basic catalysis. Chem. Rev. 5(3), 231–338. 10.1021/cr60019a001 (1928).10.1021/cr60019a001
53. Evans MG Polanyi M Inertia and driving force of chemical reactions Trans. Faraday Soc. 1938 10.1039/TF9383400011
Evans, M. G. & Polanyi, M. Inertia and driving force of chemical reactions. Trans. Faraday Soc.10.1039/TF9383400011 (1938).10.1039/TF9383400011
54. Bligaard T Nørskov JK Dahl S Matthiesen J Christensen CH Sehested J The Brønsted–Evans–Polanyi relation and the volcano curve in heterogeneous catalysis J. Catal. 2004 224 1 206 217 10.1016/j.jcat.2004.02.034
Bligaard, T. et al. The Brønsted–Evans–Polanyi relation and the volcano curve in heterogeneous catalysis. J. Catal. 224(1), 206–217. 10.1016/j.jcat.2004.02.034 (2004).10.1016/j.jcat.2004.02.034
55. Nørskov JK Trends in the exchange current for hydrogen evolution J. Electrochem. Soc. 2005 152 3 J23 10.1149/1.1856988
Nørskov, J. K. et al. Trends in the exchange current for hydrogen evolution. J. Electrochem. Soc. 152(3), J23. 10.1149/1.1856988 (2005).10.1149/1.1856988
56. Lousada CM Johansson AJ Korzhavyi PA Cover feature: Adsorption of hydrogen sulfide, hydrosulfide and sulfide at Cu(110)—Polarizability and cooperativity effects first stages of formation of a sulfide layer ChemPhysChem 2018 19 17 2120 2120 10.1002/cphc.201800756
Lousada, C. M., Johansson, A. J. & Korzhavyi, P. A. Cover feature: Adsorption of hydrogen sulfide, hydrosulfide and sulfide at Cu(110)—Polarizability and cooperativity effects first stages of formation of a sulfide layer. ChemPhysChem 19(17), 2120–2120. 10.1002/cphc.201800756 (2018).10.1002/cphc.201800756
57. Rojas KIM Villagracia ARC Narido SC Moreno JLV Arboleda NB First principles study of H2 adsorption on Ni-decorated silicene Mater. Res. Express 2019 6 5 055509 10.1088/2053-1591/ab0294
Rojas, K. I. M., Villagracia, A. R. C., Narido, S. C., Moreno, J. L. V. & Arboleda, N. B. First principles study of H2 adsorption on Ni-decorated silicene. Mater. Res. Express 6(5), 055509. 10.1088/2053-1591/ab0294 (2019).10.1088/2053-1591/ab0294
58. Watkins WL Borensztein Y Mechanism of hydrogen adsorption on gold nanoparticles and charge transfer probed by anisotropic surface plasmon resonance Phys. Chem. Chem. Phys. 2017 19 40 27397 27405 10.1039/C7CP04843B 28972603
Watkins, W. L. & Borensztein, Y. Mechanism of hydrogen adsorption on gold nanoparticles and charge transfer probed by anisotropic surface plasmon resonance. Phys. Chem. Chem. Phys. 19(40), 27397–27405. 10.1039/C7CP04843B (2017).28972603 10.1039/C7CP04843B
59. Kayode GO Zhang S Montemore MM Spivey J Han Y-F Shekhawat D Linking electronic structure to adsorption energies: metal surfaces and single-atom catalysts Catalysis 2022 The Royal Society of Chemistry 34
Kayode, G. O., Zhang, S. & Montemore, M. M. Linking electronic structure to adsorption energies: metal surfaces and single-atom catalysts. In Catalysis (eds Spivey, J. et al.) 34 (The Royal Society of Chemistry, 2022).
60. Schumann J Stamatakis M Michaelides A Réocreux R Ten-electron count rule for the binding of adsorbates on single-atom alloy catalysts Nat. Chem. 2024 16 5 749 754 10.1038/s41557-023-01424-6 38263384
Schumann, J., Stamatakis, M., Michaelides, A. & Réocreux, R. Ten-electron count rule for the binding of adsorbates on single-atom alloy catalysts. Nat. Chem. 16(5), 749–754. 10.1038/s41557-023-01424-6 (2024).38263384 10.1038/s41557-023-01424-6
61. Venkataramanan NS Suvitha A Mizuseki H Kawazoe Y A theoretical study of the effects of transition metal dopants on the adsorption and dissociation of hydrogen on nickel clusters Int. J. Quantum Chem. 2013 113 15 1940 1948 10.1002/qua.24418
Venkataramanan, N. S., Suvitha, A., Mizuseki, H. & Kawazoe, Y. A theoretical study of the effects of transition metal dopants on the adsorption and dissociation of hydrogen on nickel clusters. Int. J. Quantum Chem. 113(15), 1940–1948. 10.1002/qua.24418 (2013).10.1002/qua.24418
62. Hammer B Nørskov JK Electronic factors determining the reactivity of metal surfaces Surf. Sci. 1995 343 3 211 220 10.1016/0039-6028(96)80007-0
Hammer, B. & Nørskov, J. K. Electronic factors determining the reactivity of metal surfaces. Surf. Sci. 343(3), 211–220. 10.1016/0039-6028(96)80007-0 (1995).10.1016/0039-6028(96)80007-0
63. Wang G-C Li J Xu X-F Li R-F Nakamura J The relationship between adsorption energies of methyl on metals and the metallic electronic properties: A first-principles DFT study J. Comput. Chem. 2005 26 9 871 878 10.1002/jcc.20225 15825187
Wang, G.-C., Li, J., Xu, X.-F., Li, R.-F. & Nakamura, J. The relationship between adsorption energies of methyl on metals and the metallic electronic properties: A first-principles DFT study. J. Comput. Chem. 26(9), 871–878. 10.1002/jcc.20225 (2005).15825187 10.1002/jcc.20225
64. Bligaard T Nørskov JK Nilsson A Pettersson LGM Nørskov JK Chapter 4—heterogeneous catalysis Chemical Bonding at Surfaces and Interfaces 2008 Elsevier 255 321
Bligaard, T. & Nørskov, J. K. Chapter 4—heterogeneous catalysis. In Chemical Bonding at Surfaces and Interfaces (eds Nilsson, A. et al.) 255–321 (Elsevier, 2008).
65. Nilsson A Pettersson LGM Nilsson A Pettersson LGM Nørskov JK Chapter 2—Adsorbate electronic structure and bonding on metal surfaces Chemical Bonding at Surfaces and Interfaces 2008 Elsevier 57 142
Nilsson, A. & Pettersson, L. G. M. Chapter 2—Adsorbate electronic structure and bonding on metal surfaces. In Chemical Bonding at Surfaces and Interfaces (eds Nilsson, A. et al.) 57–142 (Elsevier, 2008).
66. Geerlings P De Proft F Langenaeker W Conceptual density functional theory Chem Rev. 2003 103 5 1793 1874 10.1021/cr990029p 12744694
Geerlings, P., De Proft, F. & Langenaeker, W. Conceptual density functional theory. Chem Rev. 103(5), 1793–1874. 10.1021/cr990029p (2003).12744694 10.1021/cr990029p
