Structural origins of the Mixed Alkali Effect in Alkali Aluminosilicate Glasses: Molecular Dynamics Study and its Assessment

Accurate and Transferable Machine Learning Potential for Molecular Dynamics Simulation of Sodium Silicate Glasses

An accurate and transferable machine learning (ML) potential for the simulation of binary sodium silicate glasses over a wide range of compositions (from 0 to 50% Na2O) was developed. The potential energy surface is approximated by the sum of atomic energy contributions mapped by a neural network algorithm from the local geometry comprising information on atomic distances and angles with neighboring atoms using the DeePMD code [Wang, H. Comput. Phys. Commun. 2018, 228, 178-184]. Our model was trained on a large data set of total energies and atomic forces computed at the density functional theory level on structures extracted from classical molecular dynamics (MD) simulations performed at several temperatures from 300 to 3000 K. This allows for the generation of a robust and transferable ML potential applicable over the full compositional range of glass formability at different temperatures that outperforms the empirical potentials available in the literature in reproducing structures and properties such as bond angle distribution, total distribution functions, and vibrational density of state. The generality of the approach enables the future training of a potential with other or more elements allowing for simulations of structures, properties, and behavior of ternary and multicomponent oxide glasses with nearly ab initio accuracy at a fraction of the computational cost.

Ionic Conductivity of K-ion Glassy Solid Electrolytes of K2S-P2S5-KOTf System

Ternary glassy electrolytes containing K2S as a glass modifier and P2S5 as a network former are synthesized by introducing a new type of complex and asymmetric salt, potassium triflate (KOTf), to obtain unprecedented K+ ion conductivity at ambient temperature. The glasses are synthesized using a conventional quenching technique at a low temperature. In general, alkali ionic glassy electrolytes of ternary systems, specifically for Li+ and Na+ ion conductivity, have been studied with the addition of halide salts or oxysalts such as M2SO4, M2SiO4, M3PO4 (M = Li or Na), etc. We introduce a distinct and complex salt, potassium triflate (KOTf) with asymmetric anion, to the conventional glass modifier and former to synthesize K+-ion-conducting glassy electrolytes. Two series of glassy electrolytes with a ternary system of (0.9-x)K2S-xP2S5-0.1KOTf (x = 0.15, 0.30, 0.45, 0.60, and 0.75) and z(K2S-2P2S5)-yKOTf (y = 0.05, 0.10, 0.15, 0.20, and 0.25) on a straight line of z(K2S-2P2S5) are studied for their K+ ionic conductivities by using electrochemical impedance spectroscopy (EIS). The composition 0.3K2S-0.6P2S5-0.1KOTf is found to have the highest conductivity among the studied glassy electrolytes at ambient temperature with the value of 1.06 × 10-7 S cm-1, which is the highest of all pure K+-ion-conducting glasses reported to date. Since the glass transition temperatures of the glasses are near 100 °C, as demonstrated by DSC, temperature-dependent conductivities are studied within the range of 25 to 100 °C to determine the activation energies. A Raman spectroscopic study shows the variation in the structural units PS43-, P2S74-, and P2S64- of the network former for different glassy electrolytes. It seems that there is a role of P2S74- and P2S64- in K+-ion conductivity in the glassy electrolytes because the spectroscopic results are compatible with the composition-dependent, room-temperature conductivity trend.

Influence of Mixed Na2O/K2O on Chemical Durability and Spectral Properties of P2O5-Al2O3-BaO-K2O-Na2O-Nd2O3 Phosphate Glasses

A series of 56P2O5-7.5Al2O3-5.9BaO-(28.56-x)K2O-xNa2O-1.51Nd2O3 phosphate glasses with different Na/(Na+K) ratios, which were specially designed for high-power laser application, were prepared by a high-temperature melting method. Except for the density, refractive index, glass transition temperature, and DC conductivity, the chemical durability and spectral properties, as emphasized by high-power and high-energy laser material, were further measured and analyzed. Regarding the chemical durability, the dissolution rates of these glasses do not show an evident mixed alkali effect with increasing the Na/(Na+K) ratio, although the effect is obvious for the glass transition temperature and DC conductivity. To better understand the nature of the dissolution mechanism, the ionic release concentrations of every element are determined. Both Na and K undergo ion exchange, but the ion exchange rate of K is much larger than that of Na. In terms of the spectral properties, the J-O parameters, emission cross-section, radiation lifetime, fluorescence lifetime, effective bandwidth, fluorescence branching ratio, and quantum efficiency are determined from absorption and emission spectra. The trend of Ω2 deviating from linearity indicates that the coordination environment symmetry of Nd3+ ions and the covalence of Nd-O also present an evident mixed alkali effect. The most important finding is that the emission cross-section and fluorescence lifetime of Nd3+ ions at 1053 nm were not affected by the change in the Na/K ratio. According to the above experimental results, the optimized value of the Na/K ratio was determined, based on which the 56P2O5-7.5Al2O3-5.9BaO-(28.56-x)K2O-xNa2O-1.51Nd2O3 glass maintains a high emission cross-section with good chemical durability.

The Structure of Gd3+-Doped Li2O and K2O Containing Aluminosilicate Glasses from Molecular Dynamics Simulations

Understanding the atomic structure of glasses is critical for developing new generations of materials with important technical applications. In particular, the local environment of rare-earth ions and their distribution and clustering is of great relevance for applications of rare earth-containing glasses in photonic devices. In this work, the structure of Gd₂O₃ doped lithium and potassium aluminosilicate glasses is investigated as a function of their network modifier oxide (NMO-Li₂O, K₂O) to aluminum oxide ratio using molecular dynamics simulations. The applied simulation procedure yields a set of configurations, the so-called inherent structures, of the liquid state slightly above the glass transition temperature. The generation of a large set of inherent structures allows a statistical sampling of the medium-range order of the Gd³⁺ ions with less computational effort compared to other simulation methods. The resulting medium-range atomic structures of network former and modifier ions are in good agreement with experimental results and simulations of similar glasses. It was found that increasing NMO/Al ratio increases the network modifier coordination number with non-bridging oxygen sites and reduces the overall stability of the network structure. The fraction of non-bridging oxygen sites in the vicinity of Gd³⁺ ions increases considerably with decreasing field strength and increasing concentration of the network modifier ions. These correlations could be confirmed even if the simulation results of alkaline earth aluminosilicate glasses are added to the analysis. In addition, the structure predictions generally indicate a low driving force for the clustering of Gd³⁺. Here, network modifier ions of large ionic radii reduce the probability of Gd-O-Gd contacts.

Disclosing crystal nucleation mechanism in lithium disilicate glass through molecular dynamics simulations and free-energy calculations

Unraveling detailed mechanism of crystal nucleation from amorphous materials is challenging for both experimental and theoretical approaches. In this study, we have examined two methods to understand the initial stage of crystal precipitation from lithium disilicate glasses using molecular dynamics simulations. One of the methods is a modified exploring method to find structurally similar crystalline clusters in the glass models, enabling us to find three different embryos, such as Li₂Si₂O₅ (LS₂), Li₂SiO₃ (LS) and Li₃PO₄ (LP), in the 33Li₂O·66SiO₂·1P₂O₅ glass (LS₂P1), in which P₂O₅ is added as a nucleating agent. Interestingly, LS₂ and LP crystals were found inside the LS₂P1 glass while LS crystal appeared on the glass surface, which agrees with experimental observations. The other method is free energy calculation using a subnano-scale spherical crystal embedded in the glass model. This method, which we called Free-Energy Seeding Method (FESM), allows us to evaluate free energy change as a function of crystal radius and to identify critical size of the crystal precipitation. The free energy profiles for LS and LS₂ crystal nuclei in the LS₂ glass models possess maximum energy at a critical radius as expected by classical nucleation theory. Furthermore, the critical radius and the energy barrier height agree well with recent experimental investigation, proving the applicability of this method to design glass–ceramics by atomistic modeling.

Ab Initio Study of Hydrolysis Effects in Single and Ion-Exchanged Alkali Aluminosilicate Glasses

Hydrolysis in alkali-doped aluminosilicate glasses is one of the most complicated mechanisms in glass science. There remain many fundamental and unresolved issues with implications on their potential applications. Herein, we address this challenge by carrying out detailed calculations on the structure and properties of both anhydrate (dry) and hydrated alkali aluminosilicate glasses on carefully constructed models. Specifically, the Na-, (Na + K)-, and K-doped aluminosilicate glasses with compositions (SiO₂)_0.6(Al₂O₃)_0.2(Na₂O)_0.2 - x(K₂O)_x (x = 0, 0.10, and 0.20) are simulated using ab initio molecular dynamics (AIMD). The local short- and intermediate-range order in these glasses is analyzed in terms of atomic pair distribution, coordination number, bond length, and bond angle distributions to delineate the subtle variations due to different alkali sizes and hydrolysis. The electronic structure, interatomic bonding, mechanical, and optical properties for these models are calculated and validated with available experimental data. We use the novel concept of total bond order density (TBOD), the quantum mechanically derived metric, to characterize the internal cohesion and strength in the simulated glasses. Detailed analysis of the hydrolysis mechanism enables us to provide information on the complex interplay of various participating elements and their interactions at the atomic level. Such detailed information provides a new platform of knowledge, which is crucial for understanding the issues related to glass corrosion and durability, and ways and means for their special applications in commercial glass products. Both undissociated molecular water and dissociated water in the form of hydroxyl groups exist in the hydrated models in the presence of alkali ions. For the first time, we observed the opposite mixed alkali effect in the Poisson's ratio for anhydrate and hydrated glasses.

Ionic self-diffusion and the glass transition anomaly in aluminosilicates

The glass transition temperature (Tg) is the temperature after which a supercooled liquid undergoes a dynamical arrest. Usually, glass network modifiers (e.g., Na2O) affect the behavior of Tg. However, in aluminosilicate glasses, the effect of different modifiers on Tg is still unclear and shows an anomalous behavior. Here, based on molecular dynamics simulations, we show that Tg decreases with increasing charge balancing cation field strength (FS) in the aluminosilicate glasses, which is an anomalous behavior as compared to other oxide glasses. The results show that the origins of this anomaly come from the dynamics of the supercooled liquid above Tg, which in turn is correlated to pair excess entropy. Our results deepen our understanding of the effect of different modifiers on the properties of the aluminosilicate glasses.