Clustering of fluoride and phosphate ions in bioactive glass from computer simulation.

In order to understand the nature of ionic clustering in bioactive glass compositions, computer simulation was used to model four different compositions of bioactive glass with various amounts of flouride and phosphate. Fluoride ions were chemically bonded only to sodium and calcium, creating regions rich in fluoride and modifiers, and fluoride clustering was seen to be present in all compositions. The majority of phosphate groups are present as orthophosphate and phosphate clustering is also seen, and shown to be stronger in compositions with a lower phosphate content. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

An investigation into the adsorption mechanism of n-butanol by ZIF-8: a combined experimental and ab initio molecular dynamics approach.

The zeolitic imidazolate framework, ZIF-8, has been shown by experimental methods to have a maximum saturation adsorption capacity of 0.36 g g-1 for n-butanol from aqueous solution, equivalent to a loading of 14 butanol molecules per unit cell or 7 molecules per sodalite ß-cage. Diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) shows the presence of hydrogen bonding between adsorbed butanol molecules within the cage; the presence of three different O-H stretching modes indicates the formation of butanol clusters of varying size. Ab initio molecular dynamics simulations show the formation of intermolecular hydrogen bonding between the butanol molecules, with an average hydrogen-bond coordination number of 0.9 after 15 ps simulation time. The simulations also uniquely demonstrate the presence of weaker interactions between the alcohol O-H group and the π-orbital of the imidazole ring on the internal surface of the cage during early stages of adsorption. The calculated adsorption energy per butanol molecule is -33.7 kJ mol-1, confirming that the butanol is only weakly bound, driven primarily by the hydrogen bonding. Solid-state MAS NMR spectra suggest that the adsorbed butanol molecules possess a reasonable degree of mobility in their adsorbed state, rather than being rigidly held in specific sites. 2D 13C-1H heteronuclear correlation (HETCOR) experiments show interactions between the butanol aliphatic chain and the ZIF-8 framework experimentally, suggesting that O-H interactions with the π-orbital are only short lived. The insight gained from these results will allow the design of more efficient ways of recovering and isolating n-butanol, an important biofuel, from low-concentration solutions.

Multi-layer Scaffolds of Poly(caprolactone), Poly(glycerol sebacate) and Bioactive Glasses Manufactured by Combined 3D Printing and Electrospinning.

Three-dimensional (3D) printing has been combined with electrospinning to manufacture multi-layered polymer/glass scaffolds that possess multi-scale porosity, are mechanically robust, release bioactive compounds, degrade at a controlled rate and are biocompatible. Fibrous mats of poly (caprolactone) (PCL) and poly (glycerol sebacate) (PGS) have been directly electrospun on one side of 3D-printed grids of PCL-PGS blends containing bioactive glasses (BGs). The excellent adhesion between layers has resulted in composite scaffolds with a Young's modulus of 240-310 MPa, higher than that of 3D-printed grids (125-280 MPa, without the electrospun layer). The scaffolds degraded in vitro by releasing PGS and BGs, reaching a weight loss of ~14% after 56 days of incubation. Although the hydrolysis of PGS resulted in the acidification of the buffer medium (to a pH of 5.3-5.4), the release of alkaline ions from the BGs balanced that out and brought the pH back to 6.0. Cytotoxicity tests performed on fibroblasts showed that the PCL-PGS-BGs constructs were biocompatible, with cell viability of above 125% at day 2. This study demonstrates the fabrication of systems with engineered properties by the synergy of diverse technologies and materials (organic and inorganic) for potential applications in tendon and ligament tissue engineering.

Atomic structure of Mg-based metallic glasses from molecular dynamics and neutron diffraction.

We use a combination of classical molecular dynamics simulation and neutron diffraction to identify the atomic structure of five different Mg-Zn-Ca bulk metallic glasses, covering a range of compositions with substantially different behaviour when implanted in vitro. There is very good agreement between the structures obtained from computer simulation and those found experimentally. Bond lengths and the total correlation function do not change significantly with composition. The zinc and calcium bonding shows differences between composition: the distribution of Zn-Ca bond lengths becomes narrower with increasing Zn content, and the preference for Zn and Ca to avoid bonding to themselves or each other becomes less strong, and, for Zn-Ca, transforms into a positive preference to bond to each other. This transition occurs at about the same Zn content at which the behaviour on implantation changes, hinting at a possible structural connection. A very broad distribution of Voronoi polyhedra are also found, and this distribution broadens with increasing Zn content. The efficient cluster packing model, which is often used to describe the structure of bulk metallic glasses, was found not to describe these systems well.

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

Phosphate-based bioactive glasses (PBGs) dissolve harmlessly in the body with a dissolution rate which depends sensitively on composition. This makes them proposed vectors for e.g. drug delivery, or other applications where an active component needs to be delivered at a therapeutically appropriate rate. Molecular dynamics (MD) simulations provide atomic-level structural information about PBG compositions. We review recent work to show that MD is an excellent tool to unravel the connections between the PBG glass composition, its atomic structure, and its dissolution rate, which can help to optimise PBGs for specific medical applications.

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 Ag's local bonding.

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 ion's higher charge and field strength. These results open the path to precise control and optimization of the PBG dissolution rate for specific biomedical applications.

Polarizable force field development and molecular dynamics study of phosphate-based glasses.

Molecular dynamics simulations of phosphate-based glasses P(2)O(5)-CaO-Na(2)O have been carried out, using an interatomic force field that has been parameterized to reproduce the structural and mechanical properties of crystalline phosphorus pentoxide, o(')(P(2)O(5))(∞) orthorhombic phase. Polarization effects have been included through the shell-model potential and formal charges have been used to aid transferability. A modification to the DL_POLY code (version 2.20) was used to model the high temperature shell dynamics. Structural characterizations of three biomedically applicative molar compositions, (P(2)O(5))(0.45)(CaO)(x)(Na(2)O)(0.55-x) (x = 0.30, 0.35, and 0.40), have been undertaken. Good agreement with available experimental and ab initio data is obtained. The simulations show that, dependent on composition, the phosphorus atoms are primarily bonded to two or three oxygens that in turn bridge to neighbouring phosphorus atoms. Na(+) and Ca(2+) modifiers are found to occupy a pseudo-octahedral bonding environment with mean oxygen coordination numbers of 6.55 and 6.85, respectively, across all compositions studied.

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.

Bioactive glasses as potential radioisotope vectors for in situ cancer therapy: investigating the structural effects of yttrium.

The incorporation of yttrium in bioactive glasses (BGs) could lead to a new generation of radionuclide vectors for cancer therapy, with high biocompatibility, controlled biodegradability and the ability to enhance the growth of new healthy tissues after the treatment with radionuclides. It is essential to assess whether and to what extent yttrium incorporation affects the favourable properties of the BG matrix: ideally, one would like to combine the high surface reactivity typical of BGs with a slow release of radioactive yttrium. Molecular Dynamics simulations show that, compared to a BG composition with the same silica fraction, incorporation of yttrium results in two opposing effects on the glass durability: a more fragmented silicate network (leading to lower durability) and a stronger yttrium-mediated association between separate silicate fragments (leading to higher durability). The simulations also highlight a high site-selectivity and some clustering of yttrium cations, which are likely linked to the observed slow rate of yttrium released from related Y-BG compositions. Optimisation of yttrium BG compositions for radiotherapy applications thus depends on the delicate balance between these effects.

Fluorine environment in bioactive glasses: ab initio molecular dynamics simulations.

Fluorinated bioactive glasses (FBGs) combine the antibacterial properties of fluorine with the biological activity of phosphosilicate glasses. Because their biomedical application depends on the release of fluorine, the detailed characterization of the fluorine environment in FBGs is the key to understand their properties. Car-Parrinello molecular dynamics (CPMD) simulations have been performed on a 45S5 Bioglass composition in which 10 mol % of the CaO has been replaced with CaF(2), and have allowed us to resolve some longstanding issues about the atomic structure of fluorinated bioglasses, with particular regard to the structural role of fluorine. F is coordinated almost entirely to the modifier ions Na and Ca, with a very small amount of residual Si-F bonds, whose fraction only becomes significant in the melt precursor. High temperature leads to Si-F bonds in both tetra- (SiO(3)F) and, less frequently, penta-coordinated (SiO(4)F and SiO(3)F(2)) complexes, showing that formation of these bonds through the expansion of the SiO(4) coordination shell is generally less favored. There is no evidence for preferential bonding of F to either modifier ion: almost all F atoms are coordinated to both calcium and sodium in a "mixed state", rather than exclusively to either, as had been conjectured. We discuss the consequences of these findings on the properties of fluorine-containing bioglasses.

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.