n-type Mg3Sb2-based materials have become a top candidate for efficient thermoelectric applications within 300–700 K, due to its high band degeneracy, inherently high carrier mobility and low lattice thermal conductivity, as well as its advantages of less toxicity and abundance. Existing works showed that Mg3Bi2-alloying largely help ensure the exceptional performance, leaving a key issue to be uncovered on the primary mechanisms favoring or limiting the thermoelectric performance of Mg3Sb2-xBix alloys. Here we focus on the alloy composition dependent transport properties at various temperatures, with a large volume of experimental data. It is revealed that, with increasing x, the reduction in both inertial mass and lattice thermal conductivity is significantly beneficial, but the closure in band gap leads to a strong compensation due to the bipolar effect. Such a compromise between band structure and phonon scattering results in optimal Mg3Bi2-alloying concentrations to be about 50%–75% at 300 K, 50%–60% at 450 K and 50% at 600 K, which successfully guiding this work to realize extraordinary thermoelectric figure of merit at these temperatures.

DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 1000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an atomate workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. To demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-dU and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-pU value results in a significantly improved agreement between the computed lattice parameters and experiment.