Analytic Model for U-Nb Liquidus and U-6Nb Melting Curve

Uranium–niobium (U-Nb) alloys, used in a variety of industrial and energy applications that require high density, ductility, and good corrosion resistance, comprise a highly complex, multiphasic system with a phase diagram well established through decades of extensive experimental and theoretical research. They are also one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. The key factor determining the performance and safety of a metallic fuel such as U-Nb is its operational limits in the application environment, which are closely related to material’s structure and thermodynamic stability. They are in turn closely related to the ambient (zero-pressure) melting point (Tm); thus, Tm is an important engineering parameter. However, the current knowledge of Tm of the U-Nb system is limited, as the only experimental study of its Nb-rich portion dates back to 1958. In addition, it has not yet been adequately modeled based on general thermodynamics principles or using an equation-of-state approach. In this study, we present a theoretical model for the melting curve (liquidus) of a mixture, and apply it to U-Nb, which is considered as a mixture of pure U and pure Nb. The model uses the known melting curves of pure constituents as an input and predicts the melting curve of their mixture. It has only one free parameter, which must be determined independently. The ambient liquidus of U-Nb predicted by the model appears to be in good agreement with the available experimental data. We calculate the melting curve (the pressure dependence of Tm) of pure U using ab initio quantum molecular dynamics (QMD), the knowledge of which is required for obtaining the model parameters for U. We also generalize the new model to nonzero pressure and consider the melting curve of U-6 wt.% Nb (U-6Nb) alloy as an example. The melting curve of U-6Nb alloy predicted by the model appears to be in good agreement with the ab initio melting curve obtained from our QMD simulations. We suggest that the U-18Nb alloy can be considered as a proxy for protactinium (Pa) and demonstrate that the melting curves of U-18Nb and Pa are in good agreement with each other.