Relatively resource-rich but property inferior-Ce-Fe-B magnet can be improved by partial replacement of Ce by Nd and/or Pr. In addition to the amount of Nd/Pr, their distribution profile in microstructure plays an important role. From our first principles density functional theory (DFT) calculation, the substitution energy of Ce by Pr/Nd is negative in Ce 2 Fe 14 B (2:14:1) while that for laves phase, CeFe 2 is positive, implying that Nd/Pr stabilize 2:14:1 and suppress the formation of the CeFe 2 phase. Micromagnetic simulation indicates that homogenized distribution of Nd/Pr improves squareness of demagnetization curve, while core (Ce-rich)-shell (Nd/Pr-rich) 2:14:1 grain structure enhances coercivity. Magnetic properties of Ce-Fe-B can be optimized by manipulating distribution profile of chemical element in microstructure based on their subtle difference in thermodynamic property, which is an effective pathway to design optimized chemical composition and processing route for high-performance magnet.