A first principles investigation of the oxygen adsorption on Zr(0001) surface using cluster expansions

Abstract

The design of corrosion resistant zircalloys is important for a variety of technological applications ranging from medicine to the nuclear industry. Since corrosion resistance is mainly attributed to the formation of a surface oxide layer, developing a detailed understanding of this process may assist in future corrosion resistance design. In this work, we conduct a systematic multi-scale investigation of the early stages of oxide formation. This was accomplished by first using a database of fully relaxed DFT calculations to build a cluster-expansion description of the potential function. The developed potential was reasonably good at predicting DFT energies as evidenced by the cross-validation score of 4.4meV/site. The effective cluster expansion parameters were indicative of repulsive adsorbate interactions in the adlayer in agreement with the literature. The potential then allowed for a systematic investigation of the oxygen configurations on the Zr(0001) surface via Monte Carlo simulations. The adsorption energy was recorded as a function of coverage and an increasing trend was observed in agreement with DFT predictions and the repulsive nature of interactions in the adlayer. The convex hull diagram was recorded indicating the most stable configuration to occur around a coverage of 0.6ML. The adsorption isotherm was also recorded and contrasted for two temperatures relevant for different applications.