This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Predictive model of surface adsorption in dissolution on transition metals and alloys

Bo Li    Wang Gao wgao@jlu.edu.cn    Qing Jiang Key Laboratory of Automobile Materials (Jilin University), Ministry of Education, Department of Materials Science and Engineering, Jilin University, Changchun 130022, China
Abstract

Surface adsorption, which is often coupled with surface dissolution, is generally unpredictable on alloys due to the complicated alloying and dissolution effects. Herein, we introduce the electronic gradient and cohesive properties of surface sites to characterize the effects of alloying and dissolution. This enables us to build a predictive model for the quantitative determination of the adsorption energy in dissolution, which holds well for transition metals, near-surface alloys, binary alloys, and high-entropy alloys. Furthermore, this model uncovers a synergistic mechanism between the dd-band upper-edge ratio, dd-band width and ss-band depth in determining the alloying and dissolution effects on adsorption. Our study not only provides fundamental mechanistic insights into surface adsorption on alloys but also offers a long-sought tool for the design of advanced alloy catalysts.

pacs:
75.80.+q, 77.65.-j
preprint: Physical Review Letters

Alloys, covering from near-surface alloys (NSAs), binary alloys (BAs), to high-entropy alloys (HEAs), are essential for heterogeneous catalysis and anticorrosion engineering Greeley and Mavrikakis (2004); Stamenkovic et al. (2007); Yeh et al. (2004); Ma et al. (2017). Their catalytic activity, which is governed by the adsorption energy in terms of the Sabatier principle Sabatier (1913) and Brønsted–Evans–Polanyi (BEP) relation Brønsted (1928); Evans and Polanyi (1938); Bligaard et al. (2004), is the footing stone for the development of electrode materials. However, when the adsorbates are continually adsorbed to the surfaces of metallics during the catalytic process, the surfaces also suffer severe corrosion under electrochemical condition with the surface-atoms dissolution Wang et al. (2011); Topalov et al. (2012); Kang et al. (2014); Lopes et al. (2016, 2020). Namely, the adsorption of surface sites is generally competitive and coupled with surface dissolution, prohibiting the establishment of principles for guiding material design.

The main challenge stems from the alloying and dissolution effects, since alloying and dissolution can generate a diverse coordination environment of surface sites. In particular, the random composition and chemical ordering of HEAs offer multiple active centers, possess significant charge redistribution, and create different crystal phases, all of which bring more opportunities to modulate the electronic and geometric structures of surface, unfortunately also making it more difficult to explicitly understand the picture of alloying effects Löffler et al. (2018); Xie et al. (2019); Batchelor et al. (2019); Lu et al. (2020); Xin et al. (2020); Yao et al. (2020). Therefore, the current understanding of adsorption achieves mainly on simple metallics instead of complex alloys Hammer and Nørskov (1995a, b); Hammer and Nørskov (2000); Abild-Pedersen et al. (2007); Vojvodic et al. (2009, 2014); Xin et al. (2014); İnoğlu and Kitchin (2010); Mpourmpakis et al. (2010); Calle-Vallejo et al. (2014, 2015); Ma and Xin (2017); Calle-Vallejo et al. (2012); Gao et al. (2020); Li et al. (2021). The adsorption energy has been correlated to the d-band center or upper edge in d-band model Hammer and Nørskov (1995a, b); Hammer and Nørskov (2000); Vojvodic et al. (2009, 2014); Xin et al. (2014); İnoğlu and Kitchin (2010), coordination number Mpourmpakis et al. (2010); Calle-Vallejo et al. (2014, 2015); Ma and Xin (2017), valence and electronegativity of surface atoms Calle-Vallejo et al. (2012); Gao et al. (2020); Li et al. (2021), which hold for transition metals (TMs), TM nanoparticles (NPs), and partial NSAs. However, it still presents a fundamental challenge to uncover the mechanism and do the quantitative prediction for the adsorption on various alloys (e.g. BAs and HEAs), especially in dissolution.

Herein we introduce the electronic gradient and cohesive properties of surface sites to characterize the effects of alloying and dissolution. Accordingly, a predictive model has been established to accurately determine the adsorption energy in dissolution on TMs, NSAs, BAs, and HEAs. This model further unravels the electronic origin of metallics in determining the adsorption energy, which outlines a novel physical picture for understanding surface adsorption on alloys.

To ensure the generality of our scheme, we have considered the widespread adsorption of CHx\rm{CH}_{\textit{x}} (x = 0-3), CCH3\rm{CCH_{3}}, CO, NHx\rm{NH}_{\textit{x}} (x = 0-2), OH, F and Cl on (100), (110), (111), (211) and (532) surfaces of TMs, NSAs, BAs, and HEAs (Fig. S1 li2 (2021)). The adsorption process in dissolution is illustrated in Note S1, the inset of Fig. 1a and b, and Figs. S2-S4 li2 (2021). All calculations were performed with Perdew-Burke-Ernzerhof (PBE) Perdew et al. (1996) and Perdew-Wang-91 (PW91) Perdew et al. (1992) exchange-correlation (XC) functionals (see more details in Note S2 li2 (2021)).

Refer to caption
Figure 1: Comparison between the electronic descriptor ψ0\psi_{0} and cohesive energy EcohE_{\rm{coh}} Zhang (2014) of adsorption-site atoms in describing the top-site adsorption energy in dissolution. (a) and (b) (111) surface of TMs. (c)-(f) (111) surface of Cu-based NSAs, (100) surface of Ag-based BAs, and (111) surface of RuRhIrPdPt-based HEAs. The symbol numbers \small{1}⃝-\small{4}⃝ of the inset of (a) and (b) denote the adsorption order in dissolution (see details in Figs. S2-S4 and Note S1 li2 (2021)). Taking HEAs as example, the red circle and black triangle of the inset of (c) denote the adsorption site and its nearest neighbors. The violet square and cyan rhombus icons for HEAs in (c)-(f) represent the cases that changing only the element of adsorption sites and that changing all the element that is the same as the adsorption site.

To understand the alloying effects on adsorption, we study the active center of adsorption, which includes the adsorption site and its nearest neighbors (see the inset of Fig. 1c and Fig. S1 li2 (2021)). The alloying effects can be roughly separated into the adsorption-site effect and the environment effect. The former refers to the change of chemical composition of adsorption sites with the fixed surrounding properties, whereas the latter does the change of the chemical properties of the nearest neighbors with the fixed adsorption-site properties. We adopt (i=1N𝑆vi)2/N/(i=1N𝜒i)1/N{\mathop{{{\left({{\mathop{\prod}\nolimits^{{N}}_{{i=1}}{\mathop{{S}}\nolimits_{{\rm{v}\textit{i}}}}}}\right)}}}\nolimits^{{2/N}}/\mathop{{{\left({{\mathop{\prod}\nolimits^{{N}}_{{i=1}}{\mathop{{\chi}}\nolimits_{{i}}}}}\right)}}}\nolimits^{{1/N}}} to calculate the electronic descriptor for the adsorption site ψ0\psi_{{0}} and for the active center ψ1\psi_{{1}}, which correspond to the local and non-local electronic effects respectively. NN is the atom number of the adsorption site or the active center, while SviS_{\rm{v}\textit{i}} and χi\chi_{i} represent the valence-electron number (including that of both dd- and ss-electrons) and Pauling electronegativity of i{i}th involved atom. Hence, we introduce the electronic gradient descriptor of surface sites as,

ψ=ψ0m/ψ1n\displaystyle\psi^{{}^{\prime}}=\psi_{0}^{m}/\psi_{1}^{n} (1)

With the change of m and n, ψ\psi^{{}^{\prime}} provides a simple way to estimate the gradient of surface electronic properties of alloys. Notably, ψ\psi^{{}^{\prime}} can (partially) characterize the adsorption-site effect or environment effect of alloying, or the coupling of them depending on the surface sites as well as m and n.

Now we try to understand the adsorption during the dissolution process. Fig. 1a and b and Fig. S5c and d li2 (2021) show the adsorption energy of \chCH2 on TM(111) surface in dissolution. Clearly, the adsorption-site descriptor ψ0\psi_{0} performs well in characterizing the trends of adsorption energy on the undissolved surfaces but gradually loses efficiency in the characterization as the dissolution proceeds. In contrast, the cohesive energy of adsorption-site atoms EcohE_{\rm{coh}} Zhang (2014) becomes increasingly accurate in describing the adsorption energy with the dissolution process going on. Based on our findings on the undissolved surfaces of TMs Gao et al. (2020), we propose an expression of the adsorption energy in dissolution as,

Ead=μ1×Dad+μ2×CN¯+θ=110×XmXXm+1×𝐷ad+15×X+1𝑋m+1×CN¯+θ\displaystyle\begin{aligned} E_{\rm{ad}}&=\mu_{1}\times{D_{\rm{ad}}}+\mu_{2}\times\overline{CN}+\theta\\ &=\frac{{1}}{{10}}\times\frac{{X_{\rm{m}}}-X}{{X_{\rm{m}}}+1}\times\mathop{{D}}\nolimits_{{\rm{ad}}}+\frac{{1}}{{5}}\times\frac{{X+1}}{{\mathop{{X}}\nolimits_{{\rm{m}}}+1}}\times\overline{CN}+\theta\end{aligned} (2)

where XmX_{\rm{m}} and XX correspond to the maximum bondable and actual bonding numbers of the central atom for a given adsorbate. The generalized coordination number CN¯\overline{CN} is obtained by dividing the summation of the usual coordination number CNCN of the adsorption sites nearest neighbors with the CNCN in bulk Calle-Vallejo et al. (2014, 2015). θ\theta is a constant for a given adsorbate and is likely controlled by the coupling between the sp states of substrates and the states of adsorbates on TMs (θ\theta can be easily determined by only doing calculations on one substrate) Gao et al. (2020). The electronic descriptor Dad{\rm\textit{D}_{ad}} for the adsorption energy on TMs in dissolution can be expressed as,

Dad=k1ψ+k2(252𝐸coh+115)=(CN¯+1CNu+1)2ψ0+[1(CN¯+1CNu+1)2](252Ecoh+115)\displaystyle\begin{aligned} &{{D_{\rm{ad}}=k_{1}\psi^{{}^{\prime}}+k_{2}}}{\left({\frac{{25}}{{2}}\mathop{{E}}\nolimits_{{\rm{coh}}}+115}\right)}\\ &={\left(\frac{\overline{CN}+1}{CN_{\rm{u}}+1}\right)}^{2}\psi_{0}+\left[{1-{\left(\frac{\overline{CN}+1}{CN_{\rm{u}}+1}\right)}^{2}}\right]\left(\frac{25}{2}E_{\rm{coh}}+115\right)\end{aligned} (3)
Refer to caption
Figure 2: Adsorption energies of the different adsorbates as a function of the proposed descriptors in dissolution. (a) and (b) CHx\rm{CH}_{\textit{x}} (x = 1-3), CO, NHx\rm{NH}_{\textit{x}} (x = 1, 2), and OH on (111) surface of TMs Roling and Abild-Pedersen (2018). (c) CH2\rm{CH}_{2} on (100) and (111) surfaces of Ag-based BAs in the adsorption-site effect of alloying (blue triangle) and the environment effect of alloying (cyan rhombus), and (111) surface of Cu-based NSAs in the adsorption-site effect of alloying. (d) CH2\rm{CH}_{2} and CO on (100), (110), (111), (211) and (532) surfaces of RuRhIrPdPt-based HEAs. The blue triangle and violet square represent the adsorption-site effect of alloying, while the rest do the environment effect of alloying. (e) CH, NH and OH on (111) surface of RuRhIrPdPt-based HEAs in the adsorption-site effect of alloying. The data in (a)-(e) are obtained at the top site. (f) CH at the hcp site and CO at the bridge site on the undissolved (211) surface of AgAu, AgPd, IrRu, and PtRh BAs Andersen et al. (2019), and OH at the bridge site of (111) surface of RuRhIrPdPt-based HEAs in both the adsorption-site and environment effects of alloying.
Table 1: Comparison between the predicted and DFT-calculated prefactors μ1=1/10×(XmX)/(Xm+1)\mu_{1}=1/10\times(X_{\rm{m}}-X)/(X_{\rm{m}}+1) and μ2=1/5×(X+1)/(Xm+1)\mu_{2}=1/5\times(X+1)/(X_{\rm{m}}+1) of Eq. (2) on TMs and HEAs. The calculated values for CH2\rm{CH_{2}} adsorption are μ1\mu_{1} = 0.035 and μ2\mu_{2} = 0.099 on NSAs and μ1\mu_{1} = 0.040 and μ2\mu_{2} = 0.120 on BAs. All data are obtained at the top site.
species XmX_{\rm{m}} XX Predicted DFT-calculated
TMs HEAs
μ1\mu_{1} μ2\mu_{2} μ1\mu_{1} μ2\mu_{2} μ1\mu_{1} μ2\mu_{2}
CH 4 1 0.060 0.080 0.062 0.081 0.054 0.080
CH2\rm{CH_{2}} 4 2 0.040 0.120 0.039 0.110 0.035 0.125
CH3\rm{CH_{3}} 4 3 0.020 0.160 0.016 0.115 0.015 0.120
CO 4 2 0.040 0.120 0.033 0.114 0.033 0.122
CCH3\rm{CCH_{3}} 4 1 0.060 0.080 0.064 0.079
NH 3 1 0.050 0.100 0.046 0.104 0.038 0.098
NH2\rm{NH_{2}} 3 2 0.025 0.150 0.024 0.148 0.025 0.145
OH 2 1 0.033 0.134 0.021 0.156 0.022 0.135

where CNu{\rm\textit{CN}_{u}} is the usual coordination number of adsorption sites on the undissolved surfaces. The coefficients 25/2 and 115 of the EcohE_{\rm{coh}} term stem from the slope and offset of the fitting relation between ψ0\psi_{{{0}}} and EcohE_{\rm{coh}} (Fig. S6 li2 (2021)). As the dissolution proceeds, there exists a gradual transition from ψ0\psi_{0} to EcohE_{\rm{coh}} for the electronic determinant of adsorption energy on TMs (Figs. S5, S7-S11, S12a and S12b li2 (2021)). Eqs. (2) and (3) describe well the adsorption energy of CHx\rm{CH}_{\textit{x}} (x = 1-3), CCH3\rm{CCH_{3}}, CO, NHx\rm{NH}_{\textit{x}} (x = 1, 2), and OH at the different adsorption coverage and sites on TMs in dissolution (Fig. 2a and b and Figs. S13-S15 and Note S3 li2 (2021)) Roling and Abild-Pedersen (2018); Andersen et al. (2019); Hoffmann et al. (2016), with the fitted mean absolute error (MAE) of 0.16 eV. The fitted prefactors μ1\mu_{1} and μ2\mu_{2} are also in good agreement with the predictions by Eq. (2) (Tables I and S1 li2 (2021)). Note that the prefactors μ1\mu_{1} and μ2\mu_{2} of Eq. (2), determined by the bonding characters of adsorbates (XmX)/(Xm+1)(X_{\rm{m}}-X)/(X_{\rm{m}}+1) and (X+1)/(Xm+1)(X+1)/(X_{\rm{m}}+1), can be deduced and rationalized from the first-order approximation of effective medium theory (EMT) and bond-order conservation framework (Note S4 li2 (2021)) Nørskov and Lang (1980); Gao et al. (2020).

We now quantify the adsorption-site effect of alloying on the adsorption energy during dissolution. Fig. 1c-f and Figs. S12c, S12d and S16-S19 li2 (2021) show that EcohE_{\rm{coh}} of adsorption-site atoms generally performs better than the adsorption-site descriptor ψ0\psi_{0} in describing the adsorption energy of CHx\rm{CH}_{\textit{x}} (x = 1-3), CO, NHx\rm{NH}_{\textit{x}} (x = 1, 2), and OH at the top site of NSAs, BAs, and HEAs in the adsorption-site effect of alloying for each dissolution step. For HEAs, changing only the element of adsorption sites or changing all the element that is the same as the adsorption site follows the same linear relation (Fig. 1c-f). By further including the coordination-chemistry effect with the electronic gradient descriptor ψ\psi^{{}^{\prime}}, we introduce a universal descriptor to quantitatively describe the electronic contribution of alloys to the adsorption energy in dissolution as,

Dad=k1ψ+k2(252Ecoh+115)=k1ψ02ψ1+k2(252Ecoh+115)\displaystyle\begin{aligned} D_{\rm{ad}}&=k_{1}\psi^{{}^{\prime}}+k_{2}\left(\frac{25}{2}E_{\rm{coh}}+115\right)\\ &=k_{1}\frac{\psi_{0}^{2}}{\psi_{1}}+k_{2}\left(\frac{25}{2}E_{\rm{coh}}+115\right)\end{aligned} (4)
Refer to caption
Figure 3: (a) The d-band center (εd\varepsilon_{\rm d}) and d-band upper edge (εu\varepsilon_{\rm u}) against the electronic descriptor ψ0\psi_{0} for atoms on the undissolved (111) surface of TMs Vojvodic et al. (2014) and RuRhIrPdPt-based HEAs. (b) The d-band width (WdW_{\rm d}) for surface sites on the undissolved (111) surface (\small{1}⃝) and that with three surface atoms dissolved (\small{4}⃝), against the cohesive energy Zhang (2014). (c) Density of states of the s band of Au atoms on the undissolved (111) surface and that with three atoms dissolved. (d) Comparison between the predicted and DFT-calculated adsorption energies of CHx\rm{CH}_{\textit{x}} (x = 0-3), CCH3\rm{CCH_{3}}, CO, NHx\rm{NH}_{\textit{x}} (x = 0-2), OH, F and Cl on TMs, NSAs, BAs, and HEAs Lu et al. (2020); Roling and Abild-Pedersen (2018); Andersen et al. (2019); Hoffmann et al. (2016); Ma et al. (2015); Calle-Vallejo et al. (2012); Xin et al. (2012a, b); Mamun et al. (2019); Ghosh et al. (2017); Gan et al. (2012).

where ψ02/ψ1{\psi_{0}^{2}}/\psi_{1} automatically transforms to ψ0\psi_{0} for TMs, implying that our electronic gradient descriptor of adsorption energy is universal for TMs and alloys. Since the EcohE_{\rm{coh}} term corresponds to the effect of adsorption-site atoms, namely the local electronic effect, while the ψ\psi^{{}^{\prime}} term considers the coupling between the adsorption site and its surrounding environments, the prefactors k1k_{1} and k2k_{2} can effectively reflect the electronic localization of alloys in adsorption. We find that k1k_{1} = 1/10 and k2k_{2} = 9/10 are the optimum values for describing the adsorption-site effect of NSAs, BAs, and HEAs on the adsorption energy. Accordingly, Eqs. (2) and (4) can describe well the adsorption energy of CHx\rm{CH}_{\textit{x}} (x = 1-3), CO, NHx\rm{NH}_{\textit{x}} (x = 1, 2), and OH at the top site of NSAs, BAs, and HEAs (Fig. 2c-e and Figs. S14b, S20 and S21 li2 (2021)), where the fitted MAEs are merely 0.18 eV for NSAs, 0.12 eV for BAs, and 0.14 eV for HEAs. Moreover, the fitted prefactors μ1\mu_{1} and μ2\mu_{2} are also consistent with the predictions by Eq. (2) (Tables I, S2 and S3 li2 (2021)). Namely, the prefactors μ1\mu_{1} and μ2\mu_{2} of alloys are also determined by the bonding characters XmX_{\rm{m}} and XX as those of TMs do, which comply with the first-order approximation of EMT Nørskov and Lang (1980) and bond-order conservation framework, indicating that the homogeneous-electron-gas model, EMT, is suitable for describing the adsorption behavior on both TMs and alloys (Note S4 li2 (2021)).

We now turn to understand the environment effect of alloying on the adsorption energy. Fig. 2c and d and Fig. S22 li2 (2021) show that the adsorption energies of CHx\rm{CH}_{\textit{x}} (x = 0-3), CO, NHx\rm{NH}_{\textit{x}} (x = 0-2), OH, F and Cl at the top site of Ag-based BAs, Pt(Pd)-based NSAs Ma et al. (2015); Calle-Vallejo et al. (2012); Xin et al. (2012a, b) and Ru(Cu)RhIrPdPt-based HEAs in the environment effect of alloying also follow the linear relations in the adsorption-site effect of alloying with the fitted MAE 0.06 eV for NSAs, BAs, and HEAs, which further demonstrates the robustness of our scheme. Notably, the prefactors k1k_{1} and k2k_{2} in Eq. (4) are different upon adsorbing on the different elements of HEAs in the environment effect of alloying, indicating the distinct localization of the different alloying elements in determining the adsorption energy. In particular, the localization of RuRhIrPdPt HEAs in adsorbing CH2\rm{CH}_{2} obeys the order of Ru >> Ir >> Pt >> Pd >> Rh (Note S5 and Table S4 li2 (2021)).

When both the adsorption site and surrounding environments change simultaneously from one system to the next, Eqs. (2) and (4) can also capture these complicated alloying effects. For instance, the adsorption of OH at the bridge site of HEAs, the adsorption of CH, CO, C and OH at the top, bridge and hollow sites of BAs and the adsorption of OH, H and OH+H at the hollow site of NSAs comply with our model by Eqs. (2) and (4), with the fitted MAE of 0.05 eV for NSAs, 0.17 eV for BAs, and 0.08 eV for HEAs (Fig. 2f and Figs. S23-S26, Tables S2-S6 and Note S6 li2 (2021)) Lu et al. (2020); Andersen et al. (2019); Mamun et al. (2019); Gan et al. (2012); Ghosh et al. (2017). Overall, our scheme is universal for capturing the adsorption energy on TMs and alloys in dissolution.

Our scheme reveals a novel physical picture for surface adsorption. The d-band model by Nø\orskov et al. Hammer and Nørskov (1995b, a); Hammer and Nørskov (2000); Vojvodic et al. (2014), which takes the d-band center or d-band upper edge as the descriptors of adsorption and assumes a constant contribution of s bands to adsorption for metals, is mainly applied into TMs but not complex alloys like BAs and HEAs Andersen et al. (2019). Our scheme, which is effective in complex alloys, can also be related to the d-band and s-band properties of TMs but in a significantly different way from the d-band model. ψ\psi^{{}^{\prime}} reflects the atoms d-band upper-edge ratio between the adsorption site and its nearest-neighbors (Fig. 3a), whereas the adsorption-site atoms cohesive energy does their own d-band width and s-band depth Turchanin and Agraval (2008). Namely, the d-band upper-edge ratio of active centers and the d-band width and s-band depth of adsorption sites together control the alloying effects in adsorption. The dissolution of a surface atom reduces significantly its neighbors d-band width and s-band depth (Fig. 3b and c and Fig. S27 li2 (2021)), making their cohesive energy crucial to determining the adsorption energy on metallics in dissolution. Notably, the exclusion of the s-band contribution makes the model lose its prediction ability by increasing the MAE to >> 0.32 eV (Fig. S28 and Table S7 li2 (2021)), demonstrating a material-dependent contribution of s bands of adsorption sites to the adsorption of complex alloys. Clearly, our scheme goes beyond the d-band model in understanding the adsorption on TMs and alloys in dissolution.

Now we access the prediction accuracy of our model. Encouragingly, our model, Eqs. (2-4), is convenient to predict the adsorption energy on metallics as all parameters are easily accessible. The MAE of the predicted adsorption energies relative to the DFT results is \sim0.14 eV for the adsorption of CHx\rm{CH}_{\textit{x}} (x = 0-3), CCH3\rm{CCH_{3}}, CO, NHx\rm{NH}_{\textit{x}} (x = 0-2), OH, F and Cl on TMs, NSAs, BAs, and HEAs (Fig. 3d) regardless of the adopted (semi-)local functionals (Note S6 li2 (2021)), which is less than the approximate error of (semi-)local functionals, ±0.2 eV. This model is also predictive for the reactivity at surfaces including the activation energy, reaction energy, and experimentally measured activity for 11 different reactions, and even including those on HEAs that can break the BEP relation (Note S7 and Figs. S29-S40 li2 (2021)). Considering the large amount (2586 adsorption energies and 539 reactivity data) and wide window (11.16 eV for adsorption energy and 6.30 eV for reactivity) of data as well as the easy accessibility of the descriptors, our model exhibits promising prediction power for surface adsorption and reaction on TMs and alloys.

In summary, we have identified the electronic gradient and cohesive energy of surface sites as new powerful descriptors for understanding surface adsorption and reaction on metallics in dissolution. With these descriptors, a predictive model has been established for the quantitative determination of the adsorption energy on TMs, NSAs, BAs, and HEAs. Our scheme is a crucial development of d-band model, revealing that the s-band contribution is material-dependent and indispensable in the accurate description of the alloying and dissolution effects on surface adsorption and reaction. The electronic gradient descriptor ψ02/ψ1{\psi_{0}^{2}}/\psi_{1}, corresponding to the d-band upper-edge ratio of active centers, is crucial for determining the adsorption energy in both the adsorption-site and environment effects of alloying, whereas the cohesive-energy descriptor, corresponding to the d-band width and s-band depth of adsorption sites, plays an increasingly important role in determining the adsorption energy as the alloying and dissolution proceed and is dominant in the adsorption-site effect of alloying. This scheme unravels the fundamental physics of surface adsorption on TMs and alloys, and provides a convenient predictive tool for the rational design of advanced alloy catalysts.

We gratefully acknowledge support from the National Natural Science Foundation of China (Nos. 22173034, 21673095, 11974128, 51631004, 52130101), the Opening Project of State Key Laboratory of High Performance Ceramics and Superfine Microstructure (SKL201910SIC), the Program of Innovative Research Team (in Science and Technology) in University of Jilin Province, the Program for JLU (Jilin University) Science and Technology Innovative Research Team (No. 2017TD-09), the Fundamental Research Funds for the Central Universities, and the computing resources of the High Performance Computing Center of Jilin University, China.

References

  • Greeley and Mavrikakis (2004) J. Greeley and M. Mavrikakis, Nat. Mater. 3, 810 (2004).
  • Stamenkovic et al. (2007) V. R. Stamenkovic, B. S. Mun, M. Arenz, K. J. Mayrhofer, C. A. Lucas, G. Wang, P. N. Ross, and N. M. Markovic, Nat. Mater. 6, 241 (2007).
  • Yeh et al. (2004) J.-W. Yeh, S.-K. Chen, S.-J. Lin, J.-Y. Gan, T.-S. Chin, T.-T. Shun, C.-H. Tsau, and S.-Y. Chang, Adv. Energy Mater. 6, 299 (2004).
  • Ma et al. (2017) H. Ma, X.-Q. Chen, R. Li, S. Wang, J. Dong, and W. Ke, Acta Mater. 130, 137 (2017).
  • Sabatier (1913) P. Sabatier, Ch. Beranger, Editeur. Paris p. 180 (1913).
  • Brønsted (1928) J. Brønsted, Chem. Rev. 5, 231 (1928).
  • Evans and Polanyi (1938) M. Evans and M. Polanyi, Trans. Faraday Soc. 34, 11 (1938).
  • Bligaard et al. (2004) T. Bligaard, J. K. Nørskov, S. Dahl, J. Matthiesen, C. H. Christensen, and J. Sehested, J. Catal. 224, 206 (2004).
  • Wang et al. (2011) C. Wang, D. Van Der Vliet, K. L. More, N. J. Zaluzec, S. Peng, S. Sun, H. Daimon, G. Wang, J. Greeley, J. Pearson, et al., Nano Lett. 11, 919 (2011).
  • Topalov et al. (2012) A. A. Topalov, I. Katsounaros, M. Auinger, S. Cherevko, J. C. Meier, S. O. Klemm, and K. J. Mayrhofer, Angew. Chem. Int. Ed. 51, 12613 (2012).
  • Kang et al. (2014) Y. Kang, J. Snyder, M. Chi, D. Li, K. L. More, N. M. Markovic, and V. R. Stamenkovic, Nano Lett. 14, 6361 (2014).
  • Lopes et al. (2016) P. P. Lopes, D. Strmcnik, D. Tripkovic, J. G. Connell, V. Stamenkovic, and N. M. Markovic, Acs Catal. 6, 2536 (2016).
  • Lopes et al. (2020) P. P. Lopes, D. Li, H. Lv, C. Wang, D. Tripkovic, Y. Zhu, R. Schimmenti, H. Daimon, Y. Kang, J. Snyder, et al., Nat. Mater. 19, 1207 (2020).
  • Löffler et al. (2018) T. Löffler, H. Meyer, A. Savan, P. Wilde, A. Garzón Manjón, Y.-T. Chen, E. Ventosa, C. Scheu, A. Ludwig, and W. Schuhmann, Adv. Energy Mater. 8, 1802269 (2018).
  • Xie et al. (2019) P. Xie, Y. Yao, Z. Huang, Z. Liu, J. Zhang, T. Li, G. Wang, R. Shahbazian-Yassar, L. Hu, and C. Wang, Nat. Commun. 10, 4011 (2019).
  • Batchelor et al. (2019) T. A. Batchelor, J. K. Pedersen, S. H. Winther, I. E. Castelli, K. W. Jacobsen, and J. Rossmeisl, Joule 3, 834 (2019).
  • Lu et al. (2020) Z. Lu, Z. W. Chen, and C. V. Singh, Matter 3, 1318 (2020).
  • Xin et al. (2020) Y. Xin, S. Li, Y. Qian, W. Zhu, H. Yuan, P. Jiang, R. Guo, and L. Wang, ACS Catal. 10, 11280 (2020).
  • Yao et al. (2020) Y. Yao, Z. Liu, P. Xie, Z. Huang, T. Li, D. Morris, Z. Finfrock, J. Zhou, M. Jiao, J. Gao, et al., Sci. Adv. 6, eaaz0510 (2020).
  • Hammer and Nørskov (1995a) B. Hammer and J. K. Nørskov, Surf. Sci. 343, 211 (1995a).
  • Hammer and Nørskov (1995b) B. Hammer and J. K. Nørskov, Nature 376, 238 (1995b).
  • Hammer and Nørskov (2000) B. Hammer and J. K. Nørskov, Adv. Catal. 45, 71 (2000).
  • Abild-Pedersen et al. (2007) F. Abild-Pedersen, J. Greeley, F. Studt, J. Rossmeisl, T. R. Munter, P. G. Moses, E. Skulason, T. Bligaard, and J. K. Nørskov, Phys. Rev. Lett. 99, 016105 (2007).
  • Vojvodic et al. (2009) A. Vojvodic, A. Hellman, C. Ruberto, and B. I. Lundqvist, Phys. Rev. Lett. 103, 146103 (2009).
  • Vojvodic et al. (2014) A. Vojvodic, J. K. Nørskov, and F. Abild-Pedersen, Top. Catal. 57, 25 (2014).
  • Xin et al. (2014) H. Xin, A. Vojvodic, J. Voss, J. K. Nørskov, and F. Abild-Pedersen, Phys. Rev. B 89, 115114 (2014).
  • İnoğlu and Kitchin (2010) N. İnoğlu and J. R. Kitchin, Phys. Rev. B 82, 045414 (2010).
  • Mpourmpakis et al. (2010) G. Mpourmpakis, A. N. Andriotis, and D. G. Vlachos, Nano Lett. 10, 1041 (2010).
  • Calle-Vallejo et al. (2014) F. Calle-Vallejo, J. I. Martínez, J. M. García-Lastra, P. Sautet, and D. Loffreda, Angew. Chem. Int. Ed. 53, 8316 (2014).
  • Calle-Vallejo et al. (2015) F. Calle-Vallejo, J. Tymoczko, V. Colic, Q. H. Vu, M. D. Pohl, K. Morgenstern, D. Loffreda, P. Sautet, W. Schuhmann, and A. S. Bandarenka, Science 350, 185 (2015).
  • Ma and Xin (2017) X. Ma and H. Xin, Phys. Rev. Lett. 118, 036101 (2017).
  • Calle-Vallejo et al. (2012) F. Calle-Vallejo, J. I. Martínez, J. M. García-Lastra, J. Rossmeisl, and M. T. M. Koper, Phys. Rev. Lett. 108, 116103 (2012).
  • Gao et al. (2020) W. Gao, Y. Chen, B. Li, S.-P. Liu, X. Liu, and Q. Jiang, Nat. Commun. 11, 1196 (2020).
  • Li et al. (2021) B. Li, W. Gao, and Q. Jiang, J. Phys. Energy 3, 022001 (2021).
  • li2 (2021) Supplemental Materials for ”Predictive model of surface adsorption in dissolution on transition metals and alloys” (2021).
  • Perdew et al. (1996) J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
  • Perdew et al. (1992) J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, and C. Fiolhais, Phys. Rev. B 46, 6671 (1992).
  • Zhang (2014) G.-X. Zhang, Ph.D. thesis, Freie Universität Berlin (2014).
  • Roling and Abild-Pedersen (2018) L. T. Roling and F. Abild-Pedersen, ChemCatChem 10, 1643 (2018).
  • Andersen et al. (2019) M. Andersen, S. V. Levchenko, M. Scheffler, and K. Reuter, ACS Catal. 9, 2752 (2019).
  • Hoffmann et al. (2016) M. J. Hoffmann, A. J. Medford, and T. Bligaard, J. Phys. Chem. C 120, 13087 (2016).
  • Nørskov and Lang (1980) J. K. Nørskov and N. D. Lang, Phys. Rev. B 21, 2131 (1980).
  • Ma et al. (2015) X. Ma, Z. Li, L. E. Achenie, and H. Xin, J. Phys. Chem. Lett. 6, 3528 (2015).
  • Xin et al. (2012a) H. Xin, A. Holewinski, N. Schweitzer, E. Nikolla, and S. Linic, Top. Catal. 55, 376 (2012a).
  • Xin et al. (2012b) H. Xin, A. Holewinski, and S. Linic, ACS Catal. 2, 12 (2012b).
  • Mamun et al. (2019) O. Mamun, K. T. Winther, J. R. Boes, and T. Bligaard, Sci. Data 6, 76 (2019).
  • Ghosh et al. (2017) S. Ghosh, S. Hariharan, and A. K. Tiwari, J. Phys. Chem. C 121, 16351 (2017).
  • Gan et al. (2012) L.-Y. Gan, R.-Y. Tian, X.-B. Yang, H.-D. Lu, and Y.-J. Zhao, J. Phys. Chem. C 116, 745 (2012).
  • Turchanin and Agraval (2008) M. Turchanin and P. Agraval, Powder Metall. Met. Ceram. 47, 26 (2008).