Extended x-ray-absorption fine-structure experiments have previously demonstrated that for each composition x, the sample average of all nearest-neighbor A-C distances in an A1-xBxC semiconductor alloy is closer to the values in the pure (x→0) AC compound than to the composition-weighted (virtual) lattice average. Such experiments do not reveal, however, the distribution of atomic positions in an alloy, so the principle displacement directions and the degrees of correlation among such atomic displacements remain unknown. Here we calculate both structural and thermodynamic properties of Ga1-xInxP alloys using an explicit occupation- and position-dependent energy functional. The latter is taken as a modified valence force field, carefully fit to structural energies determined by first-principles local-density calculations. Configurational and vibrational degrees of freedom are then treated via the continuous-space Monte Carlo approach. We find good agreement between the calculated and measured mixing enthalpy of the random alloy, nearest-neighbor bond lengths, and temperature-composition phase diagram. In addition, we predict yet unmeasured quantities such as (a) distributions, fluctuations, and moments of first- and second-neighbor bond lengths as well as bond angles, (b) radial distribution functions, (c) the dependence of short-range order on temperature, and (d) the effect of temperature on atomic displacements. Our calculations provide a detailed picture of how atoms are arranged in substitutionally random but positionally relaxed alloys, and offer an explanation for the effects of site correlations, static atomic relaxations, and dynamic vibrations on the phase-diagram and displacement maps. We find that even in a chemically random alloy (where sites are occupied by Ga or In according to a coin toss), there exists a highly correlated static position distribution whereby the P atoms are displaced deterministically in certain high-symmetry directions. © 1995 The American Physical Society.