At the end of Earth's accretion and after the core-mantle segregation, the existence of a basal magma ocean at the top of the CMB depends on the physical properties of mantle materials at relevant pressure and temperature. Present-day deep mantle structures such as ultralow-velocity zones (ULVZs) and low-shear velocity provinces (LLSVPs) might be directly linked to the still ongoing crystallization of a primordial magma ocean. We provide the first steps towards a self-consistent thermodynamic model of magma ocean crystallization at high-pressure. We build a solid-liquid thermodynamic database for silicates in the MgO-FeO-SiO2 system from 20 GPa to 140 GPa. We use already published chemical potentials for solids, liquid MgO and SiO2. We derive standard state chemical potential for liquid FeO and mixing relations from various indirect observations. Using this database, we compute the ternary phase diagram in the MgO-FeO-SiO2 system as a function of temperature and pressure. We confirm that the melt is lighter than the solid of same composition for all mantle conditions but at thermodynamic equilibrium, the iron-rich liquid is denser than the solid in the deep mantle. We compute a whole fractional crystallization sequence of the mantle and show that an iron rich and fusible layer should be left above the CMB at the end of the crystallization.