Assessing the Phosphate Distribution in Bioactive Phosphosilicate Glasses by 31P Solid-State NMR and Molecular Dynamics Simulations

Sub-Nanometer-Range Structural Effects From Mg2+ Incorporation in Na-Based Borosilicate Glasses Revealed by Heteronuclear NMR and MD Simulations

Magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) experiments and molecular dynamics (MD) simulations were employed to investigate Na2O-B2O3-SiO2 and MgO-Na2O-B2O3-SiO2 glass structures up to ≈0.3 nm. This encompassed the {Na[p]}, {Mg[p]}, and {B[3], B[4]} speciations and the {Si, B[p], M[p]}-BO and {Si, B[p], M[p]}-NBO interatomic distances to the bridging oxygen (BO) and nonbridging oxygen (NBO) species, where the superscript indicates the coordination number. The MD simulations revealed the dominance of Mg[5] coordinations, as mirrored in average Mg2+ coordination numbers in the 5.2-5.5 range, which are slightly lower than those of the larger Na+ cation but with a narrower coordination distribution stemming from the higher cation field strength (CFS) of the smaller divalent Mg2+ ion. We particularly aimed to elucidate such Na+/Mg2+ CFS effects, which primarily govern the short-range structure but also the borosilicate (BS) glass network order, where both MD simulations and heteronuclear double-resonance 11B/29Si NMR experiments revealed a reduction of B[4]-O-Si linkages relative to B[3]-O-Si upon Mg2+-for-Na+ substitution. These effects were quantified and discussed in relation to previous literature on BS glasses, encompassing the implications for simplified structural models and descriptions thereof.

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)'.

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'.

Improving Molybdenum and Sulfur Vitrification in Borosilicate Nuclear Waste Glasses Using Phosphorus: Structural Insights from NMR

Reclaiming residual fissile materials from spent nuclear fuel rods results in highly radioactive waste that must be immobilized in durable materials to protect the biosphere from harmful environmental exposure. The glasses typically used for this purpose are less effective at incorporating highly charged cations than cations in lower oxidation states. Molybdenum and sulfur solubilities in model aluminoborosilicate nuclear glasses were studied using solid-state nuclear magnetic resonance (NMR) spectroscopy. The extent of molybdenum and sulfur incorporation was significantly improved by doping the glasses with phosphorus, a high-field-strength cation capable of competing with Mo⁶⁺ and S⁶⁺ for oxygen, resulting in more oxygen sharing and a higher degree of Mo and S integration into the glass network. The addition of 3 mol % each of molybdenum and sulfur to phosphorus-free glasses led to the separation of crystalline sodium molybdates and sodium sulfates, whereas the replacement of 5 mol % SiO₂ by P₂O₅ resulted in complete Mo and S incorporation into the glass, enhancing their respective loading limits by several times. NMR spectroscopy indicates that this improvement originates from the availability of chain phosphates as binding sites. This work demonstrates the value of using phosphate additives to improve Mo and S integration in nuclear glasses, providing insight into its specific structural role and foundational understanding for further development as a material design principle.

Quantitative determination of the phosphorus environment in lithium aluminosilicate glasses using solid-state NMR techniques

The phosphorus environment is determined quantitatively using ³¹P Magic Angle Spinning NMR constrained by ³¹P–²⁷Al Multiple-Quantum Coherence-based NMR techniques.

Role of the Short Distance Order in Glass Reactivity

In 2005, our group described for the first time the structural characterization at the atomic scale of bioactive glasses and the influence of the glasses’ nanostructure in their reactivity in simulated body fluids. In that study, two bioactive sol-gel glasses with composition 80%SiO₂–20%CaO and 80%SiO₂–17%CaO–3%P₂O₅ (in mol-%) were characterized by High-Resolution Transmission Electron Microscopy (HRTEM). Such characterization revealed unknown features of the glasses’ structure at the local scale that allowed the understanding of their different in vitro behaviors as a consequence of the presence or absence of P₂O₅. Since then, the nanostructure of numerous bioactive glasses, including melt-prepared, sol-gel derived, and mesoporous glasses, was investigated by HRTEM, Nuclear Magnetic Resonance (NMR) spectroscopy, Molecular Dynamics (MD) simulations, and other experimental techniques. These studies have shown that although glasses are amorphous solids, a certain type of short distance order, which greatly influences the in vitro and in vivo reactivity, is always present. This paper reviews the most significant advances in the understanding of bioactive glasses that took place in the last years as a result of the growing knowledge of the glasses’ nanostructure.

Structure–composition trends in multicomponent borosilicate-based glasses deduced from molecular dynamics simulations with improved B–O and P–O force fields

Improved B–O and P–O force fields provide accurate molecular dynamics simulations of multicomponent glasses.

Medium-Range Structural Organization of Phosphorus-Bearing Borosilicate Glasses Revealed by Advanced Solid-State NMR Experiments and MD Simulations: Consequences of B/Si Substitutions

The short and intermediate range structures of a large series of bioactive borophosphosilicate (BPS) glasses were probed by solid-state nuclear magnetic resonance (NMR) spectroscopy and atomistic molecular dynamics (MD) simulations. Two BPS glass series were designed by gradually substituting SiO₂ by B₂O₃ in the respective phosphosilicate base compositions 24.1Na₂O-23.3CaO-48.6SiO₂-4.0P₂O₅ ("S49") and 24.6Na₂O-26.7CaO-46.1SiO₂-2.6P₂O₅ ("S46"), the latter constituting the "45S5 Bioglass" utilized for bone grafting applications. The BPS glass networks are built by interconnected SiO₄, BO₄, and BO₃ moieties, whereas P exists mainly as orthophosphate anions, except for a minor network-associated portion involving P-O-Si and P-O-B^[4] motifs, whose populations were estimated by heteronuclear ³¹P{¹¹B} NMR experimentation. The high Na⁺/Ca²⁺ contents give fragmented glass networks with large amounts of nonbridging oxygen (NBO) anions. The MD-generated glass models reveal an increasing propensity for NBO accommodation among the network units according to BO₄ < SiO₄ < BO₃ ≪ PO₄. The BO₄/BO₃ intermixing was examined by double-quantum-single-quantum correlation ¹¹B NMR experiments, which evidenced the presence of all three BO₃-BO₃, BO₃-BO₄, and BO₄-BO₄ connectivities, with B^[3]-O-B^[4] bridges dominating. Notwithstanding that B^[4]-O-B^[4] linkages are disfavored, both NMR spectroscopy and MD simulations established their presence in these modifier-rich BPS glasses, along with non-negligible B^[4]-NBO contacts, at odds with the conventional structural view of borosilicate glasses. We discuss the relative propensities for intermixing of the Si/B/P network formers. Despite the absence of pronounced preferences for Si-O-Si bond formation, the glass models manifest subtle subnanometer-sized structural inhomogeneities, where SiO₄ tetrahedra tend to self-associate into small chain/ring motifs embedded in BO₃/BO₄-dominated domains.

Proton Environments in Biomimetic Calcium Phosphates Formed from Mesoporous Bioactive CaO-SiO2-P2O5 Glasses in Vitro: Insights from Solid-State NMR

When exposed to body fluids, mesoporous bioactive glasses (MBGs) of the CaO-SiO₂-P₂O₅ system develop a bone-bonding surface layer that initially consists of amorphous calcium phosphate (ACP), which transforms into hydroxy-carbonate apatite (HCA) with a very similar composition as bone/dentin mineral. Information from various ¹H-based solid-state nuclear magnetic resonance (NMR) experiments was combined to elucidate the evolution of the proton speciations both at the MBG surface and within each ACP/HCA constituent of the biomimetic phosphate layer formed when each of three MBGs with distinct Ca, Si, and P contents was immersed in a simulated body fluid (SBF) for variable periods between 15 min and 30 days. Directly excited magic-angle-spinning (MAS) ¹H NMR spectra mainly reflect the MBG component, whose surface is rich in water and silanol (SiOH) moieties. Double-quantum-single-quantum correlation ¹H NMR experimentation at fast MAS revealed their interatomic proximities. The comparatively minor H species of each ACP and HCA component were probed selectively by heteronuclear ¹H-³¹P NMR experimentation. The initially prevailing ACP phase comprises H₂O and "nonapatitic" HPO₄^2-/PO₄^3- groups, whereas for prolonged MBG soaking over days, a well-progressed ACP → HCA transformation was evidenced by a dominating O¹H resonance from HCA. We show that ¹H-detected ¹H → ³¹P cross-polarization NMR is markedly more sensitive than utilizing powder X-ray diffraction or ³¹P NMR for detecting the onset of HCA formation, notably so for P-bearing (M)BGs. In relation to the long-standing controversy as to whether bone mineral comprises ACP and/or forms via an ACP precursor, we discuss a recently accepted structural core-shell picture of both synthetic and biological HCA, highlighting the close relationship between the disordered surface layer and ACP.

Structure–composition relationships of bioactive borophosphosilicate glasses probed by multinuclear 11B, 29Si, and 31P solid state NMR

By combining ¹¹B, ²⁹Si, and ³¹P nuclear magnetic resonance (NMR) experimental results, we present the first comprehensive structural investigation of 15 borophosphosilicate (BPS) glasses of the Na₂O–CaO–B₂O₃–SiO₂–P₂O₅ system.

Applications of DQ-DRENAR for the structural analysis of phosphate glasses

A new solid state NMR technique entitled DQ-DRENAR (Double-Quantum based Dipolar Recoupling Effects Nuclear Alignment Reduction) has been recently described for measuring homonuclear dipole-dipole interactions in multi-spin-1/2 systems under magic-angle spinning conditions. As in rotational echo double resonance (REDOR), the homonuclear dipole-dipole coupling constant can be extracted from a plot of a normalized difference signal (S0-S')/S0 versus dipolar mixing time, where S is the signal amplitude with the DQ-Hamiltonian present, and S0 is the signal amplitude in the absence of the DQ-Hamiltonian, which is used for normalization. Within the range of (S0-S)/S0≤0.3-0.5 such "homonuclear REDOR curves" can be approximated by simple parabolae, yielding effective squared dipole-dipole coupling constants ∑bjk(2) summed over all the pairwise interactions present. The effect of glassy disorder has been studied by simulations, replacing singular-valued internuclear distances by Gaussian distance distributions with the same central value. This situation results in a systematic over-estimation effect, which tends to compensate the implicit under-estimation effect caused by the parabolic fitting approach. The present contribution describes applications to a number of phosphate-based glasses and glass ceramics. The method turns out to be well suited for the differentiation of the various Q((n)) phosphate species, for characterizing the spatial distribution of isolated orthophosphate ions and for the detection of incipient nano-segregation and/or phase separation effects in glass ceramics.

Na/Ca Intermixing around Silicate and Phosphate Groups in Bioactive Phosphosilicate Glasses Revealed by Heteronuclear Solid-State NMR and Molecular Dynamics Simulations

We characterize the intermixing of network-modifying Na(+)/Ca(2+) ions around the silicate (QSi(n)) and phosphate (QP(n)) tetrahedra in a series of 16 Na2O–CaO–SiO2–P2O5 glasses, whose P content and silicate network connectivity were varied independently. The set includes both bioactive and bioinactive compositions and also encompasses two soda-lime-silicate members devoid of P, as well as two CaO–SiO2 glasses and one Na2O–SiO2–P2O5 glass. The various Si/P↔Na/Ca contacts were probed by molecular dynamics (MD) simulations together with heteronuclear magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) experimentation utilizing (23)Na{(31)P} and (23)Na{(29)Si} REDOR, as well as (31)P{ (23)Na} and (29)Si{(23)Na} REAPDOR. We introduce an approach for quantifying the extent of Na(+)/Ca(2+) ordering around a given QP(n) or QSi(n) group, encoded by the preference factor 0⩽ PM ⩽ 1 conveying the relative weights of a random cation intermixing (PM = 0) and complete preference/ordering (PM = 1) for one of the species M, which represents either Na(+) or Ca(2+). The MD-derived preference factors reveal phosphate and silicate species surrounded by Na(+)/Ca(2+) ions intermixed nearly randomly (PM ≲ 0.15), except for the QSi(4) and QSi(1) groups, which manifest more significant cation ordering with preference for Na+ and Ca2+, respectively. The overall weak preferences are essentially independent of the Si and P contents of the glass, whereas PM primarily correlates with the total amount of network modifiers: as the latter is increased, the Na/Ca distribution around the {QP(0), QSi(1), QSi(2)} groups with preference for Ca2(+ )tend to randomize (i.e., PCa decreases), while the PNa-values grow slightly for the {QP(1), QSi(3), QSi(4)} species already preferring coordination of Na. The set of experimental preference factors {PCa} for the orthophosphate (QP(0)) groups extracted from (31)P{(23)Na} REAPDOR NMR-derived M2(P–Na) dipolar second moments agrees well with the MD-generated counterparts. Our results on the Na/Ca intermixing in soda-lime-silicate glasses are discussed in relation to previous reports, highlighting the dependence of the conclusion on the approach to data evaluation.

Two heteronuclear dipolar results at the price of one: quantifying Na/P contacts in phosphosilicate glasses and biomimetic hydroxy-apatite

The analysis of S{I} recoupling experiments applied to amorphous solids yields a heteronuclear second moment M(2)(S-I) that represents the effective through-space dipolar interaction between the detected S spins and the neighboring I-spin species. We show that both M(2)(S-I) and M(2)(I-S) values are readily accessible from a sole S{I} or I{S} experiment, which may involve either S or I detection, and is naturally selected as the most favorable option under the given experimental conditions. For the common case where I has half-integer spin, an I{S} REDOR implementation is preferred to the S{I} REAPDOR counterpart. We verify the procedure by (23)Na{(31)P} REDOR and (31)P{(23)Na} REAPDOR NMR applied to Na(2)O-CaO-SiO(2)-P(2)O(5) glasses and biomimetic hydroxyapatite, where the M(2)(P-Na) values directly determined by REAPDOR agree very well with those derived from the corresponding M(2)(Na-P) results measured by REDOR. Moreover, we show that dipolar second moments are readily extracted from the REAPDOR NMR protocol by a straightforward numerical fitting of the initial dephasing data, in direct analogy with the well-established procedure to determine M(2)(S-I) values from REDOR NMR experiments applied to amorphous materials; this avoids the problems with time-consuming numerically exact simulations whose accuracy is limited for describing the dynamics of a priori unknown multi-spin systems in disordered structures.

Composition-dependent in vitro apatite formation at mesoporous bioactive glass-surfaces quantified by solid-state NMR and powder XRD

Solid-state NMR and powder XRD are employed to quantify the ACP (amorphous calcium phosphate) and HCA (hydroxy-carbonate apatite) components grown from three mesoporous bioactive glasses with variable compositions.