The specific exposed surfaces of monoclinic LaPO4 was revealed with HRTEM. DFT calculations were used to calculate the surface structure and energies. The equilibrium shapes of monoclinic LaPO4 were predicted. The mechanism for the morphology evolution of LaPO4 was identified. Provide a theoretical method to calculate the morphologies of phosphates.