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 OSTI2026

Dramatic Slowdown of Bimolecular Reactions of Criegee Intermediates in Organic Aerosols

Although Criegee intermediates (CIs) play a central role in driving the transformation of organic molecules, their condensed-phase reaction kinetics remains poorly understood. While CIs undergo unimolecular decomposition to produce OH radicals, stabilized CIs participate in a host of other bimolecular reactions forming larger molecular weight products, including secondary ozonides (kSOZ), alkoxyalkyl hydroperoxides (kAlAHP), and acyloxyalkyl hydroperoxides (kAcAHP). Here, we demonstrate, using prior literature and new measurements of heterogeneous organic aerosol reactions using O3 or OH radicals, that these condensed-phase bimolecular reaction rate coefficients kSOZ (10-20 cm3·molecule-1·s-1), kAlAHP, and kAcAHP (10-19 cm3·molecule-1·s-1) are 7-10 orders of magnitude slower than analogous gas-phase reactions (∼ 10-10-10-13 cm3·molecule-1·s-1). This contrasts with current evidence that the rate coefficients for unimolecular steps in the condensed phase are similar in magnitude to those in the gas phase. Kinetic models using gas-phase-like bimolecular coefficients consistently overestimate the importance of bimolecular relative to unimolecular pathways while underpredicting the aerosol reactivity by up to 30-fold. The slowing of bimolecular reactivity in aerosols is likely due to solvation effects and the increased steric hindrance of larger molecular weight CIs and their scavengers in the condensed phase. These results challenge the conventional view that the high reactant density of an organic aerosol should favor bimolecular over unimolecular reactions, underscoring the need for more accurate CI condensed-phase rate constants for reliable modeling of oxidative processes across a diverse set of atmospheric and environmental systems.

Zeng, Meirong↗
DOE OSTI2025

Ultrathin liquid sheets: water gets in shape for VUV absorption

We present absorption spectra of thin, free-flowing liquid sheets in the vacuum ultraviolet energy range using a gas-squeezed liquid jet. Compared to liquid flow cells, operation without transmission windows eliminates restrictions on the energy range. The temperature of the water sheet is estimated at 0 ± 3 °C, at the verge of the supercooled regime. By adjusting flow conditions in situ, we recorded absorption spectra at water sheet thicknesses ranging from 20 to 50 nm. We show that the absorption spectra of thin jets contain significant contributions from interference effects that need to be deconvoluted from spectral contributions due to the electronic structure. We employ a Fresnel propagation model to model the spectral changes and understand the impact of thickness variations and thin film interference. This opens the door for the investigation of solvation, interface, and similar effects by recording valence band spectra.

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

Influence of solution phase environmental heterogeneity and fluctuations on vibronic spectra: Perylene diimide molecular chromophore complexes in solution

Ensembles of ab initio parameterized Frenkel-exciton model Hamiltonians for different perylene diimide dimer systems are used, together with various dissipative quantum dynamics approaches, to study the influence of the solvation environment and fluctuations in chromophore relative orientation and packing on the vibronic spectra of two different dimer systems: a π-stacked dimer in aqueous solution in which the relative chromophore geometry is strongly confined by a phosphate bridge and a side-by-side dimer in dichloromethane involving a more flexible alkyne bridge that allows quasi-free rotation of the chromophores relative to one another. These entirely first-principles calculations are found to accurately reproduce the main features of the experimental absorption spectra, providing a detailed mechanistic understanding of how the structural fluctuations and environmental interactions influence the vibronic dynamics and spectroscopy of solutions of these multi-chromophore complexes.

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

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗
DOE OSTI2024

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

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

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

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

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