Search NASA⌕ Search

SEARCH · Search NASA

Results for “Silicates”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

29 Si solid state MAS NMR study on leaching behaviors and chemical stability of different Mg-silicate structures for CO 2 sequestration

Silicon is one of the most earth abundant elements, and thus, the fate and reactivity of silicate materials are often important for various energy and environmental technologies including carbon sequestration, where CO 2 is captured and stored as a thermodynamically stable solid carbonate phase. Thus, understanding the structures and chemistries of different silicate phases has become an important research aim. Here in this study, the changes in the silicate structures (Q 0 –Q 4 ) of heat-treated Mg-bearing mineral (serpentine) exposed to a CO 2 -water system (carbonic acid) was investigated using 29 Si MAS NMR, XRPD and ICP-OES and the identified structures were employed to explain complex leaching behaviors of silicate materials. The 29 Si MAS NMR and XRPD analysis indicated that the heat-treated serpentine is a mixture of amorphous (Q 1 : dehydroxylate I, Q 2 : enstatite, Q 4 : silica) and crystalline (Q 0 : forsterite, Q 3 : dehydroxylate II and serpentine) phase, while natural serpentine mineral has single crystalline Q 3 silicate structure. The leaching experiments showed that both Mg and Si in the amorphous silicate structures (Q 1 : dehydroxylate I, Q 2 : enstatite) are more soluble than those in crystalline phase (Q 0 : forsterite, Q 3 : dehydroxylate II and serpentine). Therefore, tuning the silicate structure towards Q 1 and Q 2 would significantly improve carbon sequestration potential of silicate minerals, whereas silicate materials with Q 3 structure would provide great chemical stabilities in acidic conditions. The solubilities of silicate structures were in the order of Q 1 (dehydroxylate I) > Q2 (enstatite) >> Q0 (forsterite) > Q3 (dehydroxylate II) > Q 3 (serpentine) and this finding can be used to better design a wide range of energy and environmental materials and reaction systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Probing Carbon Mineralization Mechanisms in Pore and Bulk Fluids by Harnessing Architected Calcium Silicates

The ability to synthesize materials with well-controlled pore structures gives us unprecedented control over probing fluid interactions with reactive interfaces and advancing calibrated insights into coupled chemo-morphological interactions. One of the primary challenges in developing crystalline silicate materials lies in achieving ordered pore structures. Existing approaches of producing amorphous mesoporous metal silicates via sol–gel methods and heat-treatment of these materials to produce crystalline phases cause the pore structures in the amorphous phases to collapse. To overcome this challenge, carbon coating of amorphous mesoporous calcium silicate particles is carried out to retain the pore structure, while the material is heated to produce crystalline calcium silicate with calcium sulfate inclusions. The pore diameter in these materials is about 3.9 nm, with a surface area and a pore volume of 28.75 m 2 /g and 0.092 cm 2 /g, respectively. The mechanisms of carbon mineralization are investigated by reacting architected calcium silicate with 1 M Na 2 CO 3 and monitoring the evolution in the structural phases using operando wide-angle X-ray scattering (WAXS) measurements. Formation of stable calcium carbonate polymorph or calcite and metastable calcium carbonate polymorph or vaterite in pore and bulk fluids, respectively, resulting from the reaction between Na 2 CO 3 and CaSiO 3 , are noted. The mechanisms associated with carbon mineralization are delineated using ReaxFF molecular dynamics (MD) simulations. The surface dissolution reaction is initiated by 2H + ions that replace a Ca 2+ ion in the Ca–silicate matrix. Ca 2+ ions in the solution initially react with water to form calcium hydroxide and eventually form calcium (bi)carbonate. A slow and gradual increase in the formation of sodium silicate in the solution resulting from the reactions of silicic acid or the silicon dioxide reaction with sodium hydroxide is noted. When carbon mineralization occurs in environments bearing interfacial fluids, calcite is the dominant calcium carbonate polymorph, as determined using experiments with pore fluids and molecular-scale simulations. In conclusion, these studies provide fundamental insights into the mechanisms underlying the carbon mineralization of calcium silicate informed by experiments and molecular-scale simulations.

Calcium↗
DOE OSTI2020

Influences of Mg Ions and Al Ions on Deposition Behavior of Supersaturated Silicic Acid - 20261

Silicic acid can change into various species such as deposit, colloidal, and soluble form under certain conditions: for example, gradual pH changes around geological disposal repositories. In particular, silicic acid deposits narrow groundwater flow paths, thereby restricting the flow of groundwater. Such effects may contribute to a delay in radionuclide migration. Previous studies have reported that silicic acid supersaturated in sodium and calcium ions can narrow the flow paths. This study examined the effect of magnesium (Mg) and aluminum (Al) ions on the deposition of supersaturated silicic acid on the solid phase. The deposition behavior of supersaturated silicic acid on the solid phase was investigated in a batch system while changing the initial supersaturation concentration, the coexisting ion concentration, and the amount of added solid phase. In the presence of Mg ions, silicic acid precipitation increased with increasing Mg ion concentration due to the influence of the electrolyte, and crosslinking reactions. The coexistence of Al ions significantly suppressed silicic acid deposition due to the formation of a covering on the solid surface during Al ion hydrolysis under conditions of high Al concentration (1.4 mM). On the basis of these experimental results, the apparent deposition rate constant of silicic acid was estimated to be in the range from 10{sup -11} to 10{sup -10} m/s. The Damkoehler number (the dimensionless number expressed by the ratio of the characteristic time of mass transfer and the chemical reaction constant) was sufficient to cause the clogging of the host rock by deposition. This suggests that supersaturated silicic acid deposition, regardless of the presence of Mg ions (and Al ions), contributes to improving barrier performance against radionuclide migration in the host rock. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Estimation of Sorption Behavior of Europium to Calcium Silicate Hydrate by Thermal Analysis - 20172

High-level radioactive waste generated from reprocessing processes will be disposed in a geological disposal that is deep underground. Calcium silicate hydrate is formed as a secondary mineral in the repository via the reaction of calcium ions leaching from the cement and the silicic acid dissolved from the host rocks. In this study, the sorption behavior of europium, which is also known as a chemical analog of americium, on calcium silicate hydrate was examined. The solid phase was analyzed via thermogravimetry- differential thermal analysis, differential scanning calorimetry, Raman spectroscopy, fluorescence spectroscopy, and fluorescence lifetime analysis. To understand the behavior of calcium silicate hydrate under conditions saturated with groundwater, the study samples were analyzed using the above methods, except for deferential scanning calorimetry, without a drying process. The results of the differential scanning calorimetry analyses showed that the temperatures for the endothermic peaks of calcium silicate hydrate for the lower Ca/Si ratios were larger than those for the higher Ca/Si ratios for coexistent conditions of europium ions. Hence, the relationships of the temperatures of the endothermic peaks and the Ca/Si ratios were reversed by the addition of europium ions. These study results suggest that calcium silicate hydrate can immobilize europium ions, thus preventing the production of colloidal forms. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Evidence of formation of an amorphous magnesium silicate (AMS) phase during alkali activation of (Na-Mg) aluminosilicate glasses

Highlights: • No hydrotalcite-group phase is formed as Al preferentially incorporates into zeolites. • Mg forms an amorphous Mg silicate (AMS) phase. • AMS is a more depolymerized phase when compared to Mg silicate hydrate (M-S-H). • High charge density of Mg effectively stabilizes depolymerized silicates in AMS phase. • Existence of AMS phase should be considered as possibility in Mg-rich AAMs. There is some ambiguity regarding the fate of Mg during the alkali activation of Mg-rich precursors within the broader field of alkali activated materials (AAMs). The present work addresses this issue by studying the reaction products in AAMs synthesized from (Na-Mg) aluminosilicate glasses. Here, instead of magnesium silicate hydrate (M-S-H) phase, Mg exclusively forms an amorphous magnesium silicate (AMS) phase. Compared to M-S-H, AMS is a more depolymerized phase, which has not previously been well documented. The formation of AMS seems to be driven by the high charge density of the Mg cation which effectively stabilizes the depolymerized silicate species. We also show that the lack of hydrotalcite-group phases is due to aluminum depletion by zeolite formation. This work highlights the need to consider the existence of the AMS phase in Mg-containing AAMs, especially in complex systems, where its identification may be difficult.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A numerical examination of the effect of sulfide dissolution on silicate weathering

Silicate weathering regulates climate on geological time scales as a net, climate-sensitive sink of atmospheric CO 2 . On the other hand, sulfuric acid produced through sulfide dissolution affects silicate weathering, diminishing its effectiveness as a climate regulator based on evidence from river chemistry. This study takes a theoretical approach to quantitatively examine the effect of sulfide dissolution, coupling a one-dimensional model of pyrite weathering to that of albite transformation to kaolinite. The coupled model reveals that when the reaction front of sulfide is deeper than that of silicates reacting with CO 2 , the silicate-weathering feedback on climate is not directly affected by sulfide dissolution, but only indirectly through oceanic processes such as reverse weathering and carbonate deposition. In turn, when sulfide dissolution occurs within zones of CO 2 -silicate reactions close to the surface, the feedback between climate and silicate weathering can be significantly weakened because sulfuric acid competes effectively with carbonic acid in mineral dissolution. As the depth of the sulfide reaction front can be changed by variations in water table depth, the latter feedback weakening could have occurred in concert with climate-change-driven water table fluctuations, amplifying atmospheric CO 2 variations during the Phanerozoic.

58 GEOSCIENCES↗
DOE OSTI2020

In situ X-ray diffraction of silicate liquids and glasses under dynamic and static compression to megabar pressures

Properties of liquid silicates under high-pressure and high-temperature conditions are critical for modeling the dynamics and solidification mechanisms of the magma ocean in the early Earth, as well as for constraining entrainment of melts in the mantle and in the present-day core–mantle boundary. In this paper, we present in situ structural measurements by X-ray diffraction of selected amorphous silicates compressed statically in diamond anvil cells (up to 157 GPa at room temperature) or dynamically by laser-generated shock compression (up to 130 GPa and 6,000 K along the MgSiO 3 glass Hugoniot). The X-ray diffraction patterns of silicate glasses and liquids reveal similar characteristics over a wide pressure and temperature range. Beyond the increase in Si coordination observed at 20 GPa, we find no evidence for major structural changes occurring in the silicate melts studied up to pressures and temperatures exceeding Earth’s core mantle boundary conditions. This result is supported by molecular dynamics calculations. Our findings reinforce the widely used assumption that the silicate glasses studies are appropriate structural analogs for understanding the atomic arrangement of silicate liquids at these high pressures.

58 GEOSCIENCES↗
DOE OSTI2024

Unlocking High Capacity and Reversible Alkaline Iron Redox Using Silicate-Sodium Hydroxide Hybrid Electrolytes

Alkaline iron (Fe) batteries are attractive due to the high abundance, low cost, and multiple valent states of Fe but show limited columbic efficiency and storage capacity when forming electrochemically inert Fe 3 O 4 on discharging and parasitic H 2 on charging. Herein, sodium silicate is found to promote Fe(OH) 2 /FeOOH against Fe(OH) 2 /Fe 3 O 4 conversions. Electrochemical experiments, operando X-ray characterization, and atomistic simulations reveal that improved Fe(OH) 2 /FeOOH conversion originates from (i) strong interaction between sodium silicate and iron oxide and (ii) silicate-induced strengthening of hydrogen-bond networks in electrolytes that inhibits water transport. Furthermore, the silicate additive suppresses hydrogen evolution by impairing energetics of water dissociation and hydroxyl de-sorption on iron surfaces. In conclusion, this new silicate-assisted redox chemistry mitigates H 2 and Fe 3 O 4 formation, improving storage capacity (199 mAh g -1 in half-cells) and coulombic efficiency (94 % after 400 full-cell cycles), paving a path to realizing green battery systems built from earth-abundant materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molten Flux Synthesis of Plutonium (IV) Silicates K 2 PuSi 2 O 7 and Rb 2 PuSi 6 O 15

Few plutonium containing silicates have been reported in the literature and only one other tetravalent plutonium silicate structure has been reported using the flux growth method. With the goal of understanding actinide incorporation into silicates, we describe herein two new plutonium (IV) silicates, Rb 2 PuSi 6 O 15 and K 2 PuSi 2 O 7 , synthesized using a mixed alkali chloride/alkali fluoride flux growth. Both structures crystallize in the C2/c space group with lattice parameters a=24.2980(10) Å, b=7.1466(3) Å, c=17.2025(8) Å and $\blacksquare$96.6560(10)° for Rb 2 PuSi 6 O 15 , and a=9.9478(2) Å, b=5.59820 (10) Å, c=13.1848(2) Å, and $\blacksquare$105.7400(10)° for K 2 PuSi 2 O 7 . Comparisons to the synthesis of other reported alkali actinide silicates (Ce, U, and Th) are made. Raman spectra and DFT calculations of formation enthalpies and vibrational modes provide confirmation of structure refinements and enhanced product characterization for comparison with other materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

New interaction potentials for alkaline earth silicate and borate glasses

New interaction potentials were developed for molecular dynamics simulations to study the role of Mg and Ca in modifying the structure and properties of alkaline earth silicates and borates. Competition between the depolymerization of the silica network and the formation of new bonds between oxygen atoms and modifiers leads to the enhancement of the elastic moduli with increasing modifier content in alkaline earth silicate glasses. Compared with calcium silicate, the higher elastic moduli of magnesium silicate result from a higher connectivity of the overall glass network due to the incorporation of fourfold coordinated magnesium and a more rigid connection between oxygen atoms and modifiers. In contrast to the silicates, the effect of modifier on the elastic moduli of alkaline earth borates is dominated by the formation of fourfold coordinated boron (N 4 ). Calcium borate with higher N 4 shows a more rigid network structure and higher elastic moduli.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Investigation of N in Ammonium-bearing Silicates with Electron Probe Microanalysis (EPMA)

Nitrogen is a primary constituent of Earth’s atmosphere and an essential component of many biological processes. Recent developments in understanding N incorporation into silicate minerals suggest that the solid earth and lithosphere could provide a substantial reservoir for the planetary N budget. Furthermore, determining the concentration and distribution of low abundances of N in silicates presents a significant challenge in microanalysis. We tested capabilities for high spatial resolution quantification of low levels of N using EPMA with the goal of developing standards to be used in coupled N and noble gas analyses of silicate minerals. We investigated the concentration and distribution of N in a suite of ammonium-bearing silicate minerals that reflect the dominant N-bearing phases of the lithosphere: tobelite, buddingtonite, ammonioleucite, tsaregorodtsevite, and hyalophane. These minerals form during NH4-rich alteration of K-bearing silicate minerals, and the bonding environment of N in each mineral is likely to differ based on crystal structure. Using an optimized analytical routine for N quantification enables sensitivity approaching 0.02 weight% N and high-resolution quantitative mapping of sub-micron variations in N contents.

58 GEOSCIENCES↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
DOE OSTI2020

Reversible Formation of Silanol Groups in Two-Dimensional Siliceous Nanomaterials under Mild Hydrothermal Conditions

Monitoring the effects of mild hydrothermal conditions, in situ, on siliceous materials remains challenging using surface science techniques, which often require electrically conductive substrates. The emergence of two-dimensional (2-D) siliceous nanomaterials deposited on metal single crystals overcomes this limitation. Here, we use infrared reflection absorption spectroscopy (IRRAS) to study the effects of mild hydrothermal conditions, in situ, on 2-D model systems, namely, all-Si MFI nanosheets supported on Au(111) and a polymorphous bilayer silicate supported on Ru(0001). We find that the formation of silanol groups (SiOH) occurs at 473 and 573 K under a H 2 O pressure of 3 mbar in the MFI nanosheets, but not in the polymorphous bilayer silicate. The effects of mild hydrothermal conditions are reversible in the MFI nanosheets and do not result in framework degradation. Implications shown here provide a fundamental understanding of the impact of mild hydrothermal conditions on the 2-D siliceous nanomaterials and serve as a starting point when considering these effects on three-dimensional (3-D) ones.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Hydration and reaction mechanisms on sodium silicate glass surfaces from molecular dynamics simulations with reactive force fields

Recent development of reactive force fields have enabled molecular dynamics simulations of interactions between silicate glasses and water at the atomistic scale. While multicomponent silicate glasses encompass a wide variety of compositions and properties, one common structural feature in these glasses is the combination of the network structure that is made up of silica tetrahedra linked through corner sharing interspersed with network modifiers like alkali and alkaline-earth ions that break up the Si–O–Si linkages by forming nonbridging oxygen. In reactions with water, ion exchange between alkali ions in the glass and proton or hydronium in the solution, as well as hydrolysis reaction of the Si–O–Si linkages and subsequent silanol formation, is observed and well documented. Furthermore, we have used a set of recently developed reactive force field to investigate the reactions between water and the surfaces of silica and sodium silicate glasses of different compositions for reactions up to 8 nanoseconds. Our results indicate sodium leaching into water and diffusion of water molecules up to 25 Å into the glass surface. We examined the structural and compositional changes inside the glass and around the diffused ions and use these to explain the rates of silanol formation at the surface. We also observed proton transport in the glass which has an indirect influence on the silanol formation rates. While the surface of the glass was rough to start with, it undergoes further modification into a hydrated gel-like structure in the glass for up to 5 Å in the higher alkali containing glasses. It was found that the leached sodium ions remain close to the interface and that fragments of silicate network from the surface is capable of dislodging from the bulk glass and enter the aqueous solution. These simulations thus provide insights into the formation and structure of an alteration layers commonly observed in multicomponent silicate glasses corroded in aqueous solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LionGlass™: A low‐melting, carbonate‐free alternative to soda lime silicate glass

LionGlass™ is Penn State University's patent-pending alumino-silicophosphate glass compositional family that offers, for the first time, a viable low-melting, carbonate-free alternative to soda lime silicate glass by (a) lowering the melting temperature by up to 400°C compared to soda lime silicate glass and (b) eliminating the use of carbonate batch materials. The viscosity curve of LionGlass is shifted downward in temperature compared to soda lime silicate glass, enabling melting at much lower temperatures. LionGlass also has significantly improved resistance to crack initiation compared to soda lime silicate glass. Other properties, such as thermal expansion coefficient and optical dispersion curve, are comparable to those of traditional soda lime silicate glass.

Lion-Glass↗
DOE OSTI2021

Initial hydration process of calcium silicates in Portland cement: A comprehensive comparison from molecular dynamics simulations

As the main components of Portland cement, calcium silicates show substantial differences in their hydration reactivity which have not been fully explained. A comprehensive comparison of the initial hydration process of calcium silicates, namely dicalcium silicate (C{sub 2}S) and tricalcium silicate (C{sub 3}S), was conducted using molecular dynamics simulations. The initial hydration process was divided into three stages using cut-off times of 0.001 ns and 3 ns. The hydration of M{sub 3}-C{sub 3}S (010) was more evident than that of β-C{sub 2}S (100), supported by the hydroxylation degree, radial distribution function, atomic density profile, etc. The coordination number of the surface Ca atoms might be the underlying reason for such a difference. Interactions between cement surfaces and water molecules were mainly characterised by solid OH bonding and Ca-water O bonding. Dissolution of Ca atoms was observed, although quite scarce, while no dissolution of Si atoms was observed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Searching for correlations between vibrational spectral features and structural parameters of silicate glass network

Infrared (IR) and Raman spectroscopic features of silicate glasses are usually 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 [Na 2 O] x [Al 2 O 3 ] 2 [SiO 2 ] 98- x with x = 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. Hence, the bond length and angle distributions as well as Si– Q n species 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 the Q n speciation; and (c) assignment of the Raman bands in the 420-600 cm -1 region to the bending modes of (SiO) n rings 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 result 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.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deciphering the non-linear impact of Al on chemical durability of silicate glass

The role of Al in aluminosilicate glasses remains somewhat a mystery: at low concentrations, it increases the resistance to hydrolysis of the glass, whereas at high concentrations an opposite effect is observed. To understand the origin of the phenomenon on a fundamental atomistic scale, we performed 577 MD simulations and applied potential mean force (PMF) calculations to estimate the activation barriers for hydrolysis and to statistically correlate them with local structural features of the glass. Models of pure silicate and aluminosilicate glasses are constructed and investigated. PMF simulation results are further validated by the experimental measurements and revealed that Al is very easy to dissociate, but it also increases the glass chemical durability through significantly increasing both the strength of Si and network connectivity of the glass. In contrast, at high Al concentration, preferential dissolution of Al weakens the silicate network, which it supposes to strengthen, and so the glass resistance becomes poor. Through PMF calculations, we evaluated the activation barriers for dissociating bonds around Al as 0.49 eV, which is less than a half of the energy to dissociate bonds around Si in pure silicate (1.22 eV) and around Si in aluminosilicate glass (1.34 eV), all these energy differences being statistically significant. Molecular structural level investigation revealed that Si with Al as a second neighbor in the glass network has a significantly higher activation energy for dissociation than Si in pure silicate glass. Furthermore, the proposed approach opens the way to the development of quantitative predictive models of glass durability.

36 MATERIALS SCIENCE↗