Development of a ReaxFF reactive force field for ternary phosphate-based bioactive glasses

Phosphate-based glasses (PBGs) in the CaO-Na2O-P2O5 system have diverse applications as biomaterials due to their unique dissolution properties. A reactive force field (ReaxFF) has been developed to simulate these materials at the atomic level. The ReaxFF parameters of PBGs, including the interaction between phosphorus and calcium atoms, have been developed using a published code based on genetic algorithms. The training data, including the atomic charges, atomic forces, bond and angle parameters, and different differential energies, are chosen and measured from static quantum-mechanical calculations and ab initio molecular dynamics of PBGs. We did different short- and medium-range structural analyses on the bulk simulated PBGs with different compositions to validate the developed potential. Radial and angular distribution functions and coordination numbers of network formers and modifiers, as well as the network connectivity of the glass, are in agreement with experimental and previous simulations using both shell-model classical force fields and ab initio simulated data; for example, the coordination number of phosphorus is 4.0. This successful development of ReaxFF parameters being able to describe the bulk PBGs enables us to work on the dissolution behavior of the glasses, including the interaction of water molecules with PBGs, in future works.

Effect of strontium inclusion on the bioactivity of phosphate-based glasses

We have conducted first-principles and classical molecular dynamics simulations of various compositions of strontium-containing phosphate glasses, to understand how strontium incorporation will change the glasses' activity when implanted into the body (bioactivity). To perform the classical simulations, we have developed a new interatomic potential, which takes account of the polarizability of the oxygen ions. The Sr-O bond length is ∼2.44-2.49 Å, and the coordination number is 7.5-7.8. The Q n distribution and network connectivity were roughly constant for these compositions. Sr bonds to a similar number of phosphate chains as Ca does; based on our previous work (Christie et al. in J Phys Chem B 117:10652, 2013), this implies that SrO ↔ CaO substitution will barely change the dissolution rate of these glasses and that the bioactivity will remain essentially constant. Strontium could therefore be incorporated into phosphate glass for biomedical applications.

Structures and properties of phosphate-based bioactive glasses from computer simulation: a review

Computer simulations have enabled breakthroughs in understanding the connections between the atomic structure and properties of bioactive phosphate glasses.

Atomic structure of biodegradable Mg-based bulk metallic glass

We have used highly accurate first-principles molecular dynamics simulations to elucidate the structure of Mg60Zn35Ca5 and Mg72Zn23Ca5 bulk metallic glasses, which are candidate materials for biomedical implants; these two compositions exhibit different behaviours when implanted. The environments of each species are different, and average coordination numbers are ∼13 for Mg, ∼11 for Zn and ∼18-19 for Ca. A wide range of local environments were found and icosahedral motifs, often seen in bulk metallic glasses, were among the most common for both Mg and Zn. Through the computation of a chemical short-range order parameter, a moderate avoidance of Zn-Zn bonding over Zn-Mg or Zn-Ca was observed. No statistically significant difference in structure was observed between the two compositions.