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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Probing Condensed-Phase Structure and Dynamics in Hierarchical Zeolites and Nanosheets for Catalytic Upgradation of Biomass (Final Report)

Understanding complex reaction pathways in systems governed by multi-scale collective interactions across time and length scales remains a central scientific challenge. This project was guided by the hypothesis that the interplay among oligomers, solvents, and active sites can be tuned by a suitable choice of solvation environment and pore architecture in solid-acid catalysts to direct chemical transformations relevant to biomass conversion. Zeolites and zeolite nanosheets were used as model platforms, allowing for the interaction of macromolecules with the surface of the zeolite nanosheets and with smaller pores that host catalytically active sites. To investigate these coupled phenomena, we employ a multi-scale computational framework that integrates molecular-level descriptions with advanced sampling approaches to capture key physical and chemical interactions. Our work through this project improved fundamental understanding of how reactants and solid-acid catalysts interact in solvent-rich environments, thereby enabling the rational design of catalytic systems that upgrade biomass with enhanced selectivity and energy efficiency. In addition, the project developed advanced sampling methodologies critical for disentangling complex, reactive processes in multi-component catalytic environments.

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

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗
DOE OSTI2026

Cooperative Lanthanide Solvation in an Ionic Liquid for Critical Materials Separations

Effectively separating critical materials, including intra-lanthanide separations, is crucial for meeting growing application demands. Liquid–liquid extraction (LLE) is the industry standard for lanthanide separations, where the selectivity can depend on small changes in metal coordination. In this work, we investigate representative lanthanide Eu coordination with a neutral malonamide extractant in an imidazolium bistriflimide ionic liquid (IL) solvent. Through systematic titrations of the water and extractant and under extraction conditions, we observe surprising cooperative Eu solvation with the IL anion and extractant. Time-resolved fluorescence spectroscopy measurements show strong extractant coordination in water-saturated IL. Lifetime measurements show no water coordination, and extended X-ray fine structure spectroscopy data provide a coordination number of 10. Molecular dynamics simulations confirm this coordination number and reveal IL anion coordination in the final Eu complex, even though it is ordinarily a significantly weaker ligand compared to water. The lack of water in the final extracted complex and IL anion coordination potentially explain the increased extraction in LLE systems using ILs, as evidenced by higher distribution ratios for cation exchange extraction, despite the energetic cost of cation transfer to the aqueous phase. These results highlight the opportunities for tuning metal coordination to drive extraction in unique solvent systems.

Carr, Amanda J.↗
DOE OSTI2022

Dissolving salt is not equivalent to applying a pressure on water

Abstract Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water’s structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

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

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗
DOE OSTI2020

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

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

The Influence of Alkyl Spacers and Molecular Weight on the Charge Transport and Storage Properties of Oxy‐Bithiophene‐Based Conjugated Polymers

Conjugated polymers (CPs) with polar side chains can conduct electronic and ionic charges simultaneously, making them promising for bioelectronics, electrocatalysis and energy storage. Recent work showed that adding alkyl spacers between CP backbones and polar side chains improved electronic charge carrier mobility, reduced swelling and enhanced stability, without compromising ion transport. However, how alkyl spacers impact polymer backbone conformation and, subsequently, electronic properties remain unclear. In this work, we design two oxy-bithiophene-based CP series, each featuring progressively extended alkyl spacer lengths and two distinct molecular weight (MW) distributions. Using operando characterisations, we evaluate the (spectro)electrochemical and swelling properties of the polymer thin films, and their performance in organic field-effect transistors and organic electrochemical transistors. Surprisingly, alkyl spacers negatively impact the hole mobility of our polymers, with higher MW amplifying this effect. Using molecular dynamics simulations, we show that it is thermodynamically favourable for adjacent non-polar alkyl spacers to aggregate in polar electrolytes, leading to backbone twisting. Further spectroscopic measurements corroborate this prediction. Our findings demonstrate the active interactions between side chain structure, MW and electrolyte/solvent polarity in influencing polymer performance, underscoring the importance of considering solvation environment effects on polymer conformation when designing new mixed conducting CPs for electrochemical applications.

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

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗
DOE OSTI2026

Nanoconfined Grain Boundaries Increase the Conductivity of Polycrystalline Molecular Crystals

Soft-solid molecular crystals consist of crystalline grains and fluid grain boundaries (GBs) that enhance the grain binding and transport of Li + ions between the grains. The total ionic conductivity consists of ion migration in both the grains and GBs. To unravel these contributions in adiponitrile (Adpn):LiPF 6 molecular crystals, the GB volume fraction was varied by changing the size of the crystals and the Adpn:LiPF 6 molar ratio. Molecular dynamics (MD) simulations indicate that ion motion was subdiffusive in the grains and “well-diffusive” in the GBs, with GBs characterized as disordered nanoconfined regions of higher charge carrier concentration (∼1 M) than in saturated Adpn:LiPF 6 solutions (0.04 M), and Li + ions predominantly solvated by cyano groups with few contact ion pairs. The diffusivity in the GBs is at least an order of magnitude higher than that in the crystalline grains. The emergent picture is the grains as a reservoir of ions that migrate to faster-conducting GBs.

Crystals↗
DOE OSTI2026

Strong Effect of Nonpolar Solvent Molecular Structure on CdSe Nanoplatelet Stacking

We report a drastic difference in stacking behavior of oleic acid-stabilized 4-monolayer (4 ML) CdSe nanoplatelets (NPLs) in toluene and methylcyclohexane (MCH), two nonpolar solvents that differ in the conformational flexibility of their molecules. Using liquid cell transmission electron microscopy (TEM) and small angle scattering (SAXS) techniques, we show that NPLs form microns-long ribbons consisting of 4 ML CdSe NPLs in toluene, the solvent widely used to form stable colloidal solutions of a broad range of quasi-spherical nanoparticles. In contrast, 4 ML CdSe NPLs are well dispersed in MCH. The difference in stacking behavior of NPLs in toluene and MCH suggests that the conformational flexibility of the solvent molecules, such as the ability to adopt multiple chair conformations, modulates nanoplatelet interactions. Molecular dynamics (MD) simulations reveal that solvent molecules subtly alter the structure of the organic ligand shell. These solvent-dependent changes propagate to the inorganic core, modulating the degree of CdSe nanoplatelet (NPL) twisting and, consequently, the properties of the nanoparticles. We show that toluene better solvates oleate ligands while MCH induces a bimodal oleate span distribution, which can lead to increased solubility of CdSe NPLs. In addition, the solvent can also influence the inorganic core, which, in turn, can modify the nanoparticle properties. We demonstrate that destabilization of toluene solution containing ribbons of 4 ML CdSe NPLs without CdS shells results in the formation of NPL assemblies with amplified spontaneous emission (ASE) with a low threshold of 14 µJ cm−2 that is comparable with that of CdSe/CdS core/shell NPLs. Our results emphasize that the solvent plays a major role in mediating interactions between NPLs and hence their processability for fabrication of functional structures.

CdSe↗
DOE OSTI2021

Exploring the Ion Solvation Environments in Solid-State Polymer Electrolytes through Free-Energy Sampling

The success of poly(ethylene oxide) (PEO) in solid-state polymer electrolytes for lithium-ion batteries is well established. Recently, in order to understand this success and to explore possible alternatives, we studied polyacetal electrolytes to deepen the understanding of the effect of the local chemical structure on ion transport. Advanced molecular dynamics techniques using newly developed, tailored interaction potentials have helped elucidate the various coordination environments of ions in these systems. In particular, the competition between cation-Anion pairing and coordination by the polymer has been explored using free-energy sampling (metadynamics). At equivalent reduced temperatures, with respect to the polymer-specific glass-Transition temperature, two-dimensional free-energy plots reveal the existence of multiple coordination environments for the lithium (Li) ions in these systems and their relative stabilities. Furthermore, we observe that the Li-ion movement in PEO follows a serial, stepwise pathway when moving from one coordination state to another, whereas this happens in a more continuous and concerted fashion in a polyacetal such as poly(1,3-dioxalane) [P(EO-MO)]. The implication is that interconversion between coordination states of the Li ions may be easier in P(EO-MO). However, the overarching observation from our free-energy analysis is that Li-ion coordination is dominated by the polymer (in either case) and contact-ion pairs are rare. We rationalize the observed higher increase in glass-Transition temperature (Tg) with salt loading in polyacetals as due to intermolecular Li-ion coordination involving multiple polymer chains, rather than just one chain for PEO-based electrolytes. This interchain coupling in the polyacetals, resulting in the higher Tg, works against any gains due to variations in Li-ion coordination that might enhance transport processes over PEO. Further research is required to overcome the interdependence between local coordination and macroscopic properties to compete with PEO electrolytes at the same absolute working temperature.

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

Partial Solvation of Lithium Ions Enhances Conductivity in a Nanophase-Separated Polymer Electrolyte

In this study, we demonstrate that a multiblock lithium-ion-conducting polymer can be swollen with ethylene carbonate solvent to increase the conductivity relative to the dry polymer material by nearly 4 orders of magnitude. This increase is due to the partial solvation of lithium ions by ethylene carbonate, which leads to Li + diffusion along the solvent–polymer interface. This differs from the vehicular transport mechanism for lithium ions in pure solvent. We use a combination of broadband dielectric spectroscopy, X-ray scattering, and all-atom molecular dynamics simulations to probe the effect of the solvent on the polymer morphology and to elucidate the mechanism of lithium ion transport.

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

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Stabilization of Hydroxide Ions at the Interface of a Hydrophobic Monolayer on Water via Reduced Proton Transfer

We report a joint study using surface-specific sum-frequency vibrational spectroscopy and ab initio molecular dynamics simulations, respectively, on a pristine hydrophobic (sub)monolayer hexane-water interface, namely, the hexane/water interface with varied vapor pressures of hexane and different pHs in water. We show clear evidence that hexane on water revises the interfacial water structure in a way that stabilizes the hypercoordinated solvation structure and slows down the migration of hydroxide ion (OH – ) relative to that in bulk water. Furthermore, this mechanism effectively attracts the OH – to the water-hydrophobic interface with respect to its counterion. The result illustrates the striking difference of proton transfer of hydrated OH – at the interface and in the bulk, which is responsible for the intrinsic charging effect at the hydrophobic interface.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Behavior of Water Near Multimodal Chromatography Ligands and Its Consequences for Modulating Protein–Ligand Interactions

Multimodal chromatography is a powerful approach for purifying proteins that uses ligands containing multiple modes of interaction. Recent studies have shown that selectivity in multimodal chromatographic separations is a function of the ligand structure and geometry. Here, we performed molecular dynamics simulations to explore how the ligand structure and geometry affect ligand–water interactions and how these differences in solution affect the nature of protein–ligand interactions. Our investigation focused on three chromatography ligands: Capto MMC, Nuvia cPrime, and Prototype 4, a structural variant of Nuvia cPrime. First, the solvation characteristics of each ligand were quantified via three metrics: average water density, fluctuations, and residence time. We then explored how solvation was perturbed when the ligand was bound to the protein surface and found that the probability of the phenyl ring dewetting followed the order: Capto MMC > Prototype 4 > Nuvia cPrime. To explore how these differences in dewetting affect protein–ligand interactions, we calculated the probability of each ligand binding to different types of residues on the protein surface and found that the probability of binding to a hydrophobic residue followed the same order as the dewetting behavior. This study illustrates the role that wetting and dewetting play in modulating protein–ligand interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data provided valuable insights into another type of non-DLVO interaction, namely the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. In conclusion, our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

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

Compact Peptoid Molecular Brushes for Nanoparticle Stabilization

Controlling the interfaces and interactions of colloidal nanoparticles (NPs) via tethered molecular moieties is crucial for NP applications in engineered nanomaterials, optics, catalysis, and nanomedicine. Despite a broad range of molecular types explored, there is a need for a flexible approach to rationally vary the chemistry and structure of these interfacial molecules for controlling NP stability in diverse environments, while maintaining a small size of the NP molecular shell. Here, we demonstrate that low-molecular-weight, bifunctional comb-shaped, and sequence-defined peptoids can effectively stabilize gold NPs (AuNPs). The generality of this robust functionalization strategy was also demonstrated by coating of silver, platinum, and iron oxide NPs with designed peptoids. Each peptoid (PE) is designed with varied arrangements of a multivalent AuNP-binding domain and a solvation domain consisting of oligo-ethylene glycol (EG) branches. Among designs, a peptoid (PE5) with a diblock structure is demonstrated to provide a superior nanocolloidal stability in diverse aqueous solutions while forming a compact shell (similar to 1.5 nm) on the AuNP surface. We demonstrate by experiments and molecular dynamics simulations that PE5-coated AuNPs (PE5/AuNPs) are stable in select organic solvents owing to the strong PE5 (amine)-Au binding and solubility of the oligo-EG motifs. At the vapor-aqueous interface, we show that PE5/AuNPs remain stable and can self-assemble into ordered 2D lattices. The NP films exhibit strong near-field plasmonic coupling when transferred to solid substrates.

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

Alkali Cation Inhibition of Imidazolium-Mediated Electrochemical CO 2 Reduction on Silver

Imidazolium-based ionic liquids have led to enhanced CO 2 electroreduction activity due to cation effects at the cathode surface, stabilizing the reaction intermediates and decreasing the activation energy. In aqueous media, alkali cations are also known to improve CO 2 reduction activity on metals such as Ag, with the enhancement attributed to electrical double layer effects and trending with the size of the alkali cation. However, the effect of a mixed catholyte solution of alkali cations in the presence of an imidazolium-based ionic liquid has not been well-explored. Herein, 1-ethyl-3-methylimidazolium tetrafluoroborate, [EMIM][BF 4 ], in water was investigated with alkali salts to unravel the interaction effects for CO 2 electroreduction on Ag. Although both [EMIM] + and alkali cations have individually improved CO 2 to CO conversion on Ag in water, electrochemical results showed that alkali cations hindered imidazolium-mediated CO 2 electroreduction in most conditions. Li + , in particular, was sharply inhibitory compared to other alkali cations and strongly redirected the selectivity to hydrogen evolution. The nature of the alkali cation inhibition was investigated with spectroscopic techniques, including in situ surface-enhanced Raman spectroscopy (SERS) and dynamic electrochemical impedance spectroscopy (DEIS). Along with computational insights from density functional theory (DFT), the electrochemical and spectroscopic data suggest that alkali cations inhibit [EMIM]-mediated CO 2 reduction by competing for surface adsorption sites, preventing the potential-dependent structural reorientation of imidazolium, and promoting hydrogen evolution by bringing solvated water to the cathode surface.

cations↗