Search NASA⌕ Search

SEARCH · Search NASA

Results for “Contact thermodynamics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Thermodynamic driving forces in contact electrification between polymeric materials

Abstract Contact electrification, or contact charging, refers to the process of static charge accumulation after rubbing, or even simple touching, of two materials. Despite its relevance in static electricity, various natural phenomena, and numerous technologies, contact charging remains poorly understood. For insulating materials, even the species of charge carrier may be unknown, and the direction of charge-transfer lacks firm molecular-level explanation. Here, we use all-atom molecular dynamics simulations to investigate whether thermodynamics can explain contact charging between insulating polymers. Based on prior work suggesting that water-ions, such as hydronium and hydroxide ions, are potential charge carriers, we predict preferred directions of charge-transfer between polymer surfaces according to the free energy of water-ions within water droplets on such surfaces. Broad agreement between our predictions and experimental triboelectric series indicate that thermodynamically driven ion-transfer likely influences contact charging of polymers. Furthermore, simulation analyses reveal how specific interactions of water and water-ions proximate to the polymer-water interface explain observed trends. This study establishes relevance of thermodynamic driving forces in contact charging of insulators with new evidence informed by molecular-level interactions. These insights have direct implications for future mechanistic studies and applications of contact charging involving polymeric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Kinetic versus Thermodynamic Stability of LLZO in Contact with Lithium Metal

Li 7 La 3 Zr 2 O 12 (LLZO) garnet-based oxides are a promising class of solid electrolytes used as the separator in all-solid-state batteries (ASSBs). While LLZO is considered to have a wide electrochemical stability window, its intrinsic stability in contact with lithium metal is not sufficiently well understood, and there is still a debate on the key question of whether LLZO does or does not form passivation layers before and during cycling. Utilizing both in situ and operando X-ray photoelectron spectroscopy techniques, in this work we reveal the presence of a kinetic barrier to the reduction of LLZO by Li metal, with the extent of oxygen-deficient interphase (ODI) formation depending sensitively on the energetics of Li metal arriving at the Li vertical bar LLZO interface. Despite the clear presence of a kinetic barrier to reduction, the electrochemical response of the Li vertical bar LLZO interface is unchanged by the presence of the ODI, indicating that ODI formation during electrochemical cycling does not hinder charge transfer across the Li vertical bar LLZO interface. Overall, these results reveal that the reactivity of LLZO with Li metal depends not only on the material properties of the adjoining phases (i.e., surface purity and active contact) and their resulting thermodynamic stability but also on the energy input at the interface and the resulting reaction kinetics. Furthermore, the presence of a kinetic barrier to reduction highlights the additional complexities governing the reactivity of solid-state interfaces in ASSBs and underscores the importance of operando characterization of interfacial stability to design more robust, high-performance protection strategies for solid electrolytes in contact with reactive electrodes.

25 ENERGY STORAGE↗
DOE OSTI2023

Ultrathin stable Ohmic contacts for high-temperature operation of β -Ga2O3 devices

Beta gallium oxide (β-Ga2O3) shows significant promise in high-temperature, high-power, and sensing electronics applications. However, long-term stable metallization layers for Ohmic contacts at high temperatures present unique thermodynamic challenges. The current most common Ohmic contact design based on 20 nm of Ti has been repeatedly demonstrated to fail at even moderately elevated temperatures (300–400 °C) due to a combination of nonstoichiometric Ti/Ga2O3 interfacial reactions and kinetically favored Ti diffusion processes. Here, we demonstrate stable Ohmic contacts for Ga2O3 devices operating up to 500–600 °C using ultrathin Ti layers with a self-limiting interfacial reaction. The ultrathin Ti layer in the 5 nm Ti/100 nm Au contact stack is designed to fully oxidize while forming an Ohmic contact, thereby limiting both thermodynamic and kinetic instability. This novel contact design strategy results in an epitaxial conductive anatase titanium oxide interface layer that enables low-resistance Ohmic contacts that are stable both under long-term continuous operation (>500 h) at 600 °C in vacuum (≤10−4 Torr), as well as after repeated thermal cycling (15 times) between room temperature and 550 °C in flowing N2. This stable Ohmic contact design will accelerate the development of high-temperature devices by enabling research focus to shift toward rectifying interfaces and other interfacial layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Reliable operation of Cr 2 O 3 :Mg/ $β$-Ga 2 O 3 p–n heterojunction diodes at 600 °C

Beta gallium oxide (β-Ga 2 O 3 )-based semiconductor heterojunctions have recently demonstrated improved performance at high voltages and elevated temperatures and are, thus, promising for applications in power electronic devices and harsh environment sensors. However, the long-term reliability of these ultra-wideband gap (UWBG) semiconductor devices remains barely addressed and may be strongly influenced by chemical reactions at the p–n heterojunction interface. Here, we experimentally demonstrate operation and evaluate the reliability of Cr 2 O 3 :Mg/β-Ga 2 O 3 p–n heterojunction diodes during extended operation at 600 °C, as well as after 30 repeated cycles between 25 and 550 °C. The calculated pO 2 -temperature phase stability diagram of the Ga-Cr-O material system predicts that Ga 2 O 3 and Cr 2 O 3 should remain thermodynamically stable in contact with each other over a wide range of oxygen pressures and operating temperatures. The fabricated Cr 2 O 3 :Mg/β-Ga 2 O 3 p–n heterojunction diodes show room-temperature on/off ratios >104 at ±5 V and a breakdown voltage (V Br ) of -390 V. The leakage current increases with increasing temperature up to 600 °C, which is attributed to Poole–Frenkel emission with a trap barrier height of 0.19 eV. Over the course of a 140-h thermal soak at 600 °C, both the device turn-on voltage and on-state resistance increase from 1.08 V and 5.34 mΩ cm 2 to 1.59 V and 7.1 mΩ cm 2 , respectively. This increase is attributed to the accumulation of Mg and MgO at the Cr 2 O 3 /Ga 2 O 3 interface as observed from the time-of-flight secondary ion mass spectrometry analysis. These findings inform future design strategies of UWBG semiconductor devices for harsh environment operation and underscore the need for further reliability assessments for β-Ga 2 O 3 -based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Role of water model on ion dissociation at ambient conditions

In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.

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

Degradation of Lithium Iron Phosphate Sulfide Solid-State Batteries by Conductive Interfaces

The superionic solid-state argyrodite electrolyte Li 6 PS 5 Br can improve lithium and lithium-ion batteries' safety and energy density. Despite many reports validating the conductivity of this electrolyte, it still suffers from passivating electrode degradation mechanisms. At first analysis, lithium iron phosphate (LFP) should be more thermodynamically stable in contact with sulfide electrolytes. However, without substantial improvements to interfacial engineering, we find that LFP is not inherently stable against Li 6 PS 5 Br. We hypothesize argyrodite oxidation favorably competes with LFPAs delithiation, insulating the electrolyte-electrode interface and causing large overpotential losses. We show that compared to LiNi x Mn y Co z O 2 , LFP has no actual electrochemical stability advantage despite operating at a lower voltage. We utilize tender energy XAS and XPS to show that chemical reactions occur between LFP and the Li 6 PS 5 Br solid electrolyte and these reactions are exacerbated by cycling. We also show that electrochemical degradation occurs at the interface between the solid electrolyte ion conductor and any electron conductor, namely the active material and carbon additives. We further demonstrate that LiNbO 3 cathode coatings on LFP can delay electrochemical degradation by electronically insulating the LFP-sulfide electrolyte interface but not prevent its occurrence at the carbon-electrolyte interface.

25 ENERGY STORAGE↗
DOE OSTI2021

Modeling the electrical double layer at solid-state electrochemical interfaces

Models of the electrical double layer (EDL) at electrode/liquid-electrolyte interfaces no longer hold for all-solid-state electrochemistry. Here we show a more general model for the EDL at a solid-state electrochemical interface based on the Poisson–Fermi–Dirac equation. By combining this model with density functional theory predictions, the interconnected electronic and ionic degrees of freedom in all-solid-state batteries, including the electronic band bending and defect concentration variation in the space-charge layer, are captured self-consistently. Along with a general mathematical solution, the EDL structure is presented in various materials that are thermodynamically stable in contact with a lithium metal anode: the solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) and the solid interlayer materials LiF, Li 2 O and Li 2 CO 3 . Here, the model further allows design of the optimum interlayer thicknesses to minimize the electrostatic barrier for lithium ion transport at relevant solid-state battery interfaces.

25 ENERGY STORAGE↗
DOE OSTI2021

Molecular basis for higher affinity of SARS-CoV-2 spike RBD for human ACE2 receptor

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused substantially more infections, deaths, and economic disruptions than the 2002-2003 SARSCoV. The key to understanding SARS-CoV-2's higher infectivity lies partly in its host receptor recognition mechanism. Experiments show that the human angiotensin converting enzyme 2 (ACE2) protein, which serves as the primary receptor for both CoVs, binds to the receptor binding domain (RBD) of CoV-2's spike protein stronger than SARS-CoV's spike RBD. The molecular basis for this difference in binding affinity, however, remains unexplained from X-ray structures. To go beyond insights gained from X-ray structures and investigate the role of thermal fluctuations in structure, we employ all-atom molecular dynamics simulations. Microseconds-long simulations reveal that while CoV and CoV-2 spike-ACE2 interfaces have similar conformational binding modes, CoV-2 spike interacts with ACE2 via a larger combinatorics of polar contacts, and on average, makes 45% more polar contacts. Correlation analysis and thermodynamic calculations indicate that these differences in the density and dynamics of polar contacts arise from differences in spatial arrangements of interfacial residues, and dynamical coupling between interfacial and non-interfacial residues. Furthermore, these results recommend that ongoing efforts to design spike-ACE2 peptide blockers will benefit from incorporating dynamical information as well as allosteric coupling effects.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. Here, we also present results for the energy as a function of temperature in the continuum limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1998

Dynamics of Liquids in Edges and Corners (DYLCO): IML-2 Experiment for the BDPU

Knowledge of the behavior of fluids possessing free surfaces is important to many fluid systems, particularly in space, where the normally subtle effects of surface wettability play a more dramatic and often surprising role. DYLCO for the IML-2 mission was proposed as a simple experiment to probe the particular behavior of capillary surfaces in containers of irregular cross section. Temperature control was utilized to vary the fluid-solid contact angle, a questionable thermodynamic parameter of the system, small changes in which can dramatically influence the configuration, stability, and flow of a capillary surface. Container shapes, test fluid, and temperature ranges were selected for observing both local changes in interface curvature as well as a global change in fluid orientation due to a critical wetting phenomenon. The experiment hardware performed beyond what was expected and fluid interfaces could be readily digitized post flight to show the dependence of the interface curvature on temperature. For each of the containers tested surfaces were observed which did not satisfy the classic equations for the prediction of interface shape with constant contact angle boundary condition. This is explained by the presence of contact angle hysteresis arising from expansion and contraction of the liquid during the heating and cooling steps of the test procedure. More importantly, surfaces exceeding the critical surface curvature required for critical wetting were measured, yet no wetting was observed. These findings are indeed curious and pose key questions concerning the role of hysteresis for this critical wetting phenomena. The stability of such surfaces was determined numerically and it is shown that stability is enhance (reduced) when a surface is in its 'advancing' ('receding') state, The analysis shows complete instability as the critical wetting condition is reached. The case of ideal dynamic wetting is addressed analytically in detail with results of significant flow characteristics presented in closed form. The solutions indicate a square root of T dependence of the capillary 'rise' rate which is corroborated by drop tower tests. The analysis clearly shows that infinite time is necessary for surfaces to reorient at the critical wetting transition.

Langbein, D.↗
DOE OSTI2025

Explaining Snowball-in-Hell Phenomena in Heavy-Ion Collisions Using a Novel Thermodynamic Variable

A loosely bound hadronic molecule produced by a relativistic heavy-ion collision has been described as a “snowball in hell” since it emerges from a hadron resonance gas whose temperature is orders of magnitude larger than the binding energy of the molecule. This remarkable phenomenon can be explained in terms of a novel thermodynamic variable called the “contact” that is conjugate to the binding momentum of the molecule. The production rate of the molecule can be expressed in terms of the contact density at the kinetic freeze-out of the hadron resonance gas. It approaches a nonzero limit as the binding energy goes to 0.

Relativistic heavy-ion collisions↗
NASA NTRS2004

Unfolding of Proteins: Thermal and Mechanical Unfolding

We have employed a Hamiltonian model based on a self-consistent Gaussian appoximation to examine the unfolding process of proteins in external - both mechanical and thermal - force elds. The motivation was to investigate the unfolding pathways of proteins by including only the essence of the important interactions of the native-state topology. Furthermore, if such a model can indeed correctly predict the physics of protein unfolding, it can complement more computationally expensive simulations and theoretical work. The self-consistent Gaussian approximation by Micheletti et al. has been incorporated in our model to make the model mathematically tractable by signi cantly reducing the computational cost. All thermodynamic properties and pair contact probabilities are calculated by simply evaluating the values of a series of Incomplete Gamma functions in an iterative manner. We have compared our results to previous molecular dynamics simulation and experimental data for the mechanical unfolding of the giant muscle protein Titin (1TIT). Our model, especially in light of its simplicity and excellent agreement with experiment and simulation, demonstrates the basic physical elements necessary to capture the mechanism of protein unfolding in an external force field.

Hur, Joe S.↗
DOE OSTI2026

Mechanism of ion pairing–induced co-ion penetration into the stern layer

In classical electric double layer (EDL) theory, co-ions—ions carrying the same charge as the electrode—are assumed to be excluded from the Stern layer because of electrostatic repulsion. However, ion pairing with adsorbed counterions may enable co-ions to access the EDL, a phenomenon that remains underexplored. Key questions therefore remain unresolved, including how deeply co-ions penetrate, how ion-specific properties govern them, and how ion pairing within the EDL differs from bulk solution. Here, in this study, we employ molecular dynamics simulations combined with well-tempered metadynamics to quantify co-ion penetration and ion-pairing thermodynamics at positively charged Mg-Al layered double hydroxide interfaces (σ = +20.6 μC/cm 2 ) in aqueous chloride electrolytes across 0.005–3 M range. We reveal a concentration-driven crossover: at dilute conditions, counterion–co-ion contact ion pairing (CIP) is thermodynamically unfavorable or only weakly favorable (ΔG CIP = +0.12 kcal/mol for Li + to −0.42 kcal/mol for K + ), leading to co-ion exclusion. As concentration increases (≥0.5 M), enhanced screening reduces surface repulsion and stabilizes CIP thermodynamics (ΔG CIP = −0.89 kcal/mol for Li + to −1.05 kcal/mol for K + at 3 M), enabling monovalent co-ions to penetrate the Stern layer following the hydration hierarchy K + > Na + > Li + , while divalent co-ions (Mg 2+ , Ca 2+ ) remain in the diffuse layer, forming only solvent-separated pairs. Comparison with bulk solution shows interfacial ion pairing is suppressed at low concentrations but converges to bulk-like thermodynamics at high concentrations, indicating surface effects vanish under strong screening. These results establish a quantitative framework linking ion-pairing thermodynamics to co-ion penetration and providing strategies to design electrochemical interfaces.

Co-ion penetration↗
DOE OSTI2025

Thermally induced PtOx interfacial layer enhances stability of Pt/β-Ga2O3 vertical Schottky diodes

This study investigates the long-term stability of Pt/β-Ga2O3 field-plated Schottky barrier diodes at high temperatures, with extended thermal cycling and soaking stress. The device characteristics reveal a strong dependence on operating temperature, leading to an ON/OFF current ratio of ∼10× at ±3 V and 600 °C. Extensive thermal cycling (25 cycles over 100 h) between room temperature and 350 °C shows that leakage currents stabilize after the initial ten cycles. This stabilization is likely due to the formation of a thin PtOx layer at the Pt/β-Ga2O3 interface, as confirmed by high-resolution scanning transmission electron microscopy. Further exposure to a constant 350 °C environment for 10 h reveals a slight increase in the leakage current from 10−7 to 10−6 A at −5 V and degradation of the forward current, due to potential thermodynamic instability of the Ohmic contacts. This study provides insights into the changes in the Pt/β-Ga2O3 diode performance with thermal stress, offering a basis for predicting the device’s operational reliability at high temperatures for sensing and power applications.

08 HYDROGEN↗
DOE OSTI2026

Nonlinear thermodynamic computing out of equilibrium

We present the design for a thermodynamic computer that can perform arbitrary nonlinear calculations in or out of equilibrium. Simple thermodynamic circuits, fluctuating degrees of freedom in contact with a thermal bath and confined by a quartic potential, display an activity that is a nonlinear function of their input. Such circuits can therefore be regarded as thermodynamic neurons, and can serve as the building blocks of networked structures that act as thermodynamic neural networks, universal function approximators whose operation is powered by thermal fluctuations. We simulate a digital model of a thermodynamic neural network, and show that its parameters can be adjusted by genetic algorithm to perform nonlinear calculations at specified observation times, regardless of whether the system has attained thermal equilibrium. This work expands the field of thermodynamic computing beyond the regime of thermal equilibrium, enabling fully nonlinear computations, analogous to those performed by classical neural networks, at specified observation times.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
DOE OSTI2026

Molecular Dynamics Investigation of Oil Wetting on Synthetic Polymer Substrates

Understanding the interaction of polymer surfaces with nonpolar, low surface tension liquids, or whether a substrate is oleophobic versus oleophilic, is critical for applications ranging from antifingerprint coatings to oil–water separation membranes and oil spill remediation. Despite its technological importance, the molecular mechanisms underpinning polymer–oil wetting are not well characterized. Here, we employ molecular dynamics simulations to investigate the behavior of n-hexadecane in contact with chemically distinct polymer surfaces, spanning eight constitutional unit chemistries as well as amorphous and crystalline morphologies. This permits a critical examination of both thermodynamic and dynamic descriptors of oleophobicity, including oil contact angle, dewetting free energy, interfacial diffusivity, and a proposed “ghost probe energy,” which does not require explicit simulation of oil–polymer interactions. We find that the oil contact angle does not reliably distinguish oleophobic behavior across polymer chemistries, whereas other metrics provide clearer and more consistent differentiation. Analysis of the results reveals that polymer–oil wetting behavior is primarily governed by interfacial van der Waals interactions and modulated by surface flexibility and morphology. Collectively, this work establishes a computational framework for characterizing oil wetting and provides additional insight into what molecular-level factors dictate trends in oleophobicity.

lipids↗