Abstract
In contrast to pseudobinary alloys, the relative number of bonds in quaternary alloys cannot be determined uniquely from the composition. Indeed, we do not know if the Ga0.5In0.5As0.5N0.5 alloy should be thought of as InAs+GaN or as InN+GaAs. We study the distribution of bonds using Monte Carlo simulation and find that the number of In-N and Ga-As bonds increases relative to random alloys. This quaternary-unique short range order affects the band structure: we calculate a blueshift of the band gap and predict the emergence of a broadband tail of localized states around the conduction band minimum.