Adsorption of Oxygen to High Entropy Alloy Surfaces for up to 2 ML Coverage Using Density Functional Theory and Monte Carlo Calculations

Abstract

High-entropy alloys (HEAs) manufactured with refractory elements are candidates for high-temperature structural applications. To enhance our understanding of oxidation in these complex systems, the early stages of oxidation on the surface of Al 10 Nb 15 Ta 5 Ti 30 Zr 40 were studied using density functional theory and thermodynamic modeling. Surface slabs were generated from bulk configurations sampled from equilibrium using a multicell Monte Carlo method for phase prediction. The bulk structure was found to be a single BCC phase in good agreement with experimental observations. The oxygen adsorbed with a strong preference for sites with Ti and Zr and avoided sites with Nb−Al and Nb−Ta. The surface was shown to be highly reactive to oxygen, yielding a dominating oxygen coverage of two monolayers over the temperature range of 100 to 2600 K and an oxygen pressure range of 10 −30 to 10 5 bar. Recovering a clean surface slab was not achieved until pressures approached vacuum conditions and temperature exceeded 1900 K, demonstrating the difficulty of oxygen removal from the surface. Grand canonical Monte Carlo simulations showed that high Nb content in the top surface layer reduced the surface reactivity to incoming oxygen. Inward oxygen diffusion at low coverage was preferred in regions rich with Zr but slowed with the addition of Ti and Al. Diffusion rates drastically decreased at 1 ML, especially in the region rich with Ti and Zr, where strong metal−oxygen bonds were reported. Our results indicated that a high content of Ti and Zr increased the reactivity of the HEA surface to oxygen. The presence of Nb also enhanced the resistance to oxygen adsorption, especially when partnered with Al and Ta. Inward oxygen diffusion was likely to occur at low coverage in regions rich with Zr but could be protected against with the addition of Al and Ti. The limitations of the present work are discussed. This study may provide insights that assist with devising short-and long-term mitigation strategies against material degradation related to high-temperature oxidation.

A computational approach for mapping electrochemical activity of multi-principal element alloys

“…A state-of-the-art hybrid MC and DFT approach was employed in order to develop MPEA models sampled from equilibrium 30 . Herein, the ability to predict electrochemical activity of an MPEA using computer simulation was investigated.…”

The Design of Aluminum-Matrix Composites Reinforced with AlCoCrFeNi High-Entropy Alloy Nanoparticles by First-Principles Studies on the Properties of Interfaces

“…Al and E bulk HEA are the total energies of bulk Al and bulk HEA, γ Al and γ HEA are surface energies of Al and HEA layers with free surfaces, respectively. To figure out the interfacial energy, the surface energy γ sur , i.e., either γ Al or γ HEA , is calculated by using the equation below [40],…”

Rapid mapping of alloy surface phase diagrams via Bayesian evolutionary multitasking

“…Despite being challenging, a handful of endeavors have been made in the past to search for the most stable adsorbate-alloy configuration under a given condition. The majority of theoretical studies rely on Monte Carlo (MC) simulations to search for the stationary states of alloy catalysts in the presence of adsorbates [32][33][34][35][36][37][38][39][40] . Recently, evolutionary algorithm (EA) emerges as an alternative solution to search for the most stable adsorbate-alloy configuration [41][42][43] .…”