Ab initio molecular dynamics simulations of structural changes associated with the incorporation of fluorine in bioactive phosphate glasses

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.

Review: understanding the properties of amorphous materials with high-performance computing methods

Amorphous materials have no long-range order in their atomic structure. This makes much of the formalism for the study of crystalline materials irrelevant, and so elucidating their structure and properties is challenging. The use of computational methods is a powerful complement to experimental studies, and in this paper we review the use of high-performance computing methods in the simulation of amorphous materials. Five case studies are presented to showcase the wide range of materials and computational methods available to practitioners in this field. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Unveiling the Effect of CaF2 on the Microstructure and Transport Properties of Phosphosilicate Systems

As an effective flux, CaF2 is beneficial in improving the fluidity of slag in the steel-making process, which is crucial for dephosphorization. To reveal the existence form and functional mechanism of CaF2 in phosphosilicate systems, the microstructures and transport properties of CaO-SiO2-CaF2-P2O5 quaternary slag systems are investigated by molecular dynamics simulations (MD) combined with experiments. The results demonstrate that the Si-O coordination number does not vary significantly with the increasing CaF2 content, but the P-O coordination number dramatically decreases. CaF2 has a minor effect on the single [SiO4] but makes the structure of the silicate system simple. On the contrary, F− ions could reduce the stability of P-O bonds and promoted the transformation of [PO4] to [PO3F], which is beneficial for making the P element-enriched phosphate network structure more aggregated. However, the introduction of CaF2 does not alter the tetrahedral character of the original fundamental structural unit. In addition, the results of the investigation of the transport properties show that the self-diffusion coefficients of each ion are positively correlated with CaF2 content and arranged in the order of F− > Ca2+ > O2− ≈ P5+ > Si4+. Due to CaF2 reducing the degree of polymerization of the whole melts, the viscosity decreases from 0.39 to 0.13 Pa·s as the CaF2 content increases from 0% to 20%. Moreover, the viscosity of the melt shows an excellent linear dependence on the structural parameters.

Phosphate/oxyfluorophosphate glass crystallization and its impact on dissolution and cytotoxicity

The role of fluorine in bioactive glasses is of interest due to the potential of precipitating fluorapatite, a phase with higher chemical resistance than the typical hydroxyapatite precipitated from oxide bioactive glasses. However, the introduction of fluorine in silicate bioactive glasses was found deleterious to the bioactivity of the glass. Here, phosphate glasses with the composition 75NaPO₃-(25-x) CaO-xCaF₂ (in mol%), with x = 0-20 and glass-ceramics were investigated to evaluate their potential as substitutes to the traditional silicate bioactive glass. An increase in CaF₂ substitution for CaO led to an increase in the glass solubility, due to an increase in highly soluble F(M)n species (where M is a cation) and to an increased polymerization of the phosphate network. Structural analysis reveals the formation of FP bonds, in addition to the F(M)n species, in the glass with the higher CaF₂ content. Furthermore, with heat treatment, CaF₂ crystals precipitate within the bulk in the newly developed glass, when x = 20. This bulk crystallization reduces the glass dissolution without compromising the precipitation of a reactive layer at the glass surface. Finally, in vitro cell tests were performed using MC3T3 pre-osteoblastic cells. While the substitution of CaF₂ for CaO led to an increased cytotoxicity, the controlled crystallization of the fluorine containing glasses decreased such cytotoxicity to similar values than traditional bioactive phosphate glass (x0). This study reports on new oxyfluorophosphate glass and glass-ceramics able, not only, to precipitate a Ca-P reactive layer but also to be processed into glass-ceramics with controlled crystal size, density and cellular activity. STATEMENT OF SIGNIFICANCE: Uncontrolled crystallization of bioactive glasses has negative effect on the materials' bioactivity. While in silicate glass the bioactivity is solely reduced, in phosphate glasses it is often completely suppressed. Furthermore, the need for fluorine containing bioactive glasses, not only for use in bone reconstruction but also in toothpaste as emerged. The addition of F in both silicate and phosphate has led to challenges due the lack of Si-F or P-F bonds, generally leading to a decrease in bioactivity. Here, we developed a bioactive invert phosphate glass where up to 20 mol% of CaO was replaced with CaF₂. In the new developed glasses, NMR demonstrated formation of P-F bonds. The content of fluorine was tailored to induce CaF₂ bulk crystallization. Overall an increase in F was associated with an increase network connectivity. In turns it led to an increased dissolution rate which was linked to a higher cytotoxicity. Upon (partial to full) surface crystallization of the F-free glass, the bioactivity (ability to form a reactive layer) was loss and the cytotoxicity again increased due to the rapid dissolution of one crystal phase and of the remaining amorphous phase. On another hand, the controlled bulk precipitation of CaF₂ crystals, in the F-containing glass, was associated with a reduced cytotoxicity. The new oxyfluorophosphate glass-ceramic developed is promising for application in the biomedical field.

Recent advances and future perspectives of sol–gel derived porous bioactive glasses: a review

Sol–gel derived bioactive glasses have been extensively explored as a promising and highly porous scaffold materials for bone tissue regeneration applications owing to their exceptional osteoconductivity, osteostimulation and degradation rates.

Strong-coupling induced damping of spin-echo modulations in magic-angle-spinning NMR: Implications for J coupling measurements in disordered solids

In the context of improving J coupling measurements in disordered solids, strong coupling effects have been investigated in the spin-echo and refocused INADEQUATE spin-echo (REINE) modulations of three- and four-spin systems under magic-angle-spinning (MAS), using density matrix simulations and solid-state NMR experiments on a cadmium phosphate glass. Analytical models are developed for the different modulation regimes, which are shown to be distinguishable in practice using Akaike's information criterion. REINE modulations are shown to be free of the damping that occurs for spin-echo modulations when the observed spin has the same isotropic chemical shift as its neighbour. Damping also occurs when the observed spin is bonded to a strongly-coupled pair. For mid-chain units, the presence of both direct and relayed damping makes both REINE and spin-echo modulations impossible to interpret quantitatively. We nonetheless outline how a qualitative comparison of the modulation curves can provide valuable information on disordered networks, possibly also pertaining to dynamic effects therein.

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.

Modelling the local atomic structure of molybdenum in nuclear waste glasses with ab initio molecular dynamics simulations

The nature of chemical bonding of molybdenum in high level nuclear waste glasses has been elucidated by ab initio molecular dynamics simulations. Two compositions, (SiO₂)_57.5-(B₂O₃)₁₀-(Na₂O)₁₅-(CaO)₁₅-(MoO₃)_2.5 and (SiO₂)_57.3-(B₂O₃)₂₀-(Na₂O)_6.8-(Li₂O)_13.4-(MoO₃)_2.5, were considered in order to investigate the effect of ionic and covalent components on the glass structure and the formation of the crystallisation precursors (Na₂MoO₄ and CaMoO₄). The coordination environments of Mo cations and the corresponding bond lengths calculated from our model are in excellent agreement with experimental observations. The analysis of the first coordination shell reveals two different types of molybdenum host matrix bonds in the lithium sodium borosilicate glass. Based on the structural data and the bond valence model, we demonstrate that the Mo cation can be found in a redox state and the molybdate tetrahedron can be connected with the borosilicate network in a way that inhibits the formation of crystalline molybdates. These results significantly extend our understanding of bonding in Mo-containing nuclear waste glasses and demonstrate that tailoring the glass composition to specific heavy metal constituents can facilitate incorporation of heavy metals at high concentrations.

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.

Simulations reveal the role of composition into the atomic-level flexibility of bioactive glass cements

Bioactive glass ionomer cements (GICs), the reaction product of a fluoro-alumino-silicate glass and polyacrylic acid, have been in effective use in dentistry for over 40 years and more recently in orthopaedics and medical implantation. Their desirable properties have affirmed GIC's place in the medical materials community, yet are limited to non-load bearing applications due to the brittle nature of the hardened composite cement, thought to arise from the glass component and the interfaces it forms. Towards helping resolve the fundamental bases of the mechanical shortcomings of GICs, we report the 1st ever computational models of a GIC-relevant component. Ab initio molecular dynamics simulations were employed to generate and characterise three fluoro-alumino-silicate glasses of differing compositions with focus on resolving the atomic scale structural and dynamic contributions of aluminium, phosphorous and fluorine. Analyses of the glasses revealed rising F-content leading to the expansion of the glass network, compression of Al-F bonding, angular constraint at Al-pivots, localisation of alumino-phosphates and increased fluorine diffusion. Together, these changes to the structure, speciation and dynamics with raised fluorine content impart an overall rigidifying effect on the glass network, and suggest a predisposition to atomic-level inflexibility, which could manifest in the ionomer cements they form.

Structure-solubility relationships in fluoride-containing phosphate based bioactive glasses

Structural role of fluoride on chemical dissolution behavior of bioactive phosphate glasses has been studied.

On the structure of biomedical silver-doped phosphate-based glasses from molecular dynamics simulations

First-principles and classical molecular dynamics simulations have been carried out on undoped and silver-doped phosphate-based glasses with 50 mol% P₂O₅, 0–20 mol% Ag₂O, and varying amounts of Na₂O and CaO.