Search NASA⌕ Search

SEARCH · Search NASA

Results for “solvation models”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Understanding Self-Assembly and the Stabilization of Liquid/Liquid Interfaces: The Importance of Ligand Tail Branching and Oil-Phase Solvation

Organophosphorus-based ligands represent a versatile set of solvent extraction reagents whose chemical makeup plays an important role in extraction mechanism. We hypothesize that the branching of the extractant hydrophobic tail and its oil-phase solvation affect the liquid/liquid interfacial structure. Understanding the structure mediated adsorption and interfacial ordering becomes key in designing ligands with enhanced selectivity and efficiency for targeted extractions. We employed vibrational sum frequency generation spectroscopy and interfacial tension measurements to extract thermodynamic adsorption energies, map interfacial ordering, and rationalize disparate behaviors of model di-(2-ethylhexyl) phosphoric acid and dioctyl phosphoric acid ligands at the hexadecane water interface. With increased surface loading, ligands with branched hydrophobic tails formed stable interfaces at much lower concentrations than those observed for ligands with linear alkyl tails. The lack of an oil phase and associated solvation results in markedly different interfacial properties, and thus measurements made at air/liquid surfaces cannot be assumed to correlate with the processes occurring at buried liquid/liquid interfaces. Here, we attribute these differences in the surface mediated self-assembly to key variations in hydrophobic interactions and tail solvation taking place in the oil phase demonstrating that interactions in both the polar and nonpolar phases are essential to understand self-assembly and function.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Relation between the Hydrated Electron Solvation Structure and Its Partial Molar Volume

It is now generally accepted that the hydrated electron occupies a cavity in water, but the size of the cavity and the arrangements of the solvating water molecules are not fully characterized. Here, we use the Kirkwood-Buff (KB) approach to examine how the partial molar volume (V M ) provides insight into these issues. The KB method relates V M to an integral of the electron-water radial distribution function, a key measure of the hydrated electron structure. Here we have applied it to three widely-used pseudopotentials and the results show that V M is a sensitive measure of the fidelity of hydrated electron descriptions. Thus, the measured V M places constraints on the hydrated electron structure that are important in developing and evaluating model descriptions. Importantly, we find that V M does not reflect only the cavity size (and thus should not be used to infer the cavity radius), but is strongly dependent on the extended solvation structure.

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

Locking water molecules via ternary O–H⋯O intramolecular hydrogen bonds in perhydroxylated closo -dodecaborate

A multitude of applications related to perhydroxylated closo-dodecaborate B 12 (OH) 12 2− in the condensed phase are inseparable from the fundamental mechanisms underlying the high water orientation selectivity based on the base B 12 (OH) 12 2− . Herein, we directly compare the structural evolution of water clusters, ranging from monomer to hexamer, oriented by functional groups in the bases B 12 H 12 2− , B 12 H 11 OH 2− and B 12 (OH) 12 2− using multiple theoretical methods. A significant revelation is made regarding B 12 (OH) 12 2− : each additional water molecule is locked into the intramolecular hydrogen bond B–O–H ternary ring in an embedded form. This new pattern of water cluster growth suggests that B–(H–O)⋯H–O interactions prevail over the competition from water–hydrogen bonds (O⋯H–O), distinguishing it from the behavior observed in B 12 H 12 2− and B 12 H 11 OH 2− bases, in which competition arises from a mixed competing model involving dihydrogen bonds (B–H⋯H–O), conventional hydrogen bonds (B–(H–O)⋯H–O) and water hydrogen bonds (O⋯H–O). Through aqueous solvation and ab initio molecular dynamics analysis, we further demonstrate the largest water clusters in the first hydrated shell with exceptional thermodynamic stability around B 12 (OH) 12 2− . These findings provide a solid scientific foundation for the design of boron cluster chemistry incorporating hydroxyl-group-modified borate salts with potential implications for various applications.

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

Mean-field theory for the structure of strongly interacting active liquids

Active systems, which are driven out of equilibrium by local non-conservative forces, exhibit unique behaviors and structures with potential utility for the design of novel materials. An important and difficult challenge along the path toward this goal is to precisely predict how the structure of active systems is modified as their driving forces push them out of equilibrium. Here, we use tools from liquid-state theories to approach this challenge for a classic minimal active matter model. First, we construct a nonequilibrium mean-field framework that can predict the structure of systems of weakly interacting particles. Second, motivated by equilibrium solvation theories, we modify this theory to extend it with surprisingly high accuracy to systems of strongly interacting particles, distinguishing it from most existing similarly tractable approaches. Furthermore, our results provide insight into spatial organization in strongly interacting out-of-equilibrium systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
DOE OSTI2023

Uncovering the 3d and 4d Electronic Interactions in Solvated Ru Complexes with 2p3d Resonant Inelastic X-ray Scattering

The electronic structure and dynamics of ruthenium complexes are widely studied given their use in catalytic and light-harvesting materials. Here we investigate three model Ru complexes, [Ru III (NH 3 ) 6 ] 3+ , [Ru II (bpy) 3 ] 2+ , and [Ru II (CN) 6 ] 4– , with L 3 -edge 2p3d resonant inelastic X-ray scattering (RIXS) to probe unoccupied 4d valence orbitals and occupied 3d orbitals and to gain insight into the interactions between these levels. In this study, the 2p3d RIXS maps contain a higher level of spectral information than the L 3 X-ray absorption near edge structure (XANES). This study provides a direct measure of the 3d spin–orbit splittings of 4.3, 4.0, and 4.1 eV between the 3d 5/2 and 3d 3/2 orbitals of the [Ru III (NH 3 ) 6 ] 3+ , [Ru II (bpy) 3 ] 2+ , and [Ru II (CN) 6 ] 4– complexes, respectively.

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

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗
DOE OSTI2022

Interfacial Water Manipulation with Ionic Liquids for the Oxygen Reduction Reaction

The role that interfacial water plays in both promoting and inhibiting the electrochemical oxygen reduction reaction (ORR) remains somewhat unclear and controversial. Here, we use electroanalytical chemistry, spectroscopy, and microkinetic modeling to probe the impact of hydrophobic ionic liquid (IL) thin films on interfacial water structure and its role in promoting or inhibiting the ORR. Through the use of in situ ATR-SEIRAS, we find that the IL thin films limit the content of water at the interface and prevent the formation of hydrogen bond stabilized water organization. The impact of this exclusion of water on the promotion of the ORR is through the reduced coverage of OH ad spectator species. The decreased solvation of “active” OH ad species weakens its interaction with the catalyst surface, lowering the barrier to the last step in the ORR mechanism. This “destabilized” OH ad impact on ORR kinetics is confirmed through a microkinetic model. Furthermore, the results presented here highlight the mechanistic pathway through which hydrophobic ILs enhance ORR kinetics and point to pathways to both further improve this performance enhancement and the potential for integration of hydrophobic ILs into other technologically relevant elementary electrochemical reactions.

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 OSTI2022

Computational Modeling of Battery Materials

This chapter provides a broad overview of the various computational modeling techniques used to gain fundamental insights into coupled electrochemical processes that occur in battery materials at electronic-to-mesoscopic scales. Furthermore, representative successes of these techniques in modeling electrodes, electrolytes, and electrode-electrolyte interfaces are highlighted to establish the current state-of-the-art in the field.

25 ENERGY STORAGE↗
DOE OSTI2022

Imaging Wavepackets in Real-Space & Disentangling Ultrafast Solvation Dynamics with High-Energy Ultrafast X-ray Scattering

The microscopic information of solute and solute-solvent motions can be measured using ultrafast diffuse x-ray scattering, and compared to molecular dynamics simulations, as demonstrated in studies of solvated molecular systems done in X-ray free-electron lasers. However, for the typical photon energies used in these experiments, the photon momentum transfer range was limited to lower Q values where the signals arising from different parts of the studied system overlap significantly. Recently, this limitation was greatly alleviated by extending the photon energy at the Linac Coherent Light-Source (LCLS), significantly increasing the accessible momentum transfer range. This improvement is transformative for ultrafast diffuse x-ray scattering measurements, enabling for the first time to disentangle the solute scattering difference signal, that persists to higher Q values, from the bulk solvent and solute-solvent cross-terms difference signals. In addition, it enhances the prediction ability of modeling and simulation methods such as the hybrid QM/MM approach. The extended Q range also opens the way to resolve in real-space details regarding coherent wavepackets motions beyond their center-of-mass positions. Here, we present the first results on high-energy (18keV) diffuse x-ray scattering in solution, demonstrating high fidelity time-resolved scattering and analysis of the photoexcited model photocatalysts PtPOP ([Pt2(POP)4]4-), and IrDimen ([Ir2(dimen)4]2+) in several solvents. These complexes provide ideal systems for demonstrating the ability of high-Q ultrafast scattering and QM/MM simulations to decompose ultrafast solvation dynamics into specific changes in the solute-solvent pair distribution function with a particular focus on how electronically excited states change the interaction between the solvent and photo-catalytically active metal sites. We introduce an approach to further utilize the high-energy capability and demonstrate a single-shot ultra-wide-angle X-ray Scattering modality using two perpendicular detectors, spanning a scattering angle range of more than 100 degrees, allowing to extend the accessible momentum transfer range up to Q~14 Å-1. We develop a model-free method to invert the x-ray scattering signals and enable the recovery of multiple pair-density motions that happen simultaneously, allowing us to directly measure in real-space nuclear wavepacket motions. . [1] Natan, Adi. "Real-Space Inversion and Super-Resolution of Ultrafast X-ray Scattering using Natural Scattering Kernels." arXiv preprint arXiv:2107.05576 (2021)

Natan, Adi↗
DOE OSTI2021

Mechanistic Study of the Li–Air Battery with a Co 3 O 4 Cathode and Dimethyl Sulfoxide Electrolyte

The lithium–air battery, a powerful competitor to replace the traditional lithium-ion battery, has attracted increasing attention due to its extremely high theoretical energy density. However, its development is limited by the cathode and electrolyte properties, which should include high stability, conductivity, and electrocatalytic properties in oxygen-rich environments. Here, we employ a systematic first-principles study of Li–O 2 discharge and charge reactions on the Co 3 O 4 -based cathode with the assistance of dimethyl sulfoxide (DMSO) electrolyte. The structure, stability, and electronic properties of different surface reconstructions of the Co 3 O 4 (100) facet are investigated. In addition, the mechanisms and thermodynamic overpotentials of multi-step reactions between Li + /e – and O 2 are provided, where lithium suboxide products (Li 2 O 2 or Li 3 O 2 ) are formed on the different Co 3 O 4 (100) terminations. The solvation shell of Li + components in explicit DMSO solvent is investigated through ab initio molecular dynamics simulations. In general, we find that the Co 3 O 4 (100)-O (oxidized) surface is the most stable one under standard conditions, and the stable Li + solvation structure is found in a tetrahedral Li(DMSO) 4 + shell in the DMSO-based electrolyte. Moreover, in the system of the Co 3 O 4 (100)-O cathode and DMSO electrolyte, the solution model pathway is energetically favorable for the Li–O 2 discharge reaction. It provides a low constant overpotential of 0.17 V during a long-term discharging process, thus causing the final toroid Li 2 O 2 formation on the cathode. During the charging process, an overpotential of 0.36 V is required to rapidly decompose Li 2 O 2 .

25 ENERGY STORAGE↗
DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
DOE OSTI2023

A redox model for NO oxidation, NH 3 oxidation and high temperature standard SCR over Cu-SSZ-13

A kinetic model is developed to predict the influence of temperature and hydrothermal aging on the redox of active Cu sites under standard SCR, NO oxidation and NH 3 oxidation conditions over a practically relevant fully-formulated Cu-SSZ-13 catalyst. NO 2 /N 2 production during NO/NH 3 titration of Cu II sites is utilized to identify rate parameters associated with NO-only RHC (reduction half cycle) and NH 3 -only RHC respectively. Integral N 2 formation during subsequent NO + NH 3 titration is consistent with the production of one NO 2 per two CuII sites reduced during NO-only RHC and one N 2 per six Cu II sites reduced during NH 3 -only RHC. Decreased reduction of Cu II sites by NO/NH 3 upon hydrothermal aging, along with the production of one NO 2 per two CuII sites during NO-only RHC, is accordant with the involvement of proximal ZCuOH and oxygen-bridged dimeric Cu II sites. Oxidation of partially solvated and framework coordinated ZCu (Cu I ) sites occurs in presence of O 2 , does not produce N 2 and can lead to the consumption of Brønsted acid sites. A global OHC kinetic model is developed to predict transient and integral N 2 formation during exposure of CuI sites to a mixture of NO and O 2 . The resulting redox kinetic model quantitatively predicts NO and NH 3 consumption during isothermal transient response Cu redox (TRCR) protocols, along with temperature and age dependent steady-state standard SCR and oxidation conditions. The redox model presented in this work synthesizes recent kinetic, spectroscopic and computational findings to provide a foundational description of active site redox during standard SCR, NO oxidation and NH3 oxidation over Cu-SSZ-13.

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

Microscopic insights into the solvation of polyethylene glycol chains in water: A machine learning potential approach

Polyethylene glycol (PEG) is a structurally simple, nontoxic, and water-soluble polymer widely utilized in medical and pharmaceutical applications. Notably, when a PEG chain is immersed in water, the surrounding water molecules play a key role in driving conformational changes of this macromolecule. In this study, we explore the solvation behavior of PEG under mechanical strain using molecular dynamics simulations, with an interatomic potential obtained from machine learning. Our focus is on the transition from the favored coil-like conformation to an extended one under external force. Through analyses of radial distribution functions, hydrogen bonding, and solvation dynamics, we uncover how mechanical stretching influences the local hydration environment. Furthermore, we disentangle the enthalpic and entropic contributions to the conformational stability of PEG in water. Surprisingly, our neural network potential model identifies dewetting of PEG C-atoms, and not water H-bonding with PEG O-atoms, as the main enthalpic driving force for the coiling of PEG in water.

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↗
DOE OSTI2024

Zn 2+ -mediated catalysis for fast-charging aqueous Zn-ion batteries

Rechargeable aqueous zinc-ion batteries (AZIBs), renowned for their safety, high energy density and rapid charging, are prime choices for grid-scale energy storage. Historically, ion-shuttling models centring on ion-migration behaviour have dominated explanations for charge/discharge processes in aqueous batteries, like classical ion insertion/extraction and pseudocapacitance mechanisms. However, these models struggle to account for the exceptional performance of AZIBs compared to other aqueous metal-ion batteries. Here, in this study, we present a catalysis model elucidating the Zn 2+ anomaly in aqueous batteries, explaining it through the concept of adsorption in catalysis. Such behaviour can serve the charge/discharge role, predominantly dictated by solvated metal cations and cathode materials. First-principles calculations suggest optimal adsorption/desorption behaviour (water dissociation process) with the Zn 2+ -vanadium nitride (VN) combination. Experimentally, AZIBs implementing VN cathodes demonstrate fast-charging kinetics, showing a capacity of 577.1 mAh g -1 at a current density of 300,000 mA g -1 . The grasp of catalysis steps within AZIBs can drive solutions beyond state-of-the-art fast-charging batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Unexpected electrochemical behavior of an anolyte redoxmer in flow battery electrolytes: solvating cations help to fight against the thermodynamic–kinetic dilemma

Redoxmers are redox-active molecules that can store energy in electrolytes for redox flow batteries (RFBs), and their electrochemical properties are significantly affected by the choice of supporting electrolytes. Herein, we use 2,1,3-benzothiadiazole (BzNSN) as a model system to scrutinize the supporting electrolyte impact. By systemically varying the components of supporting salts, BzNSN not only shows substantial redox potential shifts but also exhibits varying electrochemical stabilities. Specifically, changing the size of cations can effectively alter the coordination between the supporting salt and BzNSN species. From Li + , Na + , K + , to NEt 4 + , the redox potential of BzNSN shifts negatively, from -1.63 V to -1.82 V vs. Ag/Ag + . Molecular dynamics and density functional theory simulations revealed that smaller cations, like Li + , are closer to the charged BzNSN when coordinated, implying stronger coordination, while larger cations, like K + and NEt 4 + , are farther away. Interestingly, the large cation electrolytes also lead to much improved electrochemical stability, evidenced by the extraordinarily enhanced kinetic lifetime from electron paramagnetic resonance measurement. This study demonstrates the first example of tuning an anolyte redoxmer toward a concurrent improvement of lowered redox potentials AND enhanced calendar lives via solvation means, which is usually constrained by the thermodynamic–kinetic relation.

2,1,3-benzothiadiazole↗