Thermodynamic optimization of the binary rare earth alloy (Ce–La, Ce–Pr, Ce–Nd and La–Nd) systems was performed in this work through the CALPHAD method based on the critical evaluation of all available phase diagram and thermodynamic data reported in the literature. During the thermodynamic modeling, the solution phases including liquid, bcc, fcc and dhcp, were treated as the substitutional solution model. Thermodynamic parameters of the stable phases in the Ce–La, Ce–Pr, Ce–Nd and La–Nd binary systems are obtained finally and would be used directly to develop the thermodynamic database of the multi-component Nd–Fe–B-based alloys, which is indispensable for designing alloy compositions and processes of Nd–Fe–B permanent magnets with highly abundant rare earth metals.