Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interface 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

On the Structure–Property Relationship of Semi‐Coherent FeCr 2 O 4 /Cr 2 O 3 Spinel/Corundum Interfaces

Abstract Oxide heterointerfaces are extremely common in both natural and artificial composite structures, including corroded structural materials. Often, key properties such as segregation and atomic transport are dictated by the structure of these interfaces. However, despite this critical link, very few heterointerfaces have been studied in any detail at the atomic scale. Here, one important oxide heterointerface is examined, between spinel and corundum, using the chemical system FeCr 2 O 4 /Cr 2 O 3 as a representative and technologically important case. Using atomistic simulation techniques, it is found that the structure, particularly the local chemistry, of the interface depends on the crystal chemistry at the interface. This atomic and chemical structure further impacts important properties such as defect segregation and mass transport. It is found that defects can nucleate at some regions of these interfaces and migrate back and forth across the corundum layer, suggesting high atomic mobility that may be important for the evolution of spinel/corundum composite structures in extreme conditions.

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 OSTI2024

Electrolyte–Electrocatalyst Interfacial Effects of Polymeric Materials for Tandem CO 2 Capture and Conversion Elucidated Using In Situ Electrochemical AFM

Integrating CO 2 capture and electrochemical conversion has been proposed as a strategy to reduce the net energy required for CO 2 regeneration in traditional CO 2 capture and conversion schemes and can be coupled with carbon-free renewable electricity. Polyethylenimine (PEI)-based materials have been previously studied as CO 2 capture materials and can be integrated in these reactive capture processes. PEI-based electrolytes have been found to significantly increase the CO 2 loading, and impact selectivity and rate of product formation when compared to the conventional aqueous electrolytes. However, the influence of these materials at the catalyst–electrode interface is currently not well understood. In this study, PEI-based electrolytes were prepared and their impact on the morphology of a silver electrode performing electrochemical CO 2 reduction (CO 2 R) was studied using in situ electrochemical atomic force microscopy (EC-AFM). The presence of PEI on the electrode surface could be distinguished based on nanomechanical properties (DMT modulus), and changes were observed as negative polarization was applied, revealing a reorganization of the PEI chains due to electrostatic interactions. These changes were impacted by the electrolyte composition, including the addition of supporting electrolyte KHCO 3 salt, as well as CO 2 captured by the PEI-based electrolyte, which minimized the change in surface mechanical properties and degree of PEI alignment on the electrode surface. The changes in surface mechanical properties were also dependent on the PEI polymer length, with higher molecular weight PEI showing different reconfiguration than the shorter polymer brushes. The study highlights that the choice of polymer material, the electrolyte composition, and CO 2 captured impact the nearelectrode environment, which has implications for CO 2 R, and presents EC-AFM as a new tool that can be used to probe the dynamic behavior of these interfaces during electrocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Oversaturating Liquid Interfaces with Nanoparticle-Surfactants

Assemblies of nanoparticles at liquid interfaces hold promise as dynamic “active” systems when there are convenient methods to drive the system out of equilibrium via crowding. To this end, we show that oversaturated assemblies of charged nanoparticles can be realized and held in that state with an external electric field. Upon removal of the field, strong interparticle repulsive forces cause a high in-plane electrostatic pressure that is released in an explosive emulsification. We quantify the packing of the assembly as it is driven into the oversaturated state under an applied electric field. Finally, physiochemical conditions substantially affect the intensity of the induced explosive emulsification, underscoring the crucial role of interparticle electrostatic repulsion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2026

Petermann Glacier on the Brink: Progress, Challenges and Insights

Petermann Glacier, the largest marine-terminating glacier in northern Greenland based on catchment area and ice discharge, plays a key role in regulating ice discharge from the Greenland Ice Sheet into the Arctic Ocean. With an upstream catchment connected to the ice sheet interior via a deep subglacial canyon, its future stability has major implications for sea level rise. In this review, we synthesize recent advances in understanding Petermann’s dynamics across three critical interfaces: the ice–ocean, ice–atmosphere, and ice–bed boundaries. At the surface, observations show that reanalysis products underestimate air temperatures and melt, while regional climate models diverge significantly in their estimates of surface mass balance, underscoring the need for improved in situ data and models. At the ocean boundary, enhanced basal melting driven by both subglacial runoff and Atlantic water intrusions is identified as the dominant driver of recent mass loss of Petermann Glacier, with continued warming posing a serious threat to the stability of the floating tongue. At the bed, new geophysical synthesis reveals complex geology, likely spatial variability in geothermal heat flux, and the influence of the megacanyon on seasonal hydrology and velocity fluctuations. Petermann’s mass balance has been negative in recent decades without corresponding flow acceleration. However, the glacier has undergone significant calving events, and the current rifting that began in September 2025 highlights its vulnerability. This upcoming calving event underscores the timeliness of this review, as it will put Petermann Glacier’s terminus at its most retreated position since records began in 1923. The anticipated retreat of the ice tongue also reduces buttressing and brings the terminus closer to the grounding zone, and modeling studies suggest that calving within 12 km of the grounding zone could potentially trigger dynamic retreat, accelerating ice discharge, and a doubling of flow speeds. We conclude that improved observations, sustained monitoring of oceanographic and atmospheric properties, and high-resolution modeling are critical to constraining projections of Petermann Glacier’s future and its role in the stability of the Greenland Ice Sheet.

Dominik Fahrner↗
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2024

Accelerated Sequence Design of Star Block Copolymers: An Unbiased Exploration Strategy via Fusion of Molecular Dynamics Simulations and Machine Learning

Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatibilization, chemical transformations, and separations. s-BCPs have chain architectures where three or more linear diblock copolymer arms comprised of two chemically distinct linear polymers, e.g., solvophobic and solvophilic chains, are covalently joined at one point. The chemical composition of each of the subunit polymer chains comprising the arms, their molecular weights, and the number of arms can be varied to tailor the surface and interfacial activity of these architecturally unique molecules. Further, this makes identification of the optimal s-BCP design nontrivial as the total number of plausible s-BCP architectures is experimentally or computationally intractable. In this work, we use molecular dynamics (MD) simulations coupled with a reinforcement learning-based Monte Carlo tree search (MCTS) to identify s-BCP designs that minimize the interfacial tension between polar and nonpolar solvents. We first validate the MCTS approach for the design of small- and medium-sized s-BCPs and then use it to efficiently identify sequences of copolymer blocks for large-sized s-BCPs. The structural origins of interfacial tension in these systems are also identified by using the configurations obtained from MD simulations. Chemical insights into the arrangement of copolymer blocks that promote lower interfacial tension were mined using machine learning (ML) techniques. Overall, this work provides an efficient approach to solve design problems via fusion of simulations and ML and provides important groundwork for future experimental investigation of s-BCPs for various applications.

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

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

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

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

From Femtoseconds to Gigaseconds: The SolDeg Platform for the Performance Degradation Analysis of Silicon Heterojunction Solar Cells

Heterojunction Si solar cells exhibit notable performance degradation. Here, we modeled this degradation by electronic defects getting generated by thermal activation across energy barriers over time. To analyze the physics of this degradation, we developed the SolDeg platform to simulate the dynamics of electronic defect generation. First, femtosecond molecular dynamics simulations were performed to create a-Si/c-Si stacks, using the machine learning-based Gaussian approximation potential. Second, we created shocked clusters by a cluster blaster method. Third, the shocked clusters were analyzed to identify which of them supported electronic defects. Fourth, the distribution of energy barriers that control the generation of these electronic defects was determined. Fifth, an accelerated Monte Carlo method was developed to simulate the thermally activated time-dependent defect generation across the barriers. Our main conclusions are as follows. (1) The degradation of a-Si/c-Si heterojunction solar cells via defect generation is controlled by a broad distribution of energy barriers. (2) We developed the SolDeg platform to track the microscopic dynamics of defect generation across this wide barrier distribution and determined the time-dependent defect density N(t) from femtoseconds to gigaseconds, over 24 orders of magnitude in time. (3) We have shown that a stretched exponential analytical form can successfully describe the defect generation N(t) over at least 10 orders of magnitude in time. (4) We found that in relative terms, Voc degrades at a rate of 0.2%/year over the first year, slowing with advancing time. (5) We developed the time correspondence curve to calibrate and validate the accelerated testing of solar cells. We found a compellingly simple scaling relationship between accelerated and normal times t normal ∝ t accel T(accel)/T(normal) . (6) We also carried out experimental studies of defect generation in a-Si:H/c-Si stacks. We found a relatively high degradation rate at early times that slowed considerably at longer time scales.

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

Free Energy and Flexibility Analysis of Autoinhibited Human BRAF

The RAF serine/threonine protein kinases function as direct effectors of RAS in the intracellular transmission of extracellular growth signals, and they are key targets for drug discovery, given the high incidence of oncogenic mutations in RAF and other components of this signaling pathway. In its inactive state, RAF is held in an autoinhibited conformation in the cytosol through a combination of intramolecular interactions and binding to a regulatory 14−3−3 protein dimer. Activation of RAF is initiated by its interaction with membrane-localized GTP-bound RAS, which induces conformational changes that release RAF from its autoinhibited state. However, the molecular mechanisms governing RAF activation remain incomplete, largely due to the challenges in experimentally capturing the intermediate conformational states in this process. To address this gap, we developed a comprehensive all-atom model of BRAF based on existing cryo-EM structures. Using this model, we performed extensive molecular dynamics simulations to evaluate the stability and free energy landscape of autoinhibited BRAF in solution. Our analysis reveals conformational flexibility within the autoinhibited complex, suggesting that this dynamic behavior may play a role in facilitating BRAF activation upon engagement with the membrane-bound RAS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

Combined Gravity Gradient and Jitter Accelerations Acting on Liquid-Vapor Interface Oscillations in Reduced Gravity

The dynamical behavior of fluids affected by the asymmetric combined gravity gradient and jitter accelerations, in particular the effect of surface tension on partially-filled rotating fluids applicable to a full-scale Gravity Probe-B Spacecraft dewar tank, have been investigated. Three different cases of accelerations, one gravity gradient-dominated, one equally weighted between gravity gradient and jitter, and the others gravity jitter-dominated are studied. Results of slosh wave excitation along the liquid-vapor interface induced by gravity gradient-dominated acceleration indicate that the gravity gradient-dominated acceleration is equivalent to the combined effect of a twisting force and torsional moment acting on the spacecraft. Results of the slosh wave excitation along the liquid vapor interface induced by gravity jitter-dominated acceleration indicate that the gravity jitter-dominated acceleration is equivalent to time-dependent oscillatory forces which push the bubble in the combined directions of down-and-up and sideward -and-middleward as the bubble is rotating with respect to rotating dewar axis. This study discloses the slosh wave excitation along the liquid-vapor interface driven by the combined effects of gravity gradient and jitter accelerations which are two major driving forces affecting the stability of the fluid system in microgravity.

Hung, R. J.↗
DOE OSTI2024

Oversaturating Liquid Interfaces with Nanoparticle‐Surfactants

Abstract Assemblies of nanoparticles at liquid interfaces hold promise as dynamic “active” systems when there are convenient methods to drive the system out of equilibrium via crowding. To this end, we show that oversaturated assemblies of charged nanoparticles can be realized and held in that state with an external electric field. Upon removal of the field, strong interparticle repulsive forces cause a high in‐plane electrostatic pressure that is released in an explosive emulsification. We quantify the packing of the assembly as it is driven into the oversaturated state under an applied electric field. Physiochemical conditions substantially affect the intensity of the induced explosive emulsification, underscoring the crucial role of interparticle electrostatic repulsion.

Wu, Xuefei↗
DOE OSTI2023

Automated extraction of interfacial dislocations and disconnections from atomistic data

Here we introduce Interfacial Line Defect Analysis (ILDA), a method for identifying and extracting interfacial dislocations and disconnections with little-to-no input from the user on the nature of the interface. By simply providing the orientations and coherency strains for the crystals in a coherent reference state of the interface, ILDA provides exact Burgers vectors. Alternatively, these orientations and strains can be estimated using local atomic deformation gradients, making ILDA fully automated and providing estimated Burgers vectors. ILDA also determines the step height associated with each defect line segment, in case the associated defect is a disconnection. The heart of the method is the identification of atoms residing at co-incidence sites between the two crystals and the construction of a surface mesh connecting these sites that is used to compose Burgers circuits and insert defect line segments. We demonstrate the performance of ILDA in two test cases: a twist grain boundary and a phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Force field for calcium sulfate minerals to predict structural, hydration, and interfacial properties

Calcium sulfates such as anhydrite, hemihydrate, and gypsum find widespread use in building materials, implants, and tissue healing. We introduce a simple and compatible atomistic force field for all calcium sulfate phases that reproduces a wide range of experimental data including lattice parameters, surface, hydration, mechanical, and thermal properties in 1% to 5% accuracy relative to experiments. The performance is several times better than prior force fields and DFT methods, which lead to errors in structures and energies up to 100%. We explain (hkl) cleavage energies, the dynamics of (hkl) water interfaces, and new insights into molecular origins of crystal-facet specific hydration and solubility. Impressive agreement of computed and experimentally measured hydration energies is shown. The models add to the Interface force field (IFF) and are compatible with multiple force fields (CHARMM, AMBER, GROMOS, CVFF, PCFF, OPLS-AA) for property predictions of sulfate-containing materials from atoms to the large nanometer scale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular aspects of temperature swing solvent extraction for brine desalination using imidazole-based solvents

Solvent extraction desalination has shown to be a promising energy-efficient technique when compared to separation processes usually employed. Despite the remarkable impact of the organic solvent’s molecular characteristics on this process, few studies have been carried out so far to explicitly ascertain the influence of these characteristics on the solubility of water in the solvent. In this report we assess the viability of four different imidazole-based solvents and compare them to an amine compound for solvent extraction desalination of high salinity brines. All-atom molecular simulations of pure solvent and water-solvent bulk mixtures are used to ascertain the influence of the solvents’ molecular structures on solvent–solvent and solvent–water interaction. Our results show that water solubilization in imidazole-based solvents is an endothermic process, while for a secondary amine (diisopropylamine) it is an exothermic process. We show that imidazole-based solvents can be designed to provide characteristics favorable to brine desalination.

42 ENGINEERING↗
DOE OSTI2022

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗