1ST-PRINCIPLES CALCULATION OF TEMPERATURE-COMPOSITION PHASE-DIAGRAMS OF SEMICONDUCTOR ALLOYS

The three-dimensional spin-1/2 Ising model with multiple-site interactions provides a natural framework for describing the temperature-composition phase diagram of substitutional binary alloys. We have carried out a first-principles approach to this problem in the following way: (i) The total energy of an A1-xBx alloy in any given substitutional arrangement of A and B on a given lattice is expanded in a series of interaction energies {Jf(V)} of figures f. (ii) The Ns functions {Jf(V)} for figures f are determined as a function of volume V by equating the total energies of a set of Ns periodic structures, calculated in the local-density formalism, to a series expansion in Jf with known coefficients. The calculation includes in a natural way atomic relaxation and self-consistent charge transfer, hence providing a link between the electronic structure and the interaction energies which decide phase stability. (iii) The number Ns and range of the interaction energies needed in such an Ising description is determined by the ability of such cluster expansions to reproduce the independently calculated total energy of other structures. We find that this requires extending the expansion for zinc-blende-based alloys up to the fourth fcc neighbors. (iv) Using such a complete set of Ns interaction energies, {Jf}, we find approximate solutions to the corresponding Ising Hamiltonian within the cluster-variation method, retaining up to four-body and fourth-neighbor terms. A renormalization procedure, whereby distant-neighbor correlations are folded into a compact set of effective near-neighbor correlations, is used and tested against Monte Carlo solutions. This yields the phase diagram and thermodynamic properties. The set {Jf} is also used to perform a ground-state search of all stable structures. This identifies stable and metastable phases. This approach has been applied to five III-V pseudobinary alloys (Al1-xGaxAs, GaAs1-xPx, In1-xGaxP, In1-xGaxAs, and GaSb1-xAsx) and three II-VI pseudobinary alloys (Cd1-xHgxTe, Hg1-xZnxTe, and Cd1-xZnxTe). We have calculated (i) excess enthalpies, (ii) phase diagrams, (iii) clustering probabilities, and (iv) equilibrium bond lengths. We discuss in detail the chemical trends in these properties and offer a simple (-G) model which reveals the underlying physical factors controlling such trends. © 1990 The American Physical Society.