The solubility limit of elements plays a crucial role in alloy designing, as well as in the understanding of kinetic processes in alloys. Using first-principles calculations, we calculated the ground states of Fe-RE compound, and predicted the stable and metastable structures of Fe-Y and Fe-Ce. The results indicate that the stable structures for Fe-Y compounds are Fe12Y.tI26, Fe17Y2.hP38 and Fe2Y.cF24, and those for Fe-Ce binary are Fe17Ce2.hP38, Fe19Ce5.hR24 and Fe2Ce.cF24. Based on the statistical-thermodynamic theory of dilute lattice gas, we presented a systematic series of calculations of the solubility limits of RE in α-Fe as a function of temperature. It is found that the solubilities are sorted in the following order: SLa>SCe>SY, corresponding to the trend in the solubility formation enthalpies: Hsol(La)< Hsol(Ce)< Hsol(Y).