In this study, we utilized machine learning interatomic potentials (MLIPs) to investigate the nucleation mechanisms of calcium phosphate, a critical component of bone and teeth. Our analysis encompassed the process from pre-nucleation stage to the growth of amorphous calcium phosphate (ACP) in solution. We observed fluctuations in free calcium ion concentration and tracked the formation of uniform clusters in the early nucleation phases, confirming the existence of pre-nucleation clusters (PNCs). The PNCs are characterized by the composition Ca