Temperature-induced structural change through the glass transition of silicate glass by neutron diffraction
Not Available
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Infrared (IR) and Raman spectroscopic features of silicate glasses are usually interpreted based on the analogy with those of smaller molecules, molecular clusters, or crystalline counterparts; this study tests the accuracy and validity of these widely cited peak assignment schemes by comparing vibrational spectral features with bond parameters of the glass network created by molecular dynamics (MD) simulations. A series of sodium silicate glasses with compositions of [Na 2 O] x [Al 2 O 3 ] 2 [SiO 2 ] 98- x with x = 7, 12, 17, and 22 were synthesized and analyzed with IR and Raman. A silica glass substrate and a crystalline quartz were also analyzed for comparison. Glass structures with the same compositions were generated with MD simulations using three types of potentials: fixed partial charge pairwise (Teter), partial diffuse charge potential (MGFF), and bond order-based charge transfer potential (ReaxFF). The comparison of simulated and experimental IR spectra showed that, among these three potentials tested, ReaxFF reproduces the concentration dependence of spectral features closest to the experimentally observed trend. Hence, the bond length and angle distributions as well as Si– Q n species and ring size distributions of silica and sodium silicate glasses were obtained from ReaxFF-MD simulations and further compared with the peak assignment or deconvolution schemes—which have been widely used since 1970s and 1980s—(a) correlation between the IR peak position in the Si–O stretch region (1050-1120 cm -1 ) and the Si–O–Si bond angle; (b) deconvolution of the Raman bands in the Si–O stretch region with the Q n speciation; and (c) assignment of the Raman bands in the 420-600 cm -1 region to the bending modes of (SiO) n rings with different sizes (typically, n = 3-6). The comparisons showed that none of these widely used methods is congruent with the bond parameters or structures of silicate glass networks produced via ReaxFF-MD simulations. This result invokes that the adequacy of these spectral interpretation methods must be questioned. Alternative interpretations are proposed, which are to be tested independently in future studies.
Reactive classical molecular dynamics simulations of sodium silicate glasses, xNa 2 O–(100 – x)SiO 2 (x = 10–30), under quasi-static loading, were performed for the analysis of molecular scale fracture mechanisms. Mechanical properties of the sodium silicate glasses were consistent with experimentally reported values, and the amount of crack propagation varied with reported fracture toughness values. The most crack propagation occurred in NS20 systems (20-mol% Na 2 O) compared with the other simulated compositions. Dissipation via two mechanisms, the first through sodium migration as a lower activation energy process and the second through structural rearrangement as a higher activation energy process, was calculated and accounted for the energy that was not stored elastically or associated with the formation of new fracture surfaces. Here, a correlation between crack propagation and energy dissipation was identified, with systems with higher crack propagation exhibiting less energy dissipation. Sodium silicate glass compositions with lower energy dissipation also exhibited the most sodium movement and structural rearrangement within 10 Å of the crack tip during loading. Therefore, high sodium mobility near the crack tip may enable energy dissipation without requiring formation of structural defects. Therefore, the varying mobilities of the network modifiers near crack tips influence the brittleness and the crack growth rate of modified amorphous oxide systems.
Prediction of brittle fracture in silicate glasses requires a fundamental understanding of slow crack growth rates, far below the critical stress intensity (K IC ) of the material. Here, in this study, atomic force microscopy (AFM) was used to characterize slow crack propagation in soda lime silicate glasses. Using large-sample AFM systems with specialized tips and scanners, fractures less than 20 nm wide were measured as a function of relative humidity (RH). The results were found to be sensitive to the type of AFM tip used and to the dryness of the surface. Continuous scanning and image capture over several days allowed the determination of the crack tip velocity. The dependance of crack tip velocity on RH was measured by enclosing the sample and AFM head in a humidity-controlled chamber. The crack tip velocity was measured as low 4.0×10 -12 m/s at 15 % RH and as high as 3.7×10 -10 m/s at 40 % RH. Overall crack growth was smooth during the AFM scan. Additionally, recession of the crack tip was noted as the sample dried, indicating crack tip closure may be occurring. Further studies on alternative compositions will provide details on subcritical fracture processes in a range of different environmental conditions.
LionGlass™ is Penn State University's patent-pending alumino-silicophosphate glass compositional family that offers, for the first time, a viable low-melting, carbonate-free alternative to soda lime silicate glass by (a) lowering the melting temperature by up to 400°C compared to soda lime silicate glass and (b) eliminating the use of carbonate batch materials. The viscosity curve of LionGlass is shifted downward in temperature compared to soda lime silicate glass, enabling melting at much lower temperatures. LionGlass also has significantly improved resistance to crack initiation compared to soda lime silicate glass. Other properties, such as thermal expansion coefficient and optical dispersion curve, are comparable to those of traditional soda lime silicate glass.
Recent development of reactive force fields have enabled molecular dynamics simulations of interactions between silicate glasses and water at the atomistic scale. While multicomponent silicate glasses encompass a wide variety of compositions and properties, one common structural feature in these glasses is the combination of the network structure that is made up of silica tetrahedra linked through corner sharing interspersed with network modifiers like alkali and alkaline-earth ions that break up the Si–O–Si linkages by forming nonbridging oxygen. In reactions with water, ion exchange between alkali ions in the glass and proton or hydronium in the solution, as well as hydrolysis reaction of the Si–O–Si linkages and subsequent silanol formation, is observed and well documented. Furthermore, we have used a set of recently developed reactive force field to investigate the reactions between water and the surfaces of silica and sodium silicate glasses of different compositions for reactions up to 8 nanoseconds. Our results indicate sodium leaching into water and diffusion of water molecules up to 25 Å into the glass surface. We examined the structural and compositional changes inside the glass and around the diffused ions and use these to explain the rates of silanol formation at the surface. We also observed proton transport in the glass which has an indirect influence on the silanol formation rates. While the surface of the glass was rough to start with, it undergoes further modification into a hydrated gel-like structure in the glass for up to 5 Å in the higher alkali containing glasses. It was found that the leached sodium ions remain close to the interface and that fragments of silicate network from the surface is capable of dislodging from the bulk glass and enter the aqueous solution. These simulations thus provide insights into the formation and structure of an alteration layers commonly observed in multicomponent silicate glasses corroded in aqueous solutions.
Decarbonizing the glass industry requires alternative melting technology, as current industrial melting practices rely heavily on fossil fuels. Hydrogen has been proposed as an alternative to carbon-based fuels, but the ensuing consequences on the mechanical behavior of the glass remain to be clarified. A critical distinction between hydrogen and carbon-based fuels is the increased generation of water during combustion, which raises the equilibrium solubility of water in the melt and alters the behavior of the resulting glass. A series of five silicate glasses with 80% silica and variable [Na 2 O]/([H 2 O] + [Na 2 O]) ratios were simulated using molecular dynamics to elucidate the effects of water on fracture. Several fracture toughness calculation methods were used in combination with atomistic fracture simulations to examine the effects of hydroxyl content on fracture behavior. Here, this study reveals that the crack propagation pathway is a key metric to understanding fracture toughness. Notably, the fracture propagation path favors hydrogen sites over sodium sites, offering a possible explanation of the experimentally observed effects of water on fracture properties.
The role of Al in aluminosilicate glasses remains somewhat a mystery: at low concentrations, it increases the resistance to hydrolysis of the glass, whereas at high concentrations an opposite effect is observed. To understand the origin of the phenomenon on a fundamental atomistic scale, we performed 577 MD simulations and applied potential mean force (PMF) calculations to estimate the activation barriers for hydrolysis and to statistically correlate them with local structural features of the glass. Models of pure silicate and aluminosilicate glasses are constructed and investigated. PMF simulation results are further validated by the experimental measurements and revealed that Al is very easy to dissociate, but it also increases the glass chemical durability through significantly increasing both the strength of Si and network connectivity of the glass. In contrast, at high Al concentration, preferential dissolution of Al weakens the silicate network, which it supposes to strengthen, and so the glass resistance becomes poor. Through PMF calculations, we evaluated the activation barriers for dissociating bonds around Al as 0.49 eV, which is less than a half of the energy to dissociate bonds around Si in pure silicate (1.22 eV) and around Si in aluminosilicate glass (1.34 eV), all these energy differences being statistically significant. Molecular structural level investigation revealed that Si with Al as a second neighbor in the glass network has a significantly higher activation energy for dissociation than Si in pure silicate glass. Furthermore, the proposed approach opens the way to the development of quantitative predictive models of glass durability.
Ring size distribution of silicate glass can be determined from the first sharp diffraction peak of neutron structure factor.
Glass alteration in the presence of microorganisms has been a topic of research for over 150 years. Researchers from a variety of disciplines, including material science, biology, chemistry, geology, physics, and cultural heritage materials preservation have conducted experiments in this area to try and understand when, how, and why microorganism may interact and subsequently influence the alteration of glass. The breadth and depth of these studies are the topic of this review. This review presents a detailed history and a comprehensive overview of this field of research, while maintaining focus on the terrestrial alteration of anthropogenic silicate glasses. Within this manuscript is a schema for bio-interaction with silicate glasses and an outline of an evidence-based hypothesis on how these interactions may influence glass alteration processes. Topics discussed include microbial colonization of glass, development, and interactions of biofilms with glass surface, abiotic vs. biotic alteration processes, and signatures of bio-alteration. Future research needs and a discussion of practical drivers for this research are summarized.
The chemical and structural stability of two commercial multicomponent silicate glasses (SCN and G6) in contact with yttria-stabilized zirconia (YSZ) was investigated after exposure times of up to 40,000 h in air at 800 °C. With exposure time, interfacial layers develop at the SCN-YSZ and G6-YSZ interfaces, which were characterized in detail using both quantitative chemical analysis and atomic-resolution imaging. At the SCN-YSZ interface, a Ca-Ba-Si-O reaction phase was found to grow by diffusion control. In G6-YSZ, Raman spectroscopy and electron microscopy revealed a disorganized interfacial reaction later between G6 and YSZ, and the occurrence of cubic to tetragonal to monoclinic phase transformations in YSZ. Finally, this microstructural evolution is discussed in terms of devitrification resistance of glass and diffusion processes at interfaces.
Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.
Vibrational spectroscopy has been widely used to investigate various structural aspects of the glass network, and there are a plethora of papers reporting subtle but consistent changes in infrared and Raman spectral features of glass upon alterations of glass compositions, thermal histories, mechanical stresses, or surface treatments. However, interpretations of such spectral features are still obscured due to the lack of well-established physical principles accurately describing vibrational modes of the non-crystalline glass network. Due to the non-equilibrium nature of the glass network, three-dimensionally connected without any long-range orders, vibrational spectral features of glass cannot be interpreted using the analogy with those of isolated molecular moieties or crystalline counterparts. Here, this feature article explains why such comparisons are outdated and describes the recent advances made from theoretical calculations of vibrational spectral features of amorphous networks or comparisons of computational results with experimental data. For the interpretation of vibrational spectral features of silica and silicate glasses, the following empirical relationships are suggested: (i) the intensity-weighted peak position of the Si-O-Si stretch mode negatively correlates with the weighted average of the Si-O bond length distribution, and (ii) the broad band of the Si-O-Si bending mode negatively correlates with the Si-O-Si bond angle distribution. Selected examples of vibrational spectroscopic imaging of surface defects are discussed to deliberate the implication of these findings in the structure-property relationship of silica and silicate glass materials. Unanswered questions and continuing research challenges are identified.
Wavelength-dependent mid-infrared (400–1,800 cm -1 ; 5.56–25 μm) optical constants (real and imaginary indices of refraction; n and k) are determined using reflectance spectra at a spectral sampling of 4 cm -1 for several silicate glasses of varying SiO 2 wt% which include (a) basaltic volcanic glass from Kilauea, Hawaii, (b) synthetic andesite, (c) two synthetic dacites, (d) obsidian volcanic glass from Mount Lassen, California, (e) synthetic rhyodacite, and (f) rhyolitic volcanic glass from Mexico. Because glasses are optically isotropic, no specific orientation was required for spectral measurements, and polished glass samples were measured at random orientations using micro-FTIR spectrometer. Lorentz-Lorenz dispersion theory and Fresnel reflectance model for high symmetry materials were used to model the optical constants by optimizing oscillator parameters in modeled spectra to match laboratory spectra. A genetic algorithm (GA) approach automatically finds the natural oscillators and their parameters for each spectrum, and then uses these parameters in a non-linear least squares optimization routine. The study compared spectral parameters such as Christiansen feature, reststrahlen bands (RBs), and the peaks centered around 860–1,100 (peak 1) and 400–480 cm -1 (peak 2) of n and k to their respective SiO 2 wt% of the glasses. CF, RBs, peak 1 of n and k, and peak 2 of n shift to higher wavenumbers with increased SiO 2 wt%, whereas peak 2 of k shifts to lower wavenumbers with increased SiO 2 wt%. Derived optical constants of these glasses will improve quantitative abundance mapping of volcanic materials on the surfaces of silicate targets in the Solar System.
Incorporation of small molecules in silicate melts may provide an important mechanism for storing noble gases in the deep Earth, yet the means by which chemically inert noble gases enter and are retained in silica-based materials is not understood. High pressure, room temperature sound velocity measurements on silica and natural basalt glasses in different pressure-transmitting media reveal that neon enters the structure of silicate glasses and enhances their elastic strengths, whereas an ethanol-methanol mixture does not. Combined with literature data, we found the incorporation of small molecules into silica and basalt glasses is controlled by the void size distribution of the glass and size of the molecules. Pressure primarily reduces the size of noble gases, thereby increasing their solubilities in silicate melts and glasses.
Despite decades of studies, the nature of the glass transition remains elusive. In particular, the sharpness of the dynamical arrest of a melt at the glass transition is captured by its fragility. Here, we reveal that fragility is governed by the medium-range order structure. Based on neutron-diffraction data for a series of aluminosilicate glasses, we propose a measurable structural parameter that features a strong inverse correlation with fragility, namely, the average medium-range distance (MRD). We use in-situ high-temperature neutron-scattering data to discuss the physical origin of this correlation. We argue that glasses exhibiting low MRD values present an excess of small network rings. Such rings are unstable and deform more readily with changes in temperature, which tends to increase fragility. These results reveal that the sharpness of the dynamical arrest experienced by a silicate glass at the glass transition is surprisingly encoded into the stability of rings in its network.
Photothermal atomic force microscopy coupled with infrared spectroscopy (AFM-IR) brings significant value as a spatially resolved surface analysis technique for disordered oxide materials such as glasses, but additional development and fundamental understanding of governing principles is needed to interpret AFM-IR spectra, since the existing theory described for organic materials does not work for materials with high extinction coefficients for infrared (IR) absorption. This paper describes theoretical calculation of a transient temperature profile inside the IR-absorbing material considering IR refraction at the interface as well as IR adsorption and heat transfer inside the sample. Furthermore, this calculation explains the differences in peak positions and amplitudes of AFM-IR spectra from those of specular reflectance and extinction coefficient spectra. It also addresses the information depth of the AFM-IR characterization of bulk materials. AFM-IR applied to silica and silicate glass surfaces has demonstrated novel capability of characterizing subsurface structural changes and surface heterogeneity due to mechanical stresses from physical contacts, as well as chemical alterations manifested in surface layers through aqueous corrosion.
Past work has shown that water or acid soaking treatments can increase the mechanical strength of soda lime silicate (SLS) glasses. In this work, we show that acid leaching treatments result in an increase in the apparent crack resistance of the acid-leached surface of SLS glass during indentation. Vickers indentation tests in controlled environments show a humidity dependence of radial cracking, suggesting that the transport of water through the leached layer plays a critical role in the propagation of cracks to the glass surface. Here, molecular dynamics simulations with reactive force fields indicate that the leached surface layer can undergo pressure-induced mechanochemical reactions during indentation, which increases the bridging oxygen connectivity in the silica network of the leached layer. Such structural changes can hinder transport of water molecules from the environment to the subsurface crack tip. Based on experimental observations and simulation results, a new hypothesis is proposed that mechanochemical restructuring in the leached layer in response to the applied load may lower the transport kinetics of molecular water to critical flaws, resulting in an apparent enhancement in the crack resistance of the acid-leached surface of SLS glass.