Unbiased density functional solutions of freezing in binary mixtures of hard or soft spheres

We investigated the freezing of equi-concentration binary hard or soft sphere mixtures for various size ratios, σ2/σ1, using density functional theory. The Grand Potential is minimized using an unbiased, discrete, real-space mesh that does not constrain the shape of the density, and, in many cases, leads to solutions qualitatively different from those using Gaussians and plane-waves. Besides the usual face-centered-cubic solid-solution phase for σ2/σ1≈1.0, we find a sublattice-melt phase for σ2/σ1=0.85–0.5 (where the small-sphere density is nonlocalized and multi-peaked) and the NaCl phase for σ2/σ1=0.45–0.35 (when the small-sphere density again sharpens). For a range of size ratios of soft sphere mixtures, we could not find stable nonuniform solutions. Preliminary calculations within a Modified-Weighted Density-Approximation suggest that such multiple-peaked solutions are not unique to a particular density functional theory.