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.