Effects of system size and cooling rate on the structure and properties of sodium borosilicate glasses from molecular dynamics simulations

Optimal Molecular Dynamics System Size for Increased Precision and Efficiency for Epoxy Materials

Molecular dynamics (MD) simulation is an important tool for predicting thermo-mechanical properties of polymer resins at the nanometer length scale, which is particularly important for efficient computationally driven design of advanced composite materials and structures. Because of the statistical nature of modeling amorphous materials on the nanometer length scale, multiple MD models (replicates) are typically built and simulated for statistical sampling of predicted properties. Larger replicates generally provide higher precision in the predictions but result in higher simulation times. Unfortunately, there is insufficient information in the literature to establish guidelines between MD model size and the resulting precision in predicted thermo-mechanical properties. The objective of this study was to determine the optimal MD model size of epoxy resin to balance efficiency and precision. The results show that an MD model size of 15,000 atoms provides for the fastest simulations without sacrificing precision in the prediction of mass density, elastic properties, strength, and thermal properties of epoxy. The results of this study are important for efficient computational process modeling and integrated computational materials engineering (ICME) for the design of next-generation composite materials for demanding applications.

Insights into Lithium Sulfide Glass Electrolyte Structures and Ionic Conductivity via Machine Learning Force Field Simulations

Sulfide-based solid electrolytes (SEs) are important for advancing all-solid-state batteries (ASSBs), primarily due to their high ionic conductivities and robust mechanical stability. Glassy SEs (GSEs) comprising mixed Si and P glass formers are particularly promising for their synthesis process and their ability to prevent lithium dendrite growth. However, to date, the complexity of their glassy structures hinders a complete understanding of the relationships between their structures and properties. This study introduces a new machine learning force field (ML-FF) tailored for lithium sulfide-based GSEs, enabling the exploration of their structural characteristics, mechanical properties, and lithium ionic conductivities. Using molecular dynamic (MD) simulations with this ML-FF, we explore the glass structures in varying compositions, including binary Li2S-SiS2 and Li2S-P2S5 as well as ternary Li2S-SiS2-P2S5. Our simulations yielded consistent results in terms of density, elastic modulus, radial distribution functions, and neutron structure factors compared to DFT and experimental work. Our findings reveal distinct local environments for Si and P within these glasses, with most Si atoms in edge-sharing configurations in Li2S-SiS2 and a mix of corner- and edge-sharing tetrahedra in the ternary Li2S-SiS2-P2S5 composition. For lithium ionic conductivity at 300 K, the 50Li2S-50SiS2 glass displayed the lowest conductivity at 2.1 mS/cm, while the 75Li2S-25P2S5 composition exhibited the highest conductivity at 3.6 mS/cm. The ternary glass showed a conductivity of 2.6 mS/cm, sitting between the two. Moreover, an in-depth analysis of lithium ion diffusion over the MD trajectory in the ternary glass demonstrated a significant correlation between diffusion pathways and the rotational dynamics of nearby SiS4 or PS4 tetrahedra. The ML-FF developed in this study provides an important tool for exploring a broad spectrum of solid-state and mixed former sulfide-based electrolytes.

Structural Causes of Brittleness Changes in Aluminosilicate Glasses with Different Cooling Rates

Numerous sources have already demonstrated that varying annealing rates can result in distinct toughness and brittleness in glass. To determine the underlying mechanisms driving this phenomenon, molecular dynamic (MD) simulations were employed to investigate the microstructure of aluminosilicate glasses under different cooling rates, and then uniaxial stretching was performed on them under controlled conditions. Results indicated that compared with short-range structure, cooling rate has a greater influence on the medium-range structure in glass, and it remarkably affects the volume of voids. Both factors play a crucial role in determining the brittleness of the glass. The former adjusts network connectivity to influence force transmission by manipulating the levels of bridging oxygen (BO) and non-bridging oxygen (NBO), and the latter accomplishes the objective of influencing brittleness by modifying the environmental conditions that affect the changes in BO and NBO content. The variation in the void environment results in differences in the strategies of the changes in BO and NBO content during glass stress. These findings stem from the excellent response of BO and NBO to the characteristic points of stress-strain curves during stretching. This paper holds importance in understanding the reasons behind the effect of cooling rates on glass brittleness and in enhancing our understanding of the ductile/brittle transition (DTB) in glass.

Hanford low-activity waste vitrification: A review

This paper summarizes the vast body of literature (over 200 documents) related to vitrification of the low-activity waste (LAW) fraction of the Hanford tank wastes. Details are provided on the origins of the Hanford tank wastes that resulted from nuclear operations conducted between 1944 and 1989 to support nuclear weapons production. Waste treatment processes are described, including the baseline process to separate the tank waste into LAW and high-level waste fractions, and the LAW vitrification facility being started at Hanford. Significant focus is placed on the glass composition development and the property-composition relationships for Hanford LAW glasses. Glass disposal plans and criteria for minimizing long-term environmental impacts are discussed along with research perspectives.

Vanadium Oxidation States and Structural Role in Aluminoborosilicate Glasses: An Integrated Experimental and Molecular Dynamics Simulation Study

Vanadium-containing glasses have aroused interest in several fields such as electrodes for energy storage, semiconducting glasses, and nuclear waste disposal. The addition of V₂O₅, even in small amounts, can greatly alter the physical properties and chemical durability of glasses; however, the structural role of vanadium in these multicomponent glasses and the structural origins of these property changes are still poorly understood. We present a comprehensive study that integrates advanced characterizations and atomistic simulations to understand the composition-structure-property relationships of a series of vanadium-containing aluminoborosilicate glasses. UV-vis spectroscopy, X-ray photoelectron spectroscopy, and X-ray absorption near-edge structure (XANES) have been used to investigate the complex distribution of vanadium oxidation states as a function of composition in a series of six-component aluminoborosilicate glasses. High-energy X-ray diffraction and molecular dynamics simulations were performed to extract the detailed short- and medium-range atomistic structural information such as bond distance, coordination number, bond angle, and network connectivity, based on recently developed vanadium potential parameters. It was found that vanadium mainly exists in two oxidation states: V⁵⁺ and V⁴⁺, with the former being dominant (∼80% from XANES) in most compositions. V⁵⁺ ions were found to exist in 4-, 5-, and 6-fold coordination, while V⁴⁺ ions were mainly in 4-fold coordination. The percentage of 4-fold-coordinated boron and network connectivity initially increased with increasing V₂O₅ up to around 5 mol % but then decreased with higher V₂O₅ contents. The structural role of vanadium and the effect on glass structure and properties are discussed, providing insights into future studies of sophisticated structural descriptors to predict glass properties from composition and/or structure and aiding the formulation of borosilicate glasses for nuclear waste disposal and other applications.

Influence of Magnesium on the Structure of Complex Multicomponent Silicates: Insights from Molecular Simulations and Neutron Scattering Experiments

A series of multicomponent glasses containing up to five oxides are studied using classical molecular dynamics simulations and neutron scattering experiments. The focus is on the role of magnesium in determining the structural properties of these glasses and the possible mixed effect during a sodium/magnesium substitution. Calculated structure functions (pair correlation function and structure factor) rather accurately reproduce their experimental counterpart, and we show that more fine structural features are qualitatively reproduced well, despite some discrepancies in the preferential spatial distribution between sodium and magnesium to aluminum and boron, as well as the nonbridging oxygen, distribution. The simulated systems offer a solid basis to support previous experimental findings on the composition-structure relationship, allowing for further analysis and property calculation. It is confirmed that the substitution of sodium by magnesium leads to the decrease of four-fold boron and a modification of the alkali coordinations with a significant change of the network structure. Specifically, magnesium coordination extracted from numerical simulations highlights a potential dissociation from penta- to tetra- and hexahedral units with increasing MgO contents along the glass series, which could not be resolved experimentally.

Recent Advances in Corrosion Science Applicable To Disposal of High-Level Nuclear Waste

High-level radioactive waste is accumulating at temporary storage locations around the world and will eventually be placed in deep geological repositories. The waste forms and containers will be constructed from glass, crystalline ceramic, and metallic materials, which will eventually come into contact with water, considering that the period of performance required to allow sufficient decay of dangerous radionuclides is on the order of 10⁵-10⁶ years. Corrosion of the containers and waste forms in the aqueous repository environment is therefore a concern. This Review describes the recent advances of the field of materials corrosion that are relevant to fundamental materials science issues associated with the long-term performance assessment and the design of materials with improved performance, where performance is defined as resistance to aqueous corrosion. Glass, crystalline ceramics, and metals are discussed separately, and the near-field interactions of these different material classes are also briefly addressed. Finally, recommendations for future directions of study are provided.

Experimental method to quantify the ring size distribution in silicate glasses and simulation validation thereof

Silicate glasses have no long-range order and exhibit a short-range order that is often fairly similar to that of their crystalline counterparts. Hence, the out-of-equilibrium nature of glasses is largely encoded in their medium-range order. However, the ring size distribution-the key feature of silicate glasses' medium-range structure-remains invisible to conventional experiments and, hence, is largely unknown. Here, by combining neutron diffraction experiments and force-enhanced atomic refinement simulations for two archetypical silicate glasses, we show that rings of different sizes exhibit a distinct contribution to the first sharp diffraction peak in the structure factor. On the basis of these results, we demonstrate that the ring size distribution of silicate glasses can be determined solely from neutron diffraction patterns, by analyzing the shape of the first sharp diffraction peak. This method makes it possible to uncover the nature of silicate glasses' medium-range order.

Molecular dynamics simulations of simplified sodium borosilicate glasses: the effect of composition on structure and dynamics

The fusion of valuable material properties has led to the acceptance of sodium borosilicate (NBS) glasses for nuclear waste immobilization. Although popular, the mechanisms associated with these properties are still only partially discovered and need further exploration. Bearing this in mind, the combination of experiments, molecular dynamics (MD) simulations and the Dell, Yuan and Bray model have been used to understand the role of composition variation for structural and physical aspects of vitrified borosilicate glasses. Experiments have been conducted to evaluate the macroscopic glass parameters of density (ρ), glass transition temperature (Tg) and thermal expansion coefficient (TEC). Experimentally observed trends for ρ, Tg and TEC with composition have been found in good agreement with the MD results. MD studies also provide a microscopic understanding of the glass structure and phenomena associated with the change in the glass composition. A detailed view of local structure and medium-range connectivity for the borosilicate glasses has been explored. Owing to a large B4 population, the results showed the abundant presence of BO4-BO4 connections, we hereby omit the generally accepted "B[4] avoidance rule" for glass. The relative propensity for connecting SiO4/BO3/BO4 structural motifs is in line with the predictions made by the Dell, Yuan and Bray model. Furthermore, the effects of composition on the mechanical integrity of NBS glasses, including the elastic nature, plastic distortion, yielding, breaking stress, and brittle fracture, have been explored by MD simulations. In addition, the glass dynamics have been evaluated by diffusion coefficient and the results suggest that Na+ is likely to be more mobile in the case of NBS1 as compared to NBS2 and NBS3 due to significant disruption in the glass network introduced by a larger amount of Na2O network modifier. Also, the diffusivity was reduced with increasing B2O3 due to the altered role of Na+ ions from network modifiers to charge compensators. The combined study of experiments, MD simulations and the Dell, Yuan and Bray model establish the correlation between the microscopic structure and macroscopic properties of NBS glasses with varied composition, which might be of great scientific use for future glasses in various applications including nuclear waste immobilization.

Statistical topology of bond networks with applications to silica

Whereas knowledge of a crystalline material's unit cell is fundamental to understanding the material's properties and behavior, there are no obvious analogs to unit cells for disordered materials despite the frequent existence of considerable medium-range order. This article views a material's structure as a collection of local atomic environments that are sampled from some underlying probability distribution of such environments, with the advantage of offering a unified description of both ordered and disordered materials. Crystalline materials can then be regarded as special cases where the underlying probability distribution is highly concentrated around the traditional unit cell. The H_{1} barcode is proposed as a descriptor of local atomic environments suitable for disordered bond networks and is applied with three other descriptors to molecular dynamics simulations of silica glasses. Each descriptor reliably distinguishes the structure of glasses produced at different cooling rates, with the H_{1} barcode and coordination profile providing the best separation. The approach is generally applicable to any system that can be represented as a sparse graph.

Ab Initio Molecular Dynamics of CdSe Quantum-Dot-Doped Glasses

We have probed the local atomic structure of the interface between a CdSe quantum dot (QD) and a sodium silicate glass matrix. Using ab initio molecular dynamics simulations, we determined the structural properties and bond lengths, in excellent agreement with previous experimental observations. On the basis of an analysis of radial distribution functions, coordination environment, and ring structures, we demonstrate that an important structural reconstruction occurs at the interface between the CdSe QD and the glass matrix. The incorporation of the CdSe QD disrupts the Na-O bonds, while stronger SiO₄ tetrahedra are reformed. The existence of the glass matrix breaks the stable 4-membered (4MR) and 6-membered (6MR) Cd-Se rings, and we observe a disassociated Cd atom migrated in the glass matrix. Besides, the formation of Se-Na and Cd-O linkages is observed at the CdSe QD/glass interface. These results significantly extend our understanding of the interfacial structure of CdSe QD-doped glasses and provide physical and chemical insight into the possible defect structure origin of CdSe QD, of interest to the fabrication of the highly luminescent CdSe QD-doped glasses.

Modifier clustering and avoidance principle in borosilicate glasses: A molecular dynamics study

Oxide glasses are typically described as having a random, disordered skeleton of network-forming polyhedra that are depolymerized by network-modifying cations. However, the existence of local heterogeneity or clustering within the network-forming and network-modifying species remains unclear. Here, based on molecular dynamics simulations, we investigate the atomic structure of a series of borosilicate glasses. We show that the network-modifying cations exhibit some level of clustering that depends on composition-in agreement with Greaves' modified random network model. In addition, we demonstrate the existence of some mutual avoidance among network-forming atoms, which echoes the Loewenstein avoidance principle typically observed in aluminosilicate phases. Importantly, we demonstrate that the degree of heterogeneity in the spatial distribution of the network modifiers is controlled by the level of ordering in the interconnectivity of the network formers. Specifically, the mutual avoidance of network formers is found to decrease the propensity for modifier clustering.

B2O3/SiO2 substitution effect on structure and properties of Na2O–CaO–SrO–P2O5–SiO2 bioactive glasses from molecular dynamics simulations

The effect of B₂O₃/SiO₂ substitution in SrO-containing 55S4.3 bioactive glasses on glass structure and properties, such as ionic diffusion and glass transition temperature, was investigated by combining experiments and molecular dynamics simulations with newly developed potentials.