A density functional model for the binary crystal of hard spheres with vacancies

We study the stability of a binary mixture of hard spheres in the crystalline state in which a small fraction of lattice sites in the solid structure are vacant. The optimum vacancy concentration is obtained by minimizing the free energy of the inhomogeneous solid state. We use the modified weighted density functional approximation to compute the free energy. The necessary input for the theory is the thermodynamic properties of the homogeneous state of the mixture and is obtained from the solutions of the corresponding Percus-Yevick integral equations for the binary system. We compute the free energy for the crystal having two kinds of ordered structures in which (i) both the species lie in a disordered manner on a single face-centered-cubic lattice and (ii) each of the two species lie on a separate cubic lattice. Our theoretical model obtains equilibrium vacancy fraction of O(10(-5)) near the freezing point in both cases. The vacancy concentration decreases with the increase of density in both cases.