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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of Salt Aggregation in Perfluoroether Electrolytes

Electrolytes comprised of polymers mixed with salts have great potential for enabling the use of Li metal anodes in batteries for increased safety. Ionic conductivity is one of the key performance metrics of these polymer electrolytes and achieving high room-temperature conductivity remains a challenge to date. For a bottom-up design of the polymer electrolytes, we must first understand how the structure of polyelectrolytes on a molecular level determines their properties. Here, we use classical molecular dynamics to study the solvation structure and ion diffusion in electrolytes composed of a short-chain perfluoroether with LiFSI or LiTFSI salts. Density functional theory is also used to provide some insights into the structures and energies of the salt interactions with the perfluoroether. We observe the formation of aggregates of salts in the fluorinated systems even at low salt concentrations. The fluorine-fluorine attraction in the solvent is the governing factor for creating the salt aggregates. The aggregates’ size and lifetime change with concentration and anion. These simulations provide an insight into the structure and dynamics of perfluoroether based electrolytes that can be used to improve Li-ion batteries.

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

Solvent structure and dynamics over Brønsted acid MWW zeolite nanosheets

In the liquid phase of heterogeneous catalysis, solvent plays an important role and governs the kinetics and thermodynamics of a reaction. Although it is often difficult to quantify the role of the solvent, it becomes particularly challenging when a zeolite is used as the catalyst. This difficulty arises from the complex nature of the liquid/zeolite interface and the different solvation environments around catalytically active sites. Here, we use ab initio molecular dynamics simulations to probe the local solvation structure and dynamics of methanol and water over MWW zeolite nanosheets with varying Brønsted acidity. We find that the zeolite framework and the number and location of the acid sites in the zeolite influence the structure and dynamics of the solvent. In particular, methanol is more likely to be in the vicinity of the aluminum (Al 3+ ) at the T4 site than at T1 due to easy accessibility. The methanol oxygen binds strongly to the Al at the T4 site, weakening the Al–O for the bridging acid site, which results in the formation of the silanol group, significantly reducing the acidity of the site. The behavior of methanol is in direct contrast to that of water, where protons can easily propagate from the zeolite to the solvent molecules regardless of the acid site location. Furthermore, our work provides molecular-level insights into how solvent interacts with zeolite surfaces, leading to an improved understanding of the catalytic site in the MWW zeolite nanosheet.

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

Dynamic Copper Site Redispersion through Atom Trapping in Zeolite Defects

Single-site copper-based catalysts have shown remarkable activity and selectivity for a variety of reactions. However, deactivation by sintering in high-temperature reducing environments remains a challenge and often limits their use due to irreversible structural changes to the catalyst. Here, we report zeolite-based copper catalysts in which copper oxide agglomerates formed after reaction can be repeatedly redispersed back to single sites using an oxidative treatment in air at 550 °C. Under different environments, single-site copper in Cu–Zn–Y/deAlBeta undergoes dynamic changes in structure and oxidation state that can be tuned to promote the formation of key active sites while minimizing deactivation through Cu sintering. For example, single-site Cu 2+ reduces to Cu 1+ after catalyst pretreatment (270 °C, 101 kPa H 2 ) and further to Cu0 nanoparticles under reaction conditions (270–350 °C, 7 kPa EtOH, 94 kPa H 2 ) or accelerated aging (400–450 °C, 101 kPa H 2 ). After regeneration at 550 °C in air, agglomerated CuO was dispersed back to single sites in the presence and absence of Zn and Y, which was verified by imaging, in situ spectroscopy, and catalytic rate measurements. Ab initio molecular dynamics simulations show that solvation of CuO monomers by water facilitates their transport through the zeolite pore, and condensation of the CuO monomer with a fully protonated silanol nest entraps copper and reforms the single-site structure. Importantly, the capability of silanol nests to trap and stabilize copper single sites under oxidizing conditions could extend the use of single-site copper catalysts to a wider variety of reactions and allows for a simple regeneration strategy for copper single-site catalysts.

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

Rational Design of Weakly‐Solvating Molecules for Salt‐In‐Pre‐Ionic‐Liquid Electrolytes for Li Metal Batteries

Lithium metal batteries (LMBs) promise step‐changes in energy densities but suffer from poor cycle life due to unstable electrolyte‐lithium interfaces. Conventional carbonate electrolytes exhibit excessive lithium‐ion solvation and low oxidative stability, leading to rapid capacity loss. Herein, we report a rationally designed weakly‐solvating cyclic sulfonamide, 1‐trifluoromethanesulfonyl)amide pyrrolidine (TFMSPyr), which integrates an electron‐withdrawing trifluoromethanesulfonyl functional group at pyrrolidinic‐N. TFMSPyr acts as a pre‐ionic‐liquid solvent that forms intrinsically localized, anion‐dominated solvation, coupling molecular architecture, solvation topology, and transport dynamics. As a result, LiFSI based salt‐in‐pre‐ionic‐liquid (SIPIL) electrolytes exhibit high lithium‐ion transference number, oxidative stability > 5 V versus Li/Li + and anion‐derived solid electrolyte interphases (SEI). Li||Cu cells with SIPIL deliver a first cycle Coulombic efficiency (CE) of ≈ 99% with average CE of 99.2% for 100 cycles, and lithium half‐cells with lithium iron phosphate (LFP) cathode exhibit 82% capacity retention after 400 cycles with CE of 99.98%. In anode‐free full cells, 95% of initial capacity is retained after 63 cycles with an average CE of 99.5%. These results demonstrate that molecular engineering of solvents offers a powerful pathway to stabilize lithium metal interfaces and enable practical Anodeless LMBs.

25 ENERGY STORAGE↗
DOE OSTI2024

On the Statistical Mechanics of Mass Accommodation at Liquid–Vapor Interfaces

Here we propose a framework for describing the dynamics associated with the adsorption of small molecules to liquid-vapor interfaces using an intermediate resolution between traditional continuum theories that are bereft of molecular detail and molecular dynamics simulations that are replete with them. In particular, we develop an effective single particle equation of motion capable of describing the physical processes that determine thermal and mass accommodation probabilities. The effective equation is parametrized with quantities that vary through space away from the liquid-vapor interface. Of particular importance in describing the early time dynamics is the spatially dependent friction, for which we propose a numerical scheme to evaluate from molecular simulation. Taken together with potentials of mean force computable with importance sampling methods, we illustrate how to compute the mass accommodation coefficient and residence time distribution. Throughout, we highlight the case of ozone adsorption in aqueous solutions and its dependence on electrolyte composition.

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

Ab Initio Molecular Dynamics Study of the Reduction of Acetone by the Hydrated Electron

We have investigated the reaction dynamics of the reduction of acetone by the solvated electron in water using ab initio molecular dynamics simulations at 298 and 373 K. The rate constants derived from the simulations are consistent with experimental observations that the reaction has a low activation energy. Detailed analyses carried out to shed light on the mechanism of the reaction show that solvent reorganization plays a key role as a reaction coordinate, as expected from the Marcus electron transfer theory. Furthermore, constrained density functional theory calculations indicate that the electronic coupling is large, placing the reaction in the adiabatic limit. Indeed, the activation energies and rate constants of the simulations are in accordance with the predictions of adiabatic Marcus theory.

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

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2021

Mass action model of solution activity via speciation by solvation and ion pairing equilibria

Abstract Solutes and their concentrations influence many natural and anthropogenic solution processes. Electrolyte and solution models are used to quantify and predict such behavior. Here we present a mechanistic solution model based on mass action equilibria. Solvation and ion pairing are used to model speciated solute and solvent concentrations such that they correlate to a solution’s vapor pressure (solvent activity) according to Raoult’s law from dilute conditions to saturation. This model introduces a hydration equilibrium constant ( K ha ) that is used with either an ion dissociation constant ( K id ) or a hydration modifier ( m ) with an experimentally determined ion dissociation constant, as adjustable parameters to fit vapor–liquid equilibrium data. The modeled solvation equilibria are accompanied by molecular dynamics (MD) studies that support a decline in the observed degree of solvation with increased concentration. MD calculations indicate this finding is a combination of a solvent that solvates multiple solutes, and changes in a solute’s solvation sphere, with the dominant factor changing with concentration. This speciation-based solution model is lateral to established electrostatics-based electrolyte theories. With its basis in mass action, the model can directly relate experimental data to the modeled solute and solvent speciated concentrations and structures.

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

Polymer principles behind solubilizing lignin with organic cosolvents for bioenergy

Lignin solubilization is key to a viable biorefinery because its removal leads to facile deconstruction of biomass and because the isolated lignin can serve to derive precursors of novel high-value materials. The mixing of organic solvents with water has been shown to improve biomass fractionation and lignin conversion reactions. However, generally-applicable solubilization strategies are lacking because of the remarkable variability of lignin across plant feedstocks. Here, to obtain a predictive understanding of lignin solvation, we perform molecular dynamics simulations of model lignin polymers in two mixtures of water with polar aprotic solvents: tetrahydrofuran (THF) : water and γ-valerolactone (GVL) : water. The model lignins include H-, G- and S-only homopolymers and lignins with an S : G 1 : 1 ratio. We find that a well-established theory of self-avoiding polymers in a “good” solvent describes accurately the physical conformations of all types of lignin in both solvents. As the degree of methoxy substitution increases in the homopolymers, the distributions of the lignin radius of gyration and the Flory exponent ν , which describes the lignin-solvent interactions, do not change in THF : water, while ν shifts to slightly higher values in GVL : water. We attribute this increase to the interaction between the methyl group of GVL with the lignin methoxy groups. We also find that the reduction in the lignin radius of gyration due to branching is accurately described by the Zimm–Stockmayer theory for both THF : water and GVL : water. The above findings validate the applicability of polymer physics concepts to lignin and suggest that GVL : water may have the most favorable interaction with S-lignin, whereas the interactions of THF : water with lignin are independent of lignin monomeric content.

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

Effect of Zwitterionic Additives on Solvation and Transport of Sodium and Potassium Cations in (Ethylene Oxide)10: A Molecular Dynamics Simulation Study

Sodium- (Na + ) and potassium- (K + ) ion batteries are cost-effective alternatives to lithium-ion (Li + ) batteries due to the abundant sodium and potassium resources. Solid polymer electrolytes (SPEs) are essential for safer and more efficient Na + and K + batteries because they often exhibit low ionic conductivity at room temperature. While zwitterionic (ZW) materials enhance Li+ battery conductivity, their potential for Na + and K + transport in batteries remains unexplored. In this study, we investigated the effect of three ZW molecules (ChoPO4, i.e., 2-methacryloyloxyethyl phosphorylcholine, ImSO3, i.e., sulfobetaine ethylimidazole, and ImCO2, i.e., carboxybetaine ethylimidazole) on the dissociation of Na + and K + coordination with ethylene oxide (EO) chains in EO-based electrolytes through molecular dynamics simulations. Our results showed that ChoPO4 possessed the highest cation–EO10 dissociation ability, while ImSO3 exhibited the lowest. Such dissociation ability correlated with the cation–ZW molecule coordination strength: ChoPO4 and ImSO3 showed the strongest and the weakest coordination with cations. However, the cation–ZW molecule coordination could slow the cationic diffusion. The competition of these effects resulted in accelerating or decelerating cationic diffusion. Our simulated results showed that ImCO2 enhanced Na + diffusion by 20%, while ChoPO4 and ImSO3 led to a 10% reduction. For K + , ChoPO4 reduced its diffusion by 40%, while ImCO2 and ImSO3 caused a similar decrease of 15%. These findings suggest that the ZW structure and the cationic size play an important role in the ionic dissociation effect of ZW materials.

25 ENERGY STORAGE↗
DOE OSTI2022

Infrared Spectroscopy of Li + Solvation in EmimBF 4 and in Propylene Carbonate: Ab Initio Molecular Dynamics and Experiment

Infrared (IR) spectra of solutions of the lithium salt LiBF 4 in the ionic liquid 1-ethyl-3-methylimidazolium tetrafluoroborate (EmimBF 4 ) and in the organic solvent propylene carbonate (PC) are studied via infrared spectroscopy and ab initio molecular dynamics (AIMD) simulations. The measurements show that the major effects of LiBF 4 in both solutions, compared to their neat counterparts, are the appearance of a new broad band in the 300-450 cm -1 frequency region and a broadening of the IR structure in the 900-1200 cm -1 region with the development of a new peak at 980 cm -1 . Computational analysis indicates that hindered translational motions of Li + in its solvation cage are mainly responsible for the former, while the latter is due to Li + -induced structural changes and accompanying vibrational frequency shifts of constituent ions and molecules of the solutions. In addition, molecular motions in these and lower frequency regions are generally correlated, disclosing the collective nature of the vibrational dynamics, which involve multiple ions/molecules. Herein, detailed analysis of these features via AIMD simulations of the spectrum and its components arising from auto- and cross-correlations of motions of constituent molecular species, combined with generalized normal modes of the solutions and normal modes of small Li + -containing clusters, is presented. Furthermore, other minor spectral changes caused by the lithium salt as well as the interaction-induced effect on IR spectra are also discussed.

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

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗
DOE OSTI2026

Isolating solvent–solute hydrogen bonding interactions via 2D IR solvation shell spectroscopy

The solvation shell around a solute is a fundamental feature of liquid-phase solutions, determining the behavior and properties of both the solute and the overall solution. Direct experimental measurements of the solvation shell properties are challenging due to the strong signals generated from the bulk solvent, which overwhelm the small contribution of the solvation shell. Here, we use ultrafast two dimensional infrared (2D IR) spectroscopy and intermolecular cross-peaks to isolate the IR absorption spectrum of methanol molecules in the solvation shell surrounding the solute N-methylacetamide. We demonstrate that the intermolecular coupling between the solvent and solute vibrations is indirectly mediated by a low-frequency hydrogen-bonding mode, suggesting an important mechanism for anharmonic coupling induced by hydrogen bonds. From the relative frequency shifts and cross-peak anisotropy, we find that methanol molecules surrounding N-methylacetamide form stronger and distinctly oriented hydrogen bonds than those in the bulk solvent. Here, we also compare these results with the solvent spectra of the solute N,N-dimethylacetamide to investigate how solute structural changes alter the solvation shell and the contribution of N–H hydrogen bond donation. Our results are supported by molecular dynamics simulations, which provide detailed insights into the hydrogen-bonding distributions. Through these results, we demonstrate 2D solvation shell spectroscopy to be a valuable method for investigating solvation structures and dynamics without interference from the bulk solvent.

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

Revealing Hot and Long-Lived Metastable Spin States in the Photoinduced Switching of Solvated Metallogrid Complexes with Femtosecond Optical and X-ray Spectroscopies

An atomistic understanding of the photoinduced spin-state switching (PSS) within polynuclear systems of d 4 –d 7 transition metal ion complexes is required for their rational integration into light-driven reactions of chemical and biological interest. However, in contrast to mononuclear systems, the multidimensional dynamics of the PSS in solvated molecular arrays have not yet been elucidated due to the expected complications associated with the connectivity between the metal centers and the strong interactions with the surroundings. Herein, the PSS in a solvated triiron(II) metallogrid complex is characterized using transient optical absorption and X-ray emission spectroscopies on the femtosecond time scale. The complementary measurements reveal the photoinduced creation of energy-rich (hot) and long-lived quintet states, whose dynamics differ critically from their mononuclear congeners. We report this finding opens major prospects for developing novel schemes in solution-phase spin chemistry that are driven by the dynamic PSS process in compact oligometallic arrays.

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

Molecular-level solvation and selectivity behavior of Na + , K + , and Li + within glycerol-derived solvents

Ether-based acyclic solvents have a tendency to confine positively charged atoms, forming structures called podands. Due to the stability of these structures, podands have been identified as promising candidates for several different electrochemical applications. In this work, we use molecular dynamics simulations to investigate the solvation and selectivity of Na + , K + , and Li + by different glycerol-derived solvents. The structural properties suggest that the solvent-cation association is stable over a very wide temperature range. Further, a strong correlation is observed between the more electronegative oxygen sites and preferential interaction with the cations, which suggests that the electrostatic oxygen-ion interactions dominate the coordination between these species. Significantly more intense interactions are observed between the solvents and Li + , versus Na + or K + . These molecular-level thermodynamic insights provide fundamental support for further exploring the performance of glycerol-derived podand complexes in Li + capture applications and/or electrolyte formulations.

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

Effect of Framework Composition and NH 3 on the Diffusion of Cu + in Cu-CHA Catalysts Predicted by Machine-Learning Accelerated Molecular Dynamics

Cu-exchanged zeolites rely on mobile solvated Cu + cations for their catalytic activity, but the role of the framework composition in transport is not fully understood. Ab initio molecular dynamics simulations can provide quantitative atomistic insight but are too computationally expensive to explore large length and time scales or diverse compositions. We report a machine-learning interatomic potential that accurately reproduces ab initio results and effectively generalizes to allow multinanosecond simulations of large supercells and diverse chemical compositions. Biased and unbiased simulations of [Cu(NH 3 ) 2 ] + mobility show that aluminum pairing in eight-membered rings accelerates local hopping and demonstrate that increased NH 3 concentration enhances long-range diffusion. The probability of finding two [Cu(NH 3 ) 2 ] + complexes in the same cage, which is key for SCR-NOx reaction, increases with Cu content and Al content but does not correlate with the long-range mobility of Cu + . Supporting experimental evidence was obtained from reactivity tests of Cu-CHA catalysts with a controlled chemical composition.

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

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗
DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗