Study of nanoscale deformation mechanisms in bulk hexagonal hydroxyapatite under uniaxial loading using molecular dynamics

Exploring the hierarchical structure of lamellar bone and its impact on fracture behaviour: A computational study using a phase field damage model

Bone is a naturally occurring composite material composed of a stiff mineral phase and a compliant organic matrix of collagen and non-collagenous proteins (NCP). While diverse mineral morphologies such as platelets and grains have been documented, the precise role of individual constituents, and their morphology, remains poorly understood. To understand the role of constituent morphology on the fracture behaviour of lamellar bone, a damage based representative volume element (RVE) was developed, which considered various mineral morphologies and mineralised collagen fibril (MCF) configurations. This model framework incorporated a novel phase-field damage model to predict the onset and evolution of damage at mineral-mineral and mineral-MCF interfaces. It was found that platelet-based mineral morphologies had superior mechanical performance over their granular counterparts, owing to their higher load-bearing capacity, resulting from a higher aspect ratio. It was also found that MCFs had a remarkable capacity for energy dissipation under axial loading, with these fibrillar structures acting as barriers to crack propagation, thereby enhancing overall elongation and toughness. Interestingly, the presence of extrafibrillar platelet-based minerals also provided an additional toughening through a similar mechanism, whereby these structures also inhibited crack propagation. These findings demonstrate that the two primary constituent materials of lamellar bone play a key role in its toughening behaviour, with combined effect by both mineral and MCFs to inhibit crack propagation at this scale. These results have provided novel insight into the fracture behaviour of lamellar bone, enhancing our understanding of microstructure-property relationships at the sub-tissue level.

Nano-friction behavior and deformation study of hydroxyapatite in ultra-precision polishing process

CONTEXT

In order to study the effect of ultra-precision machining on the surface quality of hydroxyapatite semiconductor materials as well as the material removal mechanism of hydroxyapatite, the mechanical polishing behaviors of hydroxyapatite at different polishing depths were studied by molecular dynamics method. The results show that the subsurface damage of hydroxyapatite increases with increasing polishing depth. The polishing temperature and the polishing force showed a positive correlation with the polishing depth, and the variation of the polishing force was related to the accumulation-release effect of the potential energy of hydroxyapatite material. In addition, the variation of stresses in hydroxyapatite during polishing is mainly influenced by the thermal softening effect. With a smaller polishing depth, the hydroxyapatite semiconductor material has fewer structural defects, fewer atoms undergoing phase transitions, lower surface roughness, and better surface quality. Therefore, to ensure the long-lasting service life of hydroxyapatite semiconductor materials, a small polishing depth should be used in ultra-precision machining. Additionally, this study also provides a theoretical reference for future research on the mechanical properties of hydroxyapatite-based composites.

METHODS

A Large-Scale Atomic/Molecular Parallel Simulator (LAMMPS) was utilized to perform molecular dynamics simulations. The output was visualized and analyzed by the Open Visualization Tools (OVITO) software. The intermolecular interactions were described by the polymer consistent force-field and the 12/6 Lennard-Jones potential functions. The workpiece was polished under a micro-canonical ensemble with the temperature settled at 300 K. Periodic boundary conditions were adopted and the velocity-Verlet algorithm was used to integrate the atomic motion with a timestep of 0.1 femtoseconds (fs).

Towards Polycaprolactone-Based Scaffolds for Alveolar Bone Tissue Engineering: A Biomimetic Approach in a 3D Printing Technique

The alveolar bone is a unique type of bone, and the goal of bone tissue engineering (BTE) is to develop methods to facilitate its regeneration. Currently, an emerging trend involves the fabrication of polycaprolactone (PCL)-based scaffolds using a three-dimensional (3D) printing technique to enhance an osteoconductive architecture. These scaffolds are further modified with hydroxyapatite (HA), type I collagen (CGI), or chitosan (CS) to impart high osteoinductive potential. In conjunction with cell therapy, these scaffolds may serve as an appealing alternative to bone autografts. This review discusses research gaps in the designing of 3D-printed PCL-based scaffolds from a biomimetic perspective. The article begins with a systematic analysis of biological mineralisation (biomineralisation) and ossification to optimise the scaffold's structural, mechanical, degradation, and surface properties. This scaffold-designing strategy lays the groundwork for developing a research pathway that spans fundamental principles such as molecular dynamics (MD) simulations and fabrication techniques. Ultimately, this paves the way for systematic in vitro and in vivo studies, leading to potential clinical applications.

Effect of Pore Defects on Uniaxial Mechanical Properties of Bulk Hexagonal Hydroxyapatite: A Molecular Dynamics Study

Hydroxyapatite (HAP) is a calcium apatite bioceramic used in various naturally-derived and synthetic forms for bone repair and regeneration. While useful for the regrowth of osseus tissue, the poor load-bearing capacity of this material relative to other biomaterials is worsened by the propensity for pore formation during the synthetic processing of scaffolds, blocks, and granules. Here we use molecular dynamics (MD) simulations to improve the current understanding of the defect-altered uniaxial mechanical response in hexagonal HAP single crystals relative to defect-free structures. The inclusion of a central spherical pore within a repeated lattice was found to reduce both the failure stress and failure strain in uniaxial tension and compression, with up to a 30% reduction in maximum stress at the point of failure compared to a perfect crystalline structure observed when a 30 Å diameter pore was included. The Z axis ([0 0 0 1] crystalline direction) was found to be the least susceptible to pore defects in tension but the most sensitive to pore inclusion in compression. The deformation mechanisms are discussed to explain the observed mechanical responses, for which charge imbalances and geometric stress concentration factor effects caused by pore inclusion play a significant role.

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.

In silico study on probing atomistic insights into structural stability and tensile properties of Fe-doped hydroxyapatite single crystals

Hydroxyapatite (HA, Ca10PO4(OH)2) is a widely explored material in the experimental domain of biomaterials science, because of its resemblance with natural bone minerals. Specifically, in the bioceramic community, HA doped with multivalent cations (e.g., Mg2+, Fe2+, Sr2+, etc.) has been extensively investigated in the last few decades. Experimental research largely established the critical role of dopant content on mechanical and biocompatibility properties. The plethora of experimental measurements of mechanical response on doped HA is based on compression or indentation testing of polycrystalline materials. Such measurements, and more importantly the computational predictions of mechanical properties of single crystalline (doped) HA are scarce. On that premise, the present study aims to build atomistic models of Fe2+-doped HA with varying Fe content (10, 20, 30, and 40 mol%) and to explore their uniaxial tensile response, by means of molecular dynamics (MD) simulation. In the equilibrated unit cell structures, Ca(1) sites were found to be energetically favourable for Fe2+ substitution. The local distribution of Fe2+ ions significantly affects the atomic partial charge distribution and chemical symmetry surrounding the functional groups, and such signatures are found in the MD analyzed IR spectra. The significant decrease in the intensity of the IR bands found in the Fe-doped HA together with band splitting, because of the symmetry changes in the crystal structure. Another important objective of this work is to computationally predict the mechanical response of doped HA in their single crystal format. An interesting observation is that the elastic anisotropy of undoped HA was not compromised with Fe-doping. Tensile strength (TS) is systematically reduced in doped HA with Fe2+ dopant content and a decrease in TS with temperature can be attributed to the increased thermal agitation of atoms at elevated temperatures. The physics of the tensile response was rationalized in terms of the strain dependent changes in covalent/ionic bond framework (Ca-P distance, P-O bond strain, O-P-O angular strain, O-H bond distance). Further, the dynamic changes in covalent bond network were energetically analyzed by calculating the changes in O-H and P-O bond vibrational energy. Summarizing, the current work establishes our foundational understanding of the atomistic phenomena involved in the structural stability and tensile response of Fe-doped HA single crystals.

New Prospects in Nano Phased Co-substituted Hydroxyapatite Enrolled in Polymeric Nanofiber Mats for Bone Tissue Engineering Applications

The most common forms of tissue impairment are fracture bones and significant bone disorders caused by multiple traumas or normal aging. Surgical care sometimes necessitates the placement of a temporary or permanent prosthesis, which continues to be a challenge for orthopedic surgeons, including those with large bone defects. Electrospun scaffolds made from natural and synthetic nanofiber-based polymers are studied as natural extracellular matrix (ECM)-like scaffolds for tissue engineering. Besides, nanostructured materials have properties and functions depending on the scale of natural materials such as hydroxyapatite (HAP), ranging from 1 to 100 nm, which activity was proficient upon enrolled in nanofiber mats. The use of nanofibers in combination with nano-HAP has increased the scaffold's ability to replicate the construction of natural bone tissue that is the aim of the present text. In bone engineering, nanofiber substrates facilitate cell adhesion, proliferation, and differentiation, while HAP induces cells to secrete ECM for bone mineralization and development. This review aims to draw the reader's attention to the critical issues with synthetic and natural polymers containing HAP in bone tissue engineering; co-substituted hydroxyapatite has also been mentioned.