Mg3(Sb1-xBix)2 alloy has been extensively studied in the last 5 years due to its exceptional thermoelectric (TE) performance. The absence of accurate force field for inorganic alloy compounds presents great challenges for computational studies. Here, we explore the atomic microstructure, thermal, and elastic properties of the Mg3(Sb1-xBix)2 alloy at different solution concentrations through atomic simulations with a highly accurate machine learning interatomic potential (ML-IAP). We find atomic local ordering in the optimized structure with the Bi-Bi pair inclined to join adjacent layers and Sb-Sb pair preferring to stay within the same layer. The thermal conductivity changes with the solution concentrations can be correctly predicted through ML-IAP-based molecular dynamics simulations. Spectral thermal conductance analysis shows that the continuous movement of low-frequency peak to high frequency is responsible for the reduction of the thermal conductivity upon alloying. Elastic calculations reveal that similar to the thermal conductivity, solid solution alloying can reduce the overall elastic properties at both Mg3Sb2 and Mg3Bi2 ends, while anisotropic behavior is clearly observed with linear interpolation relationship upon alloying along the interlayer direction and nonlinearity along the intralayer direction. Although the atomic local ordering shows little effects on the properties of the Mg3(Sb1-xBix)2 alloy with only two alloying elements, it possesses potential important impacts on multiprincipal element inorganic TE alloys. This work provides a recipe for computational studies on the TE alloy systems and thus can accelerate the discovery and optimization of TE materials with high TE performance.
Keywords: alloying design; atomistic simulation; machine-learning interatomic potential; thermal conductivity; zintl thermoelectric.