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

Using a Chromatographic Pseudophase Model To Elucidate the Mechanism of Olefin Separation by Silver(I) Ions in Ionic Liquids

Silver(I) ions undergo selective olefin complexation and have been utilized in various olefin/paraffin separation techniques such as argentation chromatography and facilitated transport membranes. Ionic liquids (ILs) are solvents known for their low vapor pressure, high thermal stability, low melting points, and ability to promote a favorable solvation environment for silver(I) ion–olefin interactions. To develop highly selective separation systems, a fundamental understanding of analyte partitioning to the stationary phase and the thermodynamic driving forces behind solvation is required. In this study, a chromatographic model treating silver(I) ions as a pseudophase is constructed and employed for the first time to investigate the olefin separation mechanism in silver(I) salt/IL mixtures. Stationary phases containing varying amounts of noncoordinated silver(I) salt ([Ag + ][NTf 2 – ]) dissolved in the 1-decyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C 10 MIM + ][NTf 2 – ]) IL are utilized to determine the partition coefficients of various analytes including alkanes, alkenes, alkynes, aromatics, aldehyde, esters, and ketones. As ligand coordination to silver(I) ions is known to lower its olefin complexation capability, this study also examines two different types of coordinated silver(I) ion pseudophases, namely, monocoordinated silver(I) salt ([Ag + (1-decyl-2-methylimidazole, DMIM)][NTf 2 – ]) and dicoordinated silver(I) salt ([Ag + (1-methylimidazole, MIM)(DMIM)][NTf 2 – ]). The extent of olefin partitioning to the coordinated silver(I) ion pseudophases over the carrier gas and IL decreased by up to two orders of magnitude. Values for enthalpy, entropy, and free energy of solvation were determined for the three silver(I) ion-containing systems. Olefin retention was observed to be enthalpically dominated, while ligand coordination to the silver(I) ion pseudophase resulted in variations for both enthalpic and entropic contributions to the free energy of solvation. Furthermore, the developed model can be used to study chemical changes that occur in silver(I) ions over time as well as identify optimal silver(I) salt/IL mixtures that yield high olefin selectivity.

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

The Significant Differences in Solvation Thermodynamics of C1–C3 Oxygenates in Hydrophilic versus Hydrophobic Pores of a Hydrophilic Ti-FAU Zeolite Model

The rates of catalytic reactions have been observed to be dramatically different in zeolites, depending on if they are hydrophobic or hydrophilic. Hypotheses aimed at explaining this behavior have pointed to various solvent molecule and zeolite properties as having influence on entropy. Herein, the influence of various solvent and adsorbate properties on the solvation energies, entropies, and free energies of eleven C1-C3 oxygenates in hydrophobic and hydrophilic pores within a hydrophilic model of Ti-FAU zeolite are tested. The results indicate significant variation in the calculated solvation thermodynamics depending on the adsorbate type, as well as if it is bound within a hydrophobic or hydrophilic pore. Further, while solvation energies are related to solvent-adsorbate interactions, solvation entropies have multiple contributions, and these differ depending on if the adsorbate is in a hydrophobic or hydrophilic pore. Specifically, solvation entropies in hydrophobic pores are related to solvent structural properties, whereas solvation entropies in hydrophilic pores are related to adsorbate polarity. Here, the large range of results obtained from two different pores within one zeolite model with minimal unique adsorption sites suggests that solvation behavior in zeolites is complicated and that the phenomena that control observed performance depend on the zeolite, reaction, and solvent.

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

Aqueous-phase effects on ethanol decomposition over Ru-based catalysts

The effects of an aqueous phase on ethanol decomposition for hydrogen production over a Ru(0001) catalyst surface model have been investigated from first principles. Solvent effects on the reaction mechanism and kinetic parameters have been quantified with the help of a microkinetic reactor model, density functional theory, and an implicit solvation scheme (iSMS). Our calculations indicate that in both vapor- and aqueous-phase reaction environments, the ethanol decomposition starts with acetaldehyde formation on the surface, some of which further dehydrogenates to ketenyl species (CHCO), where the C–C bond cleaves to form methylidyne (CH) and CO. In the vapor phase, adsorbed CH gets hydrogenated to methane, and CO desorbs or undergoes methanation reducing the amount of hydrogen produced. In contrast, under aqueous phase reaction conditions, the methanation is inhibited, and the water–gas shift (WGS) reaction is accelerated, leading to complete conversion of CO to CO 2 and H 2 . Furthermore, calculations indicate that the observed reaction behavior under aqueous phase reforming conditions originates primarily from the higher water chemical potential, and implicit solvent models predict only a small solvation effect.

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

Empirically Optimized One-Electron Pseudopotential for the Hydrated Electron: A Proof-of-Concept Study

Mixed quantum-classical molecular dynamics simulations have been important tools for studying the hydrated electron. They generally use a one-electron pseudopotential to describe the interactions of an electron with the water molecules. Furthermore, this approximation shows both the strength and weakness of the approach. On the one hand, it enables extensive statistical sampling and large system sizes that are not possible with more accurate ab initio molecular dynamics methods. On the other hand, there has (justifiably) been much debate about the ability of pseudopotentials to accurately and quantitatively describe the hydrated electron properties. These pseudopotentials have largely been derived by fitting them to ab initio calculations of an electron interacting with a single water molecule. In this paper, we present a proof-of-concept demonstration of an alternative approach in which the pseudopotential parameters are determined by optimizing them to reproduce key experimental properties. Specifically, we develop a new pseudopotential, using the existing TBOpt model as a starting point, which correctly describes the hydrated electron vertical detachment energy and radius of gyration. In addition to these properties, this empirically optimized model displays a significantly modified solvation structure, which improves, for example, the prediction of the partial molar volume.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

How Aqueous Solvation Impacts the Frequencies and Intensities of Infrared Absorption Bands in Flavin: The Quest for a Suitable Solvent Model

FTIR spectroscopy accompanied by quantum chemical simulations can reveal important information about molecular structure and intermolecular interactions in the condensed phase. Simulations typically account for the solvent either through cluster quantum mechanical (QM) models, polarizable continuum models (PCM), or hybrid quantum mechanical/molecular mechanical (QM/MM) models. Recently, we studied the effect of aqueous solvent interactions on the vibrational frequencies of lumiflavin, a minimal flavin model, using cluster QM and PCM models. Those models successfully reproduced the relative frequencies of four prominent stretching modes of flavin’s isoalloxazine ring in the diagnostic 1450–1750 cm−1 range but poorly reproduced the relative band intensities. Here, we extend our studies on this system and account for solvation through a series of increasingly sophisticated models. Only by combining elements of QM clusters, QM/MM, and PCM approaches do we obtain an improved agreement with the experiment. The study sheds light more generally on factors that can impact the computed frequencies and intensities of IR bands in solution.

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

Electronic Structure and Mechanical Properties of Solvated Montmorillonite Clay Using Large-Scale DFT Method

Montmorillonite clay (MMT) has been widely used in engineering and environmental applications as a landfill barrier and toxic waste repository due to its unique property as an expandable clay mineral that can absorb water easily. This absorption process rendered MMT to be highly exothermic due to electrostatic interactions among molecules and hydrogen bonds between surface atoms. A detailed study of a large supercell model of structural clay enables us to predict long-term nuclear waste storage. Herein, a large solvent MMT model with 4071 atoms is studied using ab initio density functional theory. The DFT calculation and analysis clarify the important issues, such as bond strength, solvation effect, elasticity, and seismic wave velocities. These results are compared to our previous study on crystalline MMT (dry). The solvated MMT has reduced shear modulus (G), bulk modulus (K), and Young’s modulus (E). We observe that the conduction band (CB) in the density of states (DOS) of solvated MMT model has a single, conspicuous peak at –8.5 eV. Moreover, the atom-resolved partial density of states (PDOS) summarizes the roles played by each atom in the DOS. These findings illuminate numerous potential sophisticated applications of MMT clay.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H 2 O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data., Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from experiment. Here, the An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Actinides↗
DOE OSTI2023

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

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

Cations Enhance Hydride Transfer to Noncatalytic Metals in Concentrated Alkaline Electrolytes

Alkali metal cations are known to influence the kinetics of the hydrogen evolution reaction (HER), acting as either promoters or inhibitors to the rate-determining step depending on the metal surface, local pH, and cation concentration. Despite the importance of concentrated electrolytes for commercial electrochemical cells, the impact of cations on the HER in concentrated alkaline environments (>1 M) and in mixed cation systems remains poorly understood. Here, this study quantifies the HER kinetics at polycrystalline metal surfaces (Pt, Au, Cu, and Fe) and the equilibrium solvation environment in pure and mixed alkali metal hydroxide electrolytes at concentrations up to 3.0 M. Kinetic analyses of Au, Cu, and Fe revealed a positive cation-concentration-effect that was primarily driven by changes to the charge transfer coefficient. Multinuclear NMR spectroscopy examined the solvation of H 2 O/OH – species and the alkali cations as a function of (mixed) alkali cation concentration(s), and demonstrated rapid exchange between solvent, hydroxide, and solvated cations. Together, these findings support models where HER kinetics on noncatalytic metal surfaces in strongly alkaline conditions are primarily governed by the average polarization and polarizability of the metal/solution interface and that increasing cation activity continues to increase the transfer coefficient at metal-hydroxide concentrations up to 3.0 M. Electrolytes and additives which can outcompete weakly hydrated cations and disrupt the interfacial water structure are expected to suppress parasitic HER at electrodes for energy storage and electroplating.

Cations↗
DOE OSTI2026

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2020

First-principles modeling of chemistry in mixed solvents: Where to go from here?

Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.

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

Highly Concentrated Electrolytes: Electrochemical and Physicochemical Characteristics of LiPF 6 in Propylene Carbonate Solutions

Highly concentrated electrolytes (HCEs) based on LiPF 6 in propylene carbonate (PC) have been examined as lithium-ion battery electrolytes. These HCEs have lower ionic conductivities and higher viscosities than ethylene carbonate (EC) electrolytes with 1.2 M LiPF 6 , but they have higher Li + ion transference numbers. Electrochemical cycling behaviour of LiNi 0.8 Co 0.015 Al 0.05 O 2 //graphite cells with 3.2 M LiPF 6 in PC resembles that of cells with EC-based electrolytes; the HCE cells have higher impedance, which can be lowered by increasing test temperature. By employing Raman and infrared spectroscopy, combined with density functional theory and ab initio molecular dynamics simulations, we reveal that the Li + solvation structure and speciation are key factors that determine cell performance. Two distinct regimes are observed as a function of salt concentration—in the conventional regime, the solvation number (SN) is mostly constant, while in the HCE regime it decreases linearly. Graphite exfoliation is suppressed only at very high salt concentrations ( > 2.4 M), where [ PC ] free /[Li + ] < 1 and P F 6 − f r e e > P C f r e e . Results from the Advanced Electrolyte Model indicate that Li + desolvation improves at higher LiPF 6 concentrations, thereby mitigating PC co-intercalation into the graphite. However, Li + ion transport is hindered in the HCEs, which increases impedance at both the oxide-positive and graphite-negative electrodes.

25 ENERGY STORAGE↗
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 OSTI2020

Charging-driven coarsening and melting of a colloidal nanoparticle monolayer at an ionic liquid–vacuum interface

In this work, we induce and investigate the coarsening and melting dynamics of an initially static nanoparticle colloidal monolayer at an ionic liquid–vacuum interface, driven by a focused, scanning electron beam. Coarsening occurs through grain interface migration and larger-scale motions such as grain rotations, often facilitated by sliding dislocations. The progressive decrease in area fraction that drives melting of the monolayer is explained using an electrowetting model whereby particles at the interface are solvated once their accumulating charge recruits sufficient counterions to subsume the particle. Subject to stochastic particle removal from the monolayer, melting is recapitulated in simulations with a Lennard-Jones potential. This new driving mechanism for colloidal systems, whose dynamical timescales we show can be controlled with the accelerating voltage, opens the possibility to manipulate particle interactions dynamically without need to vary particle intrinsic properties or surface treatments. Furthermore, the decrease in particle size availed by electron imaging presents opportunities to observe force and time scales in a lesser-explored regime intermediate between typical colloidal and molecular systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanical Properties of a Solvated Biomolecule: RGD (1FUV) Peptide

The mechanical properties of proteins/peptides play an essential role in their functionalities and implications, as well as their structure and dynamic properties. Understanding mechanical properties is pivotal to our knowledge of protein folding and the molecular basis of diverse cellular processes. Herein, we present a computational approach using ab initio quantum mechanical calculations to determine the mechanical properties—such as bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio—of a solvated Arg-Gly-Asp (RGD) peptide model. Since this peptide serves as the RGD-directed integrin recognition site and may participate in cellular adhesion, it is considered a promising small peptide for medicinal applications. This successful approach paves the way for investigating larger and more complex biomolecules.

Biochemistry & Molecular Biology↗