Mixed-basis cluster expansion for thermodynamics of bcc alloys

To predict the ground-state structures and finite-temperature properties of an alloy, the total energies of many different atomic configurations sigmaequivalent to{sigma(i);i=1,...,N}, with N sites i occupied by atom A (sigma(i)=-1), or B (sigma(i)=+1), must be calculated accurately and rapidly. Direct local-density approximation (LDA) calculations provide the required accuracy, but are not practical because they are limited to small cells and only a few of the 2(N) possible configurations. The "mixed-basis cluster expansion" (MBCE) method allows to parametrize LDA configurational energetics E-LDA[sigma(i);i=1,...,N] by an analytic functional E-MBCE[sigma(i);i=1,...,N]. We extend the method to bcc alloys, describing how to select N-sigma ordered structures (for which LDA total energies are calculated explicitly) and N-F pair and multibody interactions, which are fit to the N-sigma energies to obtain a deterministic MBCE mapping of LDA. We apply the method to bcc Mo-Ta. This system reveals an unexpectedly rich ground-state line, pitting Mo-rich (100) superlattices against Ta-rich complex structures. Predicted finite-T properties such as order-disorder temperatures, solid-solution short-range order and the random alloy enthalpy of mixing are consistent with experiment.