Next Article in Journal
Magnetic Particle Imaging-Guided Thermal Simulations for Magnetic Particle Hyperthermia
Next Article in Special Issue
Optimizing P3HT/PCBM-Based Organic Photodetector Performance: Insights from SCAPS 1D Simulation Studies
Previous Article in Journal
Flow Behavior of Nanoparticle Agglomerates in a Fluidized Bed Simulated with Porous-Structure-Based Drag Laws
Previous Article in Special Issue
Light Absorption Analysis and Optimization of Ag@TiO2 Core-Shell Nanospheroid and Nanorod
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Surface Oxygen Coverage on the Oxygen Evolution Reaction over a CoFeNiCr High-Entropy Alloy

1
Department of Chemical, Environmental and Materials Engineering, University of Miami, Coral Gables, FL 33146, USA
2
Department of Civil and Architectural Engineering, University of Miami, Coral Gables, FL 33146, USA
*
Author to whom correspondence should be addressed.
Nanomaterials 2024, 14(12), 1058; https://doi.org/10.3390/nano14121058
Submission received: 22 May 2024 / Revised: 14 June 2024 / Accepted: 17 June 2024 / Published: 19 June 2024
(This article belongs to the Special Issue Modeling, Simulation and Optimization of Nanomaterials)

Abstract

:
Developing cost-effective and highly active electrocatalysts for the oxygen evolution reaction (OER) is crucial for advancing sustainable energy applications. High-entropy alloys (HEAs) made from earth-abundant transition metals, thanks to their remarkable stability and electrocatalytic performance, provide a promising alternative to expensive electrocatalysts typically derived from noble metals. While pristine HEA surfaces have been theoretically investigated, and the effect of oxygen coverage on conventional metal electrocatalysts has been examined, the impact of surface oxygen coverage on the electrocatalytic performance of HEAs remains poorly understood. To bridge this gap, we employ density functional theory (DFT) calculations to reconstruct the free energy diagram of OER intermediates on CoFeNiCr HEA surfaces with varying oxygen coverages, evaluating their impact on the rate-limiting step and theoretical overpotential. Our findings reveal that increased oxygen coverage weakens the adsorption of HO* and O*, but not HOO*. As a result, the theoretical overpotential for the OER decreases with higher oxygen coverage, and the rate-limiting step shifts from the third oxidation step (HOO* formation) at low coverage to the first oxidation step (HO* formation) at higher coverage.

Graphical Abstract

1. Introduction

The oxygen evolution reaction (OER), which oxidizes H2O to O2, is crucial for various sustainable energy technologies. For instance, electrochemical water splitting is a promising and environmentally friendly method for hydrogen production [1,2,3], whose efficiency is notably hindered by the high overpotential of the OER—its anodic half-reaction [4,5]. While noble metal-based oxides like IrO2 and RuO2 are currently the most practical OER electrocatalysts [6,7,8], with experimental overpotentials typically within the range of 200–300 mV at a current density of 10 mA cm−2 [9,10], their broader utilization in water splitting is limited by the high cost and scarcity of noble metals. In the pursuit of cost-effective OER electrocatalysts, high-entropy alloys (HEAs) have gained increasing attention due to their excellent OER activity and electrochemical stability [11,12,13]. HEAs are materials composed of four or more metals in roughly equal proportions, capable of forming single-phase solid solutions stabilized by the entropy of mixing [14,15,16]. Unlike conventional electrocatalyst materials, such as noble metals or their oxides, which have limited tunability in their catalytic activity, HEAs stand out due to their vast design space, the immense diversity of local binding environments that could act as active sites, and exceptional stability [14,17,18].
Considering the prevalence of earth-abundant transition metals such as Co, Fe, and Ni in OER electrocatalysts [19,20,21,22,23,24,25,26,27], HEAs developed for this purpose commonly incorporate two or three of these metals in their chemical composition, substituting noble metals [11,12]. For instance, Jin et al. investigated a nanoporous AlNiCoIrMo HEA [28], reporting both its excellent electrochemical stability and low overpotential of 233 mV at 10 mA cm−2. Studies on noble metal-free HEAs have also revealed promising results. Qiu et al. explored the effect of alloying various non-noble metals with CoFeNi on the OER performance, highlighting the AlNiCoFeMo HEA as a particularly noteworthy composition due to its significantly enhanced electrochemical durability and the lowest measured overpotential in that study of about 240 mV at 10 mA cm−2 [29]. Waag et al. demonstrated that the CoCrFeMnNi HEA exhibits comparable performance to commercial IrO2, with an average mass-specific activity of 0.104 A/mg [30]. Huang et al. found that lattice distortion in MnFeCoNiCu enhances its catalytic activity, leading to an overpotential of 263 mV at 10 mA cm−2 in the OER [31]. In addition to quinary HEAs, Dai et al. reported the promising OER performance of a quaternary MnFeCoNi HEA-based electrocatalyst, with an overpotential of 302 mV at the current density of 10 mA cm−2, comparable to that of RuO2 [32]. These experimental findings demonstrate the promise of earth-abundant transition metal HEAs for the OER. Besides HEAs, important theoretical and experimental work has also been conducted on evaluating the electrocatalytic properties of perovskite oxides, relatively new materials known for their versatile structure and tunable properties. This includes identifying universal activity descriptors [33], or exploiting lattice strain to tailor their adsorption characteristics and enhance their catalytic properties [34].
Previous theoretical studies using density functional theory (DFT) have shed light on the OER mechanisms and catalytic performance of HEAs by analyzing the free energy of OER intermediates on pristine HEA surfaces [35,36,37,38,39]. However, these studies have focused on pristine surfaces and overlooked the effect of surface oxygen coverage expected during the OER [6], which remains unknown. Considering the effect of oxygen coverage is of great importance as it can significantly influence the stability of OER intermediates and, consequently, affect the electrocatalytic performance of the HEA. The formation of an oxygen layer on the surface of the HEA during operation is likely due to the oxophilicity of pristine HEA surfaces made of transition metals like Co, Fe, Ni, and Cr [40].
Here, we use DFT calculations to study the effect of surface oxygen coverage on the adsorption energies on a CoFeNiCr HEA of OER intermediates (hydroxyl (HO*), oxygen (O*), and hydroperoxyl (HOO*) radicals), and we construct a free energy diagram based on the computational hydrogen electrode (CHE) model [41]. Through thermodynamic analysis, we gain insight into the impact of surface oxygen coverage on the theoretical OER overpotential and the associated rate-limiting step. The composition CoFeNiCr was selected because its constituent metals are among the most frequently utilized in HEAs [42], it has desirable mechanical properties [43], and it has shown promising results for hydrogen evolution reactions (HER) [44], making it a potential candidate for exploring dual-functionality in electrocatalytic applications.

2. Materials and Methods

DFT calculations were conducted using the Quickstep module of the CP2K 9.1 simulation package [45], which employs a dual basis of atom-centered Gaussian orbitals and Plane Waves known as the mixed Gaussian and plane wave (GPW) approach. The exchange-correlation functional revPBE-D3 was used, which integrates the generalized-gradient approximation (GGA) functional revPBE (revised Perdew–Burke–Ernzerhof) [46] with Grimme’s DFT-D3 dispersion corrections with zero-damping [47]. The core electrons of all the elements were represented by Goedecker–Teter–Hutter (GTH) pseudopotentials (PPs) [48,49], which have been shown to perform well across different systems and chemical conditions while also being computationally efficient for the GPW approach implemented in CP2K. The molecularly optimized (MOLOPT) Gaussian basis sets (BBs), which were specifically developed for use with GTH-PPs, were employed [50]. Specifically, the triple-zeta valence molecularly optimized basis set (TZV2P-MOLOPT-GTH) was used for H and O, which, as our previous work shows, offers a good balance between accuracy and computational cost [40]. This basis set is also the preferred choice for simulating ab initio liquid water [51,52]. The double-zeta valence short-range molecularly optimized basis set (DZVP-MOLOPT-SR-GTH) was used for Co, Fe, Ni, and Cr, as it is the largest available for those elements. The PPs and BSs used in this study were also previously employed to simulate the oxidation of copper [53].
To facilitate the convergence of the calculations, Fermi–Dirac smearing with an electronic temperature of 500 K was applied, along with the use of Broyden mixing [54]. The Brillouin zone was sampled at the Gamma point and periodic boundary conditions were applied in all three dimensions. The plane-wave energy cutoff for all calculations was set at 400 Ry. All DFT calculations involving oxygen intermediates were performed with spin polarization, with the only unpaired electrons being on the adsorbate. The HEAs were treated as non-magnetic. Previous research has reported a Curie temperature of approximately 130 K for the CoFeNiCr HEA [55,56]. Since the electrocatalytic process typically occurs at room temperature, we hypothesize that the magnetic properties of the HEA will have a modest impact under operating conditions. Furthermore, predicting the magnetic ground state of HEAs presents significant challenges due to their local compositional heterogeneity, which is beyond the scope of this paper.
The simulations featured slabs of a face-centered cubic (FCC) CoFeNiCr HEA with an FCC (111) surface and varying surface oxygen coverages. The HEA slabs were 4 × 4 atoms wide and 6 atoms thick, totaling 96 metal atoms. We assumed a random model for the HEA, where the atoms in the lattice were randomly assigned to Co, Fe, Ni, or Cr, maintaining the equimolar stoichiometry of the HEA. While this approximation is not entirely accurate, as theoretical studies have shown a tendency for Cr to segregate [57,58], it is a commonly used approximation to model HEAs. We expect Cr segregation to affect minor quantitative details, but not the qualitative trends reported here.
The following steps were carried out to generate each system. First, a bulk 4 × 4 × 6 HEA was optimized allowing adjustments to both the simulation box and the atomic positions using the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm [59] with a convergence criterion of 10−5 Hartree/Bohr. Next, a 10 Å vacuum layer was introduced along the z-direction of the simulation box to effectively simulate a slab with an FCC (111) surface. The slab’s geometry was then optimized, using the BFGS algorithm with a convergence criterion set at 10−3 Hartree/Bohr, permitting only the top three layers of atoms to relax. After optimizing the slab’s geometry, atomic oxygens were positioned on the FCC hollow sites—widely recognized as preferred binding sites for atomic oxygen on FCC (111) surfaces [60,61,62,63]—approximately 2 Å above the HEA surface. The positions of the surface oxygen atoms were optimized keeping the metal atoms fixed. Any oxygen atom that was not chemisorbed to the surface after the geometry optimization was removed from the system.
The surface oxygen coverage, θ O , is defined as the ratio of chemisorbed oxygen atoms to the total number of FCC hollow sites on the HEA surface. Thus, it is expressed as a fraction of a monolayer (ML). Considering the slightly fluctuating nature of the number of chemisorbed oxygen atoms on the HEA surfaces after optimization, we categorized θ O into three levels, ranging from low to high coverage, with θ O values falling within the ranges [ 0.25 ,   0.5 ) , [ 0.5 ,   0.75 ) , and [ 0.75 ,   1 ) , respectively. Figure 1a displays top-view representations of optimized systems corresponding to different levels of surface oxygen coverage.
Once the oxygen-covered HEA slabs were generated, we added the respective OER intermediates as adsorbates at random positions on the surface and about 3 Å above it. We optimized the geometry of the combined systems keeping the HEA and chemisorbed oxygen atoms fixed. Finally, the adsorption energy ( E ) of each intermediate was calculated using the following equations [9,64]:
E H O * = E H E A + H O * E H E A + 1 2 E H 2 E H 2 O
E O * = E H E A + O * E H E A + E H 2 E H 2 O
E H O O * = E H E A + H O O * E H E A + 3 2 E H 2 2 E H 2 O
where E H E A + H O * , E H E A + O * , and E H E A + H O O * represent the energies of the combined systems after geometry optimization, E H E A is the energy of the pristine or oxygen-covered slab, and E H 2 and E H 2 O are the energies of optimized H2 and H2O molecules in the gas phase. All the energies were corrected for basis set superposition errors (BSSEs). The adsorption free energy ( G ) of each intermediate was obtained by correcting the adsorption energy ( E ) for the zero-point energy and entropy using the following equation [41]:
G = E + Z P E T S
where G represents the change in reaction Gibbs free energy of each intermediate; Z P E and T S denote the difference in the zero-point energy and entropy, respectively. The zero-point energies and entropic corrections are taken from reference [41].
Table S1 provides a detailed breakdown of the number of adsorption energy calculations for each intermediate at each θ O , including the number of different HEA configurations, oxygen-covered configurations, and adsorbate configurations used in this study. In total, we carried out over 1200 DFT adsorption energy calculations.

3. Results

3.1. Non-Uniform Metal Oxidation

Our simulations revealed that oxygen atoms do not chemisorb equally to the different metals. As illustrated in Figure 1b, which shows the number of occurrences where the various constituent metals were found to be nearest neighbors to chemisorbed oxygen atoms, Cr exhibits the highest occurrence, followed by Fe and Co, while Ni shows the lowest occurrence. This trend reflects the high oxidation tendency of Cr and the low oxidation tendency of Ni, aligning with recent findings on the metal oxidation propensity in the CoFeNiCr HEA [44,65]. This non-uniform oxidation tendency is also reflected on the adsorption behavior of each constituent metal toward oxygen or oxygen-containing adsorbates. Specifically, Ni exhibits the weakest adsorption among all four elements [40,41], making surface Ni atoms the least likely to be the nearest neighbors of chemisorbed oxygens.

3.2. Surface Oxygen Coverage Impact on OER Intermediate Adsorption

Due to the wide diversity of binding sites on HEA surfaces, the adsorption energy of OER intermediates cannot be represented by a single value, as is typical for conventional electrocatalysts with one or a few active sites. Instead, it must be represented by a distribution of adsorption energies. Figure 2a–c show the distributions of adsorption energies for HO*, O*, and HOO*, respectively, at each level of surface oxygen coverage (top to bottom). The distributions display two distinct modes. Mode I, indicated by orange arrows in Figure 2a, corresponds to binding configurations where the adsorbate occupies a pristine site on the HEA surface, while Mode II, indicated by green arrows in Figure 2a, involves binding configurations where the adsorbate is on or near a site already occupied by a chemisorbed oxygen atom. Figure 2d provides snapshots from the simulations, illustrating both modes.
The most straightforward effect of surface oxygen coverage on the adsorption energies is its impact on the population of both modes. On pristine surfaces ( θ O = 0 ), Mode II obviously does not exist. As θ O increases, the number of available pristine binding sites on the surface decreases, leading to fewer Mode I and more Mode II occurrences. This unsurprising effect on the distributions is readily apparent in Figure 2a–c, progressing from the bottom to the top panels.
More interestingly, θ O has a non-trivial effect, different for each adsorbate, on the average adsorption energy associated with each mode (Figure 2e–g). For both HO* (Figure 2e) and O* (Figure 2f), as θ O increases, the average adsorption energy associated with Mode I shifts significantly towards weaker adsorption, while that associated with Mode II remains largely unaffected by surface oxygen coverage. The substantial shift in the average adsorption energy associated with Mode I for both HO* and O* suggests that adjacent chemisorbed oxygen atoms weaken their adsorption at pristine sites. This phenomenon has been previously reported for single metals [66,67]. In contrast, the average adsorption energy associated with Mode II for all intermediates remains mostly constant irrespective of θ O , indicating that adsorption at oxidized sites is minimally affected by other nearby chemisorbed oxygens. For HOO*, the average adsorption energy associated with both modes remains consistent across different θ O ranges (Figure 2g). However, unlike for HO* and O*, increasing θ O significantly impacts the frequency of spontaneous HOO* dissociation. Our simulations reveal that on pristine HEA surfaces, the O-O bond in HOO* dissociates in 95% of the cases during geometry optimization. This frequency drops to approximately 50% when θ O increases to the range [ 0.25 ,   0.5 ) and further declines to around 14% when θ O is higher than 0.5 ML. This phenomenon underscores the strong tendency of HEA surfaces to oxidize [68].
In an effort to gain mechanistic understanding of the adsorption of HO*, we examined how the tilt angle, α , of HO* with respect to the surface of the HEA varies depending on the oxygen coverage and whether HO* is adsorbed to a pristine or oxidized site (Figure 3a). According to our definition, when α = 0 ° , HO* is vertical on the surface with the hydrogen pointing away from it. For pristine sites (i.e., Mode I), HO* typically adopts a mostly vertical configuration at hollow sites ( α 0 °   t o   30 ° ) and a more tilted configuration at bridge sites ( α 30 °   t o   60 ° ), as evidenced by the bimodal distributions for pristine sites shown in blue in Figure 3a. The average tilt angle at pristine hollow sites increases with oxygen coverage, from approximately 5° on pristine surfaces to about 15° for θ O [ 0.25 ,   0.5 ) . For Mode II adsorption, HO* configurations range from parallel to the surface ( α = 90 ° ) to having the hydrogen in HO* pointing towards the HEA surface ( α > 90 ° ). Interestingly, Figure 3b reveals a correlation between adsorption energy and tilt angle, with smaller tilt angles resulting in stronger adsorption energies.
A recent study by Cao and Nørskov [69] demonstrated that spin polarization results in consistently weaker adsorption energies for a wide range of adsorbates on single metal surfaces such as Co, Fe, and Ni. We hypothesize that similar trends would be observed for the CoFeNiCr HEA in this study if magnetic effects were considered, particularly for the adsorption to pristine sites. Given that this weakening effect was uniform across all adsorbates studied, it is reasonable to believe that the qualitative trends we observe here would remain unchanged.

3.3. Free Energy Diagram for the OER

To assess how changes in the adsorption energy of each intermediate as θ O varies might influence the electrocatalytic performance of HEAs, we constructed and analyzed the free energy diagrams for the OER, assuming it proceeds through the four-electron pathway mechanism:
2 H 2 O H O * + H 2 O + H + + e
H O * + H 2 O + H + + e O * + H 2 O + 2 H + + 2 e
O * + H 2 O + 2 H + + 2 e H O O * + 3 H + + 3 e
H O O * + 3 H + + 3 e O 2 + 4 H + + 4 e
The free energy difference ( G ) of each oxidation step i (i = 1, 2, 3, or 4) at pH = 0 and zero potential (U = 0 V) is calculated as follows:
G 1 = G H O *
G 2 = G O * G H O *
G 3 = G H O O * G O *
G 4 = G O 2 G H O O *
Figure 4 contains the free energy diagrams of the OER on HEAs for the three surface oxygen coverage ranges studied here. The analysis excludes pristine HEA surfaces as HOO* was found to dissociate in most cases on those surfaces. In Figure 4, we employ box plots to visually capture the spread of adsorption free energies for each intermediate, a representation unnecessary for conventional catalysts.
By analyzing the median free energy of each reaction step as a function of θ O , we observe a significant increase in G 1 , the free energy difference for the initial oxidation step from H2O to HO*, as θ O increases. This increase in G 1 is due to the overall weakening of HO* adsorption, which results from both the diminished number of available pristine binding sites on the surface and the reduced adsorption strength at the remaining pristine sites. The free energy difference G 2 , the step from HO* to O*, also rises as θ O increases, albeit not as much as the increase observed for G 1 . The simultaneous weakening of HO* and O* adsorption on surfaces with higher oxygen coverages accounts for the observed trend. Since the adsorption of HOO* remains more or less the same across different oxygen coverages, the overall reduction in O* adsorption strength results in a significant decrease in the free energy difference G 3 , the third oxidation step from O* to HOO*, on surfaces with higher oxygen coverages.

3.4. Theoretical Overpotential for OER

The theoretical overpotential ( η O E R ) is calculated as the difference between the highest free energy change among all four oxidation steps and the theoretically minimum potential difference required for the reaction to proceed (1.23 V):
η O E R = max G 1 , G 2 , G 3 , G 4 e 1.23 V
In contrast to the previous section, which focused on the median free energy values of the distributions, this section calculates G 1 , G 2 , G 3 , G 4 , and therefore η O E R , at the level of individual binding sites on each HEA. Specifically, we only include the subset of cases where all three intermediates, after optimization, are adsorbed to the same binding site on the same oxygen-covered HEA surface. We consider two intermediates to be adsorbed to the same site if the relative distance between the centers of the oxygen atoms is less than 1.5 (which is approximately the distance between an adsorbed oxygen atom at a pristine site and the nearest neighbor metallic elements). Figure 5a shows snapshots of an example where all the intermediates were found to be adsorbed to the same binding site on the HEA surface. The number of such sites that we were able to sample is quantified in Table 1.
The distribution of η O E R for each θ O range is shown in Figure 5b. We find that G 4 is never the rate-limiting step for any of the binding sites studied here, and G 2 serves as the rate-limiting step in less than around 5% of the cases. The data in Table 1 suggest a strong correlation between the type of binding site (i.e., pristine site or Mode I versus oxidized site or Mode II) and the rate-limiting step of the reaction. Pristine sites predominantly result in G 3 being the highest free energy difference, which corresponds to the step from O* to HOO*. Conversely, oxidized sites typically lead to G 1 being the highest free energy difference, corresponding to the initial oxidation step from H2O to HO*. As the number of oxidized sites increases at the expense of pristine sites with increasing θ O , the rate-limiting step shifts from G 3 to G 1 .
In Figure 5b, it is clear that the η O E R associated with G 3 (yellow), typically linked to pristine sites, displays a broader distribution than that associated with G 1 (blue). This is because the composition of the binding site has a greater impact on pristine sites than on oxidized sites. The η O E R associated with G 3 (i.e., pristine binding sites) can range from as low as 0.96 V observed at θ O [ 0.5 ,   0.75 ) to as high as 4.97 V observed at θ O [ 0.25 ,   0.5 ) . We observe a shift in the rate-limiting step as θ O increases (Figure 5c). The proportion of rate-limiting steps attributed to G 3 decreases from about 63% at the lowest oxygen coverage to about 10% for θ O [ 0.75 ,   1 ) , while that attributed to G 1 increases from about 32% to over 90% across the same ranges of oxygen coverage.
Our mechanistic understanding for the shift in the rate-limiting step is as follows: Our simulations show that the adsorption of HOO* remains largely unaffected by θ O . In contrast, the adsorption of HO* and O*, which is initially strong at low θ O , weakens as θ O increases. As mentioned earlier, this weakening is due to the decreasing number of available pristine sites and the reduced adsorption energy of the remaining pristine sites. Consequently, at low θ O , the rate-limiting step is the third oxidation step (from O* to HOO*), primarily due to the weak adsorption of HOO*. As both HO* and O* adsorption weaken with increasing θ O (while HOO* adsorption remains unchanged), there comes a point where the adsorption of HO* becomes so weak that the first oxidation step (from H2O to HO*) becomes the rate-limiting step.
When examining the average overall overpotential, η O E R ¯ , at each oxygen coverage range (Figure 5d), we observe a general decrease as θ O increases. This suggests that the catalytic performance towards the OER improves with higher oxygen coverage. This improvement is primarily driven by the decreasing average overpotential associated with G 3 , as indicated in the inset of Figure 5d. In contrast, the average overpotential associated with G 1 remains consistent, regardless of oxygen coverage.

4. Conclusions

In this study, we used DFT simulations to investigate the effect of surface oxygen coverage ( θ O ) on a CoFeNiCr HEA in the oxygen evolution reaction (OER). We found that as surface oxygen coverage increases, the adsorption of intermediates HO* and O* weakens notably. This weakening is attributed to the rise in the number of oxidized binding sites, which exhibit weaker adsorption than pristine sites, and the reduced adsorption on the remaining pristine sites due to the influence of neighboring chemisorbed oxygens. Interestingly, the adsorption energy of HOO* is mostly independent of the level of oxygen coverage. This non-uniform effect on the adsorption energies of the different intermediates causes a transition in the rate-limiting step from the third oxidation step (formation of HOO*) to the first oxidation step (formation of HO*) as surface oxygen coverage rises. Furthermore, the average theoretical overpotential decreases as θ O increases.
Our work marks a significant step towards understanding the electrocatalytic potential of HEAs. Beyond their use in the OER, our simulations reveal a wide range of adsorption energies and overpotentials associated with the diversity of binding sites, highlighting their potential as multifunctional electrocatalysts.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano14121058/s1, Table S1: Number of adsorption energies obtained from DFT calculations for each intermediate at different coverage ranges.

Author Contributions

Conceptualization, L.R.P.; methodology, L.R.P. and G.Y.; software, G.Y. and L.R.P.; formal analysis, G.Y.; investigation, L.R.P. and G.Y.; writing—original draft preparation, G.Y.; writing—review and editing, L.R.P. and G.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the U.S. National Science Foundation under award #2230165 (Environmental Engineering program in the Chemical, Bioengineering, Environmental and Transport Systems division).

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Acknowledgments

We gratefully acknowledge the support received from the high-performance computing resources provided by the supercomputer Triton, managed by the Institute of Data Science and Computing (IDSC) at the University of Miami.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Chu, S.; Majumdar, A. Opportunities and Challenges for a Sustainable Energy Future. Nature 2012, 488, 294–303. [Google Scholar] [CrossRef]
  2. Bockris, J.O.; Veziroglu, T.N. Estimates of the Price of Hydrogen as a Medium for Wind and Solar Sources. Int. J. Hydrogen Energy 2007, 32, 1605–1610. [Google Scholar] [CrossRef]
  3. Hassan, N.S.; Jalil, A.A.; Rajendran, S.; Khusnun, N.F.; Bahari, M.B.; Johari, A.; Kamaruddin, M.J.; Ismail, M. Recent Review and Evaluation of Green Hydrogen Production via Water Electrolysis for a Sustainable and Clean Energy Society. Int. J. Hydrogen Energy 2024, 52, 420–441. [Google Scholar] [CrossRef]
  4. Trotochaud, L.; Boettcher, S.W. Precise Oxygen Evolution Catalysts: Status and Opportunities. Scr. Mater. 2014, 74, 25–32. [Google Scholar] [CrossRef]
  5. Suen, N.-T.; Hung, S.-F.; Quan, Q.; Zhang, N.; Xu, Y.-J.; Ming Chen, H. Electrocatalysis for the Oxygen Evolution Reaction: Recent Development and Future Perspectives. Chem. Soc. Rev. 2017, 46, 337–365. [Google Scholar] [CrossRef]
  6. Rossmeisl, J.; Qu, Z.-W.; Zhu, H.; Kroes, G.-J.; Nørskov, J.K. Electrolysis of Water on Oxide Surfaces. J. Electroanal. Chem. 2007, 607, 83–89. [Google Scholar] [CrossRef]
  7. Lee, Y.; Suntivich, J.; May, K.J.; Perry, E.E.; Shao-Horn, Y. Synthesis and Activities of Rutile IrO2 and RuO2 Nanoparticles for Oxygen Evolution in Acid and Alkaline Solutions. J. Phys. Chem. Lett. 2012, 3, 399–404. [Google Scholar] [CrossRef]
  8. Shi, Q.; Zhu, C.; Du, D.; Lin, Y. Robust Noble Metal-Based Electrocatalysts for Oxygen Evolution Reaction. Chem. Soc. Rev. 2019, 48, 3181–3192. [Google Scholar] [CrossRef]
  9. Man, I.C.; Su, H.-Y.; Calle-Vallejo, F.; Hansen, H.A.; Martínez, J.I.; Inoglu, N.G.; Kitchin, J.; Jaramillo, T.F.; Nørskov, J.K.; Rossmeisl, J. Universality in Oxygen Evolution Electrocatalysis on Oxide Surfaces. ChemCatChem 2011, 3, 1159–1165. [Google Scholar] [CrossRef]
  10. Matsumoto, Y.; Sato, E. Electrocatalytic Properties of Transition Metal Oxides for Oxygen Evolution Reaction. Mater. Chem. Phys. 1986, 14, 397–426. [Google Scholar] [CrossRef]
  11. Xin, Y.; Li, S.; Qian, Y.; Zhu, W.; Yuan, H.; Jiang, P.; Guo, R.; Wang, L. High-Entropy Alloys as a Platform for Catalysis: Progress, Challenges, and Opportunities. ACS Catal. 2020, 10, 11280–11306. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Wang, D.; Wang, S. High-Entropy Alloys for Electrocatalysis: Design, Characterization, and Applications. Small 2022, 18, 2104339. [Google Scholar] [CrossRef]
  13. Ren, J.-T.; Chen, L.; Wang, H.-Y.; Yuan, Z.-Y. High-Entropy Alloys in Electrocatalysis: From Fundamentals to Applications. Chem. Soc. Rev. 2023, 52, 8319–8373. [Google Scholar] [CrossRef]
  14. Yeh, J.-W.; Chen, S.-K.; Lin, S.-J.; Gan, J.-Y.; Chin, T.-S.; Shun, T.-T.; Tsau, C.-H.; Chang, S.-Y. Nanostructured High-Entropy Alloys with Multiple Principal Elements: Novel Alloy Design Concepts and Outcomes. Adv. Eng. Mater. 2004, 6, 299–303. [Google Scholar] [CrossRef]
  15. Troparevsky, M.C.; Morris, J.R.; Kent, P.R.; Lupini, A.R.; Stocks, G.M. Criteria for Predicting the Formation of Single-Phase High-Entropy Alloys. Phys. Rev. X 2015, 5, 011041. [Google Scholar] [CrossRef]
  16. King, D.; Middleburgh, S.; McGregor, A.; Cortie, M. Predicting the Formation and Stability of Single Phase High-Entropy Alloys. Acta Mater. 2016, 104, 172–179. [Google Scholar] [CrossRef]
  17. George, E.P.; Raabe, D.; Ritchie, R.O. High-Entropy Alloys. Nat. Rev. Mater. 2019, 4, 515–534. [Google Scholar] [CrossRef]
  18. Löffler, T.; Ludwig, A.; Rossmeisl, J.; Schuhmann, W. What Makes High-Entropy Alloys Exceptional Electrocatalysts? Angew. Chem. Int. Ed. 2021, 60, 26894–26903. [Google Scholar] [CrossRef]
  19. McCrory, C.C.L.; Jung, S.; Peters, J.C.; Jaramillo, T.F. Benchmarking Heterogeneous Electrocatalysts for the Oxygen Evolution Reaction. J. Am. Chem. Soc. 2013, 135, 16977–16987. [Google Scholar] [CrossRef]
  20. Liu, C.; Colón, B.C.; Ziesack, M.; Silver, P.A.; Nocera, D.G. Water Splitting–Biosynthetic System with CO2 Reduction Efficiencies Exceeding Photosynthesis. Science 2016, 352, 1210–1213. [Google Scholar] [CrossRef]
  21. Nardi, K.L.; Yang, N.; Dickens, C.F.; Strickler, A.L.; Bent, S.F. Creating Highly Active Atomic Layer Deposited NiO Electrocatalysts for the Oxygen Evolution Reaction. Adv. Energy Mater. 2015, 5, 1500412. [Google Scholar] [CrossRef]
  22. Gong, M.; Li, Y.; Wang, H.; Liang, Y.; Wu, J.Z.; Zhou, J.; Wang, J.; Regier, T.; Wei, F.; Dai, H. An Advanced Ni–Fe Layered Double Hydroxide Electrocatalyst for Water Oxidation. J. Am. Chem. Soc. 2013, 135, 8452–8455. [Google Scholar] [CrossRef]
  23. Chen, S.; Qiao, S.-Z. Hierarchically Porous Nitrogen-Doped Graphene–NiCo2O4 Hybrid Paper as an Advanced Electrocatalytic Water-Splitting Material. ACS Nano 2013, 7, 10190–10196. [Google Scholar] [CrossRef]
  24. Gao, M.; Sheng, W.; Zhuang, Z.; Fang, Q.; Gu, S.; Jiang, J.; Yan, Y. Efficient Water Oxidation Using Nanostructured α-Nickel-Hydroxide as an Electrocatalyst. J. Am. Chem. Soc. 2014, 136, 7077–7084. [Google Scholar] [CrossRef]
  25. Chen, J.; Zheng, F.; Zhang, S.-J.; Fisher, A.; Zhou, Y.; Wang, Z.; Li, Y.; Xu, B.-B.; Li, J.-T.; Sun, S.-G. Interfacial Interaction between FeOOH and Ni–Fe LDH to Modulate the Local Electronic Structure for Enhanced OER Electrocatalysis. ACS Catal. 2018, 8, 11342–11351. [Google Scholar] [CrossRef]
  26. Wang, Y.; Zhang, Y.; Liu, Z.; Xie, C.; Feng, S.; Liu, D.; Shao, M.; Wang, S. Layered Double Hydroxide Nanosheets with Multiple Vacancies Obtained by Dry Exfoliation as Highly Efficient Oxygen Evolution Electrocatalysts. Angew. Chem. Int. Ed. 2017, 56, 5867–5871. [Google Scholar] [CrossRef]
  27. Liu, R.; Wang, Y.; Liu, D.; Zou, Y.; Wang, S. Water-Plasma-Enabled Exfoliation of Ultrathin Layered Double Hydroxide Nanosheets with Multivacancies for Water Oxidation. Adv. Mater. 2017, 29, 1701546. [Google Scholar] [CrossRef]
  28. Jin, Z.; Lv, J.; Jia, H.; Liu, W.; Li, H.; Chen, Z.; Lin, X.; Xie, G.; Liu, X.; Sun, S.; et al. Nanoporous Al-Ni-Co-Ir-Mo High-Entropy Alloy for Record-High Water Splitting Activity in Acidic Environments. Small 2019, 15, 1904180. [Google Scholar] [CrossRef]
  29. Qiu, H.-J.; Fang, G.; Gao, J.; Wen, Y.; Lv, J.; Li, H.; Xie, G.; Liu, X.; Sun, S. Noble Metal-Free Nanoporous High-Entropy Alloys as Highly Efficient Electrocatalysts for Oxygen Evolution Reaction. ACS Mater. Lett. 2019, 1, 526–533. [Google Scholar] [CrossRef]
  30. Waag, F.; Li, Y.; Rosa Ziefuß, A.; Bertin, E.; Kamp, M.; Duppel, V.; Marzun, G.; Kienle, L.; Barcikowski, S.; Gökce, B. Kinetically-Controlled Laser-Synthesis of Colloidal High-Entropy Alloy Nanoparticles. RSC Adv. 2019, 9, 18547–18558. [Google Scholar] [CrossRef]
  31. Huang, K.; Zhang, B.; Wu, J.; Zhang, T.; Peng, D.; Cao, X.; Zhang, Z.; Li, Z.; Huang, Y. Exploring the Impact of Atomic Lattice Deformation on Oxygen Evolution Reactions Based on a Sub-5 Nm Pure Face-Centred Cubic High-Entropy Alloy Electrocatalyst. J. Mater. Chem. A 2020, 8, 11938–11947. [Google Scholar] [CrossRef]
  32. Dai, W.; Lu, T.; Pan, Y. Novel and Promising Electrocatalyst for Oxygen Evolution Reaction Based on MnFeCoNi High Entropy Alloy. J. Power Sources 2019, 430, 104–111. [Google Scholar] [CrossRef]
  33. Guan, D.; Xu, H.; Zhang, Q.; Huang, Y.-C.; Shi, C.; Chang, Y.-C.; Xu, X.; Tang, J.; Gu, Y.; Pao, C.-W.; et al. Identifying a Universal Activity Descriptor and a Unifying Mechanism Concept on Perovskite Oxides for Green Hydrogen Production. Adv. Mater. 2023, 35, 2305074. [Google Scholar] [CrossRef]
  34. Guan, D.; Zhong, J.; Xu, H.; Huang, Y.-C.; Hu, Z.; Chen, B.; Zhang, Y.; Ni, M.; Xu, X.; Zhou, W.; et al. A Universal Chemical-Induced Tensile Strain Tuning Strategy to Boost Oxygen-Evolving Electrocatalysis on Perovskite Oxides. Appl. Phys. Rev. 2022, 9, 011422. [Google Scholar] [CrossRef]
  35. He, R.; Yang, L.; Zhang, Y.; Jiang, D.; Lee, S.; Horta, S.; Liang, Z.; Lu, X.; Ostovari Moghaddam, A.; Li, J.; et al. A 3d-4d-5d High Entropy Alloy as a Bifunctional Oxygen Catalyst for Robust Aqueous Zinc–Air Batteries. Adv. Mater. 2023, 35, 2303719. [Google Scholar] [CrossRef]
  36. Lu, Y.; Huang, K.; Cao, X.; Zhang, L.; Wang, T.; Peng, D.; Zhang, B.; Liu, Z.; Wu, J.; Zhang, Y.; et al. Atomically Dispersed Intrinsic Hollow Sites of M-M1-M (M1 = Pt, Ir; M = Fe, Co, Ni, Cu, Pt, Ir) on FeCoNiCuPtIr Nanocrystals Enabling Rapid Water Redox. Adv. Funct. Mater. 2022, 32, 2110645. [Google Scholar] [CrossRef]
  37. Zhu, H.; Zhu, Z.; Hao, J.; Sun, S.; Lu, S.; Wang, C.; Ma, P.; Dong, W.; Du, M. High-Entropy Alloy Stabilized Active Ir for Highly Efficient Acidic Oxygen Evolution. Chem. Eng. J. 2022, 431, 133251. [Google Scholar] [CrossRef]
  38. Sharma, L.; Katiyar, N.K.; Parui, A.; Das, R.; Kumar, R.; Tiwary, C.S.; Singh, A.K.; Halder, A.; Biswas, K. Low-Cost High Entropy Alloy (HEA) for High-Efficiency Oxygen Evolution Reaction (OER). Nano Res. 2022, 15, 4799–4806. [Google Scholar] [CrossRef]
  39. Zhang, Z.-J.; Guo, J.-P.; Sun, S.-H.; Sun, Q.; Zhao, Y.-W.; Zhang, Y.-F.; Yu, Z.-Y.; Li, C.-S.; Sun, Y.; Zhang, M.-M.; et al. Optimized Valence State of Co and Ni in High-Entropy Alloy for High Active-Stable OER. Rare Met. 2023, 42, 3607–3613. [Google Scholar] [CrossRef]
  40. Yuan, G.; Wu, M.; Ruiz Pestana, L. Density Functional Theory-Machine Learning Characterization of the Adsorption Energy of Oxygen Intermediates on High-Entropy Alloys Made of Earth-Abundant Metals. J. Phys. Chem. C 2023, 127, 15809–15818. [Google Scholar] [CrossRef]
  41. Nørskov, J.K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.; Kitchin, J.R.; Bligaard, T.; Jonsson, H. Origin of the Overpotential for Oxygen Reduction at a Fuel-Cell Cathode. J. Phys. Chem. B 2004, 108, 17886–17892. [Google Scholar] [CrossRef]
  42. Miracle, D.B.; Senkov, O.N. A Critical Review of High Entropy Alloys and Related Concepts. Acta Mater. 2017, 122, 448–511. [Google Scholar] [CrossRef]
  43. Liu, B.; Wang, J.; Liu, Y.; Fang, Q.; Wu, Y.; Chen, S.; Liu, C.T. Microstructure and Mechanical Properties of Equimolar FeCoCrNi High Entropy Alloy Prepared via Powder Extrusion. Intermetallics 2016, 75, 25–30. [Google Scholar] [CrossRef]
  44. McKay, F.; Fang, Y.; Kizilkaya, O.; Singh, P.; Johnson, D.D.; Roy, A.; Young, D.P.; Sprunger, P.T.; Flake, J.C.; Shelton, W.A. CoCrFeNi High-Entropy Alloy as an Enhanced Hydrogen Evolution Catalyst in an Acidic Solution. J. Phys. Chem. C 2021, 125, 17008–17018. [Google Scholar] [CrossRef]
  45. Kühne, T.D.; Iannuzzi, M.; Del Ben, M.; Rybkin, V.V.; Seewald, P.; Stein, F.; Laino, T.; Khaliullin, R.Z.; Schütt, O.; Schiffmann, F. CP2K: An Electronic Structure and Molecular Dynamics Software Package-Quickstep: Efficient and Accurate Electronic Structure Calculations. J. Chem. Phys. 2020, 152, 194103. [Google Scholar] [CrossRef]
  46. Zhang, Y.; Yang, W. Comment on “Generalized Gradient Approximation Made Simple”. Phys. Rev. Lett. 1998, 80, 890. [Google Scholar] [CrossRef]
  47. 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, 154104. [Google Scholar] [CrossRef]
  48. Goedecker, S.; Teter, M.; Hutter, J. Separable Dual-Space Gaussian Pseudopotentials. Phys. Rev. B 1996, 54, 1703. [Google Scholar] [CrossRef]
  49. Hartwigsen, C.; Gœdecker, S.; Hutter, J. Relativistic Separable Dual-Space Gaussian Pseudopotentials from H to Rn. Phys. Rev. B 1998, 58, 3641. [Google Scholar] [CrossRef]
  50. VandeVondele, J.; Hutter, J. Gaussian Basis Sets for Accurate Calculations on Molecular Systems in Gas and Condensed Phases. J. Chem. Phys. 2007, 127, 114105. [Google Scholar] [CrossRef]
  51. Ruiz Pestana, L.; Mardirossian, N.; Head-Gordon, M.; Head-Gordon, T. Ab Initio Molecular Dynamics Simulations of Liquid Water Using High Quality Meta-GGA Functionals. Chem. Sci. 2017, 8, 3554–3565. [Google Scholar] [CrossRef]
  52. Ruiz Pestana, L.; Marsalek, O.; Markland, T.E.; Head-Gordon, T. The Quest for Accurate Liquid Water Properties from First Principles. J. Phys. Chem. Lett. 2018, 9, 5009–5016. [Google Scholar] [CrossRef]
  53. Nilsson, S.; El Berch, J.N.; Albinsson, D.; Fritzsche, J.; Mpourmpakis, G.; Langhammer, C. The Role of Grain Boundary Sites for the Oxidation of Copper Catalysts during the CO Oxidation Reaction. ACS Nano 2023, 17, 20284–20298. [Google Scholar] [CrossRef]
  54. Broyden, C.G. A Class of Methods for Solving Nonlinear Simultaneous Equations. Math. Comput. 1965, 19, 577–593. [Google Scholar] [CrossRef]
  55. Lucas, M.S.; Belyea, D.; Bauer, C.; Bryant, N.; Michel, E.; Turgut, Z.; Leontsev, S.O.; Horwath, J.; Semiatin, S.L.; McHenry, M.E.; et al. Thermomagnetic Analysis of FeCoCrxNi Alloys: Magnetic Entropy of High-Entropy Alloys. J. Appl. Phys. 2013, 113, 17A923. [Google Scholar] [CrossRef]
  56. Lucas, M.S.; Mauger, L.; Muñoz, J.A.; Xiao, Y.; Sheets, A.O.; Semiatin, S.L.; Horwath, J.; Turgut, Z. Magnetic and Vibrational Properties of High-Entropy Alloys. J. Appl. Phys. 2011, 109, 07E307. [Google Scholar] [CrossRef]
  57. Wynblatt, P.; Chatain, D. Modeling Grain Boundary and Surface Segregation in Multicomponent High-Entropy Alloys. Phys. Rev. Mater. 2019, 3, 054004. [Google Scholar] [CrossRef]
  58. Ferrari, A.; Körmann, F. Surface Segregation in Cr-Mn-Fe-Co-Ni High Entropy Alloys. Appl. Surf. Sci. 2020, 533, 147471. [Google Scholar] [CrossRef]
  59. Pfrommer, B.G.; Côté, M.; Louie, S.G.; Cohen, M.L. Relaxation of Crystals with the Quasi-Newton Method. J. Comput. Phys. 1997, 131, 233–240. [Google Scholar] [CrossRef]
  60. Mavrikakis, M.; Rempel, J.; Greeley, J.; Hansen, L.B.; Nørskov, J.K. Atomic and Molecular Adsorption on Rh(111). J. Chem. Phys. 2002, 117, 6737–6744. [Google Scholar] [CrossRef]
  61. Ford, D.C.; Xu, Y.; Mavrikakis, M. Atomic and Molecular Adsorption on Pt(111). Surf. Sci. 2005, 587, 159–174. [Google Scholar] [CrossRef]
  62. Bai, Y.; Kirvassilis, D.; Xu, L.; Mavrikakis, M. Atomic and Molecular Adsorption on Ni(111). Surf. Sci. 2019, 679, 240–253. [Google Scholar] [CrossRef]
  63. Eichler, A.; Mittendorfer, F.; Hafner, J. Precursor-Mediated Adsorption of Oxygen on the (111) Surfaces of Platinum-Group Metals. Phys. Rev. B 2000, 62, 4744–4755. [Google Scholar] [CrossRef]
  64. Batchelor, T.A.; Pedersen, J.K.; Winther, S.H.; Castelli, I.E.; Jacobsen, K.W.; Rossmeisl, J. High-Entropy Alloys as a Discovery Platform for Electrocatalysis. Joule 2019, 3, 834–845. [Google Scholar] [CrossRef]
  65. McKay, F.; Ismael, T.; Robinson, A.; Kizilkaya, O.; Sprunger, P.T.; Derosa, P.A. Surface Oxidation Mechanism of CoCrFeNi High Entropy Alloy. Surf. Sci. 2022, 723, 122124. [Google Scholar] [CrossRef]
  66. Yang, H.; Whitten, J.L. Energetics of Hydroxyl and Influence of Coadsorbed Oxygen on Metal Surfaces. J. Phys. Chem. B 1997, 101, 4090–4096. [Google Scholar] [CrossRef]
  67. Verga, L.G.; Aarons, J.; Sarwar, M.; Thompsett, D.; Russell, A.E.; Skylaris, C.-K. DFT Calculation of Oxygen Adsorption on Platinum Nanoparticles: Coverage and Size Effects. Faraday Discuss. 2018, 208, 497–522. [Google Scholar] [CrossRef]
  68. Montemore, M.M.; van Spronsen, M.A.; Madix, R.J.; Friend, C.M. O2 Activation by Metal Surfaces: Implications for Bonding and Reactivity on Heterogeneous Catalysts. Chem. Rev. 2018, 118, 2816–2862. [Google Scholar] [CrossRef]
  69. Cao, A.; Nørskov, J.K. Spin Effects in Chemisorption and Catalysis. ACS Catal. 2023, 13, 3456–3462. [Google Scholar] [CrossRef]
Figure 1. (a) Top view of optimized oxygen-covered HEA configurations at different θ O ranges, including the pristine surface configuration ( θ O = 0 ) . The red atoms represent chemisorbed oxygen atoms, while pink, blue, orange, and green correspond to Co, Cr, Fe, and Ni atoms, respectively. (b) The number of occurrences where the various constituent metals were found to be nearest neighbors to chemisorbed oxygen atoms.
Figure 1. (a) Top view of optimized oxygen-covered HEA configurations at different θ O ranges, including the pristine surface configuration ( θ O = 0 ) . The red atoms represent chemisorbed oxygen atoms, while pink, blue, orange, and green correspond to Co, Cr, Fe, and Ni atoms, respectively. (b) The number of occurrences where the various constituent metals were found to be nearest neighbors to chemisorbed oxygen atoms.
Nanomaterials 14 01058 g001
Figure 2. (ac) Adsorption energy distributions for HO*, O*, and HOO*, respectively, on HEA surfaces with different θ O , including pristine surfaces (bottom, θ O = 0 ). The orange and green arrows point to the peaks corresponding to Mode I and Mode II, respectively. The data visible in some of the distributions in panel (b), around 6–7 eV, correspond to desorbed O*. (d) Snapshots illustrating side and top views of example binding configurations of HO* corresponding to Mode I and Mode II of the distributions. The oxygen atom of HO* is colored in yellow. (eg) Average adsorption energy of HO*, O*, and HOO*, respectively, for both modes as a function of the oxygen coverage. The error bars represent standard deviation.
Figure 2. (ac) Adsorption energy distributions for HO*, O*, and HOO*, respectively, on HEA surfaces with different θ O , including pristine surfaces (bottom, θ O = 0 ). The orange and green arrows point to the peaks corresponding to Mode I and Mode II, respectively. The data visible in some of the distributions in panel (b), around 6–7 eV, correspond to desorbed O*. (d) Snapshots illustrating side and top views of example binding configurations of HO* corresponding to Mode I and Mode II of the distributions. The oxygen atom of HO* is colored in yellow. (eg) Average adsorption energy of HO*, O*, and HOO*, respectively, for both modes as a function of the oxygen coverage. The error bars represent standard deviation.
Nanomaterials 14 01058 g002
Figure 3. (a) Distribution of the tilt angle, α , of HO* with respect to the HEA surface for different oxygen coverages and for HO* adsorbed to pristine (blue) or oxidized (red) sites. The inset in the top plot illustrates the definition of α . Insets in the other plots show examples of the configuration of HO* corresponding to some notable peaks in the distributions, indicated by arrows. (b) Correlation between Δ E H O * and α for different oxygen coverages, with HO* adsorbed to pristine (blue) or oxidized (red) sites.
Figure 3. (a) Distribution of the tilt angle, α , of HO* with respect to the HEA surface for different oxygen coverages and for HO* adsorbed to pristine (blue) or oxidized (red) sites. The inset in the top plot illustrates the definition of α . Insets in the other plots show examples of the configuration of HO* corresponding to some notable peaks in the distributions, indicated by arrows. (b) Correlation between Δ E H O * and α for different oxygen coverages, with HO* adsorbed to pristine (blue) or oxidized (red) sites.
Nanomaterials 14 01058 g003
Figure 4. Free energy diagrams for oxygen evolution reaction (OER) on oxygen-covered HEA surfaces at zero potential. (a) θ O [ 0.25 ,   0.5 ) , (b) θ O [ 0.5 ,   0.75 ) , and (c) θ O [ 0.75 ,   1 ) . Each box plot represents the distribution of adsorption free energies across all binding site compositions for each intermediate. The solid line within each box plot represents the median adsorption free energy.
Figure 4. Free energy diagrams for oxygen evolution reaction (OER) on oxygen-covered HEA surfaces at zero potential. (a) θ O [ 0.25 ,   0.5 ) , (b) θ O [ 0.5 ,   0.75 ) , and (c) θ O [ 0.75 ,   1 ) . Each box plot represents the distribution of adsorption free energies across all binding site compositions for each intermediate. The solid line within each box plot represents the median adsorption free energy.
Nanomaterials 14 01058 g004
Figure 5. (a) Snapshots of an example where all the intermediates were found to be adsorbed to the same binding site on the HEA surface. The central oxygen atoms of HO*, O*, and HOO* are highlighted in yellow. (b) Distributions of theoretical overpotentials ( η O E R ) for each θ O range. The η O E R values corresponding to different rate-limiting steps in the OER are depicted in blue, orange, and yellow, representing G 1 , G 2 , and G 3 , respectively. (c) Percentage of times that G 1 , G 2 , or G 3 were the rate-limiting steps as a function of θ O . (d) Average theoretical overpotential ( η O E R ¯ ) accounting for contributions from all rate-limiting steps. The error bars indicate the standard deviation of the total distributions shown in Panel (b). The inset shows the η O E R ¯ as a function of θ O corresponding to the individual rate-limiting steps G 1 and G 3 .
Figure 5. (a) Snapshots of an example where all the intermediates were found to be adsorbed to the same binding site on the HEA surface. The central oxygen atoms of HO*, O*, and HOO* are highlighted in yellow. (b) Distributions of theoretical overpotentials ( η O E R ) for each θ O range. The η O E R values corresponding to different rate-limiting steps in the OER are depicted in blue, orange, and yellow, representing G 1 , G 2 , and G 3 , respectively. (c) Percentage of times that G 1 , G 2 , or G 3 were the rate-limiting steps as a function of θ O . (d) Average theoretical overpotential ( η O E R ¯ ) accounting for contributions from all rate-limiting steps. The error bars indicate the standard deviation of the total distributions shown in Panel (b). The inset shows the η O E R ¯ as a function of θ O corresponding to the individual rate-limiting steps G 1 and G 3 .
Nanomaterials 14 01058 g005
Table 1. Number of individual binding sites studied for calculating the theoretical overpotential (ηOER) at different coverage ranges. The table provides detailed information on the number of individual binding sites corresponding to either pristine binding site (Mode I) or oxidized binding site (Mode II) with oxidation step i (i = 1, 2, or 3) as the rate-limiting step.
Table 1. Number of individual binding sites studied for calculating the theoretical overpotential (ηOER) at different coverage ranges. The table provides detailed information on the number of individual binding sites corresponding to either pristine binding site (Mode I) or oxidized binding site (Mode II) with oxidation step i (i = 1, 2, or 3) as the rate-limiting step.
θ O
(ML)
Number of
Individual Sites
Number of Pristine SitesNumber of Oxidized Sites
Rate-Limiting StepRate-Limiting Step
G 1 G 2 G 3 G 1 G 2 G 3
[ 0.25 ,   0.5 ) 19N/A1126N/AN/A
[ 0.5 ,   0.75 ) 312113141N/A
[ 0.75 ,   1 ) 42N/AN/A138N/A3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yuan, G.; Ruiz Pestana, L. The Effect of Surface Oxygen Coverage on the Oxygen Evolution Reaction over a CoFeNiCr High-Entropy Alloy. Nanomaterials 2024, 14, 1058. https://doi.org/10.3390/nano14121058

AMA Style

Yuan G, Ruiz Pestana L. The Effect of Surface Oxygen Coverage on the Oxygen Evolution Reaction over a CoFeNiCr High-Entropy Alloy. Nanomaterials. 2024; 14(12):1058. https://doi.org/10.3390/nano14121058

Chicago/Turabian Style

Yuan, Geng, and Luis Ruiz Pestana. 2024. "The Effect of Surface Oxygen Coverage on the Oxygen Evolution Reaction over a CoFeNiCr High-Entropy Alloy" Nanomaterials 14, no. 12: 1058. https://doi.org/10.3390/nano14121058

APA Style

Yuan, G., & Ruiz Pestana, L. (2024). The Effect of Surface Oxygen Coverage on the Oxygen Evolution Reaction over a CoFeNiCr High-Entropy Alloy. Nanomaterials, 14(12), 1058. https://doi.org/10.3390/nano14121058

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop