Understanding surface reactions of silicate glass under interfacial shear is critical as it can provide physical insights needed for rational design of more durable glasses. Here, we performed reactive molecular dynamics (MD) simulations with ReaxFF potentials to study the mechanochemical wear of sodium silicate glass rubbed with amorphous silica in the absence and presence of interfacial water molecules. The effect of water molecules on the shear‐induced chemical reaction at the sliding interface was investigated. The dependence of wear on the number of interfacial water molecules in ReaxFF‐MD simulations was in reasonable agreement with the experimental data. Confirming this, the ReaxFF‐MD simulation was used to find further details of atomistic reaction dynamics that cannot be obtained from experimental investigations only. The simulation showed that the severe wear in the dry condition is due to the formation of interfacial Sisubstrate–O–Sicounter_surfacebond that convey the interfacial shear stress to the subsurface and the presence of interfacial water reduces the interfacial bridging bond formation. The leachable sodium ions facilitate surface reactions with water‐producing hydroxyl groups and their key role in the hydrolysis reaction is discussed.
Infrared (IR) and Raman spectroscopic features of silicate glasses are often 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 [Na2O]
- NSF-PAR ID:
- 10373556
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Journal of the American Ceramic Society
- Volume:
- 103
- Issue:
- 6
- ISSN:
- 0002-7820
- Page Range / eLocation ID:
- p. 3575-3589
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract 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. 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.
-
more » « less
Abstract Borates and borosilicates are potential candidates for the design and development of glass formulations with important industrial and technological applications. A major challenge that retards the pace of development of borate/borosilicate based glasses using predictive modeling is the lack of reliable computational models to predict the structure‐property relationships in these glasses over a wide compositional space. A major hindrance in this pursuit has been the complexity of boron‐oxygen bonding due to which it has been difficult to develop adequate B–O interatomic potentials. In this article, we have evaluated the performance of three B–O interatomic potential models recently developed by Bauchy et al [
J .Non‐Cryst. Solids , 2018, 498, 294–304], Du et al [J. Am. Ceram. Soc. https://doi.org/10.1111/jace.16082 ] and Edèn et al [Phys. Chem. Chem. Phys. , 2018, 20, 8192–8209] aiming to reproduce the short‐to‐medium range structures of sodium borosilicate glasses in the system 25 Na2Ox B2O3(75 −x ) SiO2(x = 0‐75 mol%). To evaluate the different force fields, we have computed at the density functional theory level the NMR parameters of11B,23Na, and29Si of the models generated with the three potentials and the simulated MAS NMR spectra compared with the experimental counterparts. It was observed that the rigid ionic models proposed by Bauchy and Du can both reliably reproduce the partitioning between BO3and BO4species of the investigated glasses, along with the local environment around sodium in the glass structure. However, they do not accurately reproduce the second coordination sphere of silicon ions and the Si–O–T (T = Si, B) and B‐O‐T distribution angles in the investigated compositional space which strongly affect the NMR parameters and final spectral shape. On the other hand, the core‐shell parameterization model proposed by Edén underestimates the fraction of BO4species of the glass with composition 25Na2O 18.4B2O356.6SiO2but can accurately reproduce the shape of the11B and29Si MAS‐NMR spectra of the glasses investigations due to the narrower B–O–T and Si‐O‐T bond angle distributions. Finally, the effect of the number of boron atoms (also distinguishing the BO3and BO4units) in the second coordination sphere of the network former cations on the NMR parameters have been evaluated. -
more » « less
Abstract Molecular dynamics simulations show that the expansion of silica glass occurs by the presence of the hydroxyl (SiOH) groups present in the glass as opposed to intact water (H2O) molecules, providing an accurate molecular description of the experimentally observed volume changes in silica glass exposed to water. Using a robust and accurate reactive potential, the simulations show that the expansion is caused by the rupture of siloxane (Si–O–Si) linkages in the glass via reactions with water molecules, forming SiOHs. Such reactions remove smaller rings and form larger rings, with a decrease in the overall number of rings smaller than a prescribed large ring size in comparison to dry glasses. This change in ring structure overcomes the inherently stronger hydrogen bonding in the glasses containing SiOH in comparison to the glasses containing predominantly intact H2O molecules. This stronger H‐bonding of the SiOH also causes a shift to lower frequencies in the high‐frequency OH vibrational spectrum for the silanols, as shown in previous ab‐initio calculations. This introduces a question about assuming the lower frequency part of the high‐frequency peak is only due to intact H2O molecules. A slight decrease in volume occurred in the glasses containing the largest concentration of intact H2O molecules. There is no change in the ring size distribution between the H2O glasses and dry glasses. Rather, the slight decrease in volume in the H2O system is caused by a decrease in siloxane bond angles caused by the formation of H‐bonds between the H2O molecules and the glass O in the siloxane cages surrounding the H2O molecules.
-
more » « less
Abstract Raman scattering is a powerful probe of
local structure (LS) of glasses. In Sodium Phosphate Glasses (SPGs), we show that bothLS composed of Qnspecies andExtended Range Structures (ERS) composed of Long Chains (LCs), Large Rings (LRs), and Small Rings (SRs) can be decoded by Raman scattering. The trimodal distribution of P‐Oterminalstretch modes of Q2species and P‐Obridgingatx < 50% are manifestations of theseERS . These two pairs of triads of modes are uniquely identified with Q2units present in either LCs, or LRs, or SRs. The existence three phases of c‐NaPO3composed of 3‐membered rings, 6‐membered rings, and infinitely long chains has facilitated the identification. The Intermediate Phase (IP) in SPGs extends in the 37.5 <x < 46.0% range, the Stressed‐rigid Phase in the 46.0% <x < 50%, and the Flexible Phase in the 18% <x < 37.5% range of soda. We show the IP consists predominantly of LCs (82%), with a minority of LRs (15%) and SRs (3%). The LR‐ and SR‐fractions increase measurably in the non‐IP phases. The structural finding is in harmony with the high configurational entropy of the IP glasses that leads aging to be qualitatively suppressed.
