Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface dynamics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Performing In Situ Closed-Cell Gas Reactions in the Transmission Electron Microscope

Gas reactions studied by in situ electron microscopy can be used to capture the real-time morphological and microchemical transformations of materials at length scales down to the atomic level. In situ closed-cell gas reaction (CCGR) studies performed using (scanning) transmission electron microscopy (STEM) can separate and identify localized dynamic reactions, which are extremely challenging to capture using other characterization techniques. For these experiments, we used a CCGR holder that utilizes microelectromechanical systems (MEMS)-based heating microchips (hereafter referred to as "E-chips"). The experimental protocol described here details the method for performing in situ gas reactions in dry and wet gases in an aberration-corrected STEM. This method finds relevance in many different materials systems, such as catalysis and high-temperature oxidation of structural materials at atmospheric pressure and in the presence of various gases with or without water vapor. Here, several sample preparation methods are described for various material form factors. During the reaction, mass spectra obtained with a residual gas analyzer (RGA) system with and without water vapor further validates gas exposure conditions during reactions. Integrating an RGA with an in situ CCGR-STEM system can, therefore, provide critical insight to correlate gas composition with the dynamic surface evolution of materials during reactions. In situ/operando studies using this approach allow for detailed investigation of the fundamental reaction mechanisms and kinetics that occur at specific environmental conditions (time, temperature, gas, pressure), in real-time, and at high spatial resolution.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Use of in Situ Synchrotron Techniques to Probe the Oxidized Surface of Molybdenum Nitride Oxygen Reduction Electrocatalysis

The development of active and stable earth-abundant catalysts for the oxygen reduction reaction (ORR) is needed for widespread, economic development of fuel cell technologies. Designing and optimizing these non-platinum group metals is challenging, however, because they are susceptible to composition and structure changes, including dissolution, oxidation, and corrosion, both in air and under reaction conditions. To identify the active surface, and thus understand the properties that affect activity, the catalyst surface must be characterized in situ. Herein, we utilize a grazing incidence electrochemical cell to investigate in situ composition and morphology changes of a molybdenum nitride (Mo-N) thin film catalyst using grazing incidence x-ray absorption spectroscopy (GI-XAS) and x-ray reflectivity (XRR). In rotating ring disk electrode measurements, we find that the activity, selectivity, stability, and capacitance of the Mo-N catalyst is dependent on the maximum potential to which it has been exposed. Specifically, the overpotential required to reach -2 mA cm-2geo decreases by over 90 mV when the maximum potential is increased from 0.3 to 0.8 V vs RHE (Figure 1). Because the Mo-N oxidizes rapidly in air, ex situ characterization methods including x-ray photoelectron spectroscopy and time-of-flight secondary ion mass spectrometry can provide only limited insight into these in situ catalyst changes. Using in situ GI-XAS measurements at applied potentials between 0.3 and 0.9 V vs RHE, however, we are able to determine that the surface of the film oxidizes and becomes more amorphous when exposed to increasingly higher potentials (Figure 1). Furthermore, the surface remains oxidized on the order of several hours when returned to "ORR relevant potentials" (< 0.6 V vs RHE), indicating that this surface-oxidized nitride is the active surface for ORR. Using in situ XRR measurements, we find that there is no change in surface roughness at potentials below 0.7 V vs RHE, but the film roughens significantly at 0.8 V vs RHE, correlating with ex situ measurements of Mo dissolution at this potential (Figure 1). We therefore conclude that the intrinsic activity of the Mo-N catalyst increases when exposed to potentials up to 0.7 V vs RHE, while above that potential activity enhancements are due to the exposure of more active sites through dissolution. The in situ electrochemical surface-sensitive x-ray characterization as used here is a promising methodology for understanding and leveraging surface dynamics to improve the performance of non-traditional catalysts.

Kreider, Melissa↗
DOE OSTI2023

Water Structure and Dynamics Near the Surfaces of Silicalite-1

Zeolites are crystalline microporous aluminosilicates that are commonly used as industrial sorbents and membranes. Solvent structuring and dynamics near zeolite crystal surfaces are thought to influence interfacial transport and molecular recognition processes, as well as fundamental aspects of their crystallization from solution. Here, in this work, we use molecular dynamics (MD) simulations to investigate the behavior of interfacial water near the exposed (010), (100), and (101) crystal faces of silicalite-1, one of the most widely studied zeolites. The MD simulations reveal that water’s translational and orientational order is strongly influenced by the distinct corrugations, pore apertures, and functional group distributions presented on each surface. Specifically, we observe two distinct hydration layers near each surface. Water molecules in the contact layer are hydrogen bonded to two or three exposed surface silanol groups. The relative populations of water molecules that are doubly and triply hydrogen bonded to the surface strongly depend on the distribution of the exposed silanols. The interactions with exposed surface silanols also influence water’s local orientational order and distribution across each surface. We also show that the structuring of the solvent near the different faces of silicalite-1 strongly impacts the interfacial dynamics. The translational and orientational relaxation dynamics of water are slowest in the contact layers and correlate with structural ordering near each face.

42 ENGINEERING↗
DOE OSTI2020

A Flux Detection Probe to Quantify Dynamic Groundwater-Surface Water Exchange in the Hyporheic Zone

A new probe was designed to quantify groundwater-surface water exchange in the hyporheic zone under dynamic stage condition. Current methods focus on either vertical pore water velocity or Darcy flux measurements. Both parameters must be understood to evaluate residence time and mass flux of constituents. Furthermore, most instruments are not well suited for monitoring instantaneous velocity or flux under dynamic exchange conditions. For this reason, the flux detection probe (FDP) was designed that employs electrogeophysical measurements to estimate in situ sediment porosity, which can be used to convert pore water velocity to Darcy flux. Dynamic pore water velocity is obtained by monitoring fluid conductivity and temperature along the FDP probe. Pressure sensors deployed at the top and bottom of the probe provide the additional information necessary to estimate vertical permeability. This study focuses on the use of a geophysical method to estimate pore water velocity, porosity, and permeability within a controlled soil column where simulated river water displaces simulated groundwater. The difference between probe derived and theoretical pore water velocity using natural tracers such as electrical conductivity and temperature was -4.9 and 3.9% for downward flow and 1.1 and 12.8% for upward flow, respectively. The difference in porosity calculated from mass and volume packed in the soil column and probe measure porosity ranged between -3.2% and 1.5%. Also, the calculated hydraulic conductivity differed from probe derived values by -8.9%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Detecting axion dynamics on the surface of magnetic topological insulators

Axions, initially proposed to solve the strong CP problem, have recently gained attention in condensed matter physics, particularly in topological insulators. However, detecting axion dynamics has proven challenging, with no experimental confirmations to date. In this paper, we identify the surface of magnetic topological insulators as an ideal platform for observing axion dynamics. The vanishing bulk gap at the surface allows for order O(1) variations in the axion field, making the detection of axionlike phenomena more feasible. In contrast, these phenomena are strongly suppressed in the bulk due to the small magnetic exchange gap. We investigate two-photon decay as a signature of axion dynamics and calculate the branching ratio using a perturbative approach. Our findings reveal that the photon flux emitted from the surface is in plane and orders of magnitude larger than that from the bulk, making it detectable with modern microwave technology. We also discuss potential material platforms for detecting axion two-photon decay and strategies to enhance the signal-to-noise ratio.

Gao, Zhi-Qiang↗
DOE OSTI2020

Takagi–Taupin dynamical X-ray diffraction simulations of asymmetric X-ray diffraction from crystals: the effects of surface undulations

Dynamical X-ray diffraction simulations from crystals with surface undulations are reported. The Takagi–Taupin equations are applied and used to derive results in good agreement with experimental data reported in a separate paper [Macrander, Pereira, Huang, Kasman, Qian, Wojcik & Assoufid (2020). J. Appl. Cryst. 53, 789–792]. The development of Uragami [J. Phys. Soc. Jpn, (1969), 27, 147–154] is followed. Although previous work by Olekhnovich & Olekhnovich [Acta. Cryst. (1980), A36, 22–27] treated a crystal in the shape of a round cylinder, there do not seem to be any reports of previous dynamical X-ray diffraction treatments specifically for surface undulations. The significance of the present work is that it bridges the diffraction treatment of more classical dynamical diffraction theory, which assumes a flat surface, and the simple kinematic diffraction theory. Finally, the kinematic theory has, to date, been the primary means of simulating X-ray diffraction from surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

Complex Formation and Dissociation Dynamics on Amorphous Silica Surfaces

Benzene complex formation and dissociation dynamics with silanols on the amorphous silica surfaces of nanoporous SiO 2 , from a benzene/carbon tetrachloride solution, were measured by the growth of off-diagonal peaks in the two dimensional infrared (2D IR) chemical exchange spectrum of the isolated Si-OD stretch. The presence of two types of isolated silanols, referred to as type I and II, was revealed, with dissociation time constants of 82 ps and 4.0 ps for type I and II complexes, respectively. The type I silanols are associated with the main IR absorption feature in the Si-OD stretching region, while the type II silanols give rise to a broader shoulder to lower frequency. Polarization selective pump-probe (PSPP) measurements provided the vibrational lifetimes and orientational relaxation rates of the two silanols in the CCl 4 (free) and benzene (complex) environments. The type II silanols constitute roughly 30% of the isolated silanol population and exhibit a substantially faster rate of vibrational relaxation, making the type I dynamics the dominant contribution to the PSPP and 2D IR signals. From the measured dissociation times, the enthalpies of formation for the two surface complexes were obtained; the formation of the type I complex being significantly more exothermic. Furthermore, as the type II site is preferentially removed from the amorphous silica surface with increasing activation temperature, the results provide a reasonable explanation for the increased exothermicity of benzene adsorption on silica with increasing activation temperature in previous calorimetry experiments.

14 SOLAR ENERGY↗
DOE OSTI2022

Single molecule tracking of bacterial cell surface cytochromes reveals dynamics that impact long-distance electron transport

Using a series of multiheme cytochromes, the metal-reducing bacterium Shewanella oneidensis MR-1 can perform extracellular electron transfer (EET) to respire redox-active surfaces, including minerals and electrodes outside the cell. While the role of multiheme cytochromes in transporting electrons across the cell wall is well established, these cytochromes were also recently found to facilitate long-distance (micrometer-scale) redox conduction along outer membranes and across multiple cells bridging electrodes. Recent studies proposed that long-distance conduction arises from the interplay of electron hopping and cytochrome diffusion, which allows collisions and electron exchange between cytochromes along membranes. However, the diffusive dynamics of the multiheme cytochromes have never been observed or quantified in vivo, making it difficult to assess their hypothesized contribution to the collision-exchange mechanism. Here, we use quantum dot labeling, total internal reflection fluorescence microscopy, and single-particle tracking to quantify the lateral diffusive dynamics of the outer membrane-associated decaheme cytochromes MtrC and OmcA, two key components of EET in S. oneidensis. We observe confined diffusion behavior for both quantum dot-labeled MtrC and OmcA along cell surfaces (diffusion coefficients D MtrC = 0.0192 ± 0.0018 µm 2 /s, D OmcA = 0.0125 ± 0.0024 µm 2 /s) and the membrane extensions thought to function as bacterial nanowires. We find that these dynamics can trace a path for electron transport via overlap of cytochrome trajectories, consistent with the long-distance conduction mechanism. The measured dynamics inform kinetic Monte Carlo simulations that combine direct electron hopping and redox molecule diffusion, revealing significant electron transport rates along cells and membrane nanowires.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Effect of Surface Roughness on Dynamic Stall in Pitching Motion

Dynamic stall plays a critical role in determining the performance and stability of a wide range of fluid-dynamics systems in various engineering applications. This unsteady aerodynamic phenomenon is particularly significant for maneuvering aircraft wings, jet aircraft subjected to gust encounters, helicopter rotor blades, and wind turbine blades [1–3]. The prediction of the dynamic stall vortex (DSV) is challenging due to factors such as unsteady aerodynamics, three-dimensional (3-D) effects, turbulence and flow separation, incoming gust, and surface impact effects [4–6]. Hence, advanced computational techniques and modeling approaches in computational fluid dynamics (CFD) would be required to enhance the accuracy and reliability of DSV predictions in dynamic stall scenarios. In the past, Batther and Lee [7] employed delayed detached eddy simulations (DDES) to understand the flow physics associated with the onset of dynamic stall. Their approach demonstrated that DDES achieves results comparable to those obtained from large-eddy simulations at a reduced computational cost. In another study, Khalifa et al. [8] examined the 3-D aspects of dynamic stall on a NACA 0012 airfoil using DES solvers. In conclusion, the findings underlined the superiority of 3-D simulations over two-dimensional approaches, particularly in predicting the lift coefficient values and capturing dynamic stall stages more precisely.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Insights into Spontaneous Solid Electrolyte Interphase Formation at Magnesium Metal Anode Surface from Ab Initio Molecular Dynamics Simulations

Spontaneous chemical reactivity at multivalent (Mg, Ca, Zn, Al) electrode surfaces is critical to solid electrolyte interphase (SEI) formation, and hence, directly affects the longevity of batteries. Here, we report an investigation of the reactivity of 0.5 M Mg(TFSI)2 in 1,2-dimethoxyethane (DME) solvent at a Mg(0001) surface using ab initio molecular dynamics (AIMD) simulations and detailed Bader charge analysis. Based on the simulations, the initial degradation reactions of the electrolyte strongly depend on the structure of the Mg(TFSI)2 species near the anode surface. At the surface, the dissociation of Mg(TFSI)2 species occurs via cleavage of the N-S bond for the solvent separated ion pair (SSIP) and via cleavage of the C-S bond for the contact ion pair (CIP) configuration. In the case of the CIP, both TFSI anions undergo spontaneous bond dissociation reactions to form atomic O, C, S, F, and N species adsorbed on the surface of the Mg anode. These products indicate that the initial SEI layer formed on the surface of the pristine Mg anode consists of a complex mixture of multiple components such as oxides, carbides, sulfides, fluorides, and nitrides. We believe that the atomic level insights gained from these simulations will lay the groundwork for the rational design of tailored and functional interphases that are critical for the success of multivalent battery technology.

Agarwal, Garvit↗
DOE OSTI2022

Riverbed Temperature and 4D ERT Monitoring Reveals Heterogenous Horizontal and Vertical Groundwater-Surface Water Exchange Flows Under Dynamic Stage Conditions

Groundwater surface water exchange plays a critical role in physical, biological, and geochemical function of coastal and riverine systems. Observing exchange flow behavior in heterogeneous systems is a primary challenge, particularly when flows are governed by dynamic river stage or tidal variations. In this paper we demonstrate a novel application of time-lapse 3D electrical resistivity tomography and temperature monitoring where an array of thermistors installed beneath a riverbed double as resistivity electrodes. We use the array to monitor stage driven exchange flows over a 6-day period in a dynamic, stage-driven high order stream. We present a method for addressing the otherwise confounding effects of the moving river-surface boundary on the raw resistivity data, thereby enabling successful tomographic imaging. Temperature time-series at each thermistor location and time-lapse 3D images of changes in bulk electrical conductivity together provide a detailed description of exchange dynamics over a 10-meter by 45-meter section of the riverbed, to a depth of approximately 5 m. Results reveal highly variable flux behavior throughout the monitoring domain including both horizontal and vertical exchange flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Rewritable Surface‐Grafted Polymer Brushes with Dynamic Covalent Linkages

Abstract Surface grafting of polymer brushes drastically modifies surface properties, including wettability, compatibility, responsiveness, etc. A broad variety of functionalities can be introduced to the surface via different types of polymers. Bringing together properties of two or more types of polymer brushes to one surface opens up even more possibilities in brush‐modified materials. However, while it is generally feasible to introduce several chemical compositions along the brushes via copolymerization, it is challenging to vary the types of polymer brushes along a surface. Although previous studies have demonstrated binary brushes via orthogonal polymerization techniques or partial deactivation/regrafting, they commonly limit the number of polymer types to two. Here, we propose a strategy to introduce dynamic covalent diketoenamine linkages at the root of polymer brushes. The grafting density could be precisely tuned during surface functionalization. The surface‐anchored polymer brushes were cleaved by the addition of small molecule amines. New polymer brushes can be regrafted from the surface following refunctionalization of exposed sites. The maneuverability allows tuning of the types and densities of the polymer brushes, pointing the way to the preparation of a new generation of well‐defined brush‐modified materials with mixed grafts, with potential applications in the design of smart materials and surfaces.

Cui, Feichen↗
DOE OSTI2023

Chemical looping conversion of CH4/CO2 to syngas on 5wt.%Ni@Ce0.6Zr0.4O2 catalyst: Impact of dynamic accumulation of surface carbon and oxygen vacancies

Syngas, a combination of carbon monoxide (CO) and hydrogen (H2), is a precursor to many chemicals and fuels, contributing to the billion-dollar global hydrocarbon industry. Chemical looping reforming (CLR) of the greenhouse gases methane (CH4) and carbon dioxide (CO2) allows for energy-efficient production of syngas. 5wt.% nickel (Ni) on ceria-zirconia (5wt%Ni@Ce0.6Zr0.4O2) mixed metal oxide catalyst was investigated here to explore pathways for enhanced syngas production on sustainable earth-abundant transition metal supported catalysts. The role of reduction-oxidation (redox) state of the catalyst, and carbon formation on the catalyst surface during chemical looping is explored to drive superior reaction kinetics, conversions, selectivity, and syngas ratios (H2/CO). The bulk and surface structure of the catalyst, along with carbon deposition features were characterized by electron microscopy, X-Ray diffraction, and ex situ Raman spectroscopy. The dynamic evolution of catalysts under CLR reaction conditions and the intrinsic reaction mechanisms were probed with in situ Raman spectroscopy, in situ Fourier transform spectroscopy, dynamic oxygen storage capacity (DOSC), and the Temporal Analysis of Products (TAP) reactor studies. Intrinsic kinetics of syngas production via CLR was correlated to the redox state of catalyst and the participation of nickel-catalyzed multiwalled carbon nanotube (CNT) growth, allowing enhanced CLR reaction performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Role of heterogeneous surface wettability on dynamic immiscible displacement, capillary pressure, and relative permeability in a CO 2 -water-rock system

Numerical examination of CO 2 displacement in brine filled rock pores with varying wettability. Results show that as the CO 2 contact angles on the rock surface increased (intermediate water wetting conditions) negative entry pressures for CO 2 displacement were observed, indicating that when heterogeneous surface wettability is accounted for, greater mobility of CO 2 into rock pores is expected. In addition, simulated relative permeability CO 2 curves with heterogeneous intermediate water wetting conditions exhibited a greater amount of variation, indicative of a complex interrelation of geometric capillary restrictions and surface wettability.

58 GEOSCIENCES↗
DOE OSTI2020

Wireless Interrogation of High Temperature Surface Acoustic Wave Dynamic Strain Sensor

Dynamic strain sensing is necessary for high-temperature harsh-environment applications, including powerplants, oil wells, aerospace, and metal manufacturing. Monitoring dynamic strain is important for structural health monitoring and condition-based maintenance in order to guarantee safety, increase process efficiency, and reduce operation and maintenance costs. Sensing in high-temperature (HT), harsh-environments (HE) comes with challenges including mounting and packaging, sensor stability, and data acquisition and processing. Wireless sensor operation at HT is desirable because it reduces the complexity of the sensor connection, increases reliability, and reduces costs. Surface acoustic wave resonators (SAWRs) are compact, can operate wirelessly and battery-free, and have been shown to operate above 1000°C, making them a potential option for HT HE dynamic strain sensing. This paper presents wirelessly interrogated SAWR dynamic strain sensors operating around 288.8MHz at room temperature and tested up to 400°C. The SAWRs were calibrated with a high-temperature wired commercial strain gauge. The sensors were mounted onto a tapered-type Inconel constant stress beam and the assembly was tested inside a box furnace. The SAWR sensitivity to dynamic strain excitation at 25°C, 100°C, and 400°C was .439 μV/με, 0.363μV/με, and .136 μV/με, respectively. Furthermore, the experimental outcomes verified that inductive coupled wirelessly interrogated SAWRs can be successfully used for dynamic strain sensing up to 400°C.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Predicting lake surface water phosphorus dynamics using process-guided machine learning

Phosphorus (P) loading to lakes is degrading the quality and usability of water globally. Accurate predictions of lake P dynamics are needed to understand whole-ecosystem P budgets, as well as the consequences of changing lake P concentrations for water quality. However, complex biophysical processes within lakes, along with limited observational data, challenge our capacity to reproduce short-term lake dynamics needed for water quality predictions, as well as long-term dynamics needed to understand broad scale controls over lake P. Here we use an emerging paradigm in modeling, process-guided machine learning (PGML), to produce a phosphorus budget for Lake Mendota (Wisconsin, USA) and to accurately predict epilimnetic phosphorus over a time range of days to decades. In our implementation of PGML, which we term a Process-Guided Recurrent Neural Network (PGRNN), we combine a process-based model for lake P with a recurrent neural network, and then constrain the predictions with ecological principles. We test independently the process-based model, the recurrent neural network, and the PGRNN to evaluate the overall approach. The process-based model accounted for most of the observed pattern in lake P; however it missed the long-term trend in lake P and had the worst performance in predicting winter and summer P in surface waters. The root mean square error (RMSE) for the process-based model, the recurrent neural network, and the PGRNN was 33.0 μg P L -1 , 22.7 μg P L -1 , and 20.7 μg P L -1 , respectively. All models performed better during summer, with RMSE values for the three models (same order) equal to 14.3 μg P L -1 , 10.9 μg P L -1 , and 10.7 μg P L -1 . Although the PGRNN had only marginally better RMSE during summer, it had lower bias and reproduced long-term decreases in lake P missed by the other two models. For all seasons and all years, the recurrent neural network had better predictions than process alone, with root mean square error (RMSE) of 23.8 μg P L -1 and 28.0 μg P L -1 , respectively. The output of PGRNN indicated that new processes related to water temperature, thermal stratification, and long term changes in external loads are needed to improve the process model. By using ecological knowledge, as well as the information content of complex data, PGML shows promise as a technique for accurate prediction in messy, real-world ecological dynamics, while providing valuable information that can improve our understanding of process.

54 ENVIRONMENTAL SCIENCES↗
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↗