Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atom trapping”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Losses of CO and CO 2 upon collision-activated dissociation of substituted 2-methoxyphenoxides after methyl radical loss

Some deprotonated, substituted 2-methoxyphenols fragment by methyl radical loss followed by competing losses of CO and CO 2 upon collisional activation in linear quadrupole ion trap mass spectrometers. These reactions were examined experimentally and computationally in order to determine their mechanisms. For deprotonated vanillin, the CO loss was found to involve ring contraction, with the highest free energy barrier of 58.0 kcal mol -1 (for the syn rotamer; here, syn refers to the relationship between the aldehyde oxygen atom and the oxygen atom of the methoxy group) for the entire process. The atoms lost in this fragmentation are the oxygen atom that was bound to the first eliminated methyl group and the aromatic carbon atom bound to this oxygen. Examination of carbon-13 labeled vanillin supports these assignments. Examination of several model compounds revealed that this reaction requires the presence of an electron-withdrawing substituent in the para-position relative to the phenol moiety. In contrast, the CO 2 loss from deprotonated vanillin occurs via ring opening followed by re-cyclization and then ring contraction, leading to the loss of CO 2 in a process wherein the highest free energy barrier is 84.5 kcal mol -1 (for the syn isomer). The atoms lost in this fragmentation are the carbon and oxygen atoms from the phenoxide group and the oxygen atom that was bound to the first eliminated methyl group, which is supported by examination of carbon-13 labeled vanillin. Despite the higher total free energy requirement, the CO 2 loss is competitive with CO loss, possibly due to favorable entropy and low activation energy (free energy barrier of 48.9 kcal mol -1 ) for the first reaction step (the analogous value for CO loss is 58.0 kcal mol -1 ). The extent of CO 2 loss is strongly affected by substituents – it either is not observed or is very slow for the other compounds studied here. For example, it is substantially less favorable than CO loss (highest free energy barrier 60.9 kcal mol -1 ) for deprotonated acetovanillone for which the first reaction step for CO 2 loss has a substantially greater free energy barrier (60.6 kcal mol-1) than for vanillin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Assessment of deep levels with selenium concentration in Cd 1–x Zn x Te 1–y Se y room temperature detector materials

Incorporation of Se into Cd 1–x Zn x Te (CZT) to form the quaternary compound semiconductor Cd 1–x Zn x Te 1–y Se y (CZTS) has proven to be an effective solution for compensating the major flaws associated with CZT, including poor homogeneity and high concentrations of electronically active deep levels that limit the performance of CZT detectors. In order to investigate how deep levels are affected by the Se concentration in CZTS, we performed photoinduced current transient spectroscopy (PICTS) measurements on CZTS crystals grown by the traveling heater method (THM) with 10% atomic Zn and varying atomic percentage of Se from 1.5% to 7.0%. The PICTS scans for up to 4% Se showed an exponential reduction in the capture cross section of deep levels associated with Te secondary phases in conjunction with an increase in a deep level positioned near the mid-gap, which initially increases the electron trapping time before degrading again at higher Se concentrations. The PICTS peaks present in 7% Se were anomalous relative to the other crystals and are expected to originate from transition metal impurities found in the lower-purity CdSe precursor material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of divertor material on neutral recycling and discharge fueling in DIII-D

We report that experiments with the lower divertor of DIII-D during the Metal Rings Campaign (MRC) show that the fraction F of atomic D in the total recycling flux is material-dependent and varies through the ELM cycle, which may affect divertor fueling. Between ELMs, F C ~ 10% and FW ~ 40%, consistent with expectations if all atomic recycling is due to reflections. During ELMs, F C increases to 50% and F W to 60%. In contrast, the total D recycling coefficient including atoms and molecules R stays close to unity near the strike point where the surface is saturated with D. During ELMs, R can deviate from unity, increasing during high energy ELM-ion deposition (net D release) and decreasing at the end of the ELM which leads to ability of the target to trap the ELM-deposited D. The increase of R > 1 in response to an increase in ion impact energy Ei has been studied with small divertor target samples using Divertor Materials Evaluation System (DiMES). An electrostatic bias was applied to DiMES to change E i by 90 eV. On all studied materials including C, Mo, uncoated and W-coated TZM (>99% Mo, Ti, and Zr alloy), W, and W fuzz, an increase of E i transiently increased the D yield (and R) by ~10%. On C there was also an increase in the molecular D 2 yield, probably due to ion-induced D 2 desorption. Despite the measured increase in F on W compared to C, attached H-mode shots with OSP on W during MRC did not demonstrate a higher pedestal density. About 8% increase in the edge density could be seen only in attached L-mode scenarios. The difference can be explained by higher D trapping in the divertor and lower divertor fueling efficiency in H- versus L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Radiation damage of germanium detectors

Energetic particles can produce interstitial-vacancy pairs in a crystal by knocking the atoms from their normal positions. Detectors are unique among semiconductor devices in depending on very low concentrations of electrically active impurities, and also on efficient transport of holes and electrons over relatively large distances. Because the dense regions of damage produced by energetic particles may result in donors and/or acceptors, and also provide trapping sites for holes and electrons, detectors are very sensitive to radiation damage. In addition to these effects occurring within the detector, radiation may also change the characteristics of the exposed surfaces causing unpredictable effects on the detector leakage current. Radiation-induced surface degradation has rarely, if ever, been observed for germanium detectors. The possibility of minimizing hole trapping in charge collection by the use of a high-purity germanium coaxial detector configured with the p (+) contact on the coaxial periphery is discussed.

Pehl, R. H.↗
DOE OSTI2021

Spatial coherence of room-temperature monolayer WSe 2 exciton-polaritons in a trap

The emergence of spatial and temporal coherence of light emitted from solid-state systems is a fundamental phenomenon intrinsically aligned with the control of light-matter coupling. It is canonical for laser oscillation, emerges in the superradiance of collective emitters, and has been investigated in bosonic condensates of thermalized light, as well as exciton-polaritons. Our room temperature experiments show the strong light-matter coupling between microcavity photons and excitons in atomically thin WSe 2 . We evidence the density-dependent expansion of spatial and temporal coherence of the emitted light from the spatially confined system ground-state, which is accompanied by a threshold-like response of the emitted light intensity. Additionally, valley-physics is manifested in the presence of an external magnetic field, which allows us to manipulate K and K’ polaritons via the valley-Zeeman-effect. Our findings validate the potential of atomically thin crystals as versatile components of coherent light-sources, and in valleytronic applications at room temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Exploring the Potential of Residual Impurities in Germanium Detectors for MeV-Scale Dark Matter Detection

Abstract The direct detection of MeV-scale dark matter (DM) particles hinges on achieving an exceptionally low-energy detection threshold. Germanium (Ge) detectors, meticulously tailored with precise impurity compositions, hold the potential to enhance sensitivity to energy levels below the sub-electronvolt (sub-eV) range. This study explores the behavior of residual impurities inherent to Ge detectors at helium temperatures, unveiling a captivating freeze-out phenomenon leading to the formation of excited localized states known as dipole states. Using compelling evidence from relative capacitance measurements obtained from two detectors, we elucidate the transition of impurity atoms from free charge states to these dipole states as the temperature drops from 11 to 6.5 K. Our investigation comprehensively covers the intricate formation of these dipole states in bothn-type andp-type impurities. Furthermore, we shed light on the electric field generated by these dipole states, revealing their ability to trap charges and facilitate the creation of cluster dipole states. Confirming findings from previous measurements, we establish that these excited dipole states exhibit a binding energy of less than 10 meV, offering an exceptionally low detection threshold for MeV-scale DM. Building upon this concept, we propose the development of a 1-kg Ge detector with internal charge amplification—an innovative approach poised to surpass electrical noise and enable the detection of MeV-scale DM with unprecedented sensitivity.

Physics↗
DOE OSTI2026

Measuring X-Ray Emission Line Shapes in Neutral Species for XRISM Calibration

Space X-ray spectrometers such as the Resolve instrument on XRISM require precise calibration in order to interpret the spectra of astrophysical objects. Key components of the calibration are the energy scale and the core line spread function, both of which vary with photon energy. A major issue in the calibration of high-resolution spectrometers is locating good calibrators with well-known and stable intrinsic line shapes. Neutral fluorescence is widely used, but inner-shell transitions in neutral atoms often exhibit complex, poorly documented line shapes that vary with excitation conditions. Here, in this study, we present empirical measurements of K-shell transitions in neutral O and F below 1 keV using an engineering model XRISM calorimeter array, an electron bombardment modulated X-ray source, and an electron beam ion trap (EBIT) to provide a precise energy reference. In addition, we report measurements of the Mo Lα complex with the transition-edge microcalorimeter spectrometer (TEMS), which reveal strong satellite structure and sensitivity of the line shape to the incident exciting spectrum. Together, these results demonstrate the need for empirical line-shape models, highlight the nonstationary nature of neutral fluorescence features, and define a path toward developing transfer standards for XRISM and future precision instruments such as Athena/X-IFU.

Astronomy and AstroPhysics↗
DOE OSTI2024

Refining Absorber Shroud Geometry to Maximize Power Output and Reduce Power Peaking in ATF Test Train

In the wake of the Fukushima Daiichi nuclear power plant accident in 2011, the accident tolerant fuels (ATF) program was initiated to enhance the safety of light-water reactor fuels, placing significant emphasis on cladding. A crucial step for the broad implementation of ATF in commercial reactors involves irradiation testing of the fuel designs. ATF-2D, the latest experiment in the ATF series, is slated to undergo irradiation in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). ATF-2D is a joint effort of the INL with industry partners General Electric Global Research; Framatome; General Atomics; the Japan Atomic Energy Agency; Hitachi-GE Nuclear Energy, Ltd; Global Nuclear Fuel-Japan Co., Ltd; Nippon Nuclear Fuel Development Co., Ltd; and Mitsubishi Heavy Industries, Ltd. The ATF-2D test train design consists of four tiers, each housing six rodlets. The device will be inserted in Loop 2A within the central flux trap of the ATR, and is anticipated to undergo irradiation throughout three 60-day cycles. Typical pressurized water reactor conditions will be emulated during the irradiation. The objective of the work presented here is to dimension neutron-absorbing hafnium (Hf) components surrounding the fuel rodlets, such that the axial power profile is flattened, while simultaneously ensuring that the total fission power output of the entire test train remains below 200 kW.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Refining Absorber Shroud Geometry to Maximize Power Output and Reduce Power Peaking in ATF Test Train

In the wake of the Fukushima Daiichi nuclear power plant accident in 2011, the accident tolerant fuels (ATF) program was initiated to enhance the safety of light-water reactor fuels, placing significant emphasis on cladding. A crucial step for the broad implementation of ATF in commercial reactors involves irradiation testing of the fuel designs. ATF-2D, the latest experiment in the ATF series, is slated to undergo irradiation in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). ATF-2D is a joint effort of the INL with industry partners General Electric Global Research; Framatome; General Atomics; the Japan Atomic Energy Agency; Hitachi-GE Nuclear Energy, Ltd; Global Nuclear Fuel-Japan Co., Ltd; Nippon Nuclear Fuel Development Co., Ltd; and Mitsubishi Heavy Industries, Ltd. The ATF-2D test train design consists of four tiers, each housing six rodlets. The device will be inserted in Loop 2A within the central flux trap of the ATR, and is anticipated to undergo irradiation throughout three 60-day cycles. Typical pressurized water reactor conditions will be emulated during the irradiation. The objective of the work presented here is to dimension neutron-absorbing hafnium (Hf) components surrounding the fuel rodlets, such that the axial power profile is flattened, while simultaneously ensuring that the total fission power output of the entire test train remains below 200 kW.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Assessment of the Second-Ionization Potential of Lawrencium: Investigating the End of the Actinide Series with a One-Atom-at-a-Time Gas-Phase Ion Chemistry Technique

Experiments were performed at the Lawrence Berkeley National Laboratory 88-Inch Cyclotron facility to investigate the electron-transfer reduction reaction of dipositive Lr ( Z = 103) with O 2 gas. Ions of 255 Lr were produced in the fusion-evaporation reaction 209 Bi( 48 Ca,2n) 255 Lr and were studied with a novel gas-phase ion chemistry technique. The produced 255 Lr 2+ ions were trapped and O 2 gas was introduced, such that the charge-exchange reaction to reduce 255 Lr 2+ to 255 Lr 1+ was observed and the reaction rate constant was determined to be k = 1.5(7) × 10 -10 cm 3 /mol/s. The observation that this reaction proceeds establishes the lower limit on the second ionization potential of Lr to be 13.3(3) eV. This gives further support that the actinide series terminates with Lr. Additionally, this result can be used to better interpret the situation concerning the placement of Lu and Lr on the periodic table within the current framework of the actinide hypothesis. The success of this experimental approach now identifies unique opportunities for future gas-phase reaction studies on actinide and super heavy elements.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

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

Defect Properties, Anion Ordering, and Photochromic Mechanism in Yttrium Oxyhydride

Yttrium oxyhydride (YHO) undergoes a reversible photochromic transition when exposed to ultraviolet light. However, the mechanism for this transformation is not fully understood, and the structure and precise chemical composition of YHO remain under debate. Here, we use first-principles density functional theory calculations with a hybrid functional to study the structure, chemical stability, and point defect properties of YHO. As experiments have shown, we find that YHO prefers a cubic structure, with H and O anions present in equal concentrations and located on tetrahedral sites. Stoichiometric and ordered YHO is chemically stable, but it has a wide band gap of 5.01 eV, considerably larger than that measured in experiments (2.4–3.8 eV). On the other hand, Y4H10O has a smaller band gap of 2.97 eV and also has a region of chemical stability; thus, the actual material may include some fraction of this H-rich structure. The defect chemistry of YHO is dominated by anionic antisite species (H O and O H ), with hydrogen interstitials (H i ) and vacancies (V H ) also present in reasonably high concentrations. We show that antisite disorder lowers the band gap relative to the perfectly ordered structure, bringing the magnitude of the gap into closer agreement with experiment. Based on our calculations of defect migration and the positions of defect states relative to the band edges, we link the onset of photochromic behavior to the reaction H O – → V O 0 + H i – , which follows photoexcitation of a H O + defect. H i – can subsequently migrate away and be trapped by additional H O + defects, contributing to the persistence of the reaction, while the resultant oxygen vacancy, V O 0 , introduces an occupied defect state that leads to optical absorption at visible wavelengths. Our results can explain reported discrepancies between experimental and computational results for YHO, and they allow us to propose specific atomic-scale processes that can lead to photochromism. In conclusion, understanding these mechanisms is key for unlocking YHO’s application in devices ranging from smart windows and optoelectronics to electrochemical synapses for neural networks.

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

High-precision Penning trap mass measurements of neutron-rich chlorine isotopes at the 𝑁 = 28 shell closure

Although it is known that the N = 28 spherical shell closure erodes, the strength of the closure with decreasing proton number Z < 20 is an open question in nuclear structure. Here, in this region of interest, high-precision mass measurements of neutron-rich 43-45 Cl isotopes were performed at the Low Energy Beam and Ion Trap (LEBIT) when coupled to the National Superconducting Cyclotron Lab. The resulting mass excesses (MEs) are ME(⁴³Cl) = -24114.4(1.7) keV, ME(⁴⁴Cl) = -20450.8(10.6) keV, and ME(⁴⁵Cl) = -18240.1(3.7) keV, and improve the uncertainty of these masses by up to a factor of ∼40 compared to the previous values reported in the 2020 Atomic Mass Evaluation. Comparison to ab initio calculations using the Valence-Space In-Medium Similarity Renormalization Group (VS-IMSRG) shows good agreement up to and including the closure.

Erington, H. [Michigan State Univ., East Lansing, ↗
NASA NTRS2005

Protein Crystallization

Nucleation, growth and perfection of protein crystals will be overviewed along with crystal mechanical properties. The knowledge is based on experiments using optical and force crystals behave similar to inorganic crystals, though with a difference in orders of magnitude in growing parameters. For example, the low incorporation rate of large biomolecules requires up to 100 times larger supersaturation to grow protein, rather than inorganic crystals. Nucleation is often poorly reproducible, partly because of turbulence accompanying the mixing of precipitant with protein solution. Light scattering reveals fluctuations of molecular cluster size, its growth, surface energies and increased clustering as protein ages. Growth most often occurs layer-by-layer resulting in faceted crystals. New molecular layer on crystal face is terminated by a step where molecular incorporation occurs. Quantitative data on the incorporation rate will be discussed. Rounded crystals with molecularly disordered interfaces will be explained. Defects in crystals compromise the x-ray diffraction resolution crucially needed to find the 3D atomic structure of biomolecules. The defects are immobile so that birth defects stay forever. All lattice defects known for inorganics are revealed in protein crystals. Contribution of molecular conformations to lattice disorder is important, but not studied. This contribution may be enhanced by stress field from other defects. Homologous impurities (e.g., dimers, acetylated molecules) are trapped more willingly by a growing crystal than foreign protein impurities. The trapped impurities induce internal stress eliminated in crystals exceeding a critical size (part of mni for ferritin, lysozyme). Lesser impurities are trapped from stagnant, as compared to the flowing, solution. Freezing may induce much more defects unless quickly amorphysizing intracrystalline water.

Chernov, Alexander A.↗
DOE OSTI2026

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗
NASA NTRS2002

Assembly of Colloidal Materials Using Bioadhesive Interactions

We have pursued the use of biological crosslinking molecules of several types to make colloidal materials at relatively low volume fraction of colloidal particles. The objective is to make binary alloys of colloidal particles, made of two different colloidal particles coated with complementary biological lock-and-key binding molecules, which assemble due to the biological specificity. The long-term goal is to use low affinity lock-and-key biological interactions, so that the can anneal to form crystalline states. We have used a variety of different surface chemistries in order to make colloidal materials. Our first system involved using selectin-carbohydrate (sialyl-Lewis) interactions; this chemistry is derived from immune system. This chemical interaction is of relatively low affinity, with timescales for dissociation of several seconds. Furthermore, the adhesion mediated by these molecules can be reversed by the chelation of calcium atoms; thus assembled structures can be disassembled reversibly. Our second system employed avidin-biotin chemistry. This well-studied system is of high affinity, and is generally irreversible on a laboratory time-scale. Thus, we would expect selectin-carbohydrate interactions at high molecular density and avidin-biotin interactions to give kinetically-trapped structures; however, at low densities, we would expect significant differences in the structure and dynamics of the two materials, owing to their very different release rates. We have also begun to use a third chemistry - DNA hybridization. By attaching single stranded DNA oligonucleotide chains to beads, we can drive the assembly of colloidal materials by hybridization of complementary DNA chains. It is well known that DNA adenosine-thymine (A-T) and guanine-cytosine (G-C) bases hybridize pairwise with a Gibbs free energy change of 1.7 kcal/mol per base; thus, the energy of the assembly can be modulated by altering the number of complementary bases in the DNA chains. Using these different crosslinking molecules, we have assembled colloidal materials from different-sized colloidal particles, A and B. In the first sets of experiment, we used high densities of adhesion molecules, and 0.96 micron (A) and 5.5 micron (B) diameter particles. The high density of adhesion molecules means that the structures are kinetically trapped in nonequilibrium configurations. The structure of the suspension can be varied by changing the number ratio of the two types of colloidal particles, NA and NB, where A is the smaller particle. With carbohydrate-selectin or avidin-biotin interactions, large NA/NB leads to the formation of colloidal micelles, with the large center B particle surrounded by many smaller A particles. As the ratio NA/NB decreases, the structures become more extended, approaching the formation of macro-Rouse polymers - extended linear chains where A beads are connected with intervening small B linkers.

Hammer, Daniel A.↗
DOE OSTI2020

Electron Trapping and Ion Leaching at the Li-Modified Quartz–Water Interface

Electrons and ions at metal oxide-water interfaces have a critical role in many phenomena and applications, thus making their properties of considerable interest. We here investigate the behavior of a lithium impurity at the water interface with Li-doped SiO 2 using hybrid density functional based ab initio molecular dynamics simulations. We find that the excess electron donated by the lithium dopant localizes on a surface silicon atom and is partially solvated at the aqueous interface. While the position of the excess electron does not change during our ~60 ps simulation, the lithium cation diffuses from its initial position in the SiO 2 subsurface toward the interface and eventually leaches out of the oxide surface and becomes solvated by four water molecules, forming an aqueous Li+- electron complex at the interface. The degree of interaction between the localized electron and the lithium ion controls not only the energy level of the excess electron but also the extent of structural distortion and hydrogen-bonding at the interface. The results of our study can be relevant for understanding the initial stages of glass corrosion, where alkali ions are leached from the surface region by interactions with water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗