Molecular Dynamics Simulation of Biomimetic Biphasic Calcium Phosphate Nanoparticles

Biphasic calcium phosphate (BCP) is used as a bone substitute and bone tissue repair material due to its better control over bioactivity and biodegradability. It is crucial to stabilize the implanted biomaterial while promoting bone ingrowth. However, a lack of standard experimental and theoretical protocols to characterize the physicochemical properties of BCP limits the optimization of its composition and properties. Computational simulations can help us better to learn BCP at a nanoscale level. Here, the Voronoi tessellation method was combined with simulated annealing molecular dynamics to construct BCP nanoparticle models of different sizes, which were used to understand the physicochemical properties of BCP (e.g., melting point, infrared spectrum, and mechanical properties). We observed a similar to 20 to 30 angstrom layer of calcium-deficient hydroxyapatite at the HAP/P-TCP interface due to particle migration, which may contribute to BCP stability. The BCP model may stimulate further research into BCP ceramics and multiphasic ceramics. Moreover, our study may facilitate the optimization of compositions of BCP-based biomaterials.