We show how a configurational lattice dynamics technique, in which the free energy of a number of configurations is determined directly by means of a fully dynamic structural minimization, can be used to calculate thermodynamic properties of solid solutions and phase diagrams. No assumptions are made as to the nature of the solution and both configurational and vibrational entropy contributions are determined directly. Only a small number of configurations are required. We illustrate the method using MnO/MgO, for which our results support the recent experiments of Wood, Hackler, and Dobson ͓Contrib. Mineral. Petrol. 115, 438 ͑1994͔͒ who, in contrast to previous workers, suggest the formation of a complete solid solution at temperatures only above 1100 K.Solid solutions, alloys, and grossly nonstoichiometric compounds present considerable challenges to the theoretician, as does the calculation of absolute free energy. Ceramic solid solutions in particular are often strongly nonideal and approaches such as the cluster variation method ͑CVM͒, 1 widely used for metallic alloys, often perform poorly. In addition, despite the importance of accurate thermodynamic data for oxide solid solutions in such areas as ceramic fabrication and design and mineralogy, accurate thermodynamic data are rare.We have recently developed a highly efficient method for the fully dynamic structure optimization of large unit cells 2 that uses lattice statics and quasiharmonic lattice dynamics ͑QLD͒. We calculate the full set of free-energy first derivatives analytically, and a full minimization of the free energy with respect to all structural variables for large unit cells is possible. The accurate calculation of the free energy via QLD is quick and computationally efficient and does not resort to lengthy thermodynamic integration. Here we show this can be used for the free energies of solid solutions ͑in-cluding ⌬H mix and ⌬S mix ͒ and phase diagrams. The method is readily extended to elevated T and high P; no a priori assumptions are made regarding the configurational entropy contribution and the vibrational contribution is also evaluated straightforwardly.We illustrate our approach using MnO/MgO, for which not only are there several sets of experimental enthalpy data 3,4 but also conflicting reports of the phase diagram. As shown in Ref. 5, the experiments of Raghavan, Iyengar, and Abraham 6 suggest a consolute temperature, T c , as low as 600 K, whereas the results of Wood, Hackler, and Dobson 7 are consistent with a much larger T c ͑Ϸ1100 K͒ and a markedly asymmetric phase diagram. The data of Ref. 6 are indeed surprising given the CaO/MgO phase diagram, in which there is a large two-phase region, and the mismatch in ionic radii between Mn 2ϩ and Mg 2ϩ , which is substantial although smaller than between Ca 2ϩ and Mg 2ϩ . A further aim of this paper is to attempt to resolve this issue.In principle the solid solution can assume any state, i.e., each atom can be at any position. However, the only states of practical importance away from t...