Jamieson K. Christie

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

Phosphate-based bioactive glasses containing fluoride ions offer the potential of a biomaterial which combines the bioactive properties of the phosphate glass and the protection from dental caries by fluoride. We conduct accurate first-principles molecular dynamics simulations of two compositions of fluorinated phosphate-based glass to assess its suitability as a biomaterial. There is a substantial amount of F-P bonding and as a result the glass network will be structurally homogeneous on medium-range length scales, without the inhomogeneities which reduce the bioactivity of other fluorinated bioactive glasses. We observe a decrease in the network connectivity with increasing F content, caused by the replacement of bridging oxygen atoms by non-bridging fluorine atoms, but this decrease is small and can be opposed by an increase in the phosphate content. We conclude that the structural changes caused by the incorporation of fluoride into phosphate-based glasses will not adversely affect their bioactivity, suggesting that fluorinated phosphate glasses offer a superior alternative to their silicate-based counterparts.

Nanoscale Chains Control the Solubility of Phosphate Glasses for Biomedical Applications

Bioactive phosphate-based glasses (PBGs) have several possible biomedical applications because of the chemical reactions they undergo with their surroundings when implanted into the body. The dissolution rate of PBGs in physiological conditions is a crucial parameter for these applications, to ensure, e.g., delivery of drugs or nutrients to the body at the correct rate. While it has been well-known that increasing the CaO content of these glasses at the expense of Na2O slows the dissolution rate, this paper provides an atomistic explanation of this for the first time. In this work, molecular dynamics simulations of five ternary P2O5-CaO-Na2O glasses reveal the structural properties at the atomic level that enhance the durability of PBGs as more Ca is added: (i) Ca binds together more fragments of the phosphate glass network than Na, (ii) Ca binds together more PO4 tetrahedra than Na, and (iii) Ca has a lower concentration of intratetrahedral phosphate bonding than Na. This behavior is rooted in the calcium ions higher charge and field strength. These results open the path to precise control and optimization of the PBG dissolution rate for specific biomedical applications.

Integrating biological activity into radioisotope vectors: molecular dynamics models of yttrium-doped bioactive glasses

Yttrium-containing bioactive glasses (YBGs) have the potential to extend the function of inert glass vectors for in situ radiotherapy, integrating it with the biodegradability and the ability to stimulate growth of new tissues. Combining these properties with the required low radioisotope dissolution rates requires a detailed understanding of the effect of yttrium on the structure of these materials. Molecular dynamics simulations of a range of YBG compositions with favourable network connectivity allowed us to extract relevant structural descriptors which affect the performances of these glasses in radiotherapy. We focus in particular on the density and strength of the network of non-covalent crosslinks built by Na, Ca and Y cations and on the ion clustering behaviour; we show that it should be possible to employ YBGs containing higher amounts of yttrium as radioisotope vectors, due to the increasingly stronger glass network which will limit the overall ion dissolution rates, without inhibiting the glass bioactivity.

Molecular dynamics simulations and structural descriptors of radioisotope glass vectors for in situ radiotherapy.

The low solubility (high durability) of yttrium aluminosilicate (YAS) glass is one of its most important properties for use in in situ radiotherapy. Simple parameters, such as silica or yttria content or network connectivity, are not sufficient to rationalize the dependence of the solubility on the glass composition observed experimentally. We performed classical molecular dynamics (MD) simulations of eight different YAS glasses of known solubility and analyzed the MD trajectories to identify specific structural features that are correlated and can be used to predict the solubility. We show that the (Si-)O-Si coordination number CN(SiOSi), the yttrium-yttrium clustering ratio R(YY), and the number of intratetrahedral O-Si-O bonds per yttrium atom N(intra) can be combined into a single structural descriptor s = f(CN(SiOSi),R(YY),N(intra)) with a high correlation with the solubility. The parameter s can thus be calculated from MD simulations and used to predict the solubility of YAS compositions, allowing one to adjust them to the range required by radiotherapy applications. For instance, its trend shows that high-silica- and low-yttria-content YAS glasses should be sufficiently durable for the radiotherapy application, although additional clinical considerations may set a lower limit to the yttria content.

Modeling the atomic structure of very high-density amorphous ice

The structure of very high-density amorphous (VHDA) ice has been modelled by positionally disordering three crystalline phases, namely ice IV, VI and XII. These phases were chosen because only they are stable or metastable in the region of the ice phase diagram where VHDA ice is formed, and their densities are comparable to that of VHDA ice. An excellent fit to the medium range of the experimentally observed pair-correlation function g(r) of VHDA ice was obtained by introducing disorder into the positions of the H2O molecules, as well as small amounts of molecular rotational disorder, disorder in the O--H bond lengths and disorder in the H--O--H bond angles. The low-k behaviour of the experimental structure factor, S(k), is also very well reproduced by this disordered-crystal model. The fraction of each phase present in the best-fit disordered model is very close to that observed in the probable crystallization products of VHDA ice. In particular, only negligible amounts of ice IV are predicted, in accordance with experimental observation.

Properties of water confined in hydroxyapatite nanopores as derived from molecular dynamics simulations

Abstract Bone tissue is characterized by nanopores inside the collagen-apatite matrix where fluid can exist and flow. The description of the fluid flow within the bone has however mostly relied on a macroscopic continuum mechanical treatment of the system, and, for this reason, the role of these nanopores has been largely overlooked. However, neglecting the nanoscopic behaviour of fluid within the bone volume could result in large errors in the overall description of the dynamics of fluid. In this work, we have investigated the nanoscopic origin of fluid motion by conducting atomistic molecular dynamics simulations of water confined between two parallel surfaces of hydroxyapatite (HAP), which is the main mineral phase of mammalian bone. The polarizable core–shell interatomic potential model used in this work to simulate the HAP–water system has been extensively assessed with respect to ab initio calculations and experimental data. The structural (pair distribution functions), dynamical (self-diffusion coefficients) and transport (shear viscosity coefficients) properties of confined water have been computed as a function of the size of the nanopore and the temperature of the system. Analysis of the results shows that the dynamical and transport properties of water are significantly affected by the confinement, which is explained in terms of the layering of water on the surface of HAP as a consequence of the molecular interactions between the water molecules and the calcium and phosphate ions at the surface. Using molecular dynamics simulations, we have also computed the slip length of water on the surface of HAP, the value of which has never been reported before.

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

First-principles and classical molecular dynamics simulations of undoped and silver-doped phosphate-based glasses with 50 mol% P2O5, 0-20 mol% Ag2O, and varying amounts of Na2O and CaO have been carried out. Ag occupies a distorted local coordination with a mean Ag-O bond length of 2.5 Å and an ill-defined first coordination shell. This environment is shown to be distorted octahedral/trigonal bipyramidal. Ag-O coordination numbers of 5.42 and 5.54-5.71 are calculated for first-principles and classical methodologies respectively. A disproportionation in the medium-range phosphorus Q(n) distribution is explicitly displayed upon silver-doping via CaO substitution, approximating 2Q(2)→Q(1) + Q(3), but not on silver-doping via Na2O substitution. An accompanying increase in FWHM of the phosphorus to bridging oxygen partial pair-correlation function is strong evidence for a bulk structural mechanism associated with decreased dissolution rates with increased silver content. Experimentally, Ag2O ↔ Na2O substitution is known to decrease dissolution and we show this to be a result of Ags local bonding.

Structural characteristics of positionally disordered lattices: Relation to the first sharp diffraction peak in glasses

Positional disorder has been introduced into the atomic structure of certain crystalline lattices, and the orientationally-averaged structure factor S(k) and pair-correlation function g(r) of these disordered lattices have been studied. Analytical expressions for S(k) and g(r) for Gaussian positional disorder in 2D and 3D are confirmed with precise numerical simulations. These analytic results also have a bearing on the unsolved Gauss circle problem in mathematics. As the positional disorder increases, high-k peaks in S(k) are destroyed first, eventually leaving a single peak, that with the lowest-k value. The pair-correlation function for lattices with such high levels of positional disorder exhibits damped oscillations, with a period equal to the separation between the furthest-separated (lowest-k) lattice planes. The last surviving peak in S(k) is, for example for silicon and silica, at a wavevector nearly identical to that of the experimentally-observed first sharp diffraction peak (FSDP) in the amorphous phases of those materials. Thus, for these amorphous materials at least, the FSDP can be regarded as arising from scattering from atomic configurations equivalent to the single family of positionally-disordered local Bragg planes having the furthest separation.

Structural properties and phase transitions in a silica clathrate

Melanophlogite, a low-pressure silica polymorph, has been extensively studied at different temperatures and pressures by molecular dynamics simulations. While the high-temperature form is confirmed as cubic, the low-temperature phase is found to be slightly distorted, in agreement with experiments. With increasing pressure, the crystalline character is gradually lost. At 8 GPa, the radial distribution function is consistent with an amorphous state. Like pristine glass, the topology changes, plastic behavior, and permanent densification appear above ∼12 GPa, triggered by Si coordination number changes. We predict that a partial crystalline and amorphous sample can be obtained by recovering the sample from a pressure of ∼12-16 GPa.

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.