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

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2020

Moving beyond the Solvent-Tip Approximation to Determine Site-Specific Variations of Interfacial Water Structure through 3D Force Microscopy

Although interfacial solution structure impacts environmental, biological, and technological phenomena, including colloidal stability, protein assembly, heterogeneous nucleation, and water desalination, its molecular details remain poorly understood. Here, we visualize the three-dimensional (3D) hydration structure at the boehmite(010)-water interface using fast force mapping (FFM). Using a self-consistent scheme to decouple long-range tip-sample interactions from short-range solvation forces, we obtain the solution structure with lattice resolution. The results are benchmarked against molecular dynamics simulations that explicitly include the effects of the tip with different levels of approximation and systematically account for tip size, chemistry, and confinement effects. We find four laterally structured water layers within 1 nm of the surface, with the highest water densities at sites adjacent to hydroxyl groups. The key features beyond the first two layers can only be predicted using a full-scale simulation of the boehmite-water-silica system. Our findings further reveal a complex relationship between site-specific chemistry, water density, and long-range particle interactions; and present important advances toward quantitative data interpretation in 3D FFM.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reorganization Energy of Electron Transfer in Ionic Liquids

We report bandshape analysis of charge-transfer optical bands in room-temperature ionic liquids (ILs) is performed to extract reorganization energy of electron transfer. Remarkably, the reorganization energies in ILs are close to those in cyclohexane. This result runs against common wisdom in the field since conducting ILs, characterized by an infinite static dielectric constant, and nonpolar cyclohexane fall to the opposite ends of the polarity scale based on their dielectric constants. Theoretical calculations employing structure factors of ILs from molecular dynamics simulations support low values of the reorganization energy. Standard dielectric arguments do not apply to solvation in ILs and nonergodic reorganization energies are required for a quantitative analysis.

14 SOLAR ENERGY↗
DOE OSTI2025

Hydration and Restructuring of Polar Polymer Interfaces: Implications in Antifouling and Responsive Materials

Manipulating polymer interfaces is crucial for understanding how structure influences function in applications spanning biofouling prevention to energy storage. Moreover, observing how polymers adapt their microscopic structure to changes in their local environment can reveal essential properties that govern their performance in such applications, providing key insights into how to design more effective interfaces. Here, in this study, a series of “grafting-from” polymer brushes with side chains varying in charge, are probed by sum frequency generation (SFG) and modeled using all-atom molecular dynamics (MD) simulations to elucidate how chemical makeup and charge mediate interfacial restructuring in dry versus hydrated states. Results show that charge, in progressing from nonpolar to cationic to zwitterionic, results in dramatic changes in interfacial structure and overall hydration. While net neutral systems, regardless of bulk phase polarity, show minimal interfacial water structuring, the cationic species exhibits strong bulk water signals from the surface potential. Meanwhile, the polymer brushes themselves restructure in water differently independent of charge, impacting the functional groups that are presented to the aqueous phase. Nonpolar and cationic species for instance undergo a change in alkyl group orientations to accommodate hydrating water molecules, whereas the zwitterionic polymer becomes completely disordered in water. Overall, the structure-based behavior trends presented herein have implications in antifouling applications and responsive material interfaces.

biointerface↗
DOE OSTI2020

Exciton energy transfer reveals spectral signatures of excited states in clusters

Electronic excitation and concomitant energy transfer leading to Penning ionization in argon–acetylene clusters generated in a supersonic expansion are investigated with synchrotron-based photoionization mass spectrometry and electronic structure calculations. Spectral features in the photoionization efficiency of the mixed argon–acetylene clusters reveal a blue shift from the 2 P 1/2 and 2 P 3/2 excited states of atomic argon. Analysis of this feature suggests that excited states of argon clusters transfer energy to acetylene, resulting in its ionization and successive evaporation of argon. Theoretically calculated Ar n (n = 2–6) cluster spectra are in excellent agreement with experimental observations, and provide insight into the structure and ionization dynamics of the clusters. A comparison between argon–acetylene and argon–water clusters reveals that argon solvates water better, allowing for higher-order excitons and Rydberg states to be populated. Furthermore, these results are explained by theoretical calculations of respective binding energies and structures.

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

Hydration structure of flat and stepped MgO surfaces

We investigate the solvation structure of flat and stepped MgO(001) in neutral liquid water using ab initio molecular dynamics based on a hybrid density functional with dispersion corrections. Our simulations show that the MgO surface is covered by a densely packed layer of mixed intact and dissociated adsorbed water molecules in a planar arrangement with strong intermolecular H-bonds. The water dissociation fractions in this layer are >20% and >30% on the flat and stepped surfaces, respectively. Slightly above the first water layer, we observe metastable OH groups perpendicular to the interface, similar to those reported in low temperature studies of water monolayers on MgO. Here, these species receive hydrogen bonds from four nearby water molecules in the first layer and have their hydrophobic H end directed toward bulk water, while their associated protons are bound to surface oxygens. The formation of these OH species is attributed to the strong basicity of the MgO surface and can be relevant for understanding various phenomena from morphology evolution and growth of (nano)crystalline MgO particles to heterogeneous catalysis.

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

Direct Correlation of the Salt-Reduced Diffusivities of Organic Solvents with the Solvent’s Mole Fraction

We report that lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in organic solvents (especially propylene carbonate) has demonstrated extraordinary pseudocapacitive performance as an electrolyte in the supercapacitor configuration. However, the influence of the solvated ions on the diffusivity of the solvent molecules is yet to be understood. We examine the impact of LiTFSI on the diffusivity in five organic solvents: acetonitrile (ACN), tetrahydrofuran (THF), methanol (MeOH), dimethyl sulfoxide (DMSO), and propylene carbonate (PC) using a combination of neutron scattering, conductivity measurements, and molecular dynamics simulations. The extent of the diffusivity reduction in the concentration regime of ≤1 M directly correlates with the solvent mole fraction at which the solvation shells around Li + ions are of similar size in all the solvents, resulting in a universal ~50% reduction in the solvent diffusivity. These results provide guidance for formulation of the new electrolytes to enhance the performance of energy storage devices.

conductivity↗
DOE OSTI2023

Changes in nitrate binding with lanthanides in BLPhen complexes

Preorganized ligands such as bis-lactam-1,10-phenanthroline (BLPhen) offer unique selectivity trends in solvent extraction of rare-earth ions, particularly from aqueous nitrate solutions. However, due to the ligand’s lipophilicity it is experimentally challenging to obtain crystal structures of the corresponding lanthanide (Ln):ligand complexes to describe the local structure of the first coordination shell and the role of the nitrate ions in the organic phase. In this work, we report first principles molecular dynamics (FPMD) simulations of the complexation of trivalent lanthanide ions (La, Nd, and Eu) with the BLPhen ligand in the presence of nitrate ions in the dichloroethane (DCE) solvent. We find that two nitrate anions are present in the first solvation shell and the third nitrate anion is far from the first solvation shell in all three [Ln(BLPhen) 2 ] 3+ complexes examined. Moving along the lanthanide series from La to Nd to Eu, the binding pocket formed by the two BLPhen ligands shrinks in size; as a result, the two nitrate configuration changes from one-monodentate/one-bidentate binding mode for La and Nd to all monodentate for Eu. Meanwhile, the total coordination number drops from 11 for La and Nd to 10 for Eu. More interestingly, the bidentate binding in [La(BLPhen) 2 ] 3+ and [Nd(BLPhen) 2 ] 3+ complexes is highly dynamic, frequently switching to monodentate and back. The FPMD insights into the differing modes and dynamics of nitrates in the first coordination shell will be useful for further atomistic understanding of the complex structure, formation, and stability in the organic phase.

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

Adsorption and Solvation Modulate Rhodamine B Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Confinement of solvents and solutes within nanoporous materials frequently leads to the emergence of unique mass transport behaviors that, once fully understood, may lead to improved chemical separations. Here, the diffusion of Rhodamine B (RhB) dye within 10 and 20 nm diameter anodic aluminum oxide (AAO) nanopores filled with binary ethanol/water mixtures is investigated. Mixture compositions spanning from pure ethanol to pure water are employed. The results of confocal fluorescence correlation spectroscopy studies reveal that RhB diffusion occurs by a two-component mechanism comprising composition-dependent fast and slow motions, characterized by diffusion coefficients D f and D s . The results are consistent with those of previous studies performed under more limited conditions [J. Phys. Chem. C, 2023, 127, 411-420]. The fast component scales with mixture viscosity and is assigned to hindered bulk-like diffusion in central pore regions. Slow diffusion likely involves adsorption of RhB to the pore surface and may be described by a desorption mediated mechanism. The occurrence of RhB adsorption to the AAO surface is verified at the single-molecule level by wide-field fluorescence imaging of membrane cross-sectional surfaces. Unique composition dependent trends in the autocorrelation amplitude and in D s that mimic bulk RhB solubility are revealed. D s is found to be smallest in pure ethanol and pure water and largest in intermediate mixtures. These results suggest that RhB surface adsorption is strongest in the pure liquids and weakest in mixtures of intermediate composition, where the dye is least soluble, and most soluble, respectively. As a result, molecular dynamics simulations reveal that a water layer appears on the pore surface under most conditions, while RhB is solvated primarily by ethanol. The composition dependence of RhB diffusion is concluded to reflect its solvation dependent interactions with the pore walls.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding extreme fast charge limitations in carbonate mixtures

Charging large format lithium ion batteries within ten to fifteen minutes requires changes to the electrolyte composition in addition to modification of electrode and cell architectures. Several approaches to address this need have been proposed; but there is not a lot of clarity on understanding which factors limit the performance of existing electrolytes. Furthermore, this work takes a closer look at the solvated components of a mixture of LiPF 6 in ethylene carbonate (EC) and ethyl methyl carbonate (EMC) in the context of extreme fast charging and relates these findings to cell-level requirements. Molecular dynamics studies of the Gen-2 electrolyte compositions with increasing salt concentrations, have been performed to estimate transport properties like diffusivity, transference number and conductivity. Molecular-level differences in the structure of solvation shells under extreme LiPF 6 concentrations are probed here and some key aspects on solvent structure that help overcome barriers to Li + transport under extreme fast charge are discussed.

25 ENERGY STORAGE↗
DOE OSTI2024

High‐Rate Polymeric Redox in MXene‐Based Superlattice‐Like Heterostructure for Ammonium Ion Storage

Abstract Achieving both high redox activity and rapid ion transport is a critical and pervasive challenge in electrochemical energy storage applications. This challenge is significantly magnified when using large‐sized charge carriers, such as the sustainable ammonium ion (NH 4 + ). A self‐assembled MXene/n‐type conjugated polyelectrolyte (CPE) superlattice‐like heterostructure that enables redox‐active, fast, and reversible ammonium storage is reported. The superlattice‐like structure persists as the CPE:MXene ratio increases, accompanied by a linear increase in the interlayer spacing of MXene flakes and a greater overlap of CPEs. Concurrently, the redox activity per unit of CPE unexpectedly intensifies, a phenomenon that can be explained by the enhanced de‐solvation of ammonium due to the increased volume of 3 Å‐sized pores, as indicated by molecular dynamic simulations. At the maximum CPE mass loading (MXene:CPE ratio = 2:1), the heterostructure demonstrates the strongest polymeric redox activity with a high ammonium storage capacity of 126.1 C g −1 and a superior rate capability at 10 A g −1 . This work unveils an effective strategy for designing tunable superlattice‐like heterostructures to enhance redox activity and achieve rapid charge transfer for ions beyond lithium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tuning conformational structures of imidazolium ionenes with 1-ethyl-3-methylimidazolium ionic liquid solvents

In this report the structures and conformations of imidazolium ionenes solvated in 1-ethyl-3-methylimidazolium ionic liquids (ILs) ([C 2 mim + ][X - ]) are investigated using molecular dynamics simulations. Four different ionenes poly(decylimidazolium) (PD 10 ), poly(tetraethyleneglycolimidazolium) (PE 10 ), alternating copolymer (P(ED) 5 ), and a block copolymer PE 5 D 5 are immersed in [C 2 mim + ][X - ], where [X - ] = thiocyanate [SCN - ], tetrafluoroborate [BF 4 - ], and trifluoromethanesulfonate [TfO - ]. The radius of gyration indicates that the ionene chains are more extended when immersed in [TfO - ] compared to [BF 4 - ] and [SCN - ], while the solvent accessible surface area shows that PE 10 is the most contracted ionene among all systems. The electrostatic interactions between the ionenes and the ILs are very consistent with the inherent electrostatic characteristics of the different anions, which can be quantified using the Ionic Polarity Index (IPI). The ionene configurations and dynamics strongly depend on the backbone functionality and architecture of the chain, due to the balance between the inter- and intra-molecular interactions in these systems.

42 ENGINEERING↗
NASA NTRS2006

Hybrid MD-Nernst Planck Model of Alpha-hemolysin Conductance Properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an alpha-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson-Nemst-Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K(+) and Cl(-)) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1 M KCI solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5-7 times in comparison to bulk values. Significant statistical variations (17-45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius approx. 9A with two constriction blocks where the radius is reduced to approx. 6A. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the a-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.

Cozmuta, Ioana↗
DOE OSTI2026

Operando Depth-Resolved Measurement of Solvation Entropy, Interfacial Transport, and Charge-Transfer Kinetics in Lithium-Ion Batteries

Understanding and improving the performance and longevity of lithium-ion batteries critically depends on insight into the dynamic processes occurring at buried electrode-electrolyte interfaces. However, direct, depth-resolved, and operando diagnosis of these interfaces remains a longstanding challenge due to their inaccessibility beneath bulk materials, the limitations of conventional surface- and bulk-sensitive characterization tools, and the difficulty of maintaining realistic cell environments during measurement. These challenges have made it nearly impossible to uniquely resolve important interfacial properties such as charge transfer resistance, SEI (solid electrolyte interphase) resistance, and solvation entropy at the individual electrode interfaces within a working cell, information that is essential for mechanistic insight and accelerated battery design. Here, in this study, we report the development of multiharmonic electro-thermal spectroscopy (METS), an operando technique that enables depth-resolved measurement of solvation entropy, interfacial transport resistance, charge-transfer resistance, and SEI resistance at individual electrode-electrolyte interfaces within practical lithium-ion batteries. By leveraging frequency-dependent, thermal-wave sensing and interface-specific modeling, METS uniquely attributes interfacial properties to specific electrodes, as validated by comparison with traditional electrochemical impedance spectroscopy (EIS). The ability to spatially and temporally resolve interfacial processes in real time provides new diagnostic capabilities that are crucial for mechanistic studies of battery degradation and for the rapid development of next-generation energy storage systems.

Chalise, Divya [University of California, Berkeley↗
DOE OSTI2025

Solvation directed morphological control in metal oxide nanostructures

The development of structural hierarchy on various length scales during the crystallization process is ubiquitous in biological systems and minerals and is common in synthetic nanomaterials. The driving forces for the formation of complex architectures range from local interfacial interactions, that modify interfacial speciation, local supersaturation, and nucleation barriers, to macroscopic interparticle forces. Although it is enticing to interpret the formation of hierarchical architectures as the assembly of independently nucleated building blocks, crystallization pathways often follow monomer-by-monomer addition with structural complexity arising from interfacial chemical coupling and strongly correlated fluctuation dynamics in the electric double layers. Here, we show that the development of structural hierarchy through heterogeneous nucleation is driven by dipolar and solvation forces. Specifically, coupled simulations and experimental studies revealed that dipole build-up along the slow growth direction can trigger twinning and the development of branched architectures. Enthalpic solvation interactions were shown to either enhance or reduce the dipole moment of the nanoparticles and, thereby, control crystal morphology and architecture. The systematic studies of chemical coupling between different solvents and undercoordinated surface atoms of the growing nanocrystals revealed the mechanism of dimensionality control and the development of structural hierarchy without ligands or structure-directing agents.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Solution Structures of Europium Terpyridyl Complexes with Nitrate and Triflate Counterions in Acetonitrile

Lanthanide–ligand complexes are key components of technological applications, and their properties depend on their structures in the solution phase, which are challenging to resolve experimentally or computationally. The coordination structure of the Eu 3+ ion in different coordination environments in acetonitrile is examined using ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy. AIMD simulations are conducted for the solvated Eu 3+ ion in acetonitrile, both with or without a terpyridyl ligand, and in the presence of either triflate or nitrate counterions. EXAFS spectra are calculated directly from AIMD simulations and then compared to experimentally measured EXAFS spectra. In acetonitrile solution, both nitrate and triflate anions are shown to coordinate directly to the Eu 3+ ion forming either ten- or eight-coordinate solvent complexes where the counterions are binding as bidentate or monodentate structures, respectively. Coordination of a terpyridyl ligand to the Eu 3+ ion limits the available binding sites for the solvent and anions. In certain cases, the terpyridyl ligand excludes any solvent binding and limits the number of coordinated anions. The solution structure of the Eu-terpyridyl complex with nitrate counterions is shown to have a similar arrangement of Eu 3+ coordinating molecules as the crystal structure. Finally, this study illustrates how a combination of AIMD and EXAFS can be used to determine how ligands, solvent, and counterions coordinate with the lanthanide ions in solution.

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

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗