The density functional theory for the freezing hard spheres is studied. We use a variety of the hard sphere direct correlation functions (DCFs) such as the one introduced by Roth et al. [J. Phys. Condens. Matter14, 12063 (2002)]; we call it RELK DCF, a new hard sphere DCF developed here by a combination of the RELK and the Percus–Yevick DCFs, and finally the generalized mean spherical approximation (GMSA). The structure factor, the freezing, and order parameters are calculated using these DCFs. The structure factor obtained by the new DCF is in good agreement with the Monte Carlo simulation. The best result for the freezing parameters in comparison with the Monte Carlo simulations is obtained by using our new expression for the DCF. Finally we obtain the Helmholtz free energy of the hard sphere FCC crystals using modified weighted density approximation (MWDA), and again the best results are obtained by using the new expression for the hard sphere DCF.