Search NASA⌕ Search

SEARCH · Search NASA

Results for “solvation 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗
DOE OSTI2020

Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale surface features of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. Additionally, a simple analytical model for considering the average effect of such surface features producing nanoscale roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to unambiguously determine the zeta potential in the presence of nanoscale surface features.

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

Influence of solution phase environmental heterogeneity and fluctuations on vibronic spectra: Perylene diimide molecular chromophore complexes in solution

Ensembles of ab initio parameterized Frenkel-exciton model Hamiltonians for different perylene diimide dimer systems are used, together with various dissipative quantum dynamics approaches, to study the influence of the solvation environment and fluctuations in chromophore relative orientation and packing on the vibronic spectra of two different dimer systems: a π-stacked dimer in aqueous solution in which the relative chromophore geometry is strongly confined by a phosphate bridge and a side-by-side dimer in dichloromethane involving a more flexible alkyne bridge that allows quasi-free rotation of the chromophores relative to one another. These entirely first-principles calculations are found to accurately reproduce the main features of the experimental absorption spectra, providing a detailed mechanistic understanding of how the structural fluctuations and environmental interactions influence the vibronic dynamics and spectroscopy of solutions of these multi-chromophore complexes.

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

Charge accumulation and solvation in $β$-NiOOH: Surface chemistry of an OER catalyst from ML-aided simulations

Electrochemical water splitting is a key technology for a sustainable energy transition, providing a route to store surplus electricity from renewable sources. A central bottleneck is the sluggish oxygen evolution reaction (OER), which drives the search for catalysts that are active, stable, and inexpensive enough for large-scale deployment. Within this context, pure and doped NiO x H y combine high activity with low cost, making them prime candidates for alkaline OER. Yet, despite extensive study, the atomistic structure of NiOOH under operando conditions and the associated reaction mechanisms remain debated. Here, we investigate the structural complexity of pure β-NiOOH, the scaffold for its doped derivatives. We systematically investigate the oxidation of the surface adsorbates via proton-coupled electron transfer steps across relevant facets and sites, identifying the most probable sequence of deprotonation events. Our results reveal asymmetric charge accumulation on Wulff-relevant surfaces and show how applied potential can promote morphological restructuring. Explicit solvation is included through machine-learning interatomic potential molecular dynamics of the NiOOH/water interface, which allows us to resolve the hydrophobic and hydrophilic character of different surfaces and the associated interfacial water structure. Together, these insights demonstrate how surface chemistry and solvation jointly govern the stability of NiOOH and the accumulation of surface charge, with possible implications for catalytic performance.

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

Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning

Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.

Adsorption↗
DOE OSTI2021

Effect of water models on structure and dynamics of lignin in solution

Lignin, a major biomass component, can be an excellent source for different monomers in the polymer industry. However, the complex and heterogeneous structure of lignin poses a significant challenge for designing energy-efficient processes for depolymerization. As many proposed depolymerization processes are solvothermal, it is essential to understand the structure and dynamics of lignin in solution, in particular aqueous solution. Here, we utilize molecular dynamics simulations to understand the effect of water models on the structure and dynamics of different model lignin oligomers (softwood and hardwood) as a function of temperature. We have examined three different water models: TIP3P, TIP4P/Ew, and flexible SPC/Fw. We find that the diffusion constant of lignin oligomers in an aqueous solution differs significantly depending on the water model used. The diffusion constant of lignin in the TIP3P water model is almost twice as large as that in SPC/Fw and TIP4P/Ew. The softwood and hardwood oligomers adopt an extended structure in TIP3P water compared to SPC/Fw and TIP4P/Ew. Given the different levels of sensitivity of transport and structural properties of aqueous lignin on water models, it is important to take these into account when discussing results from a specific water model.

09 BIOMASS FUELS↗
DOE OSTI2025

A Perspective on Multiscale Modeling of Explicit Solvation-Enabled Simulations of Catalysis at Liquid–Solid Interfaces

Catalysis at liquid-solid interfaces is profoundly influenced by the interfacial solvent structure, which affects catalytic activity, selectivity, and reaction pathways. This perspective discusses state-of-the-art multiscale modeling methods that integrate quantum mechanics and molecular mechanics approaches to apply explicit solvent molecules to capture these interfacial phenomena. Specifically, the construction of multiscale models, the importance of capturing the interfacial solvent structure, and the computational strategies used to achieve this are explored, and the challenges in balancing chemical accuracy with computational expense are highlighted. Additionally, this perspective addresses the limitations of current methods. Opportunities for integrating machine learning are proposed. Here, by advancing the efficiency and user friendliness of multiscale modeling, it is argued that deeper insights into heterogeneous catalysis in liquid phases can be provided, which will ultimately contribute to the development of more efficient catalytic processes.

Ab initio molecular dynamics↗
DOE OSTI2025

Water-mediated ion transport in an anion exchange membrane

Water is a critical component in polyelectrolyte anion exchange membranes (AEMs). It plays a central role in ion transport in electrochemical systems. Gaining a better understanding of molecular transport and conductivity in AEMs has been challenged by the lack of a general methodology capable of capturing and connecting water dynamics, water structure, and ionic transport over time and length scales ranging from those associated with individual bond vibrations and molecular reorientations to those pertaining to macroscopic AEM performance. In this work, we use two-dimensional infrared spectroscopy and semiclassical simulations to examine how water molecules are arranged into successive solvation shells, and we explain how that structure influences the dynamics of bromide ion transport processes in polynorbornene-based materials. We find that the transition to the faster transport mechanism occurs when the reorientation of water molecules in the second solvation shell is fast, allowing a robust hydrogen bond network to form. Our findings provide molecular-level insights into AEMs with inherent transport of halide ions, and help pave the way towards a comprehensive understanding of hydroxide ion transport in AEMs.

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

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

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

Structural Changes to the Gd‐DTPA Complex at Varying Ligand Protonation State

Abstract Diethylenetriaminepentaacetic acid (DTPA) is a chelating agent whose complex with the Gd 3+ ion is used in medical imaging. DTPA is also used in lanthanide‐actinide separation processes. As protonation of the DTPA ligand can facilitate dissociation of the Gd 3+ ion from the Gd‐DTPA complex, this work investigates the coordination structures of the aqueous Gd 3+ ion and its environment when chelated by DTPA in eight different DTPA protonation states. Both classical and ab initio molecular dynamics (MD) simulations are conducted to model the solvated complexes. Extended X‐ray absorption fine structure (EXAFS) measurements of the Gd 3+ aqua ion, and the Gd‐DTPA complex at pH 1 and 11, are compared to EXAFS spectra predicted from the MD simulations to verify the accuracy of the MD structures. The findings of this work provide atomic‐level details into the fluctuating Gd‐DTPA complex environment as the DTPA ligand gradually detaches from the Gd 3+ ion with increased protonation.

Chemistry↗
DOE OSTI2024

Weakly solvating ester electrolyte for high voltage sodium-ion batteries

Ethyl acetate (EA) was identified as a promising electrolyte solvent for sodium-ion batteries (SIBs), exhibiting low viscosity, cost-effectiveness, and low toxicity. Despite a significant portion of aggregation being linked to the weak solvation of Na + /EA as revealed by molecular dynamics (MD) simulations, pulsed-field gradient nuclear magnetic resonance (pfg-NMR) analysis identified a noteworthy Na+ diffusion coefficient of 3.95×10 -10 m 2 s -1 at 25°C in the presence of 1 m NaPF 6 salt. Employing fluoroethylene carbonate (FEC) as a film-forming additive to create electrode-electrolyte interphase, this electrolyte surprisingly made ~210 mAh Na 0.97 Ca 0.03 [Mn 0.39 Fe 0.31 Ni 0.22 Zn 0.08 ]O 2 (NCMFNZO)/hard carbon (HC) pouch cells achieve a lengthy cycling lifetime of 250 cycles with ~80 % capacity retention, cycled up to 4.0 V at 40°C. X-ray photoelectron spectroscopy (XPS) revealed increasing interphasial organic species over cycling, augmenting charge transfer resistance on both cathode and anode, particularly during fast charging or low temperatures (<10°C), promoting Na plating. Finally, gas chromatography-mass spectrometry combined with density functional theory identified CO 2 as the major gas generated from charged cathode/electrolyte interactions, exhibiting temperature/voltage dependence.

25 ENERGY STORAGE↗
DOE OSTI2022

Molecular engineering of fluoroether electrolytes for lithium metal batteries

Fluoroether solvents are promising electrolyte candidates for high-energy-density lithium metal batteries, where high ionic conductivity and oxidative stability are important metrics for design of new systems. Recent experiments have shown that these performance metrics, particularly stability, can be tuned by changing the fraction of ether and fluorine content. However, little is known about how different molecular architectures influence the underlying ion transport mechanisms and conductivity. Here, we use all-atom molecular dynamics simulations to elucidate the ion transport and solvation characteristics of fluoroether chains of varying length, and having different ether segment and fluorine terminal group contents. The design rules that emerge from this effort are that solvent size determines lithium-ion transport kinetics, solvation structure, and solvation energy. In particular, the mechanism for lithium-ion transport is found to shift from ion hopping between solvation sites located in different fluoroether chains in short-chain solvents, to ion–solvent co-diffusion in long-chain solvents, indicating that an optimum exists for molecules of intermediate length, where hopping is possible but solvent diffusion is fast. Consistent with these findings, our experimental measurements reveal a non-monotonic behavior of the effects of solvent size on lithium-ion conductivity, with a maximum occurring for medium-length solvent chains. Finally, a key design principle for achieving high ionic conductivity is that a trade-off is required between relying on shorter fluoroether chains having high self-diffusivity, and relying on longer chains that increase the stability of local solvation shells.

25 ENERGY STORAGE↗
DOE OSTI2021

First-Principles Grand-Canonical Simulations of Water Adsorption in Proton-Exchanged Zeolites

Water appears by design or as impurities in many important reactive systems. For those catalyzed by porous solid acids, such as widely used zeolites, experimentally quantifying the amount or elucidating the structure of adsorbed water clusters at reaction conditions is challenging, while computational studies (e.g., first-principles molecular dynamics simulations) often need to assume the loading to examine solvation effects. Furthermore, we perform first-principles grand-canonical simulations to predict water adsorption to H-ZSM-5 zeolites under specified experimental conditions. Presampling with inexpensive force fields and a pool-based parallelization algorithm are used to improve simulation efficiency, while molecular dynamics is used to sample configurations involving hydronium species. We observe that H + exchange dramatically increases the hydrophilicity of zeolite MFI and an appreciable amount of water is present at very low relative humidities. At all conditions examined, the zeolitic protons are found to dissociate readily from the surface basic sites and become mobile by participating in the hydrogen-bonded chains of adsorbed water molecules.

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

The Polarity of Co-solvents Regulates the Charge Storage Mechanisms in Supercapacitors with Concentrated Electrolytes

Developing better energy storage devices depends on comprehending the underlying mechanisms involved in charge storage. With the continuous conception of new electrolytes, this task becomes progressively more urgent and complex. An example is the utilization of co-solvated concentrated solutions. While these show promising electrochemical responses, their dynamic properties (especially under confinement) and their relationships with performance are not fully understood. Here, we combined modified step potential electrochemical spectroscopy and quasielastic neutron scattering to investigate systems composed of activated mesoporous carbon (AMC) and concentrated solutions of lithium bis(trifluoromethanesulfonyl)imide in acetonitrile co-solvated with either toluene or acetone. Further, we report that acetone does not impair surface-controlled mechanisms, contrary to the case with toluene, which competes with charged species to populate the AMC’s pores without contributing to charge storage. In turn, toluene promotes a greater overall capacitance owing to Faradaic processes, which may be related to changes in the solvation structures under confinement.

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

High-entropy electrolytes for practical lithium metal batteries

Electrolyte engineering is crucial for improving battery performance, particularly for lithium metal batteries. Recent advances in electrolytes have greatly improved cyclability by enhancing electrochemical stability at the electrode interfaces, but concurrently achieving high ionic conductivity has remained challenging. Here we report an electrolyte design strategy for enhanced lithium metal batteries by increasing the molecular diversity in electrolytes, which essentially leads to high entropy electrolytes (HEEs). We find that in weakly solvating electrolytes, the entropy effect reduces ion clustering while preserving the characteristic anion-rich solvation structures, which is characterized by synchrotron-based X-ray scattering and molecular dynamics simulations. Electrolytes with smaller- sized clusters exhibit a 2-fold improvement in ionic conductivity compared to conventional weakly- solvating electrolytes, enabling stable cycling at high current densities up to 2C (6.2 mA cm -2 ) in anode- free LiNi 0.6 Mn 0.2 Co 0.2 (NMC622)||Cu pouch cells. In conclusion, the efficacy of the design strategy is verified by performance improvements in three disparate weakly solvating electrolyte systems.

25 ENERGY STORAGE↗
DOE OSTI2021

Ab-Initio Simulation Studies of Cr Solvation in Fluoride Molten Salts

Understanding molten salt chemistry is essential in ongoing research of the molten salt nuclear reactor (MSR). In this context, detailed understanding of the mechanisms underlying selective oxidation of metal species, such as Cr, is required to guide the design of effective corrosion mitigation strategies in molten salts. An important starting point for such mechanistic understanding is knowledge of the solvation structure and its role in controlling metal speciation. In this work, we use ab initio molecular dynamics simulations to study the short-range (on the scale of the nearest-neighbor bond lengths) and medium-range (over length scales of several neighbor spacings) structure in three different fluoride melts with and without Cr addition; namely, 2KF-NaF, 2LiF-BeF 2 , and 3LiF-AlF3. We find that Cr 0 ,Cr 2+ ,Cr 3+ can each be coordinated by different numbers of F-, with the variance in coordination number decreasing as oxidation state increases, and that these coordination geometries are largely independent of solvent. The manner by which Cr changes the medium-range structure, however, is found to be solvent-dependent. While 2KF-NaF melts show short and medium range order that is highly dynamic, 2LiF-BeF 2 and 3LiF-AlF 3 are characterized by molecular associates that are relatively long-lived that organize into oligomer structures on larger length scales. Rather than being solvated by F- ions alone, we find that Cr can incorporate into and be solvated within this oligomer structure. Fluoroacidity, alone, may therefore prove too simple a metric for assessing the corrosivity of molten fluorides. As our work suggests, the ability of Cr to solvate must be understood in the context of the short- and medium-range structure of the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗