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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗
DOE OSTI2023

Solvation of furfural at metal–water interfaces: Implications for aqueous phase hydrogenation reactions

Metal–water interfaces are central to understanding aqueous-phase heterogeneous catalytic processes. However, the explicit modeling of the interface is still challenging as it necessitates extensive sampling of the interfaces’ degrees of freedom. Herein, we use ab initio molecular dynamics (AIMD) simulations to study the adsorption of furfural, a platform biomass chemical on several catalytically relevant metal–water interfaces (Pt, Rh, Pd, Cu, and Au) at low coverages. We find that furfural adsorption is destabilized on all the metal–water interfaces compared to the metal–gas interfaces considered in this work. This destabilization is a result of the energetic penalty associated with the displacement of water molecules near the surface upon adsorption of furfural, further evidenced by a linear correlation between solvation energy and the change in surface water coverage. To predict solvation energies without the need for computationally expensive AIMD simulations, we demonstrate OH binding energy as a good descriptor to estimate the solvation energies of furfural. Using microkinetic modeling, we further explain the origin of the activity for furfural hydrogenation on intrinsically strong-binding metals under aqueous conditions, i.e., the endothermic solvation energies for furfural adsorption prevent surface poisoning. In conclusion, our work sheds light on the development of active aqueous-phase catalytic systems via rationally tuning the solvation energies of reaction intermediates.

aqueous-phase hydrogenation↗
DOE OSTI2021

An Environmentally Benign Electrolyte for High Energy Lithium Metal Batteries

In this work, a hybrid electrolyte comprising a high content of H 2 O for a lithium metal cell is reported. At high LiFSI salt concentration, the N-methyl-N-propyl-piperidinium bis(fluorosulfonyl) imide (PMpipFSI) electrolyte can tolerate up to 1 M H2O addition without sacrificing its redox stability on both lithium nickel manganese cobalt oxide (NMC) cathode and lithium metal anode. Molecular dynamics simulations revealed the underpinned mechanism that, at high salt concentrations, H 2 O molecules are embedded in the Li + , PMpip + , and FSI - bulk as a structural material with a strong solvation with Li + and are orderly distributed at the surface of both electrodes. This electrolyte eliminates the critical moisture controls required for the state-of-the-art (SOA) carbonate/LiPF 6 electrolyte, electrode, separator and cell assembly, thus significantly reducing the cost of the mass production of the batteries.

molecular dynamics simulation↗
DOE OSTI2020

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Alchemical free energy simulations have long been utilized to predict free energy changes for binding affinity and solubility of small molecules. However, while the theoretical foundation of these methods is well established, seamlessly handling many of the practical aspects regarding the preparation of the different thermodynamic end states of complex molecular systems and the numerous processing scripts often remains a burden for successful applications. Here, we present CHARMM-GUI Free Energy Calculator (http://www.charmm-gui.org/input/fec) that provides various alchemical free energy perturbation molecular dynamics (FEP/MD) systems with input and post-processing scripts for NAMD and GENESIS. Four submodules are available: Absolute Ligand Binder (for absolute ligand binding FEP/MD), Relative Ligand Binder (for relative ligand binding FEP/MD), Absolute Ligand Solvator (for absolute ligand solvation FEP/MD), and Relative Ligand Solvator (for relative ligand solvation FEP/MD). Each module is designed to build multiple systems of a set of selected ligands at once for high-throughput FEP/MD simulations. The capability of Free Energy Calculator is illustrated by absolute and relative solvation FEP/MD of a set of ligands and absolute and relative binding FEP/MD of a set of ligands for T4-lysozyme in solution and the adenosine A 2A receptor in a membrane. The calculated free energy values are overall consistent with the experimental and published free energy results (within ~1 kcal/mol). We hope that Free Energy Calculator is useful to carry out high-throughput FEP/MD simulations in the field of biomolecular sciences and drug discovery.

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

Dynamic Molecular Investigation of the Solid-Electrolyte Interphase of an Anode-Free Lithium Metal Battery Using In Situ Liquid SIMS and Cryo-TEM

Solid electrolyte interphase (SEI) has been widely perceived to play a critical role in the stable cycling of rechargeable batteries. However, associated with the fragile and air-sensitive nature of the SEI layer, delineation of the formation process and the nature of SEI remains a big challenge. Here, we use in situ liquid time-of-flight secondary ion mass spectroscopy (TOF-SIMS), cryo- transmission electron microscope (TEM) and density functional theory (DFT) calculation to delineate molecular process on the formation of SEI layer under the dynamic operating condition. We discover that the onset potential for SEI layer formation and the thickness of the SEI show dependence on the solvation shell structure. Using LiCoO 2 as a cathode and Cu film as an anode, the SEI is noticed to start to form at around 2.0 V and reach its final thickness (irreversible part, ~ 40-50 nm) at about 3.0 V in the 1 M LiPF 6 –EC/DMC electrolyte, while for the case of 1 M LiFSI–DME, the SEI starts to form at around 1.5 V and reaches its final thickness (~ 20 nm) at about 2.0 V. The in situ TOF-SIMS clearly indicates the outer SEI layer formation and dissipation upon charging and discharging, implying a continued evolution of electrolyte structure with extended cycling. In conclusion, the present work establishes a direct correlation between the molecular signature of SEI layer with solvation feature of electrolytes in lithium batteries, providing insights for tailoring SEI layer toward improved electrochemical properties of lithium batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

A compartmentalized model of multiphase chemical kinetics

There are significant challenges in predicting multiphase chemical kinetics due to the complex coupling of reaction and mass transport across a phase boundary (i.e., interface). Here, we describe a framework for predicting multiphase kinetics that embeds the elementary kinetic steps of reaction, solvation, and diffusion into a coarse grain spatial description of two phases. The model is constructed to bridge the short-timescale interfacial dynamics observed in molecular simulations with the longer timescales observed in kinetic experiments. A simple set of governing differential equations is derived, which, when solved numerically or analytically, yield accurate predictions of multiphase kinetics in microdroplets. Although the equations are formulated for gas-liquid reactions, the underlying conceptual framework is general and can be applied to transformations in other two-phase systems (solid-liquid, liquid-liquid, etc.).

Chemical kinetics and dynamics↗
DOE OSTI2021

Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2

The association of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with human angiotensin-converting enzyme 2 (hACE2) represents the first required step for cellular entry. SARS-CoV-2 has continued to evolve with the emergence of several novel variants, and amino acid changes in the RBD have been implicated with increased fitness and potential for immune evasion. Reliably predicting the effect of amino acid changes on the ability of the RBD to interact more strongly with the hACE2 can help assess the implications for public health and the potential for spillover and adaptation into other animals. Here, we introduce a two-step framework that first relies on 48 independent 4-ns molecular dynamics (MD) trajectories of RBD-hACE2 variants to collect binding energy terms decomposed into Coulombic, covalent, van der Waals, lipophilic, generalized Born solvation, hydrogen bonding, π-π packing, and self-contact correction terms. The second step implements a neural network to classify and quantitatively predict binding affinity changes using the decomposed energy terms as descriptors. The computational base achieves a validation accuracy of 82.8% for classifying single–amino acid substitution variants of the RBD as worsening or improving binding affinity for hACE2 and a correlation coefficient of 0.73 between predicted and experimentally calculated changes in binding affinities. Both metrics are calculated using a fivefold cross-validation test. Our method thus sets up a framework for screening binding affinity changes caused by unknown single– and multiple–amino acid changes offering a valuable tool to predict host adaptation of SARS-CoV-2 variants toward tighter hACE2 binding.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

The Interplay of Solvation and Polarization Effects on Ion Pairing in Nanoconfined Electrolytes

The nature of ion–ion interactions in electrolytes confined to nanoscale pores has important implications for energy storage and separation technologies. However, the physical effects dictating the structure of nanoconfined electrolytes remain debated. Here we employ machine-learning-based molecular dynamics simulations to investigate ion–ion interactions with density functional theory level accuracy in a prototypical confined electrolyte, aqueous NaCl within graphene slit pores. We find that the free energy of ion pairing in highly confined electrolytes deviates substantially from that in bulk solutions, observing a decrease in contact ion pairing but an increase in solvent-separated ion pairing. These changes arise from an interplay of ion solvation effects and graphene’s electronic structure. Notably, the behavior observed from our first-principles-level simulations is not reproduced even qualitatively with the classical force fields conventionally used to model these systems. The insight provided in this work opens new avenues for predicting and controlling the structure of nanoconfined electrolytes.

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

Dynamic Interconversion of Metal Active Site Ensembles in Zeolite Catalysis

Catalysis science is founded on understanding the structure, number, and reactivity of active sites. Kinetic models that consider active sites to be static and noninteracting entities are routinely successful in describing the behavior of heterogeneous catalysts. Yet, active site ensembles often restructure in response to their external environment and even during steady-state catalytic turnover, sometimes requiring non-mean-field kinetic treatments to describe distance-dependent interactions among sites. Such behavior is being recognized more frequently in modern catalysis research, with the advent of experimental methods to quantify turnover rates with increasing precision, an expanding arsenal of operando characterization tools, and computational descriptions of atomic structure and motion at chemical potentials and timescales increasingly relevant to reaction conditions. This review focuses on dynamic changes to metal active site ensembles on zeolite supports, which are silica-based crystalline materials substituted with Al that generate binding sites for isolated and low-nuclearity metal site ensembles. Metal sites can become solvated and mobilized during reaction, facilitating interactions among sites that change their nuclearity and function. Such intersite communication can be regulated by the zeolite support, resulting in non-single-site and potentially non-mean-field kinetic behavior arising from mechanisms of catalytic action that combine elements of those canonically associated with homogeneous and heterogeneous catalysis. Here, we discuss recent literature examples that document dynamic active site behavior in metal-zeolites and outline methodologies to identify and interpret such behavior. We conclude with our outlook on future research directions to develop this evolving branch of catalysis science and harness it for practical applications.

operando characterization↗
DOE OSTI2021

Confinement effects and acid strength in zeolites

Chemical reactivity and sorption in zeolites are coupled to confinement and—to a lesser extent—to the acid strength of Brønsted acid sites (BAS). In presence of water the zeolite Brønsted acid sites eventually convert into hydronium ions. The gradual transition from zeolite Brønsted acid sites to hydronium ions in zeolites of varying pore size is examined by ab initio molecular dynamics combined with enhanced sampling based on Well-Tempered Metadynamics and a recently developed set of collective variables. While at low water content (1–2 water/BAS) the acidic protons prefer to be shared between zeolites and water, higher water contents (n > 2) invariably lead to solvation of the protons within a localized water cluster adjacent to the BAS. At low water loadings the standard free energy of the formed complexes is dominated by enthalpy and is associated with the acid strength of the BAS and the space around the site. Conversely, the entropy increases linearly with the concentration of waters in the pores, favors proton solvation and is independent of the pore size/shape.

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

Localized high-concentration electrolytes get more localized through micelle-like structures

Liquid electrolytes in batteries are typically treated as macroscopically homogeneous ionic transport media despite having a complex chemical composition and atomistic solvation structures, leaving a knowledge gap of the microstructural characteristics. Here, we reveal a unique micelle-like structure in a localized high-concentration electrolyte, in which the solvent acts as a surfactant between an insoluble salt in a diluent. The miscibility of the solvent with the diluent and simultaneous solubility of the salt results in a micelle-like structure with a smeared interface and an increased salt concentration at the centre of the salt–solvent clusters that extends the salt solubility. These intermingling miscibility effects have temperature dependencies, wherein a typical localized high-concentration electrolyte peaks in localized cluster salt concentration near room temperature and is used to form a stable solid–electrolyte interphase on a Li metal anode. Importantly, these findings serve as a guide to predicting a stable ternary phase diagram and connecting the electrolyte microstructure with electrolyte formulation and formation protocols of solid–electrolyte interphases for enhanced battery cyclability.

25 ENERGY STORAGE↗
DOE OSTI2020

Charging-driven coarsening and melting of a colloidal nanoparticle monolayer at an ionic liquid–vacuum interface

In this work, we induce and investigate the coarsening and melting dynamics of an initially static nanoparticle colloidal monolayer at an ionic liquid–vacuum interface, driven by a focused, scanning electron beam. Coarsening occurs through grain interface migration and larger-scale motions such as grain rotations, often facilitated by sliding dislocations. The progressive decrease in area fraction that drives melting of the monolayer is explained using an electrowetting model whereby particles at the interface are solvated once their accumulating charge recruits sufficient counterions to subsume the particle. Subject to stochastic particle removal from the monolayer, melting is recapitulated in simulations with a Lennard-Jones potential. This new driving mechanism for colloidal systems, whose dynamical timescales we show can be controlled with the accelerating voltage, opens the possibility to manipulate particle interactions dynamically without need to vary particle intrinsic properties or surface treatments. Furthermore, the decrease in particle size availed by electron imaging presents opportunities to observe force and time scales in a lesser-explored regime intermediate between typical colloidal and molecular systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ion Pairing, Clustering and Transport in a LiFSI-TMP Electrolyte as Functions of Salt Concentration using Molecular Dynamics Simulations

Battery capacity is highly related to ion-pairing mechanisms in electrolytes, since a cluster formation can lead to dead Li formation, reducing the number of charge carriers and leading to capacity fading. We use molecular dynamics simulations to model an electrolyte comprising trimethyl phosphate (TMP) solvent and a lithium bis(fluorosulfonyl)imide (LiFSI) salt, exploring effects of salt concentration on solvation and ion-transport. We simulate the LiFSI-TMP electrolyte for salt concentrations of 0.7, 1.43 and 3.82 molar. A statistical analysis was performed to study ion-pairing, clustering, diffusivity, conductivity, and coordination of Li-ions, providing insights into relations between molecular structures and transport properties. Molecular structure of ionic components changes as concentration increases, from a predominant solvent separated ion pair (SSIP) and contact ion pair (CIP) to aggregate salt (AGG) and ionic cluster formation. Given the formation of the ionic cluster, the diffusion mechanism followed by Li-ions changes from a hopping/exchange to a vehicular mechanism as concentration increases; this is reflected in a decrease of ionic conductivities. Ionicity was also calculated to reveal how the ionic motion changes from an uncorrelated to a correlated one as the salt concentration increases. Furthermore, we also compared our results with experimental calculations performed for similar electrolyte systems

25 ENERGY STORAGE↗
DOE OSTI2023

Probing pH-Dependent Dehydration Dynamics of Mg and Ca Cations in Aqueous Solutions with Multi-Level Quantum Mechanics/Molecular Dynamics Simulations

The dehydration of aqueous calcium and magnesium cations is the most fundamental process controlling their reactivity in chemical and biological phenomena, such as the formation of ionic solids or passing through ion channels. It holds particular relevance in light of recent advancements in the development of carbon capture techniques that rely on mineralization for long-term carbon storage. Specifically, dehydration of Ca 2+ and Mg 2+ is a key step in proposed carbon capture processes aiming to exploit the relatively high concentration of dissolved carbon dioxide in seawater via the formation of carbonate minerals from solvated Ca 2+ and Mg 2+ cations for sequestration and storage. Nevertheless, atomic-scale understanding of the dehydration of aqueous Ca 2+ and Mg 2+ cations remains limited. Here, we utilize rare event sampling via density functional theory molecular dynamics and embedded wavefunction theory calculations to elucidate the dehydration dynamics of aqueous Ca 2+ and Mg 2+ . Emphasis is placed on the investigation of the effect pH has on the stability of the different coordination environments. Here, our results reveal significant differences in the dehydration dynamics of the two cations and provide insight into how they may be modulated by pH changes.

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

Tailoring Na⁺ Chelation Dynamics for Expedient Sulfur Redox Kinetics in Low‐Temperature Sodium–Sulfur Batteries

Sodium–sulfur (Na–S) batteries have attracted considerable attention due to their high theoretical energy density and the abundant natural availability of sodium and sulfur. However, sluggish kinetics of sulfur conversion, slow Na⁺ transport, and interfacial instability at low temperatures pose significant challenges for their operation and limit their practical application. Herein, three solvents with well-designed molecular configurations are examined. We systematically investigate the impact of chelation effect on the desolvation behavior, sulfur conversion process, ion dynamics, and Na⁺ plating/stripping behavior. Compared with conventional linear ether solvents, the incorporation of methyl groups not only weaken the chelation capability and tailors the inner solvation sheath, but also reduces the energy barrier for Na⁺ transport, thus promoting enhanced sulfur conversion kinetics under low temperature conditions. This work elucidates the relationship among solvent molecules, Na⁺ desolvation behavior, and sulfur reaction kinetics, and offers a strategy for rational design of electrolytes for low-temperature metal–sulfur batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

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

Persistence of Methionine Side Chain Mobility at Low Temperatures in a Nine‐Residue Low Complexity Peptide, as Probed by 2 H Solid‐State NMR

Abstract Methionine side chains are flexible entities which play important roles in defining hydrophobic interfaces. We utilize deuterium static solid‐state NMR to assess rotameric inter‐conversions and other dynamic modes of the methionine in the context of a nine‐residue random‐coil peptide (RC9) with the low‐complexity sequence GGKGMGFGL. The measurements in the temperature range of 313 to 161 K demonstrate that the rotameric interconversions in the hydrated solid powder state persist to temperatures below 200 K. Removal of solvation significantly reduces the rate of the rotameric motions. We employed 2 H NMR line shape analysis, longitudinal and rotation frame relaxation, and chemical exchange saturation transfer methods and found that the combination of multiple techniques creates a significantly more refined model in comparison with a single technique. Further, we compare the most essential features of the dynamics in RC9 to two different methionine‐containing systems, characterized previously. Namely, the M35 of hydrated amyloid‐β 1–40 in the three‐fold symmetric polymorph as well as Fluorenylmethyloxycarbonyl (FMOC)‐methionine amino acid with the bulky hydrophobic group. The comparison suggests that the driving force for the enhanced methionine side chain mobility in RC9 is the thermodynamic factor stemming from distributions of rotameric populations, rather than the increase in the rate constant.

Vugmeyster, Liliya↗