Thermodynamic modeling of the Si-P and Si-Fe-P systems was performed using the CALculation of PHAse Diagram (CALPHAD) method based on critical evaluation of available experimental data in the literature. The liquid and solid solutions were described using the Modified Quasichemical Model accounting for the short-range ordering and Compound Energy Formalism considering the crystallographic structure, respectively. In the present study, the phase boundaries for the liquidus and solid Si phases of the Si-P system were reoptimized. Furthermore, the Gibbs energies of the liquid solution, (Fe) (P,Si) , (Fe) (P,Si) , and (Fe) (P,Si) solid solutions and FeSi P compound were carefully determined to resolve the discrepancies in previously assessed vertical sections, isothermal sections of phase diagrams, and liquid surface projection of the Si-Fe-P system. These thermodynamic data are of great necessity for a sound description of the entire Si-Fe-P system. The optimized model parameters from the present study can be used to predict any unexplored phase diagrams and thermodynamic properties within the Si-Fe-P alloys.