Estimation of Activity and Molar Excess Gibbs Energy of Binary Liquid Alloys Al-Cu, Al-Ni, and Al-Fe from the Partial Radial Distribution Function Simulated by Ab Initio Molecular Dynamics

Abstract

To accurately and conveniently obtain the thermodynamic data of binary liquid alloys, a new method is proposed in this study. It combines ab initio molecular dynamics (AIMD) simulation with a thermodynamic model to estimate the activity and molar excess Gibbs energy of binary liquid alloys. Additionally, two methods of grouping the partial radial distribution function (PRDF) of 5000 steps obtained by simulation are proposed for the first time. The PRDF of Al50Cu50, Al50Ni50, and Al50Fe50 is obtained by AIMD simulation. These PRDF are combined with four thermodynamic models to estimate the activity and molar excess Gibbs energy. Furthermore, the estimation results of the four models are compared with those of the Miedema model. The results show that when the first peak of the PRDF is obtained by the symmetric method, the average relative deviation (ARD) of the activity and molar excess Gibbs energy of the four models are, respectively: 28% and 32% for Molecular Interaction Volume Model (MIVM); 162% and 38% for Regular Solution Model (RSM); 508% and 65% for Wilson model; 562% and 67% for Non-Random Two-Liquid (NRTL). When the first peak of PRDF is obtained by non-symmetric method, the average ARD of the activity and molar excess Gibbs energy of the four models are, respectively: 64% and 20% for MIVM; 115% and 26% for RSM; 661% and 70% for Wilson; 727% and 72% for NRTL. In addition, the average ARD of the activity and molar excess Gibbs energy of the Miedema model are 113% and 33%. These data indicate that the estimation performance of the MIVM model is superior to the other four models, and the symmetric method performs better than the non-symmetric method. The grouping treatment of PRDF data effectively improves estimation performance.