9 Search Results

Predicting chemical durability of glass materials is important for various applications from daily life such as drink glass and kitchen ware to advanced technologies such as nuclear waste disposal and biomedicine. In this work, we explored prediction of initial dissolution rate through structural features from molecular dynamics (MD) simulations for a wide range of glass compositions (total 28) including borosilicate and aluminosilicate glasses, ZrO₂-containing and V₂O₅-containing boroaluminosilicate glasses. The initial dissolution rates (r₀) measured experimentally at 90 °C with varying solution conditions were correlated with structural features (e.g., polyhedral linkages and non-bridging oxygen species) obtained from MD simulations, either frommore » this study or from literature. Since hydrolysis of the glass network through breaking of the network former linkages (e.g., Si-O-Si, Si-O-Al, etc.) is a critical step of network glass dissolution, the statistics of these linkages obtained from MD were correlated to r0 through linear regression, where the coefficient of determination (R²) and root mean square error are found to be 0.949 and 0.681, respectively. This model was compared and discussed with existing models developed by various approaches including machine learning, the kinetic rate equation, topological constraint theory, and other descriptors from MD simulations. The discussion provides insights on future model improvements to predict glass dissolution. In addition, the impact of V₂O₅ on the glass dissolution was examined in detail, implicating that the impact is not the same across all glass compositions and test conditions.« less

Abstract Quantitative Structure Property Relationship (QSPR) analysis based on molecular dynamics (MD) simulations is a promising approach for establishing the composition‐property relationships of glasses and other materials with complex structures. A series of 20 borosilicate, aluminosilicate, and boroaluminosilicate glasses have been modeled using MD simulations with recently developed effective potentials. Short‐ and medium‐range structures of these glasses were analyzed and, based on the structural information, QSPR analysis of the initial dissolution rates ( r ₀ ) of these glasses that were measured at 90°C and pH 9 by using various structural descriptors such as percentage of bridging oxygen species, networkmore » connectivity, and average ring size. The structural descriptor, F _net , which contains both energetic information such as single bond strength and structural information such as cation coordination number and Q _n distribution, was also used. It was found that while the overall network connectivity, average ring size and F _net provide reasonable correlations with r ₀ of studied glasses, F _net gives the best correlation among the descriptors. For glasses that show incongruent dissolution, it was found that modification of glass compositions to account for preferential release of modifier cations is necessary to achieve best correlations. The findings were discussed with results of recent studies on evaluating the compositional dependence of glass dissolution behavior using the topological‐constraints‐based models.« less

Abstract High‐alumina containing high‐level waste (HLW) will be vitrified at the Waste Treatment Plant at the Hanford Site. The resulting glasses, high in alumina, will have distinct composition‐structure‐property (C‐S‐P) relationships compared to previously studied HLW glasses. These C‐S‐P relationships determine the processability and product durability of glasses and therefore must be understood. The main purpose of this study is to understand the detailed structural changes caused by Al:Si and (Al + Na):Si substitutions in a simplified nuclear waste model glass (ISG, international simple glass) by combining experimental structural characterizations and molecular dynamics (MD) simulations. The structures of these two series of glassesmore » were characterized by neutron total scattering and ²⁷ Al, ²³ Na, ²⁹ Si, and ¹¹ B solid‐state nuclear magnetic resonance (NMR) spectroscopy. Additionally, MD simulations were used to generate atomistic structural models of the borosilicate glasses and simulation results were validated by the experimental structural data. Short‐range (eg, bond distance, coordination number, etc) and medium‐range (eg, oxygen speciation, network connectivity, polyhedral linkages) structural features of the borosilicate glasses were systematically investigated as a function of the degree of substitution. The results show that bond distance and coordination number of the cation‐oxygen pairs are relatively insensitive to Al:Si and (Al + Na):Si substitutions with the exception of the B‐O pair. Additionally, the Al:Si substitution results in an increase in tri‐bridging oxygen species, whereas (Al + Na):Si substitution creates nonbridging oxygen species. Charge compensator preferences were found for Si‐[NBO] (Na ⁺ ), ^[3] B‐[NBO] (Na ⁺ ), ^[4] B (mostly Ca ²⁺ ), ^[4] Al (nearly equally split Na ⁺ and Ca ²⁺ ), and ^[6] Zr (mostly Ca ²⁺ ). The network former‐BO‐network former linkages preferences were also tabulated; Si‐O‐Al and Al‐O‐Al were preferred at the expense of lower Si‐O‐ ^[3] B and ^[3] B‐O‐ ^[3] B linkages. These results provide insights on the structural origins of property changes such as glass‐transition temperature caused by the substitutions, providing a basis for future improvements of theoretical and computer simulation models.« less

Despite their practical importance, atomistic modeling of B₂O₃ containing glasses have been challenging due to the lack of reliable empirical potentials. Fortunately, a few recent developments have shown promises to simulate these glasses where the boron coordination has complex and non-linear dependence on glass composition. This work aims to provide an evaluation of the effectiveness of three recently developed potentials by a systematic study of a series (~20) of sodium borosilicate glasses with constant K (ratio of [SiO₂]/[B₂O₃]) =2 and varying R (ratio of [Na₂O]/[B₂O₃] ranging from 0.1 to 4) values and several sodium boroaluminosilicate glasses. A comparative assessment wasmore » established on the basis of the short- and medium-range structure features, such as boron N₄ values, total correlation functions, bond angle distribution, oxygen speciation, and mechanical properties using experimental or well-established models as criteria. Furthermore, this work provides insights on the choices of empirical potentials for MD simulations of borosilicate glasses and gives directions of future potential development and refining.« less

Accurate prediction of boron coordination in multicomponent glasses is critical in glass science and technology as it strongly affects the properties of borate and borosilicate glasses. We have collected a dataset containing 657 glasses from literature with boron coordination values and developed models using analytical functions based on the well accepted Dell, Xiao and Bray model. Good prediction of boron coordination with a R² value higher than 0.8 was obtained. The large variation of boron coordination from experiments, originated from sample preparations and characterizations, led to difficulties in obtaining models with better prediction performance. Various machine learning (ML) algorithms weremore » evaluated and slightly better prediction performance was observed; however, interpretation of the ML models is less straight forward. In conclusion, this study developed various models capable of providing quantitative boron coordination predictions, providing insights into its structural roles in multi-component glasses, and suggesting fruitful areas for future research.« less

Atomistic computer simulations can offer insights into silicate glass-environment interactions with the recent development of reactive potentials. Yet, the accuracy of generated glass structures with these potential was usually not fully examined. In this paper, the capability of the reactive force field (ReaxFF) to describe the short and medium range structure features of sodium silicate glasses in molecular dynamics simulations is investigated by comparing a widely used partial charge pairwise potential and available experimental data. Glass structure information such as pair distribution function (PDF), coordination number, Qn species, neutron broadened structure factor, and X-ray broadened structure factor of the glassmore » structures from ReaxFF simulations were calculated and compared to evaluate the generated glass structure. Advantages and limitations of the potentials and glass forming procedures, as well as areas of further improvement, were discussed. The results show that the recently refined ReaxFF parameters through the proposed procedure enable the simulations of sodium silicate glass structures with minimal defects, which paves the way to investigate water-glass interaction mechanisms with the reactive enabled potentials.« less

Abstract Molecular dynamics and related atomistic computer simulations are effective ways in studying the structures and structure–property relations of glass materials. However, simulations of boron oxide (B ₂ O ₃ )‐containing oxide glasses pose a challenge due to the lack of reliable empirical potentials. This paper reports development of a set of partial charge pairwise composition‐dependent potentials for boron‐related interactions that enable simulations of multicomponent borosilicate glasses, together with some of the existing parameters. This set of potentials was tested in sodium borate glasses and sodium borosilicate glasses and it is shown capable to describe boron coordination change with glassmore » composition in wide composition ranges. Structure features such as boron N ₄ value, density, Q _n species distribution, fraction of non‐bridging oxygen around boron and silicon, total correlation function, and bond angle distribution function were calculated and compared with available experimental data. Mechanical properties of the simulated glasses calculated with the new potential also show good agreement with experiments. Therefore, this new set of potential can be used to simulate boron oxide‐containing multicomponent glasses including those with wide industrial and technology applications.« less

Vitrification is a widely accepted method to immobilize nuclear waste. Detailed structural information is critical to understand the physical and chemical behaviors, including the long-term chemical durability of these glasses that are to be stored in geological storage sites. High-energy X-ray diffraction studies are used to obtain accurate structural information of the International Simple Glass (ISG), a six-component borosilicate model nuclear waste glass with a composition (mol%) of 60.2SiO₂-16.0B₂O₃-12.6Na₂O-3.8Al₂O₃-5.7CaO-1.7ZrO₂. Classical molecular dynamics (MD) simulations with recently developed effective potentials are utilized to generate structural models of the ISG glasses to provide interpretation of the high-energy X-ray structural data and additionalmore » short and medium range structural details. The atomic structure model generated from MD simulations shows an excellent agreement with high-energy X-ray diffraction, measured for the first time for ISG, with an R_x value of 5.5%. Systematic investigations of the medium-range structural features of ISG are performed, built upon recent study of the short-range structure. In particular, the glass former cation oxygen polyhedral analysis show that [SiO₄]-[SiO₄] and [BO₄]-[SiO₄] connections are the most abundant while [ZrO₆] octohedra having higher probability to connect to another [SiO₄], [BO₄] and [BO₃]. Moreover, we have studied the ISG surface structure using MD simulations. Noticeable differences in terms of both chemical compositions and structural features between the bulk and surface are observed. There exist increased small-sized rings, under coordinated Si species, high 3-coordinated boron content (86.5%) on the surface. Very importantly, sodium ions are found to be enriched on the glass surface, consistent with recent experimental results. Furthermore, both bulk and surface structural information are discussed in terms of the dissolution and corrosion mechanisms of ISG and related nuclear waste glasses.« less

We present that borosilicate glasses form an important glass forming system in both glass science and technologies. The structure and property changes of borosilicate glasses as a function of thermal history in terms of cooling rate during glass formation and simulation system sizes used in classical molecular dynamics (MD) simulation were investigated with recently developed composition dependent partial charge potentials. Short and medium range structural features such as boron coordination, Si and B Qn distributions, and ring size distributions were analyzed to elucidate the effects of cooling rate and simulation system size on these structure features and selected glass propertiesmore » such as glass transition temperature, vibration density of states, and mechanical properties. Neutron structure factors, neutron broadened pair distribution functions, and vibrational density of states were calculated and compared with results from experiments as well as ab initio calculations to validate the structure models. The results clearly indicate that both cooling rate and system size play an important role on the structures of these glasses, mainly by affecting the 3B and 4B distributions and consequently properties of the glasses. It was also found that different structure features and properties converge at different sizes or cooling rates; thus convergence tests are needed in simulations of the borosilicate glasses depending on the targeted properties. Lastly, the results shed light on the complex thermal history dependence on structure and properties in borosilicate glasses and the protocols in MD simulations of these and other glass materials.« less