Structural differences and phase separation in alkali silicate glasses

A modified random network model for P2O5-Na2O-Al2O3-SiO2 glass studied by molecular dynamics simulations

We investigated the short- and medium-range structural features of sodium aluminosilicate glasses with various P₂O₅ (0–7 mol%) content and Al/Na ratios ranging from 0.667 to 2.000 by using molecular dynamics simulations. The local environment evolution of network former cations (Si, Al, P) and the extent of clustering behavior of modifiers (Na⁺) is determined through pair distribution function (PDF), total correlation function (TDF), coordination number (CN), Q_xⁿ distribution and oxygen speciation analysis. We show that Al–O–P and Si–O–Al linkage is preferred over other connections as compared to a random model and that Si–O–Si linkage is promoted by the P₂O₅ addition, which is related to structural heterogeneity and generates well-separated silicon-rich and aluminum–phosphorus-rich regions. Meanwhile, due to the relatively high propensity of Al to both Si and P, heterogeneity can be partly overcome with high Al content. A small amount of Si–O–P linkages have been detected at the interface of separated regions. Clustering of Na⁺ is also observed and intensified with the addition of P₂O₅. Based on the simulated structural information, a modified random network model for P₂O₅-bearing sodium aluminosilicate glass has been proposed, which could be useful to optimize the mobility of sodium ions and design novel functional glass compositions.

(A) A modified structural model proposed for P₂O₅-bearing sodium aluminosilicate glasses. (B) Degree of preferred connection (DPC) of different T–O–T network linkage for LAP, MAP and HAP glass compositions with various P₂O₅ content.

Structure-Elasticity Relationship of Potassium Silicate Glasses from Brillouin Light Scattering Spectroscopy and Molecular Dynamics Simulations

Brillouin light scattering (BLS) spectroscopy and molecular dynamic (MD) simulations allowed the identification of a relationship between the elastic properties and the structure of K-containing glasses of formula (K₂O)_x-(SiO₂)_1-x, having different K₂O concentrations. Excellent agreement was observed between experimental data and simulations. The peculiar elastic properties observed for these potassium silicate glasses have been extensively discussed in terms of structural and energetic features of the materials. Elastic properties were shown to be strongly dependent on the asymmetry of potential energy in the K-BO interactions and the K-NBO interactions. A low K₂O content (below 10-15% K₂O) appeared to be in favor of K⁺-BO interactions and high asymmetry of potential energy, whereas a high K₂O content (from 10 to 15% K₂O) was in favor of K⁺-NBO interactions with lower asymmetry. Our results suggest a possible explanation to the observed anomalous dependence of elastic properties of potassium silicate glasses with K₂O amount.

Topological Origins of the Mixed Alkali Effect in Glass

The mixed alkali effect, the deviation from expected linear property changes when alkali ions are mixed in a glass, remains a point of contention in the glass community. While several earlier models have been proposed to explain mixed alkali effects on ionic motion, models based on or containing discussion of structural aspects of mixed-alkali glasses remain rare by comparison. However, the transition-range viscosity depression effect is many orders in magnitude for mixed-alkali glasses, and the original observation of the effect (then known as the Thermometer Effect) concerned the highly anomalous temperature dependence of stress and structural relaxation time constants. With this in mind, a new structural model based on topological constraint theory is proposed herein which elucidates the origin of the mixed alkali effect as a consequence of network strain due to differing cation radii. Discussion of literature models and data alongside new molecular dynamics simulations and experimental data are presented in support of the model, with good agreement.

The role of different network modifying cations on the speciation of the Co2+ complex in silicates and implication in the investigation of historical glasses

In the last decades the speciation of the cobalt complex in a glass matrix has been extensively studied. Bivalent cobalt ions in glasses of different composition commonly adopt a tetrahedral coordination, though hexa- or penta-coordinated species are also possible. Changes in the absorbance spectrum of Co-doped glasses were attested in previous studies according to the introduction of different modifying cations. A shifting of the first sub-band characterizing the typical triplets of tetrahedral Co²⁺ ions in both the visible and near infrared regions was observed, but discrepancies in literature suggested a relevant role of glass composition on the definition of the optical signature of cobalt. Co-doped glasses with different composition (soda-lime, potash-lime, mixed alkali and ZnO-Na₂O-CaO-SiO₂) were studied via Fiber Optic Reflectance Spectroscopy (FORS). Pseudo-Voigt functions were used for the deconvolution of the absorbance spectra and the features of the bands characteristic of each cobalt complex were investigated. The structural role played by each modifying cation and the fundamental implications of glass basicity on the speciation of different Co-complexes were stressed. Changes in glass structure resulted in different equilibria between the three absorbing species whose specific optical signatures in the 480-530nm region interact to determine the resulting absorbance spectrum.

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.

Probing the Nanodomain Origin and Phase Transition Mechanisms in (Un)Poled PMN-PT Single Crystals and Textured Ceramics

Outstanding electrical properties of solids are often due to the composition heterogeneity and/or the competition between two or more sublattices. This is true for superionic and superprotonic conductors and supraconductors, as well as for many ferroelectric materials. As in PLZT ferroelectric materials, the exceptional ferro- and piezoelectric properties of the PMN-PT ((1−x)PbMg_1/3Nb_2/3O₃−xPbTiO₃) solid solutions arise from the coexistence of different symmetries with long and short scales in the morphotropic phase boundary (MPB) region. This complex physical behavior requires the use of experimental techniques able to probe the local structure at the nanoregion scale. Since both Raman signature and thermal expansion behavior depend on the chemical bond anharmonicity, these techniques are very efficient to detect and then to analyze the subtitle structural modifications with an efficiency comparable to neutron scattering. Using the example of poled (field cooling or room temperature) and unpoled PMN-PT single crystal and textured ceramic, we show how the competition between the different sublattices with competing degrees of freedom, namely the Pb-Pb dominated by the Coulombian interactions and those built of covalent bonded entities (NbO₆ and TiO₆), determine the short range arrangement and the outstanding ferro- and piezoelectric properties.

Structural models of bioactive glasses from molecular dynamics simulations

The bioactive mechanism, by which living tissues attach to and integrate with an artificial implant through stable chemical bonds, is at the core of many current medical applications of biomaterials, as well as of novel promising applications in tissue engineering. Having been employed in these applications for almost 40 years, soda-lime phosphosilicate glasses such as 45S5 represent today the paradigm of bioactive materials. Despite their strategical importance in the field, the relationship between the structure and the activity of a glass composition in a biological environment has not been studied in detail. This fundamental gap negatively affects further progress, for instance, to improve the chemical durability and tailor the biodegradability of these materials for specific applications. This paper reviews recent advances in computer modelling of bioactive glasses based on molecular dynamics simulations, which are starting to unveil key structural features of these materials, thus contributing to improve our fundamental understanding of how bioactive materials work.

Correlation between floppy to rigid transitions and non-Arrhenius conductivity in glasses

Non-Arrhenius behavior is observed for a number of virgin potassium silicate glasses (1-x)SiO2-xK2O with a potassium oxide concentration larger than a certain value x = x(c). Recovering of Arrhenius behavior is provided by the annealing that enhances densification. These various results are the manifestation of the floppy or rigid nature of the network. Compositional effects and saturation with temperature can be analyzed with a combination of constraint theory and a point defect model. They underscore the key role played by network rigidity for the understanding of conduction and saturation effects in glassy electrolytes.

Structural Change in Li and Na Aluminophosphate Glasses: Evidence of a “Structural Mixed Alkali Effect”

The short- and long-range structure of a series of single and mixed aluminophosphate glasses with the general composition [xNa(2)O (46 - x)Li(2)O], [yAl(2)O(3) (54 - y)P(2)O(5)] is analyzed using (31)P and (27)Al magic-angle spinning (MAS) NMR as well as small-angle X-ray scattering. These series of glasses allow analyzing both the effect of alumina incorporation in these glasses, for small alumina content (y = 0, 4, 8), and the structural changes associated with the so-called mixed alkali effect (x = 0, 11.5, 23, 34.5, 46). Our results indicate that aluminum is mainly octahedrally coordinated in these glasses and that there is most likely some segregation of the Al(OP)(6) species. In the pure phosphate glasses, we observe a "classical" continuous variation of the structural properties with the relative alkali content, but in the aluminophosphate, both local and long-range structural results reveal for the first time some nonlinear change as a function of the relative alkali content.

Molecular Dynamics Simulations of the Effect of the Composition of Calcium Alumino-Silicate Intergranular Films on Alumina Grain Growth

Molecular dynamics simulations were performed to study the effect of the composition of the intergranular film (IGF) on anisotropic and isotropic grain growth in alpha-Al2O3. In the simulations, the IGF is formed while in contact with two differently oriented alumina crystals, with the alumina (0001) basal plane on one side and the (110) prism plane on the other. Five different compositions in the IGFs were studied. Results show preferential growth along the [110] of the (110) surface in comparison to growth along the [0001] direction on the (0001) surface for compositions near a Ca/Al ratio of 0.5. Such preferential growth is consistent with anisotropic grain growth in alumina, where platelets form because of faster growth of the prism orientations than the basal orientation. The simulations also show the mechanism by which Ca ions in the IGF inhibit growth on the basal surface. At compositions with high or low Ca/Al ratios, growth along each surface normal is equivalent, indicating isotropic grain growth, although the attachment rates are quite different, which may indicate differences between normal grain growth and abnormal, but isotropic, grain growth. The simulations provide an atomistic view of attachment onto crystal surfaces, affecting grain growth in alumina.

Molecular dynamics study of cage decay, near constant loss, and crossover to cooperative ion hopping in lithium metasilicate

Molecular dynamics (MD) simulations of lithium metasilicate (Li2SiO3) in the glassy and supercooled liquid states have been performed to illustrate the decay with time of the cages that confine individual Li+ ions before they hop out to diffuse cooperatively with each other. The self-part of the van Hove function of Li+ ions, G(s)(r,t), is used as an indicator of the cage decay. At 700 K, in the early time regime t<t(x1), when the cage decays very slowly, the mean square displacement <r(2)> of Li+ ions also increases very slowly with time approximately as t(0.1) and has weak temperature dependence. Such <r(2)> can be identified with the near constant loss (NCL) observed in the dielectric response of ionic conductors. At longer times, when the cage decays more rapidly as indicated by the increasing buildup of the intensity of G(s)(r,t) at the distance between Li+ ion sites, <r(2)> broadly crosses over from the NCL regime to another power law t(beta) with beta approximately 0.64 and eventually it becomes t(1.0), corresponding to long-range diffusion. Both t(beta) and t(1.0) terms have strong temperature dependence and they are the analogs of the ac conductivity [sigma(omega) proportional, variant omega(1-beta)] and dc conductivity of hopping ions. The MD results in conjunction with the coupling model support the following proposed interpretation for conductivity relaxation of ionic conductors: (1) the NCL originates from very slow initial decay of the cage with time caused by few independent hops of the ions because t(x1)<<tau(o), where tau(o) is the independent hop relaxation time; (2) the broad crossover from the NCL to the cooperative ion hopping conductivity sigma(omega) proportional, variant omega(1-beta) occurs when the cage decays more rapidly starting at t(x1); (3) sigma(omega) proportional, variant omega(1-beta) is fully established at a time t(x2) comparable to tau(o) when the cage has decayed to such an extent that thereafter all ions participate in the slowed dynamics of cooperative jump motion; and (4) finally, at long times sigma(omega) becomes frequency independent, i.e., the dc conductivity. MD simulations show the non-Gaussian parameter peaks at approximately t(x2) and the motion of the Li+ ions is dynamically heterogeneous. Roughly divided into two categories of slow (A) and fast (B) moving ions, their mean square displacements <r(2)(A)> and <r(2)(B)> are about the same for t<t(x2), but <r(2)(B)> of the fast ions increases much more rapidly for t>t(x2). The self-part of the van Hove function of Li+ reveals that first jumps for some Li+ ions, which are apparently independent free jumps, have taken place before t(x2). While after t(x2) the angle between the first jump and the next is affected by the other ions, again indicating cooperative jump motion. The dynamic properties are analogous to those found in supercooled colloidal particle suspension by confocal microscopy.

Deconvolution of 29Si magic-angle spinning nuclear magnetic resonance spectra of silicate glasses revisited--some critical comments

It is common practice to quantify 29Si magic-angle spinning nuclear magnetic resonance (MAS NMR) spectra of silicate glasses by using lineshape fitting routines based on Gaussian distribution functions. This procedure depends on several -implicit -assumptions. In this contribution we have addressed and questioned these underlying assumptions. Also "goodness-of-fit" criteria and experimental considerations ("goodness-of-data") have been discussed. Some illustrative practical examples (silica glass, binary Na2O-SiO2 glasses with different Na:Si ratio) are given.

Quantification of the disorder in network-modified silicate glasses

Local order in silicate glasses has been observed by many experimental techniques to be similar to that in crystalline materials. Details of the intermediate-range order are more elusive, but essential for understanding the lack of long-range symmetry in glasses and the effect of composition on glass structure. Two-dimensional 17O dynamic-angle-spinning nuclear magnetic resonance experiments reveal intermediate-range order in the distribution of inter-tetrahedral (Si-O-Si) bond angles and a high degree of order in the disposition of oxygen atoms around the network-modifying cations.