The fracture mechanical behavior simulation of calcium-deficient hydroxyapatite crystals by molecular dynamics and first-principles calculation

Natural hydroxyapatite provides certain strength and stiffness to biological bones, and most of the studies on the strength of bone tissues have been carried out on hydroxyapatite (HAP). However, the Ca/P ratio of hydroxyapatite in bones is actually about 1.50, and the natural hydroxyapatite belongs to calcium-deficient hydroxyapatite (CDHA) with Ca vacancy defects. Therefore, this work focused on the effect of Ca vacancy defects on CDHA crystals through investigating the generation and expansion of microcracks under uniaxial tensile loading by combining molecular dynamics and first principles method. A series of crystal models with different Ca vacancy ratios are constructed and find that Ca vacancies degrade the mechanical properties of hydroxyapatite. Meanwhile the fracture behavior of crystals is detailed and find that the cracks arise at vacancies and extend along the direction of vacancies. Also, first-principles calculation is performed to reveal the mechanism of crack formation in MD simulations. It is found that the decrease of Ca–O bonding of CDHA causes the decrease of the stability of the crystal structure by analyzing the DOS of HAP and CDHA, and the cracks originate from Ca vacancies. This work performs more realistic simulations of CDHA with calcium vacancy defects in actual bone tissue and directly reveals the development and progression of bone fragility at the nanometer scale.