skip to main content


Title: Assessment of interatomic parameters for the reproduction of borosilicate glass structures via DFT‐GIPAW calculations
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 Na2OxB2O3(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
NSF-PAR ID:
10447469
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Journal of the American Ceramic Society
Volume:
102
Issue:
12
ISSN:
0002-7820
Page Range / eLocation ID:
p. 7225-7243
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The structures of glasses in the lithium–bismuth orthoborate composition range deviate significantly from the short‐range order structure of the two crystalline end‐members. Although binary Li3BO3and BiBO3are solely of comprised trigonal orthoborate anions, all glasses formed by their combination contain four‐coordinated borate tetrahedra. We analyze the structure of (75−1.5x)Li2O–xBi2O3–(25+0.5x)B2O3glasses in increments ofx = 5, with11B magic‐angle spinning nuclear magnetic resonance (NMR), infrared (IR), and Raman spectroscopy. For the full series, the oxygen‐to‐boron ratio remains constant at O/B = 3:1. NMR quantifies an increase in the fraction of tetrahedral boron with increasing bismuth oxide content. Evolution of the mid‐IR profile suggests multiple types of tetrahedral boron sites. Raman spectroscopy reveals that Bi2O3tends to cluster within the lithium borate matrix when initially introduced and that this behavior transforms into a bismuthate network with increasing bismuth oxide content. In all cases, mixed Bi–O–B linkages are observed. The dual role of bismuth as network modifier and network former is likewise observed in the far IR. The glass transition temperature continuously increases with bismuth oxide content; however, the glass stability displays a maximum in the multicomponent glass ofx = 40.

     
    more » « less
  2. Abstract

    The modifier field strength (FS) is believed to play an important role in determining the elastic–plastic responses of aluminoborosilicate (ABS) glasses, but its effect is not well understood. Three novel alkali and three alkaline earth (AE) ABS compositions were created for this study which is the first part of two studies that explored the elastoplastic responses of these glasses. Six compositions were designed using different network modifiers (NWMs) to cover a range of cation FS. The glasses were also designed such that the concentrations of NWM and Al2O3were similar, which maximized the three‐coordinated boron fraction in the network. It is well known that modifier FS can affect the coordination number (CN) of Al and B in an ABS glass structure, for example, a higher FS modifier can promote B3 → B4and higher [Al5,6], but the degree of this depends on network former (NWF) ratios. Previous work used solid‐state NMR spectroscopic analysis on the current glasses to find that there was variation between [B4] and [Al4] between the two glass series (alkali vs. AE) but that was attributed to synthesis factors and no trend with FS was associated with the varying NWF CN. Further,29Si ssNMR showed no evidence of NBOs which made sense based on composition. The conclusion, therefore, was that there was a far greater correlation with modifier FS for the increased mechanical and physical properties rather than the CN of Al and B. Part I of the current work focused on the elastic moduli, Poisson's ratio, the indentation size effect (ISE), and the bow‐in parameter. This part laid out the foundation to investigate the intersection of these elastoplastic properties with hardness and crack resistance as a function of NWM FS. Results showed that: (i) the Young's, bulk, and shear moduli increased with modifier FS, whereas Poisson's ratio did not trend with FS; (ii) the alkali glasses had a significantly higher magnitudes of ISE compared to the AE glasses; and (iii) the bow‐in parameter was load dependent and decreased with modifier FS at the highest indentation load.

     
    more » « less
  3. Abstract

    Iron oxide is commonly found in natural or industrial glass compositions and can exist as ferrous (Fe2+) or ferric (Fe3+) species, with their ratios depending on glass composition, temperature, pressure and the redox reactions during the glass forming process. The iron redox ratio plays an important role on silicate glass structures and consequently various properties. This work aims to study the effect of iron oxide, and particularly the iron redox ratio, on the structures of borosilicate and boroaluminosilicate glasses using molecular dynamics simulations with newly developed iron potential parameters that are compatible with the borosilicate potentials. The results provide detailed cation coordination states of both iron species and the effect of redox ratio on boron coordination and other structural features. Particularly, competition for charge compensating modifier cations (such as Na+) among the fourfold‐coordinated cations such as B3+, Al3+, and Fe3+is investigated by calculating the cation–cation pair distribution functions and coordination preferential ratios. The results show that the trivalent ferric ions, with a shorter Fe–O bond distance and better defined first coordiation shell with mainly four‐fold coordination, act as a glass former whereas the divalent ferrous ions mainly play the role of glass modifier. The ferrous/ferric ratio (Fe2+/Fe3+) was found to affect the glass chemistry and hence glass properties by regulating the amount of four‐coordinated boron, the fraction of non‐briding oxygen and other features. The results are compared with available experimental data to gain insights of the complex structures and charge compensation schemes of the glass system.

     
    more » « less
  4. Abstract

    Glass for pharmaceutical packaging requires high chemical durability for the safe storage and distribution of newly developed medicines. In borosilicate pharmaceutical glasses which typically contain a mixture of different modifier ions (alkali or alkaline earth), the dependence of the chemical durability on alkaline earth oxide concentrations is not well understood. Here, we have designed a series of borosilicate glasses with systematic substitutions of CaO with MgO while keeping their total concentrations at 13 mol% and a fixed Na2O concentration of 12.7 mol%. We used these glasses to investigate the influence ofR = [MgO]/([MgO] + [CaO]) on the resistance to aqueous corrosion at 80°C for 40 days. It was found that this type of borosilicate glass undergoes both leaching of modifier ions through an ion exchange process and etching of the glass network, leading to dissolution of the glass surface. Based on the concentration analysis of the Si and B species dissolved into the solution phase, the dissolved layer thickness was found to increase from ~100 to ~170 nm asRincreases from 0 to 1. The depth profiling analysis of the glasses retrieved from the solution showed that the concentration of modifier ions (Na+, Ca2+, and Mg2+) at the interface between the solution and the corroded glass surface decreased to around 40%–60% of the corresponding bulk concentrations, regardless ofRand the leaching of modifier cations resulted in a silica‐rich layer in the surface. The leaching of Ca2+and Mg2+ions occurred within ~50 and <25 nm, respectively, from the glass surface and this thickness was not a strong function ofR. The leaching of Na+ions varied monotonically; the thickness of the Na+depletion layer increased from ~100 nm atR = 0 to ~200 nm atR = 1. Vibrational spectroscopy analysis suggested that the partial depletion of the ions may have caused some degree of the network re‐arrangement or re‐polymerization in the corroded layer. Overall, these results suggested that for the borosilicate glass, replacing [CaO] with [MgO] deteriorates the chemical durability in aqueous solution.

     
    more » « less
  5. Abstract

    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]x[Al2O3]2[SiO2]98−xwithx = 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. Thus, the bond length and angle distributions as well as Si–Qnspecies 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 theQnspeciation; and (c) assignment of the Raman bands in the 420‐600 cm−1region to the bending modes of (SiO)nrings 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 finding 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.

     
    more » « less