Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interface”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Irradiation-induced formation of a spinel phase at the FeCr/MgO interface

Oxide dispersion strengthened ferritic/martensitic alloys have attracted significant attention for their potential uses in future nuclear reactors and storage vessels, as the metal/oxide interfaces act as stable high-strength sinks for point defects while also dispersing helium. Here, in order to unravel the evolution and interplay of interface structure and chemistry upon irradiation in these types of materials, an atomically sharp FeCr/MgO interface was synthesized at 500 °C and separately annealed and irradiated with Ni 3+ ions at 500 °C. After annealing, a slight enrichment of Cr atoms was observed at the interface, but no other structural changes were found. However, under irradiation, sufficient Cr diffuses across the interface into the MgO to form a Cr-enriched transition layer that contains spinel precipitates. First-principles calculations indicate that it is energetically favorable to incorporate Cr, but not Fe, substitutionally into MgO. Furthermore, our results indicate that irradiation can be used to form new phases and complexions at interfaces, which may have different radiation tolerance than the pristine structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of three-dimensional Cu/Nb interfaces on strengthening and shear banding in nanoscale metallic multilayers

Manipulation of the atomic-scale structure of two-dimensional (2D) interfaces have been shown to provide nanocomposites with enhanced strength, deformability, and radiation damage resistance. In comparison with 2D interfaces, here we investigate the mechanical response of nanocomposites containing three-dimensional (3D) Cu/Nb interfaces consisting of a chemical/structural gradient separating pure Cu and Nb layers, through which the lattice mismatch between face-centered cubic Cu and body-centered cubic Nb is accommodated over a distance of several nanometers. Furthermore, it is demonstrated that 3D interfaces increase the yield and flow strength by 50% and 22%, respectively, over composites containing 2D interfaces at similar layer thicknesses. After 14% compressive strain, the onset of shear banding results in co-deformation of both Cu and Nb phases within and outside of the shear band. We conclude with a discussion of the role of interface structure in shear band formation and growth in 3D Cu/Nb.

3-dimensional interface↗
DOE OSTI2022

Characterization of the atomic-level structure of γ-alumina and (111) Pt/γ-alumina interfaces

The atomic-level structure of platinum/γ-alumina interfaces is characterized in a model system of dense γ-alumina embedded with faceted Pt NPs produced by implantation of platinum ions into sapphire followed by thermal annealing in air at 800 °C. Aberration-corrected scanning transmission electron microscopy (STEM) was used to collect atomic-resolution images, which are compared to STEM image simulations of two experimentally-based bulk models of γ-alumina by Smrčok et al. and Zhou and Snyder. A density functional theory (DFT) based model of (111) interfaces with different chemical terminations (O, Al 1 , Al 2 ) of the γ-alumina developed by Oware Sarfo et al. is also compared to experimental STEM data from Pt/γ-alumina interfaces. Further, the Smrčok γ-alumina model provides a better fit than the Zhou structure to the bulk of the γ-alumina. The oxygen-terminated Oware Sarfo model best fits the experimental data and is a very good model close to the interface. However, the fit of the interface model to the experimental data is poorer beyond the third atomic layer in the γ-alumina. This is attributed to compromises required in the design of the model to limit the cell size and computational time for DFT calculations. Understanding the accuracy and limits of the structural models of γ-alumina and Pt/γ-alumina interfaces is important to further the understanding of the structure/property relationships in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-material hydrodynamics with algebraic sharp interface capturing

A finite volume method for Eulerian multi-material hydrodynamics with sharp interface capturing is presented here. The pressure-temperature non-equilibrium multi-material system with finite-rate pressure relaxation in mixed-cells is considered here. This pressure closure facilitates material-property-dependent pressure relaxation, rather than instantaneous pressure equilibration, which in turn allows the use of unsplit high-order time-integrators. A modified tangent of hyperbola for interface capturing (THINC) method is used to reconstruct multi-material (> 2) interfaces, on three-dimensional unstructured meshes. A simple modification which extends the THINC reconstruction to interfaces between more than two materials is proposed. It is demonstrated that the modified THINC can capture multi-material interfaces within 2–4 tetrahedral cells. Since no geometric reconstructions are required by the THINC method, the presented multi-material method is algorithmically simple, and computationally efficient. Consistent reconstructions of conserved quantities at material interfaces ensure that conservation and closure laws are satisfied at the discrete level. Through a suite of test problems solved on unstructured meshes, it is demonstrated that the presented method is a promising candidate for accurate and efficient multi-material hydrodynamics computations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Influence of processing on secondary phase formation and microstructural evolution at U-10Mo alloy and Zr interlayer interfaces

Thermo-mechanical processing of uranium-10 wt. % molybdenum (U-10Mo) fuel plates leads to microstructure changes at the U-10Mo/Zr interfaces. Secondary phases formed at this interface are particularly important to interfacial bond strength, process optimization, and maintaining structural integrity of the U-10Mo fuel plates during irradiation. In this work, we determined the phases and phase transformation products occurring at the interface of the U-10Mo fuel and Zr interlayer when the fuel plate is subjected to short and long hot isostatic pressure times. Interfacial morphology, structure and composition of phases formed, and relative hardness across the U-10Mo/Zr interfaces were studied using a multi-length scale, multi-modal characterization approach involving electron microscopy, atom probe tomography, and atomic force microscopy. Here, results highlight that the extent of phase transformations, secondary phase formation, and hardness variability across interfaces can be controlled by modifying processing parameters. Phase diagram construction and thermodynamic calculations were performed using the Thermocalc software to identify expected phases formed at interfaces during the maximum hold temperature of 560 °C experienced during HIP.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In conclusion, these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II

Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

An axial temperature profile curvature criterion for the engineering of convex crystal growth interfaces in Bridgman systems

A unifying idea is presented for the engineering of convex melt-solid interface shapes in Bridgman crystal growth systems. Previous approaches to interface control are discussed with particular attention paid to the idea of a “booster” heater. Proceeding from the idea that a booster heater promotes a converging heat flux geometry and from the energy conservation equation, we show that a convex interface shape will naturally result when the interface is located in regions of the furnace where the axial thermal profile exhibits negative curvature, i.e., where d 2 T/dz 2 < 0. This criterion is effective in explaining prior literature results on interface control and promising for the evaluation of new furnace designs. Here we posit that the negative curvature criterion may be applicable to the characterization of growth systems via temperature measurements in an empty furnace, providing insight about the potential for achieving a convex interface shape, without growing a crystal or conducting simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Conversion of ethanol to 1,3–butadiene over Ag–ZrO 2 /SiO 2 catalysts: The role of surface interfaces

A series of Ag–ZrO 2 /SiO 2 catalysts with different metal–support interfaces were synthesized in an effort to elucidate the roles of specific interfaces in controlling the ethanol to 1,3–butadiene conversion and selectivity. According to the results of detailed characterizations (e.g. CO/pyridine–DRIFTS, XPS, TEM, NH 3 –TPD, and 1 H MAS NMR), it was found that the Ag–O–Si interfaces significantly enhanced the dehydrogenation of ethanol while the presence of ZrO 2 improved the interaction between Ag and ZrO 2 /SiO 2 , creating more Ag δ+ active sites. In this research, the high dispersion of ZrO 2 on SiO 2 generated abundant Zr–O–Si interfaces with medium and weak Lewis acidity, promoting the condensation of acetaldehyde to crotonaldehyde. These Zr–O–Si interfaces in close interaction with Ag δ+ species played a critical role in the enhanced H transfer during the MPV reduction of crotonaldehyde to crotyl alcohol. The synergies among the interfaces resulted in retarded ethanol dehydration reactivity, balanced ethanol dehydrogenation and condensation reactions, and a subsequent high 1,3–butadiene yield.

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

From coherent to incoherent mismatched interfaces. A generalized continuum formulation of surface stresses

The equilibrium of coherent and incoherent mismatched interfaces is reformulated in the context of continuum mechanics based on the Gibbs dividing surface concept. Two surface stresses are introduced: a coherent surface stress and an incoherent surface stress, as well as a transverse excess strain. Additionally, the coherent surface stress and the transverse excess strain represent the thermodynamic driving forces of stretching the interface while the incoherent surface stress represents the driving force of stretching one crystal while holding the other fixed and thereby altering the structure of the interface. These three quantities fully characterize the elastic behavior of coherent and incoherent interfaces as a function of the in-plane strain, the transverse stress and the mismatch strain. The isotropic case is developed in detail and particular attention is paid to the case of interfacial thermo-elasticity. This exercise provides an insight on the physical significance of the interfacial elastic constants introduced in the formulation and illustrates the obvious coupling between the interface structure and its associated thermodynamics quantities. Finally, an example based on atomistic simulations of Cu/Cu 2 O interfaces is given to demonstrate the relevance of the generalized interfacial formulation and to emphasize the dependence of the interfacial thermodynamic quantities on the incoherency strain with an actual material system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗
DOE OSTI2016

Structural ordering at solid-liquid interfaces in Al-Sm system: A molecular-dynamics study

The structural ordering at solid-liquid interfaces far from equilibrium is studied with molecular dynamics simulations for the Al-Sm system. Using the van-Hove self-correlation function as the criterion to identify attachment/detachment events that occur at the interface, we are able to determine the time-dependent interface position, and characterize the detailed interfacial structure ordering surrounding the attached atoms. For the interface between an undercooled Al90Sm10 liquid and a metastable cubic structure, the solid induces the crystalline order of the cubic phase in the liquid layers, promoting the continuous growth of the crystal phase. When the same liquid is put in contact with f.c.c. Al, Sm from the liquid can still attach to the solid interface despite its insolubility in the Al lattice. Non-f.c.c. order is revealed surrounding the attached Sm atoms. Lastly, we show that the local structure ordering at interface is highly correlated to solid packing and liquid ordering.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Modular multi-interface nanocrystals for enhanced ethanol oxidation electrocatalysis

Electrochemical processes that utilize biomass-derived ethanol as a source of electrons and protons offer a sustainable energy strategy, yet their practical implementation is limited by sluggish ethanol oxidation reaction (EOR) kinetics and catalyst poisoning. Here, in this study, we report a modular multi-interface nanocrystal catalyst comprising core/shell Co 2 P/Pd and Pd-Au heterostructured interfaces that exhibit complementary functions for the enhanced EOR catalysis. The Co 2 P/Pd interface boosts Pd atom utilization and lowers the kinetic barriers for ethanol-to-acetate conversion, while the Pd-Au interface effectively alleviates CO poisoning caused by C–C bond cleavage of ethanol. In-depth analyses using in situ attenuated total reflectance-surface-enhanced infrared absorption spectroscopy, differential electrochemical mass spectrometry, and density functional theory calculations elucidate the mechanistic roles of these interfaces. The optimized Co 2 P/Pd-Au 0.08 nanorods achieve an excellent mass activity, underscoring the potential of modular, multi-interface nanocrystals for advancing EOR catalysis and offering a generalizable strategy for broader catalytic innovations.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Interface reduction for Hurty/Craig-Bampton substructured models: Review and improvements

The Hurty/Craig-Bampton method in structural dynamics represents the interior dynamics of each subcomponent in a substructured system with a truncated set of normal modes and retains all of the physical degrees of freedom at the substructure interfaces. This makes the assembly of substructures into a reduced-order system model relatively simple, but means that the reduced-order assembly will have as many interface degrees of freedom as the full model. When the full-model mesh is highly refined, and/or when the system is divided into many subcomponents, this can lead to an unacceptably large system of equations of motion. To overcome this, interface reduction methods aim to reduce the size of the Hurty/Craig-Bampton model by reducing the number of interface degrees of freedom. This research presents a survey of interface reduction methods for Hurty/Craig-Bampton models, and proposes improvements and generalizations to some of the methods. Some of these interface reductions operate on the assembled system-level matrices while others perform reduction locally by considering the uncoupled substructures. The advantages and disadvantages of these methods are highlighted and assessed through comparisons of results obtained from a variety of representative linear FE models.

42 ENGINEERING↗
DOE OSTI2021

Tailoring Interfaces in Solid-State Batteries Using Interfacial Thermochemistry and Band Alignment

Solid-state lithium ion batteries have enhanced thermal stability compared to batteries using liquid electrolytes. Although their practical implementation is limited by undesired side reactions between the electrode and electrolyte, different modification strategies have been proposed to stabilize these reactive interfaces. These approaches have been primarily based upon bulk materials properties, however, which may not necessarily be representative of the chemistry at the solid/solid interface. In this work, first principles calculations and X-ray scattering experiments are used to elucidate molecular-level reactivity and develop design principles for tailored interfaces based on surface and interfacial properties. Because of its well-known instability, the interface between the Li metal anode and the lithium lanthanum titanate (LLTO) solid electrolyte is applied as a model system. Density functional theory calculations are used to describe bulk, surface, and interfacial thermochemistry of the Li/LLTO system, and interfacial reconstruction is probed using ab initio molecular dynamics. These simulations of the Li/LLTO interface demonstrate facile decomposition concomitant with reduction of Ti 4+ ions. Based on further insights from computational analysis of the surface band edge positions, La 2 O 3 is proposed as an interlayer coating and is shown to provide an energetic barrier for interfacial charge transfer and decomposition reactions from X-ray reflectivity measurements and theoretical calculations. The findings in this work suggest that design strategies based on surface and interfacial properties can be used to kinetically stabilize electrode/electrolyte interfaces in solid-state batteries.

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

Self-Assembly of Silica Nanoparticles at Water–Hydrocarbon Interfaces: Insights from In Operando Small-Angle X-ray Scattering Measurements and Molecular Dynamics Simulations

To achieve tunable controls on the interactions of siliceous materials in subsurface environments bearing water and hydrocarbons, it is essential to determine the influence of hydrocarbon–water interfaces on the self-assembly of the silica particles. The hydrophilic silica particle have a tendency to aggregate on the water front of the interface. Self-assembly of silica nanoparticles proceeds via the migration of these nanoparticles to the water–hydrocarbon interface followed by aggregation at the interface. Fractal-like morphologies of assembled silica nanoparticles at water–hydrocarbon interfaces are observed. Rapid assembly of the hydrophilic silica nanoparticles at water–hydrocarbon interfaces corresponds to an overall reduction in the surface tension of water–toluene and water–heptane systems. These studies demonstrate the silica aggregation that is the precursor of silica polymerization, or nucleation and growth is influenced by the presence of hydrocarbons in subsurface geologic environments. These insights were derived from in operando ultrasmall and small-angle X-ray scattering (USAXS/SAXS) measurements, cryo-scanning electron microscopy (Cryo-SEM) imaging, and classical molecular dynamics simulations. As a result, these studies are intended to inform current and future efforts aimed at tuning silica reactivity in subsurface geologic environments to enhance permeability, the design and use of silica-based proppants to enhance fractures, and the development of effective strategies to control silica-based scaling behavior in subsurface reservoirs.

42 ENGINEERING↗
DOE OSTI2025

Reorganization of Water at Aqueous Aluminum Chloride (AlCl 3 ) Interfaces: Vibrational Sum Frequency Generation and Molecular Dynamics Simulations

AlCl 3 hydration states and complexation are not well understood both in solutions and at the air–aqueous interface despite their potential significance in natural waters and their industrial and energy-related applications. Here, we investigated Al 3+ and Cl – ion behaviors in an AlCl 3 aqueous bulk solution and at the air–aqueous interface using interface-selective vibrational sum frequency generation (SFG), Raman and infrared spectroscopies, molecular dynamics (MD) simulation, as well as molecular-informed reduced modeling. Our reduced modeling reveals relatively long-range effects for Al 3+ as compared to monovalent ions such as Na + indicating that the interfacial depth of trivalent ions can be significantly larger than that of monovalent ions at the air–water interface. MD simulations reveal interfacial stratification and multiple layering of the ions. Compression of the Al 3+ and Cl – distributions with increasing concentrations from 0.5 to 2.5 m is also observed in the subsurface regions. Significant SSP- and PPP-polarized SFG OH spectral intensity increases are observed from 0.5 to 1.5 m and 0.5 to 2.5 m, respectively, indicative of interfacial depth increases and a change in average orientation above 1.5 m. Extensive evaluation of SFG spectra, Fresnel-corrected using several approaches, shows the same trends. The nonmonotonic trend points to a changing structure in surface and subsurface water orientation and hydrogen bonding environment generally consistent with the MD simulation of stratification and water orientation changes. Furthermore, solvent-shared ion pairing is implicated with MD simulation radial distribution analysis and consistent with infrared spectral identification of the hexaaqua aluminum ion in the solution phase. Spectral evidence of a strong Al 3+ hydration shell and the acidic behavior of the Al 3+ ions is obvious in the Raman and infrared spectra of the bulk solution. In conclusion, we show that the MD dipole potential is directly related to the MD second-order susceptibility of the interface, χ SFG–MD (2) , both of which correlate up to ∼35 Å with the spectral observations of increasing and then saturating intensities, suggesting that both ion stratification and interfacial depth determine the water orientations at an air–water interface of 1-3 electrolyte solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗